HDU3081 Marriage Match II(最大流)

n个女孩和n个男孩要玩过家家的游戏。每个女孩要选择一个与她没有吵过架的男孩假装男朋友，如果这个女孩的朋友与某个男孩没有吵过架，那么这个女孩也可以找那个男孩假装男朋友。每次必须找不同的人，问这个游戏最多能进行多少轮。
把女孩放在一个并查集中然后连边，每个女孩向符合的男孩们连一条容量为1的边。假设这个游戏能进行k轮，那么源点s向每个女孩连容量为k的边，表示每个女孩都能找个k个男友，同理，每个男孩都向汇点t连容量为k的边。之后跑最大流，判断最大流是否等于n*k。于是乎直接二分最大轮数。
最先保存关系的数组开小了，奇葩地MLE好多次。
结果二分又写丑了，WA了，到网上去抄了抄二分的代码才过。

#include<cstdio>
#include<cstring>
#include<vector>
using namespace std;
inline int Min(int a,int b)
{return a<b?a:b;}
struct E
{
    int v,w,op;
    E(){}
    E(int a,int b,int c)
    {v = a; w = b; op = c;}
};
vector<E> g[210];
int n,m,f,fa[210],rela[10010][2];
int a,b,s,t,d[210],vd[210],flow;
bool mp[210][210];
int find(int x)
{
    if(fa[x] == x) return x;
    else return fa[x] = find(fa[x]);
}
void myunion(int x,int y)
{
    x = find(x);
    y = find(y);
    if(x != y) fa[y] = x;
}
int aug(int i,int augco)
{
    int j,mind = t-1,augc = augco,delta,sz = g[i].size();
    if(i == t) return augco;

    for(j = 0; j < sz; j++)
    {
        int v = g[i][j].v;
        if(g[i][j].w)
        {
            if(d[i] == d[v]+1)
            {
                delta = Min(augc,g[i][j].w);
                delta = aug(v,delta);
                g[i][j].w -= delta;
                g[v][g[i][j].op].w += delta;
                augc -= delta;
                if(d[s] >= t) return augco - augc;
                if(augc == 0) break;
            }
            if(d[v] < mind) mind = d[v];
        }
    }
    if(augc == augco)
    {
        vd[d[i]]--;
        if(vd[d[i]] == 0) d[s] = t;
        d[i] = mind+1;
        vd[d[i]]++;
    }
    return augco-augc;
}
void sap()
{
    flow = 0;
    memset(d,0,sizeof d);
    memset(vd,0,sizeof vd);
    vd[0] = t;
    while(d[s] < t)
        flow += aug(s,0x3f3f3f3f);
}
bool check(int val)
{
    memset(mp,0,sizeof mp);
    s = 2*n+1,t = s+1;
    for(int i = 1; i <= t; i++) g[i].clear();
    for(int i = 1; i <= n; i++)
    {
        g[s].push_back(E(i,val,g[i].size()));
        g[i].push_back(E(s,0,g[s].size()-1));
        g[i+n].push_back(E(t,val,g[t].size()));
        g[t].push_back(E(i+n,0,g[i+n].size()-1));
    }
    for(int i = 1; i <= m; i++)
    {
        int u = rela[i][0],v = rela[i][1];
        for(int j = 1; j <= n; j++)
        {
            if(find(u) == find(j)&&!mp[j][v])
            {    
                mp[j][v] = 1;
                g[j].push_back(E(v+n,1,g[v+n].size()));
                g[v+n].push_back(E(j,0,g[j].size()-1));
            }
        }
    }
    sap();
    if(flow >= n*val) return true;
    else return false;
}
void init()
{
    memset(rela,0,sizeof rela);
    memset(fa,0,sizeof fa);
    for(int i = 1; i <= n; i++)
        fa[i] = i;
}
int main()
{
    int T;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d%d%d",&n,&m,&f);
        init();
        for(int i = 1; i <= m; i++)
        {    
            scanf("%d%d",&rela[i][0],&rela[i][1]);
        }
        for(int i = 1; i <= f; i++)
        {
            scanf("%d%d",&a,&b);
            myunion(a,b);
        }
        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);
    }
}