HDU 3081 Marriage Match II [二分最大流]+并查集?807★

题意：给定二分图：问最多有多少种完美匹配，限制条件是，每个人对应的另一个人，在所有的完美匹配中，都不一样；

思路：二分匹配数n,对X集的每个人与s连一条容量为n的边，对Y集的每个人与t连一条容量为n的边，求最大流，如果所有与s,t连接的边都满流，则n可行：

为什么满流就可行，就是说是不是一定能找到n种完美匹配，以下都是个人猜想的：

（一）.如果X,Y集人数刚好为n,也就是说这个图是完全图，那么找到n中完美匹配是一定的（可以归纳出来）

（二）.X,Y集人数比n大1，另外两个点叫A,B,

               1.如果A,B直间没有边，则在（一）的匹配的基础上，用A,B,替换一些匹配，（一）的匹配仍成立，因为A,B,都满流，所以替换了n次，所以A,B也参与了n次，新匹配成立

               2.如果A,B直间有边，则在（一）的匹配的基础上，用A,B,替换一些匹配，（一）的匹配仍成立，因为A,B,都满流，所以替换了n-1次，所以A,B也参与了n-1次，

再加上A与B的直接匹配也不影响（一）的匹配，新匹配成立

（三）X,Y集人数比n大的其他情况可以归纳出来

#include <iostream>
#include <cstring>
#include <cstdio>
#define MAXN 210
#define MAXM 50000
using namespace std;
struct node
{
    int u,v,w;
};
struct node2
{
    int x,y;
};
node2 ee[MAXM];
int first[210],p[MAXN];
node e[MAXM];
int next[MAXM],cc,n,m,f;
int gap[210],d[210],curedges[210],pre[210];
int map[105][105];
inline void add_edge(int u,int v,int w)
{
    e[cc].u=u;
    e[cc].v=v;
    e[cc].w=w;
    next[cc]=first[u];
    first[u]=cc;
    cc++;

    e[cc].u=v;
    e[cc].v=u;
    e[cc].w=0;
    next[cc]=first[v];
    first[v]=cc;
    cc++;
}
void make_set(int n)
{
    int i;
    for(i=0;i<=n;i++)
        p[i]=i;
}
int find_set(int i)
{
    int j=i;
    while(j!=p[j])
        j=p[j];
    return p[i]=j;
}
int Union(int x,int y)
{
    x=find_set(x);
    y=find_set(y);
    if(x==y)
        return 0;
    else
    {
        p[x]=y;
        return 1;
    }
}
void build(int t,int mid)
{
    cc=0;
    int i,j;
    memset(map,0,sizeof(map));
    memset(first,-1,sizeof(first));
    memset(next,-1,sizeof(next));
    for(i=1;i<=n;i++)
    {
        add_edge(0,i,mid);
        add_edge(n+i,t,mid);
    }
    for(i=0;i<m;i++)
    {
        int u=ee[i].x;
        int v=ee[i].y;
        for(j=1;j<=n;j++)
        {
            if(find_set(u)==find_set(j))
                if(!map[j][v])
                {
                    map[j][v]=1;
                    add_edge(j,n+v,1);
                }
        }
    }

}
int sap(int s,int t)
{
    int cur_flow,flow_ans=0,u,tmp,neck,i;
    memset(d,0,sizeof(d));
    memset(gap,0,sizeof(gap));
    memset(pre,-1,sizeof(sizeof(pre)));
    for(i=0;i<=t;i++)
        curedges[i]=first[i];
    gap[0]=t+1;
    u=s;
    while(d[s]<t+1)
    {
        if(u==t)
        {
            cur_flow=99999999;
            for(i=s;i!=t;i=e[curedges[i]].v)
            {
                if(cur_flow>e[curedges[i]].w)
                {
                    neck=i;
                    cur_flow=e[curedges[i]].w;
                }
            }
            for(i=s;i!=t;i=e[curedges[i]].v)
            {
                tmp=curedges[i];
                e[tmp].w-=cur_flow;
                e[tmp^1].w+=cur_flow;
            }
            flow_ans+=cur_flow;
            u=neck;
        }
        for(i=curedges[u];i!=-1;i=next[i])
            if(e[i].w&&d[u]==d[e[i].v]+1)
                break;
        if(i!=-1)
        {
            curedges[u]=i;
            pre[e[i].v]=u;
            u=e[i].v;
        }
        else
        {
            if(0==--gap[d[u]])
                break;
            curedges[u]=first[u];
            for(tmp=t+1,i=first[u];i!=-1;i=next[i])
                if(e[i].w)
                    tmp=min(tmp,d[e[i].v]);
            d[u]=tmp+1;
            ++gap[d[u]];
            if(u!=s)
                u=pre[u];
        }
    }
    return flow_ans;
}
int main()
{
    int T;
    scanf("%d",&T);
    while(T--)
    {
        memset(ee,0,sizeof(ee));
        memset(e,0,sizeof(e));
        scanf("%d%d%d",&n,&m,&f);
        int i,s=0,t=n*2+1;
        for(i=0;i<m;i++)
            scanf("%d%d",&ee[i].x,&ee[i].y);
        make_set(t);
        for(i=0;i<f;i++)
        {
            int u,v;
            scanf("%d%d",&u,&v);
            Union(u,v);
        }
        int l=0,r=n,ans=0;
        while(l<=r)
        {
            int mid=(l+r)>>1;
            build(t,mid);
            int res=sap(s,t);
            if(res>=mid*n)
            {
                ans=mid;
                l=mid+1;
            }
            else
                r=mid-1;
        }
        printf("%d\n",ans);
    }
    return 0;
}