Elementary shortest path problems, with or without resource constraints
Elementary Shortest Path Problem with or without Resource Constraint
This module solves the shortest elementary path problem, by which it is meant that each node in a path can only be visited once. The edge weights do not need to be nonnegative.
A general algorithm to solve this problem works by solving the k-paths problem, i.e. by solving the nonelementary problem repeatedly and requiring that for the nodes that are present in nonelementary paths multiple (k) paths need to be found. The resulting algorithm finds the elementary shortest paths (ESPP) from a single source to all destinations . However, it can be very slow if very large graphs are considered. Especially because the algorithm in  is intrinsically flawed and the fix for this makes it much less efficient than expected: See note in ESPP.py.
A more powerful algorithm can be used when additionally there are positive resources associated to each edge in the graph, though this requirement can be relaxed to nonnegative resources if not too many edges have zero resource cost associated to them. The resulting Elementary Resource Constrained Shortest Path Problem (ESPPRC) can then be solved by an algorithm that starts by first preprocessing the graph . During preprocessing, the graph is first pruned by discarding all nodes that cannot be reached from the origin or reach the target without violating the resource limits. Subsequently, for each pair in the resulting pruned graph, the least-cost path is sought, which will be used later to incrase so-called dominance of paths by other paths.
After preprocessing, the ESPPRC is solved by introducing node-resources that count how many times a node has been used in a certain path, so that elementarity can be enforced. As this problem can quickly become exponential in nature, it is critical that the numer of resource nodes that have to be introdcuced kept as low as possible. In the algorithm this is done by using a two-fold dominance strategy: On the one hand, paths are only added to the list of possible shortest possible paths, if their resource consumption is not dominated by another path. Also, when a path is added, the other possible paths that are dominated by it, are discarded. On the other hand, strong dominance is a mathematical trick that is used to speed up convergence to the shortest path, by artificially incrementing possible node resources of nodes that cannot be reached . The resulting algorithm seeks the ESPPRC for single source to single destination problems.
pip install pylgrim
- or alternatively
- Clone this repository
pip install -e pylgrim
- matplotlib (for testing)
- random (for testing)
- timeit (for testing)
First of all, the directed graph (
G on which to use pylgrim needs to have a attribute
n_res set through, for example
G = nx.DiGraph(n_res=1)
G.graph['n_res'] = 1
if the graph already exists.
Also, each of the graph edges should have an attribute
weight, which is what will be minimized along the paths.
from pylgrim import tools
should be used to decouple the directed graph to make sure that the source node from which to calculate elementary shortest paths is certain not to have any in-edges, as this is not allowed by the algorithms:
tools.decouple_source(G, source, source_in='source_in')
source node to
source_in and moves all the in-edges to the new node.
This operation happens in situ.
Afterwards, the graph can be reverted to its original state using
tools.undecouple_source(G, source, source_in=`source_in`)
Apart from this,
from pylgrim import Path
should be used to import the
Path class that inherits from the
networkx.digraph and that represents a path that can be cyclic, so it contains the information of source node.
Among other things, the class provides an iterator for the next edge in the path, pretty printing, etc.
Furthermore it is used to represent the solution path(s) from ESPP and ESPPRC.
Elementary shortest path
Elementary shortest path can be used by first importing the module through
from pylgrim import ESPP
The Dynamic labelling algorithm is implemented :
return paths, costs = DLA(G, source, min_K=1, output_pos = False, max_path_len = -1)
G si the directed graph for which to find the shortest path from the
source to all other nodes.
min_K the minimum number of paths to be returned for every node can be changed, with
output_pos also positive costs can be returned and through
max_path_len, the maximum length of a path can be limited, so that the problem becomes easier to solve.
DLA outputs in
paths a dictionary where the keys are the target nodes and for every one of these, a list is returned with the best paths of type
Path, ordered by lowest cost, of length at least
min_K. The cost itself is returned in a matching dictionary
Note, however, that ESPP should no be used for graphs other than very small ones. The following algorithm should be much faster.
Elementary shortest path with resource constraints
Elementary shortest path with resource constraints can be used by first importing the module through
from pylgrim import ESPPRC
The General State Space Augmenting Algorithm,
GSSA is implemented :
best_path, best_path_label = GSSA(G, source, target, max_res, res_min, res_name='res_cost')
G is the directed graph for which to find the shortest path between
target subject to maximum resource usage given in
max_res, the array of maximum resources that can be consumed on any path, where a value should be given for every resource.
Note, furthermore, the graph
G must have been preprocessed first so that it is reduced and has the minimal resource information that is to be passed into
G, res_min = ESPPRC.preprocess(G_no_preprocess, source, target, max_res, res_name='res_cost')
Finally, in these functions,
res_name can be used to specify the name of the resource that is to be used internally for the edges in the graph.
The output of
GSSA is a
best_path of the
Path type and a corresponding
best_path_label that is a 2-tuple with the first element equal to the cost of the best path and the second equal to the array of the resources consumed, less than or equal to
- Run the tests in
[NAME]is the name of the test.
- Plot(s) will appear as well as textual output.
- To debug, change
logging.basicConfig(level=logging.WARNING)present at the top of the python files.
 On the shortest path problem with negative cost cycles by Di Puglia Pugliese, Luigi (DOI: 10.1007/s10589-015-9773-1)
 Accelerated label setting algorithms for the elementary resource constrained shortest path problem by Boland, Natashia (DOI: 10.1016/j.orl.2004.11.011)
Copyright 2018 Toon Weyens, Daan van Vugt