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
思路:想法是dij把最短路 在处理的过程中 把最短路的条数 路径 全部记录下来（由于小细节找了很久的bug）
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<string>
#include<vector>
#include<stack>
#include<bitset>
#include<cstdlib>
#include<cmath>
#include<set>
#include<list>
#include<deque>
#include<map>
#include<queue>
#define ll long long int
using namespace std;
inline ll gcd(ll a,ll b){return b?gcd(b,a%b):a;}
inline ll lcm(ll a,ll b){return a/gcd(a,b)*b;}
int moth[]={,,,,,,,,,,,,};
int dir[][]={, ,, ,-, ,,-};
int dirs[][]={, ,, ,-, ,,-, -,- ,-, ,,- ,,};
const int inf=0x3f3f3f3f;
const ll mod=1e9+;
int n,m,s,dd;
int G[][];
int num[];
int d[]; //最短路
int ma[]; //记录人数
int vis[];
int path[]; //路径数组
int ans[];
int ge[]; //最短路边数
int cnt=;
void output(int e){
if(e==-) return ;
output(path[e]);
ans[cnt++]=e;
}
void dij(int s){
memset(vis,,sizeof(vis));
memset(d,inf,sizeof(d));
d[s]=; ma[s]=num[s];
ge[s]=;
for(int i=;i<n;i++){
int k=-; int minn=inf;
for(int j=;j<n;j++){
if(d[j]<minn&&!vis[j]){
minn=d[j];
k=j;
}
}
if(k==-) return ;
vis[k]=;
for(int j=;j<n;j++){
if(j==k) continue;
if(d[j]>d[k]+G[k][j]){
d[j]=d[k]+G[k][j];
ma[j]=ma[k]+num[j];
path[j]=k;
ge[j]=ge[k];
}else if(d[j]==d[k]+G[k][j]){
ge[j]+=ge[k];
if(ma[j]<ma[k]+num[j]){
ma[j]=ma[k]+num[j];
path[j]=k;
}
}
}
}
}
int main(){
ios::sync_with_stdio(false);
cin>>n>>m>>s>>dd;
for(int i=;i<n;i++)
for(int j=;j<n;j++){
if(i==j)
G[i][j]=;
else
G[i][j]=inf;
}
memset(path,-,sizeof(path));
for(int i=;i<n;i++)
cin>>num[i];
for(int i=;i<=m;i++){
int a,b,c;
cin>>a>>b>>c;
G[a][b]=min(G[a][b],c);
G[b][a]=G[a][b];
}
dij(s);
cout<<ge[dd]<<" "<<ma[dd]<<endl;
output(dd);
for(int i=;i<cnt;i++){
if(i==)
cout<<ans[i];
else cout<<" "<<ans[i];
}
cout<<endl;
return ;
}