【PAT天梯赛练习集L2-1】紧急救援——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

SPFA的变形对于边 <v,u>,dis[]表示距离，num[]表示到达某一点的人的数量，Num[][2]:Num[][1]表示在更行之前的路径的数量， <script type="math/tex" id="MathJax-Element-106"> ,\quad dis[] 表示距离，num[]表示到达某一点的人的数量，Num[][2]:Num[][1] 表示在更行之前的路径的数量，</script>
$Num[][0]，表示在更新之后的路径的数量$
$如果dis[u]>dis[v]+e[i].w ,不言而喻，我们的首要条件是最短，所$
$以dis[u] = dis[v]+e[i].w,num[u] = a[u]+num[v],Num[u][1] = Num[u][0],Num[u][0] = Num[v][0]$
$如果dis[u] = dis[v]+e[i].w ,Num[u][1] = Num[u]$[0], Num[u][0]+= (Num[v][0]-Num[v][1])\quad num[u] ==max(num[u],num[v]+a[u])$

#include <iostream>
#include <cstring>
#include <cstdio>
#include <queue>
#include <stack>

using namespace std;

const int INF = 0x3f3f3f3f;

const int Max = 550;

typedef struct node
{
    int v,w,next;
}edge;

edge e[Max*Max];

int H[Max],top;


int Dis[Max];

int num[Max],Num[Max][2];

int pre[Max];

bool vis[Max];

int n,m,s,t;

int a[Max];

void AddEdge(int u,int v,int w)
{
    e[top].v = v; e[top].w = w;

    e[top].next =  H[u] ;  H[u] = top++;


    e[top].v = u; e[top].w = w;

    e[top].next =  H[v] ;  H[v] = top++;
}

void SPFA()
{
    memset(Dis,INF,sizeof(Dis));

    for(int i = 0;i<n;i++)
    {
        num[i] = a[i];

        pre[i] = -1;

        Num[i][0] = Num[i][1] = 0;

        vis[i] = false;
    }

    vis[s] = true;

    Dis[s] = 0;

    queue<int>Q;

    Num[s][0] = 1;

    Q.push(s);

    while(!Q.empty())
    {
        int u = Q.front();

        Q.pop();

        for(int i = H[u];~i;i= e[i].next)
        {
            int v = e[i].v;

            if(Dis[v]>Dis[u]+e[i].w)
            {
                Dis[v]  =Dis[u]+e[i].w;

                Num[v][1] = Num[v][0];
                Num[v][0] = Num[u][0];

                num[v] = a[v]+num[u];

                pre[v] = u;

                if(!vis[v])
                {
                    vis[v] = true;

                    Q.push(v);
                }
            }
            else if(Dis[v]==Dis[u]+e[i].w)
            {
                Num[v][1] = Num[v][0];
                Num[v][0] += (Num[u][0]-Num[u][1]);


                if(num[v]<num[u]+a[v])
                {
                    num[v] = num[u]+a[v];

                    pre[v] = u;
                }
                if(!vis[v])
                {
                    vis[v] = true;

                    Q.push(v);
                }

            }
        }

        vis[u]= false;
    }


}

int main()
{
    int u,v,w;

    while(~scanf("%d %d %d %d",&n,&m,&s,&t))
    {
        for(int i = 0;i<n;i++)
        {
            scanf("%d",&a[i]);
        }

        top =0 ;

        memset(H,-1,sizeof(H));

        for(int i = 0;i<m;i++)
        {
            scanf("%d %d %d",&u,&v,&w);

            AddEdge(u,v,w);
        }

        SPFA();

        printf("%d %d\n",Num[t][0],num[t]);

        stack<int>S;

        for(int i = t ;~i;i  = pre[i])
        {
            S.push(i);
        }
        bool flag = false;

        while(!S.empty())
        {
            if(!flag) flag= true;

            else printf(" ");

            printf("%d",S.top());

            S.pop();
        }
        printf("\n");
    }
    return 0;
}