UVA10779 Collectors Problem(最大流)

白皮书上的练习题0 0
大意：有$T(T \leq 20)$组数据。Bob在与他的$n-1(2 \leq n \leq 10)$个朋友交换糖纸，一共有$m(5 \leq m \leq 25)$种糖纸。每个朋友手里有一些糖纸，可以一对一与Bob交换自己没有的糖纸，但是每个朋友只会交出自己的重复的糖纸。朋友之间不互相交换糖纸。问Bob最多能收集到多少种糖纸。
还是按照书上的方式建模吧，自己没有经验- -。
把Bob看作源点，把他与每种糖纸之间连一条边，容量为他有的该种糖纸数量。虚拟一个汇点，每一种糖纸与汇点连一条容量为1的边，因为每种糖纸只被计算一次。对于第i朋友，假设第k种糖纸他有p张，那么i到k之间连一条容量为p-1的边(因为自己要留一张)，如果没有第k种糖纸，那么k到i之间连一条容量为1的边(每种糖纸最多接受一张)。好了，在构建好的图中跑个最大流，就是Bob能拥有的最多种数。

#include<cstdio>
#include<cstring>
using namespace std;
int c[50][50],t,k,x,a[50][60],n,m;
int flow,d[50],vd[50];
inline int Min(int a,int b)
{return a<b?a:b;}
int aug(int i,int augco)
{
    int j,augc = augco,mind = t-1,delta;
    if(i == t) return augco;
    for(j = 1; j <= t; j++)
    if(c[i][j])
    {
        if(d[i] == d[j]+1)
        {
            delta = Min(augc,c[i][j]);
            delta = aug(j,delta);
            c[i][j] -= delta;
            c[j][i] += delta;
            augc -= delta;
            if(d[1] >= t) return augco-augc;
            if(augc == 0) break;
        }
        if(mind > d[j]) mind = d[j];
    }
    if(augc == augco)
    {
        vd[d[i]]--;
        if(vd[d[i]] == 0) d[1] = n;
        d[i] = mind+1;
        vd[d[i]]++;
    }
    return augco - augc;
}
int sap()
{
    t = n+m+1;
    flow = 0;
    memset(d,0,sizeof d);
    memset(vd,0,sizeof vd);
    vd[0] = t;
    while(d[1] < t)
        flow += aug(1,0x3f3f3f3f);
    return flow;
}
void init()
{
    memset(c,0,sizeof c);
    memset(a,0,sizeof a);
}
int main()
{
    int T;
    scanf("%d",&T);
    for(int o = 1; o <= T; o++)
    {
        init();
        scanf("%d%d",&n,&m);
        for(int i = 1; i <= n; i++)
        {
            scanf("%d",&k);
            for(int j = 1; j <= k; j++)
            {
                scanf("%d",&x);
                a[i][x]++;
            }
        }
        for(int i = 1; i <= m; i++)
            if(a[1][i]) c[1][i+n] = a[1][i];
        for(int i = 1; i <= m; i++)
            c[i+n][n+m+1] = 1;
        for(int i = 2; i <= n; i++)
        {
            for(int j = 1; j <= m; j++)
                if(a[i][j] > 1)
                    c[i][j+n] = a[i][j]-1;
                else if(!a[i][j])
                    c[j+n][i] = 1;
        }
        printf("Case #%d: %d\n",o,sap());
    }
}