// C++ program for finding minimum cut using Ford-Fulkerson
#include <iostream>
#include <limits.h>
#include <string.h>
#include <queue>
using namespace std;
// Number of vertices in given graph
#define V 6
/* Returns true if there is a path from source 's' to sink 't' in
residual graph. Also fills parent[] to store the path */
int bfs(int rGraph[V][V], int s, int t, int parent[])
{
// Create a visited array and mark all vertices as not visited
bool visited[V];
memset(visited, 0, sizeof(visited));
// Create a queue, enqueue source vertex and mark source vertex
// as visited
queue <int> q;
q.push(s);
visited[s] = true;
parent[s] = -1;
// Standard BFS Loop
while (!q.empty())
{
int u = q.front();
q.pop();
for (int v=0; v<V; v++)
{
if (visited[v]==false && rGraph[u][v] > 0)
{
q.push(v);
parent[v] = u;
visited[v] = true;
}
}
}
// If we reached sink in BFS starting from source, then return
// true, else false
return (visited[t] == true);
}
// A DFS based function to find all reachable vertices from s. The function
// marks visited[i] as true if i is reachable from s. The initial values in
// visited[] must be false. We can also use BFS to find reachable vertices
void dfs(int rGraph[V][V], int s, bool visited[])
{
visited[s] = true;
for (int i = 0; i < V; i++)
if (rGraph[s][i] && !visited[i])
dfs(rGraph, i, visited);
}
// Prints the minimum s-t cut
void minCut(int graph[V][V], int s, int t)
{
int u, v;
// Create a residual graph and fill the residual graph with
// given capacities in the original graph as residual capacities
// in residual graph
int rGraph[V][V]; // rGraph[i][j] indicates residual capacity of edge i-j
for (u = 0; u < V; u++)
for (v = 0; v < V; v++)
rGraph[u][v] = graph[u][v];
int parent[V]; // This array is filled by BFS and to store path
// Augment the flow while there is a path from source to sink
while (bfs(rGraph, s, t, parent))
{
// Find minimum residual capacity of the edhes along the
// path filled by BFS. Or we can say find the maximum flow
// through the path found.
int path_flow = INT_MAX;
for (v=t; v!=s; v=parent[v])
{
u = parent[v];
path_flow = min(path_flow, rGraph[u][v]);
}
// update residual capacities of the edges and reverse edges
// along the path
for (v=t; v != s; v=parent[v])
{
u = parent[v];
rGraph[u][v] -= path_flow;
rGraph[v][u] += path_flow;
}
}
// Flow is maximum now, find vertices reachable from s
bool visited[V];
memset(visited, false, sizeof(visited));
dfs(rGraph, s, visited);
// Print all edges that are from a reachable vertex to
// non-reachable vertex in the original graph
for (int i = 0; i < V; i++)
for (int j = 0; j < V; j++)
if (visited[i] && !visited[j] && graph[i][j])
cout << i << " - " << j << endl;
return;
}
// Driver program to test above functions
int main()
{
// Let us create a graph shown in the above example
int graph[V][V] = { {0, 16, 13, 0, 0, 0},
{0, 0, 10, 12, 0, 0},
{0, 4, 0, 0, 14, 0},
{0, 0, 9, 0, 0, 20},
{0, 0, 0, 7, 0, 4},
{0, 0, 0, 0, 0, 0}
};
minCut(graph, 0, 5);
return 0;
}
Ly8gQysrIHByb2dyYW0gZm9yIGZpbmRpbmcgbWluaW11bSBjdXQgdXNpbmcgRm9yZC1GdWxrZXJzb24gCiNpbmNsdWRlIDxpb3N0cmVhbT4gCiNpbmNsdWRlIDxsaW1pdHMuaD4gCiNpbmNsdWRlIDxzdHJpbmcuaD4gCiNpbmNsdWRlIDxxdWV1ZT4gCnVzaW5nIG5hbWVzcGFjZSBzdGQ7IAoKLy8gTnVtYmVyIG9mIHZlcnRpY2VzIGluIGdpdmVuIGdyYXBoIAojZGVmaW5lIFYgNiAKCi8qIFJldHVybnMgdHJ1ZSBpZiB0aGVyZSBpcyBhIHBhdGggZnJvbSBzb3VyY2UgJ3MnIHRvIHNpbmsgJ3QnIGluIApyZXNpZHVhbCBncmFwaC4gQWxzbyBmaWxscyBwYXJlbnRbXSB0byBzdG9yZSB0aGUgcGF0aCAqLwppbnQgYmZzKGludCByR3JhcGhbVl1bVl0sIGludCBzLCBpbnQgdCwgaW50IHBhcmVudFtdKSAKeyAKCS8vIENyZWF0ZSBhIHZpc2l0ZWQgYXJyYXkgYW5kIG1hcmsgYWxsIHZlcnRpY2VzIGFzIG5vdCB2aXNpdGVkIAoJYm9vbCB2aXNpdGVkW1ZdOyAKCW1lbXNldCh2aXNpdGVkLCAwLCBzaXplb2YodmlzaXRlZCkpOyAKCgkvLyBDcmVhdGUgYSBxdWV1ZSwgZW5xdWV1ZSBzb3VyY2UgdmVydGV4IGFuZCBtYXJrIHNvdXJjZSB2ZXJ0ZXggCgkvLyBhcyB2aXNpdGVkIAoJcXVldWUgPGludD4gcTsgCglxLnB1c2gocyk7IAoJdmlzaXRlZFtzXSA9IHRydWU7IAoJcGFyZW50W3NdID0gLTE7IAoKCS8vIFN0YW5kYXJkIEJGUyBMb29wIAoJd2hpbGUgKCFxLmVtcHR5KCkpIAoJeyAKCQlpbnQgdSA9IHEuZnJvbnQoKTsgCgkJcS5wb3AoKTsgCgoJCWZvciAoaW50IHY9MDsgdjxWOyB2KyspIAoJCXsgCgkJCWlmICh2aXNpdGVkW3ZdPT1mYWxzZSAmJiByR3JhcGhbdV1bdl0gPiAwKSAKCQkJeyAKCQkJCXEucHVzaCh2KTsgCgkJCQlwYXJlbnRbdl0gPSB1OyAKCQkJCXZpc2l0ZWRbdl0gPSB0cnVlOyAKCQkJfSAKCQl9IAoJfSAKCgkvLyBJZiB3ZSByZWFjaGVkIHNpbmsgaW4gQkZTIHN0YXJ0aW5nIGZyb20gc291cmNlLCB0aGVuIHJldHVybiAKCS8vIHRydWUsIGVsc2UgZmFsc2UgCglyZXR1cm4gKHZpc2l0ZWRbdF0gPT0gdHJ1ZSk7IAp9IAoKLy8gQSBERlMgYmFzZWQgZnVuY3Rpb24gdG8gZmluZCBhbGwgcmVhY2hhYmxlIHZlcnRpY2VzIGZyb20gcy4gVGhlIGZ1bmN0aW9uIAovLyBtYXJrcyB2aXNpdGVkW2ldIGFzIHRydWUgaWYgaSBpcyByZWFjaGFibGUgZnJvbSBzLiBUaGUgaW5pdGlhbCB2YWx1ZXMgaW4gCi8vIHZpc2l0ZWRbXSBtdXN0IGJlIGZhbHNlLiBXZSBjYW4gYWxzbyB1c2UgQkZTIHRvIGZpbmQgcmVhY2hhYmxlIHZlcnRpY2VzIAp2b2lkIGRmcyhpbnQgckdyYXBoW1ZdW1ZdLCBpbnQgcywgYm9vbCB2aXNpdGVkW10pIAp7IAoJdmlzaXRlZFtzXSA9IHRydWU7IAoJZm9yIChpbnQgaSA9IDA7IGkgPCBWOyBpKyspIAoJaWYgKHJHcmFwaFtzXVtpXSAmJiAhdmlzaXRlZFtpXSkgCgkJZGZzKHJHcmFwaCwgaSwgdmlzaXRlZCk7IAp9IAoKLy8gUHJpbnRzIHRoZSBtaW5pbXVtIHMtdCBjdXQgCnZvaWQgbWluQ3V0KGludCBncmFwaFtWXVtWXSwgaW50IHMsIGludCB0KSAKeyAKCWludCB1LCB2OyAKCgkvLyBDcmVhdGUgYSByZXNpZHVhbCBncmFwaCBhbmQgZmlsbCB0aGUgcmVzaWR1YWwgZ3JhcGggd2l0aCAKCS8vIGdpdmVuIGNhcGFjaXRpZXMgaW4gdGhlIG9yaWdpbmFsIGdyYXBoIGFzIHJlc2lkdWFsIGNhcGFjaXRpZXMgCgkvLyBpbiByZXNpZHVhbCBncmFwaCAKCWludCByR3JhcGhbVl1bVl07IC8vIHJHcmFwaFtpXVtqXSBpbmRpY2F0ZXMgcmVzaWR1YWwgY2FwYWNpdHkgb2YgZWRnZSBpLWogCglmb3IgKHUgPSAwOyB1IDwgVjsgdSsrKSAKCQlmb3IgKHYgPSAwOyB2IDwgVjsgdisrKSAKCQkJckdyYXBoW3VdW3ZdID0gZ3JhcGhbdV1bdl07IAoKCWludCBwYXJlbnRbVl07IC8vIFRoaXMgYXJyYXkgaXMgZmlsbGVkIGJ5IEJGUyBhbmQgdG8gc3RvcmUgcGF0aCAKCgkvLyBBdWdtZW50IHRoZSBmbG93IHdoaWxlIHRoZXJlIGlzIGEgcGF0aCBmcm9tIHNvdXJjZSB0byBzaW5rIAoJd2hpbGUgKGJmcyhyR3JhcGgsIHMsIHQsIHBhcmVudCkpIAoJeyAKCQkvLyBGaW5kIG1pbmltdW0gcmVzaWR1YWwgY2FwYWNpdHkgb2YgdGhlIGVkaGVzIGFsb25nIHRoZSAKCQkvLyBwYXRoIGZpbGxlZCBieSBCRlMuIE9yIHdlIGNhbiBzYXkgZmluZCB0aGUgbWF4aW11bSBmbG93IAoJCS8vIHRocm91Z2ggdGhlIHBhdGggZm91bmQuIAoJCWludCBwYXRoX2Zsb3cgPSBJTlRfTUFYOyAKCQlmb3IgKHY9dDsgdiE9czsgdj1wYXJlbnRbdl0pIAoJCXsgCgkJCXUgPSBwYXJlbnRbdl07IAoJCQlwYXRoX2Zsb3cgPSBtaW4ocGF0aF9mbG93LCByR3JhcGhbdV1bdl0pOyAKCQl9IAoKCQkvLyB1cGRhdGUgcmVzaWR1YWwgY2FwYWNpdGllcyBvZiB0aGUgZWRnZXMgYW5kIHJldmVyc2UgZWRnZXMgCgkJLy8gYWxvbmcgdGhlIHBhdGggCgkJZm9yICh2PXQ7IHYgIT0gczsgdj1wYXJlbnRbdl0pIAoJCXsgCgkJCXUgPSBwYXJlbnRbdl07IAoJCQlyR3JhcGhbdV1bdl0gLT0gcGF0aF9mbG93OyAKCQkJckdyYXBoW3ZdW3VdICs9IHBhdGhfZmxvdzsgCgkJfSAKCX0gCgoJLy8gRmxvdyBpcyBtYXhpbXVtIG5vdywgZmluZCB2ZXJ0aWNlcyByZWFjaGFibGUgZnJvbSBzIAoJYm9vbCB2aXNpdGVkW1ZdOyAKCW1lbXNldCh2aXNpdGVkLCBmYWxzZSwgc2l6ZW9mKHZpc2l0ZWQpKTsgCglkZnMockdyYXBoLCBzLCB2aXNpdGVkKTsgCgoJLy8gUHJpbnQgYWxsIGVkZ2VzIHRoYXQgYXJlIGZyb20gYSByZWFjaGFibGUgdmVydGV4IHRvIAoJLy8gbm9uLXJlYWNoYWJsZSB2ZXJ0ZXggaW4gdGhlIG9yaWdpbmFsIGdyYXBoIAoJZm9yIChpbnQgaSA9IDA7IGkgPCBWOyBpKyspIAoJZm9yIChpbnQgaiA9IDA7IGogPCBWOyBqKyspIAoJCWlmICh2aXNpdGVkW2ldICYmICF2aXNpdGVkW2pdICYmIGdyYXBoW2ldW2pdKSAKCQkJY291dCA8PCBpIDw8ICIgLSAiIDw8IGogPDwgZW5kbDsgCgoJcmV0dXJuOyAKfSAKCi8vIERyaXZlciBwcm9ncmFtIHRvIHRlc3QgYWJvdmUgZnVuY3Rpb25zIAppbnQgbWFpbigpIAp7IAoJLy8gTGV0IHVzIGNyZWF0ZSBhIGdyYXBoIHNob3duIGluIHRoZSBhYm92ZSBleGFtcGxlIAoJaW50IGdyYXBoW1ZdW1ZdID0geyB7MCwgMTYsIDEzLCAwLCAwLCAwfSwgCgkJCQkJCXswLCAwLCAxMCwgMTIsIDAsIDB9LCAKCQkJCQkJezAsIDQsIDAsIDAsIDE0LCAwfSwgCgkJCQkJCXswLCAwLCA5LCAwLCAwLCAyMH0sIAoJCQkJCQl7MCwgMCwgMCwgNywgMCwgNH0sIAoJCQkJCQl7MCwgMCwgMCwgMCwgMCwgMH0gCgkJCQkJfTsgCgoJbWluQ3V0KGdyYXBoLCAwLCA1KTsgCgoJcmV0dXJuIDA7IAp9IAo=