L2-001. 紧急救援（最短路dij+路径记录）（模板）

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

#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<cmath>
using namespace std;
const int maxx=505;
const int inf=0x3f3f3f3f;
int mapp[maxx][maxx],va[maxx],path[maxx],dis[maxx],re[maxx],totval[maxx],pathnum[maxx];
bool vis[maxx];
int n,m,s,d;
void disj(int v0)
{
for(int i=0;i<n;i++)//初始化
{
dis[i]=inf;
}
vis[v0]=true;//标记走过
dis[v0]=0;
totval[v0]=1;
pathnum[v0]=1;
for(int i=0;i<n;i++)
{
if(mapp[v0][i]!=inf&&i!=v0)
{
dis[i]=mapp[v0][i];
path[i]=v0;
totval[i]=va[v0]+va[i];//起始点到个个能到达城市能取得救援队的数目
pathnum[i]=1;
}
}
for(int i=0;i<n-1;i++)
{
int mi=inf,mival=0,u=v0;
for(int j=0;j<n;j++)
{
if(!vis[j]&&dis[j]<mi)//找当前情况下最短路
mi=dis[j],u=j;
}
vis[u]=1;//标记
for(int j=0;j<n;j++)
{
if(!vis[j])
{
if(dis[u]+mapp[u][j]<dis[j])
{
pathnum[j]=pathnum[u];
dis[j]=dis[u]+mapp[u][j];
totval[j]=totval[u]+va[j];//救援队伍数
path[j]=u;//保存路径
}
else if(dis[u]+mapp[u][j]==dis[j])
{
pathnum[j]+=pathnum[u];//记录条数
if(totval[j]<totval[u]+va[j])//最短路径有多个，召集的救援队伍不一定
{
totval[j]=totval[u]+va[j];
path[j]=u;
}
}
}
}
}
}
int main()
{
scanf("%d%d%d%d",&n,&m,&s,&d);//S出发地  D目的地
for(int i=0;i<n;i++)
{
scanf("%d",&va[i]);
}
for(int i=0;i<n;i++)
for(int j=0;j<n;j++)
mapp[i][j]=inf;
int x,y,l;
for(int i=0;i<m;i++)
{
scanf("%d%d%d",&x,&y,&l);
mapp[x][y]=min(mapp[x][y],l);
mapp[y][x]=mapp[x][y];
}
disj(s);
int num=0,cur=d;
while(cur!=s)//回溯 确定路径
{
re[num++]=cur;
cur=path[cur];
}
re[num++]=s;
printf("%d %d\n",pathnum[d],totval[d]);
for(int i=num-1;i>0;i--)
printf("%d ",re[i]);
printf("%d\n",re[0]);
return 0;
}