作者 陈越
单位 浙江大学
作为一个城市的应急救援队伍的负责人,你有一张特殊的全国地图。在地图上显示有多个分散的城市和一些连接城市的快速道路。每个城市的救援队数量和每一条连接两个城市的快速道路长度都标在地图上。当其他城市有紧急求助电话给你的时候,你的任务是带领你的救援队尽快赶往事发地,同时,一路上召集尽可能多的救援队。
输入格式:
输入第一行给出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
思路:
照题目描述找从S到D的最短路,记录路径数。
方法:
用Dijkstra跑最短路,拿一个数组记录路径,一个数组记录道路数。
void dj(){ memset(vis,0,sizeof vis); memset(v,0x3f,sizeof v); priority_queue<P,vector<P>,greater<P> >q; //开始点s权值为0,救援人数最大值=此点救援人数,最优路径1条 v[s]=0,help[s]=hp[s],lj[s]=1; q.push({0,s});//入队 while(!q.empty()){ ll x=q.top().second; q.pop(); if(vis[x])continue; vis[x]=1; for(ll i = 0 ; i < n ; i ++){ if(vis[i] || mp[x][i] == 0x3f3f3f3f3f3f3f3f)continue; if(v[i] > v[x] + mp[x][i]){//看走这条路是不是更优 v[i] = v[x] + mp[x][i]; q.push({v[i],i}); //链接前一个点,继承最优路径数,更新救援人数最大值 lian[i]=x,lj[i]=lj[x],help[i]=help[x]+hp[i]; }else if(v[i] == v[x] + mp[x][i]){ lj[i]+=lj[x];//增加路径 if(help[i] < help[x]+hp[i])//更新救援人数最大值 help[i] = help[x]+hp[i],lian[i]=x; } } } return; }
AC代码:
#include<bits/stdc++.h> using namespace std; typedef long long ll; #define endl "\n" #define P pair<ll,ll> const ll N = 5e2+7; ll n,m,s,d; //邻接表,路径,权值,初始救援人数,每个点救援人数最大值,每个点有几条最优路径 ll mp[N][N],lian[N],v[N],hp[N],help[N],lj[N]; bool vis[N];//标记有没有来过 stack<ll>st; void dj(){ memset(vis,0,sizeof vis); memset(v,0x3f,sizeof v); priority_queue<P,vector<P>,greater<P> >q; //开始点s权值为0,救援人数最大值=此点救援人数,最优路径1条 v[s]=0,help[s]=hp[s],lj[s]=1; q.push({0,s});//入队 while(!q.empty()){ ll x=q.top().second; q.pop(); if(vis[x])continue; vis[x]=1; for(ll i = 0 ; i < n ; i ++){ if(vis[i] || mp[x][i] == 0x3f3f3f3f3f3f3f3f)continue; if(v[i] > v[x] + mp[x][i]){//看走这条路是不是更优 v[i] = v[x] + mp[x][i]; q.push({v[i],i}); //链接前一个点,继承最优路径数,更新救援人数最大值 lian[i]=x,lj[i]=lj[x],help[i]=help[x]+hp[i]; }else if(v[i] == v[x] + mp[x][i]){ lj[i]+=lj[x];//增加路径 if(help[i] < help[x]+hp[i])//更新救援人数最大值 help[i] = help[x]+hp[i],lian[i]=x; } } } return; } void solve(){ cin >> n >> m >> s >> d; memset(mp,0x3f,sizeof mp); memset(lian,-1,sizeof lian); for(ll i = 0 ; i < n ; i ++)cin >> hp[i]; while(m--){ ll x,y,z;cin >> x >> y >> z; mp[x][y]=mp[y][x]=min(mp[x][y],z); } dj(); cout << lj[d] << " " << help[d] << endl; ll x=d,f=0; while(lian[x] != -1)st.push(x),x=lian[x]; cout << s; while(!st.empty())cout << " " << st.top(),st.pop(); return; } int main(){ ll t=1;//cin >> t; while(t --)solve(); return 0; }