[BZOJ1715][Usaco2006 Dec]Wormholes 虫洞

[Usaco2006 Dec]Wormholes 虫洞

时间限制: 1 Sec 内存限制: 128 MB

题目描述

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

输入

• 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秒前的虫洞。
  输出
• Lines 1..F: 如果John能在这个农场实现他的目标，输出”YES”,否则输出”NO”。

样例输入

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

样例输出

NO
YES

SPFA判负环即可，据说写成BFS会快不少，有时间补上…

var w:array[0..100000,1..3]of longint; dist:array[0..10000]of longint; p:array[0..100000,1..2]of longint; t:array[0..1000000]of longint; i,j,l:longint; a,b,c,x,y,tt,k,f,key,v:longint; n,m,start,finish:longint; head,tail:longint;begin readln(f); for l:=1 to f do begin for i:=1 to k do  begin   w[i,1]:=0;   w[i,2]:=0;   w[i,3]:=0;  end; readln(n,m,v); k:=n+1; for i:=1 to m do  begin   readln(a,b,c);   w[k,1]:=b;   w[k,2]:=c;   if w[a,3]=0   then w[a,3]:=k   else w[w[a,1],3]:=k;   w[a,1]:=k;   inc(k);   w[k,1]:=a;   w[k,2]:=c;   if w[b,3]=0   then w[b,3]:=k   else w[w[b,1],3]:=k;   w[b,1]:=k;   inc(k);  end; for i:=1 to v do  begin   readln(a,b,c);   w[k,1]:=b;   w[k,2]:=-c;   if w[a,3]=0   then w[a,3]:=k   else w[w[a,1],3]:=k;   w[a,1]:=k;   inc(k);  end; start:=1; finish:=n; for i:=1 to n do  dist[i]:=100000; for i:=1 to tail do  begin   p[i,1]:=0;   p[i,2]:=0;  end; dist[start]:=0; head:=1; tail:=2; t[1]:=start; p[start,1]:=1; p[start,2]:=1; key:=0; while head<tail do  begin   x:=head; y:=tail;   for i:=head to tail-1 do    begin     if p[t[head],2]>=n     then begin key:=1; break; end;     tt:=w[t[i],3];     while tt<>0 do      begin       if dist[t[i]]+w[tt,2]<dist[w[tt,1]]       then        begin         dist[w[tt,1]]:=dist[t[i]]+w[tt,2];         if p[w[tt,1],1]=0         then begin t[y]:=w[tt,1]; p[w[tt,1],1]:=1; inc(p[w[tt,1],2]); inc(y); end;        end;       tt:=w[tt,3];      end;     inc(x);     p[t[i],1]:=0;    end;   head:=x; tail:=y;   if key=1 then break;  end; if key=0 then writeln('NO') else writeln('YES'); end;end.