---
title: Near-Exact Hopsets in Graph Algorithms
url: https://www.emergentmind.com/topics/near-exact-hopsets
type: topic
---

# Near-Exact Hopsets in Graph Algorithms

Near-exact hopsets are graph augmentations that add weighted shortcut edges to a graph so that every pair of vertices admits a path with a bounded number of hops and length within a factor \(1+\epsilon\) of the true shortest-path distance. In the standard terminology, a \((\beta,\epsilon)\)-hopset enforces a hop-bound \(\beta\) and multiplicative stretch \(1+\epsilon\), with the regime of small \(\epsilon\) constituting the near-exact setting. Introduced by Cohen in the context of parallel computation of approximate shortest paths, hopsets later became a basic primitive in dynamic graph algorithms, distributed computing, and the streaming model [1605.04538].

## 1. Formal definition and variants

For a weighted undirected graph \(G=(V,E,\omega)\), let \(d_G(u,v)\) denote the shortest-path distance, and for an integer \(\beta \ge 1\), let
\[
d_{G'}^{(\beta)}(u,v)
\]
denote the length of a shortest \(u\)-to-\(v\) path in \(G'=(V,E\cup H)\) that uses at most \(\beta\) edges. A graph \(H=(V,H,\omega_H)\) is a \((\beta,\epsilon)\)-hopset for \(G\) if, for every \(u,v\in V\),
\[
d_G(u,v)\;\le\; d_{G'}^{(\beta)}(u,v)\;\le\;(1+\epsilon)\,d_G(u,v).
\]
Here \(\beta\) is the hop-bound and \(\epsilon\) is the multiplicative stretch parameter [1605.04538].

The same basic idea extends to directed graphs, but the formulation is less uniform. One directed definition requires exact distance preservation in the augmented graph together with a low-hop approximate path: for every \(u,v\in V\),
\[
\dist_{G\cup H}(u,v)=\dist_G(u,v),
\]
and there exists a \(u\!\to\!v\) path in \(G\cup H\) of at most \(\beta\) edges and weight at most \((1+\epsilon)\dist_G(u,v)\) [2207.04507]. A broader directed formulation fixes \(\alpha \ge 1\) and asks that every reachable ordered pair \((s,t)\) admit an \(s\to t\) path in \(G\cup H\) of at most \(\beta\) edges and total weight at most \(\alpha\,\mathrm{dist}_G(s,t)\); the case \(\alpha=1+\epsilon\) is the directed near-exact regime [2411.08151].

Two parameters dominate the theory. The first is sparsity, usually expressed as \(|H|=O(n^{1+1/\kappa})\) or by specialized linear-size bounds. The second is the hop-bound \(\beta\), whose dependence on \(\epsilon\), the sparsity exponent, and in directed settings the aspect ratio, is the central structural quantity throughout the literature.

## 2. Undirected trade-offs and existential bounds

The undirected theory developed through a sequence of progressively sharper size-hopbound trade-offs. An early milestone gave linear-size hopsets with substantially improved hopbound: setting \(\kappa=\log n\) in the Elkin–Neiman construction yields \(|H|=O(n)\) and
\[
\beta=(\log n)^{\log^{(3)}n+O(1)},
\]
while the general bound is \(|H|=O(n^{1+1/\kappa})\) with
\[
\beta=O\!\Bigl(\tfrac{\log\kappa+1/\rho}{\epsilon}\Bigr)^{\log\kappa+1/\rho+1}
\]
for a runtime parameter \(\rho\) [1704.08468].

A subsequent construction exhibited the first sparse hopsets with a constant number of hops. For parameters \(\epsilon\in(0,1)\), \(\kappa\in[2,\Theta(\log n)]\), and \(\rho\in[1/\kappa,1/2)\), one can construct a \((\beta,\epsilon)\)-hopset of expected size
\[
|H|=O\bigl(n^{1+1/\kappa}\log n\bigr),
\]
where
\[
\beta=O\!\Bigl(\tfrac{\log\kappa+1/\rho}{\epsilon}\Bigr)^{\log\kappa+O(1/\rho)},
\]
and \(\beta\) is independent of \(n\) for fixed \(\epsilon,\kappa,\rho\); the leading \(\log n\) in the size can be removed at the expense of increasing the exponent of \(\beta\) by \(+1\) [1605.04538].

The later survey literature emphasizes a parallel line of results based on scale-free constructions. For every \(n\)-vertex weighted graph, every \(\epsilon\in(0,1/6)\), and integer \(\kappa\ge 1\), there exists a \((1+\epsilon,\beta)\)-hopset of size
\[
O\bigl(n^{1+1/\kappa}\log n\bigr)
\]
with
\[
\beta = O\!\Bigl(\tfrac{\log \kappa}{\epsilon}\Bigr)^{\log \kappa-2},
\]
and a refined analysis removes the leading \(\log n\). The same parameters can be achieved by a single, \(\epsilon\)-independent construction, so that one hopset simultaneously satisfies the \((1+\epsilon,\beta)\) guarantee for all \(\epsilon\in(0,1/6)\) [2001.07477].

| Construction | Size | Hop-bound |
|---|---:|---:|
| Linear-size Elkin–Neiman specialization [1704.08468] | \(O(n)\) | \((\log n)^{\log^{(3)}n+O(1)}\) |
| Constant-hop hopsets [1605.04538] | \(O(n^{1+1/\kappa}\log n)\) | \(O\!\bigl(\frac{\log\kappa+1/\rho}{\epsilon}\bigr)^{\log\kappa+O(1/\rho)}\) |
| Scale-free hopsets [2001.07477] | \(O(n^{1+1/\kappa})\) | \(\Theta\!\bigl(\frac{\log\kappa}{\epsilon}\bigr)^{\log\kappa-2}\) |
| Thorup–Zwick optimal trade-off [1705.00327] | \(O(n^{1+1/(2^{k+1}-1)})\) | \(O((k/\epsilon)^k)\) |

These parameterizations are not identical: \(\kappa\) and \(k\) index different sparsity curves. Nonetheless, they chart the main undirected progression from linear-size quasi-polylogarithmic hopbound, to constant-hop sparse constructions, to sparsity-optimal Thorup–Zwick-based hopsets.

A further structural synthesis appears in the unified framework of Neiman and Shabat. In the size regime \(O(n^{1+1/k})\), their sampling-and-bunch formalism recovers a simultaneous \((1+\epsilon,\beta)\)-hopset with
\[
\beta = O\bigl((\log k)/\epsilon\bigr)^{\log k},
\]
and proves a matching lower bound of a different form: any \((\alpha,\beta)\)-hopset of size \(O(n^{1+1/k})\) must satisfy
\[
\alpha\cdot \beta \ge \Omega(k)
\]
[2108.09673]. For near-exact hopsets, where \(\alpha=1+\epsilon\), this gives a coarse but general obstruction to simultaneously minimizing multiplicative stretch and hopbound.

## 3. Thorup–Zwick emulators as universally optimal hopsets

A decisive simplification came from the observation that Thorup–Zwick sublinear-additive emulators are also hopsets. In the construction, one fixes an integer \(k\ge 1\), forms nested samples
\[
V=V_0\supseteq V_1\supseteq \cdots \supseteq V_k \supseteq V_{k+1}=\emptyset,
\]
and for each \(v\in V_i\setminus V_{i+1}\) adds weighted edges from \(v\) to every \(u\in B(v)\cup\{p_{i+1}(v)\}\), where
\[
B(v)=\{u\in V_i:\dist_G(v,u)<\dist_G(v,p_{i+1}(v))\},
\]
and each edge \((v,u)\) receives weight \(w(v,u)=\dist_G(v,u)\) [1705.00327].

Following Pettie, the sampling probabilities are set to
\[
q_i=n^{-(2^i-1)/(2^{k+1}-1)} \qquad i=0,\dots,k.
\]
With this choice,
\[
|H| = O\bigl(n^{1+1/(2^{k+1}-1)}\bigr).
\]
The key analytic step is to convert Thorup–Zwick’s sublinear additive emulator guarantee into multiplicative \((1+\epsilon)\)-stretch by explicitly bounding the number of hops in \(G\cup H\). Defining \(h_0=1\) and
\[
h_i=(r+1)h_{i-1}+r \qquad i=1,\dots,k,
\]
one obtains \(h_k<2(r+1)^k\), and choosing \(r=\Theta(k/\epsilon)\) yields
\[
\beta = h_k = O((k/\epsilon)^k)
\]
[1705.00327].

This trade-off is described in the paper as universally optimal. Abboud–Bodwin–Pettie proved that any hopset of size \(n^{1+1/(2^{k+1}-1)-\delta}\) must have
\[
\beta=\Omega(1/\epsilon^{\,k+1}),
\]
whereas the Thorup–Zwick modification achieves size \(n^{1+1/(2^{k+1}-1)}\) with
\[
\beta=O((k/\epsilon)^k),
\]
matching the lower-bound curve up to \(k\)-polynomial factors [1705.00327].

The same construction also has corollaries outside the hopset setting. For unweighted graphs it gives a sublinear-additive emulator of size \(O(n^{1+1/(2^{k+1}-1)})\) and stretch
\[
f(d)=d + (4+o(1))\cdot k\cdot d^{1-1/k},
\]
and, combined with standard “stretch \(\approx 1\)” pruning, yields a spanner with stretch \((1+\epsilon)\), additive error \(O((k/\epsilon)^{k-1})\), and size
\[
O((k/\epsilon)^h \cdot n^{1+1/(2^{k+1}-1)}),
\]
where \(h<0.75\) [1705.00327]. This is one of the clearest demonstrations that near-exact hopsets, near-additive emulators, and sparse spanners often occupy the same combinatorial template.

## 4. Construction paradigms and computational models

The first constant-hop construction is also notable for its algorithmic generality. Its centralized algorithm uses a hierarchy of clusterings, interleaving super-clustering and interconnection over distance scales. For fixed \(\epsilon,\kappa,\rho\), it constructs in
\[
O((|E|+n\log n)\cdot n^\rho/\rho\cdot \log n)
\]
time a \((\beta,\epsilon)\)-hopset of expected size \(O(n^{1+1/\kappa}\log n)\), with space \(O(|E|+n)\); parallel, Congested Clique, CONGEST, and streaming implementations with explicit resource bounds are also given [1605.04538].

In the EREW PRAM model, this framework yields a construction using
\[
O((|E|+n^{1+1/\kappa}\log^2 n)n^\rho\log^2 n)
\]
processors and time
\[
O(\beta(\log \kappa+1/\rho)\log^2 n),
\]
after which \((1+\epsilon)\)-approximate shortest paths from \(s\) sources can be computed via \(\beta\)-limited Bellman–Ford in \(O(\beta\log n)\) further time [1605.04538]. In the Congested Clique, one can build the hopset in
\[
O(n^\rho/\rho\cdot \log n\cdot \beta\cdot \log n\cdot 2^t)
\]
rounds, with a choice of \(t\approx \rho\log n\) making \(\beta\) constant in \(n\) [1605.04538].

Deterministic parallelization required a separate advance. The first efficient deterministic parallel hopset algorithm builds a \((1+\epsilon,\beta)\)-hopset of size
\[
O(n^{1+\rho}\log n)
\]
with
\[
\beta=(\log n)^{O(1/\rho)},
\]
in CREW-PRAM time
\[
O((\log A)(\log_{n^\rho}n)\log^2 n)
\]
and work
\[
O((|E|+n^{1+\rho})n^\rho\log n),
\]
yielding the first deterministic polylogarithmic-time PRAM algorithm for \((1+\epsilon)\)-SSSP with slightly super-linear work [2009.14729].

Distributed settings impose different sparsity-round trade-offs. In the Congested Clique, a randomized algorithm constructs, with high probability, a \((\beta,\epsilon)\)-hopset of size
\[
\tilde O\bigl(n^{1+1/(2k)}\bigr),
\]
hopbound
\[
O\!\Bigl(\tfrac{\log^2 n}{\epsilon}\,\Bigl(\tfrac{\log n\cdot\log k}{\epsilon}\Bigr)^{\log(k+1)-1}\Bigr),
\]
and round complexity
\[
O(\beta\log^2 n),
\]
providing the first polylogarithmic-round Congested Clique algorithm for sparse hopsets with polylogarithmic hopbound [1911.07154].

Taken together, these results show that the near-exact regime is not only an existential theory of sparse augmentations. It is also a design framework for sequential, PRAM, Congested Clique, CONGEST, streaming, and deterministic shortest-path algorithms.

## 5. Directed near-exact hopsets

Directed graphs are markedly harder. The first nearly linear-work directed construction produces a hopset with \(\tilde O(n)\) edges and
\[
\beta=n^{1/2+o(1)},
\]
in sequential time \(\tilde O(m)\), and a parallel version with work \(\tilde O(m)\) and span \(n^{1/2+o(1)}\) [1912.05506]. This already established that near-exact directed hopsets can be built efficiently, but not yet with the same linear-size hopbound known for directed shortcut sets.

That gap was closed by Bernstein and Wein. For any \(n\)-vertex directed graph with integer weights in \([1,W]\), fixed \(\epsilon\in(0,1)\), and \(\beta\ge 20\log n\), there exists a \((\beta,\epsilon)\)-hopset \(H\) of size
\[
|H|
=\begin{cases}
O\!\bigl(n^2\,\log^7 n\,\log^2(nW)\bigr)\,/(\epsilon^2\beta^3) & \text{if }\beta\le n^{1/3},\\[4pt]
O\!\bigl(n^{3/2}\,\log^7 n\,\log^2(nW)\bigr)\,/(\epsilon^2\beta^{3/2}) & \text{if }\beta> n^{1/3}.
\end{cases}
\]
In particular, choosing
\[
\beta=\Theta\bigl(\epsilon^{-2}n^{1/3}\log^2 n\bigr)
\]
gives \(|H|=O(n)\) and hop-bound \(\tilde O(n^{1/3})\). More generally, for every \(\delta\in(0,1)\), one can set
\[
\beta=O\bigl(n^{1-\delta}\polylog n\bigr)
\]
to obtain a \((\beta,\epsilon)\)-hopset of size
\[
O\bigl(n^{1+\delta}\polylog n\bigr),
\]
matching, up to polylogarithmic factors, the known shortcut-set trade-off [2207.04507].

The construction combines several ingredients: a weighted transitive closure, a “nice path collection” of vertex-disjoint shortest paths of exactly \(h_0=\beta/(162\log n)\) hops, Raz–Spira forward 2-hop emulators on those paths, hierarchical vertex-path sampling, and a new sixth-hop backward-shortcutting gadget that handles the cyclic obstructions absent in DAG-like settings [2207.04507].

Lower bounds reveal that the directed near-exact regime is highly sensitive to weight structure. For arbitrary weighted directed graphs, every \(O(n)\)-size approximate hopset with any given finite stretch must have
\[
\beta=\Omega(n^{1/2}),
\]
while \(O(m)\)-size exact hopsets in unweighted directed graphs require \(\beta=\Omega(n^{2/7})\), and \(O(m)\)-size shortcut sets require \(\beta=\Omega(n^{2/9})\) [2411.08151]. The \(\Omega(n^{1/2})\) bound establishes a separation from the polynomial-aspect-ratio setting, where \(O(n)\)-size approximate hopsets with hopbound \(\widetilde O(n^{1/3})\) are known [2411.08151]. A common misunderstanding is therefore to treat the \(\widetilde O(n^{1/3})\) linear-size bound as fully general in directed graphs; the lower bound shows that unbounded aspect ratio rules this out.

## 6. Reductions, unification, and optimization

Near-exact hopsets have long been studied alongside near-additive spanners and emulators. The survey literature highlighted that the results of Elkin–Peleg for spanners and Elkin–Neiman for hopsets are “strikingly similar,” and explicitly posed the problem of a black-box reduction between near-additive emulators and near-exact hopsets [2001.07477].

One direction was clarified through a reduction from hopsets to spanners, emulators, and distance preservers. The central notion is an \(r\)-missing \(t\)-spanner: a subgraph \(S\subseteq G\) such that every pair \(u,v\) has an approximate path of length at most \(t\,d_G(u,v)\) omitting at most \(r\) edges of \(S\). Given a sequence of hopsets \(H_1,\dots,H_\ell\) with parameters \((\beta_i,\epsilon)\), one can build in polynomial time a subgraph \(S\subseteq G\) that is an
\[
r=\beta_\ell \quad\text{missing}\quad t=(1+\epsilon)^\ell
\]
spanner of size
\[
|S|\le \sum_{i=1}^{\ell}|H_i|\cdot \beta_{i-1}.
\]
The same framework yields new weighted near-additive spanners, directed weighted pairwise preservers, and lower bounds such as \(\beta=\Omega(n^{1/3})\) for linear-size exact weighted hopsets [2211.06920].

The reverse direction was established more recently. Any algorithm \(A\) that constructs a \((1+\epsilon',\beta)\)-emulator of size \(S_A(n,\epsilon',\beta)\) for undirected unweighted graphs can be used as a black box to obtain, for any \(0<\epsilon\le 1\), a \((1+\epsilon,H)\)-hopset for an undirected weighted graph with
\[
H = O\!\Bigl(\tfrac{\beta^2}{\epsilon^2}\ln\tfrac{n}{\epsilon}\Bigr)
\]
and size
\[
O\!\Bigl(\bigl(S_A(n + m\tfrac{\beta}{\epsilon^2},\,\tfrac{\epsilon}{294},\,\beta)\cdot\tfrac1\epsilon + n\bigr)\ln\tfrac{n}{\epsilon}\Bigr).
\]
This answers the reduction question for sparse graphs, while leaving open the removal of the linear dependence on \(m\) in the size bound [2607.07190].

The unification program also reaches the level of algorithmic templates. A single sampling-and-bunch framework can attain all state-of-the-art spanners and hopsets for general graphs by choosing the appropriate input parameters, and in the \(O(n^{1+1/k})\)-size regime it resolves an open problem by proving that every \((\alpha,\beta)\)-hopset must satisfy \(\alpha\cdot\beta=\Omega(k)\) [2108.09673].

A further development shifts attention from worst-case existence to graph-specific optimization. The optimization variant asks, for a given graph instance, for the minimum number of edges in a hopset with prescribed hopbound and stretch. In the near-exact regime, there is a randomized polynomial-time algorithm that, for directed graphs with \(\beta\ge 20\log n\), returns a \((1+\epsilon)\)-hopset of size
\[
\widetilde O\bigl(n^{3/4+\epsilon'}\,\epsilon^{-1/4}\bigr)\cdot OPT,
\]
and in undirected graphs with moderate \(\beta\), a randomized polynomial-time algorithm returning a \((1+\epsilon)\)-hopset of size
\[
\widetilde O\!\bigl(\sqrt{\beta}\,n^{1/2+1/(2\eta)+\epsilon'}\bigr)\cdot OPT
\]
for \(\eta=\lfloor \beta^{1/W_0(\ln\beta)}\rfloor\) [2502.06522].

Open questions remain explicit and technically sharp. The survey literature asks whether the classical near-exact hopbound can be reduced to a polynomial in \(\kappa\) and \(1/\epsilon\), describes the gap as “wide open,” and lists directed analogues, black-box reductions, universal constructions, and dynamic or streaming maintenance as central directions [2001.07477]. The emulator-to-hopset reduction leaves open a size bound independent of \(m\) [2607.07190]. In directed graphs, the recent lower bounds isolate aspect ratio as a decisive parameter and ask for a finer characterization of how the best achievable \(\beta\) depends on weight spread [2411.08151].

Source: https://www.emergentmind.com/topics/near-exact-hopsets