【最大流】bzoj1711: [Usaco2007 Open]Dining吃饭

正在网络流入门（原来这种题用网络流做）

Description

农夫JOHN为牛们做了很好的食品,但是牛吃饭很挑食. 每一头牛只喜欢吃一些食品和饮料而别的一概不吃.虽然他不一定能把所有牛喂饱,他还是想让尽可能多的牛吃到他们喜欢的食品和饮料. 农夫JOHN做了F (1 <= F <= 100) 种食品并准备了D (1 <= D <= 100) 种饮料. 他的N ( 1 <= N <= 100)头牛都以决定了是否愿意吃某种食物和喝某种饮料. 农夫JOHN想给每一头牛一种食品和一种饮料,使得尽可能多的牛得到喜欢的食物和饮料. 每一件食物和饮料只能由一头牛来用. 例如如果食物2被一头牛吃掉了,没有别的牛能吃食物2.

Input

* 第一行: 三个数: N, F, 和 D

* 第2..N+1行: 每一行由两个数开始F_i 和D_i, 分别是第i 头牛可以吃的食品数和可以喝的饮料数.下F_i个整数是第i头牛可以吃的食品号,再下面的D_i个整数是第i头牛可以喝的饮料号码.

Output

* 第一行: 一个整数,最多可以喂饱的牛数.

---

题目分析

首先将源点连向所有食物；所有饮料连向汇点。接下来是中间奶牛的部分，容易发现如果简单地食物-奶牛-饮料一连，将会导致一头奶牛可能吃了很多饮料食物。这里有一种巧妙的处理方法：将一头奶牛拆成两个点，中间连一条容量为1的边，意味着一头牛只能占有一个食物/饮料。

话说这题在bzoj为什么会莫名其妙TLE啊……

 #include<bits/stdc++.h>
 const int maxn = ;
 const int maxm = ;
 const int INF = 2e9;
 
 struct Edge
 {
     int u,v,f,c;
     Edge(int a=, int b=, int c=, int d=):u(a),v(b),f(c),c(d) {}
 }edges[maxm];
 int n,f,d,S,T,lv[maxn];
 int edgeTot,head[maxn],nxt[maxm];
 
 int read()
 {
     char ch = getchar();
     int num = , fl = ;
     for (; !isdigit(ch); ch = getchar())
         if (ch=='-') fl = -;
     for (; isdigit(ch); ch = getchar())
         num = (num<<)+(num<<)+ch-;
     return num*fl;
 }
 void addedge(int u, int v, int c)
 {
     edges[edgeTot] = Edge(u, v, , c), nxt[edgeTot] = head[u], head[u] = edgeTot++;
     edges[edgeTot] = Edge(v, u, , ), nxt[edgeTot] = head[v], head[v] = edgeTot++;
 }
 bool buildLevel()
 {
     memset(lv, , sizeof lv);
     std::queue<int> q;
     q.push(S), lv[S] = ;
     for (int tmp; q.size(); )
     {
         tmp = q.front(), q.pop();
         for (int i=head[tmp]; i!=-; i=nxt[i])
         {
             int v = edges[i].v;
             if (!lv[v]&&edges[i].f < edges[i].c){
                 lv[v] = lv[tmp]+, q.push(v);
                 if (v==T) return true;
             }
         }
     }
     return false;
 }
 int fndPath(int x, int lim)
 {
     if (x==T) return lim;
     for (int i=head[x]; i!=-; i=nxt[i])
     {
         int v = edges[i].v, val;
         if (lv[x]+==lv[v]&&edges[i].f < edges[i].c){
             if ((val = fndPath(v, std::min(lim, edges[i].c-edges[i].f)))){
                 edges[i].f += val, edges[i^].f -= val;
                 return val;
             }else lv[v] = -;
         }
     }
     return ;
 }
 int dinic()
 {
     int ret = , val;
     while (buildLevel())
         while ((val = fndPath(S, INF))) ret += val;
     return ret;
 }
 int main()
 {
     memset(head, -, sizeof head);
     n = read(), f = read(), d = read();
     S = , T = n*+f+d+;
     for (int i=; i<=n; i++)
     {
         int k1 = read(), k2 = read();
         for (; k1; --k1)
             addedge(read(), f+i, );
         for (; k2; --k2)
             addedge(f+i+n, read()+n*+f, );
         addedge(f+i, f+i+n, );
     }
     for (int i=; i<=f; i++) addedge(S, i, );
     for (int i=; i<=d; i++) addedge(f+*n+i, T, );
     printf("%d\n",dinic());
     return ;
 }

END