bzoj2200[Usaco2011 Jan]道路和航线

 [Usaco2011 Jan]道路和航线

Time Limit: 10 Sec  Memory Limit: 259 MB
Submit: 742  Solved: 265
[Submit][Status][Discuss]

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号城镇。
 
题解:
  貌似是Spfa+slf优化,就是更短放队列首那个,所以用deque。
 1 #include<cstring>
 2 #include<cmath>
 3 #include<algorithm>
 4 #include<iostream>
 5 #include<cstdio>
 6 #include<deque>
 7 
 8 #define N 25007
 9 #define M 100007
10 #define inf 2000000007
11 using namespace std;
12 inline int read()
13 {
14     int x=0,f=1;char ch=getchar();
15     while(ch>'9'||ch<'0'){if (ch=='-') f=-1;ch=getchar();}
16     while(ch<='9'&&ch>='0'){x=(x<<3)+(x<<1)+ch-'0';ch=getchar();}
17     return x*f;
18 }
19 
20 int n,m1,m2,S;
21 int cnt,hed[N],rea[M*2],nxt[M*2],val[M*2];
22 int dis[N];
23 bool flag[N];
24 deque<int>q;
25 
26 void add(int u,int v,int z)
27 {
28     nxt[++cnt]=hed[u];
29     hed[u]=cnt;
30     rea[cnt]=v;
31     val[cnt]=z;
32 }
33 void Spfa()
34 {
35     for (int i=1;i<=n;i++)
36         dis[i]=inf;
37     dis[S]=0,flag[S]=1;q.push_back(S);
38     
39     while(!q.empty())
40     {
41         int u=q.front();q.pop_front();
42         for (int i=hed[u];i!=-1;i=nxt[i])
43         {
44             int v=rea[i],fee=val[i];
45             if (dis[v]>dis[u]+fee)
46             {
47                 dis[v]=dis[u]+fee;
48                 if (!flag[v])
49                 {
50                     if (!q.empty()&&dis[q.front()]>dis[v]) q.push_front(v);
51                     else q.push_back(v);
52                     flag[v]=1;
53                 }
54             }
55         }
56         flag[u]=0;
57     }
58 }
59 int main()
60 {
61     memset(hed,-1,sizeof(hed));
62     n=read(),m1=read(),m2=read(),S=read();
63     for (int i=1;i<=m1;i++)
64     {
65         int x=read(),y=read(),z=read();
66         add(x,y,z),add(y,x,z);
67     }
68     for (int i=1;i<=m2;i++)
69     {
70         int x=read(),y=read(),z=read();
71         add(x,y,z);
72     }
73     Spfa();
74     for (int i=1;i<=n;i++)
75         if (dis[i]==inf) printf("NO PATH\n");
76         else printf("%d\n",dis[i]);
77 }

 

posted @ 2017-12-29 09:14  Kaiser-  阅读(556)  评论(0编辑  收藏  举报