洛谷P2774：方格取数问题（二分图最大权独立集，最小割模型）

题解：用最小割模型求解二分图最大权独立解。方法的可行性和证明可参考《最小割模型在信息学竞赛中的应用》论文。

建图时将二分图X部和S相连，容量为点权，Y部和T点相连，容量为点权，原图的边保留。
由于是网格图，不需要二分染色，观察可以发现下标（x,y) x 和 y同奇偶的一定是在同一个部，不同奇偶的在另一个部。

#include<bits/stdc++.h>
using namespace std;
const int maxn = 1e3 + 10;
const int maxm = 2e5 + 10;
const int inf = 0x3f3f3f3f;
struct ss{
	int u,v,w,nxt;
}edg[maxm];
int head[maxn],cnt,d[maxn],a[maxn][maxn],n,m;
int x[4] = {0,1,0,-1};
int y[4] = {1,0,-1,0};
void init() {
	cnt = 0;
	memset(head,-1,sizeof head);
}
void add(int u,int v,int w) {
	edg[cnt].v = v;
	edg[cnt].u = u;
	edg[cnt].w = w;
	edg[cnt].nxt = head[u];
	head[u] = cnt++;
}
bool bfs(int s,int t) {
	queue<int> q;
	memset(d,0,sizeof d);
	q.push(s);
	d[s] = 1;
	while(!q.empty()) {
		int top = q.front();
		q.pop();
		for(int i = head[top]; i + 1; i = edg[i].nxt) {
			int v = edg[i].v,w = edg[i].w;
			if(!d[v] && w) {
				d[v] = d[top] + 1;
				q.push(v);
			}
		}
	}
	return d[t] > 0;
}
int dfs(int s,int t,int inflow) {
	if(s == t || !inflow) return inflow;
	int flow = 0;
	for(int i = head[s]; i + 1; i = edg[i].nxt) {
		int v = edg[i].v,w = edg[i].w;
		if(w && d[v] == d[s] + 1) {
			int x = dfs(v,t,min(inflow,w));
			edg[i].w -= x;edg[i ^ 1].w += x;
			inflow -= x;flow += x;
			if(!inflow) break;
		}
	}
	if(!flow) d[s] = -2;
	return flow;
}
long long dinic(int s,int t) {
	long long ans = 0;
	while(bfs(s,t))
		ans += dfs(s,t,inf);
	return ans;
}
int main() {
	init();
	scanf("%d%d",&n,&m);
	int s = 0,t = n * m + 1;
	long long sum = 0;
	for(int i = 1; i <= n; i++)
		for(int j = 1; j <= m; j++) {
			scanf("%d",&a[i][j]);
			sum += a[i][j];
		}
	for(int i = 1; i <= n; i++) {
		for(int j = 1; j <= m; j++) {
			int p = (i - 1) * m + j;
			if((i & 1) == (j & 1)) {
				add(s,p,a[i][j]);
				add(p,s,0);
				for(int k = 0; k < 4; k++) {
					int tx = i + x[k],ty = j + y[k];
					if(tx >= 1 && tx <= n && ty >= 1 && ty <= m) {
						int q = (tx - 1) * m + ty;
						add(p,q,inf);
						add(q,p,0);
					}
				}
			}
			else {
				add(p,t,a[i][j]);
				add(t,p,0);
			}
		}
	}
	sum -= dinic(s,t);
	printf("%lld\n",sum);
	return 0;
}