usaco 4.2 Drainage Ditches 网络流

 很明显的网络流问题，点数和边数都是200，采用Edmonds_karp的复杂度上限是O(VE^2)，实际上一般情况下要比此上限小得多，所以N=200时也能很快解决

 

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