Spoj 839 Optimal Marks（网络流最小割）

2400：Spoj 839 Optimal Marks

时间限制: 1000 ms
内存限制: 131072 KB

题目描述

定义无向图中的一条边的值为：这条边连接的两个点的值的异或值。

定义一个无向图的值为：这个无向图所有边的值的和。

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

输入

第一行，两个数$n$，$m$，表示图的点数和边数。

接下来$n$行，每行一个数，按编号给出每个点的值（若为负数则表示这个点的值由你决定，值的绝对值大小不超过$10^9$）。

接下来$m$行，每行二个数$a$，$b$，表示编号为$a$与$b$的两点间连一条边。（保证无重边与自环。）

输出

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

输入样例

3 2
2
-1
0
1 2
2 3

输出样例

2
2

提示

$n≤500，m≤2000$

样例解释

2结点的值定为0即可。

解：

介绍一个网络流最小割建图方法：一个集合，每个元素可以有两种选择，一些元素选相同选择会付出代价。
具体方法：把两个选择建成两个点，必须选某个的元素向需要的选择连inf,相关联的点连他们的代价。

那么这道题就可以做了。首先异或是可以拆位的，我们单考虑一位，只有0,1两种选择，相邻的点权值相同付出1的代价。至于点权最小，即1最少。我们把选择1的原点跑一遍dfs，连通的必须选1，否则可以选0。这样就是最小的。

code：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<queue>
using namespace std;
struct lxy{
    int to,next,flow;
}eg[5005];
struct egggg{
    int x,y;
}cold[2005];

int head[505],cnt,layer[505],s,t,n,m,data[505],k;
bool vis[505];
long long ans1,ans2;

void add(int op,int ed,int len){
    //cout<<op<<" "<<ed<<" "<<len<<endl;
    eg[++cnt].to=ed;
    eg[cnt].next=head[op];
    head[op]=cnt;
    eg[cnt].flow=len;
}

bool bfs(){
    memset(layer,0,sizeof(layer));
    queue <int> d;
    d.push(s),layer[s]=1;
    while(!d.empty()){
        int now=d.front();d.pop();
        for(int i=head[now];i!=-1;i=eg[i].next)
          if(eg[i].flow!=0&&layer[eg[i].to]==0){
            d.push(eg[i].to);layer[eg[i].to]=layer[now]+1;
          }
    }//cout<<layer[1]<<endl;
    return layer[t];
}

int dfs(int u,int a){
    if(u==t||a==0) return a;
    int flow=0,f;
    for(int i=head[u];i!=-1;i=eg[i].next)
      if(layer[u]+1==layer[eg[i].to]&&eg[i].flow!=0){
        f=dfs(eg[i].to,min(a,eg[i].flow));
        eg[i].flow-=f;
        eg[i^1].flow+=f;
        a-=f;
        flow+=f;
      }
    return flow;
}

int dinic(){
    int ret=0;
    while(bfs()){
        ret+=dfs(s,0x3f3f3f3f);
    }
    return ret;
}

void findit(int u){
    vis[u]=1;if(u!=s) k++;
    for(int i=head[u];i!=-1;i=eg[i].next)
      if(vis[eg[i].to]==0&&eg[i].flow!=0)
        findit(eg[i].to);
}

int main()
{
    memset(head,-1,sizeof(head));
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++)
      scanf("%d",&data[i]);
    for(int i=1;i<=m;i++)
      scanf("%d%d",&cold[i].x,&cold[i].y);
    for(int p=0;p<=30;p++){
        s=0,t=n+1;cnt=-1;
        memset(eg,0,sizeof(eg));
        memset(head,-1,sizeof(head));
        for(int i=1;i<=m;i++)
          add(cold[i].x,cold[i].y,1),add(cold[i].y,cold[i].x,1);
        for(int i=1;i<=n;i++)
          if(data[i]>=0){
            if((data[i]>>p)&1==1) add(s,i,0x3f3f3f3f),add(i,s,0);
            else add(i,t,0x3f3f3f3f),add(t,i,0);
          }
        ans1=ans1+(1ll*dinic()<<p);
        memset(vis,0,sizeof(vis)),k=0;
        findit(s);
        ans2=ans2+(1ll*k<<p);
    }
    printf("%lld\n%lld",ans1,ans2);
}