2 3 2 1 2 1 2 3 1 3 3 1 2 1 2 3 1 1 3 1
Case 1: 1 Case 2: 2
纯求网络最大流,用EK算法就可以实现。
代码:
#include <iostream> #include <cstdio> #include <cstring> #include <queue> using namespace std; #define maxn 20 #define INF 0x3f3f3f3f int ans, s, t, n; int a[maxn], pre[maxn]; int flow[maxn][maxn]; int cap[maxn][maxn]; void Edmonds_Karp() { queue<int> q; memset(flow, 0, sizeof(flow)); ans = 0; while(1) { memset(a, 0, sizeof(a)); a[s] = INF; q.push(s); while(!q.empty()) //bfs找增光路径 { int u = q.front(); q.pop(); for(int v = 1; v <= n; v++) if(!a[v] && cap[u][v] > flow[u][v]) { pre[v] = u; q.push(v); a[v] = min(a[u], cap[u][v]-flow[u][v]); } } if(a[t] == 0) break; for(int u = t; u != s; u = pre[u]) //改进网络流 { flow[pre[u]][u] += a[t]; flow[u][pre[u]] -= a[t]; } ans += a[t]; } } int main() { //freopen("hdu_3549.txt", "r", stdin); int T, m, u, v, c; scanf("%d", &T); for(int cas = 1; cas <= T; cas++) { scanf("%d%d", &n, &m); memset(cap, 0, sizeof(cap)); while(m--) { scanf("%d%d%d", &u, &v, &c); cap[u][v] += c; } s = 1, t = n; Edmonds_Karp(); printf("Case %d: %d\n", cas, ans); } return 0; }
HDU 3549 Flow Problem 网络最大流问题 Edmonds_Karp算法
原文地址:http://blog.csdn.net/u011439796/article/details/39485465