求无向图的最小点割集?
最近有点神经质,什么都想用数学求解。。。
最优方案显然是割点集合将图分为两部分。
因此考虑为每个点标号,可以划分为S集合、T集合、以及两个都在。
两个集合都在的点即属于割点集。
标号S集合的为0、T集合的为1,但是割点不容表示,因此每个点需要两个变量来表示:
x′i
表示是否在S集合,0表示在;
x″i
表示是否在T集合,1表示在。
对于每条边,显然边的两端点必须同时属于1个集合,即都在S或都在T。
那么答案即为
min{∑max{0,x″k−x′k}⋅1+∑max{0,x′j−x″i}⋅∞∑max{0,x′i−x″j}⋅∞}(<i,j>∈E)
即如果同时属于2个集合,就贡献答案1;每条边两端点必须在同一集合。
最小割的经典形式,但是没有规定原汇点,即S与T集合的划分还是很模糊,可以考虑规定某个点 i 必须属于S集合,即(
就我把这道题的题解写的这么麻烦?
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
const int inf = 0x3f3f3f3f, N = 105, M = 10005;
int level[N], n, cnt = 0, v[M], w[M], p[M], h[N], q[M], s, t;
void add(int a, int b, int c) {
p[++cnt] = h[a]; v[cnt] = b; w[cnt] = c; h[a] = cnt;
p[++cnt] = h[b]; v[cnt] = a; w[cnt] = 0; h[b] = cnt;
}
bool bfs() {
int f = 0, r = 0, u, i;
memset(level, -1, sizeof level);
q[r++] = s; level[s] = 1;
while (f < r) {
u = q[f++];
for (i = h[u]; i; i = p[i]) {
if (w[i] && level[v[i]] == -1) {
level[v[i]] = level[u] + 1;
q[r++] = v[i]; if (v[i] == t) return 1;
}
}
}
return false;
}
int dfs(int u, int low) {
int i, tmp = 0; int res = 0;
if (u == t) return low;
for (i = h[u]; i && res < low; i = p[i]) {
if (w[i] && level[v[i]] == level[u] + 1) {
tmp = dfs(v[i], min((int) w[i], low - res));
w[i] -= tmp; w[i ^ 1] += tmp; res += tmp;
}
}
if (!res) level[u] = -1;
return res;
}
int dinic() {
int ans = 0;
while (bfs()) ans += dfs(s, inf);
return ans;
}
void reset() {
for (int i = 2; i <= cnt; i += 2)
w[i] += w[i ^ 1], w[i ^ 1] = 0;
}
int main() {
int m, a, b, c, np, nc, i, ans;
while (scanf("%d%d", &n, &m) == 2) {
if (n == 0 && m == 0) { puts("0"); continue; }
s = 0 + n; memset(h, 0, sizeof h); cnt = 1; ans = inf;
for (i = 0; i < n; ++i) add(i, i + n, 1);
for (i = 0; i < m; ++i) {
scanf(" (%d,%d)", &a, &b);
add(a + n, b, inf);
add(b + n, a, inf);
}
for (i = 1; i < n; ++i) {
reset(); t = i;
ans = min(ans, dinic());
}
printf("%d\n", ans == inf ? n : ans);
}
return 0;
}