BZOJ 2127 happiness 最小割

标签：bzoj   bzoj2127   最小割   dinic   

题目大意：给定一个座位图，相邻两人之间是朋友，每个人选择学文或学理会有相应的喜悦值，一对朋友同时选择学文/学理也会有相应的喜悦值，求喜悦值之和最大的方案

这个题的模型显然是最小割- -

看到矩阵上相邻点之间的关系 很容易想到黑白染色 随后就能想到将某种颜色的点源汇对调

但是很可惜我没建出图来- - 自己YY了一种做法 但是喜闻乐见地发现建图方法不对- -

还是OTZ HZWER一下吧- - http://hzwer.com/2422.html

我不赘述了- -

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define M 10100
#define S 0
#define T (m*n+1)
#define INF 0x3f3f3f3f
#define P(i,j) ((i)*n-n+(j))
using namespace std;
struct abcd{
	int to,f,next;
}table[1001001];
int head[M],tot=1;
int m,n,ans;
int a[110][110],b[110][100];
void Add(int x,int y,int z)
{
	table[++tot].to=y;
	table[tot].f=z;
	table[tot].next=head[x];
	head[x]=tot;
}
void Link(int x,int y,int z)
{
	Add(x,y,z);
	Add(y,x,z);
}
namespace Max_Flow{
	int dpt[M];
	bool BFS()
	{
		static int q[M];
		int i,r=0,h=0;
		memset(dpt,-1,sizeof dpt);
		dpt[S]=1;q[++r]=S;
		while(r!=h)
		{
			int x=q[++h];
			for(i=head[x];i;i=table[i].next)
				if(table[i].f&&!~dpt[table[i].to])
				{
					dpt[table[i].to]=dpt[x]+1;
					q[++r]=table[i].to;
					if(table[i].to==T)
						return true;
				}
		}
		return false;
	}
	int Dinic(int x,int flow)
	{
		int i,left=flow;
		if(x==T) return flow;
		for(i=head[x];i&&left;i=table[i].next)
			if(table[i].f&&dpt[table[i].to]==dpt[x]+1)
			{
				int temp=Dinic(table[i].to,min(left,table[i].f) );
				left-=temp;
				table[i].f-=temp;
				table[i^1].f+=temp;
			}
		if(left) dpt[x]=-1;
		return flow-left;
	}
}
int main()
{
	int i,j,x;
	cin>>m>>n;
	for(i=1;i<=m;i++)
		for(j=1;j<=n;j++)
		{
			scanf("%d",&a[i][j]);
			ans+=a[i][j];
			a[i][j]<<=1;
		}
	for(i=1;i<=m;i++)
		for(j=1;j<=n;j++)
		{
			scanf("%d",&b[i][j]);
			ans+=b[i][j];
			b[i][j]<<=1;
		}
	for(i=1;i<m;i++)
		for(j=1;j<=n;j++)
		{
			scanf("%d",&x);ans+=x;
			a[i][j]+=x;a[i+1][j]+=x;
			Link(P(i,j),P(i+1,j),x);
			
		}
	for(i=1;i<m;i++)
		for(j=1;j<=n;j++)
		{
			scanf("%d",&x);ans+=x;
			b[i][j]+=x;b[i+1][j]+=x;
			Link(P(i,j),P(i+1,j),x);
			
		}
	for(i=1;i<=m;i++)
		for(j=1;j<n;j++)
		{
			scanf("%d",&x);ans+=x;
			a[i][j]+=x;a[i][j+1]+=x;
			Link(P(i,j),P(i,j+1),x);
		}
	for(i=1;i<=m;i++)
		for(j=1;j<n;j++)
		{
			scanf("%d",&x);ans+=x;
			b[i][j]+=x;b[i][j+1]+=x;
			Link(P(i,j),P(i,j+1),x);
		}
	for(i=1;i<=m;i++)
		for(j=1;j<=n;j++)
			Link(S,P(i,j),a[i][j]),Link(P(i,j),T,b[i][j]);
	using namespace Max_Flow;
	ans<<=1;
	while( BFS() )
		ans-=Dinic(S,INF);
	cout<<(ans>>1)<<endl;
	return 0;
}

BZOJ 2127 happiness 最小割

标签：bzoj   bzoj2127   最小割   dinic   

原文地址：http://blog.csdn.net/popoqqq/article/details/42607893