usaco 4.2 The Perfect Stall 二分图最大匹配

 这里采用网络流方法解决，建立虚拟源点src,从src到每头牛添边，建立虚拟汇点des,从每个摊点向des添边，每条边的权值都是1，求从src到des的最大流即可

 

 

1. #include <iostream>
2. #include <algorithm>
3. #include <string>
4. #include <map>
5. using namespace std;
7. /*
8. PROG: stall4
9. LANG: C++
10. ID: heben991
11. */
14. const int N = 420, INF = 2000000000, MaxFlow = 10000000;
16. int g[N][N], flow[N][N];
17. int n, cow, stall;
18. int inc[N], q[N], pre[N];
21. bool bfs(int src, int des)
22. {
23.     int fp, rp, i, j, u;
25.     for(i = 0; i <= n; ++i) pre[i]=-1;
26.     inc[src] = INF;
27.     q[fp=rp=0] = src;
28.     while(fp<=rp)
29.     {
30.         u = q[fp++];
31.         for(i = 0; i <= n; ++i)
32.         if(pre[i]==-1 && flow[u][i] < g[u][i])
33.         {
34.             inc[i] = min(inc[u], g[u][i]-flow[u][i]);
35.             pre[i] = u;
36.             if(i==des) return true;
37.             q[++rp] = i;
38.         }
39.     }
40.     return false;
41. }
44. int Edmonds_Karp(int src, int des)
45. {
46.     int i, j, maxflow=0;
48.     while(bfs(src,des))
49.     {
50.         for(i = des; i != src; i = pre[i])
51.         {
52.             j = pre[i];
53.             flow[j][i] += inc[i];
54.             flow[i][j] = -flow[j][i];
55.         }
56.         maxflow += inc[des];
57.     }
58.     return maxflow;
59. }
61. int main()
62. {
63.     int i, j, k, x, y;
65.     freopen("stall4.in", "r", stdin);
66.     freopen("stall4.out","w",stdout);
68.     scanf("%d%d", &cow, &stall);
69.     for(i = 1; i <= cow; ++i)
70.     {
71.         scanf("%d", &x);
72.         for(j = 1; j <= x; ++j)
73.         {
74.             scanf("%d", &y);
75.             g[i][y+cow] = 1;
76.         }
77.     }
79.     for(i = 1; i <= cow; ++i) g[0][i] = 1;
80.     n = cow+stall+1;
81.     for(i = cow+1; i < n; ++i) g[i][n] = 1;
83.     printf("%d/n", Edmonds_Karp(0,n) );
85.     return 0;
86. }