【HDU 1569】 方格取数2 --最大点权独立集

【题目大意】

给你一个m*n的格子的棋盘，每个格子里面有一个非负数。
从中取出若干个数，使得任意的两个数所在的格子没有公共边，就是说所取数所在的2个格子不能相邻，并且取出的数的和最大。

(n,m<=50)

【分析】

将图黑白染色，然后按照矛盾关系建成二分图。

那么求一个最大点权独立集即可。

建图：

i + j 为奇数的与S节点相连，边的权值为棋盘上对应位置的值。

i + j 为偶数的与T节点相连，边的权值为棋盘上对应位置的值。

然后让棋盘上相邻之间的节点用边相连，边的权值为INF。

定理： 
1、最大点权独立集 = sum - 最小点权覆盖集。 
2、最小点权覆盖集 = 最小割 = 最大流 

【代码】

#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <iostream>
#include <algorithm>
#include <queue>
#include <deque>

using namespace std;

#define rep(i,a,b) for(int i=(a);i<=(b);++i)
#define rrep(i,b,a) for(int i=(b);i>=(a);--i)
#define pf printf
#define sf scanf
#define p_b push_back
#define INF (~0U>>3)
#define ll long long
#define MAXR 55
#define MAXN 2510
#define MAXM 50010
#define id(x,y) (((x)-1)*M+(y))

struct edge{
	int v,c;
	edge* next;
	edge* opt;
}*head[MAXN],Edge[MAXM];
int N,M;
int S,T;
int idx;
int ans;
int dis[MAXN];
int vd[MAXN];
int mat[MAXR][MAXR];

void Read(int& x){
	char tt=getchar();
	while(tt<'0'||'9'<tt) tt=getchar();
	for(x=0;'0'<=tt&&tt<='9';x=(x<<1)+(x<<3)+tt-'0',tt=getchar());
}

void Addedge(int u,int v,int c){
	idx++;
	Edge[idx].v=v;
	Edge[idx].c=c;
	Edge[idx].next=head[u];
	Edge[idx].opt=Edge+idx+1;
	head[u]=Edge+idx;
	
	idx++;
	Edge[idx].v=u;
	Edge[idx].c=0;
	Edge[idx].next=head[v];
	Edge[idx].opt=Edge+idx-1;
	head[v]=Edge+idx;
}

void Bfs(){
	queue <int> Q;
	Q.push(T);
	dis[T]=0;
	while(!Q.empty()){
		int u=Q.front();
		Q.pop();
		for(edge* p=head[u];p;p=p->next)
		    if(dis[p->v]==-1){
		    	vd[dis[p->v]=dis[u]+1]++;
		    	Q.push(p->v);
		    }
	}
}

int Sap(int u,int flow){
	if(u==T) return flow;
	int delta=0;
	for(edge* p=head[u];p;p=p->next)
	    if(p->c>0&&dis[u]==dis[p->v]+1){
	    	int tmp=Sap(p->v,min(p->c,flow-delta));
	    	p->c-=tmp;
	    	p->opt->c+=tmp;
	    	delta+=tmp;
	    	if(delta==flow) return delta;
	    }
	if(dis[S]>T) return delta;
	if(!(--vd[dis[u]])) dis[S]=T+1;
	vd[++dis[u]]++;
	return delta;
}

void Init(){
	rep(i,1,N)
	    rep(j,1,M)
	        Read(mat[i][j]);
}

void Solve(){
	rep(i,1,N)
	    rep(j,1,M){
	    	ans+=mat[i][j];
	    	int u=id(i,j);
	    	if((i+j)&1){
	    		Addedge(S,u,mat[i][j]);
	    		if(i>1) Addedge(u,u-M,INF);
	    		if(i<N) Addedge(u,u+M,INF);
	    		if(j>1) Addedge(u,u-1,INF);
	    		if(j<M) Addedge(u,u+1,INF);
	    	}
	    	else Addedge(u,T,mat[i][j]);
	    }
	Bfs();
	while(dis[S]<=T) ans-=Sap(S,INF);
	pf("%d\n",ans);
}

int main(){
	while(sf("%d%d",&N,&M)!=EOF){
	    ans=idx=0;
    	S=0;T=N*M+1;
	    memset(head,0,sizeof head);
    	memset(dis,-1,sizeof dis);
    	memset(vd,0,sizeof vd);
	    Init();
	    Solve();
	}
	return 0;
}