poj1273

最新推荐文章于 2021-03-07 07:25:47 发布

原创最新推荐文章于 2021-03-07 07:25:47 发布 · 579 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#算法 #网络 #c

acm--网络流问题专栏收录该内容

1 篇文章

订阅专栏

本文提供了一个Edmonds-Karp算法的实现模板，该算法用于解决网络流问题，并特别注意处理两点间可能存在多条边的情况。

网络流问题！！！

完全的edmonds karp 算法模板！！！！

需要注意的是要他可能在两点之间挖多条路！！！

#include<cstdio>
#include<cstring>
#include<queue>
using namespace std;
const int INF = 1000000000;
const int MAXN = 1000;

int n, m, s, t;
int a[MAXN], p[MAXN], cap[MAXN][MAXN], flow[MAXN][MAXN];

int main() {
while(scanf("%d%d", &n, &m)!=EOF){
memset(cap, 0, sizeof(cap));
for(int i = 0; i < n; i++) {
    int u, v, c;
    scanf("%d%d%d", &u, &v, &c);
    cap[u][v] += c;//关键！！！
}
s = 1;
t = m;
int ans = 0;

queue<int> q;
memset(flow, 0, sizeof(flow));
for(;;) {
    memset(a, 0, sizeof(a));
    a[s] = INF;
    q.push(s);
    while(!q.empty()) {
      int u = q.front(); q.pop();
      for(int v = 1; v <=m; v++) if(!a[v] && cap[u][v]>flow[u][v]){
        p[v] = u; q.push(v);
        //a[v] = a[u] <? cap[u][v]-flow[u][v];
        if(a[u]<cap[u][v]-flow[u][v]) a[v]=a[u];
        else a[v]=cap[u][v]-flow[u][v];
      }
    }
    if(a[t] == 0) break;
    for(int u = t; u != s; u = p[u]) {
      flow[p[u]][u] += a[t];
      flow[u][p[u]] -= a[t];
    }
    ans += a[t];
}

printf("%d/n", ans);
}
return 0;
}