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

来源:互联网 发布:电脑微信群发软件 编辑:程序博客网 时间:2024/04/26 11:34

链接

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;}
原创粉丝点击