hdu3277

hdu3277
题意：
女选男，结婚游戏，可以选自己喜欢的，也可以选朋友喜欢的，每次换男的结婚，还可以容忍没有选的男生k次 ，问最多能玩几次(sum)

分析
二分匹配的改编
二分匹配是源点到女孩是1，女孩到男孩是1,男孩到汇点是1

之前是源点到女孩是sum，女孩到男孩是1,男孩到汇点是sum（满流，每一个女孩都可以换至少sum次）

这是是源点到女孩是sum，女孩到男孩是1,女孩可以选k个不喜欢的男孩玩，新建一个源点，每个女孩到自己源点是k,自己源点到自己不喜欢男孩是1，男孩到汇点是sum（满流，每一个女孩都可以换至少sum次）
dinic(多路)

#include <bits/stdc++.h>

using namespace std;
typedef long long ll;
const ll INF=1e18;
const int maxn=805;
const int maxm=maxn*maxn;
struct edge
{
    int v,next;
    ll w;
};
edge no[maxm];
set<int >se[305];
set<int >se_[305];
int vis[305];
int fa[305];
int head[maxn];
int level[maxn];
int cnt;
queue <int> q;
int find_(int x)
{
    return x==fa[x] ?x: fa[x]=find_(fa[x]);
}
bool bfs(int s,int t,int n)
{
    for(int i=0;i<=n;i++)
        level[i]=-1;
    while(!q.empty())
    {
        q.pop();
    }
    q.push(s);
    level[s]=0;
    while(!q.empty())
    {
        int u=q.front();
        if(u==t)
            return true;
        q.pop();
        for(int i=head[u];i;i=no[i].next)
        {
            int v=no[i].v;
            ll w=no[i].w;
            if(w>0&&level[v]==-1)
            {
                level[v]=level[u]+1;
                q.push(v);
            }
        }
    }
    return false;
}
ll dfs(int u,int t,ll max_)
{
    if(u==t)
        return max_;
    ll sum=0;
    for(int i=head[u];i;i=no[i].next)
    {
        edge & e=no[i];
        if(e.w>0&&level[e.v]==level[u]+1)
        {
            ll f=dfs(e.v,t,min(max_-sum,e.w));
            sum+=f;
            no[i].w-=f;
            no[1+((i-1)^1)].w+=f;
            if(sum==max_)
                break;
        }
    }
    if(!sum)
        level[u]=-1;
    return sum;
}
ll dinic(int s,int t,int n)
{
    ll ans=0;
    while(bfs(s,t,n))
    {
        ll f;
        while((f=dfs(s,t,INF))>0)
            ans+=f;
    }
    return ans;
}
void init()
{
    for(int i=0;i<255;i++)
    {
        se[i].clear();
        se_[i].clear();
        fa[i]=i;
    }
    cnt=1;
    for(int i=0;i<maxn;i++)
    {
        head[i]=0;
    }
    //memset(vis,0,sizeof(vis));
}
void init2()
{
    cnt=1;
    for(int i=0;i<maxn;i++)
    {
        head[i]=0;
    }
}
void add(int x,int y,ll w)
{
    no[cnt].v=y;
    no[cnt].w=w;
    no[cnt].next=head[x];
    head[x]=cnt;
    cnt++;
}
void solve()
{
    int n,m,f,k;
    cin>>n>>m>>k>>f;
    int s=0,t=3*n+1;
    for(int i=1;i<=m;i++)
    {
        int x,y;
        cin>>x>>y;
        se[x].insert(y);
    }
    for(int i=1;i<=f;i++)
    {
        int x,y;
        cin>>x>>y;
        int fax=find_(x);
        int fay=find_(y);
        if(fax!=fay)
        {
            fa[fax]=fay;
            for(auto it:se[fax])
            {
                se[fay].insert(it);
            }
        }
    }
    for(int i=1;i<=n;i++)
    {
        memset(vis,0,sizeof(vis));
        int fai=find_(i);
        for(auto it:se[fai])
        {
            vis[it]=1;
        }
        for(int j=1;j<=n;j++)
        {
            if(!vis[j])
            {
                se_[i].insert(j);
            }
        }
    }
    int l=0,r=n;
    while(l<r)
    {
        ll res=0;
        int mid=(l+r+1)/2;
        for(int i=1;i<=n;i++)
        {
            int fai=find_(i);
            add(s,i,mid);
            add(i,s,0);
            res+=mid;
            for(auto it:se[fai])
            {
                add(i,it+2*n,1);
                add(it+2*n,i,0);
            }
            for(auto it:se_[i])
            {
                add(i+n,it+2*n,1);
                add(it+2*n,i+n,0);
            }
            add(i,i+n,k);
            add(i+n,i,0);
        }
        for(int i=1;i<=n;i++)
        {
            add(i+2*n,t,mid);
            add(t,i+n*2,0);
        }
        ll ans=dinic(s,t,3*n+1);
        if(ans==res)
        {
            l=mid;
        }
        else
        {
            r=mid-1;
        }
        init2();
    }
    cout<<l<<"\n";
}
int main()
{
    ios::sync_with_stdio(0);
    cin.tie(0),cout.tie(0);
    int __=1;
    cin>>__;
    while(__--)
    {
        init();
        solve();
    }
    return 0;
}

isap

#include <bits/stdc++.h>
using namespace std;

typedef long long ll;
const ll INF=1e18;
const int maxn=255;

int fa[maxn];
struct Edge
{
    int u,v;
};
Edge e[maxn*maxn];

struct edge
{
    int v,next;
    ll w;
};
int cnt,start,End,ans,head[maxn*3];
ll gap[maxn*3],dis[maxn*3];
edge no[maxn*9*maxn];
int map_[maxn][maxn];
void init()
{
    memset(head,-1,sizeof(head));
    cnt=0;
}
int find_(int x)
{
    return x==fa[x]?x:fa[x]=find_(fa[x]);
}
void add(int x,int y,ll w)
{
    no[cnt].v=y;
    no[cnt].w=w;
    no[cnt].next=head[x];
    head[x]=cnt;
    cnt++;
}
ll dfs(int u,ll minflow)
{
    if(u==End)return minflow;
    int i,v;
    ll f,min_dis=ans-1,flow=0;
    for(i=head[u];i!=-1;i=no[i].next)
    {
        v=no[i].v;
        if(no[i].w<=0)continue;
        if(dis[v]+1==dis[u])
        {
            f=dfs(v,no[i].w>minflow-flow?minflow-flow:no[i].w);
            no[i].w-=f;
            no[i^1].w+=f;
            flow+=f;
            if(flow==minflow)break;
            if(dis[start]>=ans)return flow;
        }
        min_dis=min_dis>dis[v]?dis[v]:min_dis;
    }
    if(flow==0)
    {
        if(--gap[dis[u]]==0)
            dis[start]=ans;
        dis[u]=min_dis+1;
        gap[dis[u]]++;
    }
    return flow;
}
ll isap()///ans等于End+1
{
    ll maxflow=0;
    memset(gap,0,sizeof(gap));
    memset(dis,0,sizeof(dis));
    gap[0]=ans;
    while(dis[start]<ans)
        maxflow+=dfs(start,INF);
    return maxflow;
}
void solve()
{
    int n,m,k,f;
    cin>>n>>m>>k>>f;
    for(int i=0;i<=n;i++)
    {
        fa[i]=i;
    }
    for(int i=1;i<=m;i++)
    {
        cin>>e[i].u>>e[i].v;
    }
    for(int i=1;i<=f;i++)
    {
        int x,y;
        cin>>x>>y;
        int fax=find_(x);
        int fay=find_(y);
        if(fax!=fay)
        {
            fa[fax]=fay;
        }
    }
    for(int i=1;i<=n;i++)
    {
        fa[i]=find_(i);
    }
    for(int i=1;i<=m;i++)
    {
        int u=find_(e[i].u);
        int v=e[i].v;
        map_[u][v]=1;
    }
    int l=0,r=n;
    start=0;End=3*n+1;
    ans=End+1;
    ll bns;
    while(l<=r)
    {
        init();
        int mid=(l+r)/2;
        /*for(int i=1;i<=n;i++)
        {
            add(start,i,mid);
            add(i,start,0);
        }*/
        for(int i=1;i<=n;i++)
        {
            add(start,i,mid);
            add(i,start,0);
            for(int j=1;j<=n;j++)
            {
                if(map_[fa[i]][j])
                {
                    add(i,j+2*n,1);
                    add(j+2*n,i,0);
                }
                else
                {
                    add(i+n,j+2*n,1);
                    add(j+2*n,i+n,0);
                }
            }
            add(i,i+n,k);
            add(i+n,i,0);
        }
        for(int i=1;i<=n;i++)
        {
            add(i+2*n,End,mid);
            add(End,i+2*n,0);
        }
        if(isap()==mid*n)
        {
            bns=mid;
            l=mid+1;
        }
        else
        {
            r=mid-1;
        }
    }
    cout<<bns<<"\n";
}
int main()
{
    ios::sync_with_stdio(0);
    cin.tie(0),cout.tie(0);
    int __=1;
    cin>>__;
    while(__--)
    {
        memset(map_,0,sizeof(map_));
        //init();
        solve();
    }
    return 0;
}