Description
While exploring his many farms, Farmer John has discovered a number of amazing wormholes. A wormhole is very peculiar because it is a one-way path that delivers you to its destination at a time that is BEFORE you entered the wormhole! Each of FJ's farms comprises N (1 ≤ N ≤ 500) fields conveniently numbered 1..N, M (1 ≤ M ≤ 2500) paths, and W (1 ≤ W≤ 200) wormholes.
As FJ is an avid time-traveling fan, he wants to do the following: start at some field, travel through some paths and wormholes, and return to the starting field a time before his initial departure. Perhaps he will be able to meet himself :) .
To help FJ find out whether this is possible or not, he will supply you with complete maps to F (1 ≤ F ≤ 5) of his farms. No paths will take longer than 10,000 seconds to travel and no wormhole can bring FJ back in time by more than 10,000 seconds.
Input
Line 1 of each farm: Three space-separated integers respectively: N, M, and W
Lines 2..M+1 of each farm: Three space-separated numbers (S, E, T) that describe, respectively: a bidirectional path between S and E that requires T seconds to traverse. Two fields might be connected by more than one path.
Lines M+2..M+W+1 of each farm: Three space-separated numbers (S, E, T) that describe, respectively: A one way path from S to E that also moves the traveler back T seconds.
Output
Sample Input
2 3 3 1 1 2 2 1 3 4 2 3 1 3 1 3 3 2 1 1 2 3 2 3 4 3 1 8
Sample Output
NO YES
Hint
For farm 2, FJ could travel back in time by the cycle 1->2->3->1, arriving back at his starting location 1 second before he leaves. He could start from anywhere on the cycle to accomplish this.
1 #include<iostream> 2 #include<cstdio> 3 #include<queue> 4 #include<string.h> 5 using namespace std; 6 const int maxn=505; 7 int T,x,y,z,n,m,k,u,v,w,cnt[maxn],dis[maxn];vector<int> v1[maxn],v2[maxn];bool vis[maxn];queue<int> que; 8 bool spfa(int s){ 9 while(!que.empty())que.pop(); 10 memset(vis,false,sizeof(vis)); 11 que.push(s),vis[s]=true,cnt[s]++,dis[s]=0;//注意:自己到本身的时间为dis[s]=0,同时累计每个顶点的入队次数,用于判负环 12 while(!que.empty()){ 13 u=que.front(),que.pop(),vis[u]=false; 14 for(size_t j=0;j<v1[u].size();++j){ 15 v=v1[u][j],w=v2[u][j]; 16 if(dis[u]+w<dis[v]){//松弛 17 dis[v]=dis[u]+w; 18 if(!vis[v]){ 19 que.push(v),vis[v]=true; 20 if(++cnt[v]>n)return true;//如果第n+1次仍然更新,则存在负圈 21 } 22 } 23 } 24 } 25 return false; 26 } 27 int main(){ 28 while(~scanf("%d",&T)){ 29 while(T--){ 30 scanf("%d%d%d",&n,&m,&k);memset(cnt,0,sizeof(cnt)); 31 for(int i=0;i<=n;++i)v1[i].clear(),v2[i].clear(); 32 for(int i=0;i<=n;++i)dis[i]=2e9; 33 while(m--){ 34 scanf("%d%d%d",&x,&y,&z);//建双向边 35 v1[x].push_back(y),v1[y].push_back(x); 36 v2[x].push_back(z),v2[y].push_back(z); 37 } 38 while(k--){ 39 scanf("%d%d%d",&x,&y,&z);//建单向边,权值为负数,表示时间倒流 40 v1[x].push_back(y); 41 v2[x].push_back(-z); 42 } 43 if(spfa(1))puts("YES"); 44 else puts("NO"); 45 } 46 } 47 return 0; 48 }