hdu 4560 拆点最大流 ***

题意：

2013年一开始，一档音乐节目“我是歌手”就惊艳了大家一回。闲话少说，现在，你成为了这档节目的总导演，你的任务很简单，安排每一期节目的内容。

现
在有N个歌手，M种歌曲流派（Rock，Pop之类），每个歌手都有自己擅长的流派领域，这些资料都已整理好。你的工作是，安排尽可能多场的演唱比赛。每
一场比赛所有歌手都必须上场，为了提高收视率，每个人演唱的歌曲类型不能相同，即便一些歌手要被迫选择一些他们不擅长的。同时，为了展现全面性，在不同的
演唱比赛上，每个歌手都会安排不同的歌曲流派。

但是问题是，对于任何一个歌曲流派的歌迷，如果超过K个不擅长的歌手演唱了这种歌曲，他们就会表示不满，比如，发一些宣泄不满的帖子微博，为了表示观点挑起事端等等。你当然不希望这些事情与你的节目有关，在这个前提下，你可以任意安排尽可能多的比赛场次。

链接：点我

online_judge写错了，wa了好多发

N个歌手编号为1~N

M种歌曲，拆点。编号分别为N+1~N+M,N+M+1~N+2*M

源点为0,汇点为N+2*M+1  总共N+2*M+2个点。

N+i 和  N+M+i  之间连权值为K的边，作为限制。

对于i, j  如果i对j擅长，那么连边 i -> N+M+j

否则连边 i -> N+j 这样就可以限制了。

然后是二分答案。

0~(1~N)之间，连边 mid .

(N+M+1~N+2*M)与N+2*M+1之间也连边mid

然后判断mid*N和最大流是不是相当。

然后就是最大流模板题了。

 #include<cstdio>
 #include<iostream>
 #include<algorithm>
 #include<cstring>
 #include<cmath>
 #include<queue>
 #include<map>
 using namespace std;
 #define mOD 1000000007
 #define pb(a) push_back(a)
 const int InF=0x3f3f3f3f;
 const double eps=1e-;
 typedef long long ll;
 #define cl(a) memset(a,0,sizeof(a))
 #define ts printf("*****\n");
 const int MAXN=;
 int n,m,tt,cnt;
 int maze[MAXN][MAXN];
 int gap[MAXN],dis[MAXN],pre[MAXN],cur[MAXN];
 int flow[MAXN][MAXN];//存最大流的容量
 int sap(int start,int end,int nodenum)
 {
 memset(cur,,sizeof(cur));
 memset(dis,,sizeof(dis));
 memset(gap,,sizeof(gap));
 memset(flow,,sizeof(flow));
 int u=pre[start]=start,maxflow=,aug=-;
 gap[]=nodenum;
 while(dis[start]<nodenum)
 {
 loop:
 for(int v=cur[u];v<nodenum;v++)
 if(maze[u][v]-flow[u][v] && dis[u]==dis[v]+)
 {
 if(aug==- || aug>maze[u][v]-flow[u][v])aug=maze[u][v]-flow[u][v];
 pre[v]=u;
 u=cur[u]=v;
 if(v==end)
 {
 maxflow+=aug;
 for(u=pre[u];v!=start;v=u,u=pre[u])
 {
 flow[u][v]+=aug;
 flow[v][u]-=aug;
 }
 aug=-;
 }
 goto loop;
 }
 int mindis=nodenum-;
 for(int v=;v<nodenum;v++)
 if(maze[u][v]-flow[u][v]&&mindis>dis[v])
 {
 cur[u]=v;
 mindis=dis[v];
 }
 if((--gap[dis[u]])==)break;
 gap[dis[u]=mindis+]++;
 u=pre[u];
 }
 return maxflow;
 }
 bool check(int mid)
 {
     for(int i=;i<=n;i++)
     {
         maze[][i]=mid;
     }
     for(int i=n+m+;i<=n+*m;i++)maze[i][n+*m+]=mid;
     if(sap(,n+*m+,n+*m+)==mid*n)   return ;
     else
         return ;
 }
 int used[MAXN][MAXN];
 int main()
 {
     int i,j,k;
     int ca=;
     #ifndef ONLINE_JUDGE
     freopen("1.in","r",stdin);
     #endif
     scanf("%d",&tt);
     int l;
     while(tt--)
     {
         scanf("%d%d%d%d",&n,&m,&l,&k);
         cl(used);
         cl(maze);
         int q,p;
         for(i=;i<=l;i++)
         {
             scanf("%d%d",&q,&p);
             used[q][p]=;
         }

         for(i=n+;i<=n+m;i++)maze[i][i+m]=k;
         for(int i=;i<=n;i++)
             for(int j=;j<=m;j++)
             {
                 if(used[i][j])
                 {
                     maze[i][j+n+m]=;
                 }
                 else maze[i][j+n]=;
             }
         int l=,r=m;
         int ans=;
          while(l<=r)
         {
             int mid=(l+r)>>;
             if(check(mid))
             {
                 ans=mid;
                 l=mid+;
             }
             else r=mid-;
         }
         printf("Case %d: %d\n",ca++,ans);
     }
 }