| | 1 | | using System; |
| | 2 | | using System.Collections.Generic; |
| | 3 | |
|
| | 4 | | namespace Itinero.Network.Search.Islands; |
| | 5 | |
|
| | 6 | | internal class IslandLabelGraph |
| | 7 | | { |
| | 8 | | private const uint NoEdge = uint.MaxValue; |
| | 9 | | private const uint NoVertex = uint.MaxValue - 1; |
| | 10 | |
|
| 13 | 11 | | private readonly List<uint> _vertices = new(); // Holds all vertices pointing to it's first edge. |
| 13 | 12 | | private readonly List<(uint vertex1, uint vertex2, uint pointer1, uint pointer2)> _edges = new(); |
| | 13 | |
|
| | 14 | | public uint AddVertex() |
| 42 | 15 | | { |
| 42 | 16 | | var vertex = (uint)_vertices.Count; |
| 42 | 17 | | _vertices.Add(NoEdge); |
| 42 | 18 | | return vertex; |
| 42 | 19 | | } |
| | 20 | |
|
| | 21 | | public bool HasVertex(uint vertex) |
| 399 | 22 | | { |
| 399 | 23 | | if (vertex >= _vertices.Count) return false; |
| 399 | 24 | | if (_vertices[(int)vertex] == NoVertex) return false; |
| | 25 | |
|
| 399 | 26 | | return true; |
| 399 | 27 | | } |
| | 28 | |
|
| | 29 | | public bool RemoveVertex(uint vertex) |
| 20 | 30 | | { |
| 20 | 31 | | if (!this.HasVertex(vertex)) return false; |
| | 32 | |
|
| 20 | 33 | | this.RemoveEdges(vertex); |
| | 34 | |
|
| 20 | 35 | | _vertices[(int)vertex] = NoVertex; |
| | 36 | |
|
| 40 | 37 | | while (_vertices.Count > 0 && |
| 40 | 38 | | _vertices[^1] == NoVertex) |
| 20 | 39 | | { |
| 20 | 40 | | _vertices.RemoveAt(_vertices.Count - 1); |
| 20 | 41 | | } |
| | 42 | |
|
| 20 | 43 | | return true; |
| 20 | 44 | | } |
| | 45 | |
|
| | 46 | | public void AddEdge(uint vertex1, uint vertex2) |
| 32 | 47 | | { |
| 32 | 48 | | if (vertex1 == vertex2) { throw new ArgumentException("Given vertices must be different."); } |
| 32 | 49 | | if (!this.HasVertex(vertex1)) throw new ArgumentException($"Vertex {vertex1} does not exist."); |
| 32 | 50 | | if (!this.HasVertex(vertex2)) throw new ArgumentException($"Vertex {vertex2} does not exist."); |
| | 51 | |
|
| 32 | 52 | | this.AddEdgeInternal(vertex1, vertex2); |
| 32 | 53 | | } |
| | 54 | |
|
| | 55 | | public bool HasEdge(uint vertex1, uint vertex2) |
| 33 | 56 | | { |
| 33 | 57 | | if (vertex1 == vertex2) { throw new ArgumentException("Given vertices must be different."); } |
| 33 | 58 | | if (!this.HasVertex(vertex1)) throw new ArgumentException($"Vertex {vertex1} does not exist."); |
| 33 | 59 | | if (!this.HasVertex(vertex2)) throw new ArgumentException($"Vertex {vertex2} does not exist."); |
| | 60 | |
|
| 33 | 61 | | var edgeId = _vertices[(int)vertex1]; |
| 48 | 62 | | while (edgeId != NoEdge) |
| 16 | 63 | | { |
| | 64 | | // get the edge. |
| 16 | 65 | | var edge = _edges[(int)edgeId]; |
| 16 | 66 | | if (edge.vertex1 == vertex1) |
| 1 | 67 | | { |
| 1 | 68 | | if (edge.vertex2 == vertex2) |
| 1 | 69 | | { |
| 1 | 70 | | return true; |
| | 71 | | } |
| 0 | 72 | | edgeId = edge.pointer1; |
| 0 | 73 | | } |
| 15 | 74 | | else if (edge.vertex2 == vertex1) |
| 15 | 75 | | { |
| 15 | 76 | | edgeId = edge.pointer2; |
| 15 | 77 | | } |
| | 78 | | else |
| 0 | 79 | | { |
| 0 | 80 | | throw new Exception("Edge found on vertex set that does not contain vertex"); |
| | 81 | | } |
| 15 | 82 | | } |
| | 83 | |
|
| 32 | 84 | | return false; |
| 33 | 85 | | } |
| | 86 | |
|
| | 87 | | public int RemoveEdges(uint vertex) |
| 20 | 88 | | { |
| 20 | 89 | | var removed = 0; |
| 20 | 90 | | var edges = this.GetEdgeEnumerator(); |
| 49 | 91 | | while (edges.MoveTo(vertex) && edges.MoveNext()) |
| 29 | 92 | | { |
| 29 | 93 | | if (edges.Forward) |
| 14 | 94 | | { |
| 14 | 95 | | this.RemoveEdge(vertex, edges.Head); |
| 14 | 96 | | } |
| | 97 | | else |
| 15 | 98 | | { |
| 15 | 99 | | this.RemoveEdge(edges.Head, vertex); |
| 15 | 100 | | } |
| 29 | 101 | | removed++; |
| 29 | 102 | | } |
| | 103 | |
|
| 20 | 104 | | return removed; |
| 20 | 105 | | } |
| | 106 | |
|
| | 107 | | public bool RemoveEdge(uint vertex1, uint vertex2) |
| 29 | 108 | | { |
| 29 | 109 | | if (vertex1 == vertex2) throw new ArgumentException("Given vertices must be different."); |
| 29 | 110 | | if (!this.HasVertex(vertex1)) throw new ArgumentException($"Vertex {vertex1} does not exist."); |
| 29 | 111 | | if (!this.HasVertex(vertex2)) throw new ArgumentException($"Vertex {vertex2} does not exist."); |
| | 112 | |
|
| 29 | 113 | | var edgeId = this.RemoveEdgeFromVertex1(vertex1, vertex2); |
| 29 | 114 | | if (edgeId == NoEdge) return false; |
| | 115 | |
|
| 29 | 116 | | this.RemoveEdgeFromVertex(vertex2, edgeId); |
| | 117 | |
|
| 29 | 118 | | _edges[(int)edgeId] = (NoVertex, NoVertex, NoEdge, NoEdge); |
| | 119 | |
|
| 58 | 120 | | while (_edges.Count > 0 && _edges[^1].vertex1 == NoVertex) |
| 29 | 121 | | { |
| 29 | 122 | | _edges.RemoveAt(_edges.Count - 1); |
| 29 | 123 | | } |
| 29 | 124 | | return true; |
| 29 | 125 | | } |
| | 126 | |
|
| | 127 | | private uint RemoveEdgeFromVertex1(uint vertex1, uint vertex2) |
| 29 | 128 | | { |
| | 129 | | // find the edge and keep the previous edge id. |
| 29 | 130 | | var edgeId = _vertices[(int)vertex1]; |
| 29 | 131 | | var previousEdgeId = NoEdge; |
| 29 | 132 | | var foundEdgeId = NoEdge; |
| 29 | 133 | | while (edgeId != NoEdge) |
| 29 | 134 | | { |
| | 135 | | // get the edge. |
| 29 | 136 | | var edge = _edges[(int)edgeId]; |
| 29 | 137 | | if (edge.vertex1 == vertex1) |
| 29 | 138 | | { |
| 29 | 139 | | if (edge.vertex2 == vertex2) |
| 29 | 140 | | { |
| 29 | 141 | | foundEdgeId = edgeId; |
| 29 | 142 | | edgeId = edge.pointer1; |
| 29 | 143 | | break; |
| | 144 | | } |
| 0 | 145 | | previousEdgeId = edgeId; |
| 0 | 146 | | edgeId = edge.pointer1; |
| 0 | 147 | | } |
| 0 | 148 | | else if (edge.vertex2 == vertex1) |
| 0 | 149 | | { |
| 0 | 150 | | previousEdgeId = edgeId; |
| 0 | 151 | | edgeId = edge.pointer2; |
| 0 | 152 | | } |
| | 153 | | else |
| 0 | 154 | | { |
| 0 | 155 | | throw new Exception("Edge found on vertex set that does not contain vertex"); |
| | 156 | | } |
| 0 | 157 | | } |
| | 158 | |
|
| 29 | 159 | | if (foundEdgeId == NoEdge) return NoEdge; |
| | 160 | |
|
| | 161 | | // set pointer on last edge or vertex. |
| 29 | 162 | | if (previousEdgeId == NoEdge) |
| 29 | 163 | | { |
| 29 | 164 | | _vertices[(int)vertex1] = edgeId; |
| 29 | 165 | | } |
| | 166 | | else |
| 0 | 167 | | { |
| 0 | 168 | | var edge = _edges[(int)previousEdgeId]; |
| 0 | 169 | | if (edge.vertex1 == vertex1) |
| 0 | 170 | | { |
| 0 | 171 | | edge.pointer1 = edgeId; |
| 0 | 172 | | } |
| 0 | 173 | | else if (edge.vertex2 == vertex1) |
| 0 | 174 | | { |
| 0 | 175 | | edge.pointer2 = edgeId; |
| 0 | 176 | | } |
| 0 | 177 | | _edges[(int)previousEdgeId] = edge; |
| 0 | 178 | | } |
| | 179 | |
|
| 29 | 180 | | return foundEdgeId; |
| 29 | 181 | | } |
| | 182 | |
|
| | 183 | | private void RemoveEdgeFromVertex(uint vertex, uint edgeIdToRemove) |
| 29 | 184 | | { |
| 29 | 185 | | if (edgeIdToRemove == NoEdge) throw new Exception("This edge has to exist, it existing in other vertex"); |
| | 186 | |
|
| 29 | 187 | | var edgeId = _vertices[(int)vertex]; |
| 29 | 188 | | var previousEdgeId = NoEdge; |
| 30 | 189 | | while (edgeId != NoEdge) |
| 30 | 190 | | { |
| | 191 | | // get the edge. |
| 30 | 192 | | var edge = _edges[(int)edgeId]; |
| 30 | 193 | | var currentEdgeId = edgeId; |
| | 194 | |
|
| 30 | 195 | | if (edge.vertex1 == vertex) |
| 1 | 196 | | { |
| 1 | 197 | | edgeId = edge.pointer1; |
| 1 | 198 | | } |
| 29 | 199 | | else if (edge.vertex2 == vertex) |
| 29 | 200 | | { |
| 29 | 201 | | edgeId = edge.pointer2; |
| 29 | 202 | | } |
| | 203 | | else |
| 0 | 204 | | { |
| 0 | 205 | | throw new Exception("Edge found on vertex set that does not contain vertex"); |
| | 206 | | } |
| | 207 | |
|
| 59 | 208 | | if (currentEdgeId == edgeIdToRemove) break; |
| 1 | 209 | | previousEdgeId = currentEdgeId; |
| 1 | 210 | | } |
| | 211 | |
|
| 29 | 212 | | if (previousEdgeId == NoEdge) |
| 28 | 213 | | { |
| 28 | 214 | | _vertices[(int)vertex] = edgeId; |
| 28 | 215 | | } |
| | 216 | | else |
| 1 | 217 | | { |
| 1 | 218 | | var edge = _edges[(int)previousEdgeId]; |
| 1 | 219 | | if (edge.vertex1 == vertex) |
| 1 | 220 | | { |
| 1 | 221 | | edge.pointer1 = edgeId; |
| 1 | 222 | | } |
| 0 | 223 | | else if (edge.vertex2 == vertex) |
| 0 | 224 | | { |
| 0 | 225 | | edge.pointer2 = edgeId; |
| 0 | 226 | | } |
| 1 | 227 | | _edges[(int)previousEdgeId] = edge; |
| 1 | 228 | | } |
| 29 | 229 | | } |
| | 230 | |
|
| | 231 | | public bool HasEdge(uint vertex1) |
| 32 | 232 | | { |
| 32 | 233 | | var enumerator = this.GetEdgeEnumerator(); |
| 32 | 234 | | if (!enumerator.MoveTo(vertex1)) return false; |
| | 235 | |
|
| 32 | 236 | | return enumerator.MoveNext(); |
| 32 | 237 | | } |
| | 238 | |
|
| | 239 | | public EdgeEnumerator GetEdgeEnumerator() |
| 85 | 240 | | { |
| 85 | 241 | | return new EdgeEnumerator(this); |
| 85 | 242 | | } |
| | 243 | |
|
| | 244 | | /// <summary> |
| | 245 | | /// Returns the number of vertices in this graph. |
| | 246 | | /// </summary> |
| 0 | 247 | | public int VertexCount => _vertices.Count; |
| | 248 | |
|
| | 249 | | /// <summary> |
| | 250 | | /// Returns the number of edges in this graph. |
| | 251 | | /// </summary> |
| 0 | 252 | | public long EdgeCount => _edges.Count; |
| | 253 | |
|
| | 254 | | private uint AddEdgeInternal(uint vertex1, uint vertex2) |
| 32 | 255 | | { |
| | 256 | | // this adds an edge without checking for duplicates! |
| 32 | 257 | | var vertex1EdgeId = _vertices[(int)vertex1]; |
| 32 | 258 | | var vertex2EdgeId = _vertices[(int)vertex2]; |
| | 259 | |
|
| 32 | 260 | | var newEdgeId = (uint)_edges.Count; |
| 32 | 261 | | _vertices[(int)vertex1] = newEdgeId; |
| 32 | 262 | | _vertices[(int)vertex2] = newEdgeId; |
| | 263 | |
|
| 32 | 264 | | _edges.Add((vertex1, vertex2, vertex1EdgeId, vertex2EdgeId)); |
| | 265 | |
|
| 32 | 266 | | return newEdgeId; |
| 32 | 267 | | } |
| | 268 | |
|
| | 269 | | /// <summary> |
| | 270 | | /// An edge enumerator. |
| | 271 | | /// </summary> |
| | 272 | | public class EdgeEnumerator |
| | 273 | | { |
| | 274 | | private readonly IslandLabelGraph _islandLabelGraph; |
| 85 | 275 | | private uint _nextEdgeId = NoEdge; |
| | 276 | |
|
| | 277 | | /// <summary> |
| | 278 | | /// Creates a new edge enumerator. |
| | 279 | | /// </summary> |
| 85 | 280 | | internal EdgeEnumerator(IslandLabelGraph islandLabelGraph) |
| 85 | 281 | | { |
| 85 | 282 | | _islandLabelGraph = islandLabelGraph; |
| 85 | 283 | | } |
| | 284 | |
|
| | 285 | | /// <summary> |
| | 286 | | /// Move to the next edge. |
| | 287 | | /// </summary> |
| | 288 | | /// <returns></returns> |
| | 289 | | public bool MoveNext() |
| 279 | 290 | | { |
| 279 | 291 | | if (this.Tail == NoVertex) return false; |
| 399 | 292 | | if (_nextEdgeId == NoEdge) return false; |
| | 293 | |
|
| 159 | 294 | | var edge = _islandLabelGraph._edges[(int)_nextEdgeId]; |
| 159 | 295 | | if (edge.vertex1 == this.Tail) |
| 86 | 296 | | { |
| 86 | 297 | | this.Head = edge.vertex2; |
| 86 | 298 | | _nextEdgeId = edge.pointer1; |
| 86 | 299 | | this.Forward = true; |
| 86 | 300 | | return true; |
| | 301 | | } |
| 73 | 302 | | if (edge.vertex2 == this.Tail) |
| 73 | 303 | | { |
| 73 | 304 | | this.Head = edge.vertex1; |
| 73 | 305 | | _nextEdgeId = edge.pointer2; |
| 73 | 306 | | this.Forward = false; |
| 73 | 307 | | return true; |
| | 308 | | } |
| | 309 | |
|
| 0 | 310 | | throw new Exception("Next edge does not have vertex"); |
| 279 | 311 | | } |
| | 312 | |
|
| 787 | 313 | | public uint Tail { get; private set; } = NoVertex; |
| | 314 | |
|
| 389 | 315 | | public uint Head { get; private set; } = NoVertex; |
| | 316 | |
|
| 331 | 317 | | public bool Forward { get; private set; } = true; |
| | 318 | |
|
| | 319 | | public bool MoveTo(uint vertex) |
| 191 | 320 | | { |
| 191 | 321 | | if (!_islandLabelGraph.HasVertex(vertex)) return false; |
| | 322 | |
|
| 191 | 323 | | _nextEdgeId = _islandLabelGraph._vertices[(int)vertex]; |
| 191 | 324 | | if (_nextEdgeId == NoVertex) return false; |
| | 325 | |
|
| 191 | 326 | | this.Tail = vertex; |
| 191 | 327 | | return true; |
| 191 | 328 | | } |
| | 329 | |
|
| | 330 | | public void Reset() |
| 0 | 331 | | { |
| 0 | 332 | | if (this.Tail == NoVertex) return; |
| | 333 | |
|
| 0 | 334 | | this.MoveTo(this.Tail); |
| 0 | 335 | | } |
| | 336 | | } |
| | 337 | | } |