P3225 [HNOI2012] 礦場搭建

P3225 [HNOI2012] 礦場搭建

點雙縮點+分討結論

我們可以看看這樣的點滿足什麼條件：去掉這個點以及所連出的邊，圖不連通。顯然是割點，那麼考慮將圖用點雙縮成樹，發掘性質。

縮完是個森林，考慮度數為 \(0\)、\(1\)、\(2\) 分別討論。很明顯的，樹上的葉子結點，即度數為 \(1\) 的點一定要救援出口，否則堵住父親它就會成為獨立點。度數為 \(0\) 的點為獨立點，那麼一定要兩個救援出口。度數為 \(2\) 的可以看出不需要救援出口，因為它所在的連通塊中一定有度數為 \(1\) 的點，並且只能堵住一個出口，一定有另一個出口通向其他邊雙。

那麼計算答案就簡單了，度數為 \(0\) 就是 \(num\choose 2\)，度數為一就是 \(num\)，值得注意的是後者的 \(num\) 不包括點雙中的割點。

找點雙直接遍歷節點即可，因為去到另一個點雙一定要越過一個割點。

複雜度 \(O(n)\)。

#include <iostream>
#include <vector>
#include <cstdio>
#include <cmath>
#include <string>
#include <array>
#include <map>
#include <set>
#include <queue>
#include <stack>
#include <algorithm>
#include <cstdlib>
#include <bitset>
#define pii std::pair<int, int>
#define fi first
#define se second
#define pb push_back

typedef long long i64;
const i64 iinf = 0x3f3f3f3f, linf = 0x3f3f3f3f3f3f3f3f;
const int N = 1e3 + 10;
int n, m, tot, rt, cnt, num, q, idx;
i64 ans1, ans2 = 1;
std::vector<int> e[N];
int dfn[N], low[N], vis[N], cut[N];
void tarjan(int u) {
	dfn[u] = low[u] = ++tot;
	int flg = 0;
	for(auto v : e[u]) {
		if(!dfn[v]) {
			tarjan(v);
			low[u] = std::min(low[u], low[v]);
			if(low[v] >= dfn[u]) {
				flg++;
				if(u != rt || flg > 1) cut[u] = 1;
			}
		} else low[u] = std::min(low[u], dfn[v]);
	}
}
void dfs(int u) {
	vis[u] = idx;
	if(cut[u]) return;
	cnt++;
	vis[u] = 1;
	for(auto v : e[u]) {
		if(cut[v] && vis[v] != idx) num++, vis[v] = idx; 
		if(!vis[v]) dfs(v);
	}
}
void Solve() {
	while(1) {
		q++;
		std::cin >> m;
		if(!m) break;
		for(int i = 1; i <= m; i++) {
			int u, v;
			std::cin >> u >> v;
			e[u].pb(v), e[v].pb(u);
			n = std::max({n, u, v});
		}
		for(int i = 1; i <= n; i++) {
			if(!dfn[i]) rt = i, tarjan(i);
		}
		for(int i = 1; i <= n; i++) {
			if(!cut[i] && !vis[i]) {
				idx++;
				num = cnt = 0;
				dfs(i);
				if(!num) ans1 += 2, ans2 *= cnt * (cnt - 1) / 2;
				if(num == 1) ans1++, ans2 *= cnt; 
			}
		}
		std::cout << "Case " << q << ": " << ans1 << " " << ans2 << "\n";
		for(int i = 1; i <= n; i++) e[i].clear(), dfn[i] = low[i] = cut[i] = vis[i] = 0;
		cnt = num = n = 0;
		ans1 = 0, ans2 = 1;
	}
}
int main() {
    std::ios::sync_with_stdio(false);
    std::cin.tie(nullptr);
    
	Solve();

	return 0;
}