【bzoj2400】Spoj 839 Optimal Marks 最小割dinic

Description

定义无向图中的一条边的值为：这条边连接的两个点的值的异或值。
定义一个无向图的值为：这个无向图所有边的值的和。
给你一个有n个结点m条边的无向图。其中的一些点的值是给定的，而其余的点的值由你决定（但要求均为非负数），使得这个无向图的值最小。在无向图的值最小的前提下，使得无向图中所有点的值的和最小。

Input

第一行，两个数n，m，表示图的点数和边数。
接下来n行，每行一个数，按编号给出每个点的值（若为负数则表示这个点的值由你决定，值的绝对值大小不超过10^9）。
接下来m行，每行二个数a，b，表示编号为a与b的两点间连一条边。（保证无重边与自环。）

Output

第一行，一个数，表示无向图的值。
第二行，一个数，表示无向图中所有点的值的和。

Sample Input

3 2

2

-1

0

1 2

2 3

Sample Output

2

2

HINT

数据约定

n<=500，m<=2000

样例解释

2结点的值定为0即可。

题解
http://blog.csdn.net/u012288458/article/details/50720724
Orz Oxer
Orz YZH

代码

#include<bits/stdc++.h>
#define ll long long
#define inf 1000000000
#define mod 1000000007
using namespace std;
inline int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if (ch=='-') f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
queue<int>q;
int Head[505],ret[20005],Next[20005],len[20005],tot,a[505];
int S,T,n,m,dis[505],cur[505];
ll ans1,ans2;
struct node{int x,y;}e[2005];
inline void ins(int u,int v,int l)
{
    ret[++tot]=v;len[tot]=l;
    Next[tot]=Head[u];Head[u]=tot;
}
void build(int k)
{
    tot=1;for (int i=1;i<=T;i++) Head[i]=0;
    for (int i=1;i<=n;i++)
        if (a[i]>=0)
        {
            if (a[i]&(1<<k))
            {
                ins(S,i,1);ins(i,S,0);
                ins(i,T,inf);ins(T,i,0);
            }
            else ins(S,i,inf),ins(i,S,0);
        }
        else ins(S,i,1),ins(i,S,0);
    for (int i=1;i<=m;i++)
    {
        ins(e[i].x,e[i].y,10000);ins(e[i].y,e[i].x,0);
        ins(e[i].y,e[i].x,10000);ins(e[i].x,e[i].y,0);
    }
}
bool bfs()
{
    for (int i=1;i<=T;i++) dis[i]=-1;q.push(S);dis[S]=0;
    while (!q.empty())
    {
        int now=q.front();q.pop();
        for (int i=Head[now];i;i=Next[i])
        {
            if (dis[ret[i]]==-1&&len[i])
            {
                dis[ret[i]]=dis[now]+1;
                q.push(ret[i]);
            }
        }
    }
    if (dis[T]==-1) return 0;return 1;
}
int dfs(int u,int aug)
{
    if (u==T||aug==0) return aug;
    int flow=0;
    for (int i=cur[u];i;i=Next[i])if (dis[ret[i]]==dis[u]+1)
    {
        int t=dfs(ret[i],min(aug-flow,len[i]));
        len[i]-=t;len[i^1]+=t;flow+=t;
        if (len[i])cur[u]=i;
        if (flow==aug) return flow;
    }
    if (flow==0) dis[u]=-1;
    return flow;
}
int dinic(int k)
{
    int ans=0;
    while (bfs())
    {
        for (int i=1;i<=T;i++) cur[i]=Head[i];
        ans+=dfs(S,inf);
    }
    return ans;
}
int main()
{
    n=read();m=read();S=n+1;T=S+1;
    for (int i=1;i<=n;i++) a[i]=read();
    for (int i=1;i<=m;i++){e[i].x=read();e[i].y=read();}
    for (int k=0;k<=30;k++)
    {
        build(k);
        int ans=dinic(k);
        ans1+=1LL*(ans/10000)*(1<<k);
        ans2+=1LL*(ans%10000)*(1<<k);
    }
    cout<<ans1<<"\n"<<ans2;
    return 0;
}