天梯赛L2-001 紧急救援 (25 分)

L2-001 紧急救援 (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

这一题要求我们求最短路径，同时符合之前学习的Dijkstra算法附加情况中的两种，点权 + 最短路径条数，而且没有负权边，那么我们可以在这一题中使用这个算法来解题，端点数小于1000，用邻接矩阵来存储图会比较方便。最后输出这条路径的时候，需要用一个数组存储使用DFS的思想来进行遍历和输出，实在理解不了可以记忆。。。

#include <iostream>
#include <algorithm>
using namespace std;

const int MAXV = 550;
const int INF = 1000000000;
int G[MAXV][MAXV];
int d[MAXV], w[MAXV] = { 0 }, num[MAXV] = {0};
int weight[MAXV];
bool vis[MAXV] = { false };
//存储目标路径的端点的数组pre[]
int pre[MAXV];
int n, m, st, ed;
void Dijkstra(int s)
{
	fill(d, d + MAXV, INF);
	d[s] = 0;
	num[s] = 1;
	w[s] = weight[s];
	//将所有端点的前驱都设为本身
	for (int i = 0; i < n; i++)
		pre[i] = i;
	for (int i = 0; i < n; i++)
	{
		int u = -1;
		int min = INF;
		for (int j = 0; j < n; j++)
		{
			if (vis[j] == false && d[j] < min)
			{
				u = j;
				min = d[j];
			}
		}
		if (u == -1)
			return;
		vis[u] = true;
		for (int v = 0; v < n; v++)
		{
			if (vis[v] == false && G[u][v] != INF)
			{
				if (d[u] + G[u][v] < d[v])
				{
					d[v] = d[u] + G[u][v];
					w[v] = w[u] + weight[v];
					num[v] = num[u];
					//记录v的前驱是u
					pre[v] = u;
				}
				else if (d[u] + G[u][v] == d[v])
				{
					if (w[u] + weight[v] > w[v])
					{
						w[v] = w[u] + weight[v];
						pre[v] = u;
					}
					num[v] += num[u];
				}
			}
		}
	}
}
void DFS(int s, int v)
{
	if (v == s)
	{
		cout << v;
		return;
	}
	DFS(s, pre[v]);
	cout << " " << v;
}
int main()
{
	cin >> n >> m >> st >> ed;
	for (int i = 0; i < n; i++)
		cin >> weight[i];
	int u, v;
	fill(G[0], G[0] + MAXV * MAXV, INF);
	for (int i = 0; i < m; i++)
	{
		cin >> u >> v;
		cin >> G[u][v];
		G[v][u] = G[u][v];
	}
	Dijkstra(st);
	cout << num[ed] << " " << w[ed] << endl;
	DFS(st, ed);
	return 0;
}