luogu P1935 [国家集训队]圈地计划（最小割）

luogu P1935 [国家集训队]圈地计划（最小割）

题目大意

最近房地产商GDOI(Group of Dumbbells Or Idiots)从NOI(Nuts Old Idiots)手中得到了一块开发土地。据了解，这块土地是一块矩形的区域，可以纵横划分为N×M块小区域。GDOI要求将这些区域分为商业区和工业区来开发。根据不同的地形环境，每块小区域建造商业区和工业区能取得不同的经济价值。更具体点，对于第i行第j列的区域，建造商业区将得到Aij收益，建造工业区将得到Bij收益。另外不同的区域连在一起可以得到额外的收益，即如果区域(i,j)相邻（相邻是指两个格子有公共边）有k块（显然k不超过4）类型不同于(i,j)的区域，则这块区域能增加k×Cij收益。经过Tiger.S教授的勘察，收益矩阵A,B,C都已经知道了。你能帮GDOI求出一个收益最大的方案么？

解题思路

两种点黑白染色，根据染色的不同有不同的收益。很容易就想到这个模型，本题中，a，c分别为$A_x$和$B_x$，而e在x->y时为$C_x$而在y->x时则为$C_y$,至于b，d则与a，c对称

跑出最小割再用总和减之即可

AC代码

#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
const int maxn=10005;
const int maxm=2e5+5;
const int inf=0x3f3f3f3f;
struct Edge{
    int to,nxt,cap,flow;
}edge[maxm];
int tol;
int head[maxn];
void init(){
    tol=2;
    memset(head,-1,sizeof(head));
}
void AddEdge(int u,int v,int w,int rw=0){
    edge[tol].to=v;edge[tol].cap=w;edge[tol].flow=0;
    edge[tol].nxt=head[u];head[u]=tol++;
    edge[tol].to=u;edge[tol].cap=rw;edge[tol].flow=0;
    edge[tol].nxt=head[v];head[v]=tol++;
}
int Q[maxn];
int dep[maxn],cur[maxn],sta[maxn];
bool bfs(int s,int t,int n){
    int front=0,tail=0;
    memset(dep,-1,sizeof(dep[0])*(n+1));
    dep[s]=0;
    Q[tail++]=s;
    while(front<tail){
        int u=Q[front++];
        for(int i=head[u];i!=-1;i=edge[i].nxt){
            int v=edge[i].to;
            if(edge[i].cap>edge[i].flow&&dep[v]==-1){
                dep[v]=dep[u]+1;
                if(v==t) return true;
                Q[tail++]=v;
            }
        }
    }
    return false;
}
int dinic(int s,int t,int n){
    int maxflow=0;
    while(bfs(s,t,n)){
        for(int i=0;i<n;i++) cur[i]=head[i];
        int u=s,tail=0;
        while(cur[s]!=-1){
            if(u==t){
                int tp=inf;
                for(int i=tail-1;i>=0;i--)
                {
                    tp=min(tp,edge[sta[i]].cap-edge[sta[i]].flow);
                }
                maxflow+=tp;
                for(int i=tail-1;i>=0;i--){
                    edge[sta[i]].flow+=tp;
                    edge[sta[i]^1].flow-=tp;
                    if(edge[sta[i]].cap-edge[sta[i]].flow==0) tail=i;
                }
                u=edge[sta[tail]^1].to;
            }
            else if(cur[u]!=-1&&edge[cur[u]].cap>edge[cur[u]].flow&&dep[u]+1==dep[edge[cur[u]].to]){
                sta[tail++]=cur[u];
                u=edge[cur[u]].to;
            }
            else{
                while(u!=s&&cur[u]==-1) u=edge[sta[--tail]^1].to;
                cur[u] = edge [cur[u]].nxt;
            }
        }
    }
    return maxflow;
}
int n,m;
inline int Hash(int i,int j){return (i-1)*m+j;}
int main()
{
	init();
	scanf("%d%d",&n,&m);
	int s=0,t=n*m+1;
	int ans=0;
	int v;
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<=m;j++)
		{
			scanf("%d",&v);ans+=v;
			if((i+j)&1) AddEdge(s,Hash(i,j),v);
			else AddEdge(Hash(i,j),t,v);
		}
	}
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<=m;j++)
		{
			scanf("%d",&v);ans+=v;
			if((i+j)&1) AddEdge(Hash(i,j),t,v);
			else AddEdge(s,Hash(i,j),v);
		}
	}
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<=m;j++)
		{
			scanf("%d",&v);
			if(i-1>=1) AddEdge(Hash(i,j),Hash(i-1,j),v),AddEdge(Hash(i-1,j),Hash(i,j),v),ans+=v;
			if(i+1<=n) AddEdge(Hash(i,j),Hash(i+1,j),v),AddEdge(Hash(i+1,j),Hash(i,j),v),ans+=v;
			if(j-1>=1) AddEdge(Hash(i,j),Hash(i,j-1),v),AddEdge(Hash(i,j-1),Hash(i,j),v),ans+=v;
			if(j+1<=m) AddEdge(Hash(i,j),Hash(i,j+1),v),AddEdge(Hash(i,j+1),Hash(i,j),v),ans+=v;
		}
	}
	printf("%d\n",ans-dinic(s,t,t+1));
}

) AddEdge(Hash(i,j),Hash(i,j+1),v),AddEdge(Hash(i,j+1),Hash(i,j),v),ans+=v;
}
}
printf("%d\n",ans-dinic(s,t,t+1));
}