PAT L2-001 紧急救援

L2-001. 紧急救援

作为一个城市的应急救援队伍的负责人，你有一张特殊的全国地图。在地图上显示有多个分散的城市和一些连接城市的快速道路。每个城市的救援队数量和每一条连接两个城市的快速道路长度都标在地图上。当其他城市有紧急求助电话给你的时候，你的任务是带领你的救援队尽快赶往事发地，同时，一路上召集尽可能多的救援队。
输入格式：
输入第一行给出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算法找最短路，但是在最短时间赶去的基础上又要带尽量多的救援队。所以在最短路不唯一的情况下判断哪条路上可以召集尽可能多的救援队。

代码如下：
    
#include<stdio.h>
#include<string.h>
#define INF 0x3f3f3f3f
#define N 520

int mp

, dis
, maxd
;//mp数组存储图的信息 dis数组求最短路 maxd数组存最大救援队
int path
, pathnum
, vis
, book
;//path数组标记当前点前驱结点 pathnum数组记录路径个数
int n, m, st, en; //vis存储各地救援队数量 book标记走过的结点

void print( int x )
{
if( path[x] == -1 )
return ; //递归输出路径
print( path[x] );
printf("%d ",path[x]);
}

void Dijkstra( int v0 )
{
for( int i = 0; i < n; i++ )
{
dis[i] = mp[v0][i];//初始化起点到各点距离
book[i] = 0; //初始化各点都没走过
path[i] = -1; //各点前驱结点为-1
}
pathnum[v0] = 1;
book[v0] = 1;
maxd[v0] = vis[v0];
for( int i = 0; i < n; i++ )
{
if( mp[v0][i] != INF && i != v0 )
{
path[i] = v0;
maxd[i] = vis[i] + vis[v0];
pathnum[i] = 1;
}
}
for( int v = 1; v < n; v++ )
{
int u = v0, minn = INF;
for( int i = 0; i < n; i++ )
{
if( !book[i] && dis[i] < minn )
{
minn = dis[i];
u = i;
}
}
book[u] = 1;
for( int i = 0; i < n; i++ )
{
if( !book[i] )
{
if( dis[i] > dis[u] + mp[u][i] )
{
dis[i] = dis[u] + mp[u][i];
maxd[i] = maxd[u] + vis[i];
pathnum[i] = pathnum[u];
path[i] = u;
}
else if( dis[i] == dis[u] + mp[u][i] )
{
pathnum[i] += pathnum[u];
if( maxd[i] < maxd[u] + vis[i] )
{
maxd[i] = maxd[u] + vis[i];
path[i] = u;
}

}
}
}
}
}

int main()
{
while( ~scanf("%d%d%d%d",&n,&m,&st,&en) )
{
for( int i = 0; i < n; i++ )
scanf("%d",&vis[i]);
for( int i = 0; i < n; i++ )
{
for( int j = 0; j < n; j++ )
{ //初始化图中各点间的权值
if( i == j )
mp[i][j] = 0;
else
mp[i][j] = INF;
}
}
int t1, t2, t3;
for( int i = 0; i < m; i++ )
{
scanf("%d%d%d",&t1,&t2,&t3);
if( mp[t1][t2] > t3 ) //更新最小的权值
mp[t1][t2] = mp[t2][t1] = t3;
}
Dijkstra(st);
printf("%d %d\n",pathnum[en], maxd[en]);
print( en );
printf("%d\n",en);
}
return 0;
}