1
|
|
|
<?php |
2
|
|
|
|
3
|
|
|
namespace Mouf\Database\SchemaAnalyzer; |
4
|
|
|
|
5
|
|
|
use Fhaculty\Graph\Edge; |
6
|
|
|
use Fhaculty\Graph\Exception\UnexpectedValueException; |
7
|
|
|
use Fhaculty\Graph\Vertex; |
8
|
|
|
use SplPriorityQueue; |
9
|
|
|
|
10
|
|
|
/** |
11
|
|
|
* Dijkstra's shortest path algorithm modified to measure all possible shortest paths. |
12
|
|
|
*/ |
13
|
|
|
class MultiDijkstra |
14
|
|
|
{ |
15
|
|
|
/** |
16
|
|
|
* Get all edges on shortest path for this vertex. |
17
|
|
|
* |
18
|
|
|
* @throws UnexpectedValueException when encountering an Edge with negative weight |
19
|
|
|
* @throws MultiDijkstraNoPathException |
20
|
|
|
* |
21
|
|
|
* @return array<string, Vertex[]> where key is the destination vertex name and value is an array of possible origin vertex |
|
|
|
|
22
|
|
|
*/ |
23
|
|
|
public static function findShortestPaths(Vertex $startVertex, Vertex $endVertex) |
24
|
|
|
{ |
25
|
|
|
$totalCostOfCheapestPathTo = []; |
26
|
|
|
// start node distance |
27
|
|
|
$totalCostOfCheapestPathTo[$startVertex->getId()] = 0; |
28
|
|
|
|
29
|
|
|
$endVertexId = $endVertex->getId(); |
30
|
|
|
|
31
|
|
|
// just to get the cheapest vertex in the correct order |
32
|
|
|
$cheapestVertex = new SplPriorityQueue(); |
33
|
|
|
$cheapestVertex->setExtractFlags(SplPriorityQueue::EXTR_BOTH); |
34
|
|
|
$cheapestVertex->insert($startVertex, 0); |
35
|
|
|
|
36
|
|
|
// predecessors |
37
|
|
|
$predecesEdgeOfCheapestPathTo = []; |
38
|
|
|
|
39
|
|
|
// mark vertices when their cheapest path has been found |
40
|
|
|
$usedVertices = [$startVertex->getId() => true]; |
41
|
|
|
|
42
|
|
|
$isFirst = true; |
43
|
|
|
|
44
|
|
|
// Repeat until all vertices have been marked |
45
|
|
|
$totalCountOfVertices = count($startVertex->getGraph()->getVertices()); |
46
|
|
|
for ($i = 0; $i < $totalCountOfVertices; ++$i) { |
47
|
|
|
$currentVertex = null; |
48
|
|
|
$currentVertexId = null; |
49
|
|
|
$isEmpty = false; |
50
|
|
|
if ($isFirst) { |
51
|
|
|
$currentVertex = $startVertex; |
52
|
|
|
$currentCost = 0; |
53
|
|
|
$currentVertexId = $currentVertex->getId(); |
54
|
|
|
} else { |
55
|
|
|
do { |
56
|
|
|
// if the priority queue is empty there are isolated vertices, but the algorithm visited all other vertices |
57
|
|
|
if ($cheapestVertex->isEmpty()) { |
58
|
|
|
$isEmpty = true; |
59
|
|
|
break; |
60
|
|
|
} |
61
|
|
|
// Get cheapest unmarked vertex |
62
|
|
|
$cheapestResult = $cheapestVertex->extract(); |
63
|
|
|
$currentVertex = $cheapestResult['data']; |
64
|
|
|
$currentCost = $cheapestResult['priority']; |
65
|
|
|
$currentVertexId = $currentVertex->getId(); |
66
|
|
|
// Vertices can be in the priority queue multiple times, with different path costs (if vertex is already marked, this is an old unvalid entry) |
67
|
|
|
} while (isset($usedVertices[$currentVertexId])); |
68
|
|
|
} |
69
|
|
|
|
70
|
|
|
// Check premature end condition |
71
|
|
|
// If the end vertex is marked as done and the next lowest possible weight is bigger than end vertix, |
72
|
|
|
// we are done processing. |
73
|
|
|
|
74
|
|
|
|
75
|
|
|
// catch "algorithm ends" condition |
76
|
|
|
if ($isEmpty) { |
77
|
|
|
break; |
78
|
|
|
} |
79
|
|
|
|
80
|
|
|
if ($isFirst) { |
81
|
|
|
$isFirst = false; |
82
|
|
|
} |
83
|
|
|
|
84
|
|
|
// mark this vertex |
85
|
|
|
$usedVertices[$currentVertexId] = true; |
86
|
|
|
|
87
|
|
|
if (isset($usedVertices[$endVertexId]) && $totalCostOfCheapestPathTo[$endVertexId] < -$currentCost) { |
|
|
|
|
88
|
|
|
break; |
89
|
|
|
} |
90
|
|
|
|
91
|
|
|
// check for all edges of current vertex if there is a cheaper path (or IN OTHER WORDS: Add reachable nodes from currently added node and refresh the current possible distances) |
92
|
|
|
foreach ($currentVertex->getEdgesOut() as $edge) { |
93
|
|
|
$weight = $edge->getWeight(); |
94
|
|
|
if ($weight < 0) { |
95
|
|
|
throw new UnexpectedValueException('Dijkstra not supported for negative weights - Consider using MooreBellmanFord'); |
96
|
|
|
} |
97
|
|
|
|
98
|
|
|
$targetVertex = $edge->getVertexToFrom($currentVertex); |
99
|
|
|
$targetVertexId = $targetVertex->getId(); |
100
|
|
|
|
101
|
|
|
// if the targetVertex is marked, the cheapest path for this vertex has already been found (no negative edges) { |
102
|
|
|
if (!isset($usedVertices[$targetVertexId])) { |
103
|
|
|
// calculate new cost to vertex |
104
|
|
|
$newCostsToTargetVertex = $totalCostOfCheapestPathTo[$currentVertexId] + $weight; |
105
|
|
|
|
106
|
|
|
if ((!isset($predecesEdgeOfCheapestPathTo[$targetVertexId])) |
107
|
|
|
// is the new path cheaper? |
108
|
|
|
|| $totalCostOfCheapestPathTo[$targetVertexId] > $newCostsToTargetVertex) { |
109
|
|
|
|
110
|
|
|
// Not an update, just a new insert with lower cost |
111
|
|
|
$cheapestVertex->insert($targetVertex, -$newCostsToTargetVertex); |
112
|
|
|
// so the lowest cost will be extracted first |
113
|
|
|
// and higher cost will be skipped during extraction |
114
|
|
|
|
115
|
|
|
// update/set costs found with the new connection |
116
|
|
|
$totalCostOfCheapestPathTo[$targetVertexId] = $newCostsToTargetVertex; |
117
|
|
|
// update/set predecessor vertex from the new connection |
118
|
|
|
$predecesEdgeOfCheapestPathTo[$targetVertexId] = [$edge]; |
119
|
|
|
} elseif ($totalCostOfCheapestPathTo[$targetVertexId] == $newCostsToTargetVertex) { |
120
|
|
|
// Same length paths. We need to add the predecessor to the list of possible predecessors. |
121
|
|
|
$predecesEdgeOfCheapestPathTo[$targetVertexId][] = $edge; |
122
|
|
|
} |
123
|
|
|
} |
124
|
|
|
} |
125
|
|
|
} |
126
|
|
|
|
127
|
|
|
if (!isset($totalCostOfCheapestPathTo[$endVertexId])) { |
128
|
|
|
throw new MultiDijkstraNoPathException("No path found between vertex '".$startVertex->getId()."' and vertex '".$endVertex->getId()."'"); |
129
|
|
|
} |
130
|
|
|
|
131
|
|
|
// algorithm is done, return resulting edges |
132
|
|
|
return $predecesEdgeOfCheapestPathTo; |
133
|
|
|
} |
134
|
|
|
|
135
|
|
|
/** |
136
|
|
|
* @param array<string, Vertex[]> $predecesEdgesArray key is the destination vertex name and value is an array of possible origin vertex |
137
|
|
|
* |
138
|
|
|
* @return Edge\Base[] |
139
|
|
|
*/ |
140
|
|
|
public static function getCheapestPathFromPredecesArray(Vertex $startVertex, Vertex $endVertex, array $predecesEdgesArray) |
141
|
|
|
{ |
142
|
|
|
$edges = []; |
143
|
|
|
$currentVertex = $endVertex; |
144
|
|
|
while ($currentVertex !== $startVertex) { |
145
|
|
|
$predecessorEdges = $predecesEdgesArray[$currentVertex->getId()]; |
146
|
|
|
if (count($predecessorEdges) > 1) { |
147
|
|
|
throw new MultiDijkstraAmbiguityException("There are many possible shortest paths to link vertex '".$startVertex->getId()."' to '".$endVertex->getId()."'"); |
148
|
|
|
} |
149
|
|
|
/* @var $edge \Fhaculty\Graph\Edge\Base */ |
150
|
|
|
$edge = $predecessorEdges[0]; |
151
|
|
|
$edges[] = $edge; |
152
|
|
|
if ($currentVertex === $edge->getVerticesStart()->getVertexFirst()) { |
153
|
|
|
$currentVertex = $edge->getVerticesTarget()->getVertexFirst(); |
154
|
|
|
} else { |
155
|
|
|
$currentVertex = $edge->getVerticesStart()->getVertexFirst(); |
156
|
|
|
} |
157
|
|
|
} |
158
|
|
|
|
159
|
|
|
return array_reverse($edges); |
160
|
|
|
} |
161
|
|
|
|
162
|
|
|
/** |
163
|
|
|
* @param Vertex $startVertex |
164
|
|
|
* @param Vertex $endVertex |
165
|
|
|
* @param array $predecesEdgesArray |
166
|
|
|
* |
167
|
|
|
* @return Edge\Base[][] |
168
|
|
|
*/ |
169
|
|
|
public static function getAllPossiblePathsFromPredecesArray(Vertex $startVertex, Vertex $endVertex, array $predecesEdgesArray) |
170
|
|
|
{ |
171
|
|
|
$edgesPaths = []; |
172
|
|
|
|
173
|
|
|
if ($startVertex === $endVertex) { |
174
|
|
|
return []; |
175
|
|
|
} |
176
|
|
|
|
177
|
|
|
$predecessorEdges = $predecesEdgesArray[$endVertex->getId()]; |
178
|
|
|
|
179
|
|
|
foreach ($predecessorEdges as $edge) { |
180
|
|
|
if ($endVertex === $edge->getVerticesStart()->getVertexFirst()) { |
181
|
|
|
$nextVertex = $edge->getVerticesTarget()->getVertexFirst(); |
182
|
|
|
} else { |
183
|
|
|
$nextVertex = $edge->getVerticesStart()->getVertexFirst(); |
184
|
|
|
} |
185
|
|
|
|
186
|
|
|
$edgesPaths2 = self::getAllPossiblePathsFromPredecesArray($startVertex, $nextVertex, $predecesEdgesArray); |
187
|
|
|
if ($edgesPaths2) { |
188
|
|
|
foreach ($edgesPaths2 as &$edges2) { |
189
|
|
|
$edges2[] = $edge; |
190
|
|
|
} |
191
|
|
|
} else { |
192
|
|
|
$edgesPaths2 = [[$edge]]; |
193
|
|
|
} |
194
|
|
|
|
195
|
|
|
$edgesPaths = array_merge($edgesPaths, $edgesPaths2); |
196
|
|
|
} |
197
|
|
|
|
198
|
|
|
return $edgesPaths; |
199
|
|
|
} |
200
|
|
|
} |
201
|
|
|
|
This check marks PHPDoc comments that could not be parsed by our parser. To see which comment annotations we can parse, please refer to our documentation on supported doc-types.