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和最大流是不是相当。

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

  1 #include<cstdio>

  2 #include<iostream>

  3 #include<algorithm>

  4 #include<cstring>

  5 #include<cmath>

  6 #include<queue>

  7 #include<map>

  8 using namespace std;

  9 #define mOD 1000000007

 10 #define pb(a) push_back(a)

 11 const int InF=0x3f3f3f3f;

 12 const double eps=1e-5;

 13 typedef long long ll;

 14 #define cl(a) memset(a,0,sizeof(a))

 15 #define ts printf("*****\n");

 16 const int MAXN=1010;

 17 int n,m,tt,cnt;

 18 int maze[MAXN][MAXN];

 19 int gap[MAXN],dis[MAXN],pre[MAXN],cur[MAXN];

 20 int flow[MAXN][MAXN];//存最大流的容量

 21 int sap(int start,int end,int nodenum)

 22 {

 23 memset(cur,0,sizeof(cur));

 24 memset(dis,0,sizeof(dis));

 25 memset(gap,0,sizeof(gap));

 26 memset(flow,0,sizeof(flow));

 27 int u=pre[start]=start,maxflow=0,aug=-1;

 28 gap[0]=nodenum;

 29 while(dis[start]<nodenum)

 30 {

 31 loop:

 32 for(int v=cur[u];v<nodenum;v++)

 33 if(maze[u][v]-flow[u][v] && dis[u]==dis[v]+1)

 34 {

 35 if(aug==-1 || aug>maze[u][v]-flow[u][v])aug=maze[u][v]-flow[u][v];

 36 pre[v]=u;

 37 u=cur[u]=v;

 38 if(v==end)

 39 {

 40 maxflow+=aug;

 41 for(u=pre[u];v!=start;v=u,u=pre[u])

 42 {

 43 flow[u][v]+=aug;

 44 flow[v][u]-=aug;

 45 }

 46 aug=-1;

 47 }

 48 goto loop;

 49 }

 50 int mindis=nodenum-1;

 51 for(int v=0;v<nodenum;v++)

 52 if(maze[u][v]-flow[u][v]&&mindis>dis[v])

 53 {

 54 cur[u]=v;

 55 mindis=dis[v];

 56 }

 57 if((--gap[dis[u]])==0)break;

 58 gap[dis[u]=mindis+1]++;

 59 u=pre[u];

 60 }

 61 return maxflow;

 62 }

 63 bool check(int mid)

 64 {

 65     for(int i=1;i<=n;i++)

 66     {

 67         maze[0][i]=mid;

 68     }

 69     for(int i=n+m+1;i<=n+2*m;i++)maze[i][n+2*m+1]=mid;

 70     if(sap(0,n+2*m+1,n+2*m+2)==mid*n)   return 1;

 71     else

 72         return 0;

 73 }

 74 int used[MAXN][MAXN];

 75 int main()

 76 {

 77     int i,j,k;

 78     int ca=1;

 79     #ifndef ONLINE_JUDGE

 80     freopen("1.in","r",stdin);

 81     #endif

 82     scanf("%d",&tt);

 83     int l;

 84     while(tt--)

 85     {

 86         scanf("%d%d%d%d",&n,&m,&l,&k);

 87         cl(used);

 88         cl(maze);

 89         int q,p;

 90         for(i=1;i<=l;i++)

 91         {

 92             scanf("%d%d",&q,&p);

 93             used[q][p]=1;

 94         }

 95 

 96         for(i=n+1;i<=n+m;i++)maze[i][i+m]=k;

 97         for(int i=1;i<=n;i++)

 98             for(int j=1;j<=m;j++)

 99             {

100                 if(used[i][j])

101                 {

102                     maze[i][j+n+m]=1;

103                 }

104                 else maze[i][j+n]=1;

105             }

106         int l=0,r=m;

107         int ans=0;

108          while(l<=r)

109         {

110             int mid=(l+r)>>1;

111             if(check(mid))

112             {

113                 ans=mid;

114                 l=mid+1;

115             }

116             else r=mid-1;

117         }

118         printf("Case %d: %d\n",ca++,ans);

119     }

120 }