网络流经典/Dinic（Dining POJ - 3281）

题目链接
从这个题，我知道了一个真理：网络流问题就是一个花式建图问题，考的不是你板子敲得有多熟练。而是你的建图有多骚。
对于这个题，我的理解：每个food先去匹配一头cow，cow带着food去匹配drink。。。。

food -> cow ->drink

这样的的话，从超源点s出发到food，cow和喜欢的food之间有一个food->cow的容量为1的边， 而cow-> drink同理。然后到超汇点t。可是这样有一个显而易见的问题：每个cow有一堆喜欢的drink和food。不排除一个cow匹配一堆food，drink的可能。所以我们拆cow点：把cow拆成两个点中间建立一个容量为1的边。问题得以解决。
下面是ac代码：

#include <iostream>
#include <cstdio>
#include <cstring>
#include <string>
#include <vector>
#include <cmath>
#include <algorithm>
#include <map>
#include <queue>
#define ll long long
using namespace std;
const int NN = 1e5+5;
const int INF = 0x3f3f3f3f;
typedef unsigned long long ull;
int he[NN], ne[NN], ver[NN], e[NN], l[NN];
int tot;
void init()
{
    memset(he, 0, sizeof(he));
    tot = 1;
}
void add(int x, int y, int w)
{
    e[++tot] = w;
    ver[tot] = y;
    ne[tot] = he[x];
    he[x] = tot;
    e[++tot] = 0;
    ver[tot] = x;
    ne[tot] = he[y];
    he[y] = tot;
}
bool bfs(int s, int t)
{
    memset(l, 0, sizeof(l));
    queue<int> q;
    q.push(s);
    l[s] = 1;
    while(q.size())
    {
        int x = q.front();
        q.pop();
        if (x == t) return 1;
        for (int i = he[x]; i; i = ne[i])
        {
            int y = ver[i];
            if (!l[y] && e[i])
            {
                l[y] = l[x] + 1;
                q.push(y);
            }
        }
    }
    return 0;
}
int dfs(int x, int MaxFlow, int t)
{
    if (x == t) return MaxFlow;
    int uflow = 0;
    for (int i = he[x]; i; i = ne[i])
    {
        int y = ver[i];
        if (l[y] == l[x] + 1 && e[i])
        {
            int flow = min(e[i], MaxFlow - uflow);
            flow = dfs(y, flow, t);
            e[i] -= flow;
            e[i^1] += flow;
            uflow += flow;
            if (uflow == MaxFlow)
                break;
        }
    }
    if (uflow == 0)
        l[x] = 0;
    return uflow;
}
int Dinic(int s, int t)
{
    int MaxFlow = 0;
    while(bfs(s, t))
        MaxFlow += dfs(s, INF, t);
    return MaxFlow;
}
int main()
{
    int n, f, d;
    while(~scanf("%d%d%d", &n, &f, &d))
    {
        int all = 2 * n + f + d + 2;
        int s = 2 * n + f + d, t = s + 1;
        init();
        for (int i = 0; i < f; i++)
            add(s, 2 * n + i, 1);
        for (int i = 0; i < d; i++)
            add(2 * n + f + i, t, 1);
        for (int i = 0; i < n; i++)
            add (i, n+i, 1);
        for (int i = 0; i < n; i++)
        {
            int _f, _d;
            scanf("%d%d", &_f, &_d);
            for (int j = 0; j < _f; j++)
            {
                int fi;
                scanf("%d", &fi);
                add(2*n + fi - 1, i, 1);
            }
            for (int j = 0; j < _d; j++)
            {
                int di;
                scanf("%d", &di);
                add(n+i, 2*n+f+di-1, 1);
            }
        }
        int ans = Dinic(s,t);
        printf("%d\n", ans);
    }
}