codevs1907方格取数最大流

«问题描述：
在一个有m*n 个方格的棋盘中，每个方格中有一个正整数。现要从方格中取数，使任
意2 个数所在方格没有公共边，且取出的数的总和最大。试设计一个满足要求的取数算法。
«编程任务：
对于给定的方格棋盘，按照取数要求编程找出总和最大的数。

棋盘类问题通过黑白染色转化为二分图模型，

本题求不相邻的格点取值最大值，在黑白染色的二分图中，若两点相连，则他们连边

实际上本题求建图后，二分图的最大权独立集，来表示相邻两点是不能选的

最大权独立集=总权值-最小点权覆盖集=总权值-构图后最小割=总权值-最大流

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<algorithm>
#include<iostream>
#include<queue>
#define debug(x) cout<<#x<<"="<<x<<endl
using namespace std;
const int inf=0x3f3f3f3f;
const int N=1009;
const int qx[4]={0,0,1,-1};
const int qy[4]={1,-1,0,0};
struct aa
{
	int to,pre,cap,flow;
}edge[N*10];
int n,m,id[35][35],map[35][35],s,t;
int head[N],lev[N],tot,cur[N];
void addedge(int x,int y,int z)
{
	edge[++tot].to=y;edge[tot].pre=head[x];edge[tot].cap=z;head[x]=tot;
	edge[++tot].to=x;edge[tot].pre=head[y];edge[tot].cap=0;head[y]=tot;
}
bool bfs()
{
	memset(lev,0,sizeof(lev));
	lev[s]=1;
	queue<int> q;q.push(s);
	while (!q.empty())
	{
		int u=q.front();q.pop();
		for (int i=head[u];i;i=edge[i].pre)
		if (edge[i].cap>edge[i].flow&&!lev[edge[i].to])
		{
			lev[edge[i].to]=lev[u]+1;
			if (edge[i].to==t) return true;
			q.push(edge[i].to);
		}
	}
	return false;
}
int dfs(int u,int maxflow)
{
	if (u==t||maxflow==0) return maxflow;
	int ans=0;
	for (int &i=cur[u];i;i=edge[i].pre)
	if (lev[edge[i].to]==lev[u]+1)
	{
		int flow=dfs(edge[i].to,min(maxflow,edge[i].cap-edge[i].flow));
		ans+=flow;
		maxflow-=flow;
		edge[i].flow+=flow;
		edge[((i-1)^1)+1].flow-=flow;
		if (maxflow==0) return ans;
	}
	return ans;
}
int work()
{
	int ans=0;
	while (bfs())
	{
		for (int i=s;i<=t;i++) cur[i]=head[i];
		ans+=dfs(s,inf);
	}
	return ans;
}
int main()
{
	scanf("%d%d",&n,&m);
	int bj=0,total=0;
	s=0;
	for (int i=1;i<=n;i++)
		for (int j=1;j<=m;j++) id[i][j]=++bj; 
	t=bj+1;
	for (int i=1;i<=n;i++)
		for (int j=1;j<=m;j++) 
		scanf("%d",&map[i][j]),total+=map[i][j];
	for (int i=1;i<=n;i++)
		for (int j=1;j<=m;j++)
		if ((i+j)&1)
		{
			addedge(s,id[i][j],map[i][j]);
			for (int l=0;l<4;l++) 
			if (map[i+qx[l]][j+qy[l]])
			addedge(id[i][j],id[i+qx[l]][j+qy[l]],inf);
		}
		else addedge(id[i][j],t,map[i][j]);
	printf("%d",total-work());	
	return 0;
}