/// <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); }