【PTA】7-35 城市间紧急救援（Dijkstra最短路）

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

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

思路分析：

在原Dikjstra算法里加上对救援人数的更新即可。这道题一开始只拿了18分，原因是没考虑到最短路径条数需要继承上一个顶点的条数。假设 0->1 有3条最短路径，则cnt[1] = 3. 1->3有一条路径，但是0 ->1 -> 3 也有3条路径，所以cnt[3] = cnt[1] + cnt[3]; 这个坑点比较细微，没考虑到…
以下是AC代码

#include <iostream>
#include <cstring>
#define MAX 505
using namespace std;

int map[MAX][MAX];
int dis[MAX];	//S到i最短距离 
int vis[MAX];	
int psn[MAX];	//i城的人数 
int val[MAX];	//S到i城的人数 
int cnt[MAX];	//最短路径条数 
int path[MAX];	//path[i] = x 表示i的父结点是x 
int N, M, S, D;

void Dijkstra(int S)
{
	memset(dis,0x3f,sizeof(dis));
	dis[S] = 0;
	val[S] = psn[S];
	path[S] = -1;
	while(1)
	{
		int k = -1, Min = 0x3f3f3f3f;
		for(int i = 0;i < N;i++)	//找未遍历的最小dis下标 
		{
			if(!vis[i] && dis[i] < Min)
			{
				k = i;
				Min = dis[i];
			}
		}
		if(k == -1) return;
		vis[k] = 1;				//收录顶点 
		for(int i = 0;i < N;i++)	//更新dis和val 
		{
			if(!vis[i] && dis[i] > dis[k] + map[k][i])	//找到更短路径 
			{
				dis[i] = dis[k] + map[k][i];			//更新dis 
				val[i] = val[k] + psn[i];				//更新val 
				path[i] = k;
				cnt[i] = cnt[k];	//出度的顶点cnt要继承上一个顶点的cnt 
			}
			else if(!vis[i] && dis[i] == dis[k] + map[k][i])	//有相同长度的路径 
			{
				cnt[i] += cnt[k]; 	//相同路径数+1 
				if(val[k] + psn[i] > val[i])	 //更新到最大救援数 
				{
					path[i] = k;
					val[i] = val[k] + psn[i];
				}
			}
		}
	}
}

void print(int x)
{
	if(path[x] == -1) return;
	print(path[x]);	
	cout<<" "<<x;
}

int main()
{
	int x, y, z;
	memset(map,0x3f,sizeof(map));
	cin>>N>>M>>S>>D;
	for(int i = 0;i < N;i++) cin>>psn[i], cnt[i] = 1;	//初始化城市救援人数，路径条数 
	while(M--)
	{
		cin>>x>>y>>z;
		map[x][y] = z;
		map[y][x] = z;
	}
	Dijkstra(S);
	cout<<cnt[D]<<" "<<val[D]<<endl;
	cout<<S;
	print(D);
	return 0;
}