【BZOJ 1475】方格取数 最大点权独立集

首先可以很明确的是原图一定是二分图（i+j为奇数的和偶数的格子没有边），先摆出结论：最大独立点集=sum-最小割

原因：我们这样建图，s向所有i+j为奇数的格子连边，点权是格子上面的数大小，i+j为偶数的格子向t连边边权也是格子上的权值，然后相邻的格子之间连边，边权为inf。这样的话每割掉一条边（连向s或者t的）就相当于不选这个节点，点之间的边权为inf，边两边连接的点有你无我，有我无你，也就意味着不能同时选这两个点，最后由于是最小割，所以sum-最小割就是最大独立权值。

#include<cstdio>
#include<iostream>
#include<cstring>
#define maxn 50023
using namespace std;
int head[maxn],n,mat[51][51],id[51][51],last[maxn];
int ans,tot,q[maxn],s,t,cnt,h[maxn];
struct edge{int v,next,w;}e[maxn];
void adde(int a,int b,int c){
	e[tot].v=b,e[tot].next=head[a],e[tot].w=c;head[a]=tot++;
	e[tot].v=a,e[tot].next=head[b],e[tot].w=0;head[b]=tot++;
}
bool bfs(){
	int l=0,r=1;q[l]=s;
	for(int i=s;i<=t;i++)h[i]=-1;h[s]=0;
	while(l<r){
		int u=q[l++];
		for(int v,i=head[u];i!=-1;i=e[i].next){
			if(h[v=e[i].v]==-1&&e[i].w){
				h[v]=h[u]+1;
				q[r++]=v;
			}
		}
	}
	return h[t]!=-1;
}
int dfs(int u,int f){
	if(!f||u==t)return f;
	int used=0,w;
	for(int v,i=last[u];i!=-1;i=e[i].next)
		if(h[v=e[i].v]==h[u]+1&&e[i].w){
			last[u]=i;
			w=min(f-used,e[i].w);
			w=dfs(v,w);
			used+=w,e[i].w-=w,e[i^1].w+=w;
			if(used==f)return f;
		}
	if(!used)h[u]=-1;
	return used;
}
void dinic(){
	while(bfs()){
		for(int i=s;i<=t;i++)last[i]=head[i];
		ans-=dfs(s,1e9);
	}
}

void build(){
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n;j++){
			if(i+j&1){
				adde(s,id[i][j],mat[i][j]);
				if(i>=2)adde(id[i][j],id[i-1][j],1e9);
				if(j>=2)adde(id[i][j],id[i][j-1],1e9);
			}else{
				adde(id[i][j],t,mat[i][j]);
				if(i>=2)adde(id[i-1][j],id[i][j],1e9);
				if(j>=2)adde(id[i][j-1],id[i][j],1e9);
			}
		}
}

int main(){
	memset(head,-1,sizeof(head));
	scanf("%d",&n);t=n*n+1;
	for(int i=1;i<=n;i++)
	for(int j=1;j<=n;j++){
		scanf("%d",&mat[i][j]),id[i][j]=++cnt;
		ans+=mat[i][j];
	}
	build();
	dinic();
	printf("%d",ans);
	return 0;
}