hdu1269-图论-Tarjan-迷宫城堡

迷宫城堡

Problem Description

为了训练小希的方向感，Gardon建立了一座大城堡，里面有N个房间(N<=10000)和M条通道(M<=100000)，每个通道都是单向的，就是说若称某通道连通了A房间和B房间，只说明可以通过这个通道由A房间到达B房间，但并不说明通过它可以由B房间到达A房间。Gardon需要请你写个程序确认一下是否任意两个房间都是相互连通的，即：对于任意的i和j，至少存在一条路径可以从房间i到房间j，也存在一条路径可以从房间j到房间i。

 

Input

输入包含多组数据，输入的第一行有两个数：N和M，接下来的M行每行有两个数a和b，表示了一条通道可以从A房间来到B房间。文件最后以两个0结束。

 

Output

对于输入的每组数据，如果任意两个房间都是相互连接的，输出"Yes"，否则输出"No"。

 

Sample Input

  
  

   
   3 3
1 2
2 3
3 1
3 3
1 2
2 3
3 2
0 0
  
  

 

Sample Output

  
  

   
   Yes
No
  
  

#include<stdio.h>
#include<vector>
#include<stack>
#include<string.h>
#include<iostream>
#define maxn 10010
using namespace std;
vector<int> G[maxn];        //用来存图
//pre[u]为搜素节点u的次序编号，lowlink[u]为u及其后代能追溯到的（最先被发现的）祖先v的pre（v）的值
//sccno[i]为i所在的SCC编号，scc_cnt为SCC计数器
int pre[maxn],lowlink[maxn],sccno[maxn],dfs_clock,scc_cnt;
stack<int> S;   //保存当前SCC中的结点
int N;
void dfs(int u)
{
    pre[u]=lowlink[u]=++dfs_clock;      //为结点u设定次序编号和lowlink初值
    S.push(u);      //将结点u压入栈中
    for(int i=0;i<G[u].size();i++)      //枚举每一条边
    {
        int v=G[u][i];
        if(!pre[v])             //如果结点u未被访问过
        {
            dfs(v);             //继续向下找
            lowlink[u]=min(lowlink[u],lowlink[v]);
        }
        else if(!sccno[v])      //如果v点在栈中,说明已经形成环了，那么u的lowlink值为（其本身和v点的pre[v]中较小的那个）
        {
            lowlink[u]=min(lowlink[u],pre[v]);
        }
    }
    if(lowlink[u]==pre[u])      //结点u是强连通分量的根，（找到一个SCC）
    {
        scc_cnt++;      //用于记录SCC的个数
        for(;;)         //将一个强连通分量中所有的结点从栈中取出
        {
            int x=S.top();S.pop();
            sccno[x]=scc_cnt;
            if(x==u)break;
        }
    }
}
void find_scc(int n)
{
    dfs_clock=scc_cnt=0;
    memset(sccno,0,sizeof(sccno));
    memset(pre,0,sizeof(pre));
    for(int i=1;i<=n;i++)
    {
        if(!pre[i])
            dfs(i);
    }
}
int main()
{
    int M,i,a,b,flag;
    while(scanf("%d%d",&N,&M)!=EOF&&(N||M))
    {
        flag=0;
        while(!S.empty())
        S.pop();
        for(i=1;i<=N;i++)
        G[i].clear();
        for(i=0;i<M;i++)
        {
            scanf("%d%d",&a,&b);
            G[a].push_back(b);
        }
        find_scc(N);
        int temp=sccno[1];
        for(i=2;i<=N;i++)
        {
            if(sccno[i]!=temp)
            {
                flag=1;
                break;
            }
        }
//        for(i=1;i<=N;i++)
//            printf("%d ",lowlink[i]);
//        printf("\n");
//       printf("%d\n",scc_cnt);
        if(flag==1)
            printf("No\n");
        else
            printf("Yes\n");
    }
    return 0;
}
/*
6 8
1 3
3 5
5 6
3 4
4 6
1 2
2 4
4 1
*/