【bzoj3894】文理分科 最小割

AC通道：http://www.lydsy.com/JudgeOnline/problem.php?id=3894

【题解】

一个人要么选文，要么选理，所以这个问题就是把所有的点分为2个点集，建图很容易

但是还有额外的满意值，如何处理呢？

我们可以创造一个点表示文科的额外值，从源点向此点连一条容量为满意值的边，然后从此点向与之相邻的4个点和自身连一条容量为+oo的边

理科的处理也是一样的。

这样图中割掉边就相当于不选这种方案的损失。

跑一边最小割即可。

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<ctime>
#include<cmath>
#include<algorithm>
using namespace std;
#define FILE "read"
#define INF 1000000000
#define up(i,j,n) for(int i=j;i<=n;i++)
#define dn(i,j,n) for(int i=j;i>=n;i--)
const int dx[5]={0,1,0,-1,0};
const int dy[5]={1,0,-1,0,0};
namespace INIT{
	char buf[1<<15],*fs,*ft;
	inline char getc() {return (fs==ft&&(ft=(fs=buf)+fread(buf,1,1<<15,stdin),fs==ft))?0:*fs++;}
	inline int read(){
		int x=0,f=1;  char ch=getc();
		while(!isdigit(ch))  {if(ch=='-')  f=-1;  ch=getc();}
		while(isdigit(ch))  {x=x*10+ch-'0';  ch=getc();}
		return x*f;
	}
}using namespace INIT;
struct node{int y,next,v,rel;}e[1001000];
int n,m,S,T,len,sum,A[105][105],B[105][105],C[105][105],D[105][105],Link[100100],q[100100],level[100100];
int cal(int i,int j) {return (i-1)*m+j;}
void insert(int x,int y,int v) {
	e[++len].next=Link[x];Link[x]=len;e[len].y=y;e[len].v=v;e[len].rel=len+1;
	e[++len].next=Link[y];Link[y]=len;e[len].y=x;e[len].v=0;e[len].rel=len-1;
}
bool bfs(){
	memset(level,-1,sizeof(level));
	int head=0,tail=1;  level[S]=0;  q[1]=S;
	while(++head<=tail){
		for(int i=Link[q[head]];i;i=e[i].next){
			if(e[i].v&&level[e[i].y]<0){
				q[++tail]=e[i].y;
				level[q[tail]]=level[q[head]]+1;
			}
		}
	}
	return level[T]>=0;
}
int MAXFLOW(int x,int flow){
	if(x==T)  return flow;
	int d=0,maxflow=0;
	for(int i=Link[x];i&&maxflow<flow;i=e[i].next){
		if(level[x]+1==level[e[i].y]&&e[i].v){
			if(d=MAXFLOW(e[i].y,min(e[i].v,flow-maxflow))){
				e[i].v-=d;  e[e[i].rel].v+=d;  maxflow+=d;
			}
		}
	}
	if(!maxflow) level[x]=-1;
	return maxflow;
}
void dinic(){
	int d=0,ans=0;
	while(bfs())  while(d=MAXFLOW(S,INF))  ans+=d;
	printf("%d\n",sum-ans);
}
int main(){
	freopen(FILE".in","r",stdin);
	freopen(FILE".out","w",stdout);
	n=read();  m=read();  S=0;  T=n*m*3+1;
	up(i,1,n)  up(j,1,m)  A[i][j]=read(),sum+=A[i][j];
	up(i,1,n)  up(j,1,m)  B[i][j]=read(),sum+=B[i][j];
	up(i,1,n)  up(j,1,m)  C[i][j]=read(),sum+=C[i][j];
	up(i,1,n)  up(j,1,m)  D[i][j]=read(),sum+=D[i][j];
	up(i,1,n)  up(j,1,m)  insert(S,cal(i,j),A[i][j]),insert(cal(i,j),T,B[i][j]);
	up(i,1,n)  up(j,1,m){
		insert(S,cal(i,j)+n*m,C[i][j]);
		insert(cal(i,j)+n*m*2,T,D[i][j]);
		up(k,0,4) {
			int xx=i+dx[k],yy=j+dy[k];
			if(xx>0&&xx<=n&&yy>0&&yy<=m)  {
				insert(cal(i,j)+n*m,cal(xx,yy),INF);
				insert(cal(xx,yy),cal(i,j)+n*m*2,INF);
			}
		}
	}
	dinic();
	return 0;
}