LOJ10081
USACO 201 Jan. Gold
Farmer John 正在一个新的销售区域对他的牛奶销售方案进行调查。他想把牛奶送到 T 个城镇 ,编号为 1 到 T。这些城镇之间通过 R 条道路(编号为 1 到 R)和 P 条航线(编号为 1 到 P)连接。每条道路 i 或者航线 i 连接城镇 Ai 到 Bi,花费为 Ci。
对于道路,0 \le C_i \le 10^40≤Ci≤104,然而航线的花费很神奇,花费 Ci 可能是负数。道路是双向的,可以从 Ai 到 Bi,也可以从 Bi 到 Ai,花费都是 Ci。然而航线与之不同,只可以从 Ai 到 Bi。
事实上,由于最近恐怖主义太嚣张,为了社会和谐,出台了一些政策保证:如果有一条航线可以从 Ai 到 Bi,那么保证不可能通过一些道路和航线从 Bi 回到 Ai。由于 FJ 的奶牛世界公认十分给力,他需要运送奶牛到每一个城镇。他想找到从发送中心城镇 SS 把奶牛送到每个城镇的最便宜的方案,或者知道这是不可能的。
输入格式
第一行为四个空格隔开的整数:T,R,P,S;
第二到第 R+1 行:三个空格隔开的整数(表示一条道路):Ai,Bi 和 Ci;
第 R+2 到 R+P+1 行:三个空格隔开的整数(表示一条航线):Ai,Bi 和 Ci。
输出格式
输出 T 行,第 i 行表示到达城镇 i 的最小花费,如果不存在输出 NO PATH。
样例
样例输入
6 3 3 4
1 2 5
3 4 5
5 6 10
3 5 -100
4 6 -100
1 3 -10
样例输出
NO PATH
NO PATH
5
0
-95
-100
样例说明
一共六个城镇。在 1 和 2,3 和 4,5 和 6 之间有道路,花费分别是 5,5,10。同时有三条航线:3→5,4→6 和 1→3,花费分别是 −100,−100,−10。FJ 的中心城镇在城镇 4。FJ 的奶牛从 4 号城镇开始,可以通过道路到达 3 号城镇。然后他们会通过航线达到 5 和 6 号城镇。但是不可能到达 1 和 2 号城镇。
数据范围与提示
对于全部数据,1≤T≤2.5×104,1≤R,P≤5×104,1≤Ai,Bi,S≤T。保证对于所有道路,0≤Ci≤104,对于所有航线,−104≤Ci≤104。
_______________________________________________________
有负边,明显SPFA。
超时,注意SPFA的优化。
_______________________________________________________
![](https://images.cnblogs.com/OutliningIndicators/ContractedBlock.gif)
1 #include<bits/stdc++.h> 2 using namespace std; 3 const int maxn=25010; 4 const int maxm=150010; 5 int n,r,p,s; 6 struct edge 7 { 8 int u,v,w,nxt; 9 }e[maxm]; 10 int head[maxn],js; 11 void addage(int u,int v,int w) 12 { 13 e[++js].u=u;e[js].v=v;e[js].w=w; 14 e[js].nxt=head[u];head[u]=js; 15 } 16 int dis[maxn]; 17 bool inq[maxn]; 18 void spfa(int x) 19 { 20 memset(dis,0x3f,sizeof dis); 21 inq[x]=1;dis[x]=0; 22 deque<int>q; 23 q.push_front(x); 24 while(!q.empty()) 25 { 26 int u=q.front(); 27 q.pop_front(); 28 inq[u]=0; 29 for(int i=head[u];i;i=e[i].nxt) 30 { 31 int v=e[i].v; 32 if(dis[v]>dis[u]+e[i].w) 33 { 34 dis[v]=dis[u]+e[i].w; 35 if(!inq[v]) 36 { 37 if(q.empty() || dis[v]<dis[q.front()])q.push_front(v); 38 else q.push_back(v); 39 inq[v]=1; 40 } 41 } 42 } 43 } 44 } 45 inline int read(){ 46 int ret=0;char ch=getchar();bool flg=0; 47 while(!isdigit(ch)) flg^=!(ch^'-'),ch=getchar(); 48 while(isdigit(ch)) ret=(ret<<3)+(ret<<1)+ch-48,ch=getchar(); 49 return flg?-ret:ret; 50 } 51 int main() 52 { 53 n=read();r=read();p=read();s=read(); 54 for(int u,v,w,i=0;i<r;++i) 55 { 56 u=read();v=read();w=read(); 57 addage(u,v,w);addage(v,u,w); 58 } 59 for(int u,v,w,i=0;i<p;++i) 60 { 61 u=read();v=read();w=read(); 62 addage(u,v,w); 63 } 64 spfa(s); 65 for(int i=1;i<=n;++i) 66 if(dis[i]==0x3f3f3f3f)puts("NO PATH"); 67 else printf("%d\n",dis[i]); 68 return 0; 69 }