[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.