【hdu】3277 Marriage Match III【最大流】

题意：

n男n女，女的选择男朋友，每个女的有喜欢的男的和不喜欢的男的，每个女的有朋友，朋友与男朋友满足传递关系，即a和b是朋友，b喜欢c男生，那么a也喜欢c男生，问每次女的都选择不同的男朋友，并且每个女生有k次机会选择自己不喜欢的男生，如果所有女的都选择到了男朋友算一轮游戏，问这个游戏总共可以玩多少轮

题解：

看出来大概是一个类似婚姻匹配的游戏，然后加上了朋友圈的传递以及选择不喜欢的男生，很明显可以看出是一个二分答案+最大流的问题，那么接下来的问题就是如何建图了

我们每次二分答案mid，并且将每个女生拆点，女生连他喜欢的男生，女生'连他不喜欢的男生

1.vs向所有女生连一条容量为mid的边

2.所有女生向女生'连一条容量为k的边

3.女生向自己喜欢的男生连一条容量为1的边

4.女生'向自己喜欢的男生连一条容量为1的边

5.男生向vt连一条容量为mid的边

之前我在处理朋友圈的时候出了点问题，我是在之前读入女生喜欢的男朋友的时候将男朋友存入到vector中，然后在并查集操作的时候对vector进行去重的合并操作，TLE

之后看了下别人的题解，发现只需要vis标记就可以了，先把男女生对离线，对女生进行并查集操作会，将每一对女男直接加入到女的父节点上就好了，确实省了很多时间，虽然我觉得我的也不至于TLE，不管了，，，，、

#include<vector>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int MAXN=1010;
const int MAXM=100100;
const int inf=1<<30;
struct edge
{
    int v,cap,next;
}edge[MAXM];
int head[MAXN],pre[MAXN],cur[MAXN],level[MAXN],gap[MAXN];
int NE;
void init()
{
    NE=0;
    memset(head,-1,sizeof(head));
}
void add(int u,int v,int cap)
{
    edge[NE].v=v;edge[NE].cap=cap;
    edge[NE].next=head[u];head[u]=NE++;

    edge[NE].v=u;edge[NE].cap=0;
    edge[NE].next=head[v];head[v]=NE++;
}
int SAP(int vs,int vt,int NV)
{
    memset(pre,-1,sizeof(pre));
    memset(level,0,sizeof(level));
    memset(gap,0,sizeof(gap));
    for(int i=0;i<NV;i++)cur[i]=head[i];
    int u=pre[vs]=vs,maxflow=0,aug=-1;
    gap[0]=NV;
    while(level[vs]<NV){
loop:
        for(int &i=cur[u];i!=-1;i=edge[i].next){
            int v=edge[i].v;
            if(edge[i].cap&&level[u]==level[v]+1){
                aug==-1?aug=edge[i].cap:aug=min(aug,edge[i].cap);
                pre[v]=u;
                u=v;
                if(v==vt){
                    maxflow+=aug;
                    for(u=pre[u];v!=vs;v=u,u=pre[u]){
                        edge[cur[u]].cap-=aug;
                        edge[cur[u]^1].cap+=aug;
                    }
                    aug=-1;
                }
                goto loop;
            }
        }
        int minlevel=NV;
        for(int i=head[u];i!=-1;i=edge[i].next){
            int v=edge[i].v;
            if(edge[i].cap&&minlevel>level[v]){
                cur[u]=i;
                minlevel=level[v];
            }
        }
        gap[level[u]]--;
        if(gap[level[u]]==0)break;
        level[u]=minlevel+1;
        gap[level[u]]++;
        u=pre[u];
    }
    return maxflow;
}
bool vis[260][260];
int f[260];
int Find(int x)
{
    return f[x]==x?f[x]:f[x]=Find(f[x]);
}
void merged(int a,int b)
{
    int fa=Find(a);
    int fb=Find(b);
    if(fa!=fb)     f[fa]=fb;
}
struct node
{
    int a,b;
}p[MAXM];
int main()
{//freopen("C:\\Users\\Administrator\\Desktop\\input.txt","r",stdin);
    int T,i,j,n,m,k,t,a,b;
    scanf("%d",&T);
    while(T--){
        scanf("%d%d%d%d",&n,&m,&k,&t);
        for(i=1;i<=n;i++)f[i]=i;
        for(i=0;i<m;i++){
            scanf("%d%d",&p[i].a,&p[i].b);
        }
        for(i=0;i<t;i++){
            scanf("%d%d",&a,&b);
            merged(a,b);
        }
        memset(vis,0,sizeof(vis));
        for(i=0;i<m;i++){
            int a=p[i].a;
            int b=p[i].b;
            vis[Find(a)][b]=1;
        }
        int left=0,right=n;
        int mid,ans=0,vs=0,vt=3*n+1;
        while(left<=right){
            mid=(left+right)/2;init();
            for(i=1;i<=n;i++){
                add(i,n+i,k);
                add(vs,i,mid);
                add(2*n+i,vt,mid);
                for(j=1;j<=n;j++){
                    if(vis[Find(i)][j]) add(i,j+2*n,1);
                    else add(n+i,2*n+j,1);
                }
            }
            if(SAP(vs,vt,vt+1)==mid*n)ans=mid,left=mid+1;
            else right=mid-1;
        }
        printf("%d\n",ans);
    }
    return 0;
}