L2-001. 紧急救援

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

输入格式：

输入第一行给出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<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
#define INF 0x3f3f3f

using namespace std;

int n,m,s,d;
int a[1005],b[1005],path[1005],pathnum[1005],dist[1005];
/*a用来保存每个城市原有的救援队数量，b用来保存从起始城市到当前城市最多的救援队数量，path保存从最短路径到该城市的前一个城市，pathnum用来保存最短路径的条数*/
int Map[505][505],vis[505];
void Dijkstra(int v)
{
    int mindis,u;
    b[v]=a[v];
    for(int i=0;i<n;i++)//初始化dist，path，b和pathnum数组，dist初始为当前节点到所有节点的距离，等待更新，path若能直接到达则初始化为当前节点否则初始化为-1，若能直接到达则，b为a的值，pathnum为1
    {
        dist[i]=Map[v][i];
        path[i]=-1;
        if(dist[i]<INF)
        {
            path[i]=v;
            b[i]=a[v]+a[i];
            pathnum[i]=1;
        }
    }
    vis[v]=1;
    path[v]=-1;
    for(int i=0;i<n;i++)
    {
        mindis=INF;
        for(int j=0;j<n;j++)
            if(vis[j]==0&&mindis>dist[j])//找到dist的最小的数，并记录
            {
                mindis=dist[j];
                u=j;
            }
        vis[u]=1;
        for(int j=0;j<n;j++)
        {
            if(vis[j]==0&&dist[j]>mindis+Map[u][j])//更新dist数组同时更新b，pathnum和path数组
            {
                dist[j]=mindis+Map[u][j];
                b[j]=b[u]+a[j];
                pathnum[j]=pathnum[u];
                path[j]=u;
            }
            else if(vis[j]==0&&dist[j]==mindis+Map[u][j])//若相等则最短路径数目加一，同时更新b，pathnum和path数组
            {
                pathnum[j]=pathnum[u]+pathnum[j];
                if(b[j]<b[u]+a[j])
                {
                    b[j]=b[u]+a[j];
                    path[j]=u;
                }
            }
        }
    }
}
void Dispath(int D)
{
    if(path[D]!=-1)
    {
        Dispath(path[D]);
        cout<<path[D]<<" ";
    }
}
int main()
{
    int x,y,t;
    memset(Map,INF,sizeof(Map));
    memset(vis,0,sizeof(vis));
    cin>>n>>m>>s>>d;
    for(int i=0;i<n;i++)
    cin>>a[i];
    for(int i=0;i<m;i++)
    {
        cin>>x>>y>>t;
        Map[x][y]=Map[y][x]=t;
    }
    Dijkstra(s);
    cout<<pathnum[d]<<" "<<b[d]<<endl;
    int D=d;
    Dispath(d);
    cout<<d<<endl;
    return 0;
}