Abstract 1 Introduction 2 Preliminaries 3 Sparsity and Lightness Lower Bounds for (𝟏+𝜺)-Spanners. 4 Bounded Degree Tree Cover References

Optimal Bounds for Spanners and Tree Covers in Doubling Metrics

An La ORCID University of Massachusetts Amherst, MA, USA    Hung Le ORCID University of Massachusetts Amherst, MA, USA    Shay Solomon ORCID Tel Aviv University, Israel    Cuong Than ORCID University of Massachusetts Amherst, MA, USA    Vinayak ORCID University of Massachusetts Amherst, MA, USA    Shuang Yang University of Massachusetts Amherst, MA, USA    Tianyi Zhang ORCID State Key Laboratory for Novel Software Technology, Nanjing University, China
Abstract

It is known that any n-point set in the d-dimensional Euclidean space d, for d=O(1), admits:

  1. 1.

    A (1+ε)-spanner with maximum degree O~(εd+1) and with lightness O~(εd), for any ε>0.111The lightness is a normalized notion of weight, where we divide the spanner weight by the weight of a minimum spanning tree. Here and throughout, the O~ and Ω~ notations hide polylog(ε1) terms.

  2. 2.

    A (1+ε)-tree cover with O~(nεd+1) trees and maximum degree of O(1) in each tree.

Moreover, all the parameters in these constructions are optimal: For any 2d=O(1), there exists an n-point set in d, for which any (1+ε)-spanner has Ω~(nεd+1) edges and lightness Ω~(εd).

The upper bounds for Euclidean spanners rely heavily on the spatial property of cone partitioning in d, which does not seem to extend to the wider family of doubling metrics, i.e., metric spaces of constant doubling dimension. In doubling metrics, a simple spanner construction from two decades ago, the net-tree spanner, has O~(nεd) edges, and it could be transformed into a spanner of maximum degree O~(εd) and lightness O~(nε(d+1)) by pruning redundant edges. Moreover, a careful refinement of the net-tree spanner yields a (1+ε)-tree cover with O~(εd) trees.

Despite a large body of work, the problem of obtaining tight bounds for spanners and tree covers in the wider family of doubling metrics has remained elusive. We resolve this problem by presenting:

  1. 1.

    A surprisingly simple and tight lower bound, which shows that the net-tree spanner and its pruned version are optimal with respect to all the involved parameters.

  2. 2.

    A new construction of (1+ε)-tree covers with O~(nεd) trees, with maximum degree O(1) in each tree. This construction is optimal with respect to the number of trees and maximum degree.

Keywords and phrases:
doubling metrics, doubling spanners, Euclidean spanners, tree cover
Funding:
An La: Supported by the NSF CAREER award CCF-2237288 and an NSF grant CCF-2517033.
Hung Le: Supported by the NSF CAREER award CCF-2237288 and an NSF grant CCF-2517033.
Shay Solomon: Supported by the European Union (ERC, DynOpt, 101043159). Views and opinions expressed are however those of the author(s) only and do not necessarily reflect those of the European Union or the European Research Council. Neither the European Union nor the granting authority can be held responsible for them. Shay Solomon is also funded by a grant from the United States-Israel Binational Science Foundation (BSF), Jerusalem, Israel, and the United States National Science Foundation (NSF).
Cuong Than: Supported by the NSF CAREER award CCF-2237288 and NSF grant CCF-2121952, and a Google Ph.D. Fellowship.
Vinayak: Supported by the NSF CAREER award CCF-2237288 and an NSF grant CCF-2517033.
Shuang Yang: Supported by the NSF CAREER award CCF-2237288 and an NSF grant CCF-2517033.
Tianyi Zhang: Supported by Fundamental and Interdisciplinary Disciplines Breakthrough Plan of the Ministry of Education of China (No. JYB2025XDXM118) and the “111 Center” (No. B26023).
Copyright and License:
[Uncaptioned image] © An La, Hung Le, Shay Solomon, Cuong Than, Vinayak, Shuang Yang, and Tianyi Zhang; licensed under Creative Commons License CC-BY 4.0
2012 ACM Subject Classification:
Theory of computation Sparsification and spanners
; Theory of computation Routing and network design problems ; Mathematics of computing Graph algorithms
Related Version:
Full Version: https://arxiv.org/abs/2508.11555
Editors:
Hee-Kap Ahn, Michael Hoffmann, and Amir Nayyeri

1 Introduction

1.1 Low-Dimensional Euclidean Spaces

1.1.1 Euclidean Spanners

Let P be a set of n points in the Euclidean space d,d2, and consider the complete weighted graph GP=(P,(P2),) induced by P, where the weight of any edge (x,y)(P2) is the Euclidean distance xy between its endpoints. We say that a spanning subgraph H=(P,E,) of GP (with E(P2)) is a t-spanner for P, for a parameter t1 that is called the stretch of the spanner, if dH(x,y)txy holds x,yP. Spanners for Euclidean spaces, or Euclidean spanners, were introduced in the pioneering work of Chew [14] from 1986, which gave an O(1)-spanner with O(n) edges. The first constructions of Euclidean (1+ε)-spanners, for any parameter ε>0, were given in the seminal works of [15, 34, 35]) that introduced the Θ-graph in 2 and 3-dimensional Euclidean spaces, which was generalized for any Euclidean space d in [45, 2]. The Θ-graph is a natural variant of the Yao graph, introduced by Yao [49] in 1982, and can be described as follows.

Yao graph: pP, the space d around p is partitioned into cones of angle Θ each (O(Θd+1) cones for d=O(1)), and then edges are added between point p and its closest point in each of these cones.

The Θ-graph is defined similarly to the Yao graph: instead of connecting p to its closest point in each cone, connect it to a point whose orthogonal projection to some fixed ray contained in the cone is closest to p. Taking Θ to be cε, for small enough constant c, one can show that the stretch of the Θ and Yao graphs is at most 1+ε. Since the number of cones is asymptotically εd+1 (for d=O(1)), the maximum degree is O(εd+1), leading to a total of O(nεd+1) edges.

The tradeoff between stretch 1+ε and O(nεd+1) edges (and maximum degree O(εd+1)) is also achieved by other constructions, including the greedy spanner [2, 12, 43] and the gap-greedy spanner [46, 4]. The spatial cone partitioning of d is key to attaining the size bound O(nεd+1) in these constructions, either in the constructions themselves or in their analysis. In 2019, Le and Solomon [38] showed that this stretch-size tradeoff is existentially tight: For any constant d2, there exists an n-point set in d (basically a set of evenly spaced points on the d-dimensional sphere), for which any (1+ε)-spanner has Ω(nεd+1) edges.

Light spanners.

Another basic property of spanners is lightness, defined as the ratio of a spanner’s total weight to the weight of the Minimum Spanning Tree of P. A long line of work [2, 16, 18, 44, 43, 6, 38], starting from the paper of Das et al. [16] in 1993, showed that for any point set in d, the greedy (1+ε)-spanner of [2] has constant (depending on ε and d) lightness. The exact dependencies on ε and d in the lightness bound were not explicated in [2, 16, 18]. In their seminal work on approximating TSP in d using light spanners, Rao and Smith [44] showed that the greedy spanner has lightness εO(d), and they raised the question of determining the exact constant hiding in the exponent of their upper bound. The proofs in [2, 16, 18, 44] were incomplete; the first complete proof was given in [43], where a lightness bound of O(ε2d) was established. This line of work culminated with the work of Le and Solomon [38], which improved the lightness bound to O~(εd), where we shall use the O~ and Ω~ notations to suppress polylog(ε1) terms. The exact lightness bound here is O(εdlog(1/ε)), but for the sake of brevity we will mostly disregard polylog(ε1) terms from now on. They [38] also showed that this stretch-lightness tradeoff is existentially tight (up to a log(ε1) factor): for any constant d2, there exists an n-point set in d (the same set of evenly spaced points on the d-dimensional sphere), for which any (1+ε)-spanner has lightness Ω(εd).

1.1.2 Euclidean Tree Covers

A more structured variant of spanners is the tree cover, defined as a collection 𝒯 of trees for a given graph or metric space such that, for every pair of vertices, there exists a tree in 𝒯 that preserves their distance up to a given stretch factor, without shortening distances. Due to their strong structural properties, tree covers play a pivotal role in routing, (path-reporting) distance oracles, and various other algorithmic applications. The first tree cover construction for a graph or metric family was given by Arya et al. [3] – known as the “Dumbbell Theorem” – asserting that any low-dimensional Euclidean space admits a (1+ε)-stretch tree cover of size O(εdlog(1/ε)), for any ε>0. The lower bound by Le and Solomon [37] on the size of spanners directly implies that any (1+ε)-stretch tree cover must have size at least Ω(εd+1). Recently, Chang et al. [13] closed this longstanding gap between the upper and lower bounds on the size of tree covers, up to a log(1/ε) factor, by improving the size upper bound to O(εd+1log(1/ε)). Another important quality measure of tree covers is the maximum degree of a vertex in any tree; in particular, a constant degree tree cover enables compact routing schemes with small routing tables [13]. Chang et al. constructed a constant degree tree cover without increasing the size beyond O(εd+1log(1/ε)). Whether the log(1/ε) factor separating this upper bound from the lower bound is necessary remains an open problem.

1.2 Doubling Metrics

1.2.1 Doubling Spanners

Euclidean spanners have been extensively studied over the years [34, 35, 2, 17, 3, 18, 4, 44, 26, 22, 1, 7, 24, 9, 19, 48, 20, 38], with a plethora of applications, such as in geometric approximation algorithms [44, 27, 30, 28], geometric distance oracles [27, 30, 29, 28], network design [32, 42] and machine learning [25]. (See the book [43] for an excellent account on Euclidean spanners and their applications.) There is a growing body of work on doubling spanners, i.e., spanners for the wider family of doubling metrics;222The doubling dimension of a metric space is the smallest value d such that every ball B in the metric can be covered by at most 2d balls of half the radius of B; a metric space is called doubling if its doubling dimension is constant. The doubling dimension generalizes the standard Euclidean dimension, as the doubling dimension of the Euclidean space d is Θ(d). see [22, 8, 7, 31, 23, 24, 47, 10, 20, 11, 48, 6, 41, 33, 40], and the references therein. A common theme in this line of work is to devise constructions of spanners for doubling metrics that are just as good as the analog Euclidean spanner constructions. Alas, this may not always be possible, as doubling metrics do not possess the spatial properties of Euclidean spaces – and in particular the spatial property of cone partitioning, which is key to achieving the aforementioned stretch-size and stretch-lightness upper bounds. Despite this shortcoming, the basic packing bound in doubling metrics can be used to construct, via a simple greedy procedure, a hierarchy of nets, which induces the so-called net-tree [31, 22, 8]. (See Section 2 for the packing bound and other definitions.) Equipped with such a hierarchy of nets, a (1+ε)-spanner with O~(nεd) edges is constructed as follows [31, 22, 8].

Net-tree spanner [31, 22, 8]: Let X=N0N1N2NlogΔ be a hierarchy of nets of a doubling metric (X,dX) with minimum distance 1 and maximum distance Δ, where Ni is a 2i-net for Ni1, 1ilogΔ.

Let E1= and for each 0ilogΔ:

Ei={(x,y)|x,yNi,dX(x,y)(4+32ε)2i}Ei1 (1)

Then (X,0ilogΔEi,dX) is a (1+ε)-spanner of (X,dX).

The net-tree spanner is a simple and basic spanner construction from over 20 years ago, and it provides the state-of-the-art size bound, O(nεdlog(1/ε)). Le and Solomon [39] presented a unified framework for transforming sparse spanners into light spanners by carefully pruning redundant edges. In particular, for doubling metrics, the framework of [39] provides a pruned spanner, with a lightness bound of O~(ε(d+1)), which is the state-of-the-art lightness bound for doubling metrics. Interestingly, these upper bounds on the size and lightness in doubling metrics exceed the respective Euclidean bounds by a factor of ε1. Despite a large body of work in the area, it has remained a longstanding open question whether the superior Euclidean bounds of O~(nεd+1) edges and lightness O~(εd) could be achieved also in doubling metrics.

Question 1.

Can one get a construction of (1+ε)-spanners in doubling metrics with O~(nεd+1) edges and/or lightness O~(εd)?

This question is related to a possibly deeper question, regarding the (im)possibility to generalize the spatial cone partitioning in d to arbitrary doubling metrics.

1.2.2 Doubling Tree Covers

Bartal et al. [5] constructed a tree cover of size O(εdlog(1/ε)) using the net-tree by redistributing cross edges at each level into color classes. However, in their construction, a vertex in one of the trees may have many children, resulting in unbounded degree, potentially as large as log(Δ), where Δ is the aspect ratio. Motivated by the Euclidean setting, it is natural to ask whether one can obtain a tree cover with bounded degree in doubling metrics.

Question 2.

Can one construct a (1+ε)-tree cover in doubling metrics whose maximum degree is an absolute constant, i.e., O(1), while having size O~(εd)?

1.3 Our contribution

We give a negative answer to Question 1 by a surprisingly simple and tight lower bound.

Theorem 1.

For any integer constant d1, parameter ε(0,1), and n+ satisfying εd=O(n), there exists an n-point ultrametric space (X,dX) of doubling dimension d such that any (1+ε)-spanner G=(X,E,w) for X has Ω(nεd) edges and lightness Ω~(ε(d+1)).

Since the net-tree spanner has sparsity 333The sparsity of a spanner H is the ratio between number of edges in H over n1, for n>1. O(εd) [22], our lower bound shows, in particular, that the net-tree spanner is optimal; this is our key conceptual contribution.

Second, as our main technical contribution, we show that a tree cover with constant degree and nearly optimal size is achievable, resolving Question 2 affirmatively.

Theorem 2.

For any metric (X,dX) with doubling dimension d, parameter 0<ε<1, there exists a (1+ε)-stretch tree cover 𝒯 of size O~(εd), where each tree in 𝒯 has maximum vertex degree O(1).

To complement our result, Theorem 3 shows that it is impossible to eliminate the 𝗉𝗈𝗅𝗒𝗅𝗈𝗀(1/ε) factor from the size of any constant-degree tree cover (the detailed proof is provided in Appendix A of the full version of this paper).

Theorem 3.

For any integer constant d1, any parameter ε(0,1), and any n+ satisfying εd=O(n), there exists an n-point ultrametric space (X,dX) of doubling dimension d such that any (1+ε)-spanner G=(X,E,w) for X has maximum degree Ω(log(1/ε)εd).

Theorem 3 demonstrates an inherent multiplicative gap of log(1/ε) between the maximum degree and the optimal sparsity of (1+ε)-spanners in doubling metrics, a gap that does not appear in Euclidean spaces. Although Theorem 3 is not the main focus of our work, it provides an interesting separation between Euclidean and doubling metrics. Moreover, it shows that the 𝗉𝗈𝗅𝗒𝗅𝗈𝗀(1/ε) term in Theorem 2 (hidden via the O~-notation) is unavoidable in doubling metrics, whereas its necessity in Euclidean spaces remains an open question.

1.4 Technical Overview

Our tree cover construction provides the key technical contribution of this work. In what follows, we provide a technical overview of this construction.

Our tree cover construction relies on a hierarchical structure – specifically, the standard net-tree – and builds a tree cover so that every cross edge e (i.e., every edge in Ei) is contained in some tree. For every pair of points (u,v), there exists a cross edge (u,v) such that the path from u to v that goes through (u,v) provides a good approximation of dX(u,v). The distance between u and v is then preserved by the tree T that contains (u,v); in this case, we say that (u,v) is covered by T.

The construction proceeds in two phases. In Phase 1, we partition the cross edges into congruent classes according to their levels modulo log(1/ε), and handle each class independently. For each class, we distribute the cross edges into each tree in our tree cover. Note that after Phase 1, we already obtain a tree cover of stretch (1+ε) and size O~(εd). In Phase 2, for a fixed level, we replace edges in a controlled manner to obtain the required bounded degree property, without significantly increasing the stretch. The desired (1+ε) stretch bound can then be recovered by a simple scaling argument.

Our Phase 2 is inspired by the degree-reduction procedure of [13]. To achieve constant degree, we proceed in two steps. First, using the technique in [8], we reroute the edges in each tree T in the cover produced by Phase 1, obtaining trees of maximum degree εO(d), while the stretch increases only to 1+O(ε). Second, we use the rerouting technique in [13] to obtain absolute constant maximum degree. For each vertex u in each such tree T, we replace the edges connecting u to its children by a binary tree. Suppose that l is the maximum distance from (the representative of) u to its children, this step incurs an additive distortion of Θ(log(εd)l) for every pair of points covered by T.

A key technical challenge in integrating Phase 2 into Phase 1 is that the edges replaced by binary trees must be sufficiently short relative to the cross edges at the same level. Otherwise, the additive distortion for a pair (u,v) might become comparable to d(u,v) itself, producing a stretch larger than 2. In the construction of [5], a point u in a tree T may be incident to εd cross edges of nearly equal length; applying Phase 2 on top of [5] would therefore not yield a tree cover with good stretch.

Our first idea is to ensure that, in Phase 1, for every level we partition the cross edges into O(εd) sets such that in each tree of the resulting cover, every point u is incident to at most one edge from each set. This motivates partitioning Ei into sets of matchings.

We implement Phase 1 as follows. We initialize tree cover 𝒯 with O~(εd) empty forests. For each congruent class modulo log(1/ε), we process the levels in increasing order. For each level i, we partition Ei into sets of matchings i. For each matching Mi, we assign M to forest F𝒯, and add all edges in M to F. We continue this process until the top level.

A second challenge is ensuring that the graph F remains a forest after inserting the edges of M. Specifically, we guarantee that for any two edges (u,v) and (u,v) in M, there is no tree in the forest F containing both u and u. This leads to our second idea: each matching M must be pairwise far. That is, for every Mi and every two distinct edges e,eM, the distance between (any two endpoints of) e and e must be at least Θ(2i). Under this condition, an inductive argument shows that every tree in each forest F has small diameter (<2i); thus, after connecting each tree to an edge in M, a sufficient separation still remains, ensuring that the resulting graph is indeed a forest. Such matchings, also known as red–blue matchings, were also used by [21, 36] in different contexts.

2 Preliminaries

Let (X,dX) be a metric space. The aspect ratio Δ of X is the ratio of its largest to smallest pairwise distances. The following definitions are crucial for our construction and proofs.

Definition 4 (r-net).

For a metric space (X,dX), an r-net NX is a set satisfying:

  1. 1.

    [Packing] for every distinct u,vN, dX(u,v)>r

  2. 2.

    [Covering] for every xX, there exists uN such that dX(u,x)r.

Definition 5 ((α,β)-Net Tree).

An (α,β)-net tree τ is induced by a hierarchy of nets N1=X,N0,N1,,NlogΔ at levels 1,0,1,,logΔ of τ resp., satisfying the following:

  1. 1.

    (Hierarchy) NiNi1 for every i[logΔ].

  2. 2.

    (Net) For every i{0,1,,logΔ}, Ni is an (αβi)-net of Ni1.

The following lemma gives the basic packing bound of doubling metrics.

Lemma 6 (Packing Lemma).

Metric (X,dX) has doubling dimension d. If SX is a subset with pairwise distance at least r that is contained in a ball of radius R, then |S|(4Rr)d.

Definition 7 (Ultrametric).

An ultrametric (X,dX) is a metric space that satisfies the strong triangle inequality, i.e., for any u,v,wX, we have dX(u,w)max{dX(u,v),dX(v,w)}.

Definition 8 (k-HST).

Let k>1 be a fixed constant. A metric space (X,dX) is a k-hierarchical well-separated tree (k-HST) if there exists a bijection ϕ from X to leaves of a rooted tree T in which:

  1. 1.

    Each node vT is associated with a label Γ(v) such that Γ(v)=0 if v is a leaf, and Γ(v)kΓ(u) (Γ(v)>0) if v is an internal node and u is a child of v.

  2. 2.

    dX(x,y)=Γ(LCA(ϕ(x),ϕ(y))), where LCA(u,v) is the lowest common ancestor of u,v.

Lastly, we defer some proofs of lemmas/claims (marked by ()) to the appendix of the full version.

3 Sparsity and Lightness Lower Bounds for (𝟏+𝜺)-Spanners.

This section is devoted to the proof of Theorem 1. We state it again below: See 1

We begin by constructing an ultrametric space (X,dX) for which we later prove that any (1+ε)-spanner must satisfy the size and lightness lower bound.

Without loss of generality, assume that n is a power of 2d. More specifically, if n is not a power of 2d, let n~=(2d) be the largest integer power of 2d before n, with n~n/2d, and apply the argument to n~ instead of n. The size lower bound will decrease by at most an O(1) factor 2d and the lightness lower bound will remain the same.

We define the following 2-HST (X,dX). Let T be a rooted perfect 2d-ary tree with n leaf nodes: all leaf nodes are at the same depth log2dn=log2nd and each internal node has 2d children. The height of a node is its distance to the closest leaf, so the leaf nodes have height 0, their parents have height 1, and so on, until reaching the root, at height log2nd. For each node v in the tree we assign a label, denoted by Γ(v). Each leaf node is assigned label 0, and each node of height i=1,,log2nd is assigned label 2i. For any x,yT, we denote the lowest common ancestor of x and y by LCA(x,y). Let ψ be any arbitrary bijection mapping X to the leaf nodes of T. For any u,vX, we define dX(u,v) to be Γ(LCA(ψ(u),ψ(v))).

We start with the following basic observation in Claim 9. Although to the best of our knowledge, this precise statement does not appear in previous work, weaker (though similar) versions of this statement can be found in [7, 21].

Claim 9.

(X,dX) is an ultrametric space of doubling dimension d.

Proof.

First, we argue that (X,dX) is an ultrametric space. By the definition of dX, we have dX(u,v)=0 iff u=v. Further, dX(u,v)=dX(v,u) holds for any u,vX since LCA(ψ(u),ψ(v))=LCA(ψ(v),ψ(u)). It remains to show dX(u,v)max{dX(u,w),dX(w,v)}, for any u,v,wX. Let Tu (Tv) denote the subtree rooted at the child node of LCA(ψ(v),ψ(u)) that contains ψ(u) (ψ(v)). Note that for any wX, ψ(w) cannot belong to both Tu and Tv, so either LCA(ψ(v),ψ(w)) is not contained in Tv, or LCA(ψ(u),ψ(w)) is not contained in Tu. Thus, dX(u,v)max{dX(u,w),dX(w,v)} holds.

Second, we prove that (X,dX) has doubling dimension d.

We first show that the doubling dimension is no smaller than d. Let B(u,r) be a ball centered at uX with radius r+. Let ψ(fu) be the parent of ψ(u) in T. Consider r=2. By the construction of T, every point vX corresponding to a child of ψ(fu) satisfies dX(u,v)=Γ(ψ(u),ψ(v)). Hence, B(u,2) contains all points corresponding to the children of ψ(fu), and contains at least 2d distinct points. On the other hand, for any vX, the ball B(v,1) contains only the single point v. Hence, covering B(u,2) requires at least 2d balls of radius 1. This implies that the doubling dimension is no smaller than d.

Then, we show that the doubling dimension of (X,dX) is no bigger than d. When r<2, B(u,r) contains only one point u and is covered by B(u,r2). When 2ir<2i+1 for some integer i1, let ω be the ancestor of ψ(u) of height i in T and let Tω be the subtree rooted at ω. By the construction of T, all pairwise distances between the points corresponding to the leaf nodes of Tω are no larger than the label 2i of ω. On the other hand, the distance between u and any point corresponding to a leaf node outside Tω is no smaller than the label 2i+1 of the parent node of ω, as their lowest common ancestor is outside Tω. Thus, B(u,r) consists of the (2d)i=2di points that correspond to the leaf nodes of Tω. Let Lj be the set of points corresponding to the leaves of the subtree rooted at the jth child of ω, for each j[2d], and let vj be an arbitrary point in Lj. Since all children of ω have label 2i1, all pairwise distances between the points in Lj are no larger than 2i1. So B(vj,r2)B(vj,2i1) covers all (2d)i1=2d(i1) points in Lj, for each j[2d], implying that the ball B(u,r) can be covered by the 2d balls B(vj,r2),j[2d].

Then, we demonstrate the following property of (X,dX) in Claim 10 that will be crucial for proving the size and lightness lower bounds.

Claim 10.

Let G=(X,E,w) be an arbitrary (1+ε)-spanner for (X,dX). For any point pair u,vX such that dX(u,v)<2ε, edge (u,v) must be in G.

Proof.

Consider any pair u,vX. We first argue that if edge (u,v) is not in G, then dG(u,v)dX(u,v)+2. In this case, consider a shortest path from u to v in G, and note that it traverses an intermediate point wX. By the triangle inequality, dG(u,v)dG(u,w)+dG(w,v)dX(u,w)+dX(w,v). Since the minimum pairwise distance is 2 and (X,dX) is an ultrametric, dG(u,v)dX(u,w)+dX(w,v)=max{dX(u,w),dX(w,v)}+min{dX(u,w),dX(w,v)}dX(u,v)+2.

We conclude that if dX(u,v)<2ε, then edge (u,v) must be in G, otherwise by the above assertion we would get dG(u,v)dX(u,v)+2>(1+ε)dX(u,v), which is a contradiction to G being a (1+ε)-spanner for (X,dX).

In the following subsections, we apply Claim 10 to prove the size and lightness lower bounds of any (1+ε)-spanner G=(X,E,w) for (X,dX).

Size lower bound.

Let Dmax be the largest label among nodes in T that is strictly smaller than 2ε. Since the labels in T grow geometrically by a factor of 2, there exists a label in [1ε,2ε), hence 1εDmax<2ε. Let imax=log(Dmax)log(1ε) be the height of nodes with label Dmax. Each subtree of height imax has (2d)imaxεd leaf nodes (this is where the restriction εd=O(n) comes into play) and the pairwise distances between the corresponding points are bounded by Dmax<2ε, hence Claim 10 implies that the spanner G connects all these points by a clique. Thus each point in X has a degree of at least εd1 in G, and so G contains at least n(εd1)/2=Ω(nεd) edges, which proves the size lower bound.

Lightness lower bound.

We begin by giving the lower bound on the weight of G. Consider an arbitrary point uX and denote the subtree rooted at the ancestor of ψ(u) at level imax as Timax(u). We have shown that u is connected to all (2d)imaxεd points corresponding to the leaf nodes in Timax(u). Note also that only (2d)imax1 points among those may belong to the child subtree of Timax(u) that contains ψ(u), while the remaining (2d)imax(2d)imax1εd/2 points belong to other child subtrees of Timax(u), and are therefore at distance Dmax=2imax1/ε from u. Thus the total weight of edges incident on any point uX in G is at least εd/2(1/ε)=ε(d+1)/2, and so the weight of G is at least nε(d+1)/4=Ω(nε(d+1)).

Claim 11 concludes the proof of the lightness lower bound for d2. As for d=1, Claim 11 is not enough, since we need w(MSTX) to be O(nlog(ε1)) rather than O(nlogn) to obtain the required lightness bound; we will handle this issue after the proof of Claim 11.

Claim 11.

w(MSTX)=O(n) for any d2 and w(MSTX)=O(nlogn) for d=1.

Proof.

Let PX be the Hamiltonian path of X that traverses the points of X according to the induced left-to-right ordering in T of the corresponding leaves. To prove the claim, we will show that w(PX)=O(n) for any d2 and w(PX)=O(nlogn) for d=1. (It can be shown that PX=MSTX, but there is no need for that, since we only need to upper bound w(MSTX) and we have w(PX)w(MSTX).)

Denote the root of T by r, and note that Γ(r)=2log2nd=n1/d. Observe that the number of edges in P of weight Γ(r)=n1/d is exactly one less than the number of children, i.e., 2d1. In general, note that the number of edges in P of weight Γ(r)2i=n1/d2i is 2di(2d1), for each i[0,log2nd1]. It follows that

w(PX)=i=0log2nd12di(2d1)n1/d2i=(2d1)n1/di=0log2nd1(2d1)i.

If d=1, we get that w(PX)=nlog2n.
If d2, then i=0log2nd1(2d1)i is a geometric sum with rate 2d12, hence we have

i=0log2nd1(2d1)i 22(d1)(log2nd1)= 22log2ndlog2nd+1= 4n2dn1/d,

and we get that w(PX)=(2d1)n1/di=0log2nd1(2d1)i4n. Finally, consider the case d=1. Claim 11 gives w(MSTX)=O(nlogn), but to get the required lightness bound we need an upper bound of w(MSTX)=O(nlog(ε1)).

Without loss of generality, we may assume that n=n and that n=Θ(ε1). In this case, Claim 11 implies that w(MSTX)=O(nlogn)=O(nlog(ε1)), and the lightness of G is Ω(ε2log(ε1)), as required. We introduce two different reductions justifying this assumption.

First, one can simply take the same ultrametric as before on top of n points (instead of n) as leaves, and then add nn points that are arbitrarily close to one of the n points, to get a metric space with n points and basically the same MST weight; while this reduction provides the required lightness lower bound, it does not preserve the size lower bound.

To get a single instance for which both lower bounds apply, one can do the following. Create n/n vertex-disjoint copies of the same ultrametric (X,dX) as before, each on top of n points as leaves, denoted by (X(j),dX(j)), for each j[n/n], and place these copies “on a line”, so that neighboring copies (X(j),dX(j)) and (X(j+1),dX(j+1)) are at distance say 2n from each other. Since any two copies are sufficiently spaced apart from each other (w.r.t. the maximum pairwise distance in each copy, namely n), any (1+ε)-spanner edge between two points in the same copy must be contained in that copy, hence we can apply the aforementioned size and weight spanner lower bounds on each of the copies separately, and then aggregate their size and weight bounds to achieve the same size and weight lower bounds as before, namely Ω(nε1) and Ω(nε2), respectively, for the entire spanner. Since each copy has weight O(nlogn) by Claim 11 and as neighboring copies are sufficiently close to each other, the MST weight of this metric is bounded by O(nlogn)(n/n)=O(nlog(ε1)), which proves the required lightness lower bound Ω(ε2log(ε1)).

4 Bounded Degree Tree Cover

Theorem 2. [Restated, see original statement.]

For any metric (X,dX) with doubling dimension d, parameter 0<ε<1, there exists a (1+ε)-stretch tree cover 𝒯 of size O~(εd), where each tree in 𝒯 has maximum vertex degree O(1).

4.1 Phase 1: Unbounded Degree Tree Cover

Lemma 12.

For any metric (X,dX) with doubling dimension d, and any parameter 0<ε<1, there is a (1+ε)-stretch tree cover 𝒯 of size O(εdlog(1/ε)), where the maximum degree of a vertex in each of the trees in 𝒯 is O(εdlog1/εΔ).

We first give a construction of a tree cover with maximum degree O(εdlog1/εΔ). Later, in Section 4.2, we modify this tree cover to get constant degree bound.

  1. 1.

    Creating intra-level edge matchings. Let δp=2p for p{0,,log(1/ε)}, and let τδp be a (δp,ε1)-net tree (see Definition 5) of X and Ni be the δpεi-net in τδp. For each p, construct graph Gp, where V(Gp)={uiuNi}i and E(Gp)=. Define a set of cross edges Ec={(ui,vi)dX(u,v)(δp3εi+1,4δp3εi+1)}i. Partition Ec into γ=O(εd) matchings ={M1,,Mγ}, with the property that in each M, dX(ui,vi)>14δpεi, for all matched ui, vi. Initialize Tkp=(V(Gp),Mk) for k{1,,γ}.

  2. 2.

    Creating inter-level edges. For each T=Tkp, construct edges between levels i and (i1) in the order i=0,1,. Create a relabel function ΦiT:NiNi, initialized as ΦiT(ui)=ui for all uiNi. At each level i, process vertices in the order specified below:

    1. (a)

      Matched vertices. For each matched ui in Mk, add (ui,ui1). For vBX(u,4δpεi) such that vi1 exists, add (ui,vi1) if either vi1 is unmatched, or (vi1,xi1)Mk and xi1 does not have a parent at level i.

    2. (b)

      Unmatched vertices. For each unmatched ui: enumerate its child nodes in net-tree τδp, denoted by Cui={ai10,ai11,,ai1η}. For every ai1Cui:

      • Add (ΦiT(ui),ai1) in the following cases: 1) ai1 is unmatched and does not have a parent, and 2) (ai1,bi1)Mk and both ai1, bi1 do not have parents at level i.

      • Otherwise, if ai1 has a parent xi, or (ai1,bi1)Mk and bi1 has a parent xi, merge xi with ΦiT(ui). Set ΦiT(xi)=ΦiT(ui) if ui+1 exists, otherwise set ΦiT(ui)=xi. (Note that ai1 and bi1 will not have parents at the same time.)

  3. 3.

    Contraction. Contract all vertices sharing label u. 𝒯={Tkp}p,k forms a tree cover.

Throughout the construction (even after relabeling), the weight of any edge, wT(ui,vj)=dX(u,v). In Step 2, we slightly abuse the notation, setting ΦiT(ui)=xi, instead of ΦiT(ΦiT(ui))=xi. This is fine, since a node undergoes relabeling at most once. This finishes the construction of the tree cover.

Unless stated otherwise, we refer to the intermediate trees just prior to Step 3 as the tree cover 𝒯 and prove our bounds on them. This works since the contraction of only 0-weight edges in Step 3 does not affect the stretch.

We first prove that each T in 𝒯, after contraction, is a tree (see Corollary 16). Then, we bound the number of trees in 𝒯 (see Lemma 17). Finally, we use Claim 18 and Claim 19 to bound the final stretch in Lemma 20 (proof deferred to the Appendix B of the full version).

We first require the following combinatorial lemma of [21], used in Step 1 of our construction.

Lemma 13 (Lemma 6 of [21]).

Consider a graph G=(V,EbEr) with disjoint sets of blue edges Eb and red edges Er, such that the maximum blue degree is δb1, and the maximum red degree is δr1. Then there is a set of at most γ=O(δbδr) matchings ={M1,,Mγ} of G such that:

  1. 1.

    Ebi=1γMi, and

  2. 2.

    for every matching M, there is no red edge whose both endpoints are matched by M.

 Remark 14.

Note that we apply this lemma in Step 1 of the construction. More specifically, Eb=Ec={(ui,vi)dX(u,v)(δp3εi+1,4δp3εi+1)}i, and Er={(ui,vi)dX(u,v)14δpεi}. In each M, dX(ui,vi)>14δpεi holds for all matched ui, vi.

We begin by proving that each Tkp𝒯 is a tree. We rely on the following inductive invariant: after processing the ith level, every level-i subtree has a diameter at most O(ε) factor of the newly added level-i cross edge. This large separation ensures that the forest remains a forest.

Lemma 15.

Each Tkp𝒯 is a tree.

Proof.

Let T=Tkp. We will prove the following two properties. First, each vertex at level-i is connected to at most one vertex at level-(i+1). Second, for every level-i cross edge, at most one of the end points is connected to a level-(i+1) vertex. These two arguments are sufficient condition for T to be acyclic since all cross edges are mutually vertex-disjoint.

We begin by proving that these two properties hold in Step 2a. For any vertex vi1BX(u,4δpεi) and ui is a matched vertex of T, vi cannot also be a matched vertex of T since dX(u,v)>14δpεi should hold by Lemma 13. Hence, any vertex vi1 is connected to at most one vertex at level-(i+1). For any matched cross edge (ui1,vi1)T, we have dX(u,v)(δp3εi,4δp3εi), in which case ui and vi cannot both be matched vertices in level-i since this requires dX(u,v)>14δpεi. Thus, (ui,ui1) and (vi,vi1) will not be added at the same time, and (ui1,vi1) cannot have parents at both endpoints. Next, by construction (Step 2b), the two properties hold. This is because an edge between some ai+1 and bi is created only when either bi is unmatched and does not have a parent, or some (bi,ci) exists and neither of bi or ci have a parent. At each level of the net tree τδp, we go over all vertices.

Let rt be the vertex at the highest level of T. For any two vertices u,vV(T), we know that there is a path from u to rt, and from v to 𝑟𝑡 in T. This holds because each vertex is either connected to a vertex one level above, or to another vertex at the same level via a cross edge (since each such vertex is enumerated in Step 2b). Hence, there is a walk between u and v via rt in T. This implies that T is connected.

Corollary 16.

() After the contraction procedure, each Tkp in 𝒯 remains a tree.

Lemma 17.

The number of trees in 𝒯 is O(εdlog(1/ε)).

Proof.

For each Gp (p{0,1,log(1/ε)}, we take V=V(Gp), Eb=Ec. We define a set of dummy edges in Gp: Er={(ui,vi)dX(u,v)14δpεi}i. By packing lemma, the maximum degree in Eb at each level i is O((δpε(i+1)/δpεi)d)=O(εd).

Similarly, the maximum degree in Er is O(14δpεi/δpεi)d=cd for some constant c. We can now invoke Lemma 13 on G=(V,EbEr) with δb=O(εd) and δr=cd to get O((cε1)d)=O(εd) matchings (each corresponding to a tree) with the required properties. Hence, the total number of trees in 𝒯 is O(γlog(1/ε))=O(εdlog(1/ε)).

Claim 18.

() Let T=Tkp. Let (ai,bi) be a cross edge and (ci,di1) be an inter-level edge in T. For ε<1/156, we have

δp4εi+1dT(ai,bi)3δp2εi+1,anddT(ci,di1)272δpεi.
Claim 19.

() For any Tkp𝒯, the distance from a leaf node u0 to any of its ancestors at level i is at most 2δpε(i+1).

Lemma 20.

() For every u,vX, there exists a tree T in the tree cover 𝒯 such that: dT(u,v)(1+ε)dX(u,v).

4.2 Phase 2: Bounding the Degree

The algorithm in Section 4.1 builds the tree cover 𝒯 of size O(εdlog(1/ε)). In this section, we focus on adjusting edges of these trees to achieve constant degree bound. We state the main algorithm and refer to the Appendix B of the full version for the detailed analysis of this section.

Let T denote a tree in the tree cover 𝒯 constructed before Step 3 (contraction). Let l(x) be the highest level of x in the net tree τδp. We assign a direction to each edge of T based on the levels of its endpoints in τδp. An edge e=(aibj) is going out from ai to bj if l(a)l(b). If l(a)=l(b), we assign the direction arbitrarily. Similar to Phase 1, we categorize the set of directed edges into two categories:

  1. 1.

    An edge (xipi+1) is an inter-level edge at level i if xNiNi+1 and pNi+1. Let Di be the set of inter-level edges at level i. We denote the set of points pi+1 incident to an inter-level edge (xipi+1) at level i by Di(x), and the set of points wi incident to an inter-level edge (wixi+1) at level i by Di+(x).

  2. 2.

    An edge (xiyi) is a cross edge at level i if x,yNi. Let Ci be the set of all cross edges at level i. Given a point xi, we denote the set of all points yi having a cross edge (yixi) by Ci+(x), and the set of all points zi having a cross edge (xizi) by Ci(x).

By our Phase 1 construction, we observe the following bounds:

Observation 21.

Given a fixed level i, for every vertex xNi, |Di+(x)|=O(εd), |Di(x)|1, |Ci+(x)|1, |Ci(x)|1. Furthermore, for every point x, we have i|Di(x)|1, and i|Ci(x)|2.

We use the degree reduction technique of [8] to reduce degree associated with both cross edges and inter-level edges. The input and output of this technique are specified in Remark 22 and Lemma 23, respectively.

 Remark 22.

The input tree T must be a tree in the tree cover, and we consider a subset E of edges in T. Let Ei denote the edges in E at level i, and Ei+(x) (resp., Ei(x)) denote the set of points connected to x by incoming (resp., outgoing) edges in Ei. The subset E must satisfy the following degree constraints for every point x: 1) O(1) in-coming edges at every level, 2) at most one out-going edge at every level, and 3) totally O(1) out-going edges over all levels. Formally, |Ei+(x)|din, |Ei(x)|1 and i|Ei(x)|dout, where din>0 and dout>0 are constants.

Lemma 23.

Given a tree cover 𝒯, let T be a tree in 𝒯 that satisfies Remark 22. There exists an algorithm redirecting edges in E of T to obtain a tree T with a new edge set E such that:

  • (Bounded degree.) Every point has at most 2+2din+dout edges in E.

  • (Small stretch.) For every edge (y,x) in E of T, there exist a path π in T from y to x such that π contains edges in E, and dT(y,x)(1+O(ε))dX(y,x).

Algorithm.

We use two operations: redirect cross edges and redirect inter-level edges.

  • Step 1. Redirecting cross edges. We apply Lemma 23 to re-direct cross edges, where E=iCi. Ei and Ei+ in Remark 22 will respectively be Ci and Ci+, where |Ci+(x)|1,|Ci(x)|1,t|Ct(x)|2 for every point x and every level i (by Observation 21).

  • Step 2. Redirecting inter-level edges.

    • Step 2.1. We visit levels from the bottom to top. Let Qi be the set of new edges at level i created by this step. Initialize Qi=. Let Qi(p) be the set of points q having an edge (pq) in Qi. Let Qi+(p) be the set of points q having an edge (qp) in Qi.

      Given a fixed level i, for every vertex xi+1Ni+1, recall that the set of in-coming inter-level edges at level i of x is Di+(x)={wi(1),wi(2),}, we know |Di+(x)|O(εd). We sort Di+ by increasing order of distances to xi+1. Consider wi(j)Di+(x):

      • *

        If j2: keep (wi(j)x).

      • *

        If j>2: let j=(j1)/2; remove (wi(j)xi+1) and create edge (wi(j)wi(j)) with weight dX(wi(j),wi(j)). We add (wi(j)wi(j)) to Qi.

    • Step 2.2. We apply Lemma 23 to re-direct inter-level edges remaining after Step 2.1, where E=iDi. Ei and Ei+ in Remark 22 will be Di and Di+ respectively, where |Dk+(p)|2,|Dk(p)|1,t|Dt(p)|1 for every point p and level k. These bounds hold by Observation 21 and by Step 2.1 (when j2).

    • Step 2.3. We apply Lemma 23 to re-direct E=iQi created by Step 2.1.

      Ei and Ei+ in Remark 22 will be respectively Qi and Qi+.

This finishes the construction of the bounded degree tree cover.

References

  • [1] P. K. Agarwal, Y. Wang, and P. Yin. Lower bound for sparse Euclidean spanners. In Proc. of 16th SODA, pages 670–671, 2005.
  • [2] I. Althöfer, G. Das, D. Dobkin, D. Joseph, and J. Soares. On sparse spanners of weighted graphs. Discrete Computational Geometry, 9(1):81–100, 1993. doi:10.1007/BF02189308.
  • [3] S. Arya, G. Das, D. M. Mount, J. S. Salowe, and M. Smid. Euclidean spanners: Short, thin, and lanky. In Proceedings of the Twenty-seventh Annual ACM Symposium on Theory of Computing, STOC ’95, pages 489–498, 1995.
  • [4] S. Arya and M. H. M. Smid. Efficient construction of a bounded degree spanner with low weight. Algorithmica, 17(1):33–54, 1997. doi:10.1007/BF02523237.
  • [5] Yair Bartal, Arnold Filtser, and Ofer Neiman. On notions of distortion and an almost minimum spanning tree with constant average distortion. J. Comput. Syst. Sci., 105:116–129, 2019. preliminary version published in SODA 2016. doi:10.1016/j.jcss.2019.04.006.
  • [6] G. Borradaile, H. Le, and C. Wulff-Nilsen. Greedy spanners are optimal in doubling metrics. In Proceedings of the 30th Annual ACM-SIAM Symposium on Discrete Algorithms, SODA ‘19, pages 2371–2379, 2019. doi:10.1137/1.9781611975482.145.
  • [7] H. T.-H. Chan and A. Gupta. Small hop-diameter sparse spanners for doubling metrics. In Proc. of 17th SODA, pages 70–78, 2006.
  • [8] T-H. Hubert Chan, Anupam Gupta, Bruce M. Maggs, and Shuheng Zhou. On hierarchical routing in doubling metrics. In Proceedings of the 16th Annual ACM-SIAM Symposium on Discrete Algorithms, SODA’05, 2005.
  • [9] T.-H. Hubert Chan, Anupam Gupta, Bruce M. Maggs, and Shuheng Zhou. On hierarchical routing in doubling metrics. ACM Trans. Algorithms, 12(4):55:1–55:22, 2016. Preliminary version appeared in SODA 2005. doi:10.1145/2915183.
  • [10] T.-H. Hubert Chan, Mingfei Li, and Li Ning. Sparse fault-tolerant spanners for doubling metrics with bounded hop-diameter or degree. In ICALP (1), pages 182–193, 2012. Full version in Algorithmica’15. doi:10.1007/978-3-642-31594-7_16.
  • [11] T.-H. Hubert Chan, Mingfei Li, Li Ning, and Shay Solomon. New doubling spanners: Better and simpler. In Fedor V. Fomin, Rusins Freivalds, Marta Z. Kwiatkowska, and David Peleg, editors, Proc. 40th ICALP, volume 7965 of Lecture Notes in Computer Science, pages 315–327. Springer, 2013. Full version in SICOMP’15. doi:10.1007/978-3-642-39206-1_27.
  • [12] B. Chandra, G. Das, G. Narasimhan, and J. Soares. New sparseness results on graph spanners. In Proceedings of the Eighth Annual Symposium on Computational Geometry, 1992.
  • [13] Hsien-Chih Chang, Jonathan Conroy, Hung Le, Lazar Milenković, Shay Solomon, and Cuong Than. Optimal Euclidean Tree Covers. In Wolfgang Mulzer and Jeff M. Phillips, editors, 40th International Symposium on Computational Geometry (SoCG 2024), volume 293 of Leibniz International Proceedings in Informatics (LIPIcs), pages 37:1–37:15, Dagstuhl, Germany, 2024. Schloss Dagstuhl – Leibniz-Zentrum für Informatik. doi:10.4230/LIPIcs.SoCG.2024.37.
  • [14] L. P. Chew. There is a planar graph almost as good as the complete graph. In Proceedings of the Second Annual Symposium on Computational Geometry, SCG ‘86, pages 169–177, 1986.
  • [15] K. Clarkson. Approximation algorithms for shortest path motion planning. In Proceedings of the Nineteenth Annual ACM Symposium on Theory of Computing, STOC ‘87, pages 56–65, 1987.
  • [16] G. Das, P. Heffernan, and G. Narasimhan. Optimally sparse spanners in 3-dimensional euclidean space. In Proceedings of the 9th Annual Symposium on Computational Geometry, SCG ’93, pages 53–62, 1993.
  • [17] G. Das and G. Narasimhan. A fast algorithm for constructing sparse Euclidean spanners. In Proc. of 10th SOCG, pages 132–139, 1994.
  • [18] G. Das, G. Narasimhan, and J. Salowe. A new way to weigh malnourished euclidean graphs. In Proceedings of the 6th Annual ACM-SIAM Symposium on Discrete Algorithms, SODA ’95, pages 215–222, 1995.
  • [19] Y. Dinitz, M. Elkin, and S. Solomon. Shallow-low-light trees, and tight lower bounds for Euclidean spanners. In Proc. of 49th FOCS, pages 519–528, 2008.
  • [20] Michael Elkin and Shay Solomon. Optimal euclidean spanners: Really short, thin, and lanky. J. ACM, 62(5):35:1–35:45, 2015. doi:10.1145/2819008.
  • [21] Arnold Filtser and Hung Le. Locality-sensitive orderings and applications to reliable spanners. In Proceedings of the 54th Annual ACM SIGACT Symposium on Theory of Computing, STOC ’22, pages 1066–1079. ACM, 2022. doi:10.1145/3519935.3520042.
  • [22] Jie Gao, Leonidas J. Guibas, and An Nguyen. Deformable spanners and applications. In Proceedings of the Twentieth Annual Symposium on Computational Geometry, SCG ’04, pages 190–199, New York, NY, USA, 2004. Association for Computing Machinery. doi:10.1145/997817.997848.
  • [23] L. Gottlieb and L. Roditty. Improved algorithms for fully dynamic geometric spanners and geometric routing. In Proc. of 19th SODA, pages 591–600, 2008.
  • [24] L. Gottlieb and L. Roditty. An optimal dynamic spanner for doubling metric spaces. In Proc. of 16th ESA, pages 478–489, 2008. Another version of this paper is available via http://cs.nyu.edu/˜adi/spanner2.pdf.
  • [25] Lee-Ad Gottlieb, Aryeh Kontorovich, and Robert Krauthgamer. Efficient regression in metric spaces via approximate lipschitz extension. IEEE Trans. Inf. Theory, 63(8):4838–4849, 2017. doi:10.1109/TIT.2017.2713820.
  • [26] J. Gudmundsson, C. Levcopoulos, and G. Narasimhan. Fast greedy algorithms for constructing sparse geometric spanners. SIAM J. Comput., 31(5):1479–1500, 2002. doi:10.1137/S0097539700382947.
  • [27] J. Gudmundsson, C. Levcopoulos, G. Narasimhan, and M. H. M. Smid. Approximate distance oracles for geometric graphs. In Proc. of 13th SODA, pages 828–837, 2002.
  • [28] J. Gudmundsson, C. Levcopoulos, G. Narasimhan, and M. H. M. Smid. Approximate distance oracles for geometric spanners. ACM Transactions on Algorithms, 4(1), 2008. doi:10.1145/1328911.1328921.
  • [29] J. Gudmundsson, G. Narasimhan, and M. H. M. Smid. Fast pruning of geometric spanners. In Proc. of 22nd STACS, pages 508–520, 2005.
  • [30] Joachim Gudmundsson, Christos Levcopoulos, Giri Narasimhan, and Michiel H. M. Smid. Approximate distance oracles revisited. In Proc. of 13th ISAAC, pages 357–368, 2002. doi:10.1007/3-540-36136-7_32.
  • [31] S Har-Peled and M Mendel. Fast construction of nets in low-dimensional metrics and their applications. SIAM Journal on Computing, 35(5):1148–1184, 2006. doi:10.1137/s0097539704446281.
  • [32] Y. Hassin and D. Peleg. Sparse communication networks and efficient routing in the plane. In Proc. of 19th PODC, pages 41–50, 2000.
  • [33] Omri Kahalon, Hung Le, Lazar Milenkovic, and Shay Solomon. Can’t see the forest for the trees: Navigating metric spaces by bounded hop-diameter spanners. In Alessia Milani and Philipp Woelfel, editors, Proc. PODC, pages 151–162. ACM, 2022. doi:10.1145/3519270.3538414.
  • [34] J. M. Keil. Approximating the complete euclidean graph. In Proceedings of the first Scandinavian Workshop on Algorithm Theory, SWAT ‘88, pages 208–213, 1988.
  • [35] J. M. Keil and C. A. Gutwin. Classes of graphs which approximate the complete Euclidean graph. Discrete and Computational Geometry, 7(1):13–28, 1992. doi:10.1007/BF02187821.
  • [36] An La and Hung Le. Dynamic locality sensitive orderings in doubling metrics. In Proceedings of the 57th Annual ACM Symposium on Theory of Computing, STOC ’25, pages 1178–1189, New York, NY, USA, 2025. Association for Computing Machinery. doi:10.1145/3717823.3718209.
  • [37] H. Le and S. Solomon. Near-optimal spanners for general graphs in (nearly) linear time. In Proceedings of the 2022 Annual ACM-SIAM Symposium on Discrete Algorithms, SODA ’22, pages 3332–3361. 2022. doi:10.1137/1.9781611977073.132.
  • [38] Hung Le and Shay Solomon. Truly optimal euclidean spanners. In David Zuckerman, editor, 60th IEEE Annual Symposium on Foundations of Computer Science, FOCS 2019, Baltimore, Maryland, USA, November 9-12, 2019, pages 1078–1100. IEEE Computer Society, 2019. doi:10.1109/FOCS.2019.00069.
  • [39] Hung Le and Shay Solomon. A unified framework for light spanners. In Proceedings of the 55th Annual ACM Symposium on Theory of Computing, STOC 2023, Orlando, FL, USA, June 20-23, 2023, pages 295–308. ACM, 2023. doi:10.1145/3564246.3585185.
  • [40] Hung Le, Shay Solomon, and Cuong Than. Optimal fault-tolerant spanners in euclidean and doubling metrics: Breaking the ω (log n) lightness barrier. In 2023 IEEE 64th Annual Symposium on Foundations of Computer Science (FOCS), pages 77–97. IEEE, 2023. doi:10.1109/FOCS57990.2023.00013.
  • [41] Hung Le and Cuong Than. Greedy spanners in euclidean spaces admit sublinear separators. In Joseph (Seffi) Naor and Niv Buchbinder, editors, Proc. SODA, pages 3287–3310. SIAM, 2022. doi:10.1137/1.9781611977073.130.
  • [42] Y. Mansour and D. Peleg. An approximation algorithm for min-cost network design. DIMACS Series in Discr. Math and TCS, 53:97–106, 2000.
  • [43] G. Narasimhan and M. Smid. Geometric Spanner Networks. Cambridge University Press, 2007.
  • [44] S. B. Rao and W. D. Smith. Approximating geometrical graphs via “spanners” and “banyans”. In Proceedings of the 30th Annual ACM Symposium on Theory of Computing, STOC ’98, pages 540–550, 1998. Full version at http://graphics.stanford.edu/courses/cs468-06-winter/Papers/rs-tsp.pdf. doi:10.1145/276698.276868.
  • [45] J. Ruppert and R. Seidel. Approximating the d-dimensional complete Euclidean graph. In Proceedings of the 3rd Canadian Conference on Computational Geometry, CCCG ‘91, pages 207–210, 1991.
  • [46] Jeffrey S. Salowe. On euclidean spanner graphs with small degree. In Proceedings of the Eighth Annual Symposium on Computational Geometry, Berlin, Germany, June 10-12, 1992, pages 186–191, 1992. doi:10.1145/142675.142716.
  • [47] M. Smid. The weak gap property in metric spaces of bounded doubling dimension. In Susanne Albers, Helmut Alt, and Stefan Näher, editors, Efficient Algorithms, pages 275–289. Springer-Verlag, 2009.
  • [48] Shay Solomon. From hierarchical partitions to hierarchical covers: optimal fault-tolerant spanners for doubling metrics. In Symposium on Theory of Computing, STOC 2014, New York, NY, USA, May 31 - June 03, 2014, pages 363–372, 2014. doi:10.1145/2591796.2591864.
  • [49] A. C. Yao. On constructing minimum spanning trees in k-dimensional spaces and related problems. SIAM Journal on Computing, 11(4):721–736, 1982. doi:10.1137/0211059.