Tarjan‘s Algorithm Practice

Tarjan’s Algorithm——求强连通分量的线性时间的算法

• 如果有两个vertex可以互相到达，那么则称两个顶点强连通

• 若graph中任意两个顶点都可以互相到达，那么则称这个图为一个强连通图
  其中的一个子图称为强连通分量

• 基于对图的dfs,每个强连通分量是搜索树中的一颗子树

• 维护两个数组dfn[N],low[N],dfn[x]表示顶点x的时间戳，low[x]表示x或x的子树能够回溯的栈中最小的次序

模板题1

题目描述
给出一个 n 个点，m 条边的无向图，求图的割点。

输入格式
第一行输入两个正整数 n,m。

下面m 行每行输入两个正整数 x,y表示 x 到 y 有一条边。

输出格式
第一行输出割点个数。

第二行按照节点编号从小到大输出节点，用空格隔开。

输入输出样例

6 7
1 2
1 3
1 4
2 5
3 5
4 5
5 6

1
5

Code
Attention:“index” is a keyword int C++

#include<bits/stdc++.h>
using namespace std;

int n, m;
const int N = 2e4 + 5;
const int M = 1e5 + 5;
struct Edge {
	int end, next;
}edge[2*M];
int head[N], dfn[N], low[N];
bool vis[N];
int cnt;
void Tarjan(int x,int y) {
	dfn[x] = low[x] = ++cnt;
	int child = 0;
	for (int i = head[x]; i; i = edge[i].next) {
		int end = edge[i].end;
		if (!dfn[end]) {
			Tarjan(end, y);
			low[x] = min(low[x], low[end]);
			if (low[end] >= dfn[x] && x != y) {
				vis[x] = true;
			}
			if (x == y)child++;
		}
		/*else*/ low[x] = min(low[x], dfn[end]);
	}
	if (child >= 2 && x == y)vis[x] = true;
}
int main() {
	ios::sync_with_stdio(false);
	cin.tie(0); cout.tie(0);
	cin >> n >> m;
	int t = 0;
	for (int i = 1; i <= m; i++) {
		int u, v;
		cin >> u >> v;
		++t;
		edge[t].end = v;
		edge[t].next = head[u];
		head[u] = t;
		++t;
		edge[t].end = u;
		edge[t].next = head[v];
		head[v] = t;
	}
	for (int i = 1; i <= n; i++) {
		if (!dfn[i]) {
			Tarjan(i, i);
		}
	}
	int ans = 0;
	for (int i = 1; i <= n; i++) {
		if (vis[i])ans++;
	}
	cout << ans << endl;
	for (int i = 1; i <= n; i++) {
		if (vis[i])cout << i << " ";
	}
	return 0;
}

模板题2

题目描述
有n个城市，中间有单向道路连接，消息会沿着道路扩散，现在给出n个城市及其之间的道路，问至少需要在几个城市发布消息才能让这所有n个城市都得到消息。

输入格式
第一行两个整数n,m表示n个城市，m条单向道路。

以下m行，每行两个整数b,e表示有一条从b到e的道路，道路可以重复或存在自环。

输出格式
一行一个整数，表示至少要在几个城市中发布消息。

输入输出样例

5 4
1 2
2 1
2 3
5 1

2

Code

#define _CRT_SECURE_NO_WARNINGS
#include<bits/stdc++.h>
using namespace std;

int n, m;
const int N = 1e5 + 5;
const int M = 5e5 + 5;
struct Edge {
	int end, next;
}edge[M];
int head[N], dfn[N], low[N], Stack[N], belong[N], indeg[N];
bool vis[N];
int cnt, ind, sum;
void Tarjan(int x) {
	dfn[x] = low[x] = ++cnt;
	Stack[++ind] = x;
	vis[x] = true;
	for (int i = head[x]; i; i = edge[i].next) {
		int end = edge[i].end;
		if (!dfn[end]) {
			Tarjan(end);
			low[x] = min(low[x], low[end]);
		}
		else if(vis[end]){
			low[x] = min(low[x], dfn[end]);
		}
	}
	if (dfn[x] == low[x]) {
		++sum;
		do {
			int tmp = Stack[ind];
			--ind;
			vis[tmp] = false;
			belong[tmp] = sum;
		} while (x != Stack[ind + 1]);
	}
}
int main() {
	ios::sync_with_stdio(false);
	cin.tie(0); cout.tie(0);
	cin >> n >> m;
	for (int i = 1; i <= m; i++) {
		int u, v;
		cin >> u >> v;
		if(u == v)continue;//避免成环
		edge[i].end = v;
		edge[i].next = head[u];
		head[u] = i;
	}
	for (int i = 1; i <= n; i++) {
		if (!dfn[i])Tarjan(i);
	}
	for (int i = 1; i <= n; i++) {
		for (int j = head[i]; j; j = edge[j].next) {
			int end = edge[j].end;
			if (belong[end] != belong[i]) {
				indeg[belong[end]]++;
			}
		}
	}
	int ans = 0;
	for (int i = 1; i <= sum; i++) {
		if (!indeg[i])ans++;
	}
	cout << ans << endl;
	return 0;
}