网络流二十四题 ————(十三）、P2774 方格取数问题 二分图网络流 最小割模型

一看到网格图我们就要想到：横纵坐标相加为偶数与为奇数的点不相邻这个tips

而题目说每个数不与其上下左右共四个数相连，与其他边相连即  —>  考虑删边可能更好做。

于是自然而然想到建立网络流最小割模型，删除权值最小的边，使得剩下的点不直接相连。

连边如下:

左部图为横纵坐标和是奇数的点，右部图是横纵坐标是偶数的点。

左部图向其四周不能连边的点连一条容量为无穷大的边。

S向左部图每个点连边，容量为点权。

右部图每个点向T连边，容量为点权。

求这个网络的最小割。

我们构造的这个网络S ->  T可达就说明存在禁止的边，(我们建立模型的时候就是往这里靠近的)

而中间连边容量为无穷大肯定不能删(最小割连边的套路！！)

左部图每条边对应一个点，删边等于删对应的点(套路！！)

右部图连边同理。

这个网络的跑最小割刚好是最小的删点方式，使得剩下的点互不直接相连！

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int M = 250000+7;
const int N = 500+7;
struct Dinic{
	//N个点，M条边 
	//add建双向边，然后D.gao ,最后输出maxflow  
	#define inf 0x3f3f3f3f
	ll maxflow;int s,t,n;
	int v[N],pre[N],d[N],now[N];
	ll incf[N];
	int head[N],cnt=1;
	struct EDGE{int to,nxt;ll w;}ee[M*2];
	inline void AD(int x,int y,ll w){ee[++cnt].nxt=head[x],ee[cnt].w=w,ee[cnt].to=y,head[x]=cnt;}
	inline void add(int x,int y,ll w){AD(x,y,w);AD(y,x,0);}
	inline bool bfs()//在残量网络上构造分层图 
	{
		memset(d,0,sizeof(d));
		queue<int>q;
		q.push(s);d[s]=1;
		while(q.size())
		{
			int x=q.front();q.pop();
			for(int i=head[x];i;i=ee[i].nxt)
			{
				int y=ee[i].to;ll w=ee[i].w;
				if(w&&!d[y])
				{
					q.push(y);
					d[y]=d[x]+1;
					if(y==t)return 1;
				}
			}
		}
		return false;
	}
	inline int dinic(int x,int flow)
	{
		if(x==t)return flow;
		ll rest = flow,k;
		for(int i=now[x];i&&rest;i=ee[i].nxt)
		{
			int y=ee[i].to;ll w= ee[i].w;
			now[x]=i;
			if(w&&d[y]==d[x]+1)
			{
				k=dinic(y,min(rest,w));
				if(!k)d[y]=0;//剪枝，去掉增广完毕的点
				ee[i].w-=k;
				ee[i^1].w+=k;
				rest-=k; 
			}
		}
		return flow - rest;
	}
	inline void gao()
	{
		int flow=0;
		while(bfs())
		{
			for(int i=0;i<=n;i++)now[i]=head[i];
			while(flow=dinic(s,inf))maxflow+=flow;
		}
	}
	inline void init(int nn,int S,int T)
	{
		cnt=1;maxflow=0;
		for(int i=0;i<=n;i++)head[i]=0;
		s=S,t=T,n=nn;
	}
}D;
//求解二分图匹配问题时的复杂度是：m*sqrt(n)的
int n,m;
int id(int x,int y){
	return (x-1)*m+y;
}
int a[110][110];
int main()
{
	int s,t;
	cin>>n>>m;
	int N=n*m+2;
	s=N-1,t=N; 
	D.init(N,s,t);//初始化节点个数，起点，终点 
	//注意点边个数NM,要根据具体题目进行更改！！！！！！！！！
	int sm=0;
	for(int i=1;i<=n;i++)
	for(int j=1;j<=m;j++){
		cin>>a[i][j];
		sm+=a[i][j];
		if((i+j)&1)D.add(s,id(i,j),a[i][j]);
		else D.add(id(i,j),t,a[i][j]);
	}
	for(int i=1;i<=n;i++)
	for(int j=1;j<=m;j++){
		if((i+j)&1)//只从左部图连向右部图 
		{
			if(i>1)D.add(id(i,j),id(i-1,j),inf);
			if(i<n)D.add(id(i,j),id(i+1,j),inf);
			if(j>1)D.add(id(i,j),id(i,j-1),inf);
			if(j<m)D.add(id(i,j),id(i,j+1),inf);
		}
	}
	D.gao();
	cout<<sm-D.maxflow<<endl;
	return 0;
}