noip2014-寻找道路

题目很清楚，要满足2个要求：

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

2.最短

条件2很好做，spfa或者bfs反正怎么搞都行。

于是考虑条件1

我们当然想知道每个点能否经过。对于不能经过的点，无非是2种：

1.自身到不了终点的 2.它的出边指向的点有到不了终点的。

首先找满足条件1的，然后把所有指向它们的点都标记为0就成

对于条件1，我们建反边……从t开始搜索，所有能到达的点就是1不然就是0

大概就是这样……

注意！！！环

这里必须要建反边，不能正向搜索（如果v[i]到的了的话，x也能到），因为有环！！！一个结点是可以经过多次的！！！

而这种正向搜索对一个有向图而言只有存在拓扑序时才能成立。

emmmm差不多咯

 

后记：经过N次的摸索（https://blog.csdn.net/Mmm040403/article/details/83687775）

采用正向搜索的方法AC啦

方法：Tarjan强连通分量

先缩点，按有向无环图遍历。

对于同一个连通分量中的点，只要有一个点到达，则连通分量中所有点全可以到达。

#include <iostream> 
#include <cstdio>
#include <vector>
#include <queue>
#include <cstring>
#include <stack>

using namespace std;
stack<int> sa; 
const int maxn=200005;
const int INF=0x3f3f3f3f;
int n,m,s,t,indexs;
vector<int> v[maxn],p[maxn],vi[maxn];int ins[maxn];
int vis[maxn],can[maxn],ok[maxn],dis[maxn];
int dfn[maxn],low[maxn],belong[maxn];
queue<int> q;
void tarjan(int x)
{
    sa.push(x);ins[x]=1;
    dfn[x]=low[x]=++indexs;
    for(int i=0;i<v[x].size();i++)
    {
        if(!dfn[v[x][i]]) 
        {
            tarjan(v[x][i]);
            low[x]=min(low[x],low[v[x][i]]);
        }
        else if(ins[v[x][i]]) low[x]=min(dfn[v[x][i]],low[x]);
    }
    if(dfn[x]==low[x])
    {
        while(1)
        {
            int tmp=sa.top();sa.pop();ins[tmp]=0;
            belong[tmp]=x;
            if(tmp==x) break;
        }
    }
}
void find(int x)
{
    vis[x]=1;
    for(int i=0;i<vi[x].size();i++)
    {
        if(!vis[vi[x][i]]) find(vi[x][i]);
        if(can[vi[x][i]]) can[x]=1;
    }
}
int main()
{
    scanf("%d%d",&n,&m);
    for(int i=1;i<=m;i++)
    {
        int x,y;scanf("%d%d",&x,&y);
        if(x!=y) {v[x].push_back(y);p[y].push_back(x);}
    } 
    scanf("%d%d",&s,&t);can[t]=1;
    for(int i=1;i<=n;i++) if(!dfn[i]) tarjan(i);
    for(int i=1;i<=n;i++)
    {
        for(int j=0;j<v[i].size();j++)
        {
            if(belong[i]!=belong[v[i][j]])	
            {
                vi[belong[i]].push_back(belong[v[i][j]]);	
            }	
        }
    } 
    can[belong[t]]=1;find(belong[s]);
    for(int i=1;i<=n;i++) if(can[belong[i]]) can[i]=1;
//	for(int i=1;i<=n;i++) cout<<can[i]<<" ";
//	cout<<endl;
    for(int i=1;i<=n;i++) ok[i]=1;
    for(int i=1;i<=n;i++)
    {
        if(!can[i]) 
        {
            ok[i]=0;
            for(int j=0;j<p[i].size();j++) ok[p[i][j]]=0;
        }
    }
    if(!ok[s]) {printf("-1");return 0;}
    memset(vis,0,sizeof(vis));  
    memset(dis,INF,sizeof(dis));
    while(!q.empty()) q.pop();
    q.push(s);vis[s]=1;dis[s]=0;
    while(!q.empty())
    {
        int x=q.front();q.pop();
        for(int i=0;i<v[x].size();i++)
        {
            if(!ok[v[x][i]]||vis[v[x][i]]) continue;
            dis[v[x][i]]=dis[x]+1;
            q.push(v[x][i]);vis[v[x][i]]=1;
        } 
    }
    if(dis[t]!=INF) printf("%d",dis[t]);
    else printf("-1");
    return 0;
}