方格取数[网络流]

传送门

最大流=最小割

所以我们每个点向相邻的建边 , 用总的减去最大流就好


#include<bits/stdc++.h>
#define N 150
#define M 100050
#define inf 0x3fffffff
using namespace std;
int first[N*N],next[M],to[M],w[M],tot=1;
int n,m,val[N][N],sum,st,ed,dis[N*N];
int fx[5]={0,1,-1,0,0},fy[5]={0,0,0,1,-1};
void add(int x,int y,int z){
	next[++tot]=first[x],first[x]=tot,to[tot]=y,w[tot]=z;
	next[++tot]=first[y],first[y]=tot,to[tot]=x,w[tot]=0;
}
bool bfs(){
	memset(dis,127,sizeof(dis)); int Inf=dis[0];
	dis[st]=0; queue<int> q; q.push(st);
	while(!q.empty()){
		int x=q.front(); q.pop();
		for(int i=first[x];i;i=next[i]){
			int t=to[i]; if(dis[t]==Inf && w[i]){
				dis[t]=dis[x]+1; q.push(t);
				if(t==ed) return true;
			}
		}
	}return false;
}
int dfs(int u,int flow){
	if(u==ed) return flow;
	int ans=0;
	for(int i=first[u];i;i=next[i]){
		int t=to[i];
		if(dis[t]==dis[u]+1 && w[i]){
			int delta = dfs(t,min(flow,w[i]));
			flow -= delta , ans += delta;
			w[i] -= delta , w[i^1] += delta;
			if(!flow) break;
		}
	}return ans;
}
int dinic(){
	int ans=0;
	while(bfs()){
		if(int x=dfs(st,inf)) ans += x;
	}
	return ans;
}
int main(){
	scanf("%d%d",&n,&m); st=0,ed=n*m+1;
	for(int i=1;i<=n;i++)
		for(int j=1;j<=m;j++){ 
			scanf("%d",&val[i][j]),sum+=val[i][j];
			if((i+j)&1) add(st,(i-1)*m+j,val[i][j]); 
			else add((i-1)*m+j,ed,val[i][j]);
		} 
		
	for(int i=1;i<=n;i++)
		for(int j=1;j<=m;j++)
			if((i+j)&1)
				for(int k=1;k<=4;k++){
				  	int x=i+fx[k],y=j+fy[k];
					if(x<1||x>n||y<1||y>m) continue;
					add((i-1)*m+j,(x-1)*m+y,inf);	
				}
	printf("%d",sum-dinic()); return 0;
} 

猜你喜欢

转载自blog.csdn.net/sslz_fsy/article/details/84927388