public static float MaximumFlow(FlowNetwork fn, int source, int sink)
{
float maxFlow = 0f;
int sourceIndex = source;
int sinkIndex = sink;
foreach (FlowEdge e in fn.GetEdges())
{
e.CurrentFlow = 0f;
}
Graph g = CreateResidualNetwork(fn);
List <Edge> augmentingPath = BreadthFirstSearch(g, sourceIndex, sinkIndex);
while (augmentingPath != null)
{
float minFlow = float.MaxValue;
//find min flow
foreach (Edge e in augmentingPath)
{
if (e.Weight < minFlow)
{
minFlow = e.Weight;
}
}
maxFlow += minFlow;
foreach (Edge e in augmentingPath)
{
//Update residual graph
Edge currentForwardResid = g.FindEdge(e);
Edge currentBackwardResid = g.FindEdge(new Edge(e.GetToVertex(), e.GetFromVertex()));
float newForwardWeight = currentForwardResid.Weight - minFlow;
if (newForwardWeight <= 0f)
{
g.RemoveEdge(currentForwardResid);
}
else
{
currentForwardResid.Weight = newForwardWeight;
}
if (currentBackwardResid != null)
{
float newBackwardWeight = currentBackwardResid.Weight + minFlow;
currentBackwardResid.Weight = newBackwardWeight;
}
else
{
currentBackwardResid = new Edge(e.GetToVertex(), e.GetFromVertex(), minFlow);
g.AddEdge(currentBackwardResid);
}
//Update flow network
FlowEdge fe = fn.FindEdge(new FlowEdge(e.GetFromVertex(), e.GetToVertex(), 0, 0));
fe.CurrentFlow = fe.CurrentFlow + minFlow;
}
augmentingPath = BreadthFirstSearch(g, sourceIndex, sinkIndex);
//augmentingPath = null;
}
return(maxFlow);
}
