Introduction
We consider the problem of embedding a union of Euclidean metric spaces into Euclidean space with low distortion. Formally, let ( X , d ) (X,d) ( X , d ) be a metric space, such that X X X is a union of k k k sets U 1 , … , U k U_1,\ldots,U_k U 1 , … , U k and each U i U_i U i can be embedded isometrically into Euclidean space. Our goal is to find the smallest distortion D k D_k D k , such that the union U 1 ∪ ⋯ ∪ U k U_1\cup\cdots\cup U_k U 1 ∪ ⋯ ∪ U k embeds into ℓ 2 \ell_2 ℓ 2 with distortion D k D_k D k . When every U i U_i U i is a single point this is exactly the setting of Bourgain’s Theorem [4 ] , which states that every n n n -point metric space can be embedded into Euclidean space with distortion O ( log n ) O(\log n) O ( log n ) . More generally, our question belongs to the study of local-global properties of metric embeddings initiated by Arora, Lovasz, Newman, Rabani, Rabinovich, and Vempala [1 ] (see also [6 ] ), and asks whether the distortion of X X X is governed only by k k k .
The special case k = 2 k=2 k = 2 was asked by Assaf Naor and settled by Makarychev and Makarychev [21 ] , who showed that a union of two Euclidean sets embeds into ℓ 2 \ell_2 ℓ 2 with distortion at most 8.93. They then posed the problem for general k k k as an open question [21 Question 6] . On one hand, iterating their result gives an upper bound of D k ≤ 2 O ( k ) D_k\le2^{O(k)} D k ≤ 2 O ( k ) . On the other hand, when each U i U_i U i is a single point we have the classical lower bound of D k = Ω ( log k ) D_k=\Omega(\log k) D k = Ω ( log k ) for Bourgain’s Theorem [19 , 2 ] .
In this work, we show that a union of k k k Euclidean sets can be embedded into Euclidean space with distortion D k = O ( k 1.128 log k ) D_k=O(k^{1.128}\log k) D k = O ( k 1.128 log k ) , an exponential improvement over the previous bound, with no dependence on the number of points in X X X or any U i U_i U i . Given the additional property that every U i U_i U i stays Euclidean after adjoining any other two points from X X X , we can further establish that the distortion of embedding X X X is O ( log k ) O(\log k) O ( log k ) , thus, matching the lower bound. Finally, if each U i U_i U i remains Euclidean after adjoining any one other point from X X X , then we can obtain distortion of O ( k 0.895 log k ) O(k^{0.895}\log k) O ( k 0.895 log k ) .
Our approach contains two reductions that may be of independent interest. First, the distortion of a union of Euclidean sets is determined, up to an O ( log k ) O(\log k) O ( log k ) factor, by how well Lipschitz maps on the individual sets extend to the whole space. Second, bounding these Lipschitz extension constants is equivalent to constructing good Euclidean metric sparsifiers.
Model and Main Results
Distortion. We say that a metric space ( X , d ) (X,d) ( X , d ) embeds in ℓ 2 \ell_2 ℓ 2 with distortion D D D if there is a function f : X → ℓ 2 f:X\to\ell_2 f : X → ℓ 2 , such that for every x , y ∈ X x,y\in X x , y ∈ X :
d ( x , y ) ≤ ∥ f ( x ) − f ( y ) ∥ 2 ≤ D ⋅ d ( x , y ) . (1) d(x,y) \le\lVert f(x)-f(y)\rVert_2 \le D\cdot d(x,y). \tag*{(1)} d ( x , y ) ≤ ∥ f ( x ) − f ( y ) ∥ 2 ≤ D ⋅ d ( x , y ) . (1)
Then, the distortion c 2 ( X ) c_2(X) c 2 ( X ) of ( X , d ) (X,d) ( X , d ) is the smallest such scalar D D D .
Euclidean Classes. In a metric space ( X , d ) (X,d) ( X , d ) , a subset V ⊆ X V\subseteq X V ⊆ X is Euclidean in X X X if it embeds isometrically into ℓ 2 \ell_2 ℓ 2 . Moreover, we say a subset V ⊆ X V\subseteq X V ⊆ X is h h h -strongly Euclidean in X X X if, for every h h h -element subset S ⊆ X S\subseteq X S ⊆ X , the set V ∪ S V\cup S V ∪ S embeds isometrically into ℓ 2 \ell_2 ℓ 2 . That is, an h h h -strongly Euclidean set is a Euclidean set that remains Euclidean if we add h h h arbitrarily chosen points.
The parameter h h h interpolates between the condition that each U i U_i U i is Euclidean (h = 0 h=0 h = 0 ) and that X X X is Euclidean itself (h = ∣ X ∣ h=|X| h = ∣ X ∣ ). In the spirit of local-to-global phenomena [1 , 6 ] , one can ask how large h h h has to be before the distortion of X X X becomes the optimal Θ ( log k ) \Theta(\log k) Θ ( log k ) . Our first main result shows that h = 2 h=2 h = 2 is already sufficient.
Theorem 1.1. Let ( X , d ) (X,d) ( X , d ) be a metric space, such that X = U 1 ∪ ⋯ ∪ U k X=U_1\cup\cdots\cup U_k X = U 1 ∪ ⋯ ∪ U k and each U i U_i U i is 2-strongly Euclidean in X X X . Then X X X embeds into ℓ 2 \ell_2 ℓ 2 with distortion O ( log k ) O(\log k) O ( log k ) .
Observe that when each U i U_i U i consists of a single point, Theorem 1.1 recovers Bourgain’s Theorem [4 ] . Because an h h h -strongly Euclidean set is 2-strongly Euclidean for every h ≥ 2 h\ge2 h ≥ 2 , the same bound holds for unions of h h h -strongly Euclidean sets with h ≥ 3 h\ge3 h ≥ 3 . Therefore, what remains are, indeed, only unions of Euclidean sets and unions of 1-strongly Euclidean sets. In the sections that follow we first develop the proof of Theorem 1.1 , and later adapt what is useful to the Euclidean and 1-strongly Euclidean settings.
In the general case of a union of k k k Euclidean sets we obtain the following subquadratic upper bound on the distortion.
Theorem 1.2. Let ( X , d ) (X,d) ( X , d ) be a metric space, such that X = U 1 ∪ ⋯ ∪ U k X=U_1\cup\cdots\cup U_k X = U 1 ∪ ⋯ ∪ U k and each U i U_i U i is Euclidean in X X X . Then X X X embeds into ℓ 2 \ell_2 ℓ 2 with distortion O ( k 1.128 log k ) O(k^{1.128}\log k) O ( k 1.128 log k ) . Moreover, for the special case k = 2 k=2 k = 2 we get an improved upper bound of c 2 ( X ) < 3.997 c_2(X)<3.997 c 2 ( X ) < 3.997 .
For the case of 1-strongly Euclidean sets, we get the following sublinear distortion.
Theorem 1.3. Let ( X , d ) (X,d) ( X , d ) be a metric space, such that X = U 1 ∪ ⋯ ∪ U k X=U_1\cup\cdots\cup U_k X = U 1 ∪ ⋯ ∪ U k and each U i U_i U i is 1-strongly Euclidean in X X X . Then X X X embeds into ℓ 2 \ell_2 ℓ 2 with distortion O ( k 0.895 log k ) O(k^{0.895}\log k) O ( k 0.895 log k ) .
Technical Overview
Union of 2-strongly Euclidean sets. In § 2 , we prove Theorem 1.1 in two steps:
For each U i U_i U i we construct a D D D -terminal embedding, a 1-Lipschitz map φ i : X → ℓ 2 \varphi_i:X\to\ell_2 φ i : X → ℓ 2 that contracts distances between U i U_i U i and every point of X X X by at most a factor of D D D . This is where we use the 2-strongly Euclidean property of the sets U i U_i U i .
We combine the k k k terminal embeddings into a single embedding for X X X . In the combined embedding, each point uses only the terminal embedding of the set closest to it, and we sample subcollections of the sets as in the algorithmic proof of Bourgain’s Theorem [19 ] .
Overall, this yields an embedding of X X X into ℓ 2 \ell_2 ℓ 2 with distortion O ( log k ) O(\log k) O ( log k ) . We note that our notion of terminal embedding is slightly different from the one considered by [8 ] because they do not require the maps to be Lipschitz on the whole space, as we do here.
From distortion to Lipschitz extension constants. Only the first step above depends on the sets being 2-strongly Euclidean, while combining terminal embeddings works for any union of Euclidean sets as long as each set admits a terminal embedding. In § 3 , we show that any L L L -Lipschitz map X → ℓ 2 X\to\ell_2 X → ℓ 2 that does not contract distances on U i U_i U i can be augmented into an O ( L ) O(L) O ( L ) -terminal embedding for U i U_i U i (Lemma 3.1 ). So now the question becomes how well an isometry of U i U_i U i extends to a Lipschitz map on all of X X X . This is captured by the Lipschitz extension constant e ( U i ) e(U_i) e ( U i ) , the least Λ ≥ 1 \Lambda\ge1 Λ ≥ 1 such that every Lipschitz map f : ( U i , d ) → ℓ 2 f:(U_i,d)\to\ell_2 f : ( U i , d ) → ℓ 2 extends to a Lipschitz map F : ( X , d ) → ℓ 2 F:(X,d)\to\ell_2 F : ( X , d ) → ℓ 2 with Lip ( F ) ≤ Λ Lip ( f ) \operatorname{Lip}(F)\le\Lambda\operatorname{Lip}(f) Lip ( F ) ≤ Λ Lip ( f ) . Indeed, since U i U_i U i is Euclidean, by Kirszbraun’s Theorem [14 ] it suffices to extend an isometry of U i U_i U i in order to bound e ( U i ) e(U_i) e ( U i ) . The formal result connecting distortion to Lipschitz extension constants is the following.
Theorem 1.4. Let X = U 1 ∪ ⋯ ∪ U k X=U_1\cup\cdots\cup U_k X = U 1 ∪ ⋯ ∪ U k with each U i U_i U i Euclidean in X X X , and let e max = max i e ( U i ) e_{\max}=\max_i e(U_i) e m a x = max i e ( U i ) . Then
e max ≤ c 2 ( X ) ≤ O ( e max log ( k + 1 ) ) . e_{\max}\le c_2(X)\le O(e_{\max}\log(k+1)). e m a x ≤ c 2 ( X ) ≤ O ( e m a x log ( k + 1 )) . For instance, for a union of 2-strongly Euclidean sets, we can interpret the analysis of Theorem 1.1 as showing that e max = Θ ( 1 ) e_{\max}=\Theta(1) e m a x = Θ ( 1 ) .
Lipschitz extension constants for Euclidean and 1-strongly Euclidean sets. In § 4 , we bound the Lipschitz extension constants using metric sparsifiers [20 ] . Given a weighted graph G G G on X X X and a set of terminals A ⊆ X A\subseteq X A ⊆ X , a metric sparsifier of G G G is a weighted graph H H H on A A A alone, such that for every metric ρ \rho ρ on A A A the cost of H H H on ρ \rho ρ approximates, up to a factor Q Q Q called the quality , the least cost of G G G over all extensions of ρ \rho ρ to X X X . We take A = U i A=U_i A = U i and measure cost by the quadratic energy E G ( d 2 ) = ∑ x < y g x y d ( x , y ) 2 E_G(d^2)=\sum_{x<y}g_{xy}d(x,y)^2 E G ( d 2 ) = ∑ x < y g x y d ( x , y ) 2 , while only allowing extensions that are Euclidean on every U j U_j U j . In § 4.1 , we show that in this setting the best sparsifier quality for every G G G is exactly the square of the Lipschitz extension constant, where both quantities are taken uniformly over all metrics on X X X that are Euclidean on every U j U_j U j (Theorem 4.2 ). Bounding e max e_{\max} e m a x thus reduces to constructing a good metric sparsifier.
The following upper bounds, combined with Theorem 1.4 , give Theorems 1.2 and 1.3 for finite X X X . We take care of the finiteness assumption in Section C.
Lemma 1.5. Let ( X , d ) (X,d) ( X , d ) be a finite metric space, such that X = U 1 ∪ ⋯ ∪ U k X=U_1\cup\cdots\cup U_k X = U 1 ∪ ⋯ ∪ U k . If each U i U_i U i is Euclidean in X X X , then e max = O ( k 1.128 ) e_{\max}=O(k^{1.128}) e m a x = O ( k 1.128 ) . If each U i U_i U i is 1-strongly Euclidean in X X X , then e max = O ( k 0.895 ) e_{\max}=O(k^{0.895}) e m a x = O ( k 0.895 ) .
In the construction for Lemma 1.5 we bound the quality of the canonical sparsifier for U i U_i U i by adding the other sets back one at a time. We show that adding one Euclidean set costs a constant factor in the Lipschitz constant, and averaging the extensions over the set added last decreases the factor at the step with r r r sets down to 1 + η / r 1+\eta/r 1 + η / r (Lemma 4.8 ). The product of these factors is O ( k η ) O(k^\eta) O ( k η ) , where the exponent η \eta η is determined by the best trade-off between the growth on pairs of points inside the added set and on pairs crossing it, giving η 0 = 1.1275 … < 1.128 \eta_0=1.1275\ldots<1.128 η 0 = 1.1275 … < 1.128 for Euclidean sets and η 1 = 0.8941 … < 0.895 \eta_1 = 0.8941\ldots< 0.895 η 1 = 0.8941 … < 0.895 for 1-strongly Euclidean sets (§ 4.5 ). We leave the improved upper bound for the case k = 2 k = 2 k = 2 to Section B.
Lower bounds on Lipschitz extension constants. Theorem 1.4 bounds the distortion with the maximal Lipschitz extension constant e max e_{\max} e m a x up to a factor of O ( log k ) O(\log k) O ( log k ) . Meanwhile, the analysis of Theorem 1.1 shows that e max e_{\max} e m a x is an absolute constant for unions of 2-strongly Euclidean sets. What happens for unions of Euclidean sets, or of 1-strongly Euclidean sets? We find that the extension constant can be Ω ( log k ) \Omega(\sqrt{\log k}) Ω ( log k ) for some unions of Euclidean sets, and Ω ( log k 4 ) \Omega(\sqrt[4]{\log k}) Ω ( 4 log k ) for 1-strongly Euclidean. The lower bounds we obtain here do not improve on the general Ω ( log k ) \Omega(\log k) Ω ( log k ) lower bound for distortion. They do show, however, that the reduction to Lipschitz extensions cannot by itself give O ( log k ) O(\log k) O ( log k ) for unions of Euclidean sets.
Lemma 1.6. There are unions of k k k Euclidean sets with e max = Ω ( log k ) e_{\max} = \Omega(\sqrt{\log k}) e m a x = Ω ( log k ) , and unions of k k k 1-strongly Euclidean sets with e max = Ω ( log k 4 ) e_{\max} = \Omega(\sqrt[4]{\log k}) e m a x = Ω ( 4 log k ) .
In § 5 for unions of Euclidean sets, we build X X X from a graph with k k k vertices and girth Θ ( log k ) \Theta(\log k) Θ ( log k ) , under a specially-designed norm on R E \mathbb{R}^{E} R E , so that the cycle space of the graph is Euclidean, with one set per vertex. Then extending an isometry of one set is the same as a linear projection onto the cycle space. By a theorem of Lindenstrauss [17 ] , its Lipschitz constant is controlled by the relative projection constant of the cycle space. In § 6 , for unions of 1-strongly Euclidean sets we use a variant of the twisted cube of Naor and Rabani [24 ] .
From isometric embeddings to constant-distortion embeddings. Finally, in § 7 , we relax the assumption that each U i U_i U i is isometric to a subset of ℓ 2 \ell_2 ℓ 2 and only require that it embeds with constant distortion C i C_i C i . Similarly to the isometric setting, the distortion is naively upper bounded by O ( k ) ∏ i = 1 k C i O(k) \prod_{i=1}^{k} C_i O ( k ) ∏ i = 1 k C i , and the second term is itself exponential in k k k when all C i C_i C i equal a common constant. Fortunately, we can do much better than that.
Theorem 1.7. Let ( X , d ) (X,d) ( X , d ) be a metric space, such that X = U 1 ∪ ⋯ ∪ U k X = U_1 \cup\cdots\cup U_k X = U 1 ∪ ⋯ ∪ U k and each U i U_i U i embeds in ℓ 2 \ell_2 ℓ 2 with distortion C i ≥ 1 C_i \ge1 C i ≥ 1 , and let C max = max i C i C_{\max} = \max_i C_i C m a x = max i C i . Then X X X embeds into ℓ 2 \ell_2 ℓ 2 with distortion O ( C max k η ( C max ) log k ) O(C_{\max} k^{\eta(C_{\max})} \log k) O ( C m a x k η ( C m a x ) log k ) , where η ( c ) ≤ c − 1 + π / c \eta(c) \le c - 1 + \pi/c η ( c ) ≤ c − 1 + π / c for every c ≥ 1 c \ge1 c ≥ 1 and η ( 1 ) = η 0 \eta(1) = \eta_0 η ( 1 ) = η 0 .
In the isometric setting the 2-strongly Euclidean property is what gave the tight O ( log k ) O(\log k) O ( log k ) distortion in Theorem 1.1. We might hope for a similar effect here, by assuming that every U i ∪ { x , y } U_i \cup\{x,y\} U i ∪ { x , y } embeds with constant distortion. However, by [21 ] , Theorem 1.1, if U i U_i U i embeds with distortion C C C then U i ∪ { x , y } U_i \cup\{x,y\} U i ∪ { x , y } embeds into ℓ 2 \ell_2 ℓ 2 with distortion O ( C ) O(C) O ( C ) because { x , y } \{x,y\} { x , y } is itself Euclidean. Therefore, up to constants, this condition already follows from the assumption that each U i U_i U i embeds with constant distortion. As a result, adjoining two points provides no meaningful improvement, unlike in the isometric setting.
Related Work
Bourgain’s Theorem. While the initial result of Bourgain [4 ] was existential, Linial, London and Rabinovich [19 ] gave an algorithmic construction and initiated the use of low-distortion embeddings in approximation algorithms, with early example applications being [19 , 2 , 10 ] . Beyond approximation algorithms, Bourgain’s Theorem has been used for Internet latency estimation [15 ] , differential privacy [11 ] , and designing position-aware embeddings for graph neural networks [25 ] . We refer the reader to [12 , 18 ] for surveys on low-distortion embeddings.
Terminal embeddings. Elkin, Filtser and Neiman [8 ] considered a different generalization of Bourgain’s Theorem. Given a set of terminals T ⊆ X T \subseteq X T ⊆ X , a terminal embedding approximately preserves the distance from every point of X X X to every terminal, and they obtained distortion O ( log ∣ T ∣ ) O(\log|T|) O ( log ∣ T ∣ ) . When X X X is itself Euclidean, Mahabadi, Makarychev, Makarychev and Razenshteyn [22 ] showed that terminal embeddings can be obtained by extending an embedding of the terminals to the whole space with distortion 1 + ε 1+\varepsilon 1 + ε and O ( ε − 4 log ∣ T ∣ ) O(\varepsilon^{-4}\log|T|) O ( ε − 4 log ∣ T ∣ ) dimensions. Later Narayanan and Nelson [23 ] improved the number of dimensions to the tight O ( ε − 2 log ∣ T ∣ ) O(\varepsilon^{-2}\log|T|) O ( ε − 2 log ∣ T ∣ ) , and Cherapanamjeri and Nelson [7 ] showed that such embeddings can be computed in time sublinear in ∣ T ∣ |T| ∣ T ∣ .
Union of 2-strongly Euclidean Sets
In this section we prove Theorem 1.1 . We show how to construct terminal embeddings for 2-strongly Euclidean sets, and then combine those using a procedure inspired by Bourgain’s Theorem.
Preliminaries
Notation. For 1 ≤ p < ∞ 1 \le p < \infty 1 ≤ p < ∞ , we write A ⊕ p B A \oplus_p B A ⊕ p B for the direct sum with norm ∥ ( a , b ) ∥ = ( ∥ a ∥ p + ∥ b ∥ p ) 1 / p \|(a,b)\| = (\|a\|^p+\|b\|^p)^{1/p} ∥ ( a , b ) ∥ = ( ∥ a ∥ p + ∥ b ∥ p ) 1/ p , and ⨁ p , i = 1 r A i \bigoplus_{p,i=1}^r A_i ⨁ p , i = 1 r A i for a finite sum of this kind. For n < m n<m n < m , we will identify ℓ 2 n \ell_2^n ℓ 2 n with the subspace of ℓ 2 m \ell_2^m ℓ 2 m spanned by the first n n n coordinates. For a finite- or infinite-dimensional space ℓ 2 n \ell_2^n ℓ 2 n , the space ℓ 2 n ⊕ 2 R \ell_2^n \oplus_2 \mathbb{R} ℓ 2 n ⊕ 2 R is the Euclidean or Hilbert space obtained from ℓ 2 n \ell_2^n ℓ 2 n by adding an extra coordinate; we will identify ℓ 2 n \ell_2^n ℓ 2 n with the subspace of ℓ 2 n ⊕ 2 R \ell_2^n \oplus_2 \mathbb{R} ℓ 2 n ⊕ 2 R where the extra coordinate is 0.
Two-sided Kirszbraun theorem for combining terminal embeddings. This lemma allows us to smoothly combine terminal embeddings for different sets U i U_i U i .
Lemma 2.1 ([3 Cor. 14] ). Let S S S be a nonempty subset of ℓ 2 m \ell_2^m ℓ 2 m , where m m m may be finite or infinite. Define R ( x ) = dist ( x , S ) R(x)=\operatorname{dist}(x,S) R ( x ) = dist ( x , S ) . There exists a map ψ : ℓ 2 m → ℓ 2 \psi:\ell_2^m\to\ell_2 ψ : ℓ 2 m → ℓ 2 such that:
ψ ( x ) = 0 \psi(x)=0 ψ ( x ) = 0 for every x ∈ S x\in S x ∈ S ;
for every x , y ∈ ℓ 2 m x,y\in\ell_2^m x , y ∈ ℓ 2 m ,
∥ ψ ( x ) − ψ ( y ) ∥ = Θ ( min { ∥ x − y ∥ , R ( x ) + R ( y ) } ) . \|\psi(x)-\psi(y)\|=\Theta\left(\min\{\|x-y\|,R(x)+R(y)\}\right). ∥ ψ ( x ) − ψ ( y ) ∥ = Θ ( min { ∥ x − y ∥ , R ( x ) + R ( y )} ) .
The hidden constants in the Θ \Theta Θ notation are absolute and do not depend on S S S and m m m .
Fix a finite or infinite n n n . Let m = n + 1 m=n+1 m = n + 1 . Applying Lemma 2.1 to S = ℓ 2 n = { ( u , 0 ) : u ∈ ℓ 2 n } ⊂ ℓ 2 n ⊕ 2 R S=\ell_2^n=\{(u,0):u\in\ell_2^n\}\subset\ell_2^n\oplus_2\mathbb{R} S = ℓ 2 n = {( u , 0 ) : u ∈ ℓ 2 n } ⊂ ℓ 2 n ⊕ 2 R , we get the following corollary.
Corollary 2.2. For every n n n , there exists a map ψ : ℓ 2 n × R → ℓ 2 \psi:\ell_2^n\times\mathbb{R}\to\ell_2 ψ : ℓ 2 n × R → ℓ 2 with the following properties:
ψ ( u , 0 ) = 0 \psi(u,0)=0 ψ ( u , 0 ) = 0 for all u ∈ ℓ 2 n u\in\ell_2^n u ∈ ℓ 2 n .
∥ ψ ( u , a ) − ψ ( v , b ) ∥ = Θ ( min { ∥ u − v ∥ + ∣ a − b ∣ , a + b } ) \|\psi(u,a)-\psi(v,b)\|=\Theta(\min\{\|u-v\|+|a-b|,a+b\}) ∥ ψ ( u , a ) − ψ ( v , b ) ∥ = Θ ( min { ∥ u − v ∥ + ∣ a − b ∣ , a + b }) for all u , v ∈ ℓ 2 n u,v\in\ell_2^n u , v ∈ ℓ 2 n and a , b ≥ 0 a,b\ge0 a , b ≥ 0 , where the constants in the Θ \Theta Θ -notation are absolute.
We will construct maps ψ \psi ψ for different values of n n n , but use the same notation ψ \psi ψ for all of them, since the value of n n n will be clear from context. Moreover, all these maps can be viewed as restrictions of the map ψ \psi ψ for n = ∞ n=\infty n = ∞ to finite values of n n n .
Reduction to finite spaces. For bounding the distortion, it is enough to consider only finite metric spaces. Let X = U 1 ∪ ⋯ ∪ U k X=U_1\cup\cdots\cup U_k X = U 1 ∪ ⋯ ∪ U k be as in any of Theorems 1.1 to 1.4 and 1.7. Since each U i U_i U i embeds into ℓ 2 \ell_2 ℓ 2 with finite distortion, X X X is separable. Then the compactness argument of [21 Lemma 2.9, part I] shows that X X X embeds with distortion D D D as long as every finite subset of X X X does. Moreover, every finite S ⊆ X S\subseteq X S ⊆ X is covered by the nonempty sets S ∩ U i S\cap U_i S ∩ U i , which inherit the Euclidean property of U i U_i U i . Further details are in Section C. For the rest of the paper we assume that all metric spaces are finite.
Terminal embeddings for 2-strongly Euclidean subsets
We now show that the 2-strongly Euclidean property implies the existence of 1-terminal embeddings.
Definition 2.3. Let ( X , d ) (X,d) ( X , d ) be a metric space, and let V ⊆ X V \subseteq X V ⊆ X be a nonempty subset. For D ≥ 1 D \ge1 D ≥ 1 , a map φ : X → ℓ 2 m \varphi:X \to\ell_2^m φ : X → ℓ 2 m , where m m m may be finite or infinite, is a D D D -terminal embedding for V V V if it satisfies the following properties:
The map φ \varphi φ is 1-Lipschitz.
The map φ \varphi φ contracts distances from points in V V V by at most a factor D D D : for every x ∈ V x \in V x ∈ V and y ∈ X y \in X y ∈ X ,
∥ φ ( x ) − φ ( y ) ∥ ≥ d ( x , y ) / D . \lVert\varphi(x)-\varphi(y)\rVert\ge d(x,y)/D. ∥ φ ( x ) − φ ( y )∥ ≥ d ( x , y ) / D .
In particular, when D = 1 D=1 D = 1 , φ \varphi φ preserves all distances from points in V V V .
We use a variant of terminal embeddings introduced by [8 ] . In their terminology, we require our D D D -terminal embeddings for V V V to have terminal distortion at most D D D and to be 1-Lipschitz on all of X X X ; the latter condition is not required in their original definition.
Let us say that an isometric embedding φ : Z → ℓ 2 r \varphi:Z \to\ell_2^r φ : Z → ℓ 2 r (where r r r is finite or infinite) is tight if the closed affine span L L L of φ ( Z ) \varphi(Z) φ ( Z ) coincides with ℓ 2 r \ell_2^r ℓ 2 r . If φ \varphi φ is not tight, we can transform it to a tight isometric embedding by restricting the target space to L L L .
We will use the following simple geometric fact.
Fact 2.4 (Rigidity of ℓ 2 \ell_2 ℓ 2 embeddings). Consider a tight isometric embedding φ : Z → ℓ 2 n \varphi:Z \to\ell_2^n φ : Z → ℓ 2 n . Then every isometric embedding φ x ′ \varphi'_x φ x ′ of Z ∪ { x } Z \cup\{x\} Z ∪ { x } into ℓ 2 \ell_2 ℓ 2 can be transformed into an isometric embedding φ x \varphi_x φ x into ℓ 2 ⊕ 2 R \ell_2 \oplus_2 \mathbb{R} ℓ 2 ⊕ 2 R such that
φ x ( z ) = φ ( z ) ⊕ 2 0 \varphi_x(z)=\varphi(z)\oplus_2 0 φ x ( z ) = φ ( z ) ⊕ 2 0 for z ∈ Z z \in Z z ∈ Z (that is, φ x \varphi_x φ x is an extension of φ \varphi φ );
φ x ( x ) = u ⊕ 2 a \varphi_x(x)=u\oplus_2 a φ x ( x ) = u ⊕ 2 a for some u ∈ ℓ 2 n u \in\ell_2^n u ∈ ℓ 2 n and a ≥ 0 a \ge0 a ≥ 0 .
See Figure 1 for the geometry of the terminal embedding.
Figure 1. The terminal embedding places x x x above the subspace containing φ 0 ( V ) \varphi_0(V) φ 0 ( V ) , with projection u u u and height a a a . The projection u u u need not belong to φ 0 ( V ) \varphi_0(V) φ 0 ( V ) .
Claim 2.5. If V V V is a nonempty 2-strongly Euclidean subset of X X X , then there exists a 1-terminal embedding φ : X → ℓ 2 \varphi:X \to\ell_2 φ : X → ℓ 2 for V V V .
Proof. Since V V V is 2-strongly Euclidean in X X X , it isometrically embeds into some ℓ 2 r \ell_2^r ℓ 2 r . Let φ 0 : V → ℓ 2 r \varphi_0:V \to\ell_2^r φ 0 : V → ℓ 2 r be a tight isometric embedding. For every x ∈ X x \in X x ∈ X , we are also given that V ∪ { x } V \cup\{x\} V ∪ { x } isometrically embeds into ℓ 2 \ell_2 ℓ 2 . By Fact 2.4 , there is an embedding φ x : V ∪ { x } → ℓ 2 ⊕ 2 R \varphi_x:V \cup\{x\} \to\ell_2 \oplus_2 \mathbb{R} φ x : V ∪ { x } → ℓ 2 ⊕ 2 R extending φ 0 \varphi_0 φ 0 : φ x ( z ) = φ 0 ( z ) ⊕ 2 0 \varphi_x(z)=\varphi_0(z)\oplus_2 0 φ x ( z ) = φ 0 ( z ) ⊕ 2 0 for z ∈ V z \in V z ∈ V and φ x ( x ) = u ⊕ 2 a \varphi_x(x)=u\oplus_2 a φ x ( x ) = u ⊕ 2 a for u ∈ ℓ 2 r u \in\ell_2^r u ∈ ℓ 2 r and a ≥ 0 a \ge0 a ≥ 0 . We are ready to define φ \varphi φ :
φ ( x ) = { φ 0 ( x ) , if x ∈ V φ x ( x ) , if x ∉ V \varphi(x)=
\begin{cases}
\varphi_0(x), & \text{if } x\in V \\
\varphi_x(x), & \text{if } x\notin V
\end{cases} φ ( x ) = { φ 0 ( x ) , φ x ( x ) , if x ∈ V if x ∈ / V
It is immediate that φ \varphi φ preserves distances from points in V V V : for x ∈ V x\in V x ∈ V and y ∈ X y\in X y ∈ X ,
∥ φ ( x ) − φ ( y ) ∥ = ∥ φ y ( x ) − φ y ( y ) ∥ = d ( x , y ) . \|\varphi(x)-\varphi(y)\|=\|\varphi_y(x)-\varphi_y(y)\|=d(x,y). ∥ φ ( x ) − φ ( y ) ∥ = ∥ φ y ( x ) − φ y ( y ) ∥ = d ( x , y ) .
It remains to show that φ \varphi φ is 1-Lipschitz. That is,
∥ φ ( x ) − φ ( y ) ∥ ≤ d ( x , y ) . \|\varphi(x)-\varphi(y)\|\le d(x,y). ∥ φ ( x ) − φ ( y ) ∥ ≤ d ( x , y ) .
The claim is trivial if at least one of the points x x x and y y y is in V V V (in this case, we have an equality). So let us assume that x , y ∉ V x,y\notin V x , y ∈ / V .
Using Fact 2.4 , extend φ x \varphi_x φ x to a map φ x y : V ∪ { x , y } → ℓ 2 r ⊕ 2 R ⊕ 2 R \varphi_{xy}:V\cup\{x,y\}\to\ell_2^r\oplus_2\mathbb{R}\oplus_2\mathbb{R} φ x y : V ∪ { x , y } → ℓ 2 r ⊕ 2 R ⊕ 2 R . We have
φ x y ( x ) = ( φ x ( x ) , 0 ) = ( φ ( x ) , 0 ) = ( u , a , 0 ) , \varphi_{xy}(x)=(\varphi_x(x),0)=(\varphi(x),0)=(u,a,0), φ x y ( x ) = ( φ x ( x ) , 0 ) = ( φ ( x ) , 0 ) = ( u , a , 0 ) ,
φ x y ( y ) = ( v , b cos α , b sin α ) \varphi_{xy}(y)=(v,b\cos\alpha,b\sin\alpha) φ x y ( y ) = ( v , b cos α , b sin α )
for some u , v ∈ ℓ 2 r u,v\in\ell_2^r u , v ∈ ℓ 2 r , a , b ≥ 0 a,b\ge0 a , b ≥ 0 and α ∈ [ 0 , 2 π ) \alpha\in[0,2\pi) α ∈ [ 0 , 2 π ) . Since both φ y \varphi_y φ y and φ x y \varphi_{xy} φ x y are isometric embeddings of V ∪ { y } V\cup\{y\} V ∪ { y } , we further have
φ ( y ) = φ y ( y ) = ( v , b ) . \varphi(y)=\varphi_y(y)=(v,b). φ ( y ) = φ y ( y ) = ( v , b ) .
Now we upper bound the distance between φ ( x ) \varphi(x) φ ( x ) and φ ( y ) \varphi(y) φ ( y ) :
∥ φ ( x ) − φ ( y ) ∥ 2 = ∥ u − v ∥ 2 + a 2 + b 2 − 2 a b ≤ ∥ u − v ∥ 2 + a 2 + b 2 − 2 a b cos α = ∥ φ x y ( x ) − φ x y ( y ) ∥ 2 = d ( x , y ) 2 \|\varphi(x)-\varphi(y)\|^2=\|u-v\|^2+a^2+b^2-2ab\le\|u-v\|^2+a^2+b^2-2ab\cos\alpha=\|\varphi_{xy}(x)-\varphi_{xy}(y)\|^2=d(x,y)^2 ∥ φ ( x ) − φ ( y ) ∥ 2 = ∥ u − v ∥ 2 + a 2 + b 2 − 2 ab ≤ ∥ u − v ∥ 2 + a 2 + b 2 − 2 ab cos α = ∥ φ x y ( x ) − φ x y ( y ) ∥ 2 = d ( x , y ) 2
where we used that φ x y \varphi_{xy} φ x y is an isometry.
Claim 2.5 immediately gives an embedding with distortion at most k \sqrt{k} k under the assumptions of Theorem 1.1 . This guarantee is much weaker than the desired one, but its proof is immediate.
Construct a 1-terminal embedding φ i \varphi_i φ i for each U i U_i U i and let F = φ 1 ⊕ 2 ⋯ ⊕ 2 φ k F=\varphi_1\oplus_2\cdots\oplus_2\varphi_k F = φ 1 ⊕ 2 ⋯ ⊕ 2 φ k . Since all maps φ i \varphi_i φ i are 1-Lipschitz, Lip ( F ) ≤ k \operatorname{Lip}(F)\le\sqrt{k} Lip ( F ) ≤ k . On the other hand, for x ∈ U i x\in U_i x ∈ U i and y ∈ X y\in X y ∈ X , we have
∥ F ( x ) − F ( y ) ∥ ≥ ∥ φ i ( x ) − φ i ( y ) ∥ = d ( x , y ) . \|F(x)-F(y)\|\ge\|\varphi_i(x)-\varphi_i(y)\|=d(x,y). ∥ F ( x ) − F ( y ) ∥ ≥ ∥ φ i ( x ) − φ i ( y ) ∥ = d ( x , y ) .
Therefore, the distortion is at most k \sqrt{k} k .
The construction above concatenates terminal embeddings for the individual sets U i U_i U i . Its analogue for a k k k -point metric space concatenates the distance functions to individual points. For a metric space X = { x 1 , … , x k } X=\{x_1,\ldots,x_k\} X = { x 1 , … , x k } , there is a simple embedding F ( y ) = ( d ( x 1 , y ) , d ( x 2 , y ) , … , d ( x k , y ) ) F(y)=(d(x_1,y),d(x_2,y),\ldots,d(x_k,y)) F ( y ) = ( d ( x 1 , y ) , d ( x 2 , y ) , … , d ( x k , y )) of X X X into ℓ 2 k \ell_2^k ℓ 2 k with distortion k \sqrt{k} k . Bourgain’s proof improves on this latter construction: to get an asymptotically optimal distortion O ( log k ) O(\log k) O ( log k ) , we concatenate distance functions to different subsets A A A of X X X . Informally, the advantage of this approach is that we have much tighter control over the Lipschitz constant of the obtained map.
Now our goal is to generalize the proof of Bourgain’s theorem to our setting. Our next step is to define a map for a subcollection of U 1 , … , U k U_1,\ldots,U_k U 1 , … , U k .
Combining terminal embeddings
In this section, we construct a map f f f for a given collection of nonempty subsets V 1 , … , V r V_1,\ldots,V_r V 1 , … , V r of X X X . We do not assume that V 1 , … , V r V_1,\ldots,V_r V 1 , … , V r cover X X X . In the proof of Theorem 1.1 , we will apply this construction to a randomly chosen subcollection of U 1 , … , U k U_1,\ldots,U_k U 1 , … , U k .
Let D ≥ 1 D \ge1 D ≥ 1 , and suppose that each V i V_i V i has a D D D -terminal embedding φ i \varphi_i φ i . Under the assumptions of Theorem 1.1 , Claim 2.5 gives such embeddings with D = 1 D=1 D = 1 .
To avoid the blow-up in the Lipschitz constant caused by simply concatenating the maps φ i \varphi_i φ i , we aim to ensure that f f f is C C C -Lipschitz for some absolute constant C C C .
Let N i N_i N i be the Voronoi cell of V i V_i V i , defined by
N i = { x : d ( x , V i ) < d ( x , V j ) for all j ≠ i } . N_i=\{x:d(x,V_i)<d(x,V_j)\text{ for all }j\ne i\}. N i = { x : d ( x , V i ) < d ( x , V j ) for all j = i } .
In addition to this Lipschitz condition, we seek a useful lower bound on the distance between the images of x ∈ N 1 ∪ ⋯ ∪ N r x\in N_1\cup\cdots\cup N_r x ∈ N 1 ∪ ⋯ ∪ N r and y ∈ X y\in X y ∈ X under f f f .
Construction idea The key idea is to use the map φ i \varphi_i φ i within each cell N i N_i N i , avoiding the k \sqrt{k} k Lipschitz constant blow-up. However, this idea needs a refinement, because as-is the map will have a huge Lipschitz constant near the “boundary” between different Voronoi cells. To avoid this, we gradually reduce the map to zero at the boundary. To formalize this idea, we introduce the following parameters.
Let
R i ( x ) = d ( x , V i ) R_i(x)=d(x,V_i) R i ( x ) = d ( x , V i )
R ^ i ( x ) = min j : j ≠ i d ( x , V j ) \widehat{R}_i(x)=\min_{j:j\ne i}d(x,V_j) R i ( x ) = j : j = i min d ( x , V j )
Δ i ( x ) = { R ^ i ( x ) − R i ( x ) 2 for x ∈ N i 0 , otherwise \Delta_i(x)=
\begin{cases}
\dfrac{\widehat{R}_i(x)-R_i(x)}{2} & \text{for }x\in N_i\\
0, & \text{otherwise}
\end{cases} Δ i ( x ) = ⎩ ⎨ ⎧ 2 R i ( x ) − R i ( x ) 0 , for x ∈ N i otherwise
We let R ^ i ( x ) = Δ i ( x ) = ∞ \widehat{R}_i(x)=\Delta_i(x)=\infty R i ( x ) = Δ i ( x ) = ∞ , if i = r = 1 i=r=1 i = r = 1 (there is only one set V 1 V_1 V 1 ). By definition, R i ( x ) < R ^ i ( x ) R_i(x)<\widehat{R}_i(x) R i ( x ) < R i ( x ) inside N i N_i N i and thus Δ i ( x ) ≥ 0 \Delta_i(x)\ge0 Δ i ( x ) ≥ 0 . Note that the distance to a fixed nonempty set is a 1-Lipschitz function and therefore R i R_i R i is 1-Lipschitz. If r > 1 r>1 r > 1 , then R ^ i \widehat{R}_i R i is also 1-Lipschitz and thus R ^ i ( x ) − R i ( x ) \widehat{R}_i(x)-R_i(x) R i ( x ) − R i ( x ) is 2-Lipschitz. Therefore, Δ i ( x ) = max ( R ^ i ( x ) − R i ( x ) , 0 ) 2 \Delta_i(x)=\dfrac{\max(\widehat{R}_i(x)-R_i(x),0)}{2} Δ i ( x ) = 2 max ( R i ( x ) − R i ( x ) , 0 ) is 1-Lipschitz. Informally, Δ i ( x ) \Delta_i(x) Δ i ( x ) serves as a proxy for the distance from x ∈ N i x\in N_i x ∈ N i to the boundary of N i N_i N i , since the open ball of radius Δ i ( x ) \Delta_i(x) Δ i ( x ) around x ∈ N i x\in N_i x ∈ N i lies entirely within N i N_i N i .
Figure 2 illustrates the distances to the sets and the open ball of radius Δ i ( x ) \Delta_i(x) Δ i ( x ) contained in N i N_i N i .
Figure 2. A schematic Euclidean illustration of the Voronoi cells (not to scale). For x ∈ N i x \in N_i x ∈ N i , R i ( x ) R_i(x) R i ( x ) is the distance to V i V_i V i , while R ^ i ( x ) \widehat{R}_i(x) R i ( x ) is the distance to the nearest set V j V_j V j other than V i V_i V i .
Now we are ready to define the desired embedding f f f of X X X :
f ( x ) = ⨁ i = 1 r ψ ( φ i ( x ) , Δ i ( x ) ) , f(x)=\bigoplus_{i=1}^{r}\psi(\varphi_i(x),\Delta_i(x)), f ( x ) = i = 1 ⨁ r ψ ( φ i ( x ) , Δ i ( x )) ,
where, for i = r = 1 i=r=1 i = r = 1 , we use the convention that ψ ( φ i ( x ) , ∞ ) = φ i ( x ) \psi(\varphi_i(x),\infty)=\varphi_i(x) ψ ( φ i ( x ) , ∞ ) = φ i ( x ) . Note that in this case, f ( x ) = φ 1 ( x ) f(x)=\varphi_1(x) f ( x ) = φ 1 ( x ) is 1-Lipschitz.
Theorem 2.6. The map f f f satisfies the following properties, with constants independent of D D D .
f f f is C C C -Lipschitz for some absolute constant C C C ;
for every i ∈ [ r ] i\in[r] i ∈ [ r ] , x ∈ N i x\in N_i x ∈ N i and y ∈ X y\in X y ∈ X ,
∥ f ( x ) − f ( y ) ∥ ≥ Ω ( min ( d ( x , y ) 2 D , R ^ i ( x ) ) − R i ( x ) ) . \lVert f(x)-f(y)\rVert\ge\Omega\left(\min\left(\frac{d(x,y)}{2D},\widehat{R}_i(x)\right)-R_i(x)\right). ∥ f ( x ) − f ( y )∥ ≥ Ω ( min ( 2 D d ( x , y ) , R i ( x ) ) − R i ( x ) ) .
Proof. First, we verify that f f f is C C C -Lipschitz. If r = 1 r=1 r = 1 , f f f is 1-Lipschitz, so let us assume that r ≥ 2 r \ge2 r ≥ 2 . Since φ i \varphi_i φ i is a D D D -terminal embedding for V i V_i V i , it is 1-Lipschitz. We have observed above that Δ i \Delta_i Δ i is also 1-Lipschitz. By Corollary 2.2, item 2, ∥ ψ ( u , a ) − ψ ( v , b ) ∥ = O ( ∥ u − v ∥ + ∣ a − b ∣ ) \|\psi(u,a)-\psi(v,b)\|=O(\|u-v\|+|a-b|) ∥ ψ ( u , a ) − ψ ( v , b ) ∥ = O ( ∥ u − v ∥ + ∣ a − b ∣ ) and thus ψ \psi ψ is C ′ C' C ′ -Lipschitz for some absolute constant C ′ C' C ′ . We conclude that ψ ( φ i ( x ) , Δ i ( x ) ) \psi(\varphi_i(x),\Delta_i(x)) ψ ( φ i ( x ) , Δ i ( x )) is C ′ ′ C'' C ′′ -Lipschitz. For every x , y ∈ X x,y \in X x , y ∈ X , at most one coordinate of f ( x ) f(x) f ( x ) is non-zero and at most one of the coordinates of f ( y ) f(y) f ( y ) is non-zero. Therefore, we have that ∥ f ( x ) − f ( y ) ∥ ≤ 2 C ′ ′ d ( x , y ) \|f(x)-f(y)\| \le\sqrt{2}C''d(x,y) ∥ f ( x ) − f ( y ) ∥ ≤ 2 C ′′ d ( x , y ) . That is, f f f is ( 2 C ′ ′ ) (\sqrt{2}C'') ( 2 C ′′ ) -Lipschitz.
Now we verify the second property. We first show that for every x , y ∈ X x,y \in X x , y ∈ X , we have
∥ φ i ( x ) − φ i ( y ) ∥ ≥ d ( x , y ) / D − 2 R i ( x ) . \|\varphi_i(x)-\varphi_i(y)\| \ge d(x,y)/D-2R_i(x). ∥ φ i ( x ) − φ i ( y ) ∥ ≥ d ( x , y ) / D − 2 R i ( x ) .
Assume first that V i V_i V i is compact. Let x ′ x' x ′ be a closest point to x x x in V i V_i V i . Then d ( x , x ′ ) = R i ( x ) d(x,x')=R_i(x) d ( x , x ′ ) = R i ( x ) . Applying the triangle inequality twice, we get
∥ φ i ( x ) − φ i ( y ) ∥ ≥ ∥ φ i ( x ′ ) − φ i ( y ) ∥ − ∥ φ i ( x ) − φ i ( x ′ ) ∥ ≥ d ( x ′ , y ) / D − d ( x , x ′ ) ≥ d ( x , y ) / D − 2 d ( x , x ′ ) , (2) \begin{aligned}
\|\varphi_i(x)-\varphi_i(y)\| &\ge\|\varphi_i(x')-\varphi_i(y)\|-\|\varphi_i(x)-\varphi_i(x')\| \\
&\ge d(x',y)/D-d(x,x') \ge d(x,y)/D-2d(x,x'), \tag*{(2)}
\end{aligned} ∥ φ i ( x ) − φ i ( y ) ∥ ≥ ∥ φ i ( x ′ ) − φ i ( y ) ∥ − ∥ φ i ( x ) − φ i ( x ′ ) ∥ ≥ d ( x ′ , y ) / D − d ( x , x ′ ) ≥ d ( x , y ) / D − 2 d ( x , x ′ ) , (2)
as required. Here, we used that φ i \varphi_i φ i is 1-Lipschitz and contracts distances from points in V i V_i V i , including x ′ x' x ′ , to all points in X X X by at most a factor D D D . If V i V_i V i is not compact, choose a sequence x ′ x' x ′ in V i V_i V i such that d ( x , x ′ ) → R i ( x ) d(x,x') \to R_i(x) d ( x , x ′ ) → R i ( x ) and apply (2) . We get the desired inequality in the limit.
Let x ∈ N i x \in N_i x ∈ N i and y ∈ X y \in X y ∈ X . If r ≥ 2 r \ge2 r ≥ 2 , we have
∥ f ( x ) − f ( y ) ∥ ≥ ∥ ψ ( φ i ( x ) , Δ i ( x ) ) − ψ ( φ i ( y ) , Δ i ( y ) ) ∥ ≥ Ω ( min ( ∥ φ i ( x ) − φ i ( y ) ∥ + ∣ Δ i ( x ) − Δ i ( y ) ∣ , Δ i ( x ) + Δ i ( y ) ) ) ≥ Ω ( min ( ∥ φ i ( x ) − φ i ( y ) ∥ , Δ i ( x ) ) ) ≥ Ω ( min ( d ( x , y ) / D − 2 R i ( x ) , R ^ i ( x ) − R i ( x ) ) ) ≥ Ω ( min ( d ( x , y ) / ( 2 D ) , R ^ i ( x ) ) − R i ( x ) ) . \begin{aligned}
\|f(x)-f(y)\| &\ge\|\psi(\varphi_i(x),\Delta_i(x))-\psi(\varphi_i(y),\Delta_i(y))\| \\
&\ge\Omega(\min(\|\varphi_i(x)-\varphi_i(y)\|+|\Delta_i(x)-\Delta_i(y)|,\Delta_i(x)+\Delta_i(y))) \\
&\ge\Omega(\min(\|\varphi_i(x)-\varphi_i(y)\|,\Delta_i(x))) \\
&\ge\Omega(\min(d(x,y)/D-2R_i(x),\widehat{R}_i(x)-R_i(x))) \\
&\ge\Omega(\min(d(x,y)/(2D),\widehat{R}_i(x))-R_i(x)).
\end{aligned} ∥ f ( x ) − f ( y ) ∥ ≥ ∥ ψ ( φ i ( x ) , Δ i ( x )) − ψ ( φ i ( y ) , Δ i ( y )) ∥ ≥ Ω ( min ( ∥ φ i ( x ) − φ i ( y ) ∥ + ∣ Δ i ( x ) − Δ i ( y ) ∣ , Δ i ( x ) + Δ i ( y ))) ≥ Ω ( min ( ∥ φ i ( x ) − φ i ( y ) ∥ , Δ i ( x ))) ≥ Ω ( min ( d ( x , y ) / D − 2 R i ( x ) , R i ( x ) − R i ( x ))) ≥ Ω ( min ( d ( x , y ) / ( 2 D ) , R i ( x )) − R i ( x )) .
For r = 1 r=1 r = 1 , we have f ( x ) = φ 1 ( x ) f(x)=\varphi_1(x) f ( x ) = φ 1 ( x ) and thus ∥ f ( x ) − f ( y ) ∥ = ∥ φ 1 ( x ) − φ 1 ( y ) ∥ ≥ d ( x , y ) / D − 2 R 1 ( x ) = 2 ( min ( d ( x , y ) / ( 2 D ) , ∞ ) − R 1 ( x ) ) \|f(x)-f(y)\|=\|\varphi_1(x)-\varphi_1(y)\| \ge d(x,y)/D-2R_1(x)=2(\min(d(x,y)/(2D),\infty)-R_1(x)) ∥ f ( x ) − f ( y ) ∥ = ∥ φ 1 ( x ) − φ 1 ( y ) ∥ ≥ d ( x , y ) / D − 2 R 1 ( x ) = 2 ( min ( d ( x , y ) / ( 2 D ) , ∞ ) − R 1 ( x )) .
Proof of Theorem 1.1
We prove the following statement, which gives Theorem 1.1 when D = 1 D=1 D = 1 .
Theorem 2.7. Let ( X , d ) (X,d) ( X , d ) be a metric space covered by nonempty sets U 1 , … , U k U_1,\ldots,U_k U 1 , … , U k . Suppose that every U i U_i U i has a D D D -terminal embedding into ℓ 2 \ell_2 ℓ 2 . Then X X X embeds into ℓ 2 \ell_2 ℓ 2 with distortion O ( D log ( k + 1 ) ) O(D\log(k+1)) O ( D log ( k + 1 )) .
Proof of Theorem 2.7 . We construct a random embedding F F F of X X X into ℓ 2 \ell_2 ℓ 2 . As the space of random vectors L 2 ( Ω , ℓ 2 ) L_2(\Omega,\ell_2) L 2 ( Ω , ℓ 2 ) is isometrically isomorphic to ℓ 2 \ell_2 ℓ 2 , F F F will yield the desired embedding of X X X into ℓ 2 \ell_2 ℓ 2 . Let ℓ = ⌊ log 2 k ⌋ \ell=\lfloor\log_2 k\rfloor ℓ = ⌊ log 2 k ⌋ . Define a random subset A ⊆ [ k ] A\subseteq[k] A ⊆ [ k ] as follows.
Choose s s s uniformly at random from { 0 , 1 , … , ℓ } \{0,1,\ldots,\ell\} { 0 , 1 , … , ℓ } .
For every x x x , include x x x in A A A independently with probability 1 / 2 s + 1 1/2^{s+1} 1/ 2 s + 1 .
Construct the embedding f A f_A f A for the sets U i U_i U i with i ∈ A i\in A i ∈ A as in the previous section. If A = ∅ A=\varnothing A = ∅ , define f A f_A f A to be the identically zero map. Let F = f A F=f_A F = f A .
We now analyze the distortion of F F F . By Theorem 2.6, item 1, all maps f A f_A f A are C C C -Lipschitz for some absolute constant C C C , so the map F F F is also C C C -Lipschitz. Now, fix two points x , y ∈ X x,y\in X x , y ∈ X . We have
∥ F ( x ) − F ( y ) ∥ L 2 = E [ ∥ F ( x ) − F ( y ) ∥ 2 ] ≥ E A ∥ f A ( x ) − f A ( y ) ∥ \lVert F(x)-F(y)\rVert_{L_2}=\sqrt{\mathbb{E}\left[\lVert F(x)-F(y)\rVert^2\right]}\geq\mathbb{E}_A\lVert f_A(x)-f_A(y)\rVert ∥ F ( x ) − F ( y ) ∥ L 2 = E [ ∥ F ( x ) − F ( y ) ∥ 2 ] ≥ E A ∥ f A ( x ) − f A ( y )∥
by Jensen’s inequality. Let g i = d ( x , U i ) g_i=d(x,U_i) g i = d ( x , U i ) for i ∈ [ k ] i\in[k] i ∈ [ k ] . Note x x x belongs to some U i U_i U i , and thus one of the g i g_i g i is 0 0 0 . Without loss of generality, we may assume that g 1 ≤ g 2 ≤ ⋯ ≤ g k g_1\leq g_2\leq\cdots\leq g_k g 1 ≤ g 2 ≤ ⋯ ≤ g k . Then g 1 = 0 g_1=0 g 1 = 0 . Further, extend the definition of g i g_i g i by letting g j = ∞ g_j=\infty g j = ∞ for j > k j>k j > k .
Let s 0 s_0 s 0 be the random number chosen in the first step of the construction of A A A . Let E s \mathcal{E}_s E s be the event that exactly one number in [ 2 s ] [2^s] [ 2 s ] is included in A A A and no numbers from [ 2 s + 1 , 2 s + 1 ] [2^s+1,2^{s+1}] [ 2 s + 1 , 2 s + 1 ] are included in A A A . Conditioned on s = s 0 s=s_0 s = s 0 , the probability of this event is at least some absolute constant. Figure 3 illustrates the event E s \mathcal{E}_s E s .
Figure 3. The event E s \mathcal{E}_s E s : exactly one index in the first block is selected, none in the second, and selections in the remaining block are unrestricted.
Assume E s \mathcal{E}_s E s occurs. Let i ∈ [ 2 s ] i\in[2^s] i ∈ [ 2 s ] be the smallest element in A A A . Then
R i ( x ) = d ( x , U i ) = g i ≤ g 2 s . R_i(x)=d(x,U_i)=g_i\leq g_{2^s}. R i ( x ) = d ( x , U i ) = g i ≤ g 2 s .
On the other hand, the second smallest element in A A A , if present, is at least 2 s + 1 + 1 2^{s+1}+1 2 s + 1 + 1 . Therefore,
R ^ i ( x ) ≥ g 2 s + 1 + 1 ≥ g 2 s + 1 . \widehat{R}_i(x)\geq g_{2^{s+1}+1}\geq g_{2^{s+1}}. R i ( x ) ≥ g 2 s + 1 + 1 ≥ g 2 s + 1 .
If A A A contains only one element, then the expression on the left-hand side is ∞ \infty ∞ and the inequality holds trivially. Finally, let t ∈ { 0 , 1 , … , ℓ } t\in\{0,1,\ldots,\ell\} t ∈ { 0 , 1 , … , ℓ } be the largest integer such that g 2 t ≤ d ( x , y ) / ( 2 D ) g_{2^t}\leq d(x,y)/(2D) g 2 t ≤ d ( x , y ) / ( 2 D ) . The number t t t is well-defined, since g 2 0 = g 1 = 0 ≤ d ( x , y ) / ( 2 D ) g_{2^0}=g_1=0\leq d(x,y)/(2D) g 2 0 = g 1 = 0 ≤ d ( x , y ) / ( 2 D ) and g 2 ℓ + 1 = ∞ g_{2^{\ell+1}}=\infty g 2 ℓ + 1 = ∞ . For s ≤ t s\leq t s ≤ t , if min ( d ( x , y ) / ( 2 D ) , g 2 s + 1 ) > g 2 s \min(d(x,y)/(2D),g_{2^{s+1}})>g_{2^s} min ( d ( x , y ) / ( 2 D ) , g 2 s + 1 ) > g 2 s , then R ^ i ( x ) > R i ( x ) \widehat{R}_i(x)>R_i(x) R i ( x ) > R i ( x ) and hence x ∈ N i x\in N_i x ∈ N i ; otherwise, the corresponding summand below is zero. Applying Theorem 2.6, item 2, and averaging over the scales, we obtain a telescoping sum:
∥ F ( x ) − F ( y ) ∥ L 2 ≥ E A ∥ f A ( x ) − f A ( y ) ∥ ≥ Ω ( 1 ) ℓ + 1 ∑ s = 0 ℓ E A [ ∥ f A ( x ) − f A ( y ) ∥ ∣ s 0 = s , E s ] ≥ Ω ( 1 ) ℓ + 1 ∑ s = 0 t ( min ( d ( x , y ) 2 D , g 2 s + 1 ) − g 2 s ) = Ω ( 1 ) ℓ + 1 ( ( g 2 − g 1 ) + ( g 4 − g 2 ) + ⋯ + ( g 2 t − g 2 t − 1 ) + ( d ( x , y ) 2 D − g 2 t ) ) = Ω ( d ( x , y ) ) D ( ℓ + 1 ) . \begin{aligned}
\lVert F(x)-F(y)\rVert_{L_2}\geq\mathbb{E}_A\lVert f_A(x)-f_A(y)\rVert\\
&\ge\frac{\Omega(1)}{\ell+1}\sum_{s=0}^{\ell}\mathbb{E}_{A}\left[\left\|f_A(x)-f_A(y)\right\| \mid s_0=s,\mathcal{E}_s\right] \\
&\ge\frac{\Omega(1)}{\ell+1}\sum_{s=0}^{t}\left(\min\left(\frac{d(x,y)}{2D},g_{2s+1}\right)-g_{2s}\right) \\
&= \frac{\Omega(1)}{\ell+1}\left((g_2-g_1)+(g_4-g_2)+\cdots+(g_{2t}-g_{2t-1})+\left(\frac{d(x,y)}{2D}-g_{2t}\right)\right) \\
&= \frac{\Omega(d(x,y))}{D(\ell+1)}.
\end{aligned} ∥ F ( x ) − F ( y ) ∥ L 2 ≥ E A ∥ f A ( x ) − f A ( y )∥ ≥ ℓ + 1 Ω ( 1 ) s = 0 ∑ ℓ E A [ ∥ f A ( x ) − f A ( y ) ∥ ∣ s 0 = s , E s ] ≥ ℓ + 1 Ω ( 1 ) s = 0 ∑ t ( min ( 2 D d ( x , y ) , g 2 s + 1 ) − g 2 s ) = ℓ + 1 Ω ( 1 ) ( ( g 2 − g 1 ) + ( g 4 − g 2 ) + ⋯ + ( g 2 t − g 2 t − 1 ) + ( 2 D d ( x , y ) − g 2 t ) ) = D ( ℓ + 1 ) Ω ( d ( x , y )) .
We conclude that F F F has distortion O ( D ( ℓ + 1 ) ) = O ( D log ( k + 1 ) ) O(D(\ell+1))=O(D\log(k+1)) O ( D ( ℓ + 1 )) = O ( D log ( k + 1 )) , as required.
Proof of Theorem 1.1 . By Claim 2.5 , every U i U_i U i has a 1-terminal embedding. Apply Theorem 2.7 with D = 1 D=1 D = 1 .
Relationship between Distortion and Lipschitz Extension Constants
In this section we establish Theorem 1.4 , which wants us to show e max ≤ c 2 ( X ) ≤ O ( e max log ( k + 1 ) ) e_{\max}\le c_2(X)\le O(e_{\max}\log(k+1)) e m a x ≤ c 2 ( X ) ≤ O ( e m a x log ( k + 1 )) , where e max = max i e ( U i ) e_{\max}=\max_i e(U_i) e m a x = max i e ( U i ) is the maximal Lipschitz extension constant for any U i U_i U i . We first show how to turn a Lipschitz extension into a terminal embedding. Then Theorem 2.7 gives an upper bound on the distortion. The lower bound relies only on Kirszbraun’s Theorem [14 ] .
Upper Bound via Terminal Embeddings
Lemma 3.1 (From Lipschitz extension to terminal embedding). Let V V V be a nonempty subset of a finite metric space ( X , d ) (X,d) ( X , d ) , and let D ≥ 1 D\ge1 D ≥ 1 . If a non-contracting embedding φ : V → ℓ 2 \varphi:V\to\ell_2 φ : V → ℓ 2 extends to a map φ ~ : X → ℓ 2 \widetilde{\varphi}:X\to\ell_2 φ : X → ℓ 2 with Lip ( φ ~ ) ≤ D \operatorname{Lip}(\widetilde{\varphi})\le D Lip ( φ ) ≤ D , then V V V has a 3 2 D 3\sqrt{2}D 3 2 D -terminal embedding.
Proof. Set R ( x ) = d ( x , V ) R(x)=d(x,V) R ( x ) = d ( x , V ) and
φ ^ ( x ) = 1 D φ ~ ( x ) ⊕ 2 R ( x ) . \widehat{\varphi}(x)=\frac{1}{D}\widetilde{\varphi}(x)\oplus_2 R(x). φ ( x ) = D 1 φ ( x ) ⊕ 2 R ( x ) .
Since both φ ~ / D \widetilde{\varphi}/D φ / D and R R R are 1-Lipschitz, Lip ( φ ^ ) ≤ 2 \operatorname{Lip}(\widehat{\varphi})\le\sqrt{2} Lip ( φ ) ≤ 2 . For u ∈ V u\in V u ∈ V and y ∈ X y\in X y ∈ X , choose v ∈ V v\in V v ∈ V closest to y y y . Then
d ( u , y ) ≤ d ( u , v ) + R ( y ) ≤ ∥ φ ~ ( u ) − φ ~ ( v ) ∥ + R ( y ) because φ is non-contracting on V ≤ ∥ φ ~ ( u ) − φ ~ ( y ) ∥ + ∥ φ ~ ( y ) − φ ~ ( v ) ∥ + R ( y ) ≤ ∥ φ ~ ( u ) − φ ~ ( y ) ∥ + ( D + 1 ) R ( y ) because Lip ( φ ~ ) ≤ D ≤ 3 D ∥ φ ^ ( u ) − φ ^ ( y ) ∥ because R ( u ) = 0. \begin{aligned}
d(u,y)&\le d(u,v)+R(y)\\
&\le\left\|\widetilde{\varphi}(u)-\widetilde{\varphi}(v)\right\|+R(y) && \text{because }\varphi\text{ is non-contracting on }V\\
&\le\left\|\widetilde{\varphi}(u)-\widetilde{\varphi}(y)\right\|+\left\|\widetilde{\varphi}(y)-\widetilde{\varphi}(v)\right\|+R(y)\\
&\le\left\|\widetilde{\varphi}(u)-\widetilde{\varphi}(y)\right\|+(D+1)R(y) && \text{because }\operatorname{Lip}(\widetilde{\varphi})\le D\\
&\le3D\left\|\widehat{\varphi}(u)-\widehat{\varphi}(y)\right\| && \text{because }R(u)=0.
\end{aligned} d ( u , y ) ≤ d ( u , v ) + R ( y ) ≤ ∥ φ ( u ) − φ ( v ) ∥ + R ( y ) ≤ ∥ φ ( u ) − φ ( y ) ∥ + ∥ φ ( y ) − φ ( v ) ∥ + R ( y ) ≤ ∥ φ ( u ) − φ ( y ) ∥ + ( D + 1 ) R ( y ) ≤ 3 D ∥ φ ( u ) − φ ( y ) ∥ because φ is non-contracting on V because Lip ( φ ) ≤ D because R ( u ) = 0.
Therefore, φ ^ / 2 \widehat{\varphi}/\sqrt{2} φ / 2 is 1-Lipschitz and contracts distances from V V V by at most a factor of 3 2 D 3\sqrt{2}D 3 2 D , as required.
Lemma 3.2 (Lipschitz extension constant upper bounds distortion). Let ( X , d ) (X,d) ( X , d ) be a finite metric space covered by nonempty sets U 1 , … , U k U_1,\ldots,U_k U 1 , … , U k , so that each U i U_i U i is Euclidean in X X X . Then c 2 ( X ) ≤ O ( e max log ( k + 1 ) ) c_2(X)\le O(e_{\max}\log(k+1)) c 2 ( X ) ≤ O ( e m a x log ( k + 1 )) .
Proof. Fix i ∈ [ k ] i \in[k] i ∈ [ k ] and an isometric embedding φ i : U i → ℓ 2 \varphi_i: U_i \to\ell_2 φ i : U i → ℓ 2 . It is 1-Lipschitz, so it extends to φ ~ i : X → ℓ 2 \widetilde{\varphi}_i: X \to\ell_2 φ i : X → ℓ 2 with Lip ( φ ~ i ) ≤ e ( U i ) ≤ e max \operatorname{Lip}(\widetilde{\varphi}_i) \le e(U_i) \le e_{\max} Lip ( φ i ) ≤ e ( U i ) ≤ e m a x , and Lemma 3.1 gives a 3 2 e max 3\sqrt{2}\,e_{\max} 3 2 e m a x -terminal embedding for U i U_i U i . Now Theorem 2.7 gives the claim.
Lower Bound
Lemma 3.3 (Lipschitz extension constant lower bounds distortion). Let ( X , d ) (X,d) ( X , d ) be a metric space that embeds into ℓ 2 \ell_2 ℓ 2 with distortion D D D , and let V ⊆ X V \subseteq X V ⊆ X be a nonempty subset. Then every Lipschitz map f : V → ℓ 2 f: V \to\ell_2 f : V → ℓ 2 extends to a map f ~ : X → ℓ 2 \widetilde{f}: X \to\ell_2 f : X → ℓ 2 with Lip ( f ~ ) ≤ D Lip ( f ) \operatorname{Lip}(\widetilde{f}) \le D\operatorname{Lip}(f) Lip ( f ) ≤ D Lip ( f ) . In particular, e ( V ) ≤ c 2 ( X ) e(V) \le c_2(X) e ( V ) ≤ c 2 ( X ) .
Proof. Let F : X → ℓ 2 F: X \to\ell_2 F : X → ℓ 2 be an embedding with distortion D D D , so that d ( x , y ) ≤ ∥ F ( x ) − F ( y ) ∥ ≤ D d ( x , y ) d(x,y) \le\lVert F(x)-F(y)\rVert\le Dd(x,y) d ( x , y ) ≤ ∥ F ( x ) − F ( y )∥ ≤ D d ( x , y ) for all x , y ∈ X x,y \in X x , y ∈ X . Since F F F does not contract distances, it is injective, and the map g = f ∘ F − 1 g=f\circ F^{-1} g = f ∘ F − 1 is well-defined on F ( V ) ⊆ ℓ 2 F(V) \subseteq\ell_2 F ( V ) ⊆ ℓ 2 . For u , v ∈ V u,v \in V u , v ∈ V ,
∥ g ( F ( u ) ) − g ( F ( v ) ) ∥ = ∥ f ( u ) − f ( v ) ∥ ≤ Lip ( f ) d ( u , v ) ≤ Lip ( f ) ∥ F ( u ) − F ( v ) ∥ , \lVert g(F(u))-g(F(v))\rVert= \lVert f(u)-f(v)\rVert\le\operatorname{Lip}(f)d(u,v) \le\operatorname{Lip}(f)\lVert F(u)-F(v)\rVert, ∥ g ( F ( u )) − g ( F ( v ))∥ = ∥ f ( u ) − f ( v )∥ ≤ Lip ( f ) d ( u , v ) ≤ Lip ( f ) ∥ F ( u ) − F ( v )∥ ,
so g g g is Lip ( f ) \operatorname{Lip}(f) Lip ( f ) -Lipschitz on F ( V ) F(V) F ( V ) . By Kirszbraun’s Theorem [14 ] , the map g g g extends to a map G : ℓ 2 → ℓ 2 G: \ell_2 \to\ell_2 G : ℓ 2 → ℓ 2 with Lip ( G ) ≤ Lip ( f ) \operatorname{Lip}(G) \le\operatorname{Lip}(f) Lip ( G ) ≤ Lip ( f ) . Let f ~ = G ∘ F \widetilde{f}=G\circ F f = G ∘ F . We show that this is a Lipschitz extension of f f f . For u ∈ V u \in V u ∈ V we have f ~ ( u ) = G ( F ( u ) ) = g ( F ( u ) ) = f ( u ) \widetilde{f}(u)=G(F(u))=g(F(u))=f(u) f ( u ) = G ( F ( u )) = g ( F ( u )) = f ( u ) , so, indeed, f ~ \widetilde{f} f extends f f f . For x , y ∈ X x,y \in X x , y ∈ X ,
∥ f ~ ( x ) − f ~ ( y ) ∥ ≤ Lip ( f ) ∥ F ( x ) − F ( y ) ∥ ≤ D Lip ( f ) d ( x , y ) , \lVert\widetilde{f}(x)-\widetilde{f}(y)\rVert\le\operatorname{Lip}(f)\lVert F(x)-F(y)\rVert\le D\operatorname{Lip}(f)d(x,y), ∥ f ( x ) − f ( y )∥ ≤ Lip ( f ) ∥ F ( x ) − F ( y )∥ ≤ D Lip ( f ) d ( x , y ) ,
which gives Lip ( f ~ ) ≤ D Lip ( f ) \operatorname{Lip}(\widetilde{f}) \le D\operatorname{Lip}(f) Lip ( f ) ≤ D Lip ( f ) . The map f f f was an arbitrary Lipschitz map for V V V , so D D D must be at least e ( V ) e(V) e ( V ) . Taking D = c 2 ( X ) D=c_2(X) D = c 2 ( X ) gives e ( V ) ≤ c 2 ( X ) e(V) \le c_2(X) e ( V ) ≤ c 2 ( X ) .
Proof of Theorem 1.4 . Discard empty sets among { U i } i = 1 k \{U_i\}_{i=1}^k { U i } i = 1 k . The lower bound is Lemma 3.3 applied to each U i U_i U i , and the upper bound is Lemma 3.2 .
Lipschitz Extension Constants of Euclidean and 1-Strongly Euclidean Sets
In this section we prove Lemma 1.5. Throughout, X X X is a finite pseudometric space, so distinct points may be at distance zero, and U 1 , … , U k U_1,\ldots,U_k U 1 , … , U k are nonempty sets covering X X X that are Euclidean, or 1-strongly Euclidean, in X X X . We fix A = U i A=U_i A = U i and a 1-Lipschitz map f : A → ℓ 2 f:A \to\ell_2 f : A → ℓ 2 , and we construct an extension F : X → ℓ 2 F:X \to\ell_2 F : X → ℓ 2 of f f f with Lip ( F ) = O ( k η ) \operatorname{Lip}(F)=O(k^\eta) Lip ( F ) = O ( k η ) , where η = η 0 = 1.1275 … \eta=\eta_0=1.1275\ldots η = η 0 = 1.1275 … for Euclidean sets and η = η 1 = 0.8941 … \eta=\eta_1=0.8941\ldots η = η 1 = 0.8941 … for 1-strongly Euclidean sets.
We bound Lip ( F ) \operatorname{Lip}(F) Lip ( F ) through a metric extension problem. Given a weighted graph G G G on X X X , removing the vertices outside A A A gives a graph G A G^A G A on A A A whose energy on f f f is the least energy of an extension of f f f (Lemma 4.4 ), and the best quality it guarantees as a metric sparsifier, over all graphs G G G , is the square of the Lipschitz extension constant taken over all pseudometrics that are Euclidean on every U j U_j U j (Theorem 4.2 ). We bound the quality by adding the sets U j U_j U j back one at a time (Lemmas 4.6 and 4.9), turn the energy bounds into extensions (Lemma 4.7 ), and average over the order in which the sets are added (Lemma 4.8 ).
Euclidean metric sparsifiers and Lipschitz extensions
We obtain the Lipschitz extension bounds needed above by studying metric extension costs and their sparsifiers. All sets in this section are finite, and we allow pseudometrics: distinct points may have distance zero.
In metric extension problems , we are given a weighted graph G G G on X X X with nonnegative weights g x y = g y x g_{xy}=g_{yx} g x y = g y x and g x x = 0 g_{xx}=0 g xx = 0 , a nonempty set of terminals A ⊆ X A\subseteq X A ⊆ X , and a metric ρ \rho ρ on A A A . We seek a metric d d d on X X X extending ρ \rho ρ that minimizes
∑ x < y g x y d ( x , y ) p \sum_{x<y} g_{xy}d(x,y)^p x < y ∑ g x y d ( x , y ) p
for a fixed exponent p ≥ 1 p\geq1 p ≥ 1 , subject to problem-specific constraints on d d d . Here x < y x<y x < y means that each unordered pair is included once.
For example, in the 0-extension problem, studied in [13 ] and [5 ] , we take p = 1 p=1 p = 1 and require
d ( x , y ) = ρ ( r ( x ) , r ( y ) ) , d(x,y)=\rho(r(x),r(y)), d ( x , y ) = ρ ( r ( x ) , r ( y )) ,
where r : X → A r:X\to A r : X → A fixes every terminal. Equivalently, every vertex must be at distance zero from some terminal. The standard LP relaxation drops this requirement and allows any pseudometric extension of ρ \rho ρ , with no further constraints.
In the metric labeling problem of [16 ] , we map the vertices into a prescribed ambient metric space ( Y , d Y ) (Y,d_Y) ( Y , d Y ) and take d ( x , y ) = d Y ( r ( x ) , r ( y ) ) d(x,y)=d_Y(r(x),r(y)) d ( x , y ) = d Y ( r ( x ) , r ( y )) , where r : X → Y r:X\to Y r : X → Y . In the extension setting, ( Y , d Y ) (Y,d_Y) ( Y , d Y ) contains ( A , ρ ) (A,\rho) ( A , ρ ) as a subspace and r r r fixes A A A , while the remaining vertices may be mapped anywhere in Y Y Y . With p = 1 p=1 p = 1 , the objective adds a nonnegative assignment cost c ( x , r ( x ) ) c(x,r(x)) c ( x , r ( x )) for each nonterminal vertex to the pairwise distance cost: ∑ x < y g x y d ( x , y ) + ∑ x ∈ X ∖ A c ( x , r ( x ) ) \sum_{x<y}g_{xy}d(x,y)+\sum_{x\in X\setminus A}c(x,r(x)) ∑ x < y g x y d ( x , y ) + ∑ x ∈ X ∖ A c ( x , r ( x )) .
For a problem-specific class D \mathcal{D} D of metrics on X X X , define the optimal extension cost by
min-ext D ( ρ , G ) = inf d ∈ D d ∣ A × A = ρ ∑ x < y g x y d ( x , y ) p . \operatorname{min-ext}_{\mathcal{D}}(\rho,G)=\inf_{\substack{d\in\mathcal{D}\\d|_{A\times A}=\rho}}\sum_{x<y}g_{xy}d(x,y)^p. min-ext D ( ρ , G ) = d ∈ D d ∣ A × A = ρ inf x < y ∑ g x y d ( x , y ) p .
A metric sparsifier lets us approximate this cost using only the distances on A A A , without computing an extension to X X X . Following the formulation from [20 ] for p = 1 p=1 p = 1 , we call a weighted graph H H H on A A A , with nonnegative weights h a b h_{ab} h ab , a metric sparsifier of quality Q ≥ 1 Q\geq1 Q ≥ 1 if
1 Q ∑ a < b h a b ρ ( a , b ) p ≤ min-ext D ( ρ , G ) ≤ ∑ a < b h a b ρ ( a , b ) p \frac{1}{Q}\sum_{a<b}h_{ab}\rho(a,b)^p\leq\operatorname{min-ext}_{\mathcal{D}}(\rho,G)\leq\sum_{a<b}h_{ab}\rho(a,b)^p Q 1 a < b ∑ h ab ρ ( a , b ) p ≤ min-ext D ( ρ , G ) ≤ a < b ∑ h ab ρ ( a , b ) p
for every terminal metric ρ \rho ρ admitting an extension in D \mathcal{D} D . The same graph H H H provides this approximation for every such ρ \rho ρ .
For our purposes, we fix a partition of X X X into nonempty sets U 1 , … , U k U_1,\ldots,U_k U 1 , … , U k , take A = ⋃ i U i A=\bigcup_i U_i A = ⋃ i U i and p = 2 p=2 p = 2 , and let D \mathcal{D} D be the class of pseudometrics d d d on X X X that are Euclidean on every U j U_j U j . For any pseudometric d d d on X X X , define its energy by
E G ( d 2 ) = ∑ x < y g x y d ( x , y ) 2 . \mathcal{E}_G(d^2)=\sum_{x<y}g_{xy}d(x,y)^2. E G ( d 2 ) = x < y ∑ g x y d ( x , y ) 2 .
This definition does not require a Euclidean representation of d d d . For a map f : X → ℓ 2 f:X\to\ell_2 f : X → ℓ 2 , we also define
E G ( f ) = ∑ x < y g x y ∥ f ( x ) − f ( y ) ∥ 2 . \mathcal{E}_G(f)=\sum_{x<y}g_{xy}\lVert f(x)-f(y)\rVert^2. E G ( f ) = x < y ∑ g x y ∥ f ( x ) − f ( y ) ∥ 2 .
This is the energy of the pseudometric induced by f f f . Thus, E G ( f ) = E G ( d 2 ) \mathcal{E}_G(f)=\mathcal{E}_G(d^2) E G ( f ) = E G ( d 2 ) when d ( x , y ) = ∥ f ( x ) − f ( y ) ∥ d(x,y)=\lVert f(x)-f(y)\rVert d ( x , y ) = ∥ f ( x ) − f ( y )∥ on all of X X X . A general d ∈ D d\in\mathcal{D} d ∈ D need not admit such a map: only its restrictions to the individual sets U j U_j U j are required to be Euclidean. For real-valued f f f , E G ( f ) \mathcal{E}_G(f) E G ( f ) is the energy dissipated per unit time by the electrical network with conductances g x y g_{xy} g x y and vertex potentials f ( x ) f(x) f ( x ) . For ℓ 2 \ell_2 ℓ 2 -valued f f f , it is the sum of the coordinate energies.
Specializing the definitions above, for a Euclidean metric ρ \rho ρ on A A A let
min-ext D ( ρ , G ) = inf d ∈ D d ∣ A × A = ρ E G ( d 2 ) . (3) \operatorname{min-ext}_{\mathcal{D}}(\rho,G)=\inf_{\substack{d\in\mathcal{D}\\d|_{A\times A}=\rho}}\mathcal{E}_G(d^2). \tag*{(3)} min-ext D ( ρ , G ) = d ∈ D d ∣ A × A = ρ inf E G ( d 2 ) . (3)
Definition 4.1. A weighted graph H H H on A A A , with nonnegative weights, is a Euclidean metric sparsifier of G G G for D \mathcal{D} D with quality Q ≥ 1 Q \ge1 Q ≥ 1 if, for every Euclidean ρ \rho ρ on A A A ,
1 Q E H ( ρ 2 ) ≤ min - ext D ( ρ , G ) ≤ E H ( ρ 2 ) . (4) \frac{1}{Q}\mathcal{E}_H(\rho^2) \le\min\text{-}\operatorname{ext}_{\mathcal{D}}(\rho,G) \le\mathcal{E}_H(\rho^2). \tag*{(4)} Q 1 E H ( ρ 2 ) ≤ min - ext D ( ρ , G ) ≤ E H ( ρ 2 ) . (4)
The following theorem makes the connection with Lipschitz extension constants precise. Let Q ( D , A ) Q(\mathcal{D}, A) Q ( D , A ) be the least quality guaranteed for every weighted graph G G G on X X X . Let e ( D , A ) e(\mathcal{D}, A) e ( D , A ) be the least Λ ≥ 1 \Lambda\ge1 Λ ≥ 1 such that, for every d ∈ D d \in\mathcal{D} d ∈ D , every Lipschitz map f : ( A , d ) → ℓ 2 f : (A,d) \to\ell_2 f : ( A , d ) → ℓ 2 extends to F : ( X , d ) → ℓ 2 F : (X,d) \to\ell_2 F : ( X , d ) → ℓ 2 with Lip ( F ) ≤ Λ Lip ( f ) \operatorname{Lip}(F) \le\Lambda\operatorname{Lip}(f) Lip ( F ) ≤ Λ Lip ( f ) . The graph and the metric may vary, while the class D \mathcal{D} D and the distinguished set A A A are fixed. Let Q k Q_k Q k and e k e_k e k be the suprema of these constants over all finite partitions into at most k k k sets and all choices of A A A . For the pseudometric d d d of X X X itself, which belongs to D \mathcal{D} D , the constant e ( A ) e(A) e ( A ) of § 1 is the value at d d d of the supremum defining e ( D , A ) e(\mathcal{D}, A) e ( D , A ) , so
e ( A ) ≤ e ( D , A ) ≤ e k . (5) e(A) \le e(\mathcal{D}, A) \le e_k. \tag*{(5)} e ( A ) ≤ e ( D , A ) ≤ e k . (5)
Hence, bounding e k e_k e k allows us to bound e ( A ) e(A) e ( A ) . We have the following result connecting the quality Q k Q_k Q k and the supremum e k e_k e k . The proof is deferred to Section A.
Theorem 4.2. For every partition U 1 , … , U k U_1,\ldots,U_k U 1 , … , U k of X X X , with D \mathcal{D} D as above, and every A = U i A=U_i A = U i , we have Q ( D , A ) = e ( D , A ) 2 Q(\mathcal{D},A)=e(\mathcal{D},A)^2 Q ( D , A ) = e ( D , A ) 2 . As a consequence, Q k = e k 2 Q_k=e_k^2 Q k = e k 2 .
Remark 4.3. The same connection between sparsifier quality and Lipschitz extension constants holds for other metric extension problems. For example, [20 Theorems 5.1–5.2] showed that the optimal qualities of cut and flow sparsifiers equal the corresponding Lipschitz extension constants between suitable Banach spaces. We omit the details here.
The canonical Euclidean metric sparsifier
We first consider the extension problem in which all of X X X must be Euclidean. Fix a nonempty set A ⊆ X A \subseteq X A ⊆ X . We define G A G^A G A to be the graph on A A A with nonnegative weights whose energy gives the minimum Euclidean extension cost:
E G A ( ρ 2 ) = min d ∣ A × A = ρ d Euclidean on X E G ( d 2 ) \mathcal{E}_{G^A}(\rho^2)=\min_{\substack{d|_{A\times A}=\rho\\ d\ \text{Euclidean on }X}}\mathcal{E}_G(d^2) E G A ( ρ 2 ) = d ∣ A × A = ρ d Euclidean on X min E G ( d 2 )
for every Euclidean pseudometric ρ \rho ρ on A A A . Only the distances on A A A are fixed. All other distances can be chosen to minimize the cost. We will show that one graph G A G^A G A works for every ρ \rho ρ , and that its weights are uniquely determined. The sets U j U_j U j play no role here. We return to them when comparing this minimum with the cost of extensions in D \mathcal{D} D .
The weights of G A G^A G A can be found by removing one vertex at a time. To see how, suppose we first remove a nonisolated vertex v v v with neighbors x 1 , … , x t x_1,\ldots,x_t x 1 , … , x t . Fix a Euclidean pseudometric d d d on X ∖ { v } X \setminus\{v\} X ∖ { v } . The distances among the neighbors can be written as d ( x a , x b ) = ∥ u a − u b ∥ d(x_a,x_b)=\lVert u_a-u_b\rVert d ( x a , x b ) = ∥ u a − u b ∥ for suitable vectors u 1 , … , u t u_1,\ldots,u_t u 1 , … , u t . Choosing a position z z z for v v v adds ∑ r = 1 t g v x r ∥ z − u r ∥ 2 \sum_{r=1}^{t}g_{vx_r}\lVert z-u_r\rVert^2 ∑ r = 1 t g v x r ∥ z − u r ∥ 2 to the energy. Recall that this is the objective for a single cluster in the weighted k k k -means problem. The optimal choice of the center is simply the weighted mean
u ˉ = 1 δ v ∑ r = 1 t g v x r u r , where δ v = ∑ r = 1 t g v x r > 0. \bar{u}=\frac{1}{\delta_v}\sum_{r=1}^{t}g_{vx_r}u_r,\qquad\text{where }\delta_v=\sum_{r=1}^{t}g_{vx_r}>0. u ˉ = δ v 1 r = 1 ∑ t g v x r u r , where δ v = r = 1 ∑ t g v x r > 0.
Using the standard alternative formula for the cost of a cluster under the k k k -means objective, we express the contribution of v v v to the energy without referring to z z z explicitly as
1 δ v ∑ 1 ≤ r < s ≤ t g v x r g v x s ∥ u r − u s ∥ 2 = 1 δ v ∑ 1 ≤ r < s ≤ t g v x r g v x s d ( x r , x s ) 2 . \frac{1}{\delta_v}\sum_{1\le r<s\le t}g_{v x_r}g_{v x_s}\lVert u_r-u_s\rVert^2=\frac{1}{\delta_v}\sum_{1\le r<s\le t}g_{v x_r}g_{v x_s}d(x_r,x_s)^2. δ v 1 1 ≤ r < s ≤ t ∑ g v x r g v x s ∥ u r − u s ∥ 2 = δ v 1 1 ≤ r < s ≤ t ∑ g v x r g v x s d ( x r , x s ) 2 .
To capture the contribution of v v v in G ′ = G X ∖ { v } G'=G^{X\setminus\{v\}} G ′ = G X ∖ { v } , we add g x v g y v / δ v g_{xv}g_{yv}/\delta_v g xv g y v / δ v to each weight g x y g_{xy} g x y :
G x y ′ = g x y + g x v g y v δ v , G'_{xy}=g_{xy}+\frac{g_{xv}g_{yv}}{\delta_v}, G x y ′ = g x y + δ v g xv g y v ,
where x , y ∈ X ∖ { v } x,y\in X\setminus\{v\} x , y ∈ X ∖ { v } , x ≠ y x\ne y x = y . The energy of G ′ G' G ′ therefore includes the cost of the best position for v v v , even though v v v is no longer in the graph. We repeat this step until only A A A remains, using the current weights and the current value of δ v \delta_v δ v at each step. An isolated vertex is simply deleted. This gives the graph G A G^A G A described above, as we now prove.
Lemma 4.4. For every weighted graph G G G on X X X and nonempty A ⊆ X A\subseteq X A ⊆ X , the graph G A G^A G A defined above exists and is unique. The update rule above constructs it by removing the vertices of X ∖ A X\setminus A X ∖ A in any order. For every map f : A → ℓ 2 f:A\to\ell_2 f : A → ℓ 2 , it satisfies
E G A ( f ) = min F : X → ℓ 2 F ∣ A = f E G ( F ) . (6) \mathcal{E}_{G^A}(f)=\min_{\substack{F:X\to\ell_2\\F|_A=f}}\mathcal{E}_G(F). \tag*{(6)} E G A ( f ) = F : X → ℓ 2 F ∣ A = f min E G ( F ) . (6) Moreover, for A ⊆ A ′ ⊆ X A\subseteq A'\subseteq X A ⊆ A ′ ⊆ X we have:
( G A ′ ) A = G A . (7) (G^{A'})^A=G^A. \tag*{(7)} ( G A ′ ) A = G A . (7)
For A = U i A=U_i A = U i , allowing all extensions in D \mathcal{D} D gives min - ext D ( ρ , G ) ≤ E G A ( ρ 2 ) \min\text{-}\operatorname{ext}_{\mathcal{D}}(\rho,G)\leq\mathcal{E}_{G^A}(\rho^2) min - ext D ( ρ , G ) ≤ E G A ( ρ 2 ) for every Euclidean ρ \rho ρ on A A A , which is the right inequality in (4) .
Proof. We begin by turning the weighted k k k -means discussion into a proof that the update rule gives the graph we want. Fix positions f ( a ) ∈ ℓ 2 f(a)\in\ell_2 f ( a ) ∈ ℓ 2 for the vertices a ∈ A a\in A a ∈ A , and remove the other vertices one at a time. With the positions of the remaining vertices fixed, the calculation above replaces the cost of the best position for v v v by the added edge weights. Thus, after each removal, the current graph records the minimum energy over the positions of all vertices removed so far. When only A A A remains, the resulting graph H H H satisfies
E H ( f ) = min F : X → ℓ 2 F ∣ A = f E G ( F ) . \mathcal{E}_H(f)=\min_{\substack{F:X\to\ell_2\\F|_A=f}}\mathcal{E}_G(F). E H ( f ) = F : X → ℓ 2 F ∣ A = f min E G ( F ) .
We can obtain a minimizing extension by restoring the removed vertices in reverse order, placing each at the weighted mean of its neighbors with the weights used when it was removed. An isolated vertex can be placed anywhere, since its position does not affect the energy.
Now let ρ \rho ρ be a Euclidean pseudometric on A A A , and choose f f f with ρ ( a , b ) = ∥ f ( a ) − f ( b ) ∥ \rho(a,b)=\lVert f(a)-f(b)\rVert ρ ( a , b ) = ∥ f ( a ) − f ( b )∥ . Every Euclidean extension of ρ \rho ρ can be represented by vectors that agree with f f f on A A A , so the minimum above is exactly the minimum Euclidean extension cost of ρ \rho ρ . Thus, H H H has the defining property of G A G^A G A , and (6) holds. Its weights are nonnegative, and the same graph works for every ρ \rho ρ , since the updates depend only on the graph weights.
These energies also determine the graph uniquely. Indeed, for real-valued f f f , the coefficient of f ( x ) f ( y ) f(x)f(y) f ( x ) f ( y ) in E H ( f ) \mathcal{E}_H(f) E H ( f ) is − 2 h x y -2h_{xy} − 2 h x y for x ≠ y x\ne y x = y . Hence, every removal order gives the same graph G A G^A G A . In particular, we may first remove all vertices outside A ′ A' A ′ and then those in A ′ ∖ A A'\setminus A A ′ ∖ A . This proves ( G A ′ ) A = G A (G^{A'})^A=G^A ( G A ′ ) A = G A .
Finally, when A = U i A=U_i A = U i , every globally Euclidean extension of ρ \rho ρ belongs to D \mathcal{D} D . Allowing all extensions in D \mathcal{D} D can only lower the minimum cost, so min - ext D ( ρ , G ) ≤ E G A ( ρ 2 ) \min\text{-}\operatorname{ext}_{\mathcal{D}}(\rho,G)\leq\mathcal{E}_{G^A}(\rho^2) min - ext D ( ρ , G ) ≤ E G A ( ρ 2 ) , as claimed. ▫
It remains to compare E G A ( ρ 2 ) \mathcal{E}_{G_A}(\rho^2) E G A ( ρ 2 ) with the energy of extensions of ρ \rho ρ that are Euclidean only on the individual sets U j U_j U j . We eliminate the sets U j U_j U j with j ≠ i j \ne i j = i one at a time, and first bound the energy added when a single Euclidean set is eliminated. We use G [ B ] G[B] G [ B ] to denote the graph induced by B B B , with the original weights.
A matrix estimate
The next claim helps bound the cost of the new edges created when we remove a Euclidean set B B B . It controls weighted averages of squared distances within B B B using the energy on its internal edges. In the application, the vertex weight s b s_b s b will be the total weight of edges from b ∈ B b \in B b ∈ B to the vertices we keep. We let L L L denote the non-normalized Laplacian of a weighted graph with edge weights w b b ′ w_{bb'} w b b ′ . For vectors z b z_b z b in a Hilbert space, one for each vertex, we have
∑ b ⟨ z b , ( L z ) b ⟩ = ∑ b < b ′ w b b ′ ∥ z b − z b ′ ∥ 2 . \sum_b \langle z_b,(Lz)_b\rangle= \sum_{b<b'} w_{bb'}\lVert z_b-z_{b'}\rVert^2. b ∑ ⟨ z b , ( L z ) b ⟩ = b < b ′ ∑ w b b ′ ∥ z b − z b ′ ∥ 2 .
Claim 4.5. Let L L L be the Laplacian of a weighted graph on B B B . Assign a positive weight s b s_b s b to each vertex b ∈ B b \in B b ∈ B , and let C = diag ( s b ) C=\operatorname{diag}(s_b) C = diag ( s b ) , K = C − 1 L K=C^{-1}L K = C − 1 L , and P ( t ) = e − t K P(t)=e^{-tK} P ( t ) = e − t K (t ≥ 0 t \ge0 t ≥ 0 ). Each P ( t ) P(t) P ( t ) has nonnegative entries and satisfies P ( t ) 1 = 1 P(t)\mathbf{1}=\mathbf{1} P ( t ) 1 = 1 , C P ( t ) = P ( t ) T C CP(t)=P(t)^{\mathsf{T}}C C P ( t ) = P ( t ) T C and P ( t ) T s = s P(t)^{\mathsf{T}}s=s P ( t ) T s = s , where 1 \mathbf{1} 1 is the all-ones vector and s = ( s b ) b ∈ B s=(s_b)_{b\in B} s = ( s b ) b ∈ B is the vector of vertex weights. The first two identities say that the rows of P ( t ) P(t) P ( t ) sum to one and that C P ( t ) CP(t) C P ( t ) is symmetric. For vectors z b z_b z b in a Hilbert space,
∑ b , b ′ s b P b b ′ ( t ) ∥ z b − z b ′ ∥ 2 ≤ 2 t ∑ b < b ′ w b b ′ ∥ z b − z b ′ ∥ 2 . (8) \sum_{b,b'}s_bP_{bb'}(t)\lVert z_b-z_{b'}\rVert^2 \le2t\sum_{b<b'}w_{bb'}\lVert z_b-z_{b'}\rVert^2. \tag*{(8)} b , b ′ ∑ s b P b b ′ ( t ) ∥ z b − z b ′ ∥ 2 ≤ 2 t b < b ′ ∑ w b b ′ ∥ z b − z b ′ ∥ 2 . (8)
Proof. The nonnegative entries and the identity P ( t ) 1 = 1 P(t)\mathbf{1}=\mathbf{1} P ( t ) 1 = 1 follow by writing P ( t ) P(t) P ( t ) as the limit of ( Id − t K / N ) N (\operatorname{Id}-tK/N)^N ( Id − t K / N ) N : for sufficiently large N N N , the matrix Id − t K / N \operatorname{Id}-tK/N Id − t K / N has nonnegative entries and fixes 1 \mathbf{1} 1 . With respect to the inner product
⟨ u , v ⟩ C = ∑ b s b ⟨ u b , v b ⟩ , \langle u,v\rangle_C=\sum_b s_b\langle u_b,v_b\rangle, ⟨ u , v ⟩ C = b ∑ s b ⟨ u b , v b ⟩ ,
the matrix K K K is self-adjoint and nonnegative, since C K = L CK=L C K = L . Thus, P ( t ) P(t) P ( t ) is also self-adjoint, so C P ( t ) = P ( t ) T C CP(t)=P(t)^{\mathsf{T}}C C P ( t ) = P ( t ) T C . Since s = C 1 s=C\mathbf{1} s = C 1 , we also get
P ( t ) T s = P ( t ) T C 1 = C P ( t ) 1 = C 1 = s . P(t)^{\mathsf{T}}s=P(t)^{\mathsf{T}}C\mathbf{1}=CP(t)\mathbf{1}=C\mathbf{1}=s. P ( t ) T s = P ( t ) T C 1 = C P ( t ) 1 = C 1 = s .
Expanding the squared distances now gives
∑ b , b ′ s b P b b ′ ( t ) ∥ z b − z b ′ ∥ 2 = 2 ⟨ z , ( Id − P ( t ) ) z ⟩ C \sum_{b,b'}s_bP_{bb'}(t)\lVert z_b-z_{b'}\rVert^2=2\langle z,(\operatorname{Id}-P(t))z\rangle_C b , b ′ ∑ s b P b b ′ ( t ) ∥ z b − z b ′ ∥ 2 = 2 ⟨ z , ( Id − P ( t )) z ⟩ C
≤ 2 t ⟨ z , K z ⟩ C = 2 t ∑ b < b ′ w b b ′ ∥ z b − z b ′ ∥ 2 . \le2t\langle z,Kz\rangle_C=2t\sum_{b<b'}w_{bb'}\lVert z_b-z_{b'}\rVert^2. ≤ 2 t ⟨ z , K z ⟩ C = 2 t b < b ′ ∑ w b b ′ ∥ z b − z b ′ ∥ 2 .
The inequality follows by diagonalizing K K K and applying 1 − e − t λ ≤ t λ 1-e^{-t\lambda}\le t\lambda 1 − e − t λ ≤ t λ to each eigenvalue λ ≥ 0 \lambda\ge0 λ ≥ 0 . □ \square □
Eliminating one Euclidean set
Lemma 4.6. Let Y = A ⊔ B Y=A\sqcup B Y = A ⊔ B be a finite pseudometric space, where A ≠ ∅ A\ne\varnothing A = ∅ and B B B is Euclidean. For a weighted graph G G G on Y Y Y , let I I I be the energy on edges within B B B and J J J the energy on edges between A A A and B B B , so I = E G [ B ] ( d 2 ) I=\mathcal{E}_{G[B]}(d^2) I = E G [ B ] ( d 2 ) and J = ∑ a ∈ A , b ∈ B g a b d ( a , b ) 2 J=\sum_{a\in A,b\in B}g_{ab}d(a,b)^2 J = ∑ a ∈ A , b ∈ B g ab d ( a , b ) 2 . Then
E G A ( d 2 ) ≤ E G ( d 2 ) + J + 2 π I J . (9) \mathcal{E}_{G_A}(d^2)\le\mathcal{E}_G(d^2)+J+\sqrt{2\pi IJ}. \tag*{(9)} E G A ( d 2 ) ≤ E G ( d 2 ) + J + 2 π I J . (9)
Proof. The idea is to bound the distance across each new edge on A A A by three distances: from A A A into B B B , across B B B , and back to A A A . The two outside parts will be controlled by J J J , and Claim 4.5 will control the middle part using I I I .
For b ∈ B b \in B b ∈ B , take s b = ∑ a ∈ A g a b s_b = \sum_{a\in A} g_{ab} s b = ∑ a ∈ A g ab . We first assume that B B B is nonempty and s b > 0 s_b > 0 s b > 0 for every b ∈ B b \in B b ∈ B , so each vertex of B B B has an edge to A A A . We remove this assumption at the end. Since B B B is Euclidean, choose vectors z b z_b z b with d ( b , b ′ ) = ∥ z b − z b ′ ∥ d(b,b') = \lVert z_b-z_{b'}\rVert d ( b , b ′ ) = ∥ z b − z b ′ ∥ .
We need to account for the weights added on A A A when B B B is removed. To compute them, we give each edge from A A A to B B B its own copy of its endpoint in A A A . Each copy then has just one neighbor in B B B , which makes the new weights easy to compute. More precisely, for each positive-weight edge with endpoints a p ∈ A a_p \in A a p ∈ A and b p ∈ B b_p \in B b p ∈ B , create a vertex p p p and join it to b p b_p b p with the original edge weight c p c_p c p . Our auxiliary graph consists of G [ B ] G[B] G [ B ] and these new edges. Remove B B B from this graph, and call the resulting weights between the copies ω p q \omega_{pq} ω pq .
To recover the weights added on A A A , give all copies of each vertex the same real value and use (6) . The weight added between two vertices of A A A is the sum of ω p q \omega_{pq} ω pq over their copies. Consequently,
E G A ( d 2 ) − E G [ A ] ( d 2 ) = ∑ p < q ω p q d ( a p , a q ) 2 . (10) \mathcal{E}_{G_A}(d^2)-\mathcal{E}_{G[A]}(d^2)=\sum_{p<q}\omega_{pq}d(a_p,a_q)^2. \tag*{(10)} E G A ( d 2 ) − E G [ A ] ( d 2 ) = p < q ∑ ω pq d ( a p , a q ) 2 . (10)
This is an identity between edge weights, so it holds even though A A A need not be Euclidean.
Let L L L be the Laplacian of G [ B ] G[B] G [ B ] , and let C = diag ( s b ) C=\operatorname{diag}(s_b) C = diag ( s b ) and M = L + C M=L+C M = L + C . Since all s b s_b s b are positive, M M M is positive definite. Minimizing over the values at B B B in the auxiliary graph gives for p ≠ q p\ne q p = q
ω p q = c p c q ( M − 1 ) b p b q . (11) \omega_{pq}=c_pc_q(M^{-1})_{b_pb_q}. \tag*{(11)} ω pq = c p c q ( M − 1 ) b p b q . (11)
Indeed, if the copies have real values η p \eta_p η p and the vertices of B B B have real values ξ b \xi_b ξ b , the energy is
ξ T M ξ − 2 ∑ p c p η p ξ b p + ∑ p c p η p 2 . \xi^{T}M\xi-2\sum_p c_p\eta_p\xi_{b_p}+\sum_p c_p\eta_p^2. ξ T M ξ − 2 p ∑ c p η p ξ b p + p ∑ c p η p 2 .
Completing the square and comparing the coefficients of η p η q \eta_p\eta_q η p η q gives (11) .
We next express these weights using the matrices P ( t ) P(t) P ( t ) from Claim 4.5 , applied to G [ B ] G[B] G [ B ] with vertex weights s b s_b s b . With K = C − 1 L K=C^{-1}L K = C − 1 L and P ( t ) = e − t K P(t)=e^{-tK} P ( t ) = e − t K , we have
M − 1 = ( Id + K ) − 1 C − 1 = ∫ 0 ∞ e − t P ( t ) C − 1 d t . M^{-1}=(\operatorname{Id}+K)^{-1}C^{-1}=\int_0^\infty e^{-t}P(t)C^{-1}\,dt. M − 1 = ( Id + K ) − 1 C − 1 = ∫ 0 ∞ e − t P ( t ) C − 1 d t .
Thus, (11) gives ω p q = ∫ 0 ∞ e − t μ p q ( t ) d t \omega_{pq}=\int_0^\infty e^{-t}\mu_{pq}(t)\,dt ω pq = ∫ 0 ∞ e − t μ pq ( t ) d t for p ≠ q p\ne q p = q , where
μ p q ( t ) = c p P b p b q ( t ) c q s b q . \mu_{pq}(t)=c_pP_{b_pb_q}(t)\frac{c_q}{s_{b_q}}. μ pq ( t ) = c p P b p b q ( t ) s b q c q .
These weights are nonnegative and symmetric, and satisfy
∑ q μ p q ( t ) = c p , ∑ p μ p q ( t ) = c q . (12) \sum_q\mu_{pq}(t)=c_p,\qquad\sum_p\mu_{pq}(t)=c_q. \tag*{(12)} q ∑ μ pq ( t ) = c p , p ∑ μ pq ( t ) = c q . (12)
For the first equality, group the copies q q q by their neighbor b q b_q b q , use ∑ q : b q = b ′ c q = s b ′ \sum_{q:b_q=b'}c_q=s_{b'} ∑ q : b q = b ′ c q = s b ′ , and then use the row sums of P ( t ) P(t) P ( t ) . The second follows by symmetry, which comes from s b P b b ′ ( t ) = s b ′ P b ′ b ( t ) s_bP_{bb'}(t)=s_{b'}P_{b'b}(t) s b P b b ′ ( t ) = s b ′ P b ′ b ( t ) .
For each t ≥ 0 t\ge0 t ≥ 0 , let
S t = ∑ p , q μ p q ( t ) d ( a p , a q ) 2 . S_t=\sum_{p,q}\mu_{pq}(t)d(a_p,a_q)^2. S t = p , q ∑ μ pq ( t ) d ( a p , a q ) 2 .
By symmetry, we may replace the sum over p < q p<q p < q in (10) by half the sum over all p , q p,q p , q . The terms with p = q p=q p = q contribute zero. Hence,
E G A ( d 2 ) − E G [ A ] ( d 2 ) = 1 2 ∫ 0 ∞ e − t S t d t . (13) \mathcal{E}_{G^A}(d^2)-\mathcal{E}_{G[A]}(d^2)=\frac{1}{2}\int_0^\infty e^{-t}S_t\,\mathrm{d}t. \tag*{(13)} E G A ( d 2 ) − E G [ A ] ( d 2 ) = 2 1 ∫ 0 ∞ e − t S t d t . (13)
We can now use the three distances along the route a p → b p → b q → a q a_p\to b_p\to b_q\to a_q a p → b p → b q → a q . Let u p = d ( a p , b p ) u_p=d(a_p,b_p) u p = d ( a p , b p ) . The triangle inequality gives
d ( a p , a q ) ≤ u p + ∥ z b p − z b q ∥ + u q . d(a_p,a_q)\leq u_p+\lVert z_{b_p}-z_{b_q}\rVert+u_q. d ( a p , a q ) ≤ u p + ∥ z b p − z b q ∥ + u q .
By (12) , the two outside distances satisfy
∑ p , q μ p q ( t ) u p 2 = ∑ p , q μ p q ( t ) u q 2 = ∑ p c p u p 2 = J . \sum_{p,q}\mu_{pq}(t)u_p^2=\sum_{p,q}\mu_{pq}(t)u_q^2=\sum_p c_pu_p^2=J. p , q ∑ μ pq ( t ) u p 2 = p , q ∑ μ pq ( t ) u q 2 = p ∑ c p u p 2 = J .
For the middle distance, group the copies by their neighbors in B B B and apply (8) :
∑ p , q μ p q ( t ) ∥ z b p − z b q ∥ 2 = ∑ b , b ′ s b P b b ′ ( t ) ∥ z b − z b ′ ∥ 2 ≤ 2 t I . \sum_{p,q}\mu_{pq}(t)\lVert z_{b_p}-z_{b_q}\rVert^2=\sum_{b,b'}s_bP_{bb'}(t)\lVert z_b-z_{b'}\rVert^2\leq2tI. p , q ∑ μ pq ( t ) ∥ z b p − z b q ∥ 2 = b , b ′ ∑ s b P b b ′ ( t ) ∥ z b − z b ′ ∥ 2 ≤ 2 t I .
The triangle inequality for the Euclidean norm with weights μ p q ( t ) \mu_{pq}(t) μ pq ( t ) therefore gives
S t ≤ J + 2 t I + J . \sqrt{S_t}\leq\sqrt{J}+\sqrt{2tI}+\sqrt{J}. S t ≤ J + 2 t I + J .
Substituting into (13), we obtain
E G A ( d 2 ) − E G [ A ] ( d 2 ) ≤ 1 2 ∫ 0 ∞ e − t ( 4 J + 4 2 t I J + 2 t I ) d t = 2 J + 2 π I J + I , \begin{aligned}
\mathcal{E}_{G^A}(d^2)-\mathcal{E}_{G[A]}(d^2)\leq\frac{1}{2}\int_0^\infty e^{-t}(4J+4\sqrt{2tIJ}+2tI)\,\mathrm{d}t \\
&=2J+\sqrt{2\pi IJ}+I,
\end{aligned} E G A ( d 2 ) − E G [ A ] ( d 2 ) ≤ 2 1 ∫ 0 ∞ e − t ( 4 J + 4 2 t I J + 2 t I ) d t = 2 J + 2 π I J + I ,
where we used ∫ 0 ∞ e − t t d t = π / 2 \int_0^\infty e^{-t}\sqrt{t}\,\mathrm{d}t=\sqrt{\pi}/2 ∫ 0 ∞ e − t t d t = π /2 . Since E G ( d 2 ) = E G [ A ] ( d 2 ) + I + J \mathcal{E}_G(d^2)=\mathcal{E}_{G[A]}(d^2)+I+J E G ( d 2 ) = E G [ A ] ( d 2 ) + I + J , this proves (9) under our assumption that all s b s_b s b are positive.
To remove that assumption, let B 0 = { b ∈ B : s b > 0 } B_0=\{b\in B:s_b>0\} B 0 = { b ∈ B : s b > 0 } . If B 0 = ∅ B_0=\varnothing B 0 = ∅ , no edges join A A A to B B B , so G A = G [ A ] G^A=G[A] G A = G [ A ] and the bound is immediate. Otherwise, first remove the vertices of B ∖ B 0 B\setminus B_0 B ∖ B 0 . Throughout this process their neighbors lie in B B B , so the edges within A A A and the edges from A A A to B 0 B_0 B 0 keep their weights. Let G ′ G' G ′ be the resulting graph on A ⊔ B 0 A\sqcup B_0 A ⊔ B 0 . By (7) , ( G ′ ) A = G A (G')^A=G^A ( G ′ ) A = G A , and the graph on B 0 B_0 B 0 is ( G [ B ] ) B 0 (G[B])^{B_0} ( G [ B ] ) B 0 . Applying (6) to vectors representing the distances on B B B , we find that the internal and crossing energies of G ′ G' G ′ satisfy I ′ = E G ′ [ B 0 ] ( d 2 ) ≤ I I'=\mathcal{E}_{G'[B_0]}(d^2)\leq I I ′ = E G ′ [ B 0 ] ( d 2 ) ≤ I and J ′ = J J'=J J ′ = J . Thus, E G ′ ( d 2 ) ≤ E G ( d 2 ) \mathcal{E}_{G'}(d^2)\leq\mathcal{E}_G(d^2) E G ′ ( d 2 ) ≤ E G ( d 2 ) , and every vertex of B 0 B_0 B 0 still has an edge to A A A . The bound already proved for G ′ G' G ′ therefore gives (9) for G G G as well.
From energy bounds to Lipschitz extensions
We now turn the one-set estimate into a bound on the Lipschitz constant of an extension. Two factors matter separately: the growth on pairs of points inside the eliminated set, and the growth on pairs crossing it. We first show that an energy bound with separate coefficients for the internal and the crossing energy yields an extension with separate Lipschitz factors on the two kinds of pairs ((Lemma 4.7)). Averaging these extensions over the choice of the set added last then reduces the growth at the step with r r r sets to a factor of 1 + η / r 1+\eta/r 1 + η / r ((Lemma 4.8)). Finally, we choose the two factors to minimize the exponent η \eta η .
We write H 0 H_0 H 0 for the assumption that every U i U_i U i is Euclidean in X X X , and H 1 H_1 H 1 for the assumption that every U i U_i U i is 1-strongly Euclidean in X X X , that is, U i ∪ { x } U_i\cup\{x\} U i ∪ { x } is Euclidean for every x ∈ X x\in X x ∈ X . Both conditions are preserved when a set or the ambient space is restricted. Under H 1 H_1 H 1 , we also sharpen the one-set estimate (9) .
Lemma 4.7. Let Y = A ⊔ B Y=A\sqcup B Y = A ⊔ B be a finite pseudometric space with A ≠ ∅ A\ne\varnothing A = ∅ , and let t , s ≥ 1 t,s\ge1 t , s ≥ 1 . Suppose that every nonnegatively weighted graph G G G on Y Y Y with no edges inside A A A satisfies
E G A ( d 2 ) ≤ t 2 I + s 2 J , \mathcal{E}_{G^A}(d^2)\le t^2I+s^2J, E G A ( d 2 ) ≤ t 2 I + s 2 J , where I I I and J J J are as in Lemma 4.6 . Then every 1 1 1 -Lipschitz map f : A → ℓ 2 f:A\to\ell_2 f : A → ℓ 2 has an extension F : Y → ℓ 2 F:Y\to\ell_2 F : Y → ℓ 2 such that
∥ F ( b ) − F ( b ′ ) ∥ ≤ t d ( b , b ′ ) , ∥ F ( a ) − F ( b ) ∥ ≤ s d ( a , b ) , \lVert F(b)-F(b')\rVert\le td(b,b'),\qquad\lVert F(a)-F(b)\rVert\le sd(a,b), ∥ F ( b ) − F ( b ′ )∥ ≤ t d ( b , b ′ ) , ∥ F ( a ) − F ( b )∥ ≤ s d ( a , b ) , where a ∈ A a\in A a ∈ A , b , b ′ ∈ B b,b'\in B b , b ′ ∈ B . The extension can be chosen with values in conv ( f ( A ) ) \operatorname{conv}(f(A)) conv ( f ( A )) .
Proof. The case B = ∅ B=\varnothing B = ∅ is immediate. Otherwise, by Lemma 4.4 and the nonnegative weights of G A G^A G A ,
min F ∣ A = f E G ( F ) = E G A ( f ) ≤ E G A ( d 2 ) ≤ t 2 I + s 2 J . \min_{F\mid_A=f}\mathcal{E}_G(F)=\mathcal{E}_{G^A}(f)\le\mathcal{E}_{G^A}(d^2)\le t^2I+s^2J. F ∣ A = f min E G ( F ) = E G A ( f ) ≤ E G A ( d 2 ) ≤ t 2 I + s 2 J . Projection onto the compact convex set conv ( f ( A ) ) \operatorname{conv}(f(A)) conv ( f ( A )) fixes f ( A ) f(A) f ( A ) and does not increase distances, so this inequality remains valid when all extension values are restricted to that set. Let F \mathcal{F} F be the resulting compact convex set of extensions, and let W \mathcal{W} W be the probability simplex on unordered pairs with at least one endpoint in B B B . For such a pair, let b x y = t 2 b_{xy}=t^2 b x y = t 2 if x , y ∈ B x,y\in B x , y ∈ B and b x y = s 2 b_{xy}=s^2 b x y = s 2 otherwise, and define
Φ ( F , w ) = ∑ { x , y } w x y ( ∥ F ( x ) − F ( y ) ∥ 2 − b x y d ( x , y ) 2 ) . \Phi(F,w)=\sum_{\{x,y\}}w_{xy}\left(\lVert F(x)-F(y)\rVert^2-b_{xy}d(x,y)^2\right). Φ ( F , w ) = { x , y } ∑ w x y ( ∥ F ( x ) − F ( y ) ∥ 2 − b x y d ( x , y ) 2 ) . The energy inequality gives min F ∈ F Φ ( F , w ) ≤ 0 \min_{F\in\mathcal{F}}\Phi(F,w)\le0 min F ∈ F Φ ( F , w ) ≤ 0 for every w ∈ W w\in\mathcal{W} w ∈ W . As in Section A, convexity in F F F and affinity in w w w give
min F ∈ F max w ∈ W Φ ( F , w ) = max w ∈ W min F ∈ F Φ ( F , w ) ≤ 0. \min_{F\in\mathcal{F}}\max_{w\in\mathcal{W}}\Phi(F,w)=\max_{w\in\mathcal{W}}\min_{F\in\mathcal{F}}\Phi(F,w)\le0. F ∈ F min w ∈ W max Φ ( F , w ) = w ∈ W max F ∈ F min Φ ( F , w ) ≤ 0. Taking w w w supported on each individual pair proves the two bounds, including pairs at pseudodistance zero.
We now average such extensions over the choice of the set added last.
Lemma 4.8. Suppose that an added set satisfying a given geometric condition always admits the extension in Lemma 4.7 with fixed factors t , s ≥ 1 t,s\ge1 t , s ≥ 1 , and that the condition is preserved by restriction. Let η = max { t − 1 , 2 ( s − 1 ) } \eta=\max\{t-1,2(s-1)\} η = max { t − 1 , 2 ( s − 1 )} . If A ≠ ∅ A\ne\varnothing A = ∅ and X ∖ A X\setminus A X ∖ A is partitioned into m m m such sets, every 1 1 1 -Lipschitz map f : A → ℓ 2 f:A\to\ell_2 f : A → ℓ 2 extends to X X X with Lipschitz constant at most
P m = ∏ r = 1 m ( 1 + η r ) . (14) P_m=\prod_{r=1}^{m}\left(1+\frac{\eta}{r}\right). \tag*{(14)} P m = r = 1 ∏ m ( 1 + r η ) . (14) For fixed η \eta η , this is O ( ( m + 1 ) η ) O((m+1)^\eta) O (( m + 1 ) η ) . If η ≥ 1 \eta\ge1 η ≥ 1 , it is at most ( m + 1 ) η (m+1)^\eta ( m + 1 ) η .
Proof. We induct on m m m , starting with P 0 = 1 P_0=1 P 0 = 1 . Let X ∖ A = B 1 ⊔ ⋯ ⊔ B m X\setminus A=B_1\sqcup\cdots\sqcup B_m X ∖ A = B 1 ⊔ ⋯ ⊔ B m . For each i ∈ [ m ] i\in[m] i ∈ [ m ] , first extend the prescribed map to X ∖ B i X\setminus B_i X ∖ B i with Lipschitz constant at most P m − 1 P_{m-1} P m − 1 , and then add B i B_i B i using Lemma 4.7 , after scaling by P m − 1 P_{m-1} P m − 1 . Call the resulting map F i F_i F i . All these maps take values in the same Hilbert space and agree with the prescribed map on A A A , so their arithmetic average
F ( x ) = 1 m ∑ i = 1 m F i ( x ) F(x)=\frac{1}{m}\sum_{i=1}^{m}F_i(x) F ( x ) = m 1 i = 1 ∑ m F i ( x )
also extends that map. For a pair inside one B j B_j B j , only F j F_j F j uses the factor t t t . For a pair in two distinct complementary sets, exactly two summands use the factor s s s . A pair joining A A A to B j B_j B j has only one summand using s s s , and pairs inside A A A are unchanged. The triangle inequality therefore gives
Lip ( F ) ≤ P m − 1 max { 1 + t − 1 m , 1 + 2 ( s − 1 ) m } = P m . \operatorname{Lip}(F) \le P_{m-1}\max\left\{1+\frac{t-1}{m},1+\frac{2(s-1)}{m}\right\}=P_m. Lip ( F ) ≤ P m − 1 max { 1 + m t − 1 , 1 + m 2 ( s − 1 ) } = P m .
Finally, for m ≥ 1 m\ge1 m ≥ 1 ,
P m ≤ exp ( η ∑ r = 1 m 1 r ) ≤ e η ( m + 1 ) η . P_m\le\exp\left(\eta\sum_{r=1}^{m}\frac{1}{r}\right)\le e^\eta(m+1)^\eta. P m ≤ exp ( η r = 1 ∑ m r 1 ) ≤ e η ( m + 1 ) η .
When η ≥ 1 \eta\ge1 η ≥ 1 , the sharper estimate follows by multiplying 1 + η / r ≤ ( 1 + 1 / r ) η 1+\eta/r\le(1+1/r)^\eta 1 + η / r ≤ ( 1 + 1/ r ) η . For 0 < η < 1 0<\eta<1 0 < η < 1 , we use the preceding exponential estimate instead.
The exponent for Euclidean sets
For a graph with no edges inside A A A , the estimate (9) is
E G A ( d 2 ) ≤ I + 2 J + 2 π I J . \mathcal{E}_{GA}(d^2)\le I+2J+\sqrt{2\pi IJ}. E G A ( d 2 ) ≤ I + 2 J + 2 π I J .
Consequently, Lemma 4.7 applies whenever
t ≥ 1 , s ≥ 2 , ( t 2 − 1 ) ( s 2 − 2 ) ≥ π 2 . (15) t\ge1,\qquad s\ge\sqrt{2},\qquad(t^2-1)(s^2-2)\ge\frac{\pi}{2}. \tag*{(15)} t ≥ 1 , s ≥ 2 , ( t 2 − 1 ) ( s 2 − 2 ) ≥ 2 π . (15)
since 2 π I J ≤ ( t 2 − 1 ) I + ( s 2 − 2 ) J \sqrt{2\pi IJ}\le(t^2-1)I+(s^2-2)J 2 π I J ≤ ( t 2 − 1 ) I + ( s 2 − 2 ) J under this condition. To minimize the exponent supplied by Lemma 4.8 , take t = 1 + η t=1+\eta t = 1 + η and s = 1 + η / 2 s=1+\eta/2 s = 1 + η /2 : for a fixed permitted exponent η \eta η , these are the largest available factors.
Let η 0 > 1 \eta_0>1 η 0 > 1 be the unique root of
η 0 ( η 0 + 2 ) ( η 0 2 + 4 η 0 − 4 ) = 2 π . (16) \eta_0(\eta_0+2)(\eta_0^2+4\eta_0-4)=2\pi. \tag*{(16)} η 0 ( η 0 + 2 ) ( η 0 2 + 4 η 0 − 4 ) = 2 π . (16)
The left side is strictly increasing for η 0 ≥ 1 \eta_0\ge1 η 0 ≥ 1 , equals 3 at 1, and tends to infinity. Numerically, η 0 = 1.1275681325 … \eta_0=1.1275681325\ldots η 0 = 1.1275681325 … . The choices t = 1 + η 0 t=1+\eta_0 t = 1 + η 0 and s = 1 + η 0 / 2 s=1+\eta_0/2 s = 1 + η 0 /2 give equality in (15) . Thus, a complement consisting of m m m Euclidean sets admits extension factor at most ( m + 1 ) η 0 (m+1)^{\eta_0} ( m + 1 ) η 0 . The simpler exponent 8 / 7 8/7 8/7 is also valid, since t = 15 / 7 t=15/7 t = 15/7 and s = 11 / 7 s=11/7 s = 11/7 give ( t 2 − 1 ) ( s 2 − 2 ) = 4048 / 2401 > π / 2 (t^2-1)(s^2-2)=4048/2401>\pi/2 ( t 2 − 1 ) ( s 2 − 2 ) = 4048/2401 > π /2 .
Eliminating one 1-strongly Euclidean set
Lemma 4.9. In the setting of Lemma 4.6 , suppose that B ∪ { a } B\cup\{a\} B ∪ { a } is Euclidean for every a ∈ A a\in A a ∈ A . Then
E G A ( d 2 ) ≤ E G [ A ] ( d 2 ) + 1 2 ( J + J + I + I ( I + J ) ) 2 . (17) \mathcal{E}_{GA}(d^2)\le\mathcal{E}_{G[A]}(d^2)+\frac{1}{2}\left(\sqrt{J}+\sqrt{J+I+\sqrt{I(I+J)}}\right)^2. \tag*{(17)} E G A ( d 2 ) ≤ E G [ A ] ( d 2 ) + 2 1 ( J + J + I + I ( I + J ) ) 2 . (17)
Proof. Use the notation s b s_b s b , c p c_p c p , a p a_p a p , b p b_p b p , C C C , K K K , P ( t ) P(t) P ( t ) , μ p q ( t ) \mu_{pq}(t) μ pq ( t ) from the proof of Lemma 4.6 , and first assume that s b > 0 s_b>0 s b > 0 for every b ∈ B b\in B b ∈ B . Let
R = ( Id + K ) − 1 = ∫ 0 ∞ e − t P ( t ) d t , R=(\operatorname{Id}+K)^{-1}=\int_0^\infty e^{-t}P(t)\,dt, R = ( Id + K ) − 1 = ∫ 0 ∞ e − t P ( t ) d t ,
and
ν p q = ∫ 0 ∞ e − t μ p q ( t ) d t = c p R b p b q c q s b q . \nu_{pq}=\int_0^\infty e^{-t}\mu_{pq}(t)\,dt=c_pR b_pb_q\frac{c_q}{s_{b_q}}. ν pq = ∫ 0 ∞ e − t μ pq ( t ) d t = c p R b p b q s b q c q .
These weights are symmetric and have marginals ∑ q ν p q = c p \sum_q \nu_{pq}=c_p ∑ q ν pq = c p and ∑ p ν p q = c q \sum_p \nu_{pq}=c_q ∑ p ν pq = c q . Equation (13) gives
2 ( E G A ( d 2 ) − E G [ A ] ( d 2 ) ) = ∑ p , q ν p q d ( a p , a q ) 2 . (18) 2\left(\mathcal{E}_{G_A}(d^2)-\mathcal{E}_{G[A]}(d^2)\right)=\sum_{p,q}\nu_{pq}d(a_p,a_q)^2. \tag*{(18)} 2 ( E G A ( d 2 ) − E G [ A ] ( d 2 ) ) = p , q ∑ ν pq d ( a p , a q ) 2 . (18)
Choose a tight Euclidean realization z b z_b z b of B B B . By the one-point condition and Fact 2.4, for each a ∈ A a\in A a ∈ A there are a vector u a u_a u a and a height h a ≥ 0 h_a\ge0 h a ≥ 0 such that for every b ∈ B b\in B b ∈ B we have
d ( a , b ) 2 = ∥ u a − z b ∥ 2 + h a 2 . d(a,b)^2=\lVert u_a-z_b\rVert^2+h_a^2. d ( a , b ) 2 = ∥ u a − z b ∥ 2 + h a 2 .
Set m = ⟨ z , ( Id − R ) z ⟩ C m=\langle z,(\operatorname{Id}-R)z\rangle_C m = ⟨ z , ( Id − R ) z ⟩ C , v = ⟨ ( Id − R ) z , ( Id − R ) z ⟩ C v=\langle(\operatorname{Id}-R)z,(\operatorname{Id}-R)z\rangle_C v = ⟨( Id − R ) z , ( Id − R ) z ⟩ C , and T = ∑ p , q ν p q d ( a p , b q ) 2 T=\sum_{p,q}\nu_{pq}d(a_p,b_q)^2 T = ∑ p , q ν pq d ( a p , b q ) 2 . Expanding the squared distances and using the marginals gives
T − J = 2 ∑ p c p ⟨ u a p , ( ( Id − R ) z ) b p ⟩ = 2 m + 2 ∑ p c p ⟨ u a p − z b p , ( ( Id − R ) z ) b p ⟩ ≤ 2 m + 2 J v . \begin{aligned}
T-J&=2\sum_p c_p\langle u_{a_p},((\operatorname{Id}-R)z)_{b_p}\rangle\\
&=2m+2\sum_p c_p\langle u_{a_p}-z_{b_p},((\operatorname{Id}-R)z)_{b_p}\rangle\\
&\le2m+2\sqrt{Jv}.
\end{aligned} T − J = 2 p ∑ c p ⟨ u a p , (( Id − R ) z ) b p ⟩ = 2 m + 2 p ∑ c p ⟨ u a p − z b p , (( Id − R ) z ) b p ⟩ ≤ 2 m + 2 J v .
where the averaged squared norms of z b z_b z b and the height terms cancel, and ∑ p c p ∥ u a p − z b p ∥ 2 ≤ J \sum_p c_p\lVert u_{a_p}-z_{b_p}\rVert^2\le J ∑ p c p ∥ u a p − z b p ∥ 2 ≤ J .
Let w λ w_\lambda w λ be the squared norms of the spectral components of z z z for K K K in the C C C inner product. Then
I = ∑ λ λ w λ , m = ∑ λ λ 1 + λ w λ , v = ∑ λ λ 2 ( 1 + λ ) 2 w λ . I=\sum_\lambda\lambda w_\lambda,\qquad m=\sum_\lambda\frac{\lambda}{1+\lambda}w_\lambda,\qquad v=\sum_\lambda\frac{\lambda^2}{(1+\lambda)^2}w_\lambda. I = λ ∑ λ w λ , m = λ ∑ 1 + λ λ w λ , v = λ ∑ ( 1 + λ ) 2 λ 2 w λ .
Cauchy–Schwarz gives m 2 ≤ I ( m − v ) m^2\le I(m-v) m 2 ≤ I ( m − v ) , and 0 ≤ m ≤ I 0\le m\le I 0 ≤ m ≤ I . If I > 0 I>0 I > 0 , then m = I ( 1 + q ) / 2 m=I(1+q)/2 m = I ( 1 + q ) /2 with − 1 ≤ q ≤ 1 -1\le q\le1 − 1 ≤ q ≤ 1 . Then
2 m + 2 J v ≤ I + I q + I J 1 − q 2 ≤ I + I ( I + J ) . 2m+2\sqrt{Jv}\le I+Iq+\sqrt{IJ}\sqrt{1-q^2}\le I+\sqrt{I(I+J)}. 2 m + 2 J v ≤ I + I q + I J 1 − q 2 ≤ I + I ( I + J ) .
If I = 0 I=0 I = 0 , the spectral expressions give m = v = 0 m=v=0 m = v = 0 , so in either case T ≤ J + I + I ( I + J ) T\le J+I+\sqrt{I(I+J)} T ≤ J + I + I ( I + J ) . Finally, the triangle inequality along a p a_p a p , b q b_q b q , a q a_q a q and the Euclidean norm with weights ν p q \nu_{pq} ν pq give
( ∑ p , q ν p q d ( a p , a q ) 2 ) 1 / 2 ≤ T + ( ∑ p , q ν p q d ( a q , b q ) 2 ) 1 / 2 = T + J . \left(\sum_{p,q}\nu_{pq}d(a_p,a_q)^2\right)^{1/2}\le\sqrt{T}+\left(\sum_{p,q}\nu_{pq}d(a_q,b_q)^2\right)^{1/2}=\sqrt{T}+\sqrt{J}. ( p , q ∑ ν pq d ( a p , a q ) 2 ) 1/2 ≤ T + ( p , q ∑ ν pq d ( a q , b q ) 2 ) 1/2 = T + J .
This proves (17) . If some s b s_b s b vanish, use the reduction at the end of the proof of Lemma 4.6. It preserves E G [ A ] ( d 2 ) \mathcal{E}_{G[A]}(d^2) E G [ A ] ( d 2 ) and J J J , decreases I I I , and preserves the one-point condition, and the right side of (17) is increasing in I I I . If all s b s_b s b vanish, G A = G [ A ] G^A=G[A] G A = G [ A ] and the conclusion is immediate. □
The exponent for 1-strongly Euclidean sets
To use Lemma 4.8 with exponent η \eta η , we again take t = 1 + η t=1+\eta t = 1 + η and s = 1 + η / 2 s=1+\eta/2 s = 1 + η /2 . For I + J > 0 I+J>0 I + J > 0 , set z = 1 − I / ( I + J ) z=\sqrt{1-\sqrt{I/(I+J)}} z = 1 − I / ( I + J ) . Then
I I + J = ( 1 − z 2 ) 2 , J I + J = z 2 ( 2 − z 2 ) , 0 ≤ z ≤ 1 , \frac{I}{I+J}=(1-z^2)^2,\qquad\frac{J}{I+J}=z^2(2-z^2),\qquad0\le z\le1, I + J I = ( 1 − z 2 ) 2 , I + J J = z 2 ( 2 − z 2 ) , 0 ≤ z ≤ 1 ,
and the added-energy bound in (17) , divided by I + J I+J I + J , equals 1 2 ( z + 1 ) 2 ( 2 − z 2 ) \frac{1}{2}(z+1)^2(2-z^2) 2 1 ( z + 1 ) 2 ( 2 − z 2 ) . Thus, that bound is at most t 2 I + s 2 J t^2I+s^2J t 2 I + s 2 J for all $I,J\ge 0 precisely when
D η ( z ) = ( 1 + η ) 2 ( 1 − z 2 ) 2 + ( 1 + η / 2 ) 2 z 2 ( 2 − z 2 ) − 1 2 ( z + 1 ) 2 ( 2 − z 2 ) ≥ 0 ( 0 ≤ z ≤ 1 ) . (19) D_\eta(z)=(1+\eta)^2(1-z^2)^2+(1+\eta/2)^2z^2(2-z^2)-\frac{1}{2}(z+1)^2(2-z^2)\ge0\qquad(0\le z\le1). \tag*{(19)} D η ( z ) = ( 1 + η ) 2 ( 1 − z 2 ) 2 + ( 1 + η /2 ) 2 z 2 ( 2 − z 2 ) − 2 1 ( z + 1 ) 2 ( 2 − z 2 ) ≥ 0 ( 0 ≤ z ≤ 1 ) . (19)
Define the exponent exactly by
η 1 = min { η ≥ 0 : D η ( z ) ≥ 0 for every z ∈ [ 0 , 1 ] } . (20) \eta_{1}=\min\{\eta\ge0:D_{\eta}(z)\ge0\text{ for every }z\in[0,1]\}. \tag*{(20)} η 1 = min { η ≥ 0 : D η ( z ) ≥ 0 for every z ∈ [ 0 , 1 ]} . (20)
The minimum exists: D 2 ( z ) ≥ 0 D_{2}(z)\ge0 D 2 ( z ) ≥ 0 on [ 0 , 1 ] [0,1] [ 0 , 1 ] , the condition is closed, and D η ( z ) D_{\eta}(z) D η ( z ) is strictly increasing in η ≥ 0 \eta\ge0 η ≥ 0 for every z z z . This is the same optimization of the internal and crossing factors as under H 0 H_{0} H 0 .
For completeness, we certify the useful numerical interval
0.8941 < η 1 < 0.89417. 0.8941<\eta_{1}<0.89417. 0.8941 < η 1 < 0.89417.
Take a = 1 2 + η + 3 4 η 2 a=\frac{1}{2}+\eta+\frac{3}{4}\eta^{2} a = 2 1 + η + 4 3 η 2 . Expanding the gap gives
D η ( z ) = η 2 + 2 η − 2 z + z 2 2 + z 3 + a ( z 4 − 2 z 2 ) D_{\eta}(z)=\eta^{2}+2\eta-2z+\frac{z^{2}}{2}+z^{3}+a(z^{4}-2z^{2}) D η ( z ) = η 2 + 2 η − 2 z + 2 z 2 + z 3 + a ( z 4 − 2 z 2 )
D η ′ ( z ) = ( z + 1 ) ( 4 a z 2 + ( 3 − 4 a ) z − 2 ) . D_{\eta}'(z)=(z+1)(4az^{2}+(3-4a)z-2). D η ′ ( z ) = ( z + 1 ) ( 4 a z 2 + ( 3 − 4 a ) z − 2 ) .
Take η = 89417 / 100000 \eta=89417/100000 η = 89417/100000 and z 0 = 90189211 / 10 8 z_{0}=90189211/10^{8} z 0 = 90189211/1 0 8 . On [ 0 , 3 / 4 ] [0,3/4] [ 0 , 3/4 ] , the bounds z 3 ≤ 9 z / 16 z^{3}\le9z/16 z 3 ≤ 9 z /16 and z 2 ≤ 3 z / 4 z^{2}\le3z/4 z 2 ≤ 3 z /4 give D η ′ ( z ) ≤ − 2 + ( 13 / 4 − 7 a / 4 ) z < 0 D_{\eta}'(z)\le-2+(13/4-7a/4)z<0 D η ′ ( z ) ≤ − 2 + ( 13/4 − 7 a /4 ) z < 0 . On [ 3 / 4 , 1 ] [3/4,1] [ 3/4 , 1 ] ,
D η ′ ′ ( z ) = 1 + 6 z + 12 a z 2 − 4 a ≥ 11 2 + 11 a 4 > 10. D_{\eta}''(z)=1+6z+12az^{2}-4a\ge\frac{11}{2}+\frac{11a}{4}>10. D η ′′ ( z ) = 1 + 6 z + 12 a z 2 − 4 a ≥ 2 11 + 4 11 a > 10.
Exact rational substitution gives D η ( z 0 ) > 9 ⋅ 10 − 8 D_{\eta}(z_{0})>9\cdot10^{-8} D η ( z 0 ) > 9 ⋅ 1 0 − 8 and ∣ D η ′ ( z 0 ) ∣ < 10 − 7 \lvert D_{\eta}'(z_{0})\rvert<10^{-7} ∣ D η ′ ( z 0 )∣ < 1 0 − 7 . Consequently, throughout this latter interval,
D η ( z ) ≥ D η ( z 0 ) + D η ′ ( z 0 ) ( z − z 0 ) + 5 ( z − z 0 ) 2 ≥ D η ( z 0 ) − D η ′ ( z 0 ) 2 20 > 0. D_{\eta}(z)\ge D_{\eta}(z_{0})+D_{\eta}'(z_{0})(z-z_{0})+5(z-z_{0})^{2}\ge D_{\eta}(z_{0})-\frac{D_{\eta}'(z_{0})^{2}}{20}>0. D η ( z ) ≥ D η ( z 0 ) + D η ′ ( z 0 ) ( z − z 0 ) + 5 ( z − z 0 ) 2 ≥ D η ( z 0 ) − 20 D η ′ ( z 0 ) 2 > 0.
Monotonicity handles [ 0 , 3 / 4 ] [0,3/4] [ 0 , 3/4 ] , and strict positivity on the compact interval gives η 1 < 0.89417 \eta_{1}<0.89417 η 1 < 0.89417 . For η = 8941 / 10000 \eta=8941/10000 η = 8941/10000 , direct substitution yields
2 D η ( 9 / 10 ) = − 300526077 2000000000000 < 0 , 2D_{\eta}(9/10)=-\frac{300526077}{2000000000000}<0, 2 D η ( 9/10 ) = − 2000000000000 300526077 < 0 ,
so η 1 > 0.8941 \eta_{1}>0.8941 η 1 > 0.8941 .
Proof of Lemma 1.5
Proof of Lemma 1.5. Fix i ∈ [ k ] i\in[k] i ∈ [ k ] and a Lipschitz map f : U i → ℓ 2 f:U_{i}\to\ell_{2} f : U i → ℓ 2 . A constant map extends constantly, and scaling handles other Lipschitz constants, so we may assume that Lip ( f ) = 1 \operatorname{Lip}(f)=1 Lip ( f ) = 1 . Keep U i U_{i} U i whole and partition its complement into the sets B j = U j ∖ ( U i ∪ U 1 ∪ ⋯ ∪ U j − 1 ) B_{j}=U_{j}\setminus(U_{i}\cup U_{1}\cup\cdots\cup U_{j-1}) B j = U j ∖ ( U i ∪ U 1 ∪ ⋯ ∪ U j − 1 ) for j ≠ i j\ne i j = i , discarding empty ones. This gives m ≤ k − 1 m\le k-1 m ≤ k − 1 sets, and each B j B_{j} B j is Euclidean, respectively 1-strongly Euclidean, in X X X , because both conditions are preserved by restriction. For Euclidean sets, (9) and (15) show that Lemma 4.7 applies with t = 1 + η 0 t=1+\eta_{0} t = 1 + η 0 and s = 1 + η 0 / 2 s=1+\eta_{0}/2 s = 1 + η 0 /2 . For 1-strongly Euclidean sets, Lemma 4.9 and (20) show that it applies with t = 1 + η 1 t=1+\eta_{1} t = 1 + η 1 and s = 1 + η 1 / 2 s=1+\eta_{1}/2 s = 1 + η 1 /2 . In either case max { t − 1 , 2 ( s − 1 ) } = η j \max\{t-1,2(s-1)\}=\eta_{j} max { t − 1 , 2 ( s − 1 )} = η j , so Lemma 4.8 extends f f f to X X X with Lipschitz constant at most P m ≤ P k − 1 P_{m}\le P_{k-1} P m ≤ P k − 1 , with P m P_{m} P m as in (14) . Hence,
e ( U i ) ≤ ∏ r = 1 k − 1 ( 1 + η j r ) . e(U_{i})\le\prod_{r=1}^{k-1}\left(1+\frac{\eta_{j}}{r}\right). e ( U i ) ≤ r = 1 ∏ k − 1 ( 1 + r η j ) .
Since η 0 > 1 \eta_{0}>1 η 0 > 1 , this is at most k η 0 k^{\eta_{0}} k η 0 under H 0 H_{0} H 0 , and since η 1 < 1 \eta_{1}<1 η 1 < 1 , it is at most e η 1 k η 1 e^{\eta_{1}}k^{\eta_{1}} e η 1 k η 1 under H 1 H_{1} H 1 .
The proof applies to every d ∈ D d \in\mathcal{D} d ∈ D and every partition into at most k k k sets, with the same product, so it also bounds the uniform constants of § 4.1 .
Corollary 4.10. For every k ≥ 1 k \ge1 k ≥ 1 , we have e k ≤ ∏ r = 0 k − 1 ( 1 + η 0 r ) ≤ k η 0 e_k \le\prod_{r=0}^{k-1}\left(1+\frac{\eta_0}{r}\right) \le k^{\eta_0} e k ≤ ∏ r = 0 k − 1 ( 1 + r η 0 ) ≤ k η 0 , and Q k = e k 2 ≤ k 2 η 0 Q_k=e_k^2 \le k^{2\eta_0} Q k = e k 2 ≤ k 2 η 0 . In particular, for every partition of X X X into at most k k k Euclidean sets and every weighted graph G G G on X X X , the canonical graph G U i G^{U_i} G U i is a Euclidean metric sparsifier for D \mathcal{D} D of quality at most k 2 η 0 k^{2\eta_0} k 2 η 0 .
Proof. The first bound is the proof of Lemma 1.5, and Theorem 4.2 gives Q k = e k 2 Q_k=e_k^2 Q k = e k 2 . The proof of Theorem 4.2 in Section A shows that G U i G^{U_i} G U i attains quality e ( D , U i ) 2 ≤ e k 2 e(\mathcal{D},U_i)^2 \le e_k^2 e ( D , U i ) 2 ≤ e k 2 .
**Proof of Theorems 1.2 and 1.3 .** Discard empty sets among the $U_i$. Lemma 1.5 gives $e_{\max}=O(k^{\eta_0})$, respectively $O(k^{\eta_1})$, and {xref}Theorem 1.4` gives the distortion bound.
Lower Bound for Euclidean Sets
Let H ≠ { 0 } H \ne\{0\} H = { 0 } be a linear subspace of a finite-dimensional normed space W W W that inherits the norm of W W W . The relative projection constant λ ( H , W ) \lambda(H,W) λ ( H , W ) denotes the least norm of a linear projection P : W → H P:W \to H P : W → H that fixes H H H . Formally, this is
λ ( H , W ) = inf { ∥ P ∥ W → H ∣ P : W → H linear , P ∣ H = I H } . \lambda(H,W)=\inf\{\lVert P\rVert_{W\to H}\mid P:W\to H\text{ linear},\ P|_H=I_H\}. λ ( H , W ) = inf {∥ P ∥ W → H ∣ P : W → H linear , P ∣ H = I H } .
where ∥ P ∥ W → H = sup w ∈ W , w ≠ 0 ∥ P w ∥ / ∥ w ∥ \lVert P\rVert_{W\to H}=\sup_{w\in W,w\ne0}\lVert Pw\rVert/\lVert w\rVert ∥ P ∥ W → H = sup w ∈ W , w = 0 ∥ P w ∥ / ∥ w ∥ . Theorem 5.1 gives a classical connection between linear extensions and Lipschitz extensions.
Theorem 5.1 ([17 ] ). For finite-dimensional H ⊂ W H \subset W H ⊂ W , we have
inf F : W → H F ∣ H = I H Lip ( F ) = λ ( H , W ) . \inf_{\substack{F:W\to H\\F|_H=I_H}}\operatorname{Lip}(F)=\lambda(H,W). F : W → H F ∣ H = I H inf Lip ( F ) = λ ( H , W ) .
The cycle space of a graph
Let G = ( V , E ) G=(V,E) G = ( V , E ) be a connected graph with n n n vertices, m m m edges and girth g ≥ 3 g \ge3 g ≥ 3 . Fix an arbitrary orientation of the edges, so that each edge e = ( u , v ) e=(u,v) e = ( u , v ) goes from u u u to v v v . We view x ∈ R E x\in\mathbb{R}^{E} x ∈ R E as a flow that sends x e x_e x e units along edge e e e . When x e < 0 x_e<0 x e < 0 , the flow goes backwards along the edge. We identify each edge with its indicator vector in R E \mathbb{R}^{E} R E . The net inflow at a vertex w w w is ( ∂ x ) w = ∑ e ∈ δ − ( w ) x e − ∑ e ∈ δ + ( w ) x e (\partial x)_w=\sum_{e\in\delta^-(w)}x_e-\sum_{e\in\delta^+(w)}x_e ( ∂ x ) w = ∑ e ∈ δ − ( w ) x e − ∑ e ∈ δ + ( w ) x e . Moreover, the net inflow is a linear map ∂ : R E → R V \partial:\mathbb{R}^{E}\to\mathbb{R}^{V} ∂ : R E → R V . The cycle space Z = ker ∂ ⊂ R E Z=\ker\partial\subset\mathbb{R}^{E} Z = ker ∂ ⊂ R E consists of the flows that are conserved at every vertex. For instance, it contains every cycle, which is + 1 +1 + 1 on the edges the cycle traverses along their orientation, − 1 -1 − 1 on those it traverses against it, and 0 0 0 elsewhere. The dimension of the cycle space is r = m − n + 1 r=m-n+1 r = m − n + 1 , and we assume r > 0 r>0 r > 0 .
Next, choose a root o o o and a spanning tree of G G G of depth D D D . Let p v ∈ R E p_v\in\mathbb{R}^{E} p v ∈ R E be the signed indicator of the tree path from o o o to v v v , that is, the unit flow from o o o to v v v along this path, so that ∂ p v = 1 v − 1 o \partial p_v=\mathbf{1}_v-\mathbf{1}_o ∂ p v = 1 v − 1 o , where 1 w ∈ R V \mathbf{1}_w\in\mathbb{R}^{V} 1 w ∈ R V is the indicator of a vertex w w w . For the non-tree edges e j = ( u j , v j ) e_j=(u_j,v_j) e j = ( u j , v j ) , let c j = p u j + e j − p v j c_j=p_{u_j}+e_j-p_{v_j} c j = p u j + e j − p v j for j ∈ [ r ] j\in[r] j ∈ [ r ] . Notice that { c j } j = 1 r \{c_j\}_{j=1}^{r} { c j } j = 1 r are linearly independent. Indeed, suppose that ∑ j = 1 r a j c j = 0 \sum_{j=1}^{r}a_jc_j=0 ∑ j = 1 r a j c j = 0 where the coefficients a j a_j a j are not all zero. Only c j c_j c j has 1 1 1 at the coordinate of the non-tree edge e j e_j e j , and all other vectors have 0 0 0 , and so we must have a j = 0 a_j=0 a j = 0 for each j ∈ [ r ] j\in[r] j ∈ [ r ] . Therefore, the set { c j } j = 1 r \{c_j\}_{j=1}^{r} { c j } j = 1 r forms a basis of Z Z Z because c j ∈ Z c_j\in Z c j ∈ Z and dim Z = r \dim Z=r dim Z = r . Let C a = ∑ j a j c j C_a=\sum_j a_jc_j C a = ∑ j a j c j . The vectors p v p_v p v , v ≠ o v\ne o v = o , complement Z Z Z , since their images 1 v − 1 o \mathbf{1}_v-\mathbf{1}_o 1 v − 1 o under ∂ \partial ∂ form a basis of the zero-sum vectors in R V \mathbb{R}^{V} R V . Hence, a linear map on R E \mathbb{R}^{E} R E is determined by its values on Z Z Z and on the p v p_v p v .
A norm in which the cycle space is Euclidean
We define a new norm on W = R E W = \mathbb{R}^{E} W = R E that interpolates between the ℓ 2 \ell_2 ℓ 2 and the ℓ 1 \ell_1 ℓ 1 norm:
∥ x ∥ W = inf z ∈ Z ( ∥ z ∥ 2 + ∥ x − z ∥ 1 g ) . (21) \|x\|_{W} = \inf_{z \in Z}\left(\|z\|_2 + \frac{\|x-z\|_1}{\sqrt{g}}\right). \tag*{(21)} ∥ x ∥ W = z ∈ Z inf ( ∥ z ∥ 2 + g ∥ x − z ∥ 1 ) . (21)
Its unit ball is B W = conv ( B Z , 2 ∪ g B ℓ 1 E ) B_W = \operatorname{conv}(B_{Z,2} \cup\sqrt{g}B_{\ell_1}^{E}) B W = conv ( B Z , 2 ∪ g B ℓ 1 E ) . In particular, z = 0 z = 0 z = 0 gives ∥ e ∥ W ≤ 1 / g \|e\|_W \le1/\sqrt{g} ∥ e ∥ W ≤ 1/ g for every edge. Hence, ∥ p v ∥ W ≤ D / g \|p_v\|_W \le D/\sqrt{g} ∥ p v ∥ W ≤ D / g .
Each cycle has at least g g g edges, so ∥ z ∥ 1 ≥ g ∥ z ∥ 2 \|z\|_1 \ge\sqrt{g}\|z\|_2 ∥ z ∥ 1 ≥ g ∥ z ∥ 2 for z ∈ Z z \in Z z ∈ Z . Suppose z = ∑ i a i γ i z = \sum_i a_i\gamma_i z = ∑ i a i γ i with a i > 0 a_i > 0 a i > 0 is a decomposition into signed cycles, such that each γ i \gamma_i γ i agrees in sign with z z z on its support. Then
∥ z ∥ 2 ≤ ∑ i a i ∣ γ i ∣ ≤ g − 1 / 2 ∑ i a i ∣ γ i ∣ = g − 1 / 2 ∥ z ∥ 1 . \|z\|_2 \le\sum_i a_i\sqrt{|\gamma_i|} \le g^{-1/2}\sum_i a_i|\gamma_i| = g^{-1/2}\|z\|_1. ∥ z ∥ 2 ≤ i ∑ a i ∣ γ i ∣ ≤ g − 1/2 i ∑ a i ∣ γ i ∣ = g − 1/2 ∥ z ∥ 1 .
For x ∈ Z x \in Z x ∈ Z , applying this to x − z x-z x − z in (21) gives ∥ x ∥ W ≥ ∥ x ∥ 2 \|x\|_W \ge\|x\|_2 ∥ x ∥ W ≥ ∥ x ∥ 2 , and z = x z=x z = x gives equality. Thus, Z Z Z inherits the Euclidean norm from W W W .
The projection constant in the graph construction
Let P \mathcal{P} P be the family of linear maps P : R E → Z P:\mathbb{R}^{E}\to Z P : R E → Z that fix Z Z Z , so that λ ( Z , W ) = min P ∈ P ∥ P ∥ W → Z \lambda(Z,W)=\min_{P\in\mathcal{P}}\|P\|_{W\to Z} λ ( Z , W ) = min P ∈ P ∥ P ∥ W → Z . The norm of each P ∈ P P\in\mathcal{P} P ∈ P is determined by the images of the edges.
Lemma 5.2. For every P ∈ P P\in\mathcal{P} P ∈ P , we have ∥ P ∥ W → Z = max { 1 , g max e ∈ E ∥ P e ∥ 2 } \|P\|_{W\to Z}=\max\{1,\sqrt{g}\max_{e\in E}\|Pe\|_2\} ∥ P ∥ W → Z = max { 1 , g max e ∈ E ∥ P e ∥ 2 } . As a consequence, max { 1 , g r / m } ≤ λ ( Z , W ) ≤ g \max\{1,\sqrt{gr/m}\}\le\lambda(Z,W)\le\sqrt{g} max { 1 , g r / m } ≤ λ ( Z , W ) ≤ g .
Proof. Because P w ∈ Z Pw\in Z P w ∈ Z for w ∈ W w\in W w ∈ W , we have ∥ P w ∥ W = ∥ P w ∥ 2 \|Pw\|_W=\|Pw\|_2 ∥ P w ∥ W = ∥ P w ∥ 2 , and it suffices to maximize ∥ P w ∥ 2 \|Pw\|_2 ∥ P w ∥ 2 over the two sets whose convex hull is B W B_W B W . On the ball B Z , 2 B_{Z,2} B Z , 2 the maximum is one, since P P P fixes Z Z Z . On the ball g B ℓ 1 E \sqrt{g}B_{\ell_1}^{E} g B ℓ 1 E it is g max e ∈ E ∥ P e ∥ 2 \sqrt{g}\max_{e\in E}\|Pe\|_2 g max e ∈ E ∥ P e ∥ 2 , attained at a vertex ± g e \pm\sqrt{g}e ± g e . Combining these yields the first equality.
Next, P P P fixes an orthonormal basis of Z Z Z . Extending it to an orthonormal basis of R E \mathbb{R}^{E} R E gives r ≤ ∥ P ∥ H S 2 = ∑ e ∈ E ∥ P e ∥ 2 2 ≤ m max e ∥ P e ∥ 2 2 r\le\|P\|_{\mathrm{HS}}^2=\sum_{e\in E}\|Pe\|_2^2\le m\max_e\|Pe\|_2^2 r ≤ ∥ P ∥ HS 2 = ∑ e ∈ E ∥ P e ∥ 2 2 ≤ m max e ∥ P e ∥ 2 2 . Together with the first equality, this proves the lower bound. For the upper bound, the Euclidean orthogonal projection Q Q Q onto Z Z Z satisfies ∥ Q e ∥ 2 ≤ ∥ e ∥ 2 = 1 \|Qe\|_2\le\|e\|_2=1 ∥ Q e ∥ 2 ≤ ∥ e ∥ 2 = 1 for every edge, so ∥ Q ∥ W → Z ≤ g \|Q\|_{W\to Z}\le\sqrt{g} ∥ Q ∥ W → Z ≤ g .
:::
From Projection to Lipschitz Extension
For an integer N ≥ 1 N\ge1 N ≥ 1 , let B N = { − N , … , N } r B_N=\{-N,\ldots,N\}^{r} B N = { − N , … , N } r . The cell of a vertex v v v is U v , N = { p v + C a : a ∈ B N } ⊂ W U_{v,N}=\{p_v+Ca:a\in B_N\}\subset W U v , N = { p v + C a : a ∈ B N } ⊂ W , and their union is X N = ⋃ v ∈ V U v , N X_N=\bigcup_{v\in V}U_{v,N} X N = ⋃ v ∈ V U v , N . The cells are disjoint because ∂ p v = 1 v − 1 o \partial p_v=1_v-1_o ∂ p v = 1 v − 1 o , and they are Euclidean because their differences lie in Z Z Z . Let A N = U o , N A_N=U_{o,N} A N = U o , N and ϕ ( C a ) = C a \phi(Ca)=Ca ϕ ( C a ) = C a . Let e ( A N ) e(A_N) e ( A N ) be the optimal Lipschitz constant of an extension of ϕ \phi ϕ to X N X_N X N .
Theorem 5.3 (From projection constant to Lipschitz constant). Let β N = 2 D / ( 2 N − 1 ) \beta_N=2D/(2N-1) β N = 2 D / ( 2 N − 1 ) . Then
max { 1 , λ ( Z , W ) 1 + β N } ≤ e ( A N ) ≤ λ ( Z , W ) . (22) \max\left\{1,\frac{\lambda(Z,W)}{1+\beta_N}\right\}\le e(A_N)\le\lambda(Z,W). \tag*{(22)} max { 1 , 1 + β N λ ( Z , W ) } ≤ e ( A N ) ≤ λ ( Z , W ) . (22) In particular, N ≥ D + 1 N\ge D+1 N ≥ D + 1 gives e ( A N ) ≥ 1 2 g r / m e(A_N)\ge\frac{1}{2}\sqrt{gr/m} e ( A N ) ≥ 2 1 g r / m .
Proof. Any linear projection P ∈ P P\in\mathcal{P} P ∈ P is a ∥ P ∥ W → Z \|P\|_{W\to Z} ∥ P ∥ W → Z -Lipschitz extension of ϕ \phi ϕ . Taking the infimum gives the upper bound e ( A N ) ≤ λ ( Z , W ) e(A_N)\le\lambda(Z,W) e ( A N ) ≤ λ ( Z , W ) .
Now, take an L L L -Lipschitz extension F : X N → Z F:X_N\to Z F : X N → Z that coincides with ϕ \phi ϕ on A N A_N A N . For a ∈ B N a\in B_N a ∈ B N , let h v ( a ) = F ( p v + C a ) − F ( C a ) = F ( p v + C a ) − C a h_v(a)=F(p_v+Ca)-F(Ca)=F(p_v+Ca)-Ca h v ( a ) = F ( p v + C a ) − F ( C a ) = F ( p v + C a ) − C a be the offset of the cell U v , N U_{v,N} U v , N . Let
b v = ∣ J N ∣ − 1 ∑ a ∈ J N h v ( a ) b_v=|J_N|^{-1}\sum_{a\in J_N}h_v(a) b v = ∣ J N ∣ − 1 a ∈ J N ∑ h v ( a )
be the average of h v h_v h v over the inner box J N = { − N + 1 , … , N − 1 } r J_N=\{-N+1,\ldots,N-1\}^{r} J N = { − N + 1 , … , N − 1 } r . Since F ( C a ) = C a F(Ca)=Ca F ( C a ) = C a and ∥ p v ∥ W ≤ D / g \|p_v\|_W\le D/\sqrt{g} ∥ p v ∥ W ≤ D / g , notice that h o = 0 h_o=0 h o = 0 and ∥ h v ( a ) ∥ 2 ≤ L D / g \|h_v(a)\|_2\le LD/\sqrt{g} ∥ h v ( a ) ∥ 2 ≤ L D / g . For an edge e = ( u , v ) e=(u,v) e = ( u , v ) let q e = p u + e − p v q_e=p_u+e-p_v q e = p u + e − p v . It is zero on tree edges and equals c j c_j c j on e j e_j e j . When e e e is a tree edge, let η e = 0 \eta_e=0 η e = 0 . When e = e j e=e_j e = e j is a non-tree edge, let η e ∈ Z r \eta_e\in\mathbb{Z}^r η e ∈ Z r be the indicator vector with 1 at position j j j . This way we have C η e = q e C\eta_e=q_e C η e = q e .
Now because the p v p_v p v for v ≠ o v\ne o v = o complement Z Z Z , we can define the linear map P P P , such that P ∣ Z = I Z P|_Z=I_Z P ∣ Z = I Z and P p v = b v Pp_v=b_v P p v = b v . Then P e = q e + b v − b u Pe=q_e+b_v-b_u P e = q e + b v − b u for the edge e = ( u , v ) e=(u,v) e = ( u , v ) . For a ∈ J N a\in J_N a ∈ J N , the two points p u + C a p_u+Ca p u + C a and p v + C ( a + η e ) p_v+C(a+\eta_e) p v + C ( a + η e ) differ by e e e . Then,
F ( p v + C ( a + η e ) ) − F ( p u + C a ) = h v ( a + η e ) − h u ( a ) + C η e , F(p_v+C(a+\eta_e))-F(p_u+Ca)=h_v(a+\eta_e)-h_u(a)+C\eta_e, F ( p v + C ( a + η e )) − F ( p u + C a ) = h v ( a + η e ) − h u ( a ) + C η e ,
and
∥ F ( p v + C ( a + η e ) ) − F ( p u + C a ) ∥ 2 ≤ L ∥ p v + C η e − p u ∥ W = L ∥ p v + q e − p u ∥ W = L ∥ e ∥ W . \begin{aligned}
\|F(p_v+C(a+\eta_e))-F(p_u+Ca)\|_2&\le L\|p_v+C\eta_e-p_u\|_W\\
&=L\|p_v+q_e-p_u\|_W\\
&=L\|e\|_W.
\end{aligned} ∥ F ( p v + C ( a + η e )) − F ( p u + C a ) ∥ 2 ≤ L ∥ p v + C η e − p u ∥ W = L ∥ p v + q e − p u ∥ W = L ∥ e ∥ W .
Therefore, ∥ q e + h v ( a + η e ) − h u ( a ) ∥ 2 ≤ L / g \|q_e+h_v(a+\eta_e)-h_u(a)\|_2\le L/\sqrt{g} ∥ q e + h v ( a + η e ) − h u ( a ) ∥ 2 ≤ L / g for every a ∈ J N a\in J_N a ∈ J N because ∥ e ∥ W ≤ 1 / g \|e\|_W\le1/\sqrt{g} ∥ e ∥ W ≤ 1/ g . Averaging this inequality over a ∈ J N a\in J_N a ∈ J N for fixed edge e e e , and using the triangle inequality gives ∥ q e + b ~ v − b u ∥ 2 ≤ L / g \|q_e+\widetilde b_v-b_u\|_2\le L/\sqrt{g} ∥ q e + b v − b u ∥ 2 ≤ L / g , where b ~ v = ∣ J N ∣ − 1 ∑ a ∈ J N h v ( a + η e ) \widetilde b_v=|J_N|^{-1}\sum_{a\in J_N}h_v(a+\eta_e) b v = ∣ J N ∣ − 1 ∑ a ∈ J N h v ( a + η e ) is the average of h v h_v h v over the shifted box J N + η e J_N+\eta_e J N + η e . If e e e is a tree edge, then η e = 0 \eta_e=0 η e = 0 and b ~ v = b v \widetilde b_v=b_v b v = b v . Otherwise, η e \eta_e η e is the indicator of position j j j , and J N + η e J_N+\eta_e J N + η e gains { a ∈ J N + η e ∣ a j = N } \{a\in J_N+\eta_e\mid a_j=N\} { a ∈ J N + η e ∣ a j = N } and loses { a ∈ J N ∣ a j = − N + 1 } \{a\in J_N\mid a_j=-N+1\} { a ∈ J N ∣ a j = − N + 1 } . All other terms of the two sums cancel out, and each set has ( 2 N − 1 ) r − 1 = ∣ J N ∣ / ( 2 N − 1 ) (2N-1)^{r-1}=|J_N|/(2N-1) ( 2 N − 1 ) r − 1 = ∣ J N ∣/ ( 2 N − 1 ) points, so the offset bound gives ∥ b ~ v − b v ∥ 2 ≤ 2 2 N − 1 ⋅ L g = β N L g \|\widetilde b_v-b_v\|_2\le\frac{2}{2N-1}\cdot\frac{L}{\sqrt{g}}=\frac{\beta_NL}{\sqrt{g}} ∥ b v − b v ∥ 2 ≤ 2 N − 1 2 ⋅ g L = g β N L . Hence, we get
∥ P e ∥ 2 = ∥ q e + b v − b u ∥ 2 ≤ ∥ q e + b ~ v − b u ∥ 2 + ∥ b ~ v − b v ∥ 2 ≤ L g ( 1 + β N ) . \|Pe\|_2=\|q_e+b_v-b_u\|_2\le\|q_e+\widetilde b_v-b_u\|_2+\|\widetilde b_v-b_v\|_2\le\frac{L}{\sqrt{g}}(1+\beta_N). ∥ P e ∥ 2 = ∥ q e + b v − b u ∥ 2 ≤ ∥ q e + b v − b u ∥ 2 + ∥ b v − b v ∥ 2 ≤ g L ( 1 + β N ) .
We have L ≥ 1 L\ge1 L ≥ 1 because F F F is the identity on A N A_N A N , so Lemma 5.2 yields λ ( Z , W ) ≤ ∥ P ∥ W → Z ≤ ( 1 + β N ) L \lambda(Z,W)\le\|P\|_{W\to Z}\le(1+\beta_N)L λ ( Z , W ) ≤ ∥ P ∥ W → Z ≤ ( 1 + β N ) L . Taking the infimum over F F F gives the lower bound. □
Corollary 5.4. For every k ≥ 2 k\ge2 k ≥ 2 there is a metric space that is a union of k k k Euclidean sets with e max = Ω ( log k ) e_{\max}=\Omega(\sqrt{\log k}) e m a x = Ω ( log k ) .
Proof. By [9 ] , there is a 3-regular graph G G G with n ≤ k n\le k n ≤ k vertices and girth g = Ω ( log k ) g=\Omega(\log k) g = Ω ( log k ) . Then m = 3 n / 2 m=3n/2 m = 3 n /2 and r = m − n + 1 > m / 3 r=m-n+1>m/3 r = m − n + 1 > m /3 , so g r / m = Ω ( log k ) gr/m=\Omega(\log k) g r / m = Ω ( log k ) . Let X = X N X=X_N X = X N for N ≥ D + 1 N\ge D+1 N ≥ D + 1 , which is a union of the n n n Euclidean sets U v , N U_{v,N} U v , N , and adjoin k − n k-n k − n further points, each forming its own set, at distance Δ \Delta Δ from every other point, where Δ \Delta Δ is the diameter of X N X_N X N . Restricting an extension from X X X to X N X_N X N does not increase its Lipschitz constant, so Theorem 5.3 gives e max ≥ e ( A N ) ≥ 1 2 g r / m = Ω ( log k ) e_{\max}\ge e(A_N)\ge\frac{1}{2}\sqrt{gr/m}=\Omega(\sqrt{\log k}) e m a x ≥ e ( A N ) ≥ 2 1 g r / m = Ω ( log k ) . □
Lower Bound for 1-Strongly Euclidean Sets
We show that the Lipschitz extension constant can grow even when every set in the union is 1-strongly Euclidean. We use a variant of the twisted cube of Naor and Rabani [24 Section 4.1] , with cross distances chosen so that one layer remains Euclidean after adjoining any one point.
Theorem 6.1. For every integer k ≥ 2 k\ge2 k ≥ 2 , there is a finite metric space X = U 1 ∪ ⋯ ∪ U k X=U_1\cup\cdots\cup U_k X = U 1 ∪ ⋯ ∪ U k such that each U i U_i U i is 1-strongly Euclidean in X X X and
e max ≥ 1 2 ( log 2 k ) 1 / 4 = Ω ( log k 4 ) . e_{\max}\ge\frac{1}{2}(\log_2 k)^{1/4}=\Omega(\sqrt[4]{\log k}). e m a x ≥ 2 1 ( log 2 k ) 1/4 = Ω ( 4 log k ) .
The two-layer cube
Let Q m = { 0 , 1 } m Q_m=\{0,1\}^m Q m = { 0 , 1 } m with Hamming distance Ham \operatorname{Ham} Ham , and take two copies A = { a x : x ∈ Q m } A=\{a_x:x\in Q_m\} A = { a x : x ∈ Q m } and B = { b x : x ∈ Q m } B=\{b_x:x\in Q_m\} B = { b x : x ∈ Q m } . For r > 0 r>0 r > 0 , let s = r 2 + 1 − r = 1 / ( r 2 + 1 + r ) s=\sqrt{r^2+1}-r=1/(\sqrt{r^2+1}+r) s = r 2 + 1 − r = 1/ ( r 2 + 1 + r ) . On X m = A ∪ B X_m=A\cup B X m = A ∪ B , define
d ( a x , a y ) = Ham ( x , y ) , d ( a x , b y ) = r 2 + Ham ( x , y ) , d ( b x , b y ) = min { s Ham ( x , y ) , Ham ( x , y ) } . (23) \begin{aligned}
d(a_x,a_y)&=\sqrt{\operatorname{Ham}(x,y)},\\
d(a_x,b_y)&=\sqrt{r^2+\operatorname{Ham}(x,y)},\\
d(b_x,b_y)&=\min\{s\operatorname{Ham}(x,y),\sqrt{\operatorname{Ham}(x,y)}\}. \tag*{(23)}
\end{aligned} d ( a x , a y ) d ( a x , b y ) d ( b x , b y ) = Ham ( x , y ) , = r 2 + Ham ( x , y ) , = min { s Ham ( x , y ) , Ham ( x , y ) } . (23)
Thus, A A A is the usual Euclidean cube, while distances in B B B are shortened at small Hamming distances.
Lemma 6.2. For every integer m ≥ 1 m\ge1 m ≥ 1 and r > 0 r>0 r > 0 , (23) defines a metric on X m X_m X m , and A A A is 1-strongly Euclidean in X m X_m X m .
Proof. For each fixed y ∈ Q m y\in Q_m y ∈ Q m , the map a x ↦ ( x , 0 ) a_x\mapsto(x,0) a x ↦ ( x , 0 ) and b y ↦ ( y , r ) b_y\mapsto(y,r) b y ↦ ( y , r ) realizes A ∪ { b y } A\cup\{b_y\} A ∪ { b y } isometrically in R m + 1 \mathbb{R}^{m+1} R m + 1 . This proves the one-point Euclidean property and all triangle inequalities with at least two vertices in A A A .
The function t ↦ min { s t , t } t\mapsto\min\{st,\sqrt{t}\} t ↦ min { s t , t } is increasing, concave and vanishes at zero, so it is subadditive and its composition with Ham \operatorname{Ham} Ham is a metric on B B B . It remains to check triangles with one vertex in A A A and two in B B B .
For every integer q ≥ 0 q\ge0 q ≥ 0 , concavity of t ↦ r 2 + t − r t\mapsto\sqrt{r^2+t}-r t ↦ r 2 + t − r and the definition of s s s give
r 2 + q − r ≤ min { s q , q } . \sqrt{r^2+q}-r\le\min\{sq,\sqrt{q}\}. r 2 + q − r ≤ min { s q , q } .
Fix a z , b x , b y a_z,b_x,b_y a z , b x , b y , and write p = Ham ( z , x ) p=\operatorname{Ham}(z,x) p = Ham ( z , x ) , p ′ = Ham ( z , y ) p'=\operatorname{Ham}(z,y) p ′ = Ham ( z , y ) and q = Ham ( x , y ) q=\operatorname{Ham}(x,y) q = Ham ( x , y ) . Since ∣ p − p ′ ∣ ≤ q |p-p'|\le q ∣ p − p ′ ∣ ≤ q , subadditivity of the same function gives
∣ d ( a z , b x ) − d ( a z , b y ) ∣ ≤ r 2 + ∣ p − p ′ ∣ − r ≤ r 2 + q − r ≤ d ( b x , b y ) . \left|d(a_z,b_x)-d(a_z,b_y)\right|
\le\sqrt{r^2+|p-p'|}-r
\le\sqrt{r^2+q}-r
\le d(b_x,b_y). ∣ d ( a z , b x ) − d ( a z , b y ) ∣ ≤ r 2 + ∣ p − p ′ ∣ − r ≤ r 2 + q − r ≤ d ( b x , b y ) .
The remaining inequality follows from
d ( b x , b y ) ≤ q ≤ p + p ′ ≤ r 2 + p + r 2 + p ′ = d ( a z , b x ) + d ( a z , b y ) . d(b_x,b_y)\le\sqrt{q}\le\sqrt{p}+\sqrt{p'}\le\sqrt{r^2+p}+\sqrt{r^2+p'}=d(a_z,b_x)+d(a_z,b_y). d ( b x , b y ) ≤ q ≤ p + p ′ ≤ r 2 + p + r 2 + p ′ = d ( a z , b x ) + d ( a z , b y ) .
Extending the cube
An extension must keep each b x b_x b x close to the prescribed image of a x a_x a x , while neighboring points in B B B must have nearby images. Averaging over the cube shows that these two requirements force a large Lipschitz constant.
Lemma 6.3. For the metric in (23) , we have
e ( A ) ≥ m 2 r + s m . e(A)\ge\frac{\sqrt{m}}{2r+s\sqrt{m}}. e ( A ) ≥ 2 r + s m m .
Proof. Let F : X m → ℓ 2 F:X_m\to\ell_2 F : X m → ℓ 2 be an L L L -Lipschitz extension of the isometry f ( a x ) = x f(a_x)=x f ( a x ) = x , and let u x = F ( b x ) u_x=F(b_x) u x = F ( b x ) and c = ( 1 / 2 , … , 1 / 2 ) ∈ R m c=(1/2,\ldots,1/2)\in\mathbb{R}^m c = ( 1/2 , … , 1/2 ) ∈ R m . If x ( j ) x^{(j)} x ( j ) is obtained from x x x by flipping coordinate j j j , then ∥ u x − x ∥ ≤ L r \|u_x-x\|\le Lr ∥ u x − x ∥ ≤ L r and ∥ u x − u x ( j ) ∥ ≤ L s \|u_x-u_{x^{(j)}}\|\le Ls ∥ u x − u x ( j ) ∥ ≤ L s .
Let x x x be uniform on Q m Q_m Q m . Pairing x x x with x ( j ) x^{(j)} x ( j ) gives, for every j ∈ [ m ] j\in[m] j ∈ [ m ] ,
E x [ ( x j − 1 / 2 ) ( ( u x ) j − 1 / 2 ) ] = 1 2 E x [ ( x j − 1 / 2 ) ( u x − u x ( j ) ) j ] ≤ L s 4 . \mathbb{E}_x\left[(x_j-1/2)((u_x)_j-1/2)\right]
=\frac{1}{2}\mathbb{E}_x\left[(x_j-1/2)(u_x-u_{x^{(j)}})_j\right]
\le\frac{Ls}{4}. E x [ ( x j − 1/2 ) (( u x ) j − 1/2 ) ] = 2 1 E x [ ( x j − 1/2 ) ( u x − u x ( j ) ) j ] ≤ 4 L s .
Summing over j j j and using ∥ x − c ∥ = m / 2 \|x-c\|=\sqrt{m}/2 ∥ x − c ∥ = m /2 , we obtain
m 4 = E x [ ⟨ x − u x , x − c ⟩ ] + E x [ ⟨ u x − c , x − c ⟩ ] ≤ L r m 2 + m L s 4 . \frac{m}{4}
=\mathbb{E}_x[\langle x-u_x,x-c\rangle]
+\mathbb{E}_x[\langle u_x-c,x-c\rangle]
\le\frac{Lr\sqrt{m}}{2}+\frac{mLs}{4}. 4 m = E x [⟨ x − u x , x − c ⟩] + E x [⟨ u x − c , x − c ⟩] ≤ 2 L r m + 4 m L s .
Rearranging gives L ≥ m / ( 2 r + s m ) L\ge\sqrt{m}/(2r+s\sqrt{m}) L ≥ m / ( 2 r + s m ) . Since f f f is an isometry, the same lower bound holds for e ( A ) e(A) e ( A ) .
Proof of Theorem 6.1
Proof. Let m = ⌈ log 2 k ⌉ m=\lceil\log_{2} k\rceil m = ⌈ log 2 k ⌉ and r = m 1 / 4 / 2 r=m^{1/4}/2 r = m 1/4 /2 . Then s ≤ 1 / ( 2 r ) = m − 1 / 4 s\le1/(2r)=m^{-1/4} s ≤ 1/ ( 2 r ) = m − 1/4 , so Lemma 6.3 gives
e ( A ) ≥ m 2 r + s m ≥ m 1 / 4 2 ≥ 1 2 ( log 2 k ) 1 / 4 . e(A)\ge\frac{\sqrt{m}}{2r+s\sqrt{m}}\ge\frac{m^{1/4}}{2}\ge\frac{1}{2}(\log_{2}k)^{1/4}. e ( A ) ≥ 2 r + s m m ≥ 2 m 1/4 ≥ 2 1 ( log 2 k ) 1/4 .
Keep A A A as one set and partition B B B into pairs. By Lemma 6.2, A A A is 1-strongly Euclidean in X m X_m X m . Each pair is also 1-strongly Euclidean, because adjoining one point gives a metric space on at most three points, which is a Euclidean triangle.
This partition has 1 + 2 m − 1 ≤ k 1+2^{m-1}\le k 1 + 2 m − 1 ≤ k sets. Since k ≤ 2 m k\le2^m k ≤ 2 m , splitting some pairs into singletons gives exactly k k k sets while keeping A A A whole. All sets remain 1-strongly Euclidean, and e max ≥ e ( A ) e_{\max}\ge e(A) e m a x ≥ e ( A ) proves the theorem.
Unions of sets with constant distortion
We now assume that each U i U_i U i embeds into ℓ 2 \ell_2 ℓ 2 with distortion C i ≥ 1 C_i\ge1 C i ≥ 1 , so for each U i U_i U i there is a map f i : U i → ℓ 2 f_i:U_i\to\ell_2 f i : U i → ℓ 2 , such that d ( x , y ) ≤ ∥ f i ( x ) − f i ( y ) ∥ ≤ C i d ( x , y ) d(x,y)\le\lVert f_i(x)-f_i(y)\rVert\le C_i d(x,y) d ( x , y ) ≤ ∥ f i ( x ) − f i ( y )∥ ≤ C i d ( x , y ) for every x , y ∈ U i x,y\in U_i x , y ∈ U i . Directly applying [21 Theorem 1.1] for k = 2 k=2 k = 2 gives distortion at most 2 O ( k ) ∏ i = 1 k C i 2^{O(k)}\prod_{i=1}^{k}C_i 2 O ( k ) ∏ i = 1 k C i , which is exponential in k k k even when all C i C_i C i are equal to one.
For c ≥ 1 c\ge1 c ≥ 1 , let η ( c ) \eta(c) η ( c ) be the unique root η ≥ c − 1 \eta\ge c-1 η ≥ c − 1 of ( ( 1 + η ) 2 − c 2 ) ( ( 1 + η / 2 ) 2 − 2 ) = ( π c 2 ) / 2 ((1+\eta)^2-c^2)((1+\eta/2)^2-2)=(\pi c^2)/2 (( 1 + η ) 2 − c 2 ) (( 1 + η /2 ) 2 − 2 ) = ( π c 2 ) /2 . For η ≥ c − 1 \eta\ge c-1 η ≥ c − 1 , the left side is at most 0 while ( 1 + η / 2 ) 2 ≤ 2 (1+\eta/2)^2\le2 ( 1 + η /2 ) 2 ≤ 2 and afterwards increases to infinity, so the root exists and is unique. It is increasing in c c c , and η ( 1 ) = η 0 \eta(1)=\eta_0 η ( 1 ) = η 0 by (16) . Moreover, η ( c ) ≤ c − 1 + π / c \eta(c)\le c-1+\pi/c η ( c ) ≤ c − 1 + π / c , because at this value of η \eta η the two factors on the left are at least 2 π 2\pi 2 π and c 2 / 4 c^2/4 c 2 /4 .
Next, we establish the following variant of Lemma 4.6 , in which the removed set is only required to embed with distortion c c c .
Lemma 7.1. Let Y = A ⊔ B Y=A\mathbin{\sqcup}B Y = A ⊔ B be a finite pseudometric space with A ≠ ∅ A\ne\varnothing A = ∅ , where B B B embeds into ℓ 2 \ell_2 ℓ 2 with distortion c ≥ 1 c\ge1 c ≥ 1 . For a weighted graph G G G on Y Y Y , let I I I and J J J be the internal and crossing energies of B B B as in Lemma 4.6 . Then
E G A ( d 2 ) ≤ E G ( d 2 ) + ( c 2 − 1 ) I + J + c 2 π I J . (24) \mathcal{E}_{G_A}(d^2)\le\mathcal{E}_{G}(d^2)+(c^2-1)I+J+c\sqrt{2\pi IJ}. \tag*{(24)} E G A ( d 2 ) ≤ E G ( d 2 ) + ( c 2 − 1 ) I + J + c 2 π I J . (24)
Proof. Choose z b ∈ ℓ 2 z_b\in\ell_2 z b ∈ ℓ 2 for b ∈ B b\in B b ∈ B with d ( b , b ′ ) ≤ ∥ z b − z b ′ ∥ ≤ c d ( b , b ′ ) d(b,b')\le\lVert z_b-z_{b'}\rVert\le c\,d(b,b') d ( b , b ′ ) ≤ ∥ z b − z b ′ ∥ ≤ c d ( b , b ′ ) , and let I z = E G [ B ] ( z ) ≤ c 2 I I_z=\mathcal{E}_{G[B]}(z)\le c^2I I z = E G [ B ] ( z ) ≤ c 2 I . The proof of Lemma 4.6 applies with I z I_z I z in place of I I I . Indeed, the route a p → b p → b q → a q a_p\to b_p\to b_q\to a_q a p → b p → b q → a q still gives d ( a p , a q ) ≤ u p + ∥ z b p − z b q ∥ + u q d(a_p,a_q)\le u_p+\lVert z_{b_p}-z_{b_q}\rVert+u_q d ( a p , a q ) ≤ u p + ∥ z b p − z b q ∥ + u q because z z z does not contract distances, the middle term is bounded by 2 t I z 2tI_z 2 t I z in (8) , and in the reduction to the case s b > 0 s_b>0 s b > 0 the identity (6) applied to z z z gives I z ′ ≤ I z I'_z\le I_z I z ′ ≤ I z . Hence,
E G A ( d 2 ) ≤ E G [ A ] ( d 2 ) + 2 J + 2 π I z J + I z ≤ E G ( d 2 ) + ( c 2 − 1 ) I + J + c 2 π I J . \mathcal{E}_{G_A}(d^2)\le\mathcal{E}_{G[A]}(d^2)+2J+\sqrt{2\pi I_zJ}+I_z\le\mathcal{E}_{G}(d^2)+(c^2-1)I+J+c\sqrt{2\pi IJ}. E G A ( d 2 ) ≤ E G [ A ] ( d 2 ) + 2 J + 2 π I z J + I z ≤ E G ( d 2 ) + ( c 2 − 1 ) I + J + c 2 π I J .
Proof of Theorem 1.7. Fix i i i , keep U i U_i U i , and partition its complement into at most k − 1 k-1 k − 1 subsets of the other sets, as in the proof of Lemma 1.5, so that subsets inherit the distortion bounds. For a graph with no edges inside A A A , Lemma 7.1 reads E G A ( d 2 ) ≤ c 2 I + 2 J + c 2 π I J \mathcal{E}_{G_A}(d^2)\le c^2I+2J+c\sqrt{2\pi IJ} E G A ( d 2 ) ≤ c 2 I + 2 J + c 2 π I J . Hence, Lemma 4.7 applies with t = 1 + η t=1+\eta t = 1 + η and s = 1 + η / 2 s=1+\eta/2 s = 1 + η /2 whenever t ≥ c t\ge c t ≥ c , s ≥ 2 s\ge\sqrt{2} s ≥ 2 and ( t 2 − c 2 ) ( s 2 − 2 ) ≥ π c 2 / 2 (t^2-c^2)(s^2-2)\ge\pi c^2/2 ( t 2 − c 2 ) ( s 2 − 2 ) ≥ π c 2 /2 , that is, whenever η ≥ η ( c ) \eta\ge\eta(c) η ≥ η ( c ) . Embedding with distortion at most C max C_{\max} C m a x is preserved by restriction, so Lemma 4.8 with η = η ( C max ) > 1 \eta=\eta(C_{\max})>1 η = η ( C m a x ) > 1 extends every 1-Lipschitz map on U i U_i U i to X X X with Lipschitz constant at most k η ( C max ) k^{\eta(C_{\max})} k η ( C m a x ) . Applied to f i / C i f_i/C_i f i / C i , this extends f i f_i f i to f ~ i : X → ℓ 2 \widetilde f_i:X\to\ell_2 f i : X → ℓ 2 with Lip ( f ~ i ) ≤ C i k η ( C max ) \operatorname{Lip}(\widetilde f_i)\le C_i k^{\eta(C_{\max})} Lip ( f i ) ≤ C i k η ( C m a x ) . Since f i f_i f i does not contract distances, Lemma 3.1 gives a 3 2 C max k η ( C max ) 3\sqrt{2}C_{\max}k^{\eta(C_{\max})} 3 2 C m a x k η ( C m a x ) -terminal embedding for every U i U_i U i , and Theorem 2.7 gives the claimed distortion.
Generative AI Statement. The generalization of Bourgain’s Theorem for 2-strongly Euclidean sets (Theorem 1.1 ) was obtained by the authors, without any AI usage. After prompting with the proof of Theorem 1.1 , GPT-6 Astra was able to obtain the presented results for unions of Euclidean (Theorem 1.2 ) and 1-strongly Euclidean sets (Theorem 1.3 ). The exposition was significantly rewritten by the authors. In particular, the authors uncovered the connections to Lipschitz extendability and Euclidean metric sparsifiers. The authors take full responsibility for the presented results.
Proof of the Euclidean sparsifier–extension equivalence
We prove Theorem 4.2 . Fix a partition U 1 , … , U k U_1,\ldots,U_k U 1 , … , U k of X X X , let D \mathcal{D} D be the class of pseudometrics on X X X that are Euclidean on every U j U_j U j , and set A = ⋃ i U i A=\bigcup_i U_i A = ⋃ i U i . Let Q = Q ( D , A ) Q=Q(\mathcal{D},A) Q = Q ( D , A ) .
From sparsifiers to extensions. Fix Q ′ > Q Q'>Q Q ′ > Q . Let d ∈ D d\in\mathcal{D} d ∈ D , and let f : A → ℓ 2 f:A\to\ell_2 f : A → ℓ 2 be 1-Lipschitz. The extension claim is immediate when X X X has one point, so assume ∣ X ∣ ≥ 2 |X|\ge2 ∣ X ∣ ≥ 2 . Since A A A is finite, we may take f f f to have values in R m \mathbb{R}^m R m .
Fix a weighted graph G G G on X X X , and let H H H be a Euclidean metric sparsifier for G G G of quality at most Q ′ Q' Q ′ . For a scalar map u : A → R u:A\to\mathbb{R} u : A → R , put ρ u ( a , b ) = ∣ u ( a ) − u ( b ) ∣ \rho_u(a,b)=|u(a)-u(b)| ρ u ( a , b ) = ∣ u ( a ) − u ( b ) ∣ . For every admissible pseudometric δ \delta δ extending ρ u \rho_u ρ u , the map x ↦ min a ∈ A { u ( a ) + δ ( x , a ) } x\mapsto\min_{a\in A}\{u(a)+\delta(x,a)\} x ↦ min a ∈ A { u ( a ) + δ ( x , a )} is a 1-Lipschitz extension of u u u . Thus, (6) gives E G A ( u ) ≤ E G ( δ 2 ) \mathcal{E}_{G^A}(u)\le\mathcal{E}_G(\delta^2) E G A ( u ) ≤ E G ( δ 2 ) . Taking the infimum over δ \delta δ and using the reverse inequality from Lemma 4.4, we get
E G A ( u ) = min - ext D ( ρ u , G ) ≤ E H ( u ) . \mathcal{E}_{G^A}(u)=\min\text{-}\operatorname{ext}_{\mathcal{D}}(\rho_u,G)\le\mathcal{E}_H(u). E G A ( u ) = min - ext D ( ρ u , G ) ≤ E H ( u ) .
Apply this inequality to each coordinate of f f f and add. We obtain E G A ( f ) ≤ E H ( f ) \mathcal{E}_{G^A}(f)\le\mathcal{E}_H(f) E G A ( f ) ≤ E H ( f ) , and hence
min F ∣ A = f E G ( F ) = E G A ( f ) ≤ E H ( f ) ≤ E H ( d 2 ∣ A × A ) ≤ Q ′ min - ext D ( d ∣ A × A , G ) ≤ Q ′ E G ( d 2 ) . \begin{aligned}
\min_{F|_A=f}\mathcal{E}_G(F)=\mathcal{E}_{G^A}(f)\le\mathcal{E}_H(f) \\
&\le\mathcal{E}_H(d^2|_{A\times A}) \\
&\le Q'\min\text{-}\operatorname{ext}_{\mathcal{D}}(d|_{A\times A},G)\le Q'\mathcal{E}_G(d^2).
\end{aligned} F ∣ A = f min E G ( F ) = E G A ( f ) ≤ E H ( f ) ≤ E H ( d 2 ∣ A × A ) ≤ Q ′ min - ext D ( d ∣ A × A , G ) ≤ Q ′ E G ( d 2 ) .
Here we used that f f f is 1-Lipschitz and the weights of H H H are nonnegative. This energy bound holds for every choice of nonnegative edge weights. Let K = conv ( f ( A ) ) \mathcal{K}=\operatorname{conv}(f(A)) K = conv ( f ( A )) , and let F \mathcal{F} F be the compact convex set of extensions F : X → K F:X\to\mathcal{K} F : X → K of f f f . Projection onto K \mathcal{K} K fixes f ( A ) f(A) f ( A ) and cannot increase distances, so the bound also holds with F ∈ F F\in\mathcal{F} F ∈ F . Let W \mathcal{W} W be the probability simplex on unordered pairs of distinct points of X X X , and put
Φ ( F , w ) = ∑ x < y w x y ( ∥ F ( x ) − F ( y ) ∥ 2 − Q ′ d ( x , y ) 2 ) . \Phi(F,w)=\sum_{x<y}w_{xy}\left(\lVert F(x)-F(y)\rVert^2-Q'd(x,y)^2\right). Φ ( F , w ) = x < y ∑ w x y ( ∥ F ( x ) − F ( y ) ∥ 2 − Q ′ d ( x , y ) 2 ) .
The function Φ \Phi Φ is convex in F F F and affine in w w w . For fixed w w w , the energy bound with g x y = w x y g_{xy}=w_{xy} g x y = w x y gives min F ∈ F Φ ( F , w ) ≤ 0 \min_{F\in\mathcal{F}}\Phi(F,w)\le0 min F ∈ F Φ ( F , w ) ≤ 0 . Sion’s minimax theorem therefore gives
min F ∈ F max w ∈ W Φ ( F , w ) = max w ∈ W min F ∈ F Φ ( F , w ) ≤ 0. \min_{F\in\mathcal{F}}\max_{w\in\mathcal{W}}\Phi(F,w)=\max_{w\in\mathcal{W}}\min_{F\in\mathcal{F}}\Phi(F,w)\le0. F ∈ F min w ∈ W max Φ ( F , w ) = w ∈ W max F ∈ F min Φ ( F , w ) ≤ 0.
The maximum over w w w is the largest of the pairwise differences ∥ F ( x ) − F ( y ) ∥ 2 − Q ′ d ( x , y ) 2 \lVert F(x)-F(y)\rVert^2-Q'd(x,y)^2 ∥ F ( x ) − F ( y ) ∥ 2 − Q ′ d ( x , y ) 2 . Thus, some F ∈ F F\in\mathcal{F} F ∈ F satisfies
∥ F ( x ) − F ( y ) ∥ 2 ≤ Q ′ d ( x , y ) 2 ( x , y ∈ X ) , \lVert F(x)-F(y)\rVert^2\le Q'd(x,y)^2\qquad(x,y\in X), ∥ F ( x ) − F ( y ) ∥ 2 ≤ Q ′ d ( x , y ) 2 ( x , y ∈ X ) ,
including pairs with d ( x , y ) = 0 d(x,y)=0 d ( x , y ) = 0 . Hence, e ( D , A ) ≤ Q ′ e(\mathcal{D},A)\le\sqrt{Q'} e ( D , A ) ≤ Q ′ . Letting Q ′ ↓ Q Q'\downarrow Q Q ′ ↓ Q gives
e ( D , A ) 2 ≤ Q . e(\mathcal{D},A)^2\le Q. e ( D , A ) 2 ≤ Q .
From extensions to sparsifiers. Fix Λ > e ( D , A ) \Lambda>e(\mathcal{D},A) Λ > e ( D , A ) and a weighted graph G G G on X X X . We show that H = G A H=G^A H = G A is a Euclidean metric sparsifier of quality Λ 2 \Lambda^2 Λ 2 .
Let ρ \rho ρ be a Euclidean pseudometric on A A A , and choose a Euclidean realization f : A → ℓ 2 f:A\to\ell_2 f : A → ℓ 2 . For every admissible pseudometric d d d extending ρ \rho ρ , the map f f f is 1-Lipschitz and therefore has an extension to ( X , d ) (X,d) ( X , d ) with Lipschitz constant at most Λ \Lambda Λ . By (6) ,
E G A ( ρ 2 ) = min F ∣ A = f E G ( F ) ≤ Λ 2 E G ( d 2 ) . \mathcal{E}_{G^A}(\rho^2)=\min_{F|_A=f}\mathcal{E}_G(F)\le\Lambda^2\mathcal{E}_G(d^2). E G A ( ρ 2 ) = F ∣ A = f min E G ( F ) ≤ Λ 2 E G ( d 2 ) .
Taking the infimum over d d d gives the left inequality in (4) ; the right inequality follows from Lemma 4.4 . Thus, Q ≤ Λ 2 Q \le\Lambda^2 Q ≤ Λ 2 . Letting Λ ↓ e ( D , A ) \Lambda\downarrow e(D,A) Λ ↓ e ( D , A ) gives
Q ≤ e ( D , A ) 2 . Q \le e(D,A)^2. Q ≤ e ( D , A ) 2 .
Together with the first inequality, this proves the first part of Theorem 4.2 . Taking suprema over partitions gives Q k = e k 2 Q_k=e_k^2 Q k = e k 2 .
Explicit bounds for a union of two sets
For k = 2 k=2 k = 2 we replace the composition of Theorem 2.7 by a direct sum of the two extensions, which avoids the unspecified absolute constants of Lemma 2.1 and gives explicit distortion bounds. The direct sum is the same as in [21 Theorem 2.7] and the improvement over their bound of 8.93 comes from the extension step, where Lemma 4.7 replaces the extension through an α \alpha α -cover of [21 Lemma 2.6] .
Lemma B.1. Let X = A ⊔ B X=A\sqcup B X = A ⊔ B be a finite pseudometric space with A , B ≠ ∅ A,B\ne\varnothing A , B = ∅ , let f A : A → ℓ 2 f_A:A\to\ell_2 f A : A → ℓ 2 and f B : B → ℓ 2 f_B:B\to\ell_2 f B : B → ℓ 2 be isometric embeddings, and let t , s ≥ 1 t,s\ge1 t , s ≥ 1 and 0 < v ≤ ( s + 1 ) 2 0<v\le(s+1)^2 0 < v ≤ ( s + 1 ) 2 . Suppose that f A f_A f A extends to F A : X → ℓ 2 F_A:X\to\ell_2 F A : X → ℓ 2 with ∥ F A ( b ) − F A ( b ′ ) ∥ ≤ t d ( b , b ′ ) \lVert F_A(b)-F_A(b')\rVert\le t d(b,b') ∥ F A ( b ) − F A ( b ′ )∥ ≤ t d ( b , b ′ ) and ∥ F A ( a ) − F A ( b ) ∥ ≤ s d ( a , b ) \lVert F_A(a)-F_A(b)\rVert\le s d(a,b) ∥ F A ( a ) − F A ( b )∥ ≤ s d ( a , b ) for a ∈ A a\in A a ∈ A and b , b ′ ∈ B b,b'\in B b , b ′ ∈ B , and that f B f_B f B extends to F B : X → ℓ 2 F_B:X\to\ell_2 F B : X → ℓ 2 with the same factors on pairs inside A A A and on crossing pairs. Then
c 2 ( X ) 2 ≤ max { 1 + t 2 + v , 2 s 2 + 2 v } ⋅ ( s + 1 ) 2 + v 2 v . c_2(X)^2\le\max\{1+t^2+v,\;2s^2+2v\}\cdot\frac{(s+1)^2+v}{2v}. c 2 ( X ) 2 ≤ max { 1 + t 2 + v , 2 s 2 + 2 v } ⋅ 2 v ( s + 1 ) 2 + v . Proof. Define F = F A ⊕ 2 F B ⊕ 2 v d ( ⋅ , A ) ⊕ 2 v d ( ⋅ , B ) F=F_A\oplus_2F_B\oplus_2\sqrt{v}\,d(\cdot,A)\oplus_2\sqrt{v}\,d(\cdot,B) F = F A ⊕ 2 F B ⊕ 2 v d ( ⋅ , A ) ⊕ 2 v d ( ⋅ , B ) . For a , a ′ ∈ A a,a'\in A a , a ′ ∈ A , the first summand is an isometry, the second has factor t t t , the third vanishes, and ∣ d ( a , B ) − d ( a ′ , B ) ∣ ≤ d ( a , a ′ ) \lvert d(a,B)-d(a',B)\rvert\le d(a,a') ∣ d ( a , B ) − d ( a ′ , B )∣ ≤ d ( a , a ′ ) , so
d ( a , a ′ ) 2 ≤ ∥ F ( a ) − F ( a ′ ) ∥ 2 ≤ ( 1 + t 2 + v ) d ( a , a ′ ) 2 , d(a,a')^2\le\lVert F(a)-F(a')\rVert^2\le(1+t^2+v)d(a,a')^2, d ( a , a ′ ) 2 ≤ ∥ F ( a ) − F ( a ′ ) ∥ 2 ≤ ( 1 + t 2 + v ) d ( a , a ′ ) 2 , while pairs inside B B B are symmetric. For a ∈ A a\in A a ∈ A and b ∈ B b\in B b ∈ B , the first two summands have factor s s s and d ( b , A ) , d ( a , B ) ≤ d ( a , b ) d(b,A),d(a,B)\le d(a,b) d ( b , A ) , d ( a , B ) ≤ d ( a , b ) , so ∥ F ( a ) − F ( b ) ∥ 2 ≤ ( 2 s 2 + 2 v ) d ( a , b ) 2 \lVert F(a)-F(b)\rVert^2\le(2s^2+2v)d(a,b)^2 ∥ F ( a ) − F ( b ) ∥ 2 ≤ ( 2 s 2 + 2 v ) d ( a , b ) 2 .
For the lower bound on crossing pairs, put d = d ( a , b ) d=d(a,b) d = d ( a , b ) , r A = d ( b , A ) r_A=d(b,A) r A = d ( b , A ) , r B = d ( a , B ) r_B=d(a,B) r B = d ( a , B ) and C = s + 1 C=s+1 C = s + 1 , and let a ′ ∈ A a'\in A a ′ ∈ A satisfy d ( a ′ , b ) = r A d(a',b)=r_A d ( a ′ , b ) = r A . Then
d ≤ d ( a , a ′ ) + r A = ∥ F A ( a ) − F A ( a ′ ) ∥ + r A ≤ ∥ F A ( a ) − F A ( b ) ∥ + s r A + r A , d\le d(a,a')+r_A=\lVert F_A(a)-F_A(a')\rVert+r_A\le\lVert F_A(a)-F_A(b)\rVert+sr_A+r_A, d ≤ d ( a , a ′ ) + r A = ∥ F A ( a ) − F A ( a ′ )∥ + r A ≤ ∥ F A ( a ) − F A ( b )∥ + s r A + r A , so ∥ F A ( a ) − F A ( b ) ∥ ≥ ( d − C r A ) + \lVert F_A(a)-F_A(b)\rVert\ge(d-Cr_A)_+ ∥ F A ( a ) − F A ( b )∥ ≥ ( d − C r A ) + , and symmetrically ∥ F B ( a ) − F B ( b ) ∥ ≥ ( d − C r B ) + \lVert F_B(a)-F_B(b)\rVert\ge(d-Cr_B)_+ ∥ F B ( a ) − F B ( b )∥ ≥ ( d − C r B ) + . Over r ≥ 0 r\ge0 r ≥ 0 , the function ( d − C r ) + 2 + v r 2 (d-Cr)_+^2+vr^2 ( d − C r ) + 2 + v r 2 is minimized at r = C d / ( C 2 + v ) r=Cd/(C^2+v) r = C d / ( C 2 + v ) with value v d 2 / ( C 2 + v ) vd^2/(C^2+v) v d 2 / ( C 2 + v ) , hence
∥ F ( a ) − F ( b ) ∥ 2 ≥ ( d − C r A ) + 2 + v r A 2 + ( d − C r B ) + 2 + v r B 2 ≥ 2 v C 2 + v d 2 . \lVert F(a)-F(b)\rVert^2\ge(d-Cr_A)_+^2+vr_A^2+(d-Cr_B)_+^2+vr_B^2\ge\frac{2v}{C^2+v}d^2. ∥ F ( a ) − F ( b ) ∥ 2 ≥ ( d − C r A ) + 2 + v r A 2 + ( d − C r B ) + 2 + v r B 2 ≥ C 2 + v 2 v d 2 . Since v ≤ C 2 v\le C^2 v ≤ C 2 , this factor is at most one, so F F F contracts no distance by more than 2 v / ( C 2 + v ) \sqrt{2v/(C^2+v)} 2 v / ( C 2 + v ) , and the ratio of the squared expansion and contraction bounds proves the claim.
Theorem B.2. Let ( X , d ) (X,d) ( X , d ) be a metric space with X = U 1 ∪ U 2 X=U_1\cup U_2 X = U 1 ∪ U 2 . If U 1 U_1 U 1 and U 2 U_2 U 2 are Euclidean in X X X , then c 2 ( X ) < 3.997 c_2(X)<3.997 c 2 ( X ) < 3.997 . If U 1 U_1 U 1 and U 2 U_2 U 2 are 1-strongly Euclidean in X X X , then c 2 ( X ) ≤ 3.86 c_2(X)\le3.86 c 2 ( X ) ≤ 3.86 .
Together with the examples of [21 Theorem 1.2] , which require distortion arbitrarily close to 3, the first bound gives 3 ≤ D 2 < 3.997 3\le D_2<3.997 3 ≤ D 2 < 3.997 .
Proof. We may pass to a partition X = A ⊔ B X=A\sqcup B X = A ⊔ B , since both conditions are preserved by restriction. If one part is empty, X X X is Euclidean. Fix isometric embeddings f A f_A f A and f B f_B f B .
Suppose first that A A A and B B B are Euclidean. For a graph on X X X with no edges inside A A A , (9) reads E G A ( d 2 ) ≤ I + 2 J + 2 π I J \mathcal{E}_{G_A}(d^2)\le I+2J+\sqrt{2\pi IJ} E G A ( d 2 ) ≤ I + 2 J + 2 π I J , so by (15) , Lemma 4.7 extends f A f_A f A with factor t t t on B B B and s s s on crossing pairs whenever s ≥ 2 s\ge\sqrt{2} s ≥ 2 and ( t 2 − 1 ) ( s 2 − 2 ) ≥ π / 2 (t^2-1)(s^2-2)\ge\pi/2 ( t 2 − 1 ) ( s 2 − 2 ) ≥ π /2 , and exchanging the roles of A A A and B B B does the same for f B f_B f B . Take s = 187 / 125 s=187/125 s = 187/125 , v = 413 / 100 v=413/100 v = 413/100 and t 2 = 2 s 2 + v − 1 = 7.606032 t^2=2s^2+v-1=7.606032 t 2 = 2 s 2 + v − 1 = 7.606032 . Then ( t 2 − 1 ) ( s 2 − 2 ) = 1.572341312512 > 11 / 7 > π / 2 (t^2-1)(s^2-2)=1.572341312512>11/7>\pi/2 ( t 2 − 1 ) ( s 2 − 2 ) = 1.572341312512 > 11/7 > π /2 and v ≤ ( s + 1 ) 2 v\le(s+1)^2 v ≤ ( s + 1 ) 2 , and the two terms in the maximum of Lemma B.1 coincide, so
c 2 ( X ) 2 ≤ ( s 2 + v ) ( ( s + 1 ) 2 + v ) v < 15.975 < 3.997 2 . c_2(X)^2\le\frac{(s^2+v)((s+1)^2+v)}{v}<15.975<3.997^2. c 2 ( X ) 2 ≤ v ( s 2 + v ) (( s + 1 ) 2 + v ) < 15.975 < 3.99 7 2 .
Suppose now that A A A and B B B are 1-strongly Euclidean, so B ∪ { a } B\cup\{a\} B ∪ { a } is Euclidean for every a ∈ A a\in A a ∈ A , and symmetrically. For a graph with no edges inside A A A , Lemma 4.9 gives
E G A ( d 2 ) ≤ 1 2 ( J + J + I + I ( I + J ) ) 2 , \mathcal{E}_{G_A}(d^2)\le\frac{1}{2}\left(\sqrt{J}+\sqrt{J+I+\sqrt{I(I+J)}}\right)^2, E G A ( d 2 ) ≤ 2 1 ( J + J + I + I ( I + J ) ) 2 ,
and the substitution z = 1 − I / ( I + J ) z=\sqrt{1-\sqrt{I/(I+J)}} z = 1 − I / ( I + J ) preceding (19) shows that this is at most t 2 I + s 2 J t^2I+s^2J t 2 I + s 2 J for all I , J ≥ 0 I,J\ge0 I , J ≥ 0 precisely when
D ( z ) = t 2 ( 1 − z 2 ) 2 + s 2 z 2 ( 2 − z 2 ) − 1 2 ( z + 1 ) 2 ( 2 − z 2 ) ≥ 0 ( 0 ≤ z ≤ 1 ) . D(z)=t^2(1-z^2)^2+s^2z^2(2-z^2)-\frac{1}{2}(z+1)^2(2-z^2)\ge0\qquad(0\le z\le1). D ( z ) = t 2 ( 1 − z 2 ) 2 + s 2 z 2 ( 2 − z 2 ) − 2 1 ( z + 1 ) 2 ( 2 − z 2 ) ≥ 0 ( 0 ≤ z ≤ 1 ) .
Take s = 143 / 100 s=143/100 s = 143/100 , t 2 = 3 s 2 + s − 1 = 6.5647 t^2=3s^2+s-1=6.5647 t 2 = 3 s 2 + s − 1 = 6.5647 and v = s ( s + 1 ) v=s(s+1) v = s ( s + 1 ) . Then
D ( z ) = 5.5647 − 2 z − 9.5396 z 2 + z 3 + 5.0198 z 4 D(z)=5.5647-2z-9.5396z^2+z^3+5.0198z^4 D ( z ) = 5.5647 − 2 z − 9.5396 z 2 + z 3 + 5.0198 z 4
D ′ ( z ) = − 2 − 19.0792 z + 3 z 2 + 20.0792 z 3 . D'(z)=-2-19.0792z+3z^2+20.0792z^3. D ′ ( z ) = − 2 − 19.0792 z + 3 z 2 + 20.0792 z 3 .
On [ 0 , 9 / 10 ] [0,9/10] [ 0 , 9/10 ] , the bounds z 2 ≤ 9 z / 10 z^2\le9z/10 z 2 ≤ 9 z /10 and z 3 ≤ 81 z / 100 z^3\le81z/100 z 3 ≤ 81 z /100 give D ′ ( z ) ≤ − 2 − 0.115048 z < 0 D'(z)\le-2-0.115048z<0 D ′ ( z ) ≤ − 2 − 0.115048 z < 0 , so D D D is decreasing there. On [ 9 / 10 , 1 ] [9/10,1] [ 9/10 , 1 ] , the increasing function D ′ ′ ( z ) = − 19.0792 + 6 z + 60.2376 z 2 D''(z)=-19.0792+6z+60.2376z^2 D ′′ ( z ) = − 19.0792 + 6 z + 60.2376 z 2 is at least D ′ ′ ( 9 / 10 ) > 35 D''(9/10)>35 D ′′ ( 9/10 ) > 35 , while exact substitution at z 0 = 191 / 200 z_0=191/200 z 0 = 191/200 gives D ( z 0 ) > 7.476 ⋅ 10 − 4 D(z_0)>7.476\cdot10^{-4} D ( z 0 ) > 7.476 ⋅ 1 0 − 4 and 0 < D ′ ( z 0 ) < 4.1 ⋅ 10 − 3 0<D'(z_0)<4.1\cdot10^{-3} 0 < D ′ ( z 0 ) < 4.1 ⋅ 1 0 − 3 , so
D ( z ) ≥ D ( z 0 ) + D ′ ( z 0 ) ( z − z 0 ) + 35 2 ( z − z 0 ) 2 ≥ D ( z 0 ) − D ′ ( z 0 ) 2 70 > 0. D(z)\ge D(z_0)+D'(z_0)(z-z_0)+\frac{35}{2}(z-z_0)^2\ge D(z_0)-\frac{D'(z_0)^2}{70}>0. D ( z ) ≥ D ( z 0 ) + D ′ ( z 0 ) ( z − z 0 ) + 2 35 ( z − z 0 ) 2 ≥ D ( z 0 ) − 70 D ′ ( z 0 ) 2 > 0.
Hence, D ≥ 0 D\ge0 D ≥ 0 on [ 0 , 1 ] [0,1] [ 0 , 1 ] , and Lemma 4.7 provides both extensions with these factors. Since v ≤ ( s + 1 ) 2 v\le(s+1)^2 v ≤ ( s + 1 ) 2 and 1 + t 2 + v = 2 s 2 + 2 v = 4 s 2 + 2 s 1+t^2+v=2s^2+2v=4s^2+2s 1 + t 2 + v = 2 s 2 + 2 v = 4 s 2 + 2 s , Lemma B.1 gives
c 2 ( X ) 2 ≤ ( 4 s 2 + 2 s ) ⋅ ( s + 1 ) ( 2 s + 1 ) 2 s ( s + 1 ) = ( 2 s + 1 ) 2 , c_2(X)^2\le(4s^2+2s)\cdot\frac{(s+1)(2s+1)}{2s(s+1)}=(2s+1)^2, c 2 ( X ) 2 ≤ ( 4 s 2 + 2 s ) ⋅ 2 s ( s + 1 ) ( s + 1 ) ( 2 s + 1 ) = ( 2 s + 1 ) 2 ,
that is, c 2 ( X ) ≤ 2 s + 1 = 3.86 c_2(X)\le2s+1=3.86 c 2 ( X ) ≤ 2 s + 1 = 3.86 .
Reduction to finite metric spaces
For completeness, we include the reduction from separable to finite metric spaces used in § 2 . The argument follows [21 Lemma 2.9, part I] .
Lemma C.1 (Finite subsets determine distortion). Let ( X , d ) (X,d) ( X , d ) be a separable metric space, and let D ≥ 1 D\ge1 D ≥ 1 . If every finite subset of X X X embeds into ℓ 2 \ell_2 ℓ 2 with distortion at most D D D , then so does $X.
Proof. Choose a countable dense set { x 0 , x 1 , … } \{x_0,x_1,\ldots\} { x 0 , x 1 , … } in X X X . For every N N N , the finite set { x 0 , … , x N } \{x_0,\ldots,x_N\} { x 0 , … , x N } has an embedding F N F_N F N with distortion at most D D D , normalized by F N ( x 0 ) = 0 F_N(x_0)=0 F N ( x 0 ) = 0 and d ( x , y ) ≤ ∥ F N ( x ) − F N ( y ) ∥ ≤ D d ( x , y ) d(x,y) \le\lVert F_N(x)-F_N(y)\rVert\le Dd(x,y) d ( x , y ) ≤ ∥ F N ( x ) − F N ( y )∥ ≤ D d ( x , y ) .
The inner products ⟨ F N ( x m ) , F N ( x n ) ⟩ \langle F_N(x_m),F_N(x_n)\rangle ⟨ F N ( x m ) , F N ( x n )⟩ are bounded for every fixed pair, so a diagonal subsequence makes all of them converge. Every finite matrix of limiting inner products is positive semidefinite, so the limits are realized by vectors F ( x n ) F(x_n) F ( x n ) in a separable Hilbert space, and the distance bounds persist. Completeness extends F F F from the dense set to all of X X X with the same bounds.
Let X = U 1 ∪ ⋯ ∪ U k X=U_1\cup\cdots\cup U_k X = U 1 ∪ ⋯ ∪ U k be as in any of Theorems 1.1 to 1.4 and 1.7, and let S ⊆ X S\subseteq X S ⊆ X be finite. The nonempty sets S ∩ U i S\cap U_i S ∩ U i cover S S S , and they inherit each property we impose on the U i U_i U i : being Euclidean or h h h -strongly Euclidean in X X X , and embedding into ℓ 2 \ell_2 ℓ 2 with a given distortion, are preserved by restriction. So is the bound e ( U i ) ≤ Λ e(U_i)\le\Lambda e ( U i ) ≤ Λ . By Kirszbraun’s Theorem [14 ] , a Lipschitz map on S ∩ U i S\cap U_i S ∩ U i extends to the Euclidean set U i U_i U i with the same Lipschitz constant, then to X X X with constant multiplied by Λ \Lambda Λ , and restricting the result to S S S does not increase it. Moreover, X X X is separable, since each U i U_i U i embeds into ℓ 2 \ell_2 ℓ 2 with finite distortion. Hence, by Lemma C.1 , each upper bound on c 2 ( X ) c_2(X) c 2 ( X ) in these theorems follows from its finite case.