【题目大意】
给你一个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;
}