CODEVS 3731寻找道路 (重修版)

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

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

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

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

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

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

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

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

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

输出描述 Output Description
输出文件名为road.out。

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

路径上的所有点的出边所指向的点都直接或间接与终点连通。
对于这句话，我们该如何处理呢？
先反向建边，从终点跑一边SPFA，记录每个点能不能被终点跑到。
然后，我们需要对这些跑不到的点进行处理，跑反向图，从这些点只进行一层BFS，把他们连接的点都标为不可走。
从起点再跑SPFA，遇到不能走的点就continue掉。

#include<cstdio>
#include<queue>
using namespace std;
const int inf=0x7ffffff;
const int maxn=400000+500;
int n,m;
struct Edge
{
    int f;
    int to;
    int d;
    int next;
}edge1[maxn],edge[maxn];

int head[maxn],head1[maxn];
bool vis[maxn],alive[maxn];
bool kd[maxn];
int dist[maxn];
int ss,tt;
int tot;
void add(int f,int t,int d)
{
    edge[++tot]=(Edge){f,t,d,head[f]};
    head[f]=tot;    
}
int tot1;
void add1(int f,int t,int d)
{
    edge1[++tot1]=(Edge){f,t,d,head1[f]};
    head1[f]=tot1;  
}
queue<int>q1;
void spfa()
{
    alive[tt]=1;
    q1.push(tt);
    while(!q1.empty())
    {
        int x=q1.front();
        q1.pop();
        for(int i=head1[x];i;i=edge1[i].next)
        {
            Edge e=edge1[i];
            if(!vis[e.to])
            {
                vis[e.to]=1;
                alive[e.to]=1;
                kd[e.to]=1;
                q1.push(e.to);
            }
        }
    }
}
void kill(int x)
{
    for(int i=head1[x];i;i=edge1[i].next)
    {
        Edge e=edge1[i];
        if(e.to==ss)
            continue;
        if(alive[e.to])
            alive[e.to]=0;
    }
}
queue<int>q;
void spfa1()
{
    for(int i=1;i<=n;i++)
        dist[i]=inf,vis[i]=0;   
    dist[ss]=0;
    q.push(ss);
    while(!q.empty())
    {
        int x=q.front();
        q.pop();
        vis[x]=0;
        for(int i=head[x];i;i=edge[i].next)
        {
            Edge e=edge[i];
            if(alive[e.to])
            {
                if(dist[e.to]>dist[x]+e.d)
                {
                    dist[e.to]=dist[x]+e.d;
                    if(!vis[e.to])
                    {
                        vis[e.to]=1;
                        q.push(e.to); 
                    }
                }
            }   
        }
    }
}

int main()
{
    scanf("%d%d",&n,&m);
    for(int i=1;i<=m;i++)
    {
        int a,b;
        scanf("%d%d",&a,&b);
        add(a,b,1);
        add1(b,a,1);
    }
    scanf("%d%d",&ss,&tt);

    spfa();

    for(int i=1;i<=n;i++)
    {
        if(!alive[i]&&!kd[i])
        kill(i);
    }
    spfa1();

    if(dist[tt]!=inf)
        printf("%d\n",dist[tt]);
    else
        puts("-1"); 
    return 0;
}