Parallel Shortest-Paths Using Radius Stepping

TL;DR

This paper introduces Radius-Stepping, a parallel algorithm for single-source shortest path problems that achieves nearly-linear work and improved depth bounds by using variable radius steps and graph preprocessing.

Contribution

It presents Radius-Stepping, a novel $ riangle$-stepping-like algorithm with theoretical bounds on steps and work, improving parallel SSSP efficiency.

Findings

Achieves $O((m+n ext{log}n) ext{log}(n/ ho))$ work.

Bounds the number of steps by $O(rac{n}{ ho} ext{log} ho L)$.

Preprocessing can be done efficiently with bounded additional edges.

Abstract

The single-source shortest path problem (SSSP) with nonnegative edge weights is a notoriously difficult problem to solve efficiently in parallel---it is one of the graph problems said to suffer from the transitive-closure bottleneck. In practice, the $\Delta$-stepping algorithm of Meyer and Sanders (J. Algorithms, 2003) often works efficiently but has no known theoretical bounds on general graphs. The algorithm takes a sequence of steps, each increasing the radius by a user-specified value $\Delta$. Each step settles the vertices in its annulus but can take $\Theta(n)$ substeps, each requiring $\Theta(m)$ work ($n$ vertices and $m$ edges). In this paper, we describe Radius-Stepping, an algorithm with the best-known tradeoff between work and depth bounds for SSSP with nearly-linear ($\otilde(m)$) work. The algorithm is a $\Delta$-stepping-like algorithm but uses a variable instead of fixed-size increase in radii, allowing us to prove a bound on the number of steps. In particular, by using what we define as a vertex $k$-radius, each step takes at most $k+2$ substeps. Furthermore, we define a $(k, \rho)$-graph property and show that if an undirected graph has this property, then the number of steps can be bounded by $O(\frac{n}{\rho} \log \rho L)$, for a total of $O(\frac{kn}{\rho} \log \rho L)$ substeps, each parallel. We describe how to preprocess a graph to have this property. Altogether, Radius-Stepping takes $O((m+n\log n)\log \frac{n}{\rho})$ work and $O(\frac{n}{\rho}\log n \log (\rho{}L))$ depth per source after preprocessing. The preprocessing step can be done in $O(m\log n + n\rho^2)$ work and $O(\rho^2)$ depth or in $O(m\log n + n\rho^2\log n)$ work and $O(\rho\log \rho)$ depth, and adds no more than $O(n\rho)$ edges.