Алгоритм Форда-Фалкерсона: различия между версиями
Перейти к навигации
Перейти к поиску
Ctrlalt (обсуждение | вклад) (→Ссылки) |
Ctrlalt (обсуждение | вклад) Нет описания правки |
||
| Строка 1: | Строка 1: | ||
class Graph { | |||
struct Edge { | |||
int a, b, capacity, flow; | |||
Edge(int a, int b, int capacity) : | |||
a(a), b(b), capacity(capacity), flow(0) {} | |||
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 f) { | |||
flow += (v == b ? f : -f); | |||
} | |||
}; | |||
vector<Edge> edges; | vector<Edge> edges; | ||
vector<vector<int>> g; | vector<vector<int>> g; | ||
vector<bool> | vector<bool> visited; | ||
vector<int> edgeTo; | vector<int> edgeTo; | ||
void dfs(int v) { | void dfs(int v) { | ||
visited[v] = 1; | |||
for (int e : g[v]) { | for (int e : g[v]) { | ||
int to = edges[e].other(v); | int to = edges[e].other(v); | ||
if (! | if (!visited[to] && edges[e].capacityTo(to) > 0) { | ||
edgeTo[to] = e; | edgeTo[to] = e; | ||
dfs(to); | dfs(to); | ||
| Строка 34: | Строка 34: | ||
} | } | ||
} | } | ||
bool hasPath(int from, int to) { | bool hasPath(int from, int to) { | ||
fill( | fill(visited.begin(), visited.end(), 0); | ||
dfs(from); | dfs(from); | ||
return | return visited[to]; | ||
} | } | ||
int bottleneckCapacity(int from, int to) { | int bottleneckCapacity(int from, int to) { | ||
int bCapacity = 1e9; | int bCapacity = 1e9; | ||
| Строка 45: | Строка 47: | ||
return bCapacity; | return bCapacity; | ||
} | } | ||
void addFlow(int from, int to, int flow) { | void addFlow(int from, int to, int flow) { | ||
for (int v = to; v != from; v = edges[edgeTo[v]].other(v)) | for (int v = to; v != from; v = edges[edgeTo[v]].other(v)) | ||
edges[edgeTo[v]].addFlowTo(v, flow); | edges[edgeTo[v]].addFlowTo(v, flow); | ||
} | } | ||
public: | public: | ||
Graph(int | Graph(int vertexCount) { | ||
g.resize( | g.resize(vertexCount); | ||
visited.resize(vertexCount); | |||
edgeTo.resize( | edgeTo.resize(vertexCount); | ||
} | } | ||
void addEdge(int from, int to, int capacity) { | void addEdge(int from, int to, int capacity) { | ||
edges.push_back(Edge(from, to, capacity)); | edges.push_back(Edge(from, to, capacity)); | ||
g[from].push_back(edges.size() - 1); | g[from].push_back(edges.size() - 1); | ||
g[ to ].push_back(edges.size() - 1); | g[to].push_back(edges.size() - 1); | ||
} | } | ||
long long maxFlow(int from, int to) { | long long maxFlow(int from, int to) { | ||
long long flow = 0; | long long flow = 0; | ||
| Строка 70: | Строка 76: | ||
} | } | ||
}; | }; | ||
== Ссылки == | == Ссылки == | ||
Версия от 02:52, 26 декабря 2021
class Graph {
struct Edge {
int a, b, capacity, flow;
Edge(int a, int b, int capacity) :
a(a), b(b), capacity(capacity), flow(0) {}
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 f) {
flow += (v == b ? f : -f);
}
};
vector<Edge> edges;
vector<vector<int>> g;
vector<bool> visited;
vector<int> edgeTo;
void dfs(int v) {
visited[v] = 1;
for (int e : g[v]) {
int to = edges[e].other(v);
if (!visited[to] && edges[e].capacityTo(to) > 0) {
edgeTo[to] = e;
dfs(to);
}
}
}
bool hasPath(int from, int to) {
fill(visited.begin(), visited.end(), 0);
dfs(from);
return visited[to];
}
int bottleneckCapacity(int from, int to) {
int bCapacity = 1e9;
for (int v = to; v != from; v = edges[edgeTo[v]].other(v))
bCapacity = min(bCapacity, edges[edgeTo[v]].capacityTo(v));
return bCapacity;
}
void addFlow(int from, int to, int flow) {
for (int v = to; v != from; v = edges[edgeTo[v]].other(v))
edges[edgeTo[v]].addFlowTo(v, flow);
}
public:
Graph(int vertexCount) {
g.resize(vertexCount);
visited.resize(vertexCount);
edgeTo.resize(vertexCount);
}
void addEdge(int from, int to, int capacity) {
edges.push_back(Edge(from, to, capacity));
g[from].push_back(edges.size() - 1);
g[to].push_back(edges.size() - 1);
}
long long maxFlow(int from, int to) {
long long flow = 0;
while (hasPath(from, to)) {
int deltaFlow = bottleneckCapacity(from, to);
addFlow(from, to, 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 (несмотря на название, используется алгоритм Эдмондса-Карпа)
Задачи: