BZOJ2132 圈地计划

传送~

来自y_immortal学长的网络流推荐2333

建模比较神仙qwq

首先可以建立最小割模型 ->二选一

就是我们处理不同的贡献不好处理 那么我们考虑怎么把它变成相同的贡献

对于原图我们进行黑白染色 然后如果一个点属于T集我们让它反色就可以做到相同的有贡献

那么很明显这个贡献是双向的 对于两个点之间的贡献是相加的

所以我们这样做最小割就可以辣~

代码。

//Love and Freedom.
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<queue>
#define inf 20021225
#define ll long long
#define M 200000
#define N 30000
using namespace std;

struct edge{int to,lt,f;}e[M];
int in[N],cnt=1,dep[N];
int s,t;
void add(int x,int y,int f)
{
	e[++cnt].to = y; e[cnt].lt = in[x]; e[cnt].f = f; in[x] = cnt;
	e[++cnt].to = x; e[cnt].lt = in[y]; e[cnt].f = 0; in[y] = cnt;
}
queue<int> q;
bool bfs()
{
	while(!q.empty())	q.pop();
	memset(dep,0,sizeof(dep));
	dep[s] = 1; q.push(s);
	while(!q.empty())
	{
		int x = q.front(); q.pop();
		for(int i=in[x];i;i=e[i].lt)
		{
			int y = e[i].to;
			if(!dep[y] && e[i].f)
			{
				dep[y] = dep[x]+1;	q.push(y);
				if(dep[t])	return true;
			}
		}
	}
	return false;
}

int dfs(int x,int f)
{
	if(x==t||!f)	return f;
	int cur = f;
	for(int i=in[x];i;i=e[i].lt)
	{
		int y = e[i].to;
		if(dep[y] == dep[x]+1 && e[i].f)
		{
			int flow = dfs(y,min(cur,e[i].f));
			e[i].f -= flow; e[i^1].f += flow;
			cur -= flow; if(cur==0)	return f;
		}
	}
	dep[x] = -1;
	return f-cur;
}

int dinic()
{
	int ans = 0;
	while(bfs())	ans+=dfs(s,inf);
	return ans;
}
int n,m;
int id(int x,int y)
{
	return (x-1)*m + y;
}
int a[101][101];
int b[101][101];
int c[101][101];
int xx[4] = {0,0,1,-1};
int yy[4] = {1,-1,0,0};
int main()
{
	int fin = 0;
	scanf("%d%d",&n,&m);
	s = n*m*2; t = s+1;
	for(int i=1;i<=n;i++)
		for(int j=1;j<=m;j++)
		{
			scanf("%d",&a[i][j]);
			if((i+j)&1)	add(s,id(i,j),a[i][j]);
			else	add(id(i,j),t,a[i][j]);
			fin+=a[i][j];
		}
	for(int i=1;i<=n;i++)
		for(int j=1;j<=m;j++)
		{
			scanf("%d",&b[i][j]);
			if((i+j)&1)	add(id(i,j),t,b[i][j]);
			else	add(s,id(i,j),b[i][j]);
			fin+=b[i][j];
		}
	for(int i=1;i<=n;i++)
		for(int j=1;j<=m;j++)
			scanf("%d",&c[i][j]);
	for(int i=1;i<=n;i++)
		for(int j=1;j<=m;j++)
			for(int k=0;k<4;k++)
			{
				int x = i+xx[k], y = j+yy[k];
				if(x<1||y<1||x>n||y>m)	continue;
				int val = c[i][j] + c[x][y];
				add(id(i,j),id(x,y),val);
				fin+=c[i][j];
			}
	printf("%d\n",fin-dinic());
	return 0;
}