diff options
| author | pjungeblut <paul.jungeblut@gmail.com> | 2014-10-30 19:26:30 +0100 |
|---|---|---|
| committer | pjungeblut <paul.jungeblut@gmail.com> | 2014-10-30 19:26:30 +0100 |
| commit | 07a1c4f87dcccbfe2ca4fbbc100a07c9be801502 (patch) | |
| tree | b69d372db88dadb7cb755ca4cdaccc7000d2c9c3 /graph/edmondsKarp.cpp | |
| parent | 0550627fade1c01d6b90907de198ca13d34db9d9 (diff) | |
SCCS added, seperated from 2-SAT
Diffstat (limited to 'graph/edmondsKarp.cpp')
| -rw-r--r-- | graph/edmondsKarp.cpp | 34 |
1 files changed, 34 insertions, 0 deletions
diff --git a/graph/edmondsKarp.cpp b/graph/edmondsKarp.cpp new file mode 100644 index 0000000..26c5b0d --- /dev/null +++ b/graph/edmondsKarp.cpp @@ -0,0 +1,34 @@ +int s, t, f; //source, target, single flow +int res[MAX_V][MAX_V]; //adj-matrix +vector< vector<int> > adjList; +int p[MAX_V]; //bfs spanning tree + +void augment(int v, int minEdge) { + if (v == s) { f = minEdge; return; } + else if (p[v] != -1) { + augment(p[v], min(minEdge, res[p[v]][v])); + res[p[v]][v] -= f; res[v][p[v]] += f; +}} + +int maxFlow() { //first inititalize res, adjList, s and t + int mf = 0; + while (true) { + f = 0; + bitset<MAX_V> vis; vis[s] = true; + queue<int> q; q.push(s); + memset(p, -1, sizeof(p)); + while (!q.empty()) { //BFS + int u = q.front(); q.pop(); + if (u == t) break; + for (int j = 0; j < (int)adjList[u].size(); j++) { + int v = adjList[u][j]; + if (res[u][v] > 0 && !vis[v]) { + vis[v] = true; q.push(v); p[v] = u; + }}} + + augment(t, INF); //add found path to max flow + if (f == 0) break; + mf += f; + } + return mf; +}
\ No newline at end of file |
