using System; using System.Collections.Generic; using System.Linq; using Collections.Pooled; using MoreLinq; namespace MoonTools.Core.Graph { public class DirectedWeightedMultiGraph : MultiGraph where TNode : IEquatable { protected Dictionary weights = new Dictionary(); public Guid AddEdge(TNode v, TNode u, int weight, TEdgeData data) { var id = BaseAddEdge(v, u, data); weights.Add(id, weight); return id; } public void AddEdges(params (TNode, TNode, int, TEdgeData)[] edges) { foreach (var edge in edges) { AddEdge(edge.Item1, edge.Item2, edge.Item3, edge.Item4); } } public override void Clear() { base.Clear(); weights.Clear(); } public IEnumerable Weights(TNode v, TNode u) { CheckNodes(v, u); return edges[(v, u)].Select(id => weights[id]); } private IEnumerable ReconstructPath(PooledDictionary cameFrom, TNode currentNode) { while (cameFrom.ContainsKey(currentNode)) { var edgeID = cameFrom[currentNode]; var edge = IDToEdge[edgeID]; currentNode = edge.Item1; yield return edgeID; } } public IEnumerable AStarShortestPath(TNode start, TNode end, Func heuristic) { CheckNodes(start, end); var openSet = new PooledSet(ClearMode.Always); var closedSet = new PooledSet(ClearMode.Always); var gScore = new PooledDictionary(ClearMode.Always); var fScore = new PooledDictionary(ClearMode.Always); var cameFrom = new PooledDictionary(ClearMode.Always); openSet.Add(start); gScore[start] = 0; fScore[start] = heuristic(start, end); while (openSet.Any()) { var currentNode = openSet.MinBy(node => fScore[node]).First(); if (currentNode.Equals(end)) { openSet.Dispose(); closedSet.Dispose(); gScore.Dispose(); fScore.Dispose(); foreach (var edgeID in ReconstructPath(cameFrom, currentNode).Reverse()) { yield return edgeID; } cameFrom.Dispose(); yield break; } openSet.Remove(currentNode); closedSet.Add(currentNode); foreach (var neighbor in Neighbors(currentNode)) { if (!closedSet.Contains(neighbor)) { var lowestEdgeID = EdgeIDs(currentNode, neighbor).MinBy(id => weights[id]).First(); var weight = weights[lowestEdgeID]; var tentativeGScore = gScore.ContainsKey(currentNode) ? gScore[currentNode] + weight : int.MaxValue; if (!openSet.Contains(neighbor) || tentativeGScore < gScore[neighbor]) { cameFrom[neighbor] = lowestEdgeID; gScore[neighbor] = tentativeGScore; fScore[neighbor] = tentativeGScore + heuristic(neighbor, end); openSet.Add(neighbor); } } } } openSet.Dispose(); closedSet.Dispose(); gScore.Dispose(); fScore.Dispose(); cameFrom.Dispose(); yield break; } } }