【bzoj2400】Spoj 839 Optimal Marks 最小割

版权声明：本文为博主原创文章，未经博主允许必须转载。 https://blog.csdn.net/qq_35950004/article/details/88429679

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

麻烦的是第二问，其实也不烦，把第一问的权值乘10000再在%10000的意义下找点最小值。
也可以从汇点dfs，不走反向弧为0的边，找到的点就是一定在T集，反之则在S集，这是二分图最小割的求法。

AC Code:

#include<bits/stdc++.h>
#define maxn 505
#define maxm 30005
#define LL long long
#define inf 0x3f3f3f3f
using namespace std;

int n,m;
int buf[maxn],info[maxn],Prev[maxm],to[maxm],cap[maxm],cnt_e=1;
void Node(int u,int v,int c){ Prev[++cnt_e]=info[u],info[u]=cnt_e,to[cnt_e]=v,cap[cnt_e]=c; }
void Line(int u,int v,int c,int d=0){ Node(u,v,c),Node(v,u,d); }
int a[maxn],x[maxm],y[maxm],dis[maxn],S,T,tot;

int aug(int now,int Max)
{
	if(now == T) return Max;
	int inc , st = Max;
	for(int &i=info[now];i;i=Prev[i])
		if(cap[i] && dis[to[i]] == dis[now] - 1)
		{
			inc = aug(to[i] , min(cap[i] , st));
			if(inc) st-=inc,cap[i]-=inc,cap[i^1]+=inc;
			else dis[to[i]] = -1;
			if(!st) break;
		}
	return Max - st;
}

bool BFS()
{
	static queue<int>q;
	memset(dis,-1,sizeof dis);
	q.push(T),dis[T] = 0;
	for(int now;!q.empty();)
	{
		now = q.front() , q.pop();
		for(int i=info[now];i;i=Prev[i])
			if(cap[i^1] && dis[to[i]]==-1)
			{
				dis[to[i]] = dis[now] + 1;
				q.push(to[i]);
			}
	}
	return dis[S] != -1;
}

int main()
{
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;i++) scanf("%d",&a[i]);
	for(int i=1;i<=m;i++)
	{
		scanf("%d%d",&x[i],&y[i]);
	}
	S = n + 1 , T = n + 2;
	LL ans = 0 , ans2 = 0;
	for(int bit=0;bit<30;bit++)
	{
		memset(info,0,sizeof info),cnt_e=1;
		for(int i=1;i<=m;i++)
			Line(x[i],y[i],10000,10000);
		for(int i=1;i<=n;i++)
			if(a[i]>=0)
			{
				if(a[i]>>bit&1)
					Line(S,i,inf),
					Line(i,T,1);
				else 
					Line(i,T,inf);
			}
			else
			{
				Line(i,T,1);
			}
		int sum = 0;
		memcpy(buf,info,sizeof info);
		for(;BFS();)
			sum += aug(S,inf),
			memcpy(info,buf,sizeof buf);
		ans += (1ll<<bit) * (sum / 10000);//!!!
		ans2 += sum % 10000 * (1ll<<bit);
	}
	printf("%lld\n%lld\n",ans,ans2);
}