[NOIP2014][建图]寻找道路

题目描述：
题目背景: NOIP2014 提高组 Day2 T2。
题目链接： UOJ19 http://uoj.ac/problem/19
在有向图 G 中，每条边的长度均为 1，现给定起点和终点，请你在图中找一条从起点到终点的路径，该路径满足以下条件：
1．路径上的所有点的出边所指向的点都直接或间接与终点连通。
2．在满足条件 1 的情况下使路径最短。
注意：图 G 中可能存在重边和自环，题目保证终点没有出边。
请你输出符合条件的路径的长度。
输入格式：
第一行有两个用一个空格隔开的整数 n 和 m，表示图有 n 个点和 m 条边。
接下来的 m 行每行 2 个整数 x、y，之间用一个空格隔开，表示有一条边从点 x 指向点y。
最后一行有两个用一个空格隔开的整数 s、t，表示起点为 s，终点为 t。
输出格式：
输出只有一行，包含一个整数，表示满足题目描述的最短路径的长度。如果这样的路径不存在，输出 -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。$
题目分析：
不是很难，想想就可以知道。我们考虑，题目要求路径上的点不能连向那些与终点不联通的点。于是反向建图，从终点dfs或bfs,没有搜到的点就为不合法的点。然后以这些点为基础，它们在反向图中所指向的点都是不合法的点。最后以合法的点正向建图，跑一边最短路就可以得到答案了。需要注意的是，两种不合法的点标记不能一样，我是将没有搜的点的标记为0，搜过的标记为1，标记为0的点所连向的点标记为2（如果为0，则在枚举时你可能又会把它们所连向的点标记为不合法的点，而这显然是不正确的）。
PS：我附的是UOJ的链接，通常都会有一些hack数据，可以更大程度上检验自己的程序。
附代码：

#include<iostream>
#include<cstring>
#include<string>
#include<cstdlib>
#include<cstdio>
#include<ctime>
#include<queue>
#include<set>
#include<map>
#include<cmath>
#include<cctype>
#include<iomanip>
#include<algorithm>
using namespace std;

const int maxn=1e4+10;
const int maxm=2e5+10;
const int INF=0x3f3f3f3f;
int n,m,tot,first[maxn],nxt[maxm],to[maxm],x[maxm],y[maxm];
int dis[maxn],vis[maxn],s,t;
priority_queue<pair<int,int> > q;

int readint()
{
    char ch;int i=0,f=1;
    for(ch=getchar();(ch<'0'||ch>'9')&&ch!='-';ch=getchar());
    if(ch=='-') {ch=getchar();f=-1;}
    for(;ch>='0'&&ch<='9';ch=getchar()) i=(i<<3)+(i<<1)+ch-'0';
    return i*f;
}

void create(int x,int y)
{
    tot++;
    nxt[tot]=first[x];
    first[x]=tot;
    to[tot]=y;
}

void dfs(int u)
{
    for(int e=first[u];e;e=nxt[e])
        if(vis[to[e]]==0)
        {
            vis[to[e]]=1;
            dfs(to[e]);
        }
}

void dijkstra()//dijkstra+优先队列跑最短路
{
    memset(dis,INF,sizeof(dis));
    dis[s]=0;
    q.push(make_pair(0,s));
    while(!q.empty())
    {
        int u=q.top().second;
        q.pop();
        for(int e=first[u];e;e=nxt[e])
        {
            int v=to[e];
            if(dis[u]+1<dis[v])
            {
                dis[v]=dis[u]+1;
                if(v==t) return;
                q.push(make_pair(-dis[v],v));
            }
        } 
    }
}

int main()
{
    //freopen("lx.in","r",stdin);

    n=readint();m=readint(); 
    for(int i=1;i<=m;i++)//反向建图
    {
        x[i]=readint();y[i]=readint();
        create(y[i],x[i]);
    }
    s=readint();t=readint();
    vis[t]=1;
    dfs(t);
    for(int i=1;i<=n;i++)
        if(vis[i]==0)
        {
            for(int e=first[i];e;e=nxt[e])
                if(vis[to[e]]!=0)//注意要加这个判断，不能把那些原先为0的点修改了
                    vis[to[e]]=2;
        }
    if(vis[s]!=1||vis[t]!=1)
    {
        printf("-1");
        return 0;
    }
    memset(first,0,sizeof(first));
    memset(nxt,0,sizeof(nxt));
    memset(to,0,sizeof(to));
    tot=0;
    for(int i=1;i<=m;i++)//合法点正向建图
        if(vis[x[i]]==1&&vis[y[i]]==1)
        create(x[i],y[i]);
    dijkstra();
    if(dis[t]==INF) printf("-1");
    else printf("%d",dis[t]);

    return 0;
}