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.

巴特西

1715: [Usaco2006 Dec]Wormholes 虫洞