LightOJ 1286 - Space Shuttle Experiments 最大流

题目：http://www.lightoj.com/volume_showproblem.php?problem=1286

题意：有n个实验，做每个实验可以得到对应的收入，每个实验需要一些仪器，搬运这些仪器需要对应的花费（一个仪器被好几个实验使用，也只搬运一次），求最后的盈利，即收入减去搬运仪器的费用。输入格式：首先一个t，代表测试数据组数，然后是n m，实验数和仪器数，之后是n个数，代表做相应实验的收入，之后m个数，代表搬运相应仪器的花费。最后n行，每行开始一个数q，然后q个数，代表做第i个实验需要的仪器编号。

思路：首先，我们在源点和实验之间连边，容量为相应的收入，在仪器和汇点之间连边，容量为搬运的费用，然后在相应的实验和仪器之间连边，容量为实验收入（极大值也可以，只要保证大于等于实验收入即可），此时最大流的意义就是搬运仪器能够付出的花费。什么叫能够付出的花费呢？就是如果用到某些仪器可以盈利，那就付出搬运这些仪器的成本，如果不能盈利，最多就付出相应实验的全部收入，不会再多（保证不亏本啊）。对应盈利的项目，收入减去成本，自然就是利润啦，对于不盈利的项目，收入减去成本等于0（因为之前说过付出全部收入），其意义就是不做这个实验（亏本做毛线啊）。因此，结果就是全部收入减去最大流

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <queue>
using namespace std;

const int N = 210;
const int INF = 0x3f3f3f3f;
struct edge
{
    int to, cap, next;
}g[N*N*2];
int head[N], iter[N], level[N];
int n, m, cnt, _case = 0;
void add_edge(int v, int u, int cap)
{
    g[cnt].to = u, g[cnt].cap = cap, g[cnt].next = head[v], head[v] = cnt++;
    g[cnt].to = v, g[cnt].cap = 0, g[cnt].next = head[u], head[u] = cnt++;
}
bool bfs(int s, int t)
{
    memset(level, -1, sizeof level);
    level[s] = 0;
    queue<int> que;
    que.push(s);
    while(! que.empty())
    {
        int v = que.front(); que.pop();
        for(int i = head[v]; i != -1; i = g[i].next)
        {
            int u = g[i].to;
            if(g[i].cap > 0 && level[u] < 0)
            {
                level[u] = level[v] + 1;
                que.push(u);
            }
        }
    }
    return level[t] == -1;
}
int dfs(int v, int t, int f)
{
    if(v == t) return f;
    for(int &i = iter[v]; i != -1; i = g[i].next)
    {
        int u = g[i].to;
        if(g[i].cap > 0 && level[v] < level[u])
        {
            int d = dfs(u, t, min(g[i].cap, f));
            if(d > 0)
            {
                g[i].cap -= d, g[i^1].cap += d;
                return d;
            }
        }
    }
    return 0;
}
int dinic(int s, int t)
{
    int flow = 0, f;
    while(true)
    {
        if(bfs(s, t)) return flow;
        memcpy(iter, head, sizeof head);
        while(f = dfs(s, t, INF),f > 0)
            flow += f;
    }
}
int main()
{
    int t, q, a;
    scanf("%d", &t);
    while(t--)
    {
        cnt = 0;
        memset(head, -1, sizeof head);
        scanf("%d%d", &n, &m);
        int t1[N], t2[N];
        int sum = 0; //sum统计总收入
        for(int i = 1; i <= n; i++)
        {
            scanf("%d", t1 + i);
            sum += t1[i];
            add_edge(0, i, t1[i]); //源点和实验连边，容量为收入
        }
        for(int i = 1; i <= m; i++)
        {
            scanf("%d", t2 + i);
            add_edge(n + i, n + m + 1, t2[i]); //仪器和汇点连边，容量为花费
        }
        for(int i = 1; i <= n; i++)
        {
            scanf("%d", &q);
            for(int j = 1; j <= q; j++)
            {
                scanf("%d", &a);
                add_edge(i, n + a, t1[i]); //实验和仪器连边
            }
        }
        printf("Case %d: %d\n", ++_case, sum - dinic(0, n + m + 1));
    }

    return 0;
}