BZOJ2127/LG1646 happiness 新建点最小割

问题描述

BZOJ2127

LG1646

---

题解

和文理分科差不多

收益最大 -> 损失最小 -> 最小割

分别新建点表示互相关系就行了

---

(mathrm{Code})

#include<bits/stdc++.h>
using namespace std;

const int INF=0x3f3f3f3f;

int n,m,S,T;
int Head[100*100*6+100],to[5000007],Next[5000007],w[5000007],tot=1;
int d[100*100*6+100];

void addedge(int x,int y,int z){
	to[++tot]=y,Next[tot]=Head[x],Head[x]=tot,w[tot]=z;
}
void add(int x,int y,int z){
	addedge(x,y,z);addedge(y,x,0);
}

bool bfs(void){
	memset(d,0,sizeof(d));
	queue<int>q;q.push(S);d[S]=1;
	while(!q.empty()){
		int x=q.front();q.pop();
		for(int i=Head[x];i;i=Next[i]){
			int y=to[i];
			if(d[y]||!w[i]) continue;
			d[y]=d[x]+1;q.push(y);
			if(y==T) return true;
		}
	}
	return false;
}

int dfs(int x,int flow){
	if(x==T) return flow;
	int rest=flow;
	for(int i=Head[x];i&&rest;i=Next[i]){
		int y=to[i];
		if(d[y]!=d[x]+1||!w[i]) continue;
		int k=dfs(y,min(rest,w[i]));
		if(!k) d[y]=0;
		else w[i]-=k,w[i^1]+=k,rest-=k;
	}
	return flow-rest;
}

int Dinic(void){
	int res(0),t;
	while(bfs()){
		while(t=dfs(S,INF)) res+=t;
	}
	return res;
}

int id(int x,int y,int p){
	return (p-1)*n*m+(x-1)*m+y;
}

int total;

void Init(void){
	scanf("%d%d",&n,&m);
	S=n*m*5+1,T=S+1;
	for(int i=1;i<=n;i++) for(int j=1,x;j<=m;j++){
		scanf("%d",&x);add(S,id(i,j,1),x);
		total+=x;
	}
	for(int i=1;i<=n;i++) for(int j=1,x;j<=m;j++){
		scanf("%d",&x);add(id(i,j,1),T,x);
		total+=x;
	}
	for(int i=1;i<n;i++) for(int j=1,x;j<=m;j++){
		scanf("%d",&x);total+=x;
		add(id(i,j,2),id(i,j,1),INF);
		add(id(i,j,2),id(i+1,j,1),INF);
		add(S,id(i,j,2),x);
	}
	for(int i=1;i<n;i++) for(int j=1,x;j<=m;j++){
		scanf("%d",&x);total+=x;
		add(id(i,j,1),id(i,j,3),INF);
		add(id(i+1,j,1),id(i,j,3),INF);
		add(id(i,j,3),T,x);
	}
	for(int i=1;i<=n;i++) for(int j=1,x;j<m;j++){
		scanf("%d",&x);total+=x;
		add(id(i,j,4),id(i,j+1,1),INF);
		add(id(i,j,4),id(i,j,1),INF);
		add(S,id(i,j,4),x);
	}
	for(int i=1;i<=n;i++) for(int j=1,x;j<m;j++){
		scanf("%d",&x);total+=x;
		add(id(i,j,1),id(i,j,5),INF);
		add(id(i,j+1,1),id(i,j,5),INF);
		add(id(i,j,5),T,x);
	}
}

void Work(void){
	int maxflow=Dinic();
	printf("%d
",total-maxflow);
}

int main(){
	Init();
	Work();
	return 0;
}