vijos1909 寻找道路 [SPFA][dfs][bfs]

描述
在有向图 G 中，每条边的长度均为 1，现给定起点和终点，请你在图中找一条从起点到 终点的路径，该路径满足以下条件：
路径上的所有点的出边所指向的点都直接或间接与终点连通。
在满足条件 1 的情况下使路径最短。
注意：图 G 中可能存在重边和自环，题目保证终点没有出边。 请你输出符合条件的路径的长度。
格式
输入格式

第一行有两个用一个空格隔开的整数 n 和 m，表示图有 n 个点和 m 条边。
接下来的 m 行每行 2 个整数 x、y，之间用一个空格隔开，表示有一条边从点 x 指向点y。
最后一行有两个用一个空格隔开的整数 s、t，表示起点为 s，终点为 t。
输出格式

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

3 2
1 2
2 1
1 3
样例输出1[复制]

-1
样例2
样例输入2[复制]

6 6
1 2
1 3
2 6
2 5
4 5
3 4
1 5
样例输出2[复制]

3
限制
对于 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。
提示
【输入输出样例1说明】

如上图所示，箭头表示有向道路，圆点表示城市。起点 1 与终点 3 不连通，所以满足题目描述的路径不存在，故输出-1。
【输入输出样例2说明】

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

题意：
对于给定的有向图和开始结束节点走最短路，对于经过的点有如下要求：从该点出发所有能到的点都要直接或间接与终点相连；
分析：
显然走一个dfs判断可行点，再走一次SPFA就可以了。
当然各位大佬们都是建反图从终点bfs的再正向bfs的，个人觉得作为一个有向图，直接从起点dfs不就行了？如果不行的话求大佬指点。
这个程序我写过，这个博客是改题记录。
为什么要改，因为vijos是个好东西，我随手看到我居然WA了这道题，找了后来改的随手一交，居然T了。
当初把head2设为1，没想到tail2还可以为0,不严密，导致一直没有跳出循环，vijos卡得好QAQ。

#include<cstdio>
#include<iostream>
#include<cstring>
using namespace std;
const int MAXN=200005;
int tov[MAXN],tov_[MAXN],nex[MAXN],nex_[MAXN],fr[MAXN],fr_[MAXN]; 
int z[MAXN],head,tail,dis[10005],n,e,x,y,ca[10005],vis[10005];
int z2[10005],head2,tail2,sta,des;
void dfs(int k)
{
    int flg=fr_[k];ca[k]=1;vis[k]=1;
    while(flg!=0)
    {
        if(!vis[tov_[flg]])dfs(tov_[flg]);
        flg=nex_[flg];
    }
}
int main()
{
    freopen("road.in","r",stdin);
    freopen("road.out","w",stdout);
    scanf("%d %d",&n,&e);
    for(int i=1;i<=e;i++)
    {
        scanf("%d %d",&x,&y);
        nex[i]=fr[x];tov[i]=y;fr[x]=i;
        nex_[i]=fr_[y];tov_[i]=x;fr_[y]=i;
    }
    scanf("%d %d",&sta,&des);
    dfs(des);
    head++;tail++;z[tail]=sta;
    memset(dis,0x7f,sizeof(dis));
    dis[sta]=0;
    for(int i=1;i<=n;i++)
        if(!ca[i])z2[++tail2]=i;
    head2=1;
    while(head2<=tail2)
    {int flg=fr_[z2[head2]];
    while(flg!=0)
    {
        if(ca[tov_[flg]])
        {
        if(tail2>=10004){
        tail2%=10004;tail2-=1;}
        ca[tov_[flg]]=0;
        }
        flg=nex_[flg];
    }
    head2++;if(head2>=10005)head2%=10005;   
    }   
    while(head!=tail+1)
    {int flg=fr[z[head]];
    while(flg!=0)
    {
        if(ca[tov[flg]])
        if(dis[tov[flg]]>dis[z[head]]+1)
        {
            z[++tail]=tov[flg];
            if(tail>=MAXN-1){
            tail%=(MAXN-1);tail-=1;}
            dis[tov[flg]]=dis[z[head]]+1;
        }
        flg=nex[flg];
    }
    head++; 
    if(head>=MAXN)head%=MAXN;
    }
    if(dis[des]<=2000000)printf("%d",dis[des]);
    else printf("-1");
    return 0;
}