POJ 1966 Cable TV Network 最小割

求无向图的最小点割集？

最近有点神经质，什么都想用数学求解。。。

最优方案显然是割点集合将图分为两部分。
因此考虑为每个点标号，可以划分为S集合、T集合、以及两个都在。
两个集合都在的点即属于割点集。
标号S集合的为0、T集合的为1，但是割点不容表示，因此每个点需要两个变量来表示：$x_i'$表示是否在S集合，0表示在；$x_i''$表示是否在T集合，1表示在。
对于每条边，显然边的两端点必须同时属于1个集合，即都在S或都在T。
那么答案即为

$$\min \{\sum \max\{0,x_k''-x_k'\}\cdot1+\sum \max \{0,x_j'-x_i''\}\cdot \infty \sum \max \{0,x_i'-x_j''\}\cdot \infty \} (<i,j>\in E)$$
即如果同时属于2个集合，就贡献答案1；每条边两端点必须在同一集合。
最小割的经典形式，但是没有规定原汇点，即S与T集合的划分还是很模糊，可以考虑规定某个点 $i$必须属于S集合，即（$x_i'=x_i''=0$）、另外一个点 $j$必须属于T集合（$x_j'=x_j''=1$），答案就很显然了。

就我把这道题的题解写的这么麻烦？

#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;
}