SPOJ Optimal Marks 最小割

题目链接：https://vjudge.net/problem/SPOJ-OPTM

题意：给定一个无向图，每个顶点都有对应值，每条边的权值为val[u] ^ val[u]，同时已知其中某些点的值，求其他点的权值使得最终的图的边权和最小？

解法：首先将mark拆成31位来看待，对每一位都进行一次求解。于是问题就变成了顶点的标号可以为1也可以为0，边权和最小的问题，于是按标号为1或者为0就可以把顶点划分成两个集合，而连通这两个集合的边的数量就是我们需要累加的边权，为了让数量最少，于是就变成了最小割问题。 这样我们将S和所有mark为1的点连起来，容量为INF，将所有mark为0的点和T连起来，容量为INF，至于那些还没标号的点就不和源点汇点连了，他们的标号取决于做完最小割之后属于哪个集合。接着将其他边建好，容量为1，之后就可以做最大流了，做完最大流后从S出发沿容量不为0的边走所能到达的点都是属于mark为1的这个集合的。

#include <bits/stdc++.h>
using namespace std;
typedef long long LL;
const int maxn = 510;
const int maxm = 100010;
const int inf = 0x3f3f3f3f;
struct G
{
    int v, cap, next;
    G() {}
    G(int v, int cap, int next) : v(v), cap(cap), next(next) {}
} E[maxm];
struct Edge{
    int x, y;
}edge[maxm];
struct List{
    int id, mark;
}lable[maxm];
int p[maxn], T;
int n, m, k, source, sink, mark[maxn], vis[maxn];
int d[maxn], temp_p[maxn], qw[maxn]; //d顶点到源点的距离标号,temp_p当前狐优化，qw队列
void init()
{
    memset(p, -1, sizeof(p));
    T = 0;
}
void add(int u, int v, int cap)
{
    E[T] = G(v, cap, p[u]);
    p[u] = T++;
    E[T] = G(u, 0, p[v]);
    p[v] = T++;
}
bool bfs(int st, int en, int nn)
{
    int i, u, v, head, tail;
    for(i = 0; i <= nn; i++) d[i] = -1;
    head = tail = 0;
    d[st] = 0;
    qw[tail] = st;
    while(head <= tail)
    {
        u = qw[head++];
        for(i = p[u]; i + 1; i = E[i].next)
        {
            v = E[i].v;
            if(d[v] == -1 && E[i].cap > 0)
            {
                d[v] = d[u] + 1;
                qw[++tail] = v;
            }
        }
    }
    return (d[en] != -1);
}
int dfs(int u, int en, int f)
{
    if(u == en || f == 0) return f;
    int flow = 0, temp;
    for(; temp_p[u] + 1; temp_p[u] = E[temp_p[u]].next)
    {
        G& e = E[temp_p[u]];
        if(d[u] + 1 == d[e.v])
        {
            temp = dfs(e.v, en, min(f, e.cap));
            if(temp > 0)
            {
                e.cap -= temp;
                E[temp_p[u] ^ 1].cap += temp;
                flow += temp;
                f -= temp;
                if(f == 0)  break;
            }
        }
    }
    return flow;
}
int dinic(int st, int en, int nn)
{
    int i, ans = 0;
    while(bfs(st, en, nn))
    {
        for(i = 0; i <= n; i++) temp_p[i] = p[i];
        ans += dfs(st, en, inf);
    }
    return ans;
}
void build(int cur){
    init();
    for(int i=1; i<=m; i++) add(edge[i].x,edge[i].y,1), add(edge[i].y,edge[i].x,1);
    for(int i=1; i<=k; i++){
        if((lable[i].mark>>cur)&1) add(source, lable[i].id, inf);
        else add(lable[i].id, sink, inf);
    }
}
void DFS(int cur, int k){
    vis[cur] = 1;
    mark[cur] |= (1<<k);
    for(int i = p[cur]; ~i; i=E[i].next){
        if(!vis[E[i].v]&&E[i].cap)
            DFS(E[i].v, k);
    }
}
int main(){
    int T;
    scanf("%d", &T);
    while(T--){
        scanf("%d %d", &n,&m);
        for(int i=1; i<=m; i++) scanf("%d %d", &edge[i].x,&edge[i].y);
        scanf("%d", &k);
        for(int i=1; i<=k; i++) scanf("%d %d", &lable[i].id, &lable[i].mark);
        source = 0, sink = n+1;
        LL ans = 0;
        memset(mark, 0, sizeof(mark));
        for(int i=0; i<31; i++){
            build(i);
            ans += dinic(source, sink, sink+1) * (1LL<<i);
            memset(vis, 0, sizeof(vis));
            DFS(source, i);
        }
        for(int i=1; i<=n; i++) printf("%d\n", mark[i]);
    }
}