[国家集训队]圈地计划

一、题目

点此看题

二、解法

这道题跟happiness很像，图长这样：

SYDevil想出来的方法%%%。
我们把矩形上的点红黑染色，如果是黑色，就把商业区(0)连$s$，工业区(1)连$t$，如果是红色的话就反过来，然后再在中间建双向边，容量为不同的花费之和。
可以分类讨论这个图的情况，发现答案为$sum-flow$(总和-最小割)，跑一遍最大流即可。

#include <cstdio>
#include <iostream>
#include <cstring>
#include <queue>
#define inf 2e9
const int MAXN = 10005;
using namespace std;
int read()
{
    int num=0,flag=1;char c;
    while((c=getchar())<'0'||c>'9')if(c=='-')flag=-1;
    while(c>='0'&&c<='9')num=(num<<3)+(num<<1)+(c^48),c=getchar();
    return num*flag;
}
int n,m,S,T,tot=1,sum,ans,f[MAXN],cur[MAXN],dis[MAXN];
int dx[4]={1,-1},dy[4]={0,0,1,-1},c[105][105];
queue<int> q;
struct edge
{
	int v,c,next;
}e[MAXN*100];
int cal(int x,int y)
{
	return (x-1)*m+y;
}
void add_edge(int u,int v,int c)
{
	e[++tot]=edge{v,c,f[u]},f[u]=tot;
	e[++tot]=edge{u,0,f[v]},f[v]=tot;
}
int bfs()
{
	memset(dis,0,sizeof dis);
	dis[S]=1;
	q.push(S);
	while(!q.empty())
	{
		int u=q.front();
		q.pop();
		for(int i=f[u];i;i=e[i].next)
		{
			int v=e[i].v;
			if(e[i].c>0 && !dis[v])
			{
				dis[v]=dis[u]+1;
				q.push(v);
			}
		}
	}
	if(!dis[T]) return 0;
	return 1;
}
int dfs(int u,int ept)
{
	if(u==T) return ept;
	int flow=0,tmp=0;
	for(int &i=cur[u];i;i=e[i].next)
	{
		int v=e[i].v;
		if(dis[u]+1==dis[v] && e[i].c>0)
		{
			tmp=dfs(v,min(e[i].c,ept));
			if(!tmp) continue;
			ept-=tmp;
			e[i].c-=tmp;
			e[i^1].c+=tmp;
			flow+=tmp;
			if(!ept) break;
		}
	}
	return flow;
}
int main()
{
	n=read();m=read();
	S=0;T=cal(n,m)+1;
	for(int i=1;i<=n;i++)
		for(int j=1;j<=m;j++)
		{
			int t=read();
			if((i+j)%2)
				add_edge(S,cal(i,j),t);
			else
				add_edge(cal(i,j),T,t);
			sum+=t;
		}
	for(int i=1;i<=n;i++)
		for(int j=1;j<=m;j++)
		{
			int t=read();
			if((i+j)%2==0)
				add_edge(S,cal(i,j),t);
			else
				add_edge(cal(i,j),T,t);
			sum+=t;
		}
	for(int i=1;i<=n;i++)
		for(int j=1;j<=m;j++)
			c[i][j]=read();
	for(int i=1;i<=n;i++)
		for(int j=1;j<=m;j++)
		{
			if(i<n)
			{
				sum+=c[i][j]+c[i+1][j];
				add_edge(cal(i,j),cal(i+1,j),c[i][j]+c[i+1][j]);
				add_edge(cal(i+1,j),cal(i,j),c[i][j]+c[i+1][j]);
			}
			if(j<m)
			{
				sum+=c[i][j]+c[i][j+1];
				add_edge(cal(i,j),cal(i,j+1),c[i][j]+c[i][j+1]);
				add_edge(cal(i,j+1),cal(i,j),c[i][j]+c[i][j+1]);
			}
		}
	while(bfs())
	{
		for(int i=0;i<=T;i++)
			cur[i]=f[i];
		ans+=dfs(S,inf);
	}
	printf("%d\n",sum-ans);
}