SPOJ839 Optimal Marks（最小割）

题意：
$\ \ \ \ $给出$n$个点$m$条边的无向图，其中一些点已经有标号，定义一条边的权值是该边两个端点的标号异或值，剩余没有标号的点要你去给它标号， 问这个图的所有边权的最小总和是多少。

思路：
$\ \ \ \ $胡伯涛的论文里面的题目， 就直接贴图了。



$\ \ \ \ $然后是一个位一个位的建图， 已知的标号在第$i$位是$1$的时候，和源点连一条容量无穷大的边，是$0$的时候和汇点连一条无穷大的边（代表改边是定值，不可分割），每条边之间容量为$1$建边，代表如果分为两个不同的类别产生的值为 1 1 <script type="math/tex" id="MathJax-Element-503">1</script>

#include<bits/stdc++.h>
typedef long long ll;
const int maxn = 1e3 + 10;
const int INF = 1e9;
using namespace std;

struct st {
    int to, cap, re, id;
    st(int to = 0, int cap = 0, int re = 0, int id = 0) :
        to(to), cap(cap), re(re), id(id) {}
};
int n, m, T, kase = 1, s, t;
int mark[maxn], know[maxn], vis[maxn];
vector<int> g[maxn];
int lv[maxn], it[maxn];
vector<st> G[maxn];

void add(int f, int t, int c, int id) {
    G[f].push_back(st(t, c, G[t].size(), id));
    G[t].push_back(st(f, 0, G[f].size() - 1, -1));
}

void bfs() {
    memset(lv, -1, sizeof(lv));
    queue<int> q; lv[s] = 0; q.push(s);
    while(!q.empty()) {
        int u = q.front(); q.pop();
        for(int i = 0; i < G[u].size(); i++) {
            st &e = G[u][i];
            if(e.cap > 0 && lv[e.to] < 0) {
                lv[e.to] = lv[u] + 1;
                q.push(e.to);
            }
        }
     }
}

int dfs(int v, int t, int f) {
    if(v == t) return f;
    for(int &i = it[v]; i < G[v].size(); i++) {
        st &e = G[v][i];
        if(e.cap > 0 && lv[v] < lv[e.to]) {
            int d = dfs(e.to, t, min(f, e.cap));
            if(d > 0) {
                e.cap -= d;
                G[e.to][e.re].cap += d;
                return d;
            }
        }
    }
    return 0;
}

int maxflow() {
    int f = 0;
    while(1) {
        bfs();
        if(lv[t] < 0) return f;
        memset(it, 0, sizeof(it));
        int fl;
        while((fl = dfs(s, t, INF)) >0) f += fl;
    }
}

void dfs(int x) {
    vis[x] = 1;
    for(int i = 0; i < G[x].size(); i++) {
        int nxt = G[x][i].to;
        if(vis[nxt] || !G[x][i].cap) continue;
        dfs(nxt);
    }
}

void build_graph(int dig) {
    s = 0; t = 2 * n + 1;
    for(int i = 0; i <= t; i++)  G[i].clear();
    for(int i = 1; i <= n; i++) {
        if(know[i]) {
            int flag = (mark[i] >> dig) & 1;
            if(flag) add(s, i, INF, -1);
            else add(i, t, INF, -1);
        }
        for(int j = 0; j < g[i].size(); j++) {
            int to = g[i][j];
            add(i, to, 1, 1);
        }
    }
    maxflow();
    memset(vis, 0, sizeof vis);
    dfs(s);
    for(int i = 1; i <= n; i++) {
        if(know[i]) continue;
        if(vis[i]) mark[i] = mark[i] + (1 << dig);
    }
}


int main() {
    scanf("%d", &T);
    while(T--) {
        for(int i = 0; i < maxn; i++) {
            G[i].clear();
            mark[i] = know[i] = 0;
            g[i].clear();
        }
        scanf("%d %d", &n, &m);
        for(int i = 0; i < m; i++) {
            int u, v; scanf("%d %d", &u, &v);
            g[u].push_back(v);
            g[v].push_back(u);
        }
        scanf("%d", &kase);
        int si;
        while(kase--) {
            scanf("%d", &si);
            scanf("%d", &mark[si]);
            know[si] = 1;
        }
        for(int i = 0; i < 31; i++) build_graph(i);
        for(int i = 1; i <= n; i++) cout << mark[i] << endl;
    }
    return 0;
}