bzoj3894: 文理分科
Posted AKCqhzdy
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了bzoj3894: 文理分科相关的知识,希望对你有一定的参考价值。
这题感觉就是最大权闭合子图嘛。
但是是最小割?我懒得写blog 先去看看一个sb的blog http://www.cnblogs.com/CHerish_OI/p/8126311.html
然后主要要想一想为啥要用sum-ans
经过sb CHerish_OI也就是蓝蓝的讲解,终于懂了。
这个图是用割分成两个集合,然后割出来的边是不要的。
#include<cstdio> #include<iostream> #include<cstring> #include<cstdlib> #include<algorithm> #include<cmath> using namespace std; const int inf=999999999; int n,m; struct node { int x,y,c,next,other; }a[5100000];int len,last[210000]; void ins(int x,int y,int c) { int k1,k2; len++;k1=len; a[len].x=x;a[len].y=y;a[len].c=c; a[len].next=last[x];last[x]=len; len++;k2=len; a[len].x=y;a[len].y=x;a[len].c=0; a[len].next=last[y];last[y]=len; a[k1].other=k2; a[k2].other=k1; } int st,ed,h[210000],list[210000]; bool bt_h() { memset(h,0,sizeof(h));h[st]=1; int head=1,tail=2;list[1]=st; while(head!=tail) { int x=list[head]; for(int k=last[x];k;k=a[k].next) { int y=a[k].y; if(h[y]==0&&a[k].c>0) { h[y]=h[x]+1; list[tail]=y; tail++; } } head++; } if(h[ed]==0)return false; return true; } int findflow(int x,int f) { if(x==ed)return f; int s=0; for(int k=last[x];k;k=a[k].next) { int y=a[k].y; if(h[y]==h[x]+1&&a[k].c>0&&s<f) { int t=findflow(y,min(a[k].c,f-s)); s+=t;a[k].c-=t;a[a[k].other].c+=t; } } return s; } //-----------wang luo liu------------- int sum; int point(int x,int y){return (x-1)*m+y;} void compostion() { st=n*m+1,ed=n*m+2; int x; for(int i=1;i<=n;i++) for(int j=1;j<=m;j++) scanf("%d",&x), ins(st,point(i,j),x), sum+=x; for(int i=1;i<=n;i++) for(int j=1;j<=m;j++) scanf("%d",&x), ins(point(i,j),ed,x), sum+=x; int cnt=n*m+2; for(int i=1;i<=n;i++) for(int j=1;j<=m;j++) { scanf("%d",&x); cnt++; ins(st,cnt,x); ins(cnt,point(i,j),inf); if(i-1>0) ins(cnt,point(i-1,j),inf); if(i+1<=n)ins(cnt,point(i+1,j),inf); if(j-1>0) ins(cnt,point(i,j-1),inf); if(j+1<=m)ins(cnt,point(i,j+1),inf); sum+=x; } for(int i=1;i<=n;i++) for(int j=1;j<=m;j++) { scanf("%d",&x); cnt++; ins(cnt,ed,x); ins(point(i,j),cnt,inf); if(i-1>0) ins(point(i-1,j),cnt,inf); if(i+1<=n)ins(point(i+1,j),cnt,inf); if(j-1>0) ins(point(i,j-1),cnt,inf); if(j+1<=m)ins(point(i,j+1),cnt,inf); sum+=x; } } int main() { scanf("%d%d",&n,&m); sum=0; compostion(); int ans=0; while(bt_h()==true) { ans+=findflow(st,2147483647); } printf("%d\\n",sum-ans); return 0; }
以上是关于bzoj3894: 文理分科的主要内容,如果未能解决你的问题,请参考以下文章