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

慢啊,其实建边可以不用每次重建,改一改源汇点的边就好了。