http://poj.org/problem?id=3281

我的最大流的第二题，，，话说刚学网络流，，用了最常用的算法E_k算法还调试了一会，，，，杯具啊，，看了网上大牛们的解题报告，用的都是Dinci和ISPA，，，很是神奇啊，，看来是得学学了，，，这道题难在建图上，，，为了保证每种食物，和饮料只对应一头牛，用了拆点，，，建图顺序为源点-食物-牛-牛-饮料-汇点，，每个弧的容量为1，做完这题我明白了为什么最大二分匹配可以用网络流来做了，，嘎嘎，，，

#include<iostream>
#include<cstdio>
#include<limits.h>
#include<algorithm>
#include<queue>
#define M  405
#include<string.h>
using namespace std;
int map[M][M];
int N,F,D;
void E_k(int s,int t)
{  int flow[M][M],f[M],pre[M],flowmax=0;
   queue<int> Q;
    memset(flow,0,sizeof(flow));
     while(1)
     {    Q.push(s);
      memset(f,0,sizeof(f));
         f[s]=INT_MAX;
         while(!Q.empty())
         {  int u=Q.front();
             Q.pop();
             for(int i=s;i<=t;++i)
                if(!f[i]&&map[u][i]>flow[u][i])
                  {  Q.push(i);
                    f[i]=min(f[u],map[u][i]-flow[u][i]); 
                    pre[i]=u;
                   }
           }
           if(!f[t]) break;
           for(int i=t;i!=s;i=pre[i])
             { flow[pre[i]][i]+=f[t];
               flow[i][pre[i]]-=f[t];
              }
             flowmax+=f[t];
    }
             printf("%d\n",flowmax);
}
int main()
{   while(~scanf("%d%d%d",&N,&F,&D))
    {  memset(map,0,sizeof(map));
        int a,b;
       for(int i=1;i<=N;++i)
       {  scanf("%d%d",&a,&b);
          map[i][N+i]=1;
          while(a--)
          { int c;
            scanf("%d",&c);
            map[2*N+c][i]=1;
          }
          while(b--)
          { int c;
             scanf("%d",&c);
             map[N+i][2*N+F+c]=1;
           }
          }
          for(int i=2*N+1;i<=2*N+F;++i)
            map[0][i]=1;
            for(int i=2*N+F+1;i<=2*N+F+D;++i)
              map[i][2*N+F+D+1]=1;
              E_k(0,2*N+F+D+1);
     }return 0;
}