题干：

农夫约翰在探索他的许多农场，发现了一些惊人的虫洞。虫洞是很奇特的，因为它是一个单向通道，可让你进入虫洞的前达到目的地！他的N（1≤N≤500）个农场被编号为1..N，之间有M（1≤M≤2500）条路径，W（1≤W≤200）个虫洞。FJ作为一个狂热的时间旅行的爱好者，他要做到以下几点：开始在一个区域，通过一些路径和虫洞旅行，他要回到最开时出发的那个区域出发前的时间。也许他就能遇到自己了:)。为了帮助FJ找出这是否是可以或不可以，他会为你提供F个农场的完整的映射到（1≤F≤5）。所有的路径所花时间都不大于10000秒，所有的虫洞都不大于万秒的时间回溯。

Input

第1行：一个整数F表示接下来会有F个农场说明。 每个农场第一行：分别是三个空格隔开的整数：N，M和W 第2行到M+1行：三个空格分开的数字（S，E，T）描述，分别为：需要T秒走过S和E之间的双向路径。两个区域可能由一个以上的路径来连接。 第M +2到M+ W+1行：三个空格分开的数字（S，E，T）描述虫洞，描述单向路径，S到E且回溯T秒。

Output

F行，每行代表一个农场 每个农场单独的一行，” YES”表示能满足要求，”NO”表示不能满足要求。

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 1, FJ cannot travel back in time. 
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.

 

AC代码：（spfa，邻接矩阵存图）（1079ms）

#include<cstdio>
#include<algorithm>
#include<iostream>
#include<queue>
#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;
int n;
const int MAX = 505;
const int INF = 0x3f3f3f3f;
int dis[MAX],maze[MAX][MAX],cnt[MAX];
bool vis[MAX];
bool spfa(int s){
//	for(int i=0; i<=n; i++) dis[i]=99999999; //初始化每点i到s的距离dis[s]=0; vis[s]=1;  //队列初始化,s为起点int i, v;queue<int > q;q.push(s);while (!q.empty()){   //队列非空v=q.front();  //取队首元素q.pop();vis[v]=0;   //释放队首结点，因为这节点可能下次用来松弛其它节点，重新入队for(i=1; i<=n; i++)  //对所有顶点if (maze[v][i]!=INF && dis[i]>dis[v]+maze[v][i]){  dis[i] = dis[v]+maze[v][i];   //修改最短路if (vis[i]==0){  //如果扩展结点i不在队列中，入队cnt[i]++;vis[i]=1;q.push(i);if(cnt[i] >=n) return true;}}}return false;
}void init() {memset(vis,0,sizeof(vis));memset(maze,INF,sizeof(maze));memset(dis,INF,sizeof(dis));memset(cnt,0,sizeof(cnt));
}
int main()
{int t;int M,W,u,v,w;cin>>t;	while(t--) {init();scanf("%d%d%d",&n,&M,&W);for(int i = 1; i<=M; i++) {scanf("%d%d%d",&u,&v,&w);if(w<maze[u][v])maze[u][v] = maze[v][u] = w;}for(int i = 1; i<=W; i++) {scanf("%d%d%d",&u,&v,&w);maze[u][v] = -w;}cnt[1] =1;if(spfa(1)) printf("YES\n");else printf("NO\n");}return 0 ;
}

AC代码2：（Bellman_Ford算法）（125ms）

//Bellman_Ford算法试试
#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstring>
using namespace std;
#define MAXN 3000*2
#define INF 0xFFFFFFFint t , n , m, w;
int dis[MAXN];
struct Edge{int x;int y;int value;
}e[MAXN];bool judge(){for(int i = 0 ; i < m*2+w ; i++){if(dis[e[i].y] > dis[e[i].x] + e[i].value)return false;}return true;
}void Bellman_Ford(){dis[1] = 0;for(int i = 2 ; i <= n ; i++)dis[i] = INF;for(int i = 1 ; i <= n ; i++){for(int j = 0 ; j < m*2+w; j++){if(dis[e[j].y] > dis[e[j].x] + e[j].value)dis[e[j].y] = dis[e[j].x] + e[j].value;}}if(judge())printf("NO\n");elseprintf("YES\n");
}int main()
{int uu,vv,ww,i;scanf("%d",&t);while(t--){scanf("%d%d%d",&n,&m,&w);for(i=0;i<m*2;){scanf("%d%d%d",&uu,&vv,&ww);e[i].x=uu;e[i].y=vv;e[i++].value=ww;e[i].x=vv;e[i].y=uu;e[i++].value=ww;}for(;i<m*2+w;i++){scanf("%d%d%d",&uu,&vv,&ww);e[i].x=uu;e[i].y=vv;e[i].value=-ww;}Bellman_Ford();}return 0;
}

 AC代码3：（spfa，邻接表储存图）（266ms）

#include<cstdio>
#include<algorithm>
#include<iostream>
#include<queue>
#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;
int n,cnt;
const int MAX = 505;
const int INF = 0x3f3f3f3f;
int dis[MAX],maze[MAX][MAX],cntt[MAX],head[MAX];
bool vis[MAX];
struct Edge {int to,w,ne;Edge(){}//没有此构造函数不能写  node t  这样Edge(int to,int w,int ne):to(to),w(w),ne(ne){}//可以写node(pos,cost)这样} e[200000 + 5];//数组别开小了 
void add(int u,int v,int w) {e[cnt].to = v;e[cnt].w = w;e[cnt].ne = head[u];head[u] = cnt++;
}bool spfa(int s){dis[s]=0; vis[s]=1;  //队列初始化,s为起点int v;queue<int > q;q.push(s);while (!q.empty()){   //队列非空v=q.front();  //取队首元素q.pop();vis[v]=0;   //释放队首结点，因为这节点可能下次用来松弛其它节点，重新入队for(int i=head[v]; i!=-1; i=e[i].ne)  //对所有顶点{if (dis[e[i].to]>dis[v]+e[i].w) {  dis[e[i].to] = dis[v]+e[i].w;   //修改最短路if (vis[e[i].to]==0) {  //如果扩展结点i不在队列中，入队cntt[e[i].to]++;vis[e[i].to]=1;q.push(e[i].to);if(cntt[e[i].to] >=n) return true;}}}}return false;
}void init() {cnt = 0;memset(vis,0,sizeof(vis));memset(maze,0,sizeof(maze));memset(e,0,sizeof(e));memset(dis,INF,sizeof(dis));memset(cntt,0,sizeof(cntt));memset(head,-1,sizeof(head)); 
}
int main()
{int t;int M,W,u,v,w;cin>>t;	while(t--) {init();scanf("%d%d%d",&n,&M,&W);for(int i = 1; i<=M; i++) {scanf("%d%d%d",&u,&v,&w);add(u,v,w);add(v,u,w);}for(int i = 1; i<=W; i++) {scanf("%d%d%d",&u,&v,&w);add(u,v,-w);}cntt[1] =1;if(spfa(1)) printf("YES\n");else printf("NO\n");}return 0 ;
}

AC代码4：（Bellmanford，加优化）（47ms）

//bellman算法：
#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#define N 20100using namespace std;
const int INF = 0x3f3f3f3f;
int n,m,mm;
int tot;
int dis[6005];
struct Edge {int u,v,w;
} e[6005];
void add(int u,int v,int w) {e[++tot].u = u;e[tot].v = v;e[tot].w = w;
}
bool bell() {memset(dis,INF,sizeof(dis));dis[1]=0;int all = n;int flag = 0;while(all--) {flag = 0;for(int i = 1; i<=tot; i++) {if(dis[e[i].v] > dis[e[i].u] + e[i].w) {flag = 1;dis[e[i].v] = dis[e[i].u] + e[i].w;}}if(flag == 0) return false;}flag = 0;for(int i = 1; i<=tot; i++) {if(dis[e[i].v] > dis[e[i].u] + e[i].w) return 1;}return 0;
}int main()
{int f,a,b,c;cin>>f;while(f--) {tot = 0;scanf("%d%d%d",&n,&m,&mm);for(int i = 1; i<=m; i++) {scanf("%d%d%d",&a,&b,&c) ;add(a,b,c);add(b,a,c);}for(int i = 1; i<=mm; i++) {scanf("%d%d%d",&a,&b,&c);add(a,b,-c);}if(bell()) puts("YES");else puts("NO");}return 0 ;
}
//2
//3 3 1
//1 2 2
//1 3 4
//2 3 1
//3 1 3