HDU 4560 我是歌手 （网络流+拆点+二分）

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=4560

题目大意

有n个歌手，m种歌曲，每一场演唱会中
约束描述如下：每一场演唱会中每个歌手都要上台，
每个歌手都有其擅长曲目，一个歌曲全程被不擅长的歌手
演唱的次数不超过k，问最多可以半多少场演唱会。

题目分析 

不难发现题目中的答案隐含二分性质，
演唱会场数越多其出错的可能性会越大。
下面分析下约束关系，歌手与歌曲的约束关系，
歌手对其擅长的歌曲，约束应该是1，因为每一场要唱不同的歌曲，
但我们无法简单描述k的这个约束，所以我们决定用拆点，
把一个歌曲拆成：被擅长的歌手演唱，被不擅长的歌手演唱
这两种组合，即两个点，明显从被不擅长的歌手演唱的点到被擅长的歌手演唱的点其边
的容量应该是k，表示最大只有k次机会。
剩余的，建立超级源点和汇点，假设我们二分到的答案是举办p场演唱，
那么源点到歌手的边容量是p，被擅长的歌手演唱的曲目到汇点的容量也是p，如此
进行dinic即可。
如果求出来的流量等于场数*人数（总流量），那么说明这个枚举的场次可能是下界，往上再二分。
时间复杂度：O(n*n*m*logn).

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

#define debug puts("YES");
#define rep(x,y,z) for(int (x)=(y);(x)<(z);(x)++)
#define ll long long

#define lrt int l,int r,int rt
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1
#define root l,r,rt
#define mst(a,b) memset((a),(b),sizeof(a))
#define pii pair<int,int>
#define fi first
#define se second
#define mk(x,y) make_pair(x,y)
int mod=1e9+7;
const int maxn=1e4+10;
const int ub=1e6;
ll powmod(ll x,ll y){ll t; for(t=1;y;y>>=1,x=x*x%mod) if(y&1) t=t*x%mod; return t;}
ll gcd(ll x,ll y){return y?gcd(y,x%y):x;}

/*
题目大意：
有n个歌手，m种歌曲，每一场演唱会中
约束描述如下：每一场演唱会中每个歌手都要上台，
每个歌手都有其擅长曲目，一个歌曲全程被不擅长的歌手
演唱的次数不超过k，问最多可以半多少场演唱会。

题目分析：
不难发现题目中的答案隐含二分性质，
演唱会场数越多其出错的可能性会越大。
下面分析下约束关系，歌手与歌曲的约束关系，
歌手对其擅长的歌曲，约束应该是1，因为每一场要唱不同的歌曲，
但我们无法简单描述k的这个约束，所以我们决定用拆点，
把一个歌曲拆成：被擅长的歌手演唱，被不擅长的歌手演唱
这两种组合，即两个点，明显从被不擅长的歌手演唱的点到被擅长的歌手演唱的点其边
的容量应该是k，表示最大只有k次机会。
剩余的，建立超级源点和汇点，假设我们二分到的答案是举办p场演唱，
那么源点到歌手的边容量是p，被擅长的歌手演唱的曲目到汇点的容量也是p，如此
进行dinic即可。
如果求出来的流量等于场数*人数（总流量），那么说明这个枚举的场次可能是下界，往上再二分。
时间复杂度：O(n*n*m*logn).
*/
int n,m,l,k,a,b;
struct node{
    int u,nxt,flow;
}e[maxn<<1];
int head[maxn],tot=0;
void init(){
    mst(head,-1),tot=0;
}
void add(int x,int y,int v){
    e[tot]=node{y,head[x],v};
    head[x]=tot++;
    e[tot]=node{x,head[y],0};
    head[y]=tot++;
}
///
int d[maxn],S,T;
bool bfs(){
    queue<int> q;q.push(S);
    mst(d,-1),d[S]=0;
    while(!q.empty()){
        int x=q.front();q.pop();
        for(int i=head[x];~i;i=e[i].nxt){
            int y=e[i].u;
            if(e[i].flow&&d[y]==-1){
                d[y]=d[x]+1;
                q.push(y);
                if(y==T) return 1;
            }
        }
    }
    return 0;
}
int dfs(int x,int all){
    if(x==T) return all;
    int cost=0;
    ///cout<<x<<" dfs"<<endl;
    for(int i=head[x];~i;i=e[i].nxt){
        int y=e[i].u;
        if(e[i].flow&&d[y]==d[x]+1){
            int tmp=dfs(y,min(e[i].flow,all-cost));
            e[i].flow-=tmp;
            e[i^1].flow+=tmp;
            cost+=tmp;
            if(cost==all) break;
        }
    }
    if(cost==0) d[x]=-1;
    return cost;
}
int dinic(){
    int ret=0;
    while(bfs()) ret+=dfs(S,mod);
    return ret;
}
int g[105][105];
int solve(int ub){
    init();S=0,T=n+m+m+1;
    for(int i=1;i<=m;i++) add(n+i,n+m+i,k);
    rep(i,1,n+1) rep(j,1,m+1)
        if(g[i][j]) add(i,n+m+j,1);
        else add(i,n+j,1);
    rep(i,1,n+1) add(S,i,ub);
    rep(i,1,m+1) add(n+m+i,T,ub);
    return dinic();
}

int main(){
    int t;scanf("%d",&t);
    for(int ca=1;ca<=t;ca++){
        scanf("%d%d%d%d",&n,&m,&l,&k);
        mst(g,0);
        while(l--){
            int x,y;
            scanf("%d%d",&x,&y);
            g[x][y]=1;
        }
        int l=0,r=m,ans=0;
        while(l<r){
            int mid=(l+r+1)>>1;
            if(solve(mid)==mid*n) l=mid;
            else r=mid-1;
        }
        printf("Case %d: %d\n",ca,l);
    }
    return 0;
}