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;
}

 

转载于:https://www.cnblogs.com/Ishtar/p/10291854.html

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值