HDU 3277 Marriage Match III

题意

n个女孩n个男孩，每个女孩可以和一些男孩配对，此外还可以和k个任意的男孩配对，然后有些女孩是朋友，满足这个朋友圈里面的人，如果有一个能和某个男孩配对，其他就都可以，然后每轮要求每个女孩匹配到一个男孩，且每轮匹配到的都不同，问最多能匹配几轮

解

二分轮数，源向妹子连mid的边，妹子向满足条件的男孩连1的边，再连向一个点，容量为k，再连不满足的男孩连1的边，跑最大流，判断满流即可

代码

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
#include<cstdlib>
#define For(i,j,k) for(int i=(j);i<=(int)k;i++)
#define Forr(i,j,k) for(int i=(j);i>=(int)k;i--)
#define Set(a,b) memset(a,b,sizeof a)
#define Rep(i,u) for(int i=Begin[u],v=to[i];i;i=Next[i],v=to[i])
using namespace std;
const int N=760,M=63510,INF=0x3f3f3f3f;
template <class T> inline void read(T &x){
    x=0;char c=getchar();int f(0);
    while(c<'0'||c>'9')f|=(c=='-'),c=getchar();
    while(c>='0'&&c<='9')x=(x<<1)+(x<<3)+(c^48),c=getchar();
    x=f?-x:x;
}
template <class T> inline void chkmin(T &a,T b){a=a>b?b:a;}
template <class T> inline void chkmax(T &a,T b){a=a<b?b:a;}
int Begin[N],Next[M<<1],to[M<<1],e=1,n,m,f[M<<1],k,K;
inline void add(int x,int y,int z){
    to[++e]=y,Next[e]=Begin[x],Begin[x]=e,f[e]=z;
    to[++e]=x,Next[e]=Begin[y],Begin[y]=e,f[e]=0;
}
struct isap{
    int s,t,gap[N],d[N],cur[N],p[N];
    #define Retr(u) for(u=t;u^s;u=to[p[u]^1])
    inline void init(){
        queue<int>q;
        For(i,1,t)d[i]=INF,gap[i]=0,cur[i]=Begin[i];
        d[t]=0;gap[0]++;
        q.push(t);
        while(!q.empty()){
            int r=q.front();q.pop();
            Rep(i,r)
                if(f[i^1]&&d[v]>d[r]+1)
                    ++gap[d[v]=d[r]+1],q.push(v);
        }
    }
    inline int Augment(){
        int a=INF,u;
        Retr(u)chkmin(a,f[p[u]]);
        Retr(u)f[p[u]]-=a,f[p[u]^1]+=a;
        return a;
    }
    inline int Maxflow(int _s,int _t){
        int i,flow=0,u=s=_s;t=_t;
        init();
        while(d[s]<t){
            if(u==t)flow+=Augment(),u=s;
            for(i=cur[u];i;i=Next[i])
                if(f[i]&&d[to[i]]+1==d[u])break;
            if(!i){
                if(--gap[d[u]]==0)break;
                d[u]=t;cur[u]=Begin[u];
                Rep(j,u)
                    if(f[j])chkmin(d[u],d[v]+1);
                ++gap[d[u]];
                if(u^s)u=to[p[u]^1];
            }else cur[u]=i,p[u=to[i]]=i;
        }
        return flow;
    }
}F;
struct ufs{
    int f[N];
    inline void Init(int n){For(i,1,n)f[i]=i;}
    inline int Find(int x){return x==f[x]?x:f[x]=Find(f[x]);}
    inline bool Union(int x,int y){
        int a=Find(x),b=Find(y);
        return a^b?f[a]=b,1:0;
    }
}U;
int s_,t_;
int re[N][N],l[M][2];

inline bool check(int turn){
    Set(Begin,0),e=1;
    For(i,1,n)add(s_,i,turn),add(i,i+n*2,K);
    For(i,1,n)
        For(j,1,n)
            if(re[U.Find(i)][j])
                add(i,j+n,1);
            else add(i+n*2,j+n,1);
    For(i,n+1,n+n)
        add(i,t_,turn);
    return F.Maxflow(s_,t_)==n*turn;
}
inline void solve(){
    int L=0,R=n,mid,ans;
    while(L<=R){
        mid=(L+R)/2;
        if(check(mid))L=mid+1,ans=mid;
        else R=mid-1;
    }
    printf("%d\n",ans);
}
inline void work(){
    read(n),read(m),read(K),read(k);
    Set(re,0);
    s_=n*3+1,t_=s_+1;
    U.Init(n);
    For(i,1,m){
        int u,v;
        read(l[i][0]),read(l[i][1]);
    }
    For(i,1,k){
        int u,v;
        read(u),read(v);
        U.Union(u,v);
    }
    For(i,1,m)
        re[U.Find(l[i][0])][l[i][1]]=1;
    solve();
}
int main(){
    int T;
    read(T);
    while(T--)work();
    return 0;
}