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

题目的大致意思就是给一个图，判断图中任何两点能否强连通，或许你会想到用最短路的方法求解，floyd求最短路，如果两个点的距离为 inf （无穷）则无法连通，但是数据量太大了，用最短路求解会超时，本题目可以用强连通分量（tarjan算法）求解。

什么是强连通？

举个例子，在一个有向图中，a 能抵达 b，b 也能抵达 a，则称 a 和 b 强连通

那么对于一个图而言，在一个图中，有若干个子图，如果一个子图中的所有点之间都强连通，就称为强连通子图。

强连通分量则是指一个图中的极大强连通子图

很显然对于属于任何一个强连通子图的任何两个点，一定可以互相到达，这也是本题要让我们求解的东西

既然要判断图中的任意两个点是否都能强连通，只要判断图中的所有点是否在同一个强连通分量即可。（或者说判断图中强连通分量的个数，如果只有一个强连通分量，那么说明所有点属于同一个强连通子图，都可以互相到达，如果多于一个，说明有一部分点之间无法互相抵达）

如何求解强连通分量

求解强连通分量有三种方法，分别是Tarjan，Kosaraju，Garbow（蒟蒻刚学的tarjan，只能讲讲tarjan了）

tarjan的求解其实是通过求解树的 dfs 序的方法来求解的，需要两个额外数组，dfn[maxn]，low[maxn]

·dfn记录的是某一个点被访问的序号（时间戳）

·low[ now ]记录的是，通过一条不在当前搜索序列里面的边能到达的最小的dfn的值

在求解dfs序列的过程中，每一个点分别入栈，并且对应的标记设置为1，初始化low[ i ] = dfn[ i ]，扫描 i 的邻接边，扫描的过程中判断

1、如果当前的点now没有被访问过，递归访问进去，当返回的时候令low[ i ] = min(low[ now ], low[ i ])；

2、如果当前点now被访问过并且依然存在于栈中，令 low[ i ] = min(low [ i ], dfn[ now ]);

3、在回溯之前，判断是否有 low[ i ] == dfn[ i ]，如果有，则不断的从栈中弹出节点，直到 x 出栈，出栈过程中将出栈的值的标记设置为0，并且强连通分量的个数增加 1 个（因为每一个强连通分量中只有一个 low[ i ] == dfn[ i ]）

实现

void tarjan(int x){
	st[++top] = x;
	v[x] = 1;
	dfn[x] = low[x] = ++num;
	for (int i = head[x]; i; i = nex[i]){
		int y = node[i];
		if(dfn[y] == 0){
			tarjan(y);
			low[x] = min(low[x], low[y]);
		}
		else if(v[y]){
			low[x] = min(low[x], dfn[y]);
		}
	}
	if(dfn[x] == low[x]){
		++cnt;
		while (1){
			int t = st[top--];
			v[t] = 0;
			if(t == x)break;
		}
	}
}

将tarjan应用到本题中去就行了

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cstring>
using namespace std;
const int maxn = 1e5 + 10;
int st[maxn], dfn[maxn], low[maxn], v[maxn], nex[maxn], head[maxn], node[maxn], top, n, m, cnt, num, k;
void pre_work(){
	memset(st, 0, sizeof st);
	memset(dfn, 0, sizeof dfn);
	memset(low, 0, sizeof low);
	memset(v, 0, sizeof v);
	memset(nex, 0, sizeof nex);
	memset(head, 0, sizeof head);
	memset(node, 0, sizeof node);
	top = 0;
	cnt = 0;
	num = 0;
	k = 0;
}
void add(int a, int b){
	++k;
	node[k] = b;
	nex[k] = head[a];
	head[a] = k;
}
void init(){
	for (int i = 1; i <= m; i++){
		int a, b;
		scanf("%d %d", &a, &b);
		add(a, b);
	}
}
void tarjan(int x){
	st[++top] = x;
	v[x] = 1;
	dfn[x] = low[x] = ++num;
	for (int i = head[x]; i; i = nex[i]){
		int y = node[i];
		if(dfn[y] == 0){
			tarjan(y);
			low[x] = min(low[x], low[y]);
		}
		else if(v[y]){
			low[x] = min(low[x], dfn[y]);
		}
	}
	if(dfn[x] == low[x]){
		++cnt;
		while (1){
			int t = st[top--];
			v[t] = 0;
			if(t == x)break;
		}
	}
}
void solve(){
	pre_work();
	init();
	for (int i = 1; i <= n; i++){
		if(dfn[i] == 0)tarjan(i);
	}
	if(cnt == 1)printf("Yes\n");
	else printf("No\n");
}
int main()
{
	while (~scanf("%d %d", &n, &m) && (n + m)){
		solve();
	}
	return 0;
}