题目地址:
https://www.luogu.com.cn/problem/P3388
题目背景:
割点
题目描述:
给出一个
n
n
n个点,
m
m
m条边的无向图,求图的割点。
输入格式:
第一行输入两个正整数
n
,
m
n,m
n,m。下面
m
m
m行每行输入两个正整数
x
,
y
x,y
x,y表示
x
x
x到
y
y
y有一条边。
输出格式:
第一行输出割点个数。第二行按照节点编号从小到大输出节点,用空格隔开。对于全部数据,
1
≤
n
≤
2
×
1
0
4
1\leq n \le 2\times 10^4
1≤n≤2×104,
1
≤
m
≤
1
×
1
0
5
1\leq m \le 1 \times 10^5
1≤m≤1×105。点的编号均大于
0
0
0小于等于
n
n
n。
割点指的是这样的一种点,将其删去后整个图的连通块个数会增加。可以用Tarjan算法。与求桥类似,可以从任意一个点
u
u
u开始做DFS,接着开两个数组dfn
和low
,dfn[u]
代表第一次访问
u
u
u时的时间戳,而low[u]
代表从
u
u
u不经由其DFS树的祖宗节点而能走到的时间戳最小的点的时间戳(其祖宗节点可以计入,但是如果经由祖宗节点继续向上能走到时间戳更小的点,这些点就不能计入了。也就是说,这里的low[u]
表示从
u
u
u能“绕上去”的时间戳最小的节点的时间戳)。首先如果
u
u
u是DFS树根,且其有两个孩子互相不可达(经过
u
u
u可达不算),那么
u
u
u显然是割点,因为去掉其之后就会产生至少两个连通块;如果
u
u
u不是DFS树根,那么要看一下
u
u
u是否有个孩子
v
v
v满足low[v] >= dfn[u]
,如果成立,则说明从
v
v
v无法不沿着DFS树而绕到
u
u
u或
u
u
u的上方,则说明
u
u
u是割点。注意这里不能写成low[v] >= low[u]
,否则有可能
v
v
v是经过
u
u
u再向上走到时间戳更小的点的,从而有可能即使
u
u
u是割点但是漏解了。同样,必须写if (v != from) low[u] = min(low[u], dfn[v]);
,不能用low[v]
来更新,因为有可能是通过祖宗向上走的。看下图,
1
→
2
→
3
→
4
→
5
1\to 2 \to 3\to 4\to 5
1→2→3→4→5,如果用low[v]
来更新low[u]
的话,则low[5] = 1
,从而low[4] = low[3] = 1
,回溯到
3
3
3的时候,发现low[4] < dfn[3]
从而漏掉
3
3
3这个割点。
代码如下:
#include <iostream>
#include <cstring>
using namespace std;
// 边开两倍
const int N = 2e4 + 10, M = 2e5 + 10;
int n, m;
int h[N], e[M], ne[M], idx;
int dfn[N], low[N], timestamp;
bool cut[N];
int tot;
void add(int a, int b) {
e[idx] = b, ne[idx] = h[a], h[a] = idx++;
}
void tarjan(int u, int from) {
dfn[u] = low[u] = ++timestamp;
int child = 0;
for (int i = h[u]; ~i; i = ne[i]) {
int v = e[i];
if (!dfn[v]) {
tarjan(v, u);
low[u] = min(low[u], low[v]);
// 累加能割出的子树数
if (low[v] >= dfn[u]) child++;
} else if (v != from)
// 这里需要注意,要用dfn[v]来更新low[u],
// 防止v能经过DFS树上的某个祖宗走到了时间戳更小的点,从而漏解
low[u] = min(low[u], dfn[v]);
}
// 如果u是树根,且能割出两个子树,则u是割点;若u不是树根,则只需割出一个子树就是割点
if (u == from && child >= 2 || u != from && child) {
cut[u] = true;
tot++;
}
}
int main() {
scanf("%d%d", &n, &m);
memset(h, -1, sizeof h);
while (m--) {
int a, b;
scanf("%d%d", &a, &b);
add(a, b), add(b, a);
}
for (int i = 1; i <= n; i++)
if (!dfn[i])
tarjan(i, i);
printf("%d\n", tot);
for (int i = 1; i <= n; i++)
if (cut[i])
printf("%d ", i);
}
时间复杂度 O ( n + m ) O(n+m) O(n+m),空间 O ( n ) O(n) O(n)。