HDU 1269 迷宫城堡（tarjan求连通分量)

题目描述

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

输入描述:

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

输出描述:

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

输入例子:

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

输出例子:

Yes
No

思路来自：http://www.cnblogs.com/shadowland/p/5872257.html

---

//运行时间：78ms 占用内存：3252k
#include <iostream>
#include <stdio.h>
#include <vector>
#include <stack>
#define MAXN 10000
#define MAXM 100000
using namespace std;
int n;
stack<int> sta;         //深搜访问栈
vector<int> road[MAXN]; //road[i][j] 从第i个房间可以直接到第i个房间
vector<int> scon[MAXN]; //scon[i][j] 编号为j的点属于第i个连通分量
int dfn[MAXN];          //dfn[i] 的值表示第i个点的访问时间
int low[MAXN];          //形成的回路中的所有的点中最早被访问到得点的访问时间
int is_vis[MAXN];      //is_vis[i] 节点i是否在栈中, 2在， 1已访问， 0未访问
int time, scon_num;     //时间,强连通子图数
int min(int a, int b) {
    return a > b ? b : a;
}
void tarjan(int x) {
    int i, size, d;
    sta.push(x);
    is_vis[x] = 2;
    dfn[x] = ++time;
    low[x] = time;
    size = road[x].size();
    for (i = 0; i < size; i++) {
        d = road[x][i];
        if (!dfn[d]) {
            tarjan(d);
            low[x] = min(low[x], low[d]);
        }
        else if (is_vis[d] == 2) {
            low[x] = min(low[x], dfn[d]);
        }
    }
    //如果dfn 和 low相等，可以构成一个连通分量
    if (dfn[x] == low[x]) {
        while (!sta.empty()) {
            d = sta.top();
            sta.pop();
            scon[scon_num].push_back(d);
            is_vis[d] = 1;
            if (d == x)
                break;
        }
        ++scon_num;
    }
}
void solve()
{
    int i;
    //如果没有访问到，从该点出发深搜
    for (i = 1; i <= n; i++) {
        if (!dfn[i])
            tarjan(i);
    }
    if (scon_num > 1)
        printf("No\n");
    else
        printf("Yes\n");
}
int main()
{
    int m, a, b, i;
    while (scanf("%d%d", &n, &m), n) {
        //初始化操作
        scon_num = 0;
        time = 0;
        for (i = 1; i <= n; i++) {
            road[i].clear();
            scon[i].clear();
        }
        memset(low, 0, sizeof low);
        memset(dfn, 0, sizeof dfn);
        memset(is_vis, 0, sizeof is_vis);
        while (!sta.empty())
            sta.pop();
        //输入
        for (i = 0; i < m; i++) {
            scanf("%d%d", &a, &b);
            road[a].push_back(b);
        }
        solve();
    }
    return 0;
}