L2-001. 紧急救援（PAT 最短路+记录路径）

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

输入格式：

输入第一行给出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

思路：最短路的变形题,a数组记录的是每个点上消防队员的多少，b数组记录的是这一条路径上总共的消防队员数，path记录的每条最短路径经过的节点，pathnum记录的是最短路径有多少条，之所以用数组来存的原因是当走完完整的一条之后最短路才增加一，因此拿数组来记录状态。

#include<stdio.h>
#include<stdlib.h>
#include<string.h>
#define INF 0x3f3f3f3f
int Map[550][550],path[1010],a[1010],b[1010],pathnum[1010],dis[1010],vis[550];
int n,m,d,s;
void dijstra(int v)
{
    int i,j;
    b[v] = a[v];
    for(i=0; i<n; i++)
    {
        dis[i] = Map[v][i];
        path[i] = -1;
        if(dis[i]<INF)
        {
            b[i] = a[v]+a[i];
            pathnum[i] = 1;
            path[i] = v;
        }
    }
    vis[v] = 1;
    path[v] = -1;
    for(i=0; i<n; i++)
    {
        int min = INF,pos;
        for(j=0; j<n; j++)
        {
            if(!vis[j]&&dis[j]<min)
            {
                min = dis[j];
                pos = j;
            }
        }
        vis[pos]  =1;
        for(j=0; j<n; j++)
        {
            if(!vis[j]&&dis[j]>min+Map[pos][j])
            {
                dis[j] = min+Map[pos][j];
                b[j] = b[pos] + a[j];
                path[j] = pos;
                pathnum[j] = pathnum[pos];
            }
            else if(!vis[j]&&dis[j] == min+Map[pos][j])
            {
                pathnum[j] += pathnum[pos];
                if(b[j] < b[pos] + a[j])
                {
                    b[j] = b[pos] + a[j];
                    path[j] = pos;
                }
            }
        }

    }
}
void Print(int d)
{
    if(path[d]!=-1)
    {
        Print(path[d]);
        printf("%d ",path[d]);
    }
}
int main()
{

    memset(Map,INF,sizeof(Map));
    memset(vis,0,sizeof(vis));
    scanf("%d %d %d %d",&n,&m,&s,&d);
    for(int i=0; i<n; i++)
        scanf("%d",&a[i]);
    int x,y,t;
    for(int i=0; i<m; i++)
    {
        scanf("%d %d %d",&x,&y,&t);
        Map[x][y] = Map[y][x] = t;
    }
    dijstra(s);
    printf("%d %d\n",pathnum[d],b[d]);
    Print(d);
    printf("%d\n",d);
    return 0;
}