L2-001. 紧急救援 SPFA+记录路径，统计最短路

作为一个城市的应急救援队伍的负责人，你有一张特殊的全国地图。在地图上显示有多个分散的城市和一些连接城市的快速道路。每个城市的救援队数量和每一条连接两个城市的快速道路长度都标在地图上。当其他城市有紧急求助电话给你的时候，你的任务是带领你的救援队尽快赶往事发地，同时，一路上召集尽可能多的救援队。

输入格式：

输入第一行给出4个正整数N、M、S、D，其中N（2<=N<=500）是城市的个数，顺便假设城市的编号为0~(N-1)；M是快速道路的条数；S是出发地的城市编号；D是目的地的城市编号。第二行给出N个正整数，其中第i个数是第i个城市的救援队的数目，数字间以空格分隔。随后的M行中，每行给出一条快速道路的信息，分别是：城市1、城市2、快速道路的长度，中间用空格分开，数字均为整数且不超过500。输入保证救援可行且最优解唯一。

输出格式：

第一行输出不同的最短路径的条数和能够召集的最多的救援队数量。第二行输出从S到D的路径中经过的城市编号。数字间以空格分隔，输出首尾不能有多余空格。

输入样例：

4 5 0 3
20 30 40 10
0 1 1
1 3 2
0 3 3
0 2 2
2 3 2

输出样例：

2 60
0 1 3

#include<cstdio>
#include<queue>
#include<iostream>
#include<cstring>
#include<algorithm>
using namespace std;
#define inf 0x3f3f3f3f
const int M =1e4+7;
int n,m,s,d,cnt;
struct edge
{
    int u,v,w;
    int ne;
}e[M*100];
int head[M*100],pre[M],dp[M][2];
int cost[M],dis[M],ans[M];
bool vis[M];
void init(int u,int v,int w)
{
    e[cnt].u=u;
    e[cnt].v=v;
    e[cnt].w=w;
    e[cnt].ne=head[u];
    head[u]=cnt++;
}
void solve()
{
    queue<int>q;
    for(int i=0;i<=n;i++)
    {
        dis[i]=inf;
        vis[i]=0;
    }
    dis[s]=0;
    vis[s]=1;
    q.push(s);
    dp[s][0]=1;
    while(q.size())
    {
        s=q.front();
        q.pop(); vis[s]=0;
        for(int i=head[s];i!=-1;i=e[i].ne)
        {
            int v=e[i].v;
            if(v==s)
                continue;
            if(dis[v]>dis[s]+e[i].w)
            {
                ans[v]=ans[s]+cost[v];
                dis[v]=dis[s]+e[i].w;
                if(!vis[v])
                {
                    q.push(v);
                    vis[v]=1;
                }
                pre[v]=s;
                dp[v][1]=dp[v][0];
                dp[v][0]=dp[s][0];
            }
            else
            if(dis[v]==dis[s]+e[i].w)
            {
               if(!vis[v])
                {
                    q.push(v);
                    vis[v]=1;
                }
                if(ans[v]<ans[s]+cost[v])
                {
                    ans[v]=ans[s]+cost[v];
                    pre[v]=s;
                }
                dp[v][1]=dp[v][0];
                dp[v][0]+=(dp[s][0]-dp[s][1]);
            }
        }
    }
}
void dfs(int x)
{
    if(pre[x]!=-1){
        dfs(pre[x]);
    printf("%d ",pre[x]);
    }
}
int main()
{
    cnt=0;
    memset(ans,0,sizeof(ans));
    memset(dp,0,sizeof(dp));
    memset(pre,-1,sizeof(pre));
    memset(head,-1,sizeof(head));
    cin>>n>>m>>s>>d;
    for(int i=0;i<n;i++){
        cin>>cost[i];
        ans[i]=cost[i];
    }
    while(m--)
    {
        int u,v,w;
        cin>>u>>v>>w;
        init(u,v,w);
        init(v,u,w);
    }
    solve();
    cout<<dp[d][0]<<' '<<ans[d]<<endl;
    dfs(d);
    printf("%d\n",d);
    return 0;
}