题目地址:
https://www.acwing.com/problem/content/397/
为了从 F F F个草场中的一个走到另一个,奶牛们有时不得不路过一些她们讨厌的可怕的树。奶牛们已经厌倦了被迫走某一条路,所以她们想建一些新路,使每一对草场之间都会至少有两条相互分离的路径,这样她们就有多一些选择。每对草场之间已经有至少一条路径。给出所有 R R R条双向路的描述,每条路连接了两个不同的草场,请计算最少的新建道路的数量,路径由若干道路首尾相连而成。两条路径相互分离,是指两条路径没有一条重合的道路。但是,两条分离的路径上可以有一些相同的草场。对于同一对草场之间,可能已经有两条不同的道路,你也可以在它们之间再建一条道路,作为另一条不同的道路。
输入格式:
第
1
1
1行输入
F
F
F和
R
R
R。接下来
R
R
R行,每行输入两个整数,表示两个草场,它们之间有一条道路。
输出格式:
输出一个整数,表示最少的需要新建的道路数。
数据范围:
1
≤
F
≤
5000
1≤F≤5000
1≤F≤5000
F
−
1
≤
R
≤
10000
F−1≤R≤10000
F−1≤R≤10000
首先介绍边双连通分量的概念。如果一个连通块里面没有桥,则称这个连通块是个边双连通分量。桥是指的这样的一种边,删去它之后,连通块个数会变化(可以证明个数如果变化了,那一定是增了 1 1 1,而不会是其它的变化。考虑桥所在的连通块,如果删掉它,那么该连通块不会凭空消失,所以个数不可能是减少,因为减少了的话连通块个数就变成 0 0 0了,这是不可能的。如果是增加了 2 2 2,那么总共就是 3 3 3个连通块,而一条边不可能将三个连通块同时连接起来的)。所以桥其实就是删掉它,连通块个数会增 1 1 1的那种边。那么原题其实就是问,一个连通块至少增加多少条边,就能变成一个双连通分量。回忆强连通分量的概念,一个有向连通图求完强连通分量后缩点之后,会得到一个拓扑图;而一个无向连通图求完强连通分量后缩点之后,则会得到一棵树,每条边是原图的桥。具体结论是,要加的边数恰好等于树的叶子数 c c c加上 1 1 1后除以 2 2 2下取整,即 ⌊ c + 1 2 ⌋ \lfloor\frac{c+1}{2}\rfloor ⌊2c+1⌋(证明略繁琐,大致意思是可以将叶子配对连起来)。而求无向图的双连通分量可以用Tarjan算法,和求有向图的强连通分量差不多,只不过要注意这是无向图,DFS的时候不能走回头路。代码如下:
#include <iostream>
#include <cstring>
using namespace std;
const int N = 5000, M = 20010;
int n, m;
int h[N], e[M], ne[M], idx;
int dfn[N], low[N], timestamp;
int stk[N], top;
int id[N], dcc_cnt;
bool is_bridge[M];
int d[N];
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;
stk[top++] = u;
for (int i = h[u]; ~i; i = ne[i]) {
int j = e[i];
if (!dfn[j]) {
tarjan(j, i);
low[u] = min(low[u], low[j]);
if (dfn[u] < low[j])
is_bridge[i] = is_bridge[i ^ 1] = true;
// 不走回头路,from ^ 1是from的反向边
} else if (i != (from ^ 1)) low[u] = min(low[u], dfn[j]);
}
if (dfn[u] == low[u]) {
dcc_cnt++;
int y;
do {
y = stk[--top];
id[y] = dcc_cnt;
} while (y != u);
}
}
int main() {
cin >> n >> m;
memset(h, -1, sizeof h);
while (m--) {
int a, b;
cin >> a >> b;
add(a, b), add(b, a);
}
tarjan(1, -1);
// 缩点之后,每两个点之间就是被桥连起来的,所以累加度数就是在把桥的两个端点累加度数
for (int i = 0; i < idx; i++)
if (is_bridge[i])
d[id[e[i]]]++;
// 数叶子
int cnt = 0;
for (int i = 1; i <= dcc_cnt; i++)
if (d[i] == 1)
cnt++;
cout << (cnt + 1) / 2 << endl;
return 0;
}
时间复杂度 O ( n + m ) O(n+m) O(n+m),空间 O ( n ) O(n) O(n), n n n和 m m m分别是图的点数和边数。