[POJ3352]Road Construction

题目

传送门 to POJ

题意概要
求最少加入多少条边，使得图中不存在桥（即，图变为一个边双连通图）。

数据范围与提示
点数 $n$ 与边数 $m$ 均不超过 $10^3$ 。虽然事实上这个限制过于宽松了……

思路

首先求出所有桥。将边双连通分量缩点。图变为一棵树。

考虑此时加入一条边的作用：使得两个端点路径上经过的桥都不再是桥。

于是问题转化为，找最少的路径，使得覆盖所有边至少一次。这个听上去像《欧拉图相关的生成与计数问题探究》里的东西，但是好像又不是。

事实上，我们在洛谷月赛「PMOI-4」可怜的团主里面见过它。可以给出一种 神奇的构造：按照 $\mathrm{d}\mathrm{f}\mathrm{s}$ 序把叶子节点排序得到 $a_1,a_2,\dots ,a_k(2\mid k)$，则可以 $a_i(i⩽\frac{k}{2})$ 与 $a_{i+\frac{k}{2}}$ 之间作为一条路径。此时恰好将每条边、每个点都经过了至少一次。

构造是困难的，证明是容易的。譬如一条边被经过，等价于它的下端点的子树内的某个点走向了子树外的一点。而任意一个子树（非根）的 $\mathrm{d}\mathrm{f}\mathrm{s}$ 区间长度如果超过 $\frac{k}{2}$，那么区间外的叶节点肯定与区间内的叶节点匹配，得证；如果区间长度不超过 $\frac{k}{2}$，那么区间内的叶节点肯定与区间外的叶节点匹配，亦得证。

容易注意到这也是下界。因为 $k$ 个叶节点到其父节点的边，一次最多覆盖 $2$ 条。太离谱了，构造一下子就得到了下界……

那么 $2\nmid k$ 咋办呢？随便补一个节点，反正都能达到下界。但在此之前先判断——当根节点只有一个子节点时，要把根节点算作叶节点，否则可能不被经过。根节点可以作为 $a$ 序列的首项，也可以是末项。它不会影响正确性。

所以这道题就简单极了——统计出度为 $1$ 的边双连通分量即可。时间复杂度 $\mathcal{O}(n+m)$ 。

代码

#include <cstdio>
#include <iostream>
#include <vector>
#include <cstring>
#include <algorithm>
using namespace std;
# define rep(i,a,b) for(int i=(a); i<=(b); ++i)
# define drep(i,a,b) for(int i=(a); i>=(b); --i)
typedef long long int_;
inline int readint(){
	int a = 0, c = getchar(), f = 1;
	for(; '0'>c||c>'9'; c=getchar())
		if(c == '-') f = -f;
	for(; '0'<=c&&c<='9'; c=getchar())
		a = (a<<3)+(a<<1)+(c^48);
	return a*f;
}
void writeUnsigned(const unsigned &x){
	if(x > 9) writeUnsigned(x/10);
	putchar((x-x/10*10)^48);
}
inline void writeint(const int &x){
	if(x < 0){
		putchar('-');
		writeUnsigned(-x);
	}
	else writeUnsigned(x);
}

const int MAXN = 1005;
struct Edge{
	int to, nxt;
	Edge() = default;
	Edge(int _to,int _nxt){
		to = _to, nxt = _nxt;
	}
};
Edge e[MAXN<<1];
int head[MAXN], cntEdge;
void addEdge(int a,int b){
	e[cntEdge] = Edge(b,head[a]);
	head[a] = cntEdge ++;
}

int dfn[MAXN], low[MAXN], dfsClock;
bool cut[MAXN]; /// if it's cut edge
void tarjan(int x,int pre){
	dfn[x] = low[x] = ++ dfsClock;
	for(int i=head[x]; ~i; i=e[i].nxt){
		if(e[i].to == pre) continue;
		if(dfn[e[i].to] == 0){
			tarjan(e[i].to,x);
			low[x] = min(low[x],low[e[i].to]);
			if(low[e[i].to] > dfn[x])
				cut[(i>>1)+1] = true;
		}
		else low[x] = min(low[x],dfn[e[i].to]);
	}
}

bool vis[MAXN];
int dfs(int x){
	vis[x] = true; int res = 0;
	for(int i=head[x]; ~i; i=e[i].nxt)
		if(cut[(i>>1)+1]) ++ res;
		else if(!vis[e[i].to])
			res += dfs(e[i].to);
	return res;
}

int main(){
	int n = readint(), m = readint();
	memset(head+1,-1,n<<2);
	for(int a,b; m; --m){
		a = readint(), b = readint();
		addEdge(a,b), addEdge(b,a);
	}
	tarjan(1,0); int cnt = 0;
	for(int i=1; i<=n; ++i)
		if(!vis[i] && dfs(i) == 1)
			++ cnt; // leaf
	writeint((cnt+1)>>1);
	putchar('\n');
	return 0;
}