这道题卡dij？就卡一个点

#include<cstdio>
#include<queue>
#include<cstring>
using namespace std;
const int maxn=6000;
const int inf=0x7f7f7f7f;
const int mod=1000000007;

struct node
{
    int to,next,w,id;
}edge[maxn],edge1[maxn];
struct point
{
    int id,w;
    point(int a,int b)
    {
        id=a;
        w=b;
    }
    bool friend operator < (point a,point b)
    {
        return a.w>b.w;
    }
};
priority_queue<point > q;
int k,n,m,top,k1;
int dis[maxn],head[maxn],head1[maxn],vis[maxn],topz[maxn],sta[maxn],topf[maxn],ans[maxn];
void add(int u,int v,int w)
{
    edge[++k].next=head[u];
    edge[k].to=v;
    edge[k].w=w;
    head[u]=k;
}
void add1(int u,int v,int w,int c)
{
    edge1[++k1].next=head1[u];
    edge1[k1].to=v;
    edge1[k1].w=w;
    head1[u]=k1;
    edge1[k1].id=c;
}
void dij(int s)
{
    top=0;
    for(int i=1;i<=n;i++) dis[i]=inf;
    memset(topz,0,sizeof(topz));
    memset(vis,0,sizeof(vis));
    dis[s]=0;
    topz[s]=1;
    q.push(point(s,0));
    while(!q.empty())
    {
        int u=q.top().id;
        q.pop();
        if(vis[u]) continue ;
        vis[u]=1;
        sta[++top]=u;
        for(int i=head[u];i;i=edge[i].next)
        {
            int v=edge[i].to;
            if(dis[v]==dis[u]+edge[i].w) topz[v]=(topz[u]+topz[v])%mod;
            if(dis[v]>dis[u]+edge[i].w)
            {
                dis[v]=dis[u]+edge[i].w;
                topz[v]=topz[u];
                q.push(point(v,dis[v]));
            }
        }
    }
}
void topu()
{
    memset(topf,0,sizeof(topf));
    while(top)
    {
        topf[sta[top]]++;
        for(int i=head1[sta[top]];i;i=edge1[i].next)
        {
            int v=edge1[i].to;
            if(dis[sta[top]]==dis[v]+edge1[i].w) topf[v]=(topf[v]+topf[sta[top]])%mod,ans[edge1[i].id]=(ans[edge1[i].id]+topz[v]*topf[sta[top]])%mod;
        }
        top--;
    }
}
int main()
{
    int u,v,w;
    scanf("%d%d",&n,&m);
    for(int i=1;i<=m;i++)
    {
        scanf("%d%d%d",&u,&v,&w);
        add(u,v,w);
        add1(v,u,w,i);
    }
    for(int i=1;i<=n;i++)
    {
        dij(i);
        topu();
    }
    for(int i=1;i<=m;i++) printf("%d\n",ans[i]);
    return 0;
}