城市间紧急救援 ,复杂点的最短路，记录路径，记录最短路个数等

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

 

• 时间限制：400ms
• 内存限制：64MB
• 代码长度限制：16kB
• 判题程序：系统默认
• 作者：陈越
• 单位：浙江大学

#include<bits/stdc++.h>
using namespace std;
#define LL long long
#define put() puts("***********")
#define INF 12111111
const int N = 3e3+10;
int a[N][N];
int sum[N];
int dis[N],pre[N], vis[N],cnt[N], val[N];
void path(int d) {
	if(pre[d] != -1) {

		path(pre[d]);printf("%d ", pre[d]);

	}
}
int main(){
    int n,m,s,d;
    cin>>n>>m>>s>>d;
    for(int i=0;i<n;i++){
        scanf("%d",&val[i]);
    }
    for(int i=0;i<n;i++){
        for(int j=0;j<n;j++){
            if(i==j)a[i][j]=0;
            else a[i][j]=INF;
        }
    }
    for(int i=0;i<m;i++){
            int u,v,w;
        scanf("%d%d%d",&u,&v,&w);
        a[u][v]=a[v][u]=w;
    }
    for(int i=0;i<n;i++){
        dis[i]=INF;
    }


    dis[s] = 0;
	vis[s] = 1;
	cnt[s] = 1;
	sum[s] = val[s];
    memset(pre,-1,sizeof(pre));
    for(int i=0;i<n;i++){
        int mmin=INF;
        int flag=s;
        for(int j=0;j<n;j++){
            if(!vis[j]&&dis[j]<mmin){
                mmin=dis[j];
                flag=j;
            }
        }

        vis[flag]=1;
        for(int j=0;j<n;j++){
            if(!vis[j]){
                if(a[j][flag]+dis[flag]<dis[j]){
                   dis[j] = dis[flag] + a[flag][j];
                   sum[j] = sum[flag] + val[j];
                   cnt[j] = cnt[flag];
                   pre[j] = flag;
                }
                else if(dis[j] == dis[flag] + a[flag][j]) {
					cnt[j] = cnt[flag] + cnt[j];
					if(sum[j] < sum[flag] + val[j]) {
						sum[j] = sum[flag] + val[j];
						pre[j] = flag;
					}
				}
            }
        }
    }

    printf("%d %d\n", cnt[d], sum[d]);
    path(d);
    printf("%d\n",d);
    return 0;
}