HDU3081Marriage Match II

HDU3081Marriage Match II
问题简述：
女生和男生配对，有些女生相互是朋友，每个女生也可以跟她朋友所配对的男生配对每次配对，每个女生都要跟不同的男生配对。问最多能配对几轮。

算法：二分+最大流，用并floyed处理女生之间的朋友关系，最少配0轮，最多配n轮，二分解之，源点向女生建边，男生向汇点建边，容量均为mid，女生跟所有能配对的男生连线，容量为1，如果最大流 = mid * n，就满足条件

#include <cstdio>
#include <cstdlib>
#include <iostream>
#include <cstring>
#include <cmath>
#include <algorithm>
#define Inf 0x7fffffff
using namespace std;

struct node
{
    int v,w,next;
}t[50010];
int  g[110][110],c[110][110],tu[110][110];
int st,ed,m,n,len,f;
int dis[220],cur[220],head[220],gap[220],pre[220];
void build()
{
    int i,j,k;
    for (k=1; k<=n; k++)
        for (i=1; i<=n; i++)
          for (j=1; j<=n; j++)
            if (!tu[i][j])
                tu[i][j]= ( tu[i][k]&&tu[k][j] );
    for (i=1; i<=n; i++)
        for (j=1; j<=n; j++)
        {
            c[i][j]=g[i][j];
            if (!c[i][j])
            for (k=1; k<=n; k++)
                if (tu[i][k] && g[k][j])
                {
                    c[i][j]=1;
                    break;
                }
        }
}

int Isap()
{
    int flow=0,i,d=Inf,u;
    bool flag=0;
    for (i=0; i<=ed; i++)
    {
        dis[i]=gap[i]=0;
        cur[i]=head[i];
    }
    gap[st]=ed+1;
    u=pre[st]=st;
    while (dis[st]<ed+1)
    {
        flag=0;// wrong
            printf("%d\n",dis[st]);
        for (int &j=cur[u]; j!=-1; j=t[j].next)
        {
            int v=t[j].v;
            if (t[j].w>0 &&  dis[v]+1==dis[u])
            {
                d=min(d,t[j].w);
                flag=true;
                pre[v]=u;
                u=v;
                if (u==ed)
                {
                    flow+=d;
                    while (u!=st)
                    {
                        u=pre[u];
                        t[cur[u]].w-=d;
                        t[cur[u]^1].w+=d;
                    }
                    d=Inf;
                }
            break;
            }
        }
        //printf("v==%d\n",u);
        if (flag)  continue;
        int minh=ed+1;
        for (int j=head[u]; j!=-1; j=t[j].next)
        {
            int v=t[j].v;
            if (t[j].w>0 && dis[v]<minh)
            {
                minh=dis[v];
                cur[u]=j;
            }
        }

        if ((--gap[dis[u]])==0) break;
        dis[u]=minh+1;
        gap[dis[u]]++;
        u=pre[u];
    }
    return flow;
}

void add(int x,int y,int w)
{
    t[len].w=w;
    t[len].v=y;
    t[len].next=head[x];
    head[x]=len++;
}

bool check(int k)
{
    int i,j;
    if (k==0)  return true;
    len=0;
    memset(head,-1,sizeof(head));
    for (i=1; i<=n; i++)
    {
        add(st,i,k);
        add(i,st,0);
    }
    for (i=1; i<=n; i++)
      for (j=1; j<=n; j++)
      if (c[i][j])
    {
        add(i,n+j,1);
        add(n+j,i,0);
    }
    for (i=n+1; i<=2*n; i++)
    {
        add(i,ed,k);
        add(ed,i,0);
    }
    int ans=Isap();
    if (ans==n*k)  return true;
     else return false;
}

void getans()
{
    int l=0,r=n;
    while (r-l>=2)
    {
        int mid=(l+r)/2;
        if (check(mid)) l=mid;
          else r=mid;
        //printf("%d %d\n",l,r);
    }
    for (int i=r; i>=l; i--)
    if (i<=n && i>=0)
    if (check(i))
    {
        printf("%d\n",i);
        return ;
    }
}

int main()
{
    int T,i,j,a,b;
    freopen("1.in","r",stdin);
    scanf("%d",&T);
    while (T--)
    {
        scanf("%d%d%d",&n,&m,&f);
        memset(g,0,sizeof(g));
        memset(tu,0,sizeof(tu));//³õʼ»¯´íÁË
        st=0,ed=2*n+1;
        len=0;
        for (i=1; i<=m; i++)
        {
            scanf("%d%d",&a,&b);
            g[a][b]=1;
        }
        for (i=1; i<=f; i++)
        {
            scanf("%d%d",&a,&b);
            tu[a][b]=1;
            tu[b][a]=1;
        }
        build();
        getans();
    }
    return 0;
}