/// <summary>
/// Znajduje najliczniejszy multizbiór tras
/// </summary>
/// <returns>Znaleziony multizbiór</returns>
/// <param name="g">Graf opisujący muzeum</param>
/// <param name="cLevel">Tablica o długości równej liczbie wierzchołków w grafie -- poziomy ciekawości wystaw</param>
/// <param name="entrances">Wejścia</param>
/// <param name="exits">Wyjścia</param>
public static MuseumRoutes FindRoutes(Graph g, int[] cLevel, int[] entrances, int[] exits)
{
Graph seekFlowGraph = BuildSeekFlowGraph(g, cLevel, entrances, exits);
//var ge = new GraphExport();
//ge.Export(g);
//ge.Export(seekFlowGraph);
var maxFlow = MaxFlowGraphExtender.FordFulkersonDinicMaxFlow(seekFlowGraph, 0, 1, MaxFlowGraphExtender.OriginalDinicBlockingFlow);
return(new MuseumRoutes((int)maxFlow.value, null));
}
/// <summary>
/// Metoda znajdująca cyrkulację w grafie, z określonymi żądaniami wierzchołków.
/// Żądania opisane są w tablicy demands. Szukamy funkcji, która dla każdego wierzchołka będzie spełniała warunek:
/// suma wartości na krawędziach wchodzących - suma wartości na krawędziach wychodzących = demands[v]
/// </summary>
/// <param name="G">Graf wejściowy, wagi krawędzi oznaczają przepustowości</param>
/// <param name="demands">Żądania wierzchołków</param>
/// <returns>Graf reprezentujący wynikową cyrkulację.
/// Reprezentacja cyrkulacji jest analogiczna, jak reprezentacja przepływu w innych funkcjach w bibliotece.
/// Należy zwrócić kopię grafu G, gdzie wagi krawędzi odpowiadają przepływom na tych krawędziach.
/// Zwróć uwagę na rozróżnienie sytuacji, kiedy mamy zerowy przeływ na krawędzi (czyli istnieje
/// krawędź z wagą 0) od sytuacji braku krawędzi.
/// Jeśli żądana cyrkulacja nie istnieje, zwróć null.
/// </returns>
/// <remarks>
/// Nie można modyfikować danych wejściowych!
/// Złożoność metody powinna być asymptotycznie równa złożoności metody znajdującej największy przeływ (z biblioteki).
/// </remarks>
public Graph FindCirculation(Graph G, double[] demands)
{
Graph Network = new AdjacencyMatrixGraph(true, G.VerticesCount + 2);
int n = G.VerticesCount;
//g.VerticesCount - zrodlo, g.VerticesCount + 1 - ujscie
//zrodlo z ujemnymi, dodatnie z ujsciem
double flows = 0;
for (int i = 0; i < demands.Length; i++)
{
flows += demands[i];
}
if (flows != 0)
{
return(null);
}
for (int i = 0; i < demands.Length; i++)
{
if (demands[i] < 0)
{
Network.AddEdge(new Edge(n, i, Math.Abs(demands[i])));
}
else
{
Network.AddEdge(new Edge(i, n + 1, demands[i]));
}
foreach (Edge e in G.OutEdges(i))
{
Network.Add(e);
}
}
(double max_flow, Graph Flow) = MaxFlowGraphExtender.FordFulkersonDinicMaxFlow(Network, n, n + 1, MaxFlowGraphExtender.DFSBlockingFlow);
Graph Circulation = new AdjacencyMatrixGraph(true, n);
double[] degrees = new double[n];
for (int i = 0; i < n; i++)
{
foreach (Edge e in Flow.OutEdges(i))
{
if (e.To != n + 1)
{
Circulation.Add(e);
degrees[i] -= e.Weight;
degrees[e.To] += e.Weight;
}
}
}
for (int i = 0; i < n; i++)
{
if (degrees[i] != demands[i])
{
return(null);
}
}
return(Circulation);
}
