作为一个城市的应急救援队伍的负责人,你有一张特殊的全国地图。在地图上显示有多个分散的城市和一些连接城市的快速道路。每个城市的救援队数量和每一条连接两个城市的快速道路长度都标在地图上。当其他城市有紧急求助电话给你的时候,你的任务是带领你的救援队尽快赶往事发地,同时,一路上召集尽可能多的救援队。
输入格式:
输入第一行给出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;
但此题还加了一些条件 要算出最短路径(关键就是0到j的距离>0到k+k到j 的距离就更新)且救援队伍数量最多和输出路径
if(dis[j]>dis[k]+mapp[k][j]){//如果距离0到j的距离>0到k+k到j的距离就更新
dis[j]=dis[k]+mapp[k][j];
maxx[j]=maxx[k]+a[j];
num[j]=num[k];
pre[j]=k;//更新j点可以到k点
}
else if(dis[j]==dis[k]+mapp[k][j]){//如果距离一样
num[j]=num[j]+num[k];
if(maxx[j]<maxx[k]+a[j]){//且救援队伍少就更新
maxx[j]=maxx[k]+a[j];
pre[j]=k;//更新路径
}
}
#include<bits/stdc++.h>
using namespace std;
int n,m,s,d;
int a[505];//输入的救援队伍
int num[505];//最短路径的条数
int dis[505];//起始点s到每个点的最短路径
int maxx[505];//每个点能召集最多的救援队伍
int vis[505];//判断此点是否可以通过0是可以通过1是不能通过
int pre[505];//输出路径上的城市
int mod=10e8+7;
int mapp[505][505];
void dijkstra(int s){
for(int i=0;i<n;i++){//先把0到每个点的路径和救援队伍更新
if(i==s){
dis[i]=0;
maxx[i]=a[i];
}else{
dis[i]=mapp[s][i];
maxx[i]=a[s]+a[i];
}
if(i!=s&&mapp[s][i]!=mod){
pre[i]=s;
}else{
pre[i]=-1;
}
}
vis[s]=1;//0点(s点)已经走过不能再走
for(int i=0;i<n;i++){
int t=mod,k;
for(int j=0;j<n;j++){
if(!vis[j]&&t>dis[j]){
t=dis[j];
k=j;
}
}
if(t==mod){
break;
}
vis[k]=1;
for(int j=0;j<n;j++){
if(!vis[j]){
if(dis[j]>dis[k]+mapp[k][j]){//如果距离0到j的距离>0到k+k到j的距离就更新
dis[j]=dis[k]+mapp[k][j];
maxx[j]=maxx[k]+a[j];
num[j]=num[k];
pre[j]=k;//更新j点可以到k点
}
else if(dis[j]==dis[k]+mapp[k][j]){
num[j]=num[j]+num[k];
if(maxx[j]<maxx[k]+a[j]){
maxx[j]=maxx[k]+a[j];
pre[j]=k;
}
}
}
}
}
return ;
}
void p(){//输出路径上的城市
stack<int> s;
for(int i=d;i!=-1;i=pre[i]) s.push(i);
int l=s.size();
l--;
while(l--){
cout<<s.top()<<" ";
s.pop();
}
cout<<s.top();
}
int main(){
cin>>n>>m>>s>>d;
for(int i=0;i<n;i++){
cin>>a[i];
num[i]=1;
}
memset(mapp,mod,sizeof(mapp));
for(int i=1;i<=m;i++){
int x,y,z;
cin>>x>>y>>z;
mapp[x][y]=min(mapp[x][y],z);
mapp[y][x]=min(mapp[y][x],z);
}
dijkstra(s);
printf("%d %d\n",num[d],maxx[d]);
p();
return 0;
}