hust 1024 dance party

题目链接：http://acm.hust.edu.cn/problem.php?id=1024

题目大意：舞会上，男孩和女孩配对，求最大完全匹配个数，要求每个人最多与k个不喜欢的人配对，且每次都和不同的人配对。

题目思路：网络流+拆点

       这题可以看成求最多次能做几次二分图完全匹配，即最大流刚好为点的倍数，于是可以二分答案，然后用最大流判断是否可行，关键在于构图：

       二分答案ans；

       将每个男孩看成左边的点X，女孩为右边的点Y，然后将X分为与自己喜欢的跳舞Xa，不喜欢的Xb，Y同分为Ya，Yb，在源点src与所以Xa之间连一条容量为ans的有向边，在Ya与汇点dest之间连一条容量为ans的有向边，在Xa与Xb，Yb与Ya之间连一条有向边，之后在枚举所以男孩i，女孩j，如果i可以和j一起，则在Xai与Yaj之间连一条容量为1的有向边，否则在Xbi与Ybj之间连一条有向边，当前最大流如果为n*ans可行。。。（摘自网上）

#include<stdio.h>#include<stdlib.h>#include<string.h>#include<string>#include<queue>#include<algorithm>#include<vector>#include<stack>#include<list>#include<iostream>#include<map>using namespace std;#define inf 0x7f3f3f3f#define Max 500int max(int a,int b){    return a>b?a:b;}int min(int a,int b){    return a<b?a:b;}int dis[Max],gap[Max],pre[Max],cur[Max],p[Max];//int d[4][2]={0,1,1,0,0,-1,-1,0};int mp[220][220];int n,m,s,t,eid;struct node {     int to,next,c; }e[4*Max*Max],ee[4*Max*Max];void addedge(int u,int v,int c){    ee[eid].to=v;    ee[eid].c=c;    ee[eid].next=p[u];    p[u]=eid++;}int  ISAP(int st,int ed,int n,int count)   ///起点，终点，顶点数{    memset(dis, 0, sizeof(dis));    memset(gap, 0, sizeof(gap));  gap[0]=n;    memcpy(cur, p, sizeof(p));      ///memcpy!    int  i,flag,v,u=pre[st]=st,maxflow=0,aug=inf; //puts("akk");    while(dis[st] < n)    {        for(flag=0,i=cur[u];i!=-1; i=e[i].next)  /// cur[u]            if(e[i].c&& dis[u] == dis[e[i].to]+1)            {                flag = 1;                break;            }        if(flag)        {            if(aug > e[i].c)                aug = e[i].c;            v = e[i].to;            pre[v] = u;            cur[u] = i;            u = v;            if(u == ed)            {                for(u=pre[u]; 1;u=pre[u])    ///notice!                {                    e[cur[u]].c -= aug;                    e[cur[u]^1].c += aug;                    if(u==st) break;                   // puts("akkk");                }                maxflow += aug;                aug = inf;            }        }        else        {            int minx = n;            for(i=p[u]; i!=-1; i=e[i].next)                if(e[i].c&& dis[e[i].to]<minx)                {                    minx = dis[e[i].to];                    cur[u] = i;                }            if(--gap[dis[u]] == 0)                break;            dis[u] = minx+1;            gap[dis[u]]++;            u = pre[u];        }    }  // printf("Case %d:\n%d\n",count,maxflow);   // printf("%d\n",maxflow);    return maxflow;}int main(){    int m,n,t,count=1;    int u,v,c,i,j,k,x,y;    scanf("%d",&t);    while(t--)    {        scanf("%d%d%d",&n,&m,&k);        memset(mp,0,sizeof(mp));        memset(p,-1,sizeof(p));        eid=0;        while(m--)        {            scanf("%d%d",&u,&v);            mp[u][v]=1;            addedge(u,v+n,1);            addedge(v+n,u,0);        }        for(i=1;i<=n;i++)        {            addedge(i,i+2*n,k);            addedge(i+2*n,i,0);            addedge(i+3*n,i+n,k);            addedge(i+n,i+3*n,0);            addedge(0,i,0);            addedge(i,0,0);            addedge(i+n,4*n+1,0);            addedge(4*n+1,i+n,0);            for(j=1;j<=n;j++)            {                if(!mp[i][j])                {                    addedge(i+2*n,j+3*n,1);                    addedge(j+3*n,i+2*n,0);                }            }        }        int l,r,mid;        l=1;r=n;      //  printf("%d\n",eid); puts("akkkk");        while(l<=r)        {            mid=(l+r)>>1;            for(i=0;i<=eid;i++)            {                e[i]=ee[i];            }            for(i=p[0];i!=-1;i=e[i].next)            {                e[i].c=mid;                e[i^1].c=0;            }            for(i=p[4*n+1];i!=-1;i=e[i].next)            {                e[i].c=0;                e[i^1].c=mid;            }            if(ISAP(0,4*n+1,4*n+2,count++)==mid*n)                l=mid+1;            else                r=mid-1;        }      //  ISAP(0,n*m+1,n*m+2,count++);        printf("%d\n",r);    }}