poj 2289 Jamie's Contact Groups（二分答案+网络流判定）

【题目大意】：n个人，告诉你他可以从属于那些小组，现在问你将n个人分成m组，最大的组最小的人数是多少。

【解题思路】：最大最小值...想到的是二分出上界，符合条件则缩小上界，反之则增大。如何判断，可以用网络流来判定，构图是由s出发，连边到人容量为1，每个人和小组之间的对应关系连边，容量为1，小组和汇点t连边，容量为二分出来的上界...然后跑一个网络流判断汇点是否是n就可以了....

热身赛，队友不在，不会图论，硬着头皮扛，-_-!!....用的还是队友的模版......除掉模版之后...没了....

【代码】：

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <vector>
#include <queue>
#include <cmath>
#include <string>
#include <cctype>
#include <map>
#include <iomanip>
                   
using namespace std;
                   
#define eps 1e-8
#define pi acos(-1.0)
#define inf 1<<30
#define linf 1LL<<60
#define pb push_back
#define lc(x) (x << 1)
#define rc(x) (x << 1 | 1)
#define lowbit(x) (x & (-x))
#define ll long long
#define maxn 2000
#define maxm 1000000

struct per{
     int st[600],cnt;
}per[2000];

int s,t;
int n,m;
int dist[maxn],low[maxn],tot,eh[maxn],pre[maxn],cnt[maxn],cur[maxn];
int maz[maxn][maxn],mapp[maxn][maxn];
 
struct Edge {
    int u,v,cap,flow,next;
}et[maxm];
 
void init() {
    tot=0;
    memset(eh,-1,sizeof(eh));
}
 
void add(int u,int v,int cap,int flow) {
    Edge E={u,v,cap,flow,eh[u]};
    et[tot]=E;
    eh[u]=tot++;
}
 
void addedge(int u,int v,int cap) {
    add(u,v,cap,0),add(v,u,0,0);
}
 
int isap(int s,int t,int n) {
    int u,v,now;
    memset(dist,0,sizeof(dist));
    memset(low,0,sizeof(low));
    for (u=0; u<=n; u++) cur[u]=eh[u];
    int maxflow=0;
    u=s;
    low[s]=inf,cnt[0]=n;
    while (dist[s]<n) {
        for (now=cur[u]; now!=-1; now=et[now].next)
            if (et[now].cap-et[now].flow && dist[u]==dist[v=et[now].v]+1) break;
        if (now!=-1) {
            cur[u]=pre[v]=now;
            low[v]=min(low[u],et[now].cap-et[now].flow);
            u=v;
            if (u==t) {
                for (; u!=s; u=et[pre[u]].u) {
                    et[pre[u]].flow+=low[t];
                    et[pre[u]^1].flow-=low[t];
                }
                low[s]=inf;
                maxflow+=low[t];
            }
        } 
        else {
            if (--cnt[dist[u]]==0) break;
            dist[u]=n;
            cur[u]=eh[u];
            for (now=eh[u]; now!=-1; now=et[now].next)
                if (et[now].cap-et[now].flow && dist[u]>dist[et[now].v]+1) 
                    dist[u]=dist[et[now].v]+1;
             cnt[dist[u]]++;
             if(u!=s) u=et[pre[u]].u;
         }
     }
     return maxflow;
 }


bool check(int cap) {    
    s=0;
    t=n+m+1;
    init();
    for (int i=1; i<=n; i++) {
        addedge(s,i,1);
        for (int j=0; j<per[i].cnt; j++)
            addedge(i,per[i].st[j]+n+1,1);
    }
    for (int i=1; i<=m; i++) addedge(i+n,t,cap);
    if (isap(s,t,t-s+1)==n) return true;
    else return false;
}
 
int main(){
    char ch;
    string name;
    while(~scanf("%d%d",&n,&m)){
        if (n==0 && m==0) break;
        for (int i=1; i<=n; i++) {
            cin >> name;
            per[i].cnt=0;
            ch=getchar();
            while (ch!='\n'){
                scanf("%d",&per[i].st[per[i].cnt]);
                per[i].cnt++;
                ch=getchar();
            }
        }
        int l=1,r=n,ans=-1;
        while (l<=r) {
            int mid=(l+r)/2;
            if (check(mid)) {
                ans=mid;
                r=mid-1;
            }
            else l=mid+1;
        }
        printf("%d\n",ans);
    }
    return 0;
}