---
title: 'BMSSP: Bounded Multi-Source Shortest Path'
url: https://www.emergentmind.com/topics/bounded-multi-source-shortest-path-bmssp
type: topic
---

# BMSSP: Bounded Multi-Source Shortest Path

Searching arXiv for BMSSP and related shortest-path papers to ground the article in the cited literature.
Bounded Multi-Source Shortest Path (BMSSP) is the problem of computing distances from a designated source set \(S \subseteq V\) to all vertices of a graph, typically requiring either exact distances \(\mathrm{dist}(s,v)\) or estimates \(\tilde d(s,v)\) satisfying \(d_G(s,v)\le \tilde d(s,v)\le (1+\epsilon)d_G(s,v)\) for all \((s,v)\in S\times V\). The problem appears under several closely related names. In the HYBRID literature, the “Bounded Multi-Source Shortest Paths” problem on \(n\) nodes with \(k\) sources is exactly \(k\)-SSP, while in centralized and parallel settings it is often presented as MSSP or as \(S\times V\) shortest paths [2304.07107]. Recent work places BMSSP at the intersection of hopsets, emulators, rectangular matrix multiplication, low-diameter decompositions, and communication-constrained distributed computation, yielding near-tight or optimal bounds in centralized, PRAM, Congested Clique, and HYBRID models [2004.07572, 2003.03058, 2606.26554].

## 1. Problem formulation and model-dependent variants

BMSSP is defined on a graph \(G=(V,E)\) with a distinguished source set \(S\). In the weighted approximation setting, the goal is to compute \((1+\epsilon)\)-approximate distances for all pairs in \(S\times V\). In the exact undirected unweighted setting, the task is to compute \(\mathrm{dist}(s,v)\) for every \(s\in S\) and \(v\in V\). The reachability analogue replaces finite distances by a Boolean distinction between finite and infinite distance. These formulations are all explicit in the recent BMSSP/MSSP literature [2004.07572, 2606.26554].

The source-set size is usually parameterized either as \(|S|=n^r\) with \(0<r\le 1\), or as \(|S|=n^\sigma\) with \(\sigma\in[0,1]\), depending on the paper. In HYBRID, the same problem is written with \(k=|S|\), and every node \(v\in V\) must learn \(d(v,s)\) for all \(s\in S\), either exactly or up to a stretch \(\alpha\ge 1\) [2304.07107].

A plausible implication is that “bounded” is not used with a single universal formal restriction across all papers. Rather, it denotes a multi-source shortest-path regime in which the number of sources is treated as a parameter of the complexity analysis, such as \(k\), \(n^r\), or \(n^\sigma\). This explains why the same underlying task is alternately discussed as BMSSP, MSSP, or \(k\)-SSP in different computational models.

## 2. Hopsets, emulators, and distance products as a general algorithmic template

A central line of work reduces BMSSP to repeated approximate distance products between an \(|S|\times n\) matrix and an \(n\times n\) adjacency matrix. In the weighted undirected setting, the framework of hopsets plus rectangular fast matrix multiplication gives a \((1+\epsilon)\)-approximation for all distances from \(|S|=n^r\) sources. For any integer \(\kappa\ge 1\), one can build a hopset \(H\) of \(O(n^{1+1/\kappa})\) edges in time \(\tilde O(|E|\cdot n^{1/\kappa})\), with hopbound \(\beta=\bigl(\frac{\kappa}{\epsilon}\bigr)^{O(\kappa)}\), and then perform \(\beta-1\) relaxation rounds based on approximate distance products. With the choice \(\kappa=(\log\log n)/\log\log\log n\), this yields total centralized running time
\[
\tilde O\Bigl(|E|\cdot n^{o(1)} + n^{\omega(r)}\Bigr),
\]
where \(n^{\omega(r)}\) is the time for multiplying an \(n^r\times n\) matrix by an \(n\times n\) matrix [2004.07572].

The same template extends to PRAM. The hopset can be built in parallel time \((\log n/\epsilon)^{O(\kappa)}\) with work \(\tilde O(|E|\cdot n^{1/\kappa})\), and each approximate distance product takes \(O(R\log n)\) parallel time with \(\tilde O(R\,n^{\omega(r)})\) work. Tuning \(\kappa=1/\rho\) yields PRAM time \((\log n)^{O(1/\rho)}\) and work \(\tilde O(|E|\,n^\rho+n^{\omega(r)})\) for any arbitrarily small constant \(\rho>0\) [2004.07572].

In the Congested Clique, the same conceptual decomposition remains useful but requires a distributed rectangular matrix-multiplication primitive. The resulting bound is
\[
T_{\rm CC}(r)=\tilde O\!\bigl(n^{\,1-\tfrac{2}{\omega(r')}}\bigr),
\quad r'=1-(1-r)\omega(r').
\]
In particular, for \(r\le 0.655\), the paper obtains polylogarithmic round complexity in weighted graphs, and for unweighted graphs with \(r\le 0.655\), it achieves \(O(\log\log n)\) rounds with high probability [2004.07572].

This family of results establishes BMSSP as an intermediate problem between single-source and all-pairs shortest paths. The performance varies smoothly with the source exponent \(r\), rather than forcing an immediate transition from repeated SSSP to full APSP.

## 3. Deterministic BMSSP in the Congested Clique

For unweighted undirected graphs and \(k=O(\sqrt n)\) sources, a separate deterministic line of work obtains a poly\((\log\log n)\)-round Congested-Clique algorithm. The target guarantee is \((1+\epsilon)\)-approximate distances \( \delta(s,v) \) satisfying
\[
d_G(s,v)\le \delta(s,v)\le (1+\epsilon)\,d_G(s,v).
\]
The algorithm separates source-target pairs into “long” distances and “short” distances around a threshold
\[
t=\Theta(\beta/\epsilon),
\qquad
\beta=O\bigl((\log\log n/\epsilon)^{\log\log n}\bigr),
\]
and handles the two regimes by different mechanisms [2003.03058].

For long distances, it constructs a sparse near-additive emulator \(H\) with \(O(n\log\log n)\) edges and \((1+\epsilon,\beta)\)-stretch. The emulator is built from a hierarchy
\[
S_0 = V \supseteq S_1 \supseteq \cdots \supseteq S_r \supseteq S_{r+1}=\emptyset
\]
with sampling probabilities \(p_{i+1}=n^{-2^i/2^r}\) for \(i<r\), together with radii \(\delta_i\) and cumulative budgets \(R_i\). The stretch analysis proves
\[
d_H(u,v)\le (1+20\epsilon i)\,d_G(u,v)+\beta_i,
\]
which after rescaling yields the stated \((1+\epsilon,\beta)\)-emulator [2003.03058].

For short distances, the algorithm builds a bounded hopset \(H'\) for radius \(t\), forms \(G''=G\cup H'\), and then runs \(k\)-source detection from the source set \(S\) with hop-bound \(\beta'\). The output is a \((1+\epsilon/2)\)-approximation for all source-target pairs of true distance at most \(t\), and the final answer takes the minimum of the long-regime emulator estimate and the short-regime hopset estimate [2003.03058].

A notable feature is derandomization. Every random hitting-set step is replaced by a soft-hitting-set routine that can be derandomized in \(O((\log\log n)^3)\) rounds, using a pseudorandom generator of seed \(O(\log N\cdot(\log\log N)^3)\) bits and the method of conditional expectations. The final outcome is a deterministic \(O(\mathrm{poly}(\log\log n))\)-round BMSSP algorithm for \(k=O(\sqrt n)\), improving on previous deterministic poly\((\log n)\)-round bounds in the Congested Clique [2003.03058].

## 4. Near-tight BMSSP in the HYBRID model

In the synchronous message-passing HYBRID model, nodes have two communication modes per round: local communication along graph edges with unbounded message size, \(\lambda=\infty\), and global communication between arbitrary pairs with total per-node bandwidth \(\gamma=\Theta(\log^2 n)\) bits. The notation is
\[
\mathrm{HYBRID}=(\lambda=\infty,\gamma=\Theta(\log^2 n)).
\]
Within this model, BMSSP is exactly \(k\)-SSP: given an undirected weighted graph \(G=(V,E)\) and a set \(S\subseteq V\) of \(|S|=k\) sources, every node must learn its distance to all sources, either exactly or approximately [2304.07107].

The main upper bound is a randomized algorithm that, with high probability, computes a \((1+\epsilon)\)-approximation to \(k\)-SSP in
\[
T=\tilde O(\sqrt k)=O(\sqrt k\cdot \mathrm{polylog}\, n)
\]
rounds for any \(k\in[n]\) and any constant \(\epsilon>0\). More generally, if the global bandwidth is \(\gamma\) bits per node per round, the complexity becomes
\[
T=\tilde O\bigl(\sqrt{k/\gamma}\bigr)=O\bigl(\sqrt{k/\gamma}\cdot \mathrm{polylog}\, n\bigr).
\]
This matches the known lower bound up to polylogarithmic factors for the full range of \(k\) [2304.07107].

The algorithm combines several ingredients. First, it simulates the Minor-Aggregation interface in HYBRID in \(O(\mathrm{polylog}\,n)\) rounds and implements a PRAM-based Eulerian-Orientation oracle via network decompositions and forest-cover reductions. This yields an \((1+\epsilon)\)-approximate SSSP routine in \(T_0=O(\epsilon^{-2}\cdot\mathrm{polylog}\,n)\) rounds. Second, it samples a skeleton set \(V_S\) by selecting each node independently with probability \(p=1/x\), where \(x\approx \sqrt{k/\gamma}\). A weighted skeleton graph \(S\) is then built on \(V_S\), inserting \((u,v)\) whenever \(\mathrm{hop\mbox{-}distance}_G(u,v)\le h=O(x)\) and giving that edge weight \(d_{h,G}(u,v)\). In \(O(h)\) rounds, the construction ensures that for all \(u,v\in V_S\), \(d_S(u,v)=d_G(u,v)\), and no shortest-path segment of at least \(h\) hops avoids the sampled skeleton [2304.07107].

Third, for each skeleton node \(u\in V_S\), the algorithm computes a helper set \(H_u\subseteq V\) of size \(\Theta(x)\) within hop-radius \(x\), such that each real node participates in \(O(1)\) sets. The \(k\) SSSP instances are distributed across helpers, with each helper simulating \(O(k/x)=O(\sqrt{k\gamma})\) instances. By the scheduling theorem in the paper, all instances on the skeleton can be completed in \(O(x\cdot T_0)\) rounds. Finally, every node \(v\in V\) fetches approximate distances \(d_S(u,s)\) from skeleton nodes \(u\) within \(h\) hops and outputs
\[
\hat d(v,s)=\min_{u\in V_S\cap B_h(v)}\bigl[d_G(v,u)+d_S(u,s)\bigr].
\]
By the sampling lemma and the skeleton stretch guarantee, this gives overall stretch \((1+\epsilon)\) for random sources and stretch \((3+\epsilon)\) for adversarial sources [2304.07107].

The lower bound, due to Kuhn–Schneider, shows that any HYBRID algorithm, even randomized and even for huge stretch, requires
\[
T=\tilde \Omega(\sqrt{k/\gamma})
\]
rounds. The reduction encodes \(k\)-bit inputs across \(k\) disjoint source sets whose pairwise distances differ by \(0\) or \(1\) hop, and yields the sketch inequality
\[
\gamma\cdot T + T^2 \ge \Omega(k/\mathrm{polylog}\,n),
\]
which implies \(T=\tilde\Omega(\sqrt{k/\gamma})\) [2304.07107].

These bounds close the earlier gap between upper bounds such as \(O(n^{1/3}+\sqrt k)\) or \(O(n^{1/2})\) and the lower bound for small \(k\). Consequences stated explicitly include \(\tilde O(\sqrt n)\)-round APSP when \(k=n\), \(\tilde O(1)\)-round SSSP when \(k=1\), exact \((1+o(1))\) solutions for \(k\) random sources, and \((3+\epsilon)\)-approximate solutions for \(k\) adversarial sources in the same time bound [2304.07107].

## 5. Exact MSSP through low-diameter small-neighborhood decomposition

A recent exact result for undirected unweighted graphs gives an almost optimal BMSSP/MSSP algorithm running in
\[
\widetilde O\bigl(n^{\omega(\sigma,1,1)}\bigr)
\]
time for \(|S|=n^\sigma\). This matches the Boolean-matrix-multiplication lower bound and interpolates smoothly between \(O(n^2)\) for single-source and \(O(n^\omega)\) for all-pairs shortest paths [2606.26554].

The key structural ingredient is a \((d,\phi)\)-decomposition, a partition
\[
V=V_1\sqcup V_2\sqcup\cdots\sqcup V_k
\]
such that \(\max_i \mathrm{diam}(G[V_i])\le d\) and \(\sum_i |N(V_i)|\le \phi n\), where \(N(V_i)=\{y\notin V_i\mid \exists x\in V_i,(x,y)\in E\}\). For any \(\psi>1\), the paper computes such a decomposition in time
\[
O\bigl(n^2\psi^4\log n\log_\psi n\bigr)
\]
with
\[
d=6\cdot 3^{\log_\psi n},
\qquad
\phi=O(\psi\log n\log_\psi n).
\]
Setting \(\psi=2^{\sqrt{\log n}}\) gives \(d,\phi=n^{o(1)}\) in time \(n^{2+o(1)}\) [2606.26554].

The BMSSP algorithm maintains, for each source \(s\), a current level \(\ell(s)\), the frontier \(\widehat L(s,\ell)\) of vertices at exact distance \(\ell\), and an active bucket of not-yet-explored frontier vertices. Each cluster \(V_i\) carries a score equal to the number of sources whose active set intersects \(V_i\). The algorithm repeatedly selects a cluster with
\[
\mathrm{score}(i)>
\frac{|S|\cdot |N(V_i)|}{2\phi n}.
\]
If \(|N(V_i)|\) is small, it performs direct relaxation. If \(|N(V_i)|\) is large, it constructs Boolean matrices
\[
X:S'\times V_i,
\qquad
Y:V_i\times N(V_i),
\]
and uses their product \(XY\) to update many source-boundary pairs simultaneously [2606.26554].

The complexity analysis relies on three combinatorial facts stated in the paper: each source touches a cluster at most \(d\) times; each cluster is selected at most \(O(d\phi n\log n/|N(V_i)|)\) times; and a convexity lemma collapses the total matrix-multiplication cost over all relevant iterations. After bucketing by \(|N(V_i)|\) and applying the rectangular matrix-multiplication exponent \(\omega(a,1,1)\), the running time becomes \(\widetilde O(n^{\omega(\sigma,1,1)})\) [2606.26554].

The paper also states several extensions. For integer weights in \([1..B]\), bounded-weight MSSP can be solved in \(\widetilde O(B\,n^{\omega(\sigma,1,1)})\) time. Multiple-source reachability in directed graphs can be solved in \(\widetilde O(T(|S|,n,n))\) time, and the same divide-and-conquer approach extends to MSSP on DAGs. Applications include construction of a linear-size \(\beta\)-hop-set of size \(O(n^2/\beta^2)\) in \(\tilde O(n^{\omega(1/2,1,1)})=O(n^{2.043})\), a \(+4\)-emulator of size \(O(n^{4/3})\) in \(O(n^{\omega(2/3,1,1)})=O(n^{2.132})\), and an \(O(n^{2.084})\)-time algorithm for a \(\widetilde O(n)\)-size shortcut set reducing diameter to \(O(n^{1/3})\) [2606.26554].

## 6. Complexity landscape and conceptual significance

The recent BMSSP literature shows that the problem does not admit a single canonical complexity profile; instead, the dominant technique depends strongly on the computational model, graph class, approximation target, and source-set regime.

| Setting | Regime | Stated result |
|---|---|---|
| Centralized / PRAM / Congested Clique | \(|S|=n^r\), weighted undirected | \((1+\epsilon)\)-approximate BMSSP via hopsets and rectangular matrix multiplication [2004.07572] |
| Congested Clique | \(k=O(\sqrt n)\), unweighted undirected | Deterministic poly\((\log\log n)\)-round \((1+\epsilon)\)-approximate BMSSP [2003.03058] |
| HYBRID | \(k\in[1..n]\), weighted undirected | \(\tilde O(\sqrt{k/\gamma})\)-round \(k\)-SSP, near-tight up to polylogarithmic factors [2304.07107] |
| Centralized exact | \(|S|=n^\sigma\), undirected unweighted | \(\widetilde O(n^{\omega(\sigma,1,1)})\)-time MSSP matching the BMM lower bound [2606.26554] |

One recurring misconception is that BMSSP is merely repeated SSSP with bookkeeping. The current results show a more intricate picture. In centralized exact computation, the trivial bounds are to run BFS/SSSP from each of the \(n^\sigma\) sources in \(O(n^2 n^\sigma)\) time or run APSP in \(\tilde O(n^\omega)\) time, but the 2026 algorithm instead matches the rectangular-BMM lower bound \(\widetilde O(n^{\omega(\sigma,1,1)})\) [2606.26554]. In distributed settings, communication structure rather than algebraic complexity becomes decisive: the HYBRID lower bound is governed by the joint effect of global bandwidth and local congestion, while Congested Clique algorithms exploit global all-to-all communication through emulators, hopsets, and distributed matrix multiplication [2304.07107, 2004.07572].

A second important distinction is between exactness and approximation. The centralized 2026 result is exact for undirected unweighted graphs, whereas the 2020 Congested-Clique and weighted centralized/PRAM results are \((1+\epsilon)\)-approximate, and the HYBRID result gives \((1+\epsilon)\)-approximation in general, together with exact \((1+o(1))\) behavior for random sources and \((3+\epsilon)\)-approximation for adversarial sources in the stated consequence for BMSSP [2606.26554, 2003.03058, 2004.07572, 2304.07107].

Taken together, these developments identify BMSSP as a boundary problem linking shortest paths, emulators, hopsets, matrix multiplication, and distributed lower bounds. The modern view is not merely that BMSSP generalizes SSSP, but that it exposes intermediate regimes in which neither repeated single-source computation nor full APSP is the right abstraction.

Source: https://www.emergentmind.com/topics/bounded-multi-source-shortest-path-bmssp