noip2014寻找道路

题目描述

在有向图G 中，每条边的长度均为1 ，现给定起点和终点，请你在图中找一条从起点到终点的路径，该路径满足以下条件：

1 ．路径上的所有点的出边所指向的点都直接或间接与终点连通。

2 ．在满足条件1 的情况下使路径最短。

注意：图G 中可能存在重边和自环，题目保证终点没有出边。

请你输出符合条件的路径的长度。

输入输出格式

输入格式：

输入文件名为road .in。

第一行有两个用一个空格隔开的整数n 和m ，表示图有n 个点和m 条边。

接下来的m 行每行2 个整数x 、y ，之间用一个空格隔开，表示有一条边从点x 指向点y 。

最后一行有两个用一个空格隔开的整数s 、t ，表示起点为s ，终点为t 。

输出格式：

输出文件名为road .out 。

输出只有一行，包含一个整数，表示满足题目᧿述的最短路径的长度。如果这样的路径不存在，输出- 1 。

输入输出样例

输入样例#1：

3 2
1 2
2 1
1 3

输出样例#1：

-1

输入样例#2：

6 6
1 2
1 3
2 6
2 5
4 5
3 4
1 5

输出样例#2：

3

说明

解释1：



如上图所示，箭头表示有向道路，圆点表示城市。起点1 与终点3 不连通，所以满足题

目᧿述的路径不存在，故输出- 1 。

解释2：



如上图所示，满足条件的路径为1 - >3- >4- >5。注意点2 不能在答案路径中，因为点2连了一条边到点6 ，而点6 不与终点5 连通。

对于30%的数据，0<n≤10，0<m≤20；

对于60%的数据，0<n≤100，0<m≤2000；

对于100%的数据，0<n≤10,000，0<m≤200,000，0<x，y，s，t≤n，x≠t。

脑残起来我自己都害怕 orz orz  O|￣|_

第一遍做得时候正着dfs判断一个点是不是所有连的点都可以到终点成功的wa9个点 orz，

因为正着dfs在你搜到x点的时候你给vis[x]=true然后它是否符合要求的标签还没有更新仍然是false

然后你去访问它所连的点，但是如果是个环的话它所连的点又连到了它然后就都被置成false QAQ

此方法扑街 



然后就需要换个思路，既然要求都到终点，那就先只建反边，从终点往回搜，记录一个点访问的次数，

如果一个点的出度（正边的时候）等于反向dfs搜索时的访问次数说明满足条件，那么f[i]=true;

同时我们要删去出度为0的点,此时会把终点也删去，所以更新完f[i]数组后单独的 f[终点]=true

那么此时我们有哪些点可以选入最短路就已经打好标签了

然后我们建上正边，跑spfa就可以了

注意：

一定要在把反边都删了以后再建正边，不然spfa会挂（因为会走原本不存在的反边！！！）

邻接表删边的方法就是清零last数组并把边数l清0

输出时注意对-1的判断就ok了=w=

因为没删反边给跪了 QAQ

var
n,m,l,x,y,s,t :longint;
i :longint;
que,dis :array[0..10010] of longint;
last,outdu,time :array[0..10010] of longint;
xx,yy :array[0..200010] of longint;
vis,flag :array[0..10010] of boolean;
pre,other :array[0..400010] of longint;
procedure connect(x,y:longint);
begin
inc(l);
pre[l]:=last[x];
last[x]:=l;
other[l]:=y;
end;

procedure dfs(x:longint);
var
p,q:longint;
begin
q:=last[x];
vis[x]:=true;
while (q<>0) do
begin
p:=other[q];
inc(time[p]);
if not vis[p] then dfs(p);
q:=pre[q];
end;
end;

procedure spfa;
var
p,q,cur,h,tl:longint;
begin
fillchar(vis,sizeof(vis),false);
fillchar(dis,sizeof(dis),127);
h:=0;tl:=1;que[1]:=s;dis[s]:=0;
while (h<>tl) do
begin
h:=h mod 10000+1;
cur:=que[h];
vis[cur]:=false;
q:=last[cur];
while (q<>0) do
begin
p:=other[q];
if (dis[p]>dis[cur]+1) and flag[p] then
begin
dis[p]:=dis[cur]+1;
if not vis[p] then
begin
vis[p]:=true;
tl:=tl mod 10000+1;
que[tl]:=p;
end;
end;
q:=pre[q];
end;
end;
end;

begin
read(n,m);
for i:=1 to m do
begin
read(x,y);
inc(outdu[x]);
xx[i]:=x;yy[i]:=y;
connect(y,x);
end;
read(s,t);
dfs(t);
for i:=1 to n do
if (outdu[i]=time[i]) and (outdu[i]<>0) then flag[i]:=true;
flag[t]:=true;
if not flag[s] then writeln(-1) else
begin
for i:=1 to n do last[i]:=0;l:=0;
for i:=1 to m do connect(xx[i],yy[i]);
spfa;
if dis[t]=dis[0] then writeln(-1) else writeln(dis[t]);
end;
end.

——by Eirlys
转载请注明出处=w=