poj The Perfect Stall

    去这里找我拉的题目：http://blog.csdn.net/zhongshijunacm/article/details/39925067

题目:

   给出N头牛，M给位置。每个牛最一些位置有好感，如果对该位置有好感，则他愿意产奶，否则不产。现在给你N个牛对哪些位置有好感，让你安排使得最多的牛产奶。

 

算法分析：

   题目很裸的二分图。一开始用向前星存图，wrong到无语T_T。。。。

   后来一怒改成矩阵过了。再说用网络流的经历吧。。

   可以容易的想到是要拆点的，但是一开始的时候，脑残了。。。一直只对牛或者位置拆点。其实这两个都是只能使用一次的所以都要拆点。悲剧啊。好几个小时没了。好像EK网络流不用拆点，我研究研究，也把代码贴上。

/*

  去你大爷！！！！！！！！！！！！
  去去屎啊！！！！！！1

*/


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

const int INF = 1 << 20;
const int MAXN = 1000 + 10;
struct Edge{
   int from,to,cap,flow;
   Edge(){};
   Edge(int _from,int _to,int _cap,int _flow)
      :from(_from),to(_to),cap(_cap),flow(_flow){};
};
vector<Edge> edges;
vector<int> G[MAXN];
int cur[MAXN],d[MAXN];
bool vst[MAXN];
int src,sink,N,M;

void init(){
   src = N + N + M + M + 1; sink = src + 1;
   for(int i = 0;i <= sink;++i)
     G[i].clear();
   edges.clear();
}

void addEdge(int from,int to,int cap){
    edges.push_back(Edge(from,to,cap,0));
    edges.push_back(Edge(to,from,0,0));
    int sz = edges.size();
    G[from].push_back(sz - 2);
    G[to].push_back(sz - 1);
}

bool BFS(){
   queue<int> Q;
   Q.push(src);
   d[src] = 0;
   memset(vst,0,sizeof(vst));
   vst[src] = 1;

   while(!Q.empty()){
       int u = Q.front(); Q.pop();
       for(int i = 0;i < G[u].size();++i){
           Edge& e = edges[G[u][i]];
           if(!vst[e.to] && e.cap > e.flow){
               vst[e.to] = 1;
               d[e.to] = d[u] + 1;
               Q.push(e.to);
           }
       }

   }
   return vst[sink];
}

int DFS(int x,int a){
   if(x == sink || a == 0)
     return a;
   int flow = 0,f;
   for(int& i = cur[x];i < G[x].size();++i){
      Edge& e = edges[G[x][i]];
      if(d[e.to] == d[x] + 1 && (f = DFS(e.to,min(a,e.cap - e.flow))) > 0){
          e.flow += f;
          edges[G[x][i]^1].flow -= f;
          flow += f;
          a -= f;
          if(a == 0) break;
      }
   }
   return flow;
}

int maxFlow(){
    int flow = 0;
    while(BFS()){
        memset(cur,0,sizeof(cur));
        flow += DFS(src,INF);
    }
    return flow;
}

int main()
{
   while(~scanf("%d%d",&N,&M)){
        init();
        for(int i = 1;i <= N;++i){
            addEdge(src,i,INF);
        }

        for(int i = 1;i <= N;++i){ //前拆点
            addEdge(i,i + N,1);
        }

        int cas,x,t = N + N;
        for(int i = 1;i <= N;++i){
            scanf("%d",&cas);
            while(cas-- > 0){
                scanf("%d",&x);
                addEdge(i + N,x + t,1);
            }
        }

        for(int i = 1;i <= M;++i){     //后拆点
           addEdge(i + t,i + t + M,1);
        }

        t += M;
        for(int i = 1;i <= M;++i){
            addEdge(t + i,sink,INF);
        }

        printf("%d\n",maxFlow());
   }
   return 0;
}