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=