zoj 3408 dp

做法：先记录最短路径,然后两次dfs求出到当前点的方法总数,和可以到的总数,相乘得出答案.

由于其实最后得出的路径是不存在循环的，可以dfs

(真是气炸了 函数返回写int ，错了一下午 (\ _/)！！ ) 

(学到 布思乘法处理 10^10)

#include <iostream>
#include <stdio.h>
#include <string.h>
#include <algorithm>
#include <queue>
using namespace std;
#define MAXA 25005
#define MAXN 50005
#define LL long long
#define MOD 10000000000

LL dis[MAXA];

LL dp1[MAXA];
LL dp2[MAXA];

LL ans[MAXA];
bool vis[MAXA];

struct node
{
    int u;
    int v;
    int c;
    int p;
};
node r[MAXN];
node r0[MAXN];

int head[MAXA];
int head0[MAXA];
int cnt;
int cnt0;

void addedge(int u,int v,int c)
{
    r[cnt].p=head[u];
    r[cnt].u=u;
    r[cnt].v=v;
    r[cnt].c=c;
    head[u]=cnt++;

    r0[cnt0].p=head0[v];
    r0[cnt0].u=v;
    r0[cnt0].v=u;
    r0[cnt0].c=c;
    head0[v]=cnt0++;
}

void spfa()
{
    memset(dis,0x3f,sizeof(dis));
    memset(vis,0,sizeof(vis));
    dis[0]=0;

    queue<int> q;
    q.push(0);

    while(!q.empty())
    {
        int f=q.front();
        q.pop();
        vis[f]=0;
        for(int i=head[f];i!=-1;i=r[i].p)
        {
            if(dis[f]+r[i].c<dis[r[i].v])
            {
                dis[r[i].v]=dis[f]+r[i].c;
                if(!vis[r[i].v])
                {
                    vis[r[i].v]=1;
                    q.push(r[i].v);
                }
            }
        }
    }
}

LL dfs1(int u)
{
    if(vis[u])
        return dp1[u];
    vis[u]=true;
    for(int i=head0[u];i!=-1;i=r0[i].p)
    {
        if(dis[u]-r0[i].c==dis[r0[i].v])
        {
            dp1[u]+=dfs1(r0[i].v);
            dp1[u]%=MOD;
        }
    }
    return dp1[u];
}

LL dfs2(int u)
{
    if(dp2[u])
        return dp2[u];
    dp2[u]=1;
    for(int i=head[u];i!=-1;i=r[i].p)
    {
        if(dis[u]+r[i].c==dis[r[i].v])
        {
            dp2[u]+=dfs2(r[i].v);
            dp2[u]%=MOD;
        }
    }
    return dp2[u];
}

LL Mul(LL lhs,LL rhs)
{
    LL lhs2=lhs%100000;
    LL rhs2=rhs%100000;
    return ((lhs/100000*rhs2+rhs/100000*lhs2)*100000+lhs2*rhs2)%MOD;
}

int main()
{
    int n,m,q;
    int u,v,c;
    while(scanf("%d%d%d",&n,&m,&q)!=EOF)
    {
        cnt=0;
        cnt0=0;

        memset(head,-1,sizeof(head));
        memset(head0,-1,sizeof(head0));
        memset(dp1,0,sizeof(dp1));
        memset(dp2,0,sizeof(dp2));
        for(int i=0;i<m;i++)
        {
            scanf("%d%d%d",&u,&v,&c);
            addedge(u,v,c);
        }
        spfa();
        dfs2(0);
        dp1[0]=1;
        memset(vis,0,sizeof(vis));
        for(int i=0;i<n;i++)
                dfs1(i);
        for(int i=0;i<q;i++)
        {
            scanf("%d",&c);
            LL ans=Mul(dp1[c],dp2[c]);
            printf("%010lld\n",ans);
        }
    }
    return 0;
}