HDU 3277 Marriage Match III

http://acm.hdu.edu.cn/showproblem.php?pid=3277

题意：有2N个孩子，其中有N个女生，N个男生，每一个女生可以找一个没有争吵过得男生组成一个家庭，并且可以和与她关系好的女生互换男生。与

HDU 3081 Marriage Match II  

不同的是，女生交换朋友的时候也不能和争吵过得男生组成家庭。

问能交换多少次。

题解：与3081那道题主要不同的地方是多了拆点。

　　　建边：将女生拆点为女生1，女生2。源点和女生1建边，边权依旧是mid；男生和汇点建边，边权是mid；女生1和女生2建边，边权是K；如果女生和男生可以组成家庭，女生1和男生建　　　　　　　边，边权为1，否则女生2和男生建边，边权为2。

T了很多发……优化了并查集的判断。抄了别人的ISAP的板子http://www.350351.com/bianchengyuyan/Cyuyan/318117.html

#include <iostream>
#include <cstring>
#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <string>
#include <vector>
#include <list>
#include <map>
#include <queue>
#include <stack>
#include <bitset>
#include <algorithm>
#include <numeric>
#include <functional>
#include <set>
#include <fstream>

using namespace std;

const int INF=0xfffffff;
const int maxm=10000010;
const int maxn=260;

struct edge{
    int u,v,cap,next;
}G[maxm];
struct node{
    int girl,boy;
}friend1[maxm];
int judge[maxn][maxn];
int idx;
int head[maxn*maxn],level[maxn*maxn];
int num[maxn*maxn];
int cur[maxn*maxn];
int pre[maxn*maxn];
int s,t,tt;
int par[maxn*maxn];
int rankh[maxn*maxn];
int N,M,K,F;

void init(int n)
{
    for(int i=0;i<=n;i++)
    {
        par[i]=i;
    }
}

int find(int x)
{
    if(par[x]==x){
        return x;
    }else{
        return par[x]=find(par[x]);
    }
}

void unite(int x,int y)
{
    x=find(x);
    y=find(y);
    if(x==y) return;
    else{
        par[x]=y;
    }
}

void build(int u,int v,int cap)
{
    G[idx].v=v;
    G[idx].cap=cap;
    G[idx].next=head[u];
    head[u]=idx++;
}

void add_edge(int u,int v,int cap)
{
    build(u,v,cap);
    build(v,u,0);
}


void bfs()
{
    memset(level,-1,sizeof(level));
    memset(num,0,sizeof(num));
    queue<int>q;
    q.push(t);
    level[t]=0;
    num[0]=1;
    while(!q.empty())
    {
        int u=q.front();
        q.pop();
        for(int i=head[u];i!=-1;i=G[i].next)
        {
            int v=G[i].v;
            if(level[v]==-1)
            {
                level[v]=level[u]+1;
                num[level[v]]++;
                q.push(v);
            }
        }
    }
}

int ISAP()
{
    memcpy(cur,head,sizeof(cur));
    bfs();
    int flow=0;
    int u=pre[s]=s;
    while(level[s]<tt)
    {
        if(u==t){
            int f=INF,pos;
            for(int i=s;i!=t;i=G[cur[i]].v)
            {
                if(f>G[cur[i]].cap){
                    f=G[cur[i]].cap;
                    pos=i;
                }
            }
            for(int i=s;i!=t;i=G[cur[i]].v)
            {
                G[cur[i]].cap-=f;
                G[cur[i]^1].cap+=f;
            }
            flow+=f;
            u=pos;
        }
        int i;
        for(i=cur[u];i!=-1;i=G[i].next)
        {
            if(level[G[i].v]+1==level[u]&&G[i].cap) break;
        }
        if(i!=-1){
            cur[u]=i;
            pre[G[i].v]=u;
            u=G[i].v;
        }
        else{
            if(--num[level[u]]==0) break;
            int mind=tt;
            for(int i=head[u];i!=-1;i=G[i].next)
            {
                if(mind>level[G[i].v]&&G[i].cap){
                    mind=level[G[i].v];
                    cur[u]=i;
                }
            }
            level[u]=mind+1;
            num[level[u]]++;
            u=pre[u];
        }
    }
    return flow;
}

void build_graph(int cap)
{
    memset(head,-1,sizeof(head));
    idx=0;
    for(int i=1;i<=N;i++)
    {
        for(int j=1;j<=N;j++)
        {
            if(judge[find(i)][j]){
                add_edge(i,2*N+j,1);
            }
            else{
                add_edge(N+i,2*N+j,1);
            }
        }
    }
    for(int i=1;i<=N;i++)
    {
        add_edge(s,i,cap);//源点-女生(1~N)
        add_edge(2*N+i,t,cap);//男生-汇点(2*N+1~3*N)
        add_edge(i,N+i,K);//女生1-女生2
    }
}

int main()
{
   // freopen("/Users/apple/Desktop/暑假/27(2)/27(2)/in","r",stdin);
    int T;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d%d%d%d",&N,&M,&K,&F);
        s=0;
        t=3*N+1;
        tt=t+1;
        memset(judge,0,sizeof(judge));
        
        for(int i=0;i<M;i++)
            scanf("%d%d",&friend1[i].girl,&friend1[i].boy);
        
        init(3*N);
        
        for(int i=0;i<F;i++)
        {
            int g1,g2;
            scanf("%d%d",&g1,&g2);
            unite(g1,g2);
        }
        for(int i=0;i<M;i++)
        {
            int x=friend1[i].girl,y=friend1[i].boy;
            judge[find(x)][y]=1;
        }
        
        int low=0;
        int high=N;
        int res=0;
        while(low<=high)
        {
            int mid=(low+high)/2;
            build_graph(mid);
            if(ISAP()>=N*mid) {
                res=mid;
                low=mid+1;
            }
            else{
                high=mid-1;
            }
        }
        printf("%d\n",res);
    }
    
    return 0;
}