bzoj2127 happiness(最小割)
2017-12-25 10:16
232 查看
这题好神啊orz。
hzw题解传送门:http://hzwer.com/2422.html
利用最小割考虑。
对于原图中所有相邻的两个人A,B,我们建边:
s->A:cost[A文]+c[文][A][B]/2,s->B:cost[B文]+c[文][A][B]/2;
A->t:cost[A理]+c[理][A][B]/2,B->t:costB[理]+c[理][A][B]/2;
A<–>B:c[文][A][B]/2+c[理][A][B]/2
这样原图中的所有割,无非三种,分别对应了两人都选文(删掉了两人的理和都选理的贡献),都选理(删掉了两人的文和都选文的贡献),一个选文一个选理(删掉了另一个选文和另一个选理,都选理和都选文的贡献)。为了避免/2,给所有容量都*2,最后再/2就是答案。
hzw题解传送门:http://hzwer.com/2422.html
利用最小割考虑。
对于原图中所有相邻的两个人A,B,我们建边:
s->A:cost[A文]+c[文][A][B]/2,s->B:cost[B文]+c[文][A][B]/2;
A->t:cost[A理]+c[理][A][B]/2,B->t:costB[理]+c[理][A][B]/2;
A<–>B:c[文][A][B]/2+c[理][A][B]/2
这样原图中的所有割,无非三种,分别对应了两人都选文(删掉了两人的理和都选理的贡献),都选理(删掉了两人的文和都选文的贡献),一个选文一个选理(删掉了另一个选文和另一个选理,都选理和都选文的贡献)。为了避免/2,给所有容量都*2,最后再/2就是答案。
#include <cstdio> #include <cstring> #include <algorithm> #include <queue> using namespace std; #define ll long long #define inf 0x3f3f3f3f #define N 10010 inline int read(){ int x=0,f=1;char ch=getchar(); while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();} while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar(); return x*f; } int n,m,id[110][110],h ,num=1,T=10001,lev ,ans=0,tot=0,cur ; int owo=0,dx[]={0,0,1,-1},dy[]={1,-1,0,0},a[110][110],b[110][110],c[110][110]; struct edge{ int to,next,val; }data[100000]; inline void add(int x,int y,int val,int op){ data[++num].to=y;data[num].next=h[x];h[x]=num;data[num].val=val; data[++num].to=x;data[num].next=h[y];h[y]=num;data[num].val=val*op; } inline bool bfs(){ queue<int>q;memset(lev,0,sizeof(lev)); lev[0]=1;q.push(0); while(!q.empty()){ int x=q.front();q.pop(); for(int i=h[x];i;i=data[i].next){ int y=data[i].to;if(lev[y]||!data[i].val) continue; lev[y]=lev[x]+1;q.push(y); } }return lev[T]; } inline int dinic(int x,int low){ if(x==T) return low;int tmp=low; for(int i=h[x];i;i=data[i].next){ int y=data[i].to;if(lev[y]!=lev[x]+1||!data[i].val) continue; int res=dinic(y,min(data[i].val,tmp)); if(!res) lev[y]=0;else tmp-=res,data[i].val-=res,data[i^1].val+=res; if(!tmp) return low; }return low-tmp; } int main(){ // freopen("a.in","r",stdin); n=read();m=read(); for(int i=1;i<=n;++i) for(int j=1;j<=m;++j) id[i][j]=++owo; for(int i=1;i<=n;++i) for(int j=1;j<=m;++j) a[i][j]=read(),tot+=a[i][j],a[i][j]<<=1; for(int i=1;i<=n;++i) for(int j=1;j<=m;++j) b[i][j]=read(),tot+=b[i][j],b[i][j]<<=1; for(int i=1;i<n;++i) for(int j=1;j<=m;++j){ int x=read();tot+=x;a[i][j]+=x;a[i+1][j]+=x;c[i][j]=x; } for(int i=1;i<n;++i) for(int j=1;j<=m;++j){ int x=read();tot+=x;b[i][j]+=x;b[i+1][j]+=x;c[i][j]+=x; add(id[i][j],id[i+1][j],c[i][j],1); } for(int i=1;i<=n;++i) for(int j=1;j<m;++j){ int x=read();tot+=x;a[i][j]+=x;a[i][j+1]+=x;c[i][j]=x; } for(int i=1;i<=n;++i) for(int j=1;j<m;++j){ int x=read();tot+=x;b[i][j]+=x;b[i][j+1]+=x;c[i][j]+=x; add(id[i][j],id[i][j+1],c[i][j],1); } for(int i=1;i<=n;++i) for(int j=1;j<=m;++j) add(0,id[i][j],a[i][j],0),add(id[i][j],T,b[i][j],0); while(bfs()){memcpy(cur,h,sizeof(cur));ans+=dinic(0,inf);} printf("%d\n",tot-ans/2); return 0; }
相关文章推荐
- bzoj2127: happiness 最小割
- bzoj2127: happiness(最小割)
- bzoj2127 happiness 最小割
- BZOJ2127:happiness(最小割)
- bzoj2127: happiness 浅析一类最小割问题——二元关系的应用
- BZOJ 2127 Happiness 最小割
- bzoj2127(还没过+最小割)
- 【BZOJ2127】happiness 最小割
- 【COGS 1873】 [国家集训队2011]happiness(吴确) 最小割
- 【bzoj2127】happiness 最小割
- BZOJ2127: happiness
- [国家集训队2011]happiness(吴确) (最小割)
- bzoj2127 happiness
- 2127: happiness 最小割
- bzoj 2127 happiness(最小割)
- [bzoj2127]happiness——最小割
- [NBUT 1224 Happiness Hotel 佩尔方程最小正整数解]连分数法解Pell方程
- BZOJ_3894_文理分科&&BZOJ_2127_happiness_最小割
- BZOJ 2127 happiness 最小割
- BZOJ 2127: happiness( 最小割 )