bzoj1305 [CQOI2009]dance跳舞

题目

n<=50,不用说些什么了吧,网络流呀。

怎么跑呢?注意到,两个人之间的关系有喜欢与不喜欢,所以可以想到拆点,然后用关系建边,然后每人的两个点之间建一条容量为k的边,源点连向男生,女生连向汇点,这样就可以了。

日常语文不好,画个图吧。
3 1
YYN
NYY
YNY

这里写图片描述

图丑233
能看吧。。。

之后如何计算答案呢?

显然满足二分性质。。

源汇点的边连容量为二分的mid,然后满足这些边满流即可。

#include<bits/stdc++.h>
#define inf ~0u>>1
#define N 205
#define M 6000
using 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;
}


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

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值