PTA城市间紧急救援（Dijkstra）

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

输入格式:

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

比旅游规划那道题目复杂点，那道题只需要统计最短路径的长度和最小费用，而这道题需要记录路径，而且还需要记录最短路径长度相等的条数。
一开始我还看错题了，以为需要输出经过的点的个数…

使用path[i]表示i的前置路径，pathcnt[i]表示从s到i最短路径条数，res[i]表示从s到i的召集的救援队的数量。
参考：数据结构与算法题目集（中文） - 7-35 城市间紧急救援（25 分）

#include<iostream>
#include<cstring>
#include<vector>
#include<algorithm>
#include<stack>
#include<queue>
#include<map>
#include<list>
#include<set>
#include<cmath>
using namespace std;
typedef long long ll;
#define N 550
#define INF1 0x3f3f3f3f
#define INF2 0xc0c0c0c0
int n,m,s,d;

int dis[N],res[N],pathcnt[N],path[N],vis[N],edge[N][N],num[N];

void dijkstra(){
	
	for(int i=0;i<n;i++) dis[i] = edge[s][i];
	int e = s;path[s] = -1;dis[e] = 0;pathcnt[e]=1;
	while(e!=d){ 
		vis[e] = 1;
		for(int j=0;j<n;j++){
			if(!vis[j]){
				if(dis[j]>dis[e]+edge[e][j]){
					dis[j]=dis[e]+edge[e][j];
					res[j] = num[j] + res[e];
					pathcnt[j] = pathcnt[e];
					path[j] = e;
				}
				else if(dis[j]==dis[e]+edge[e][j]){
					pathcnt[j] += pathcnt[e];
					if(res[j]<num[j] + res[e]){
						res[j] = num[j] + res[e];
						path[j] = e;
					}
					
				}
			}	
		}
		int mindis = INF1;e = n;
		for(int j=0;j<n;j++){
			if(!vis[j]&&dis[j]<mindis){
				mindis = dis[j];
				e = j;
			}
		}
		if(e==n) break;
	}
	
}
int main() {
	cin >> n >> m >> s >> d;
	for(int i=0;i<n;i++) {
		scanf("%d",num+i); res[i] = num[i];
	}
	int x,y,z,h;
	for(int i=0;i<n;i++){
		for(int j=0;j<n;j++){
			edge[i][j]=edge[i][j]=INF1;
		}
		edge[i][i] = 0;
	}
	for(int i=0;i<m;i++){
		scanf("%d%d%d",&x,&y,&z);
		edge[x][y]=edge[y][x]= z;
	}
	dijkstra();
	vector<int> p;
	h = d;
	while(h!=-1){
		p.push_back(h);
		h = path[h];
	}
	cout << pathcnt[d] << " " << res[d] << endl; 
	for(int i=p.size()-1;i>=0;i--){
		if(i==p.size()-1) cout << p[i];
		else cout << " " << p[i];
	}
	return 0;
}