[网络流24题] 洛谷P2744 方格取数问题 最小割

给出一个$m\times n$的方格，每个方格都有一个权值，选出若干个方格并且使得权值之和最大但不能有任意两个方格是相邻的，求最大权值。
不太复杂的最小割模型题，但我对这种类型目前还是不太熟悉。
尝试把方格分成两部分，这个可以黑白染色，两种颜色都不会彼此相邻。然后对于黑格，源点连每个格子的权值的流量，对于白格子，连汇点每个格子的权值的流量。对于黑白之间只要在相邻的格子连边，这个流量是无穷大，因为相邻的关系是固定的。
然后跑最小割，答案是全部选上减去最小割。

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int inf=0x3f3f3f3f;
const int N=1e2+7;
int a[N][N]; 
int mp[N][N]; 
namespace Dinic {
	const int N=1e3+3,M=1e4+3; 
	int front[N],head[N],dep[N],tot=-1,n,S,T;
	struct edge{ int v,nxt;ll w; }e[2*M];
	void init() {
		tot=-1;
		for(int i=1;i<=n;i++) front[i]=-1;
	}
	void add(int u,int v,ll w) {
		e[++tot].v=v,e[tot].w=w,e[tot].nxt=front[u],front[u]=tot;
		e[++tot].v=u,e[tot].w=0,e[tot].nxt=front[v],front[v]=tot;
	}
	bool bfs() {
		for(int i=1;i<=n;i++) dep[i]=0;
		queue<int> q;
		q.push(S);
		dep[S]=1;
		while(!q.empty()) {
			int u=q.front();
			q.pop();
			for(int i=front[u];i!=-1;i=e[i].nxt) {
				int v=e[i].v;
				if(dep[v]==0&&e[i].w) {
					dep[v]=dep[u]+1;
					q.push(v); 
				}
			}
		}
		return dep[T];
	}
	ll dfs(int cur,ll dis=1e18) {
		if(cur==T) return dis;
		ll flow=0,sum=0;
		for(int &i=head[cur];i!=-1;i=e[i].nxt) {
			int v=e[i].v;
			if(dep[v]==dep[cur]+1&&e[i].w) {
				flow=dfs(v,min(dis,e[i].w));
				dis-=flow,sum+=flow;
				e[i].w-=flow,e[i^1].w+=flow;
				if(!dis) break;
			}
		}
		return sum;
	}
	ll solve() {
		ll ans=0;
		while(bfs()) {
			for(int i=1;i<=n;i++) head[i]=front[i];
			ans+=dfs(S);
		}
		return ans;
	}
}
int m,n;  
bool check(int x,int y) {
	if(x>=1&&x<=m&&y>=1&&y<=n) return 1;
	else return 0;
}
int main() {
	scanf("%d%d",&m,&n);
	int tot=0;
	int S=n*m+1;
	int T=n*m+2;
	Dinic::S=S;
	Dinic::T=T;
	Dinic::n=n*m+2;
	Dinic::init();
	for(int i=1;i<=m;i++) {
		for(int j=1;j<=n;j++) {
			scanf("%d",&a[i][j]);
			mp[i][j]=++tot;
			if((i+j)%2==0) Dinic::add(S,tot,a[i][j]);
			else Dinic::add(tot,T,a[i][j]);
		}
	}
	for(int i=1;i<=m;i++) {
		for(int j=1;j<=n;j++) {
			if((i+j)%2==0) {
				if(check(i-1,j)) Dinic::add(mp[i][j],mp[i-1][j],1e9); 
				if(check(i,j-1)) Dinic::add(mp[i][j],mp[i][j-1],1e9);
				if(check(i+1,j)) Dinic::add(mp[i][j],mp[i+1][j],1e9);
				if(check(i,j+1)) Dinic::add(mp[i][j],mp[i][j+1],1e9);
			} 
		}
	}
	ll sum=0;
	for(int i=1;i<=m;i++) {
		for(int j=1;j<=n;j++) {
			sum+=a[i][j];
		}
	} 
	printf("%lld\n",sum-Dinic::solve());
	return 0;
}