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: 1715
    User: HansBug
    Language: Pascal
    Result: Accepted
    Time:160 ms
    Memory:1060 kb
****************************************************************/
 
type
    point=^node;
    node=record
               g,w:longint;
               next:point;
    end;
var
   i,j,k,l,m,n,t:longint;
   a:array[0..1500] of point;
   b,c,g,e:array[0..1500] of longint;
   d:array[0..100000] 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),0);
              fillchar(b,sizeof(b),0);
              fillchar(g,sizeof(g),0);
              for i:=1 to n do c[i]:=maxlongint div 2;
              f:=1;r:=2;g[x]:=1;d[1]:=x;c[x]:=0;b[x]:=1;e[x]:=1;
              while f<>r do
                    begin
                         p:=a[d[f]];
                         while p<>nil do
                               begin
                                    if (e[p^.g]=0) or (c[p^.g]>(c[d[f]]+p^.w)) then
                                       begin
                                            e[p^.g]:=1;
                                            c[p^.g]:=c[d[f]]+p^.w;
                                            if g[p^.g]=0 then
                                               begin
                                                    inc(b[p^.g]);
                                                    if b[p^.g]>n then exit(true);
                                                    g[p^.g]:=1;
                                                    d[r]:=p^.g;
                                                    r:=r mod 100000+1;
                                               end;
                                       end;
                                    p:=p^.next;
                               end;
                         g[d[f]]:=0;f:=f mod 100000+1;
                    end;
              exit(false);
         end;
begin
     readln(t);
     while not(eof) do
         begin
              readln(n,m,t);
              for i:=1 to n do a[i]:=nil;
              for i:=1 to m do
                  begin
                       readln(j,k,l);
                       add(j,k,l);add(k,j,l);
                  end;
              for i:=1 to t do
                  begin
                       readln(j,k,l);
                       add(j,k,-l);
                  end;
              if spfa(1) then writeln('YES') else writeln('NO');
         end;
end.