bzoj1305 [CQOI2009]dance跳舞
来源:互联网 发布:hosts 一个域名多个ip 编辑:程序博客网 时间:2024/05/18 03:23
题目
n<=50,不用说些什么了吧,网络流呀。
怎么跑呢?注意到,两个人之间的关系有喜欢与不喜欢,所以可以想到拆点,然后用关系建边,然后每人的两个点之间建一条容量为k的边,源点连向男生,女生连向汇点,这样就可以了。
日常语文不好,画个图吧。
3 1
YYN
NYY
YNY
图丑233
能看吧。。。
之后如何计算答案呢?
显然满足二分性质。。
源汇点的边连容量为二分的mid,然后满足这些边满流即可。
#include<bits/stdc++.h>#define inf ~0u>>1#define N 205#define M 6000using namespace std;int n,k,mx,ans,S,T;int like[N+1][N+1];char opt;int first[N+1],nex[M+1],to[M+1],flow[M+1],siz;int dis[N+1];inline char nc(){ static char buf[100000],*p1=buf,*p2=buf; return p1==p2&&(p2=(p1=buf)+fread(buf,1,100000,stdin),p1==p2)?EOF:*p1++;}inline char getc(){ char c=nc(); for(;!(c<='Z'&&c>='A');)c=nc(); return c;}inline int read(){ int x=0,b=1; char c=nc(); for(;!(c<='9'&&c>='0');c=nc())if(c=='-')b=-1; for(;c<='9'&&c>='0';c=nc())x=x*10+c-'0'; return x*b;}inline void add(int x,int y,int z){ nex[siz]=first[x]; first[x]=siz; to[siz]=y; flow[siz]=z; siz++;}inline void insert(int x,int y,int z){ add(x,y,z),add(y,x,0);}inline void build(int val){ siz=0,memset(first,-1,sizeof(first)); for(int i=1;i<=n;i++)insert(S,i,val); for(int i=n+1;i<=2*n;i++)insert(i,T,val); for(int i=1;i<=n;i++)insert(i,i+2*n,k); for(int i=n+1;i<=2*n;i++)insert(i+2*n,i,k); for(int i=1;i<=n;i++) for(int j=1;j<=n;j++) if(like[i][j])insert(i,j+n,1); else insert(i+2*n,j+3*n,1);}inline bool BFS(int s,int t){ queue <int> Q; for(int i=S;i<=T;i++)dis[i]=inf; dis[s]=0,Q.push(s); while(!Q.empty()) { int x=Q.front();Q.pop(); for(int i=first[x];i!=-1;i=nex[i]) { int u=to[i]; if(flow[i]&&dis[x]+1<dis[u]) { dis[u]=dis[x]+1; if(u==t)return true; Q.push(u); } } } if(dis[t]==inf)return false; return true;}inline int DFS(int s,int f,int t){ if(s==t)return f; int rest=f; for(int i=first[s];i!=-1;i=nex[i]) { int u=to[i]; if(flow[i]&&dis[u]==dis[s]+1&&rest) { int tmp=DFS(u,min(rest,flow[i]),t); if(tmp==0)dis[u]=0; flow[i]-=tmp; flow[i^1]+=tmp; rest-=tmp; } } return f-rest;}inline int Dinic(int s,int t){ int ans=0; while(BFS(s,t)) ans+=DFS(s,~0u>>1,t); return ans;}int main(){ //freopen("in.txt","r",stdin); n=read(),k=read(); S=0,T=4*n+1; for(int i=1;i<=n;i++) for(int j=1;j<=n;j++) { opt=getc(); if(opt=='Y')like[i][j]=1; else like[i][j]=0; } int l=0,r=50; while(l<=r) { int mid=(l+r)/2; build(mid); ans=Dinic(S,T); if(ans>=n*mid)mx=mid,l=mid+1; else r=mid-1; } cout<<mx; return 0;}
慢啊,其实建边可以不用每次重建,改一改源汇点的边就好了。
阅读全文
0 0
- BZOJ1305 [CQOI2009]dance跳舞
- bzoj1305【CQOI2009】dance 跳舞
- bzoj1305[CQOI2009]dance跳舞
- bzoj1305: [CQOI2009]dance跳舞
- bzoj1305 [CQOI2009]dance跳舞
- 【bzoj1305】【CQOI2009】【dance】【跳舞】
- BZOJ1305: [CQOI2009]dance跳舞
- bzoj1305 [CQOI2009]dance跳舞
- bzoj1305[CQOI2009]dance跳舞
- bzoj1305: [CQOI2009]dance跳舞
- 【bzoj1305】[CQOI2009]dance跳舞
- bzoj1305: [CQOI2009]dance跳舞
- bzoj1305 [CQOI2009]dance跳舞
- BZOJ1305: [CQOI2009]dance跳舞
- BZOJ1305 [CQOI2009]dance跳舞
- bzoj1305 [CQOI2009]dance跳舞 最大流 二分
- 【bzoj1305】[CQOI2009]dance跳舞 最大流
- 【BZOJ1305】[CQOI2009]dance跳舞【最大流】【二分】
- FastReport 导出pdf时中文乱码的解决办法
- leetcode 140. Word Break II
- bzoj1214 [HNOI2004]FTP服务器
- Java小程序之阶乘
- 数据库锁总结
- bzoj1305 [CQOI2009]dance跳舞
- 半夜冒着大雨跑网吧写文章是怎样的体验?
- 二分法的时间复杂度+算法的时间复杂度计算
- bzoj2190 [SDOI2008]仪仗队
- 为什么你的linux时间总是不准
- javaSE(java基础库)私人学习笔记
- js模块化规约
- JSP 9大作用域(内置对象)
- bzoj2956 模积和