【bzoj2400】Spoj 839 Optimal Marks 最小割

定义无向图中的一条边的值为：这条边连接的两个点的值的异或值。
定义一个无向图的值为：这个无向图所有边的值的和。
给你一个有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);
}

$\mathrm{s}\mathrm{a}\mathrm{p}\mathrm{v}\mathrm{i}\mathrm{s}\mathrm{i}\mathrm{o}\mathrm{n}\mathrm{o}\mathrm{n}\mathrm{s}\mathrm{p}\mathrm{o}\mathrm{j}(\mathrm{i}\mathrm{n}\mathrm{p}\mathrm{u}\mathrm{t}\mathrm{s}\mathrm{v}\mathrm{a}\mathrm{r}\mathrm{y})$
sap真的短。
sap的当前弧优化原来一直都打错了，2hr。

#include<bits/stdc++.h>
#define maxn 505
#define maxm 8005
#define B 1000
#define inf 10000000
#define Clear(h,x) memset(h,x,sizeof h)
#define Copy(x,y) memcpy(x,y,sizeof x)
using namespace std;

int n,m,q,a[maxn],b[maxn],id[maxn],x[maxm],y[maxm];
int info[maxn],Prev[maxm],to[maxm],cap[maxm],cnt_e;
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 h[maxn],g[maxn],buf[maxn],S,T,vis[maxn];

int aug(int u,int mx){
	if(u==T) return mx;
	int st=mx,inc;
	for(int &i=info[u],v;i;i=Prev[i])
		if(cap[i] && h[v=to[i]]+1==h[u]){
			inc = aug(v,min(st,cap[i]));
			st -= inc , cap[i] -= inc , cap[i^1] += inc;
			if(!st){ info[u]=buf[u]; return mx - st; }
		}
	if(!--g[h[u]]) h[S]=T+1;
	++g[++h[u]];info[u]=buf[u];
	return mx - st;
}

void dfs(int u){
	if(vis[u]) return;vis[u]=1;
	for(int i=info[u];i;i=Prev[i]) if(cap[i]) dfs(to[i]);
}

int main(){
	int cas;
	for(scanf("%d",&cas);cas--;){
		scanf("%d%d",&n,&m);S=n+1,T=n+2;
		for(int i=1;i<=m;i++) scanf("%d%d",&x[i],&y[i]);
		scanf("%d",&q);
		Clear(a,-1),Clear(b,0);
		for(int i=1,u,p;i<=q;i++) scanf("%d%d",&u,&p),a[u]=b[u]=p;
		for(int i=30;i>=0;i--){
			cnt_e=1,Clear(info,0);
			for(int j=1;j<=m;j++) Line(x[j],y[j],B,B);
			for(int j=1;j<=n;j++)
				if(a[j]==-1) Line(j,T,1);
				else if(a[j]>>i&1) Line(S,j,inf);
				else Line(j,T,inf);
			Clear(h,0),Clear(g,0),Copy(buf,info);
			int stm = 0;
			for(;h[S]<=T;) stm+=aug(S,inf);
			Clear(vis,0);
			dfs(S);
			for(int j=1;j<=n;j++)
				if(a[j]==-1) b[j]|=(vis[j]<<i);
		}
		for(int i=1;i<=n;i++) printf("%d\n",b[i]);
	}
}