1715: [Usaco2006 Dec]Wormholes 虫洞

1715: [Usaco2006 Dec]Wormholes 虫洞

Time Limit: 5 Sec  Memory Limit: 64 MB
Submit: 501  Solved: 278
[Submit][Status][Discuss]

Description

John在他的农场中闲逛时发现了许多虫洞。虫洞可以看作一条十分奇特的有向边，并可以使你返回到过去的一个时刻（相对你进入虫洞之前）。John的每个农场有M条小路（无向边)连接着N （从1..N标号）块地，并有W个虫洞。其中1<=N<=500,1<=M<=2500,1<=W<=200。 现在John想借助这些虫洞来回到过去（出发时刻之前），请你告诉他能办到吗。 John将向你提供F(1<=F<=5)个农场的地图。没有小路会耗费你超过10000秒的时间，当然也没有虫洞回帮你回到超过10000秒以前。

Input

* Line 1: 一个整数 F, 表示农场个数。

* Line 1 of each farm: 三个整数 N, M, W。

* Lines 2..M+1 of each farm: 三个数（S, E, T）。表示在标号为S的地与标号为E的地中间有一条用时T秒的小路。

* Lines M+2..M+W+1 of each farm: 三个数（S, E, T）。表示在标号为S的地与标号为E的地中间有一条可以使John到达T秒前的虫洞。

Output

* Lines 1..F: 如果John能在这个农场实现他的目标，输出"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

 

Source

Gold

题解：很明显是找负权回路，可是我还一直没写过，直到今天第一次写spfa找负权回路——bfs的spfa只要判断下每个节点的入队次数，只要大于N，就可以判断有负环

 /**************************************************************
     Problem:
     User: HansBug
     Language: Pascal
     Result: Accepted
     Time: ms
     Memory: kb
 ****************************************************************/

 type
     point=^node;
     node=record
                g,w:longint;
                next:point;
     end;
 var
    i,j,k,l,m,n,t:longint;
    a:array[..] of point;
    b,c,g,e:array[..] of longint;
    d:array[..] of longint;
 procedure add(x,y,z:longint);
           var p:point;
           begin
                new(p);p^.g:=y;p^.w:=z;p^.next:=a[x];a[x]:=p;
           end;

 function spfa(x:longint):boolean;
          var f,r:longint;p:point;
          begin
               fillchar(e,sizeof(e),);
               fillchar(b,sizeof(b),);
               fillchar(g,sizeof(g),);
               for i:= to n do c[i]:=maxlongint div ;
               f:=;r:=;g[x]:=;d[]:=x;c[x]:=;b[x]:=;e[x]:=;
               while f<>r do
                     begin
                          p:=a[d[f]];
                          while p<>nil do
                                begin
                                     if (e[p^.g]=) or (c[p^.g]>(c[d[f]]+p^.w)) then
                                        begin
                                             e[p^.g]:=;
                                             c[p^.g]:=c[d[f]]+p^.w;
                                             if g[p^.g]= then
                                                begin
                                                     inc(b[p^.g]);
                                                     if b[p^.g]>n then exit(true);
                                                     g[p^.g]:=;
                                                     d[r]:=p^.g;
                                                     r:=r mod +;
                                                end;
                                        end;
                                     p:=p^.next;
                                end;
                          g[d[f]]:=;f:=f mod +;
                     end;
               exit(false);
          end;
 begin
      readln(t);
      while not(eof) do
          begin
               readln(n,m,t);
               for i:= to n do a[i]:=nil;
               for i:= to m do
                   begin
                        readln(j,k,l);
                        add(j,k,l);add(k,j,l);
                   end;
               for i:= to t do
                   begin
                        readln(j,k,l);
                        add(j,k,-l);
                   end;
               if spfa() then writeln('YES') else writeln('NO');
          end;
 end.