PTA 5-7 城市间紧急救援 【Dijkstra+记录路径】

5-7 城市间紧急救援   (25分)

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

输入格式:

输入第一行给出4个正整数$N$、$M$、$S$、$D$，其中$N$（$2\le N\le 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<cstdio>
#include<math.h>
#include<cstring>
#include<climits>
#include<string>
#include<queue>
#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<climits>
#include<string>
#include<queue>
#include<stack>
#include<set>
#include<map>
#include<list>
#include<sstream>
#include<algorithm>
using namespace std;
#define rep(i,j,k)for(i=j;i<k;i++)
#define per(i,j,k)for(i=j;i>k;i--)
#define MS(x,y)memset(x,y,sizeof(x))
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
#define ll long long
#define abs(x) (x>0?x:-x)
const int INF=0x7ffffff;

const int M=1e4+10;
int i,j,k,n,m;
int d,s;
int mp[M][M];
int vis[M],dis[M],val[M],cnt[M],pre[M],sum[M];

void init()
{
    for(i=0;i<=n;i++){
     pre[i]=-1;
     dis[i]=INF;
     for(j=0;j<=n;j++){
        mp[i][j]=INF;
        if(i==j)mp[i][j]=0;
     }
    }
    MS(vis,0);MS(cnt,0);MS(sum,0);
}

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

void Dijkstra()
{
    vis[s]=1;sum[s]=val[s];
    dis[s]=0;cnt[s]=1;
    for(i=0;i<n;i++){
        int mink=s,minx=INF;
        for(j=0;j<n;j++){
            if(!vis[j]&&dis[j]<minx)
                mink=j,minx=dis[j];
        }
        vis[mink]=1;
        for(j=0;j<n;j++){
            if(!vis[j]&&dis[j]>dis[mink]+mp[mink][j]){
                dis[j]=dis[mink]+mp[mink][j];
                pre[j]=mink;
                sum[j]=sum[mink]+val[j];
                cnt[j]=cnt[mink];
            }
            else if(!vis[j]&&dis[j]==dis[mink]+mp[mink][j]){
                cnt[j]=cnt[mink]+cnt[j];
                if(sum[j]<sum[mink]+val[j]){
                    sum[j]=sum[mink]+val[j];
                    pre[j]=mink;
                }
            }
        }
    }
    printf("%d %d\n",cnt[d],sum[d]);
    path(d);
    printf("%d\n",d);
}

int main()
{
    while(~scanf("%d%d%d%d",&n,&m,&s,&d)){
        for(i=0;i<n;i++)scanf("%d",&val[i]);
        init();
        for(i=0;i<m;i++){
           int a,b,c;
           scanf("%d%d%d",&a,&b,&c);
            if(c<mp[a][b])mp[a][b]=c,mp[b][a]=c;
        }
        Dijkstra();
    }
}