L2-001 紧急救援

关于SPFA,他死了

没什么好说的

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


垃圾题不给清楚数据范围

Solution

很综合的最短路,最短路计数,最短路求解,双关键字最短路,打印路径
用Djs写就行

Code

#include<iostream>
#include<cstdio>
#include<queue>
#include<cmath>
#include<cstring>
#include<algorithm>
#include<climits>
#include<string>
#define LL long long
#define REP(i, x, y) for(LL i = (x);i <= (y);i++)
using namespace std;
int RD(){
    LL out = 0,flag = 1;char c = getchar();
    while(c < '0' || c >'9'){if(c == '-')flag = -1;c = getchar();}
    while(c >= '0' && c <= '9'){out = out * 10 + c - '0';c = getchar();}
    return flag * out;
    }
const int inf = 5000100, maxn = 510;
int num, nr, S, T;
int d[maxn], val[maxn], tot[maxn], from[maxn], cnt[maxn];
//臭逼:d:最短路   val:本地志愿者人数
//tot:到某点志愿者最大人数
//from:某点前驱 用于打印路径
//cnt 最短路计数
int head[maxn], nume;
struct Edge{
	int v, dis, nxt;
	}E[maxn * 1000];
void add(int u, int v, int dis){
	E[++nume] = (Edge){v, dis, head[u]};
	head[u] = nume;
	}
bool vis[maxn];
void Djs(int s){
	REP(i, 0, num - 1)d[i] = inf;
	d[s] = 0, tot[s] = val[s], cnt[s] = 1, from[s] = -1;
	priority_queue<pair<int, int> >Q;
	Q.push(make_pair(-d[s], s));
	while(!Q.empty()){
		int u = Q.top().second;Q.pop();
		if(vis[u])continue;
		vis[u] = 1;
		for(int i = head[u];i;i = E[i].nxt){
			int v = E[i].v, dis = E[i].dis;
			if(d[u] + dis < d[v]){
				d[v] = d[u] + dis;
				tot[v] = tot[u] + val[v];
				cnt[v] = cnt[u];
				from[v] = u;
				Q.push(make_pair(-d[v], v));
				}
			else if(d[u] + dis == d[v]){
				cnt[v] += cnt[u];
				if(tot[u] + val[v] > tot[v]){
					tot[v] = tot[u] + val[v];
					from[v] = u;
					Q.push(make_pair(-d[v], v));
					}
				}
			}
		}
	}
void init(){
	num = RD(), nr = RD(), S = RD(), T = RD();
	REP(i, 0, num - 1)val[i] = RD();
	while(nr--){
		int u = RD(), v = RD(), dis = RD();
		add(u, v, dis);
		add(v, u, dis);
		}
	}
void print(int x){
	if(from[x] == -1)return ;
	print(from[x]);
	printf("%d ", from[x]);
	}
void work(){
	Djs(S);
	printf("%d %d
", cnt[T], tot[T]);
	print(T);
	printf("%d
", T);
	}
int main(){
	init();
	work();
	return 0;
	}

原文地址:https://www.cnblogs.com/Tony-Double-Sky/p/14005812.html