Алгоритм Форда-Фалкерсона
Перейти к навигации
Перейти к поиску
class Graph { struct Edge { int a, b, capacity, flow = 0; Edge(int a, int b, int capacity) : a(a), b(b), capacity(capacity) {} int other(int v) const { return v == a ? b : a; } int capacityTo(int v) const { return v == b ? capacity - flow : flow; } void addFlowTo(int v, int deltaFlow) { flow += (v == b ? deltaFlow : -deltaFlow); } }; vector<Edge> edges; vector<vector<int>> graph; vector<bool> visited; vector<int> edgeTo; void dfs(int v) { visited[v] = 1; for (int e : graph[v]) { int to = edges[e].other(v); if (!visited[to] && edges[e].capacityTo(to) > 0) { edgeTo[to] = e; dfs(to); } } } bool hasPath(int start, int finish) { visited.assign(visited.size(), 0); dfs(start); return visited[finish]; } int bottleneckCapacity(int start, int finish) { int bCapacity = 1e9; for (int v = finish; v != start; v = edges[edgeTo[v]].other(v)) bCapacity = min(bCapacity, edges[edgeTo[v]].capacityTo(v)); return bCapacity; } void addFlow(int start, int finish, int deltaFlow) { for (int v = finish; v != start; v = edges[edgeTo[v]].other(v)) edges[edgeTo[v]].addFlowTo(v, deltaFlow); } public: Graph(int vertexCount) : graph(vertexCount), visited(vertexCount), edgeTo(vertexCount) {} void addEdge(int from, int to, int capacity) { edges.push_back(Edge(from, to, capacity)); graph[from].push_back(edges.size() - 1); graph[to].push_back(edges.size() - 1); } long long maxFlow(int start, int finish) { long long flow = 0; while (hasPath(start, finish)) { int deltaFlow = bottleneckCapacity(start, finish); addFlow(start, finish, deltaFlow); flow += deltaFlow; } return flow; } };
Ссылки
Теория:
- algs4.cs.princeton.edu — 6.4 Maximum Flow
- neerc.ifmo.ru/wiki — Алгоритм Форда-Фалкерсона
- Brilliant.org — Ford-Fulkerson Algorithm
Демонстрация:
Код:
- CodeLibrary — Maximum flow. Ford-Fulkerson alogithm in O(V^2 * FLOW)
- Algos — Ford-Fulkerson maxflow
- algs4.cs.princeton.edu/code — capacitated edge with flow, capacitated network, maxflow–mincut (несмотря на название, используется алгоритм Эдмондса-Карпа)
Задачи: