poj 1149 PIGS 最大流 太神奇的建图方式！！！

有M个猪圈，每个猪圈里初始时有若干头猪。一开始所有猪圈都是关闭的。依次来了N个顾客，每个顾客分别会打开指定的几个猪圈，从中买若干头猪。每个顾客分别都有他能够买的数量的上限。每个顾客走后，他打开的那些猪圈中的猪，都可以被任意地调换到其它开着的猪圈里，然后所有猪圈重新关上。问总共最多能卖出多少头猪。（1 <= N <= 100, 1 <= M <= 1000）

这道题是通过特殊的化简规律，把一个复杂网络流化简成了特别简单的形式！！

资料：https://wenku.baidu.com/view/0ad00abec77da26925c5b01c.html  

这分资料特别棒，详细的解释了网络流的化简问题。

1. #include <cstdio>  
2. #include <cstring>  
3. #include <iostream>  
4. #include <algorithm>  
5. #include <queue>  
6. #include <vector>  
7. using namespace std;  
8.   
9. const int maxn = 1000 + 50;  
10. const int INF = 1e9;  
11.   
12. int n, m;  
13.   
14. struct Edge{  
15.     int from, to, cap, flow;  
16. };  
17.   
18. struct Dinic{  
19.     int n, m, s, t;  
20.     vector<Edge> edges;  
21.     vector<int> G[maxn];  
22.     bool vis[maxn];  
23.     int d[maxn];  
24.     int cur[maxn];  
25.   
26.     void init(int n)  
27.     {  
28.         this->n = n;  
29.         for(int i=0; i<=n; ++i) G[i].clear();  
30.         edges.clear();  
31.     }  
32.   
33.     void AddEdge(int from, int to, int cap)  
34.     {  
35.         edges.push_back((Edge){  
36.                 from, to, cap, 0  
37.                 });  
38.         edges.push_back((Edge){  
39.                 to, from, 0, 0  
40.                 });  
41.         m = edges.size();  
42.         G[from].push_back(m-2);  
43.         G[to].push_back(m-1);  
44.     }  
45.   
46.     bool BFS(){  
47.         memset(vis, 0, sizeof vis );  
48.         queue<int> Q;  
49.         Q.push(s);  
50.         vis[s] = 1;  
51.         d[s] = 0;  
52.         while(!Q.empty()){  
53.             int x = Q.front(); Q.pop();  
54.             for(int i=0; i<G[x].size(); ++i){  
55.                 Edge& e = edges[G[x][i]];  
56.                 if(!vis[e.to] && e.cap > e.flow){  
57.                     vis[e.to] = 1;  
58.                     d[e.to] = d[x] + 1;  
59.                     Q.push(e.to);  
60.                 }  
61.             }  
62.         }  
63.         return vis[t];  
64.     }  
65.   
66.     int DFS(int x, int a){  
67.         if(x==t || a==0) return a;  
68.         int flow = 0, f;  
69.         for(int& i=cur[x]; i<G[x].size(); ++i){  
70.             Edge& e = edges[G[x][i]];  
71.             if(d[x]+1==d[e.to] && (f=DFS(e.to, min(a, e.cap-e.flow)))>0){  
72.                 e.flow += f;  
73.                 edges[G[x][i]^1].flow -= f;  
74.                 flow += f;  
75.                 a -= f;  
76.                 if(a==0) break;  
77.             }  
78.         }  
79.         return flow;  
80.     }  
81.   
82.     int MaxFlow(int s, int t){  
83.         this->s = s; this->t = t;  
84.         int flow = 0;  
85.         while(BFS()){  
86.             memset(cur, 0, sizeof cur );  
87.             flow += DFS(s, INF);  
88.         }  
89.         return flow;  
90.     }  
91. };  
92.   
93. Dinic gao;  
94. int pig[maxn], pre[maxn];  
95.   
96. void fk()  
97. {  
98.     int A, B, K;  
99.     int s = 0, t = n + 1;  
100.     gao.init(t+2);    
101.     memset(pre, 0, sizeof pre );  
102.     for(int i=1; i<=m; ++i){  
103.         scanf("%d", &pig[i]);  
104.     }  
105.     for(int i=1; i<=n; ++i){  
106.         scanf("%d", &A);  
107.         for(int j=1; j<=A; ++j){  
108.             scanf("%d", &K);  
109.             if(!pre[K]){  
110.                 gao.AddEdge(s, i, pig[K]);  
111.             }else {  
112.                 gao.AddEdge(pre[K], i, INF);  
113.             }  
114.             pre[K] = i;  
115.         }  
116.         scanf("%d", &B);  
117.         gao.AddEdge(i, t, B);  
118.     }  
119.     int ans = gao.MaxFlow(s, t);  
120.     printf("%d\n", ans);  
121. }  
122.   
123. int main()  
124. {  
125. //  freopen("in.txt", "r",stdin);  
126.     while(~scanf("%d%d", &m, &n)){  
127.         fk();  
128.     }  
129.     return 0;  
130. }