[luogu3008]/[USACO 11JAN]道路和飞机Roads and Planes

##链接##
Description

Farmer John正在一个新的销售区域对他的牛奶销售方案进行调查。他想把牛奶送到T个城镇 (1 <= T <= 25,000),编号为1T。这些城镇之间通过R条道路 (1 <= R <= 50,000,编号为1到R) 和P条航线 (1 <= P <= 50,000,编号为1到P) 连接。每条道路i或者航线i连接城镇A_i (1 <= A_i <= T)到B_i (1 <= B_i <= T),花费为C_i。对于道路,0 <= C_i <= 10,000;然而航线的花费很神奇,花费C_i可能是负数(-10,000 <= C_i <= 10,000)。道路是双向的,可以从A_i到B_i,也可以从B_i到A_i,花费都是C_i。然而航线与之不同,只可以从A_i到B_i。事实上,由于最近恐怖主义太嚣张,为了社会和谐,出台 了一些政策保证:如果有一条航线可以从A_i到B_i,那么保证不可能通过一些道路和航线从B_i回到A_i。由于FJ的奶牛世界公认十分给力,他需要运送奶牛到每一个城镇。他想找到从发送中心城镇S(1 <= S <= T) 把奶牛送到每个城镇的最便宜的方案,或者知道这是不可能的。
Input

第1行:四个空格隔开的整数: T, R, P, and S * 第2到R+1行:三个空格隔开的整数(表示一条道路):A_i, B_i 和 C_i * 第R+2到R+P+1行:三个空格隔开的整数(表示一条航线):A_i, B_i 和 C_i
Output

第1到T行:从S到达城镇i的最小花费,如果不存在输出”NO PATH”。

Sample Input

6 3 3 4

1 2 5

3 4 5

5 6 10

3 5 -100

4 6 -100

1 3 -10

样例输入解释:

一共六个城镇。在1-2,3-4,5-6之间有道路,花费分别是5,5,10。同时有三条航线:3->5,

4->6和1->3,花费分别是-100,-100,-10。FJ的中心城镇在城镇4。

Sample Output

NO PATH

NO PATH

5

0

-95

-100

样例输出解释:

FJ的奶牛从4号城镇开始,可以通过道路到达3号城镇。然后他们会通过航线达到5和6号城镇。

但是不可能到达1和2号城镇。

//坑爹的USACO(竟然卡SPFA (加了SLF也不行
//正解好麻烦,以后改

SLF优化(deque双端队列
在节点入队列时进行比较
如果比队首优先就加在队首
否则加在队尾

#include<iostream>
#include<cstdio>
#include<queue>
#include<cstring>
using namespace std;

const int maxn = 100000 + 100;
const int inf = 1<<30;
int n,r,t,s;
struct edge {
	int u,v,w;
	int next;
}e[maxn * 2];
int head[maxn],tot = 0;

int read() {
	int x = 0, f = 1;
	char ch = getchar();
	while(ch < '0' || ch > '9') {
		if(ch == '-') f = -1;
		ch = getchar();
	}
	while(ch >= '0' && ch <= '9') {
		x = x * 10 + ch - '0';
		ch = getchar();
	}
	return x * f;
}

void add(int u, int v, int w) {
	e[++tot] = (edge){u,v,w,head[u]};
	head[u] = tot;
}

int vis[maxn],d[maxn];
void spfa(int x) {
	deque<int>q;
	for(int i = 1; i <= n; i++) d[i] = inf;
	q.push_back(x);
	vis[x] = 1, d[x] = 0;
	while(!q.empty()) {
		int k = q.front();
		q.pop_front();
		vis[k] = 0;
		for(int i = head[k]; i; i = e[i].next) {
			int v = e[i].v;
			if(d[v] > d[k] + e[i].w) {
				d[v] = d[k] + e[i].w;
				if(!vis[v]) {
					if(!q.empty() && d[v] >= d[q.front()]) q.push_back(v);
					else q.push_front(v);
					vis[v] = 1;
				}
			}
		}
	}
}

int main() {
	n = read(), r = read(), t = read(), s = read();
	for(int i = 1; i <= r; i++) {
		int u = read(), v = read(), w = read();
		add(u,v,w);
		add(v,u,w);
	}
	for(int i = 1; i <= t; i++) {
		int u = read(), v = read(), w = read();
		add(u,v,w);
	}
	spfa(s);
	for(int i = 1; i <= n; i++) {
		if(d[i] == inf) puts("NO PATH");
		else printf("%d\n",d[i]);
	}
	return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值