Алгоритм Форда-Фалкерсона: различия между версиями

Материал из Олимпиадное программирование в УлГТУ
Перейти к навигации Перейти к поиску
(Новая страница: « #include <stdio.h> #include <algorithm> #include <vector> using namespace std; struct Edge { int a, b, cap, flow; Edge(int a, int b, int cap) : a…»)
 
Строка 58: Строка 58:
* [http://informatics.mccme.ru/course/view.php?id=6 informatics.mccme.ru &mdash; Курс &laquo;Алгоритмы на графах&raquo; &mdash; часть 9]
* [http://informatics.mccme.ru/course/view.php?id=6 informatics.mccme.ru &mdash; Курс &laquo;Алгоритмы на графах&raquo; &mdash; часть 9]
* [http://github.com/indy256/codelibrary/blob/master/java/src/MaxFlowFordFulkerson.java CodeLibrary &mdash; Maximum flow. Ford-Fulkerson alogithm in O(V^2 * FLOW)]
* [http://github.com/indy256/codelibrary/blob/master/java/src/MaxFlowFordFulkerson.java CodeLibrary &mdash; Maximum flow. Ford-Fulkerson alogithm in O(V^2 * FLOW)]
* [http://github.com/ADJA/algos/blob/master/Graphs/FordFulkerson.cpp Algos &mdash; Ford-Fulkerson maxflow)]
* [http://github.com/ADJA/algos/blob/master/Graphs/FordFulkerson.cpp Algos &mdash; Ford-Fulkerson maxflow]
* [http://visualgo.net/maxflow.html VisuAlgo &mdash; Network Flow]
* [http://visualgo.net/maxflow.html VisuAlgo &mdash; Network Flow]

Версия от 22:10, 23 мая 2015

#include <stdio.h>
#include <algorithm>
#include <vector>
using namespace std;

struct Edge {
    int a, b, cap, flow;
    Edge(int a, int b, int cap) : a(a), b(b), cap(cap), flow(0) {}
};
vector<Edge> edges;

int n, m, u[1010], INF = 1 << 30;
vector<int> g[1010];

int findPath(int v, int bottleneckFlow) {
    u[v] = 1;
    if (v == n - 1)
        return bottleneckFlow;
    for (int i = 0; i < g[v].size(); i++) {
       Edge &e = edges[g[v][i]], &re = edges[g[v][i] ^ 1];
       if (!u[e.b] && e.cap > e.flow) {
           int deltaFlow = findPath(e.b, min(bottleneckFlow, e.cap - e.flow));
           if (deltaFlow) {
               e.flow += deltaFlow;
               re.flow -= deltaFlow;
               return deltaFlow;
           }
       }
    }
    return 0;
}

int main() {
    scanf("%d%d", &n, &m);
    int a, b, cap;
    for (int i = 0; i < m; i++) {
        scanf("%d%d%d", &a, &b, &cap);
        edges.push_back(Edge(a - 1, b - 1, cap));
        g[a - 1].push_back(edges.size() - 1);
        edges.push_back(Edge(b - 1, a - 1, cap));
        g[b - 1].push_back(edges.size() - 1);
    }
    
    int flow = 0;
    while (1) {
        fill(u, u + n, 0);
        int deltaFlow = findPath(0, INF);
        if (!deltaFlow)
            break;
        flow += deltaFlow;
    }
    printf("%d", flow);
}

Ссылки