BZOJ3774: 最优选择【网络流多元关系转化为二元关系】

题目描述：

$n,m<=50$

题目分析：

先考虑$x$和$x_1$的二元关系：

• 如果$x$不放并且$x_1$为否，则失去$v$的收益

然后考虑$x_1$与($x$的相邻格子) $y$的二元关系：

• 如果$x_1$为是，则$y$必须要选，如果$y$不选，将花费正无穷的代价。
• 值得注意的是，$x_1$选是或否，所割的边的代价都是0，故当$y$都选时，$x_1$为是必然比$x_1$为否更优，所以另外一种二元关系（$x_1$为否，则$y$都选的代价为正无穷）可以不予表示。（此条为个人见解）

依据上述二元关系解方程，再进行建图(需要用到二分图关系)即可解决。

Code：

#include<cstdio>
#include<cstring>
#include<algorithm>
#define maxn 5005
#define maxm 30005
#define Clear(a,b) memset(a,0,(b)<<2)
using namespace std;
const int inf = 0x3f3f3f3f;
int n,m,S,T,sum;
int fir[maxn],cur[maxn],nxt[maxm],to[maxm],c[maxm],tot=1;
inline void line(int x,int y,int z){
	nxt[++tot]=fir[x],fir[x]=tot,to[tot]=y,c[tot]=z;
	nxt[++tot]=fir[y],fir[y]=tot,to[tot]=x,c[tot]=0;
}
namespace Maxflow{
	int d[maxn],vd[maxn],sz;
	int aug(int u,int augco){
		if(u==T) return augco;
		int need=augco,delta;
		for(int &i=cur[u];i;i=nxt[i]) if(c[i]&&d[u]==d[to[i]]+1){
			delta=aug(to[i],min(need,c[i])),c[i]-=delta,c[i^1]+=delta;
			if(!(need-=delta)||d[S]==sz) return augco-need;
		}
		if(!--vd[d[u]]) d[S]=sz;
		else vd[++d[u]]++,cur[u]=fir[u];
		return augco-need;
	}
	int SAP(){
		sz=T+1,Clear(d,sz),Clear(vd,sz),Clear(cur,sz);
		int flow=0;while(d[S]<sz) flow+=aug(S,inf);
		return flow;
	}
}
int p[55][55],dx[4]={1,-1,0,0},dy[4]={0,0,1,-1};
int main()
{
	int x,y;
	scanf("%d%d",&n,&m),S=0,T=n*m*2+1;
	for(int i=1;i<=n;i++) for(int j=1;j<=m;j++) p[i][j]=(i-1)*m+j;
	for(int i=1;i<=n;i++) for(int j=1;j<=m;j++) scanf("%d",&x),(i+j)&1?(line(p[i][j],T,x),0):(line(S,p[i][j],x),0);
	for(int i=1;i<=n;i++) for(int j=1;j<=m;j++) scanf("%d",&x),(i+j)&1?(line(p[i][j]+n*m,p[i][j],x),0):(line(p[i][j],p[i][j]+n*m,x),0),sum+=x;
	for(int i=1;i<=n;i++) for(int j=1;j<=m;j++)
		if((i+j)&1) {for(int k=0;k<4;k++) if(p[x=i+dx[k]][y=j+dy[k]]) line(p[x][y],p[i][j]+n*m,inf);}
		else {for(int k=0;k<4;k++) if(p[x=i+dx[k]][y=j+dy[k]]) line(p[i][j]+n*m,p[x][y],inf);}
	printf("%d",sum-Maxflow::SAP());
}