最大流算法之Ford-Fulkerson

最大流算法

该算法的核心就是找建立残余网络并且找增广路，挑战程序设计上讲的很清楚，其实找增广路就是看有没有一条路径可以从起点走到终点的路径，如果没有那么肯定就不存在增广路了，那么此时找到的流就是最大流。

HDU 1532 Drainage Ditches

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
const int maxn = 205;
const int INF = 0x3f3f3f3f;
int head[maxn], tot, vis[maxn];
struct Edge{
    int to, cap, next;
}e[maxn<<2];

void addedge(int from, int to, int cap){
    e[tot].to = to;
    e[tot].cap = cap;
    e[tot].next = head[from];
    head[from] = tot++;
}

int dfs(int u, int v,int f){   //dfs找增广路
    if(u == v) return f;
    vis[u] = 1;   //对于走过的边进行标记
    for(int i = head[u]; i != -1; i = e[i].next){  //邻接表
        if(!vis[e[i].to] && e[i].cap > 0){
            int d = dfs(e[i].to, v, min(f, e[i].cap));
            if(d > 0){
                e[i].cap -= d;   //残余网络中正向边减流量，反向边加流量
                e[i^1].cap += d;
                return d;
            }
        }
    }
    return 0;
}

int max_frow(int s, int t){
    int flow = 0;
    for(;;){
        memset(vis, 0, sizeof(vis));
        int f = dfs(s, t, INF);
        flow += f;
        if(f == 0) break;   //直到找不到增广路那么输出当前最大流
    }
    return flow;
}

void init(){
    memset(head, -1, sizeof(head));
    tot = 0;
}

int main(){
    int n, m;
    while(scanf("%d%d", &n, &m) != EOF){
        int from, to, cap;
        init();
        for(int i = 0; i < n; i++){
            scanf("%d%d%d", &from, &to, &cap);
            addedge(from, to, cap);
            addedge(to, from, 0);
        }
        printf("%d\n", max_frow(1, m));
    }
    return 0;
}