BZOJ2127 happiness(最小割)

传送门

【题目分析】

Orz常数大到爆炸。。。。。。马也。。。。。

网站排名975我自豪!

其实这道题可以参考一下这道题,也是类似分两边,如果有若干个分在一起有额外的收益。

然后改了一下边权就可以了。。。。

但为什么建边的时候反向边容量为z???qwq真是令人智熄。。。求dalao讲解qwq

【附上近9kms的代码~】

#include<bits/stdc++.h>
using namespace std;
const int maxn=1e2+10;
const int MAXN=5e4+10;
const int MAXM=1e6+10;
const int INF=0x3f3f3f3f;

int n,m,cnt,s,t,x,sum;
int ma[maxn][maxn];
int head[MAXN],cur[MAXN],depth[MAXN];
int nxt[MAXM],to[MAXM],w[MAXM];

int Read(){
	int i=0,f=1;
	char c;
	for(c=getchar();(c>'9'||c<'0')&&c!='-';c=getchar());
	if(c=='-')
	  f=-1,c=getchar();
	for(;c>='0'&&c<='9';c=getchar())
	  i=(i<<3)+(i<<1)+c-'0';
	return i*f;
}

void Add(int x,int y,int z){
	nxt[cnt]=head[x];
	head[x]=cnt;
	to[cnt]=y;
	w[cnt]=z;
	cnt++;
}

void add(int x,int y,int z){
	Add(x,y,z);
	Add(y,x,z);
}

int bfs(){
	queue<int> q;
	memset(depth,0,sizeof(depth));
	depth[s]=1;
	q.push(s);
	while(!q.empty()){
		int u=q.front();
		q.pop();
		for(int i=head[u];i!=-1;i=nxt[i]){
			int v=to[i];
			if(w[i]&&!depth[v]){
				depth[v]=depth[u]+1;
				q.push(v);
			}
		}
	}
	if(depth[t]==0)
	  return 0;
	return 1;
}

int dfs(int u,int flow){
	if(u==t)
	  return flow;
	for(int &i=cur[u];i!=-1;i=nxt[i]){
		int v=to[i];
		if(depth[v]==depth[u]+1&&w[i]){
			int di=dfs(v,min(flow,w[i]));
			if(di){
				w[i]-=di;
				w[i^1]+=di;
				return di;
			}
		}
	}
	return 0;
}

int dinic(){
	int ans=0;
	while(bfs()){
		for(int i=s;i<=x;++i)
		  cur[i]=head[i];
		while(int d=dfs(s,INF))
		  ans+=d;
	}
	return ans;
}

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

int main(){
	memset(head,-1,sizeof(head));
	n=Read(),m=Read();
	s=0,t=n*m+1,x=t;
	for(int i=1;i<=n;++i){
	    for(int j=1;j<=m;++j){
	  	    int a=Read();
	  	    sum+=a;
	  	    add(s,id(i,j),a);
	    }
	}
	for(int i=1;i<=n;++i){
		for(int j=1;j<=m;++j){
			int a=Read();
			sum+=a;
			add(id(i,j),t,a);
		}
	}
	for(int i=1;i<n;++i){
		for(int j=1;j<=m;++j){
			int a=Read();
			sum+=a;
			x++;
			add(id(i,j),x,a),add(id(i+1,j),x,a);
			add(s,x,a);
		}
	}
	for(int i=1;i<n;++i){
		for(int j=1;j<=m;++j){
			int a=Read();
			sum+=a;
			x++;
			add(id(i,j),x,a),add(id(i+1,j),x,a);
			add(x,t,a);
		}
	}
	for(int i=1;i<=n;++i){
		for(int j=1;j<m;++j){
			int a=Read();
			sum+=a;
			x++;
			add(id(i,j),x,a),add(id(i,j+1),x,a);
			add(s,x,a);
		}
	}
	for(int i=1;i<=n;++i){
		for(int j=1;j<m;++j){
			int a=Read();
			sum+=a;
			x++;
			add(id(i,j),x,a),add(id(i,j+1),x,a);
			add(x,t,a);
		}
	}
	cout<<sum-dinic();
	return 0;
}

 

posted @ 2019-01-19 11:02  Ishtar~  阅读(158)  评论(0编辑  收藏  举报