bzoj2200[Usaco2011 Jan]道路和航线
[Usaco2011 Jan]道路和航线
Time Limit: 10 Sec Memory Limit: 259 MBSubmit: 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 }
