Introduction
Approximating a graph by a distribution of trees is a basic tool in the design of algorithms. A tree makes distances, cuts, and routing particularly simple, while a distribution can avoid the distortion forced on any one tree. There are, however, two different notions of a tree approximation. An abstract tree metric may introduce Steiner vertices and edges unrelated to the input graph. A spanning tree must use the original vertices and original edges, with their given lengths. The latter restriction is essential when a tree path must be a path in the network itself.
Alon, Karp, Peleg, and West [5 , 6 ] formulated the spanning-tree problem as a zero-sum game between a tree player and an edge player. Their extended abstract appeared at the 1991 DIMACS Workshop on On-Line Algorithms, with proceedings published in 1992; the full paper appeared in 1995. They proved a bound of exp ( O ( log n log log n ) ) \exp(O(\sqrt{\log n\log\log n})) exp ( O ( log n log log n )) and conjectured an O ( log n ) O(\log n) O ( log n ) bound. Their original payoff differs pointwise by at most one from d T , ℓ ( u , v ) / ℓ u v d_{T,\ell}(u,v)/\ell_{uv} d T , ℓ ( u , v ) / ℓ uv , so the two game values differ by at most one. The modified game asks for a distribution over spanning trees whose expected stretch is small on every original edge. An edge guarantee immediately extends to every pair of vertices by summing along a shortest graph path.
Subsequent work brought the upper bound close to logarithmic. Elkin, Emek, Spielman, and Teng [13 ] obtained average stretch O ( log 2 n log log n ) O(\log^{2} n \log\log n) O ( log 2 n log log n ) . Abraham, Bartal, and Neiman [2 ] obtained stochastic stretch O ( log n log log n ( log log log n ) 3 ) O(\log n\log\log n(\log\log\log n)^{3}) O ( log n log log n ( log log log n ) 3 ) . The petal decomposition of Abraham and Neiman [7 ] gives average stretch O ( log n log log n ) O(\log n\log\log n) O ( log n log log n ) and a construction with the same logarithmic factors in its nearly linear running time. Low-stretch trees also underlie preconditioners for fast solvers of symmetric diagonally dominant linear systems [23 ] . The remaining gap is a geometric one: local partitions must be chosen so that they can be reconnected using original edges without losing the distance estimate.
We close this gap by proving the logarithmic bound for spanning trees. The theorem is weighted: the demand assigned to an edge is independent of its length. This is the form needed both for the spanning-tree game and for a capacity decomposition. Our construction runs in polynomial time for arbitrary rational input. Its running time is not nearly linear; the contribution is the logarithmic approximation within the spanning-tree family and an explicit polynomial construction.
Results
Let G = ( V , E ) G=(V,E) G = ( V , E ) be a connected undirected graph, with n = ∣ V ∣ ≥ 2 n=|V|\ge2 n = ∣ V ∣ ≥ 2 and m = ∣ E ∣ m=|E| m = ∣ E ∣ . Throughout the paper, a spanning tree has vertex set V V V and uses only edges of G G G . For lengths ℓ e > 0 \ell_e>0 ℓ e > 0 , write d G , ℓ d_{G,\ell} d G , ℓ and d T , ℓ d_{T,\ell} d T , ℓ for graph and tree distances. For capacities c e > 0 c_e>0 c e > 0 and a tree edge f f f , let S T , f S_{T,f} S T , f be either component of T − f T-f T − f . Define
c T ( f ) = c G ( δ G ( S T , f ) ) , L T ( f ) = { c T ( f ) , f ∈ T , 0 , f ∉ T . (1) c_T(f)=c_G(\delta_G(S_{T,f})),\qquad L_T(f)=
\begin{cases}
c_T(f), & f\in T,\\
0, & f\notin T.
\end{cases} \tag*{(1)} c T ( f ) = c G ( δ G ( S T , f )) , L T ( f ) = { c T ( f ) , 0 , f ∈ T , f ∈ / T . (1)
Thus the tree topology is a subgraph, while its capacities are determined by its fundamental cuts. These are exactly the smallest capacities on that topology that dominate every graph cut (Proposition 9.5 ).
Theorem 1.1 (Logarithmic spanning-tree decompositions). There is an absolute integer C 0 > 0 C_0>0 C 0 > 0 with the following property. Put ρ n = C 0 ⌈ log 2 ( 2 n ) ⌉ \rho_n=C_0\lceil\log_2(2n)\rceil ρ n = C 0 ⌈ log 2 ( 2 n )⌉ and fix 0 < ϑ ≤ 1 0<\vartheta\le1 0 < ϑ ≤ 1 . Let α n = ( 1 + ϑ ) ρ n \alpha_n=(1+\vartheta)\rho_n α n = ( 1 + ϑ ) ρ n .
(i) For every positive length vector ℓ \ell ℓ , there is a distribution q d i s t q^{\mathrm{dist}} q dist over spanning trees such that, for all u , v ∈ V u,v\in V u , v ∈ V ,
d T , ℓ ( u , v ) ≥ d G , ℓ ( u , v ) for every supported T , d_{T,\ell}(u,v)\ge d_{G,\ell}(u,v)\qquad\text{for every supported }T, d T , ℓ ( u , v ) ≥ d G , ℓ ( u , v ) for every supported T , E q d i s t d T , ℓ ( u , v ) ≤ α n d G , ℓ ( u , v ) ( u , v ∈ V ) . \mathbb{E}_{q^{\mathrm{dist}}}d_{T,\ell}(u,v)\le\alpha_n d_{G,\ell}(u,v)\qquad(u,v\in V). E q dist d T , ℓ ( u , v ) ≤ α n d G , ℓ ( u , v ) ( u , v ∈ V ) . (ii) For every positive capacity vector c c c , there is a distribution q c a p q^{\mathrm{cap}} q cap over spanning trees, equipped with (1), such that
E q c a p L T ( f ) ≤ α n c f ( f ∈ E ) . \mathbb{E}_{q^{\mathrm{cap}}}L_T(f)\le\alpha_n c_f\qquad(f\in E). E q cap L T ( f ) ≤ α n c f ( f ∈ E ) . Consequently, for every A ⊆ V A\subseteq V A ⊆ V ,
c T ( δ T ( A ) ) ≥ c G ( δ G ( A ) ) for every supported T , c_T(\delta_T(A))\ge c_G(\delta_G(A))\qquad\text{for every supported }T, c T ( δ T ( A )) ≥ c G ( δ G ( A )) for every supported T , E q c a p c T ( δ T ( A ) ) ≤ α n c G ( δ G ( A ) ) ( A ⊆ V ) . \mathbb{E}_{q^{\mathrm{cap}}}c_T(\delta_T(A))\le\alpha_n c_G(\delta_G(A))\qquad(A\subseteq V). E q cap c T ( δ T ( A )) ≤ α n c G ( δ G ( A )) ( A ⊆ V ) . Splitting every pair’s demand among its tree paths with the fixed probabilities q T c a p q_T^{\mathrm{cap}} q T cap gives an α n \alpha_n α n -competitive fractional oblivious routing.
Each distribution can be supported on at most m m m trees. For rational input and rational ϑ \vartheta ϑ , it can be constructed deterministically using at most
1 + ⌈ 20 m ⌈ log 2 ( 2 m ) ⌉ ϑ 2 ⌉ 1+\left\lceil\frac{20m\lceil\log_2(2m)\rceil}{\vartheta^2}\right\rceil 1 + ⌈ ϑ 2 20 m ⌈ log 2 ( 2 m )⌉ ⌉
weighted-tree calls and polynomial additional work. The complete running time is polynomial in the binary input length and 1 / ϑ 1/\vartheta 1/ ϑ .
The capacity distribution is fixed before the cut or the demand vector is chosen. Its guarantee is max f E [ L T ( f ) / c f ] \max_f \mathbb{E}[L_T(f)/c_f] max f E [ L T ( f ) / c f ] , and the routing uses the fractional mixture of tree paths. The theorem does not require each sampled tree to have small congestion. The distance and capacity distributions may differ; this is necessary even when all graph lengths and capacities are one (Proposition 10.4 ).
The main ingredient is a deterministic weighted-tree theorem. Given arbitrary demands a e ≥ 0 a_e \ge0 a e ≥ 0 and lengths ℓ e > 0 \ell_e > 0 ℓ e > 0 , it constructs a spanning tree satisfying
∑ e = u v ∈ E a e d T , ℓ ( u , v ) ≤ ρ n ∑ e ∈ E a e ℓ e . (2) \sum_{e=uv\in E} a_e d_{T,\ell}(u,v) \le\rho_n \sum_{e\in E} a_e\ell_e. \tag*{(2)} e = uv ∈ E ∑ a e d T , ℓ ( u , v ) ≤ ρ n e ∈ E ∑ a e ℓ e . (2)
The geometric and entropy arguments proving this statement occupy § 4 –§ 6 . We use symbolic constants there and give an explicit admissible choice, including C 0 C_0 C 0 , in § 11 . The distribution construction in § 8 works with any weighted response bound ρ \rho ρ , loses only a factor 1 + ϑ 1+\vartheta 1 + ϑ , and has a call bound independent of the numerical size of ρ \rho ρ .
Tightness on diamonds. The logarithmic lower bound for probabilistic tree metrics on recursive diamonds is classical [18 ] . We give a direct spanning-tree proof in § 10.1 . For every spanning tree of the height-k k k unit diamond D k D_k D k , the sum of tree distances over the 4 k 4^k 4 k original edges is at least k 4 k / 9 k4^k/9 k 4 k /9 . The same inequality, by double counting, gives an Ω ( log n ) \Omega(\log n) Ω ( log n ) lower bound for maximum expected fundamental-cut congestion. Summing singleton tree cuts gives the lower bound for expected cut distortion as well. Thus all three logarithmic bounds in Theorem 1.1 are asymptotically tight on the same series-parallel family.
The size of an explicit decomposition. Charikar, Chekuri, Goel, Guha, and Plotkin [10 ] showed that O ( n log n ) O(n\log n) O ( n log n ) abstract tree metrics suffice for their probabilistic metric embedding. This raises the corresponding support question for spanning trees. Our support bound m m m cannot in general be replaced by O ( n log n ) O(n\log n) O ( n log n ) for fundamental-cut congestion, even when the optimal congestion approaches one. For an integer k ≥ 2 k \ge2 k ≥ 2 , let G k G_k G k consist of two k k k -vertex stars whose edges have capacity k 2 k^2 k 2 , together with all k 2 k^2 k 2 cross edges of capacity one.
Theorem 1.2 (Support required by fundamental-cut congestion). For G k G_k G k , every distribution with maximum expected congestion at most ρ \rho ρ has support at least ⌈ k 2 / ρ ⌉ \lceil k^2/\rho\rceil ⌈ k 2 / ρ ⌉ . There is a distribution with congestion at most β k = 1 + 2 ( k − 1 ) / k 2 \beta_k=1+2(k-1)/k^2 β k = 1 + 2 ( k − 1 ) / k 2 . Moreover, for every
ρ ≥ γ k = k 2 + 2 k − 2 k 2 − 1 , \rho\ge\gamma_k=\frac{k^2+2k-2}{k^2-1}, ρ ≥ γ k = k 2 − 1 k 2 + 2 k − 2 , the minimum support size is exactly ⌈ k 2 / ρ ⌉ \lceil k^2/\rho\rceil ⌈ k 2 / ρ ⌉ . On the same family, one fundamental-cut tree approximates every cut within factor three; with unit metric lengths, one spanning tree has distance distortion at most three.
The proof is in § 10.2 . Since n = 2 k n = 2k n = 2 k and m = k 2 + 2 k − 2 m = k^{2} + 2k - 2 m = k 2 + 2 k − 2 , the theorem gives Ω ( m / ρ ) \Omega(m/\rho) Ω ( m / ρ ) support. In particular, constant congestion can require Θ ( m ) \Theta(m) Θ ( m ) trees, and logarithmic congestion can require Ω ( n 2 / log n ) \Omega(n^{2}/\log n) Ω ( n 2 / log n ) trees. The lower bound concerns individual original-edge loads, the property used by the oblivious routing guarantee. It does not follow from cut approximation alone. At the attainable bound ρ = β k = 1 + O ( 1 / n ) \rho= \beta_k = 1 + O(1/n) ρ = β k = 1 + O ( 1/ n ) , every distribution needs at least m − O ( n ) = ( 1 − o ( 1 ) ) m m - O(n) = (1 - o(1))m m − O ( n ) = ( 1 − o ( 1 )) m trees (Corollary 10.3 ). Thus the dependence on m m m is asymptotically sharp for congestion tending to one.
Related work
Probabilistic metric embeddings. Bartal [8 Bar98] established probabilistic approximation by hierarchically well-separated trees (HSTs), first with distortion O ( log 2 n ) O(\log^{2} n) O ( log 2 n ) and then with distortion O ( log n log log n ) O(\log n \log\log n) O ( log n log log n ) . Fakcharoenphol, Rao, and Talwar [16 ] obtained the optimal O ( log n ) O(\log n) O ( log n ) bound: every finite metric admits a distribution of dominating tree metrics with this expected distortion for every pair. Their construction also produces HSTs, whose edge lengths decrease geometrically along root-to-leaf paths. The distinction between deterministic domination and expected distortion is central to the applications. A fixed feasible solution has a cost controlled in expectation by linearity, while domination supports transferring a tree solution back to the original metric. The allowed trees can have Steiner vertices and edges outside the input graph. Replacing such edges by graph paths can create cycles; extracting a spanning tree from their union need not preserve the distance bound. Consequently the metric theorem does not supply the original-edge guarantee studied here.
Approximation and online algorithms. Metric labeling illustrates the role of HST structure. Kleinberg and Tardos [20 ] give a constant-factor algorithm when the metric on the k k k labels is an HST with a fixed separation greater than two. Combined with FRT, this yields an O ( log k ) O(\log k) O ( log k ) approximation for arbitrary label metrics: assignment costs are unchanged, and the metric controls the pairwise separation costs. The related 0-extension problem fixes the labels of designated terminals. Calinescu, Karloff, and Rabani [11 ] obtain an O ( log k ) O(\log k) O ( log k ) approximation, where k k k counts terminals, through a metric relaxation and randomized partitions. Their partitioning method is also an ingredient in the development of FRT [16 ] ; this is a connection between the techniques, rather than an attribution of the 0-extension theorem to FRT.
For network design, Awerbuch and Azar [1 ] use tree embeddings for uniform buy-at-bulk, where cable types and their costs per unit length are common to all edges. Their constant-factor tree algorithm and the FRT distortion bound give an O ( log n ) O(\log n) O ( log n ) expected approximation, including the multiple source–sink setting. For group Steiner tree with k k k groups, Garg, Konjevod, and Ravi [17 ] obtain an O ( log n log k ) O(\log n\log k) O ( log n log k ) approximation on trees; combining it with FRT gives an O ( log 2 n log k ) O(\log^{2} n\log k) O ( log 2 n log k ) on general graphs. In online computation, Fiat and Mendel [15 ] give an O ( log n log log n ) O(\log n\log\log n) O ( log n log log n ) -competitive algorithm on HSTs. Its classical combination with FRT gives an O ( log 2 n log log n ) O(\log^{2} n\log\log n) O ( log 2 n log log n ) randomized competitive ratio for metrical task systems against oblivious adversaries [16 ] . These examples explain why both the distortion and the structure of the target trees matter: the reduction must preserve the particular objective and feasibility conditions of its application.
Spanning trees and numerical algorithms. The original-edge problem has a separate history, beginning with the tree edge game of Alon, Karp, Peleg, and West [6 ] , motivated by the k k k -server problem on road networks. Their reduction relates the game value to a randomized competitive ratio against an oblivious adversary. The subsequent star decompositions of Elkin, Emek, Spielman, and Teng [13 ] , the portal highways of Abraham, Bartal, and Neiman [2 ] , and the petal decompositions of Abraham and Neiman [7 ] successively improved the distortion while maintaining a spanning subgraph. The latter construction achieves O ( log n log log n ) O(\log n \log\log n) O ( log n log log n ) average stretch in O ( m log n log log n ) O(m\log n\log\log n) O ( m log n log log n ) time. For unweighted series-parallel graphs, Emek and Peleg [14 ] already obtain a tight logarithmic distribution over spanning trees. The diamond lower bound of Gupta, Newman, Rabinovich, and Sinclair [18 ] therefore holds even in a graph class where the spanning-tree question was previously settled.
Low-stretch spanning trees also enter combinatorial preconditioning. Spielman and Teng [23 ] construct sparse subgraph preconditioners by adding edges to low-stretch trees and use recursive preconditioning to solve symmetric diagonally dominant systems in nearly linear time. This application requires efficient construction, in addition to a stretch bound. Our polynomial-time weighted oracle establishes the logarithmic distortion guarantee; obtaining a comparable running time remains a separate algorithmic question.
Recent work also studies faster constructions. Kyng et al. [19 ] obtain efficient logarithmic metric-tree embeddings and ℓ 1 \ell_1 ℓ 1 -oblivious routing by analyzing random-shift decompositions across scales. Dorkofikis et al. [12 ] obtain spanning trees from randomly perturbed shortest-path computations. Their worst-case construction has O ( log 5 n ) O(\log^5 n) O ( log 5 n ) expected stretch and O ~ ( m ) \widetilde{O}(m) O ( m ) running time for positive, polynomially bounded integer lengths. These results address computational efficiency and alternative constructions alongside the distortion question.
Cut decompositions and minimum bisection. Räcke [22 ] constructs a convex combination of capacitated decomposition trees that supports O ( log n ) O(\log n) O ( log n ) -competitive fractional oblivious routing in undirected networks. The leaves correspond to graph vertices, while tree edges are mapped to graph paths; these trees are not required to be spanning subgraphs. The same framework gives an O ( log n ) O(\log n) O ( log n ) approximation for minimum bisection with exact cardinality balance: when n n n is even, both output parts have exactly n / 2 n/2 n /2 vertices, and the cut capacity is within the stated factor of the minimum among such bisections. Räcke computes a minimum-cost leaf bisection in each supported tree by dynamic programming and returns the best induced graph bisection [22 Section 3, Lemma 9] . The balance counts leaves, and hence original vertices; internal tree nodes do not contribute to it. The framework also yields logarithmic approximations for sparsest cut and multicut. Thus the logarithmic bisection approximation is an established application. Corollary 9.4 recovers this guarantee, and its analogue for any prescribed side cardinality, using spanning trees on the original vertex set with prescribed fundamental-cut capacities. The additional contribution concerns this representation and its construction.
Family-preserving duality and explicit distributions. Andersen and Feige [3 Theorem 6] formulate the equivalence between stretch and congestion for a fixed family of admissible mappings. Their discussion explicitly distinguishes the distributions needed for the two guarantees. For spanning trees, the mapping sends each edge to its tree path, and the induced loads are precisely fundamental-cut capacities. They already use this family for a near-logarithmic minimum-bisection approximation and explain how abstract tree metrics recover Räcke’s logarithmic bound [3 Sections 4.5 and 5] . Our weighted theorem supplies a logarithmic response within the spanning-tree family itself.
Charikar, Chekuri, Goel, Guha, and Plotkin [10 ] give a deterministic polynomial construction of a distribution on O ( n log n ) O(n\log n) O ( n log n ) metric trees with O ( log n log log n ) O(\log n\log\log n) O ( log n log log n ) distortion. Their work predates FRT and connects weighted response oracles to explicit distributions through fractional packing. Our rational master follows this framework and the multiplicative-weights method [21 , 4 ] ; it gives a factor 1 + ϑ 1+\vartheta 1 + ϑ loss, uniformly bounded query encodings, and exact compression to m m m trees. The support obstruction proved here concerns fundamental-cut congestion on original-edge trees. Its Ω ( n 2 / ρ ) \Omega(n^2/\rho) Ω ( n 2 / ρ ) bound explains why the O ( n log n ) O(n\log n) O ( n log n ) support guarantee for metric approximation cannot be transferred to that requirement, even at logarithmic congestion.
Organization. Section 2 explains the proof and construction. Section 3 fixes conventions. Sections 4–7 give the weighted theorem and its running time. Sections 8 and 9 construct the distributions and prove the capacity and routing statements. Section 10 proves the lower bounds, and Section 11 checks the constants.
Our approach and techniques
Our starting point is a weighted formulation of the spanning-tree problem. Give each graph edge a positive length ℓ e \ell_e ℓ e and an arbitrary nonnegative demand a e a_e a e . We construct one actual spanning tree T T T for which
∑ e = u v a e d T , ℓ ( u , v ) ≤ O ( log n ) ∑ e a e ℓ e . \sum_{e=uv} a_e d_{T,\ell}(u,v) \le O(\log n)\sum_e a_e\ell_e. e = uv ∑ a e d T , ℓ ( u , v ) ≤ O ( log n ) e ∑ a e ℓ e .
The demands may depend on the lengths. This freedom is what allows the same theorem to yield both distance and capacity decompositions. For distances, a weighted response controls a specified linear combination of edge stretches. For capacities, given positive edge prices p e p_e p e , we use auxiliary lengths p e / c e p_e/c_e p e / c e and demands c e c_e c e . Double counting the resulting tree paths identifies the weighted distance objective with the priced sum of fundamental-cut congestions. This is the family-preserving distance–capacity equivalence of Räcke and Andersen–Feige [22 , 3 ] . Our task is therefore to prove the weighted bound while keeping every tree edge in the input graph. The distance and capacity distributions obtained from it may be different.
A demand routable in the graph with congestion γ \gamma γ induces load at most γ c T ( f ) \gamma c_T(f) γ c T ( f ) on every tree edge, by the cut inequality for its fundamental cut. Averaging the corresponding original-edge tree paths therefore gives the stated oblivious routing.
The logarithmic weighted bound
Protecting an entire shortest path. Recursive tree constructions must control the distances created when their clusters are reconnected. Petal decompositions achieve this by coordinating cluster boundaries with paths leading toward a root [2 , 7 ] . We retain this geometric principle and change the object carried by the recursion: each instance protects an entire physical shortest path starting at its root. Every child avoids that path, which is retained verbatim in the final tree. The child receives a new protected path, inherited from the root’s shortest-path tree through the chosen source; see Figure 1 .
Figure 1. A recursive step. The child avoids the inherited protected path and is joined through the original edge a p b p a_p b_p a p b p . Its new protected path runs from a p a_p a p to the selected source p p p . Only the protected paths, the relevant root-path segment, and one branch of the child are shown; physical lengths are not to scale.
The construction uses temporary lengths that halve precisely the protected-path edges. A petal boundary is defined by a distance gauge adapted to a shortest-path tree in these temporary lengths. Removing a petal leaves an ancestor-closed residual, so its surviving root distances stay unchanged. Each child has one designated original portal edge leading to the later residual. Its witness paths lie entirely within the child and contain no edge of the old protected path. Consequently they have their original physical lengths, and the new designated path is physically shortest in the induced child. These properties are preserved for every admissible radius choice.
The resulting tree has constant distortion from every vertex of the protected path to every vertex in the instance. This gives the cost scale needed in the analysis: an edge first separated in an instance can be charged according to its endpoint’s distance to that instance’s protected path. Once this distance becomes comparable to the edge length, the edge’s entire remaining tree distance is paid immediately. The geometry also permits every surviving off-path vertex to compete as a source, which is essential for the local comparisons below.
Local scales and a carried threshold. For a vertex v v v in an instance X X X , let t X ( v ) t_X(v) t X ( v ) be its physical distance to the protected path, measured at entry to the instance. We fix an allowance R v = ε t X ( v ) R_v = \varepsilon t_X(v) R v = ε t X ( v ) , for a sufficiently small absolute constant ε \varepsilon ε . These allowances stay fixed while the parent is carved. Their Lipschitz property supplies meaningful comparisons between nearby vertices even after residual distances have increased.
We also fix a positive vertex measure once for the entire construction. Write b e = a e ℓ e b_e = a_e\ell_e b e = a e ℓ e and B = ∑ e b e B = \sum_e b_e B = ∑ e b e , and assume B > 0 B > 0 B > 0 ; the zero-demand case is immediate. Assign each vertex the baseline weights of its incident edges and an additional mass B / n B/n B / n . Its total mass is 3 B 3B 3 B , and its entropy is at most 3 B log n 3B\log n 3 B log n . This single entropy budget will pay for all charged edges throughout the recursion. In a current residual, the core of v v v is a small physical ball of radius δ R v \delta R_v δ R v ; denote its mass by m v m_v m v and the fixed parent mass by M M M .
The personal rate of a vertex is proportional to K v / R v K_v/R_v K v / R v , where K v = ⌈ log 2 ( 2 M / m v ) ⌉ K_v = \lceil\log_2(2M/m_v)\rceil K v = ⌈ log 2 ( 2 M / m v )⌉ . A vertex is eligible when this rate is at most the current threshold Λ \Lambda Λ . Among eligible vertices we maximize h ( v ) − R v h(v)-R_v h ( v ) − R v , where h ( v ) h(v) h ( v ) is its current root distance. The threshold doubles when no candidate is eligible. Each child starts with twice the threshold used at its own source selection. The priority rule ensures that an eligible captured vertex has an allowance comparable to the source’s allowance and a child root distance at most ε t X ( v ) \varepsilon t_X(v) ε t X ( v ) . Thus eligibility gives strong contraction; ineligibility gives a mass inequality that will instead pay for the capture.
Summing over physical scales. The natural charge at a capture of v v v is Λ t X ( v ) \Lambda t_X(v) Λ t X ( v ) . Although t X ( v ) t_X(v) t X ( v ) need not decrease at every capture, it always lies within a constant factor of its running minimum. The reason is that current root heights decrease along a recursive history, and every capture resets the next root height to at most a constant times the present distance to the protected path. This also controls any later increase in that distance.
We group captures into dyadic bands of the running minimum. Within one band, the physical scales are comparable while the carried thresholds grow geometrically. The total charge in that band is therefore bounded by a constant times its last charge. For these representatives, place an inner and an outer physical ball around the vertex. Their radii are fixed constant multiples of the running minimum. After a constant number of occupied bands, the later outer ball is contained in the earlier inner ball. Residuals are nested, and intrinsic distances cannot decrease under restriction. The logarithms of these local mass ratios consequently telescope, independently of the number of recursive levels or the metric aspect ratio.
Far captures and inherited rates. Call a capture far when the source allowance is at most 1 / χ 1/\chi 1/ χ times the captured vertex’s allowance, for a fixed constant χ ≥ 4 \chi\ge4 χ ≥ 4 . Such a vertex must be ineligible: the comparison for eligible captures rules out this discrepancy. Using K v K_v K v , within constant factors of log ( 2 M / m v ) \log(2M/m_v) log ( 2 M / m v ) , the two eligibility comparisons bound the charge by a difference of logarithms. Their common parent-mass term cancels, giving
Δ t X ( v ) ≤ O ( 1 ) log m p m v ≤ O ( 1 ) log μ ( O v ) μ ( I v ) . \Delta t_X(v) \le O(1)\log\frac{m_p}{m_v} \le O(1)\log\frac{\mu(O_v)}{\mu(I_v)}. Δ t X ( v ) ≤ O ( 1 ) log m v m p ≤ O ( 1 ) log μ ( I v ) μ ( O v ) .
Here the constants absorb 1 / ε 1/\varepsilon 1/ ε , and O v , I v O_v,I_v O v , I v are the outer and inner balls just described. Geometry places the source core inside O v O_v O v , while I v I_v I v lies in the endpoint core. The local entropy telescope therefore pays for every far capture, including captures deep inside a cell.
Near captures require additional care because an inherited threshold may greatly exceed a vertex’s present personal rate. We begin a new run at the first capture, whenever the vertex is ineligible, or when the threshold has increased inside its current parent. At these seed captures, the threshold is controlled by the vertex’s personal rate. The first subsequent charge is at most a constant multiple of the seed charge; later charges decrease geometrically because eligibility contracts the scale faster than threshold inheritance increases it. The run’s total charge is therefore bounded by a constant times its seed charge.
A seed bound initially contains an additive constant from log ( 2 M / m v ) \log(2M/m_v) log ( 2 M / m v ) . Paying that constant once per occupied band would destroy the desired estimate. We prove that a seed with sufficiently small entropy is deep inside its cell and cannot be charged for separating an active near edge. Its inner ball contains almost all the parent mass, while the seed bound gives the source core a constant fraction of that mass. The two balls must therefore intersect. The next distance to the protected path is then so small that a further capture requires a fresh threshold increase. Such a seed cannot sponsor a later inherited-rate charge either. The remaining seeds have enough entropy to absorb the additive constant.
Accounting for near vertices together. At a near seed, the rate estimate still involves the whole parent mass. Split its logarithm into a local inner-to-outer mass ratio and the remaining parent-to-outer ratio. The first part has already been handled by nested balls. To bound the second, consider one actual child cell and let N N N and F F F be the masses of its near and far vertices, respectively. All near vertices lie within each other’s outer balls. Writing V = N + F V=N+F V = N + F , their total remaining charge is bounded using the elementary inequality
N log ( M / N ) ≤ V log ( M / V ) + F . N\log(M/N) \le V\log(M/V)+F. N log ( M / N ) ≤ V log ( M / V ) + F .
The first term belongs to the entropy of the parent partition and telescopes over the recursive partition tree. Every far vertex has a constant lower bound on its capture charge, so the second term is paid by the far-capture estimate. This aggregate step is what removes the global parent-mass term from the local analysis.
Boundary events and the upper endpoint of a radius. For analysis, a radius is drawn from a truncated exponential law on the source’s allowance interval. The exact separation identity consists of a capture term and a truncation remainder. The capture term can be restricted to vertices whose entry thresholds are near the interval; a deep capture contributes nothing to it. For far opportunities, the remainder is paid by source-core masses: their normalized squares sum to at most a constant within a parent, and the same nested-ball argument controls their carried rates across parents. This covers earlier potential separations even when the actual capture occurs later. For near edges entering in the last exponential length of the interval, we instead charge the binary entropy of the actual cell and its complement. Both kinds of remainder therefore fit within the original entropy budget.
Exact constructions and the size of the distribution
A deterministic radius from a local inequality. The preceding accounting also gives a deterministic weighted oracle. At a fixed source step, let C ( r ) C(r) C ( r ) bound the weighted distance cost of edges separated by radius r r r . The one-step exponential analysis shows that its expectation is at most the sum of three credits: the charge of actual captures, the binary entropy of the split, and the far truncation remainder. Each credit has a bound that holds for every complete admissible execution. It therefore suffices to choose a radius satisfying this same inequality locally and then sum the resulting inequalities.
Such a radius can be found exactly. The cell, its portal, and all edge-entry events are constant between consecutive vertex gauges. Enumerating the rational midpoints of these intervals gives a finite list containing a satisfactory choice. Binary entropy is replaced by a rational upper bound within an absolute constant factor, and eligibility uses only the integers K v K_v K v . No probability needs to be evaluated. A source becomes permanently protected after it is selected, so there are at most n − 1 n-1 n − 1 source selections. Fixed denominators for physical distances and polynomially bounded threshold exponents give polynomial running time for arbitrary rational input.
From weighted responses to an explicit distribution. We use the response-oracle viewpoint of Charikar et al., Räcke, and Andersen–Feige [10 , 22 , 3 ] . For either application, a tree supplies an m m m -dimensional nonnegative payoff vector, and the weighted oracle returns a vector whose cost against any positive price distribution is at most ρ \rho ρ . An adaptive-step construction produces a mixture with every expected coordinate at most ( 1 + ϑ ) ρ (1+\vartheta)\rho ( 1 + ϑ ) ρ , using O ( m log ( 2 m ) / ϑ 2 ) O(m\log(2m)/\vartheta^2) O ( m log ( 2 m ) / ϑ 2 ) calls.
The update size is inversely proportional to the largest coordinate of the current response. Thus even a response with very large individual stretch or congestion makes only a bounded coordinate increment. Prices are integer powers of a rational base, with exponents determined by rounded accumulated coordinates. A potential bound limits these exponents uniformly. Normalizing the stopping horizon by 1 / ρ 1/\rho 1/ ρ makes the call count independent of the numerical size of ρ \rho ρ . All probabilities and queries consequently have polynomial encodings. Finally, affine dependence reduces the mixture to at most m + 1 m+1 m + 1 trees while preserving its mean. One additional rational linear-algebra step either removes an affine dependence or decreases all mean coordinates together until a probability vanishes. This gives support at most m m m without another oracle call.
Why substantially smaller support can fail. For fundamental-cut congestion, a quadratic dependence on the number of vertices can be necessary even when the best congestion tends to one. Take two sets of k k k vertices, put a capacity-k 2 k^2 k 2 star inside each (Figure 3 ), and add every cross edge with capacity one. Every spanning tree contains a cross edge, and its fundamental cut has capacity at least k 2 k^2 k 2 . A tree given probability q T q_T q T therefore contributes at least q T k 2 q_Tk^2 q T k 2 to some original edge’s expected congestion. Any distribution with congestion at most ρ \rho ρ must satisfy
Figure 3. The support lower bound. All k 2 k^2 k 2 cross edges are present; only representative leaves and cross edges are drawn. A canonical tree retains both heavy stars and one bridge. That bridge carries the capacity-k 2 k^2 k 2 fundamental cut A ∣ B A \mid B A ∣ B . A mixture spreads this load among different original cross edges.
q T ≤ ρ / k 2 , ∣ supp q ∣ ≥ ⌈ k 2 / ρ ⌉ . q_T \le\rho/k^2,\qquad|\operatorname{supp} q| \ge\lceil k^2/\rho\rceil. q T ≤ ρ / k 2 , ∣ supp q ∣ ≥ ⌈ k 2 / ρ ⌉ .
Nevertheless, retaining both heavy stars and choosing one cross edge uniformly gives congestion at most 1 + 2 ( k − 1 ) / k 2 1+2(k-1)/k^2 1 + 2 ( k − 1 ) / k 2 . Balanced sets of cross edges attain the support lower bound for every ρ ≥ ( k 2 + 2 k − 2 ) / ( k 2 − 1 ) \rho\ge(k^2+2k-2)/(k^2-1) ρ ≥ ( k 2 + 2 k − 2 ) / ( k 2 − 1 ) . Thus the support requirement reflects the large congestion of an individual sampled tree, despite the small congestion of the fixed mixture. It applies to the original-edge load guarantee used for routing. On the same family, a single fundamental-cut tree already approximates every cut within a constant factor. This explains why the O ( n log n ) O(n\log n) O ( n log n ) support bound for abstract metric-tree embeddings [10 ] does not carry over to the stronger congestion requirement.
Preliminaries and the recursive construction
We delete loops at the outset and allow parallel edges. All lengths in the metric theorem are strictly positive. Unless stated otherwise, capacities are also strictly positive; Section 9 explains the extension to nonnegative capacities. For a vertex set A A A , δ G ( A ) \delta_G(A) δ G ( A ) is the set of original edges with exactly one endpoint in A A A , and c G ( δ G ( A ) ) c_G(\delta_G(A)) c G ( δ G ( A )) is their total capacity. All graph distances are intrinsic to the specified induced graph. Balls are closed vertex balls. Logarithms without a subscript are natural. Every tie is broken by a fixed ordering.
The weighted theorem takes a positive length ℓ e \ell_e ℓ e and a nonnegative demand a e a_e a e on each edge. We call ℓ \ell ℓ the physical length vector : every edge of the output tree keeps this length. Put
b e = a e ℓ e , B = ∑ e b e , C a ( T ) = ∑ e = u v a e d T , ℓ ( u , v ) . b_e=a_e\ell_e,\qquad B=\sum_e b_e,\qquad C_a(T)=\sum_{e=uv}a_e d_{T,\ell}(u,v). b e = a e ℓ e , B = e ∑ b e , C a ( T ) = e = uv ∑ a e d T , ℓ ( u , v ) .
If B = 0 B=0 B = 0 , any spanning tree suffices. For B > 0 B>0 B > 0 , the proof fixes the positive vertex measure
μ ( v ) = ∑ e ∋ v b e + B / n . \mu(v)=\sum_{e\ni v}b_e+B/n. μ ( v ) = e ∋ v ∑ b e + B / n .
The measure is only restricted under recursion; it is never recomputed from a child graph. Its total mass is 3 B 3B 3 B , and
H 0 : = ∑ v μ ( v ) log 3 B μ ( v ) ≤ 3 B log n . H_0:=\sum_v\mu(v)\log\frac{3B}{\mu(v)}\le3B\log n. H 0 := v ∑ μ ( v ) log μ ( v ) 3 B ≤ 3 B log n .
The geometric lemmas do not depend on this particular measure. The entropy argument first applies to any fixed positive measure and is then specialized to the one above.
Each recursive instance consists of an induced connected graph X X X , a root x x x , and a physical shortest path P X P_X P X starting at x x x . We call P X P_X P X the protected path . The initial instance is G G G , with an arbitrary root and the one-vertex path at that root. Within an instance, the algorithm repeatedly removes a connected petal disjoint from P X P_X P X . A removed petal becomes a child instance and has one designated original portal edge joining it to the residual graph. Once only P X P_X P X remains, the algorithm retains its edges and recurses on the children. The union of these paths, recursive trees, and portal edges is the output tree.
The next section defines petals and proves the geometric invariants for every admissible source and radius. The source rule and inherited threshold are specified in Section 5; a child’s threshold is twice the threshold at the selection of its own source. Section 6 gives the finite deterministic radius rule. Temporary lengths used to define a petal are reset in each instance: edges on its protected path have length ℓ e / 2 \ell_e/2 ℓ e /2 , and all other edges have length ℓ e \ell_e ℓ e . In particular, lengths are never repeatedly halved along recursion.
Protected geodesics and recursive spanning trees
The decomposition preserves one physical geodesic in each recursive instance. Temporarily shortening this path makes it possible to attach the remaining vertices through original graph edges while controlling their distances to every vertex of the preserved path. Throughout this section, physical lengths are the original positive lengths ℓ \ell ℓ . The parameters used later satisfy
0 < ε ≤ 1 8 , 0 < δ ≤ 1 256 , χ ≥ 4 , A 0 ≥ 54 χ . 0 < \varepsilon\le\frac{1}{8}, \qquad0 < \delta\le\frac{1}{256}, \qquad\chi\ge4, \qquad A_0 \ge54\chi. 0 < ε ≤ 8 1 , 0 < δ ≤ 256 1 , χ ≥ 4 , A 0 ≥ 54 χ .
The structural and distance bounds hold for arbitrary admissible radii. The scale consequences use fixed allowances proportional to the parent distance from the protected path.
One recursive instance. Let X X X be a connected induced graph, let x ∈ X x \in X x ∈ X be its root, and let P P P be a physical shortest path starting at x x x . Assign current length w e = ℓ e / 2 w_e = \ell_e/2 w e = ℓ e /2 to the edges of P P P and w e = ℓ e w_e = \ell_e w e = ℓ e to all other edges. Every prefix of P P P remains current-shortest: a competing route has current length at least half its physical length. Choose a current shortest-path tree S S S rooted at x x x that contains P P P , and set h ( v ) = d X , w ( x , v ) h(v) = d_{X,w}(x,v) h ( v ) = d X , w ( x , v ) . Such a tree exists by fixing the tight predecessors on P P P and then choosing a tight predecessor for every remaining vertex. Positive lengths make this predecessor relation acyclic.
Initially the residual is Y = X Y = X Y = X . We shall maintain that Y Y Y is ancestor-closed in S S S , so that every surviving root path remains and the current root distances in Y Y Y are still h h h . To form a cell, select p ∈ Y ∖ P p \in Y \setminus P p ∈ Y ∖ P . Write Q p = P S ( x , p ) Q_p = P_S(x,p) Q p = P S ( x , p ) , let z p z_p z p be its last vertex on P P P , and set s p = h ( p ) − h ( z p ) s_p = h(p) - h(z_p) s p = h ( p ) − h ( z p ) . Choose a radius 0 < r p < s p 0 < r_p < s_p 0 < r p < s p and put c p = h ( p ) − r p c_p = h(p) - r_p c p = h ( p ) − r p . Define the entry gauge and the cell by
D p ( v ) = h ( p ) − 2 h ( v ) + min q ∈ Q p { h ( q ) + 2 d Y , w ( q , v ) } , X p = { v ∈ Y : D p ( v ) < r p } . (3) D_p(v) = h(p) - 2h(v) + \min_{q \in Q_p}\{h(q) + 2d_{Y,w}(q,v)\}, \qquad X_p = \{v \in Y : D_p(v) < r_p\}. \tag*{(3)} D p ( v ) = h ( p ) − 2 h ( v ) + q ∈ Q p min { h ( q ) + 2 d Y , w ( q , v )} , X p = { v ∈ Y : D p ( v ) < r p } . (3)
Let a p a_p a p be the first vertex of Q p Q_p Q p whose height exceeds c p c_p c p , and let b p b_p b p be its predecessor. Remove X p X_p X p from Y Y Y and retain the portal edge a p b p a_pb_p a p b p . Continue until only P P P remains.
After completing this parent decomposition, recurse on each induced child X [ X p ] X[X_p] X [ X p ] , with root a p a_p a p and protected path P S ( a p , p ) P_S(a_p,p) P S ( a p , p ) . In the child the current lengths are reset: exactly its new protected path is halved, and every other edge receives its original physical length. The output consists of the child trees, all portal edges, and the edges of P P P . In particular, no metric edge is introduced.
Lemma 4.1 (Structure of a petal). The gauge D p D_p D p is nonnegative, satisfies D p ( p ) = 0 D_p(p) = 0 D p ( p ) = 0 , and is 4-Lipschitz in physical residual distance. Each cell X p X_p X p is connected, contains a p a_p a p and p p p , and avoids P P P ; its portal endpoint b p b_p b p remains in the residual. Removing the cell preserves ancestor closure in S S S . The path P S ( a p , p ) P_S(a_p,p) P S ( a p , p ) is physically shortest in the induced child.
If R p = ε d X , ℓ ( p , P ) R_p = \varepsilon d_{X,\ell}(p,P) R p = ε d X , ℓ ( p , P ) and r p ∈ [ R p / 2 , R p ] r_p \in[R_p/2,R_p] r p ∈ [ R p /2 , R p ] , then the radius is admissible and the cell contains the entire physical core B Y , ℓ ( p , δ R p ) B_{Y,\ell}(p,\delta R_p) B Y , ℓ ( p , δ R p ) .
Proof. For every witness q ∈ Q p q \in Q_p q ∈ Q p , the triangle inequality gives
h ( p ) + h ( q ) + 2 d Y , w ( q , v ) − 2 h ( v ) ≥ h ( p ) − h ( q ) ≥ 0. h(p) + h(q) + 2d_{Y,w}(q,v) - 2h(v) \ge h(p) - h(q) \ge0. h ( p ) + h ( q ) + 2 d Y , w ( q , v ) − 2 h ( v ) ≥ h ( p ) − h ( q ) ≥ 0.
The witness q = p q = p q = p gives D p ( p ) = 0 D_p(p) = 0 D p ( p ) = 0 . Both h h h and each residual distance function are 1-Lipschitz in physical residual distance. Therefore, the expression in (3) , and hence its minimum, is therefore 4-Lipschitz. On the source path, ancestor closure and the shortest-path property give
d Y , w ( q , v ) = ∣ h ( q ) − h ( v ) ∣ ( q , v ∈ Q p ) , D p ( v ) = h ( p ) − h ( v ) ( v ∈ Q p ) . d_{Y,w}(q,v) = \lvert h(q) - h(v)\rvert\quad(q,v \in Q_p), \qquad D_p(v) = h(p) - h(v) \quad(v \in Q_p). d Y , w ( q , v ) = ∣ h ( q ) − h ( v )∣ ( q , v ∈ Q p ) , D p ( v ) = h ( p ) − h ( v ) ( v ∈ Q p ) .
Thus precisely the source-path vertices above the cutoff enter the cell. In particular, a p ∈ X p a_p \in X_p a p ∈ X p , b p ∉ X p b_p \notin X_p b p ∈ / X p , and the segment from a p a_p a p to p p p lies in X p X_p X p .
For a captured vertex v v v , choose a minimizing witness q q q and define its excess
η = h ( q ) + d Y , w ( q , v ) − h ( v ) ≥ 0. 2 η < h ( q ) − c p (4) \eta= h(q) + d_{Y,w}(q,v) - h(v) \ge0. \qquad2\eta< h(q) - c_p \tag*{(4)} η = h ( q ) + d Y , w ( q , v ) − h ( v ) ≥ 0. 2 η < h ( q ) − c p (4)
Consequently q q q lies on the captured segment a p p a_p p a p p . Every vertex of a shortest q q q –v v v route is captured as well. Indeed, moving from v v v toward q q q decreases the distance to q q q by exactly the traversed length, whereas the height can decrease by at most that length. This proves connectivity. Along an outward edge of S S S , the height increases by the edge’s current length while any witness distance increases by at most that length. The gauge cannot increase. The cell therefore contains every surviving S S S -descendant of each of its vertices, proving ancestor closure of the new residual.
To see that P P P is avoided, fix v ∈ P v \in P v ∈ P . For a witness q q q on P P P , the one-dimensional distance formula gives
h ( q ) + 2 d Y , w ( q , v ) − 2 h ( v ) ≥ − h ( q ) ≥ − h ( z p ) . h(q) + 2d_{Y,w}(q,v) - 2h(v) \ge-h(q) \ge-h(z_p). h ( q ) + 2 d Y , w ( q , v ) − 2 h ( v ) ≥ − h ( q ) ≥ − h ( z p ) .
For a witness beyond z p z_p z p , physical geodesicity of P P P and ℓ ≤ 2 w \ell\le2w ℓ ≤ 2 w give
2 ∣ h ( v ) − h ( z p ) ∣ ≤ h ( q ) − h ( z p ) + 2 d Y , w ( q , v ) . 2\lvert h(v) - h(z_p)\rvert\le h(q) - h(z_p) + 2d_{Y,w}(q,v). 2 ∣ h ( v ) − h ( z p )∣ ≤ h ( q ) − h ( z p ) + 2 d Y , w ( q , v ) .
The same lower bound − h ( z p ) -h(z_p) − h ( z p ) follows. Hence D p ( v ) ≥ h ( p ) − h ( z p ) = s p > r p D_p(v) \ge h(p) - h(z_p) = s_p > r_p D p ( v ) ≥ h ( p ) − h ( z p ) = s p > r p . All witness routes inside the cell consequently avoid the shortened edges of P P P , and their current lengths equal their physical lengths. The new protected path is a subpath of a current shortest path. Since the entire induced child contains no shortened edge, this path is physically shortest in the child.
Finally, the source tail from z p z_p z p to p p p has physical length s p s_p s p , so d X , ℓ ( p , P ) ≤ s p d_{X,\ell}(p,P) \le s_p d X , ℓ ( p , P ) ≤ s p and R p < s p R_p < s_p R p < s p . On the prescribed core the Lipschitz bound gives D p ( v ) ≤ 4 δ R p < R p / 2 D_p(v) \le4\delta R_p < R_p/2 D p ( v ) ≤ 4 δ R p < R p /2 , proving inclusion for every allowed radius.
Lemma 4.2 (Tree construction and root control). Every adaptive sequence of sources and admissible radii produces an original-edge spanning tree T T T containing P P P . For each captured vertex v ∈ X p v \in X_p v ∈ X p , there is a physical route from v v v to b p b_p b p , contained in X p X_p X p except for its portal edge, of length at most 5 3 ( h ( v ) − h ( b p ) ) \frac{5}{3}(h(v) - h(b_p)) 3 5 ( h ( v ) − h ( b p )) . The output tree satisfies
d T ( v , b p ) ≤ 2 ( h ( v ) − h ( b p ) ) , d T ( x , v ) ≤ 2 h ( v ) . (5) d_T(v,b_p) \le2(h(v) - h(b_p)), \qquad d_T(x,v) \le2h(v). \tag*{(5)} d T ( v , b p ) ≤ 2 ( h ( v ) − h ( b p )) , d T ( x , v ) ≤ 2 h ( v ) . (5) There are at most ∣ X ∣ − 1 |X| - 1 ∣ X ∣ − 1 source selections throughout the recursion from X X X . Moreover, each recursive tree T Z T_Z T Z remains a connected subtree of the final output: the final tree path between vertices of Z Z Z is exactly their path in T Z T_Z T Z .
Proof. Every child avoids the nonempty path P P P and is therefore smaller than its parent. Recursion terminates by induction on the number of vertices. At the parent level each portal points to a later cell or to P P P . The child trees, portals, and P P P consequently form a connected graph. If there are k k k cells, its number of edges is
∑ p ( ∣ X p ∣ − 1 ) + k + ( ∣ V ( P ) ∣ − 1 ) = ∣ X ∣ − 1 , \sum_p (|X_p| - 1) + k + (|V(P)| - 1) = |X| - 1, p ∑ ( ∣ X p ∣ − 1 ) + k + ( ∣ V ( P ) ∣ − 1 ) = ∣ X ∣ − 1 ,
so it is a tree.
Fix v ∈ X p v \in X_p v ∈ X p and its witness q q q from (4) . Since d Y , w ( q , v ) ≥ h ( q ) − h ( v ) d_{Y,w}(q,v) \ge h(q) - h(v) d Y , w ( q , v ) ≥ h ( q ) − h ( v ) , we have h ( q ) − h ( v ) ≤ η / 2 h(q) - h(v) \le\eta/2 h ( q ) − h ( v ) ≤ η /2 . Combining this with the capture inequality gives η < 2 3 ( h ( v ) − c p ) \eta< \frac{2}{3}(h(v) - c_p) η < 3 2 ( h ( v ) − c p ) . The route that follows the source path from b p b_p b p to q q q and then a shortest witness route to v v v is physical, by Lemma 4.1, and has length
h ( q ) − h ( b p ) + d Y , w ( q , v ) = h ( v ) − h ( b p ) + η ≤ 5 3 ( h ( v ) − h ( b p ) ) . h(q) - h(b_p) + d_{Y,w}(q,v) = h(v) - h(b_p) + \eta\le\frac{5}{3}(h(v) - h(b_p)). h ( q ) − h ( b p ) + d Y , w ( q , v ) = h ( v ) − h ( b p ) + η ≤ 3 5 ( h ( v ) − h ( b p )) .
For the tree bound, let h ^ \hat{h} h ^ be the child’s current root distance. Its new shortening gives
h ^ ( v ) ≤ 1 2 ( h ( q ) − h ( a p ) ) + d Y , w ( q , v ) . \hat{h}(v) \le\frac{1}{2}\left(h(q)-h(a_p)\right)+d_{Y,w}(q,v). h ^ ( v ) ≤ 2 1 ( h ( q ) − h ( a p ) ) + d Y , w ( q , v ) .
Inductively the child tree has root distances at most 2 h ^ 2\hat{h} 2 h ^ . Including the portal and using (4) yields
d T ( v , b p ) ≤ h ( a p ) − h ( b p ) + 2 h ^ ( v ) ≤ 2 h ( v ) − h ( q ) − h ( b p ) + 2 η < 2 h ( v ) − c p − h ( b p ) ≤ 2 ( h ( v ) − h ( b p ) ) . \begin{aligned}
d_T(v,b_p) &\le h(a_p)-h(b_p)+2\hat{h}(v)\\
&\le2h(v)-h(q)-h(b_p)+2\eta\\
&< 2h(v)-c_p-h(b_p)\le2\left(h(v)-h(b_p)\right).
\end{aligned} d T ( v , b p ) ≤ h ( a p ) − h ( b p ) + 2 h ^ ( v ) ≤ 2 h ( v ) − h ( q ) − h ( b p ) + 2 η < 2 h ( v ) − c p − h ( b p ) ≤ 2 ( h ( v ) − h ( b p ) ) .
The root bound now follows backwards through the parent cells. On P P P it holds with equality, since the retained physical lengths are twice the current lengths.
For the source count, a selected source becomes a vertex of its child’s protected path and can never again be selected in any descendant. Different branches have disjoint vertex sets, so all selected sources are distinct initially unprotected vertices. Finally, every recursive tree is retained as a connected subgraph of the final tree. Uniqueness of tree paths proves the last claim. Several outside branches may attach to this subtree; connectivity and acyclicity alone give the required path property.
Lemma 4.3 (Protection of the entire geodesic). For every q ∈ P q \in P q ∈ P and u ∈ X u \in X u ∈ X ,
d T ( q , u ) ≤ 34 d X , ℓ ( q , u ) . (6) d_T(q,u) \le34d_{X,\ell}(q,u). \tag*{(6)} d T ( q , u ) ≤ 34 d X , ℓ ( q , u ) . (6) In particular, writing t X ( u ) = d X , ℓ ( u , P ) t_X(u)=d_{X,\ell}(u,P) t X ( u ) = d X , ℓ ( u , P ) , every original edge u v uv uv of the instance satisfies
d T ( u , v ) ≤ 34 ( 2 t X ( u ) + ℓ u v ) . (7) d_T(u,v) \le34(2t_X(u)+\ell_{uv}). \tag*{(7)} d T ( u , v ) ≤ 34 ( 2 t X ( u ) + ℓ uv ) . (7) These inequalities also hold for the corresponding paths in the final tree of any larger instance containing this recursive call.
Proof. Follow the same-parent portal chain from u u u until its first vertex b ∈ P b \in P b ∈ P . Put L = h ( u ) − h ( b ) L=h(u)-h(b) L = h ( u ) − h ( b ) and D = d X , ℓ ( u , q ) D=d_{X,\ell}(u,q) D = d X , ℓ ( u , q ) . The tree and physical comparison bounds of Lemma 4.2 telescope along this chain: the tree route from u u u to b b b has length at most 2 L 2L 2 L , and there is a physical comparison route of length at most 5 L / 3 5L/3 5 L /3 . Since h h h is 1-Lipschitz in physical distance and P P P is physically geodesic,
2 ( L − D ) ≤ 2 ∣ h ( q ) − h ( b ) ∣ = ℓ ( P [ q , b ] ) ≤ D + 5 L / 3. 2(L-D) \le2\lvert h(q)-h(b)\rvert= \ell(P[q,b]) \le D+5L/3. 2 ( L − D ) ≤ 2 ∣ h ( q ) − h ( b )∣ = ℓ ( P [ q , b ]) ≤ D + 5 L /3.
Thus L ≤ 9 D L \le9D L ≤ 9 D and ℓ ( P [ q , b ] ) ≤ 16 D \ell(P[q,b]) \le16D ℓ ( P [ q , b ]) ≤ 16 D . The retained path from q q q to b b b followed by the tree route to u u u has length at most 16 D + 18 D = 34 D 16D+18D=34D 16 D + 18 D = 34 D .
For an edge u v uv uv , choose a nearest protected-path vertex q q q to u u u . Then d X , ℓ ( q , v ) ≤ t X ( u ) + ℓ u v d_{X,\ell}(q,v) \le t_X(u)+\ell_{uv} d X , ℓ ( q , v ) ≤ t X ( u ) + ℓ uv . Applying (6) to both endpoints and using the tree triangle inequality proves (7). The final assertion is the connected-subtree property of Lemma 4.2 .
Lemma 4.4 (Ambient locality and eligible contraction). Suppose p p p captures u u u . With the parent distance t X ( u ) = d X , ℓ ( u , P ) t_X(u)=d_{X,\ell}(u,P) t X ( u ) = d X , ℓ ( u , P ) fixed before any residual removals,
h c h i l d ( u ) ≤ h ( u ) − c p , (8) h_{\mathrm{child}}(u) \le h(u)-c_p, \tag*{(8)} h child ( u ) ≤ h ( u ) − c p , (8) d X p , ℓ ( p , u ) ≤ 3 t X ( u ) − s p + 4 r p , (9) d_{X_p,\ell}(p,u) \le3t_X(u)-s_p+4r_p, \tag*{(9)} d X p , ℓ ( p , u ) ≤ 3 t X ( u ) − s p + 4 r p , (9) h c h i l d ( u ) ≤ 3 t X ( u ) − s p + 3 r p . (10) h_{\mathrm{child}}(u) \le3t_X(u)-s_p+3r_p. \tag*{(10)} h child ( u ) ≤ 3 t X ( u ) − s p + 3 r p . (10) :::`
In particular, these bounds use the original parent distance to P P P , not a residual distance that may increase during the carving.
Suppose further that R v = ϵ t X ( v ) R_v=\epsilon t_X(v) R v = ϵ t X ( v ) is fixed throughout the parent, r p ∈ [ R p / 2 , R p ] r_p \in[R_p/2,R_p] r p ∈ [ R p /2 , R p ] , and p p p maximizes h ( v ) − R v h(v)-R_v h ( v ) − R v over an eligible set containing u u u . Then
h c h i l d ( u ) ≤ R u , d X p , ℓ ( p , u ) ≤ R p + R u , 1 − ϵ 1 + ϵ ≤ R p R u ≤ 1 + ϵ 1 − ϵ . (11) h_{\mathrm{child}}(u) \le R_u,\qquad d_{X_p,\ell}(p,u) \le R_p+R_u,\qquad\frac{1-\epsilon}{1+\epsilon} \le\frac{R_p}{R_u} \le\frac{1+\epsilon}{1-\epsilon}. \tag*{(11)} h child ( u ) ≤ R u , d X p , ℓ ( p , u ) ≤ R p + R u , 1 + ϵ 1 − ϵ ≤ R u R p ≤ 1 − ϵ 1 + ϵ . (11)
Proof. Put A = h ( u ) − h ( p ) A=h(u)-h(p) A = h ( u ) − h ( p ) and let q 0 q_0 q 0 be a minimizing witness for u u u , with excess η 0 \eta_0 η 0 as in (4) . The source–witness–u u u route is entirely inside the cell and has physical length
d 0 = A + 2 ( h ( p ) − h ( q 0 ) ) + η 0 ≤ A + 2 ( h ( p ) − h ( q 0 ) ) + 2 η 0 < A + 2 r p . d_0=A+2\bigl(h(p)-h(q_0)\bigr)+\eta_0 \le A+2\bigl(h(p)-h(q_0)\bigr)+2\eta_0 < A+2r_p. d 0 = A + 2 ( h ( p ) − h ( q 0 ) ) + η 0 ≤ A + 2 ( h ( p ) − h ( q 0 ) ) + 2 η 0 < A + 2 r p .
Choose a nearest vertex q ∈ P q \in P q ∈ P to u u u in the original parent metric, and write t = t X ( u ) t=t_X(u) t = t X ( u ) . Concatenate the physical source tail z p p z_p p z p p , the preceding in-cell route, and an ambient physical shortest u u u –q q q path. The latter path need not survive in Y Y Y ; it is used only to compare lengths in X X X . Physical geodesicity of P P P and h ( q ) ≥ h ( u ) − t h(q) \ge h(u)-t h ( q ) ≥ h ( u ) − t give
2 ( A + s p − t ) ≤ 2 ∣ h ( q ) − h ( z p ) ∣ ≤ s p + A + 2 r p + t . 2(A+s_p-t) \le2\lvert h(q)-h(z_p)\rvert\le s_p+A+2r_p+t. 2 ( A + s p − t ) ≤ 2 ∣ h ( q ) − h ( z p )∣ ≤ s p + A + 2 r p + t .
Therefore A + s p ≤ 3 t + 2 r p A+s_p \le3t+2r_p A + s p ≤ 3 t + 2 r p , and the bound on d 0 d_0 d 0 proves (9) .
The same witness route, with the child’s new shortening, gives
h c h i l d ( u ) ≤ h ( u ) − 1 2 h ( q 0 ) − 1 2 h ( a p ) + η 0 < h ( u ) − 1 2 c p − 1 2 h ( a p ) ≤ h ( u ) − c p . \begin{aligned}
h_{\mathrm{child}}(u) \le h(u)-\frac{1}{2}h(q_0)-\frac{1}{2}h(a_p)+\eta_0 \\
&< h(u)-\frac{1}{2}c_p-\frac{1}{2}h(a_p) \le h(u)-c_p.
\end{aligned} h child ( u ) ≤ h ( u ) − 2 1 h ( q 0 ) − 2 1 h ( a p ) + η 0 < h ( u ) − 2 1 c p − 2 1 h ( a p ) ≤ h ( u ) − c p .
This proves (8) . Substituting h ( u ) − c p = A + r p h(u)-c_p=A+r_p h ( u ) − c p = A + r p proves (10) .
For an eligible u u u , source maximality gives
c p ≥ h ( p ) − R p ≥ h ( u ) − R u , A ≤ R u − R p . c_p \ge h(p)-R_p \ge h(u)-R_u,\qquad A \le R_u-R_p. c p ≥ h ( p ) − R p ≥ h ( u ) − R u , A ≤ R u − R p .
The child-height bound follows immediately, while the canonical in-cell route has length at most A + 2 r p ≤ R u + R p A+2r_p \le R_u+R_p A + 2 r p ≤ R u + R p . Finally, t X t_X t X is 1-Lipschitz in the ambient physical metric, so
∣ R p − R u ∣ ≤ ϵ d X , ℓ ( p , u ) ≤ ϵ ( R p + R u ) . \lvert R_p-R_u\rvert\le\epsilon d_{X,\ell}(p,u) \le\epsilon(R_p+R_u). ∣ R p − R u ∣ ≤ ϵ d X , ℓ ( p , u ) ≤ ϵ ( R p + R u ) .
Rearranging yields both allowance ratios.
The scale information needed by the entropy analysis follows from these geometric bounds. In particular, the factor-three comparison below does not compound with recursion depth: current root heights decrease from one instance to the next.
Lemma 4.5 (Consequences for scales and boundary charges). Use fixed allowances R v = ϵ t X ( v ) R_v=\epsilon t_X(v) R v = ϵ t X ( v ) , allowed radii r p ∈ [ R p / 2 , R p ] r_p \in[R_p/2,R_p] r p ∈ [ R p /2 , R p ] , and the eligible priority rule of Lemma 4.4. Every capture satisfies
d Y , ℓ ( p , u ) ≤ 3 t X ( u ) , t X ( p ) < 3 t X ( u ) , h c h i l d ( u ) ≤ min { h ( u ) , 3 t X ( u ) } , t X ( u ) ≤ h ( u ) . (12) d_{Y,\ell}(p,u) \le3t_X(u),\qquad t_X(p)<3t_X(u),\qquad h_{\mathrm{child}}(u)\le\min\{h(u),3t_X(u)\},\qquad t_X(u)\le h(u). \tag*{(12)} d Y , ℓ ( p , u ) ≤ 3 t X ( u ) , t X ( p ) < 3 t X ( u ) , h child ( u ) ≤ min { h ( u ) , 3 t X ( u )} , t X ( u ) ≤ h ( u ) . (12) Along a vertex’s recursive history, let t i t_i t i be its parent distance to the protected path and σ i = min j ≤ i t j \sigma_i=\min_{j\le i}t_j σ i = min j ≤ i t j . Then σ i ≤ t i ≤ 3 σ i \sigma_i\le t_i\le3\sigma_i σ i ≤ t i ≤ 3 σ i .
A capture with R p ≤ R u / χ R_p \le R_u/\chi R p ≤ R u / χ must have an ineligible endpoint u u u . For the near set N p = { u ∈ X p : R u < χ R p } N_p = \{u \in X_p : R_u < \chi R_p\} N p = { u ∈ X p : R u < χ R p } , every two vertices are at physical distance less than 6 χ t X ( p ) 6\chi t_X(p) 6 χ t X ( p ) in the same residual Y Y Y . If σ X ( u ) \sigma_X(u) σ X ( u ) is the running minimum at this parent, then
N p ⊆ B Y , ℓ ( u , A 0 σ X ( u ) ) ( u ∈ N p ) . (13) N_p \subseteq B_{Y,\ell}(u,A_0\sigma_X(u)) \qquad(u \in N_p). \tag*{(13)} N p ⊆ B Y , ℓ ( u , A 0 σ X ( u )) ( u ∈ N p ) . (13)
Moreover, any physical shortest p p p –u u u path in Y Y Y of length less than R p / 8 R_p/8 R p /8 lies in every allowed cell, and bounds the child’s distance from u u u to its new protected path.
For an edge of length ℓ \ell ℓ , stop tracking it in an instance as soon as an endpoint has t X ≤ 16 χ ℓ / ε t_X \le16\chi\ell/\varepsilon t X ≤ 16 χ ℓ / ε . Its final tree distance is then at most 34 ( 32 χ / ε + 1 ) ℓ 34(32\chi/\varepsilon+1)\ell 34 ( 32 χ / ε + 1 ) ℓ . While both endpoints remain above this threshold, either endpoint u u u gives the bound d T ( u , v ) ≤ W t X ( u ) d_T(u,v) \le Wt_X(u) d T ( u , v ) ≤ W t X ( u ) with W = 3 ⋅ 34 = 102 W = 3\cdot34 = 102 W = 3 ⋅ 34 = 102 . If R u < χ R p R_u < \chi R_p R u < χ R p for one endpoint and the source can separate such an edge by an allowed radius, its lower entry threshold is greater than R p / 4 R_p/4 R p /4 .
Proof. The allowance satisfies r p ≤ ε t X ( p ) ≤ ε s p r_p \le\varepsilon t_X(p) \le\varepsilon s_p r p ≤ ε t X ( p ) ≤ ε s p . Since ε ≤ 1 / 8 \varepsilon\le1/8 ε ≤ 1/8 , the two locality bounds in Lemma 4.4 give d Y , ℓ ( p , u ) ≤ 3 t X ( u ) d_{Y,\ell}(p,u) \le3t_X(u) d Y , ℓ ( p , u ) ≤ 3 t X ( u ) and h c h i l d ( u ) ≤ 3 t X ( u ) h_{\mathrm{child}}(u) \le3t_X(u) h child ( u ) ≤ 3 t X ( u ) . The cutoff is positive, so (8) also gives h c h i l d ( u ) ≤ h ( u ) h_{\mathrm{child}}(u) \le h(u) h child ( u ) ≤ h ( u ) . The suffix of the S S S -root path to u u u after its last vertex on P P P is physical and has length at most h ( u ) h(u) h ( u ) , proving t X ( u ) ≤ h ( u ) t_X(u) \le h(u) t X ( u ) ≤ h ( u ) . Physical Lipschitzness and the more precise source-locality bound yield
t X ( p ) ≤ t X ( u ) + d X , ℓ ( p , u ) ≤ 4 t X ( u ) − t X ( p ) + 4 ε t X ( p ) . t_X(p) \le t_X(u) + d_{X,\ell}(p,u) \le4t_X(u) - t_X(p) + 4\varepsilon t_X(p). t X ( p ) ≤ t X ( u ) + d X , ℓ ( p , u ) ≤ 4 t X ( u ) − t X ( p ) + 4 ε t X ( p ) .
Consequently t X ( p ) ≤ 2 t X ( u ) / ( 1 − 2 ε ) ≤ 8 t X ( u ) / 3 < 3 t X ( u ) t_X(p) \le2t_X(u)/(1-2\varepsilon) \le8t_X(u)/3 < 3t_X(u) t X ( p ) ≤ 2 t X ( u ) / ( 1 − 2 ε ) ≤ 8 t X ( u ) /3 < 3 t X ( u ) . For an earlier index j < i j < i j < i , decreasing root heights give t i ≤ h i ≤ h j + 1 ≤ 3 t j t_i \le h_i \le h_{j+1} \le3t_j t i ≤ h i ≤ h j + 1 ≤ 3 t j . Taking the running minimum proves the claimed comparison with σ i \sigma_i σ i .
An eligible captured endpoint would satisfy R p / R u ≥ ( 1 − ε ) / ( 1 + ε ) ≥ 7 / 9 > 1 / χ R_p/R_u \ge(1-\varepsilon)/(1+\varepsilon) \ge7/9 > 1/\chi R p / R u ≥ ( 1 − ε ) / ( 1 + ε ) ≥ 7/9 > 1/ χ , excluding a far capture. For u ∈ N p u \in N_p u ∈ N p we have t X ( p ) / 3 < t X ( u ) < χ t X ( p ) t_X(p)/3 < t_X(u) < \chi t_X(p) t X ( p ) /3 < t X ( u ) < χ t X ( p ) , and hence σ X ( u ) > t X ( p ) / 9 \sigma_X(u) > t_X(p)/9 σ X ( u ) > t X ( p ) /9 . Joining two near vertices through the source gives residual distance less than 6 χ t X ( p ) 6\chi t_X(p) 6 χ t X ( p ) . Since A 0 ≥ 54 χ A_0 \ge54\chi A 0 ≥ 54 χ , their outer balls contain the full near set, proving (13) . For a shortest residual path of length less than R p / 8 R_p/8 R p /8 , every vertex on the path has gauge below 4 ( R p / 8 ) = R p / 2 4(R_p/8) = R_p/2 4 ( R p /8 ) = R p /2 and is therefore captured at every allowed radius. The source belongs to the new protected path, so this retained physical path bounds the child’s distance to that path.
The stopping and active bounds follow from (7) . At a stop use the endpoint with t X ≤ 16 χ ℓ / ε t_X \le16\chi\ell/\varepsilon t X ≤ 16 χ ℓ / ε ; during activity, ℓ < t X ( u ) \ell< t_X(u) ℓ < t X ( u ) and 34 ( 2 t X ( u ) + ℓ ) ≤ 102 t X ( u ) 34(2t_X(u)+\ell) \le102t_X(u) 34 ( 2 t X ( u ) + ℓ ) ≤ 102 t X ( u ) . If an active endpoint is near, then
R p > R u / χ = ε t X ( u ) / χ > 16 ℓ . R_p > R_u/\chi= \varepsilon t_X(u)/\chi> 16\ell. R p > R u / χ = ε t X ( u ) / χ > 16 ℓ .
The two entry gauges differ by at most 4 ℓ 4\ell 4 ℓ . If an allowed radius can separate the edge, the larger gauge is at least R p / 2 R_p/2 R p /2 , so the smaller is at least R p / 2 − 4 ℓ > R p / 4 R_p/2 - 4\ell> R_p/4 R p /2 − 4 ℓ > R p /4 , as claimed. □
Discrete choices and hypothetical captures. All conclusions hold without uniqueness of shortest paths, edge subdivision, or perturbation of lengths. A cutoff equal to a path-vertex height excludes that vertex and uses the first vertex of greater height as a p a_p a p ; any minimizing gauge witness gives the same inequalities. All child graphs are induced and retain every original edge whose endpoints survive. The source, the eligible set, the residual, and the parent data are fixed before the radius is chosen. The lemmas therefore also apply to any hypothetical allowed-radius capture at that state, independently of the numerical rule that defines eligibility.
Entropy bounds for the recursive construction
Fix positive rational constants ε , δ , χ , A 0 , ζ \varepsilon, \delta, \chi, A_0, \zeta ε , δ , χ , A 0 , ζ satisfying
0 < ε ≤ 1 8 , 0 < δ ≤ 1 256 , χ ≥ 4 , A 0 ≥ 54 χ , 0 < ζ ≤ 2 − 10 . (14) 0 < \varepsilon\le\frac{1}{8}, \qquad0 < \delta\le\frac{1}{256}, \qquad\chi\ge4, \qquad A_0 \ge54\chi, \qquad0 < \zeta\le2^{-10}. \tag*{(14)} 0 < ε ≤ 8 1 , 0 < δ ≤ 256 1 , χ ≥ 4 , A 0 ≥ 54 χ , 0 < ζ ≤ 2 − 10 . (14)
The parameters control the local allowance, core radius, far-source separation, outer-ball radius, and low-entropy cutoff, respectively. They are absolute constants, independent of the input. We give one concrete choice in § 11 . Throughout this section the physical lengths and the positive vertex measure μ \mu μ remain fixed. Write
M 0 = μ ( V ) , H 0 = ∑ v ∈ V μ ( v ) log M 0 μ ( v ) . M_0 = \mu(V), \qquad H_0 = \sum_{v\in V} \mu(v)\log\frac{M_0}{\mu(v)}. M 0 = μ ( V ) , H 0 = v ∈ V ∑ μ ( v ) log μ ( v ) M 0 .
All unmarked logarithms are natural.
Rates and charges
The source rule. In a parent X X X with protected path P X P_X P X , keep
t X ( v ) = d X , ℓ ( v , P X ) , R v = ε t X ( v ) , L v = h ( v ) − R v , M = μ ( X ) . t_X(v) = d_{X,\ell}(v, P_X), \qquad R_v = \varepsilon t_X(v), \qquad L_v = h(v) - R_v, \qquad M = \mu(X). t X ( v ) = d X , ℓ ( v , P X ) , R v = ε t X ( v ) , L v = h ( v ) − R v , M = μ ( X ) .
fixed throughout its carving. The candidates are the surviving vertices outside P X P_X P X . In residual Y Y Y , define
C v = B Y , ℓ ( v , δ R v ) , m v = μ ( C v ) , κ v = log 2 M m v , K v = ⌈ log 2 2 M m v ⌉ . C_v = B_{Y,\ell}(v,\delta R_v), \qquad m_v = \mu(C_v), \qquad\kappa_v = \log\frac{2M}{m_v}, \qquad K_v = \left\lceil\log_2 \frac{2M}{m_v} \right\rceil. C v = B Y , ℓ ( v , δ R v ) , m v = μ ( C v ) , κ v = log m v 2 M , K v = ⌈ log 2 m v 2 M ⌉ .
The integer K v K_v K v is computed by comparisons with powers of two. A candidate is eligible when 4 K v / R v ≤ Λ 4K_v/R_v \le\Lambda 4 K v / R v ≤ Λ ; among eligible candidates choose a source maximizing L v L_v L v , with a fixed tie rule. The root threshold is the minimum initial candidate rate. A threshold doubles only when no candidate is eligible, and a child inherits twice the threshold at its own source’s selection. The radius may be any r p ∈ ( R p / 2 , R p ) r_p \in(R_p/2,R_p) r p ∈ ( R p /2 , R p ) . Since κ v ≥ log 2 \kappa_v \ge\log2 κ v ≥ log 2 ,
κ v ≤ K v ≤ κ v log 2 + 1 ≤ 2 κ v log 2 < 3 κ v . (15) \kappa_v \le K_v \le\frac{\kappa_v}{\log2} + 1 \le\frac{2\kappa_v}{\log2} < 3\kappa_v. \tag*{(15)} κ v ≤ K v ≤ log 2 κ v + 1 ≤ log 2 2 κ v < 3 κ v . (15)
Thus source eligibility implies Λ R p ≥ 4 log ( 2 M / m p ) \Lambda R_p \ge4\log(2M/m_p) Λ R p ≥ 4 log ( 2 M / m p ) . Within one parent, residual restriction can only decrease m v m_v m v , so K v K_v K v is nondecreasing.
Geometric consequences. Lemmas 4.4 and 4.5 imply, for every capture of v v v by p p p ,
d Y , ℓ ( p , v ) ≤ 3 t X ( v ) , t X ( p ) < 3 t X ( v ) , h c h i l d ( v ) ≤ min { h ( v ) , 3 t X ( v ) } , t X ( v ) ≤ h ( v ) . (16) d_{Y,\ell}(p,v) \le3t_X(v), \qquad t_X(p) < 3t_X(v), \qquad h_{\mathrm{child}}(v) \le\min\{h(v),3t_X(v)\}, \qquad t_X(v) \le h(v). \tag*{(16)} d Y , ℓ ( p , v ) ≤ 3 t X ( v ) , t X ( p ) < 3 t X ( v ) , h child ( v ) ≤ min { h ( v ) , 3 t X ( v )} , t X ( v ) ≤ h ( v ) . (16)
Indeed ε ≤ 1 / 8 \varepsilon\le1/8 ε ≤ 1/8 gives t X ( p ) ≤ 2 t X ( v ) / ( 1 − 2 ε ) ≤ 8 t X ( v ) / 3 < 3 t X ( v ) t_X(p) \le2t_X(v)/(1-2\varepsilon) \le8t_X(v)/3 < 3t_X(v) t X ( p ) ≤ 2 t X ( v ) / ( 1 − 2 ε ) ≤ 8 t X ( v ) /3 < 3 t X ( v ) . If the captured vertex is eligible, then
h c h i l d ( v ) ≤ ε t X ( v ) , 1 − ε 1 + ε ≤ R p R v ≤ 1 + ε 1 − ε . (17) h_{\mathrm{child}}(v) \le\varepsilon t_X(v), \qquad\frac{1-\varepsilon}{1+\varepsilon} \le\frac{R_p}{R_v} \le\frac{1+\varepsilon}{1-\varepsilon}. \tag*{(17)} h child ( v ) ≤ ε t X ( v ) , 1 + ε 1 − ε ≤ R v R p ≤ 1 − ε 1 + ε . (17)
These conclusions also hold for a hypothetical allowed-radius capture at the unchanged pre-draw state. The entry gauge D p D_p D p is nonnegative, 4-Lipschitz in physical residual distance, and zero at p p p . Consequently C p C_p C p lies in every allowed cell, since 4 δ R p ≤ R p / 64 < R p / 2 4\delta R_p \le R_p/64 < R_p/2 4 δ R p ≤ R p /64 < R p /2 . If d Y , ℓ ( p , v ) < R p / 8 d_{Y,\ell}(p,v) < R_p/8 d Y , ℓ ( p , v ) < R p /8 , the entire physical shortest p p p –v v v path lies in every allowed cell. Since p p p is on the child’s protected path, this also gives t c h i l d ( v ) ≤ d Y , ℓ ( p , v ) t_{\mathrm{child}}(v) \le d_{Y,\ell}(p,v) t child ( v ) ≤ d Y , ℓ ( p , v ) .
Follow a vertex’s actual capturing cells until it reaches a protected path. At its i i i th capture, write t i , Λ i , Y i , M i t_i,\Lambda_i,Y_i,M_i t i , Λ i , Y i , M i for the parent scale, capture threshold, pre-capture residual, and parent mass. Put
σ i = min j ≤ i t j , x i = Λ i t i , O i = B Y i , ℓ ( v , A 0 σ i ) , I i = B Y i , ℓ ( v , δ ε σ i ) , B i = log M i μ ( I i ) . \sigma_i=\min_{j\le i}t_j,\qquad x_i=\Lambda_i t_i,\qquad O_i=B_{Y_i,\ell}(v,A_0\sigma_i),\qquad I_i=B_{Y_i,\ell}(v,\delta\varepsilon\sigma_i),\qquad B_i=\log\frac{M_i}{\mu(I_i)}. σ i = j ≤ i min t j , x i = Λ i t i , O i = B Y i , ℓ ( v , A 0 σ i ) , I i = B Y i , ℓ ( v , δ ε σ i ) , B i = log μ ( I i ) M i .
A capture is far if R p ≤ R v / χ R_p\le R_v/\chi R p ≤ R v / χ and flagged if D p ( v ) > R p / 4 D_p(v)>R_p/4 D p ( v ) > R p /4 . Define the total charges
Q F = ∑ v μ ( v ) ∑ i : f a r x i , Q M = ∑ v μ ( v ) ∑ i : f l a g g e d x i . Q_F=\sum_v\mu(v)\sum_{i:\,\mathrm{far}}x_i,\qquad Q_M=\sum_v\mu(v)\sum_{i:\,\mathrm{flagged}}x_i. Q F = v ∑ μ ( v ) i : far ∑ x i , Q M = v ∑ μ ( v ) i : flagged ∑ x i .
Theorem 5.1 (Main-charge bound). For every admissible sequence of radii,
Q F ≤ C F H 0 , Q M ≤ C M H 0 . Q_F\le C_FH_0,\qquad Q_M\le C_MH_0. Q F ≤ C F H 0 , Q M ≤ C M H 0 . where
L = 1 + ⌈ log 2 A 0 δ ε ⌉ , C F = 144 χ L ( χ − 3 ) ε , L=1+\left\lceil\log_2\frac{A_0}{\delta\varepsilon}\right\rceil,\qquad C_F=\frac{144\chi L}{(\chi-3)\varepsilon}, L = 1 + ⌈ log 2 δ ε A 0 ⌉ , C F = ( χ − 3 ) ε 144 χ L , C R = 1 + 6 1 − 2 ε , C S = 24 ε ( 1 + log 2 ζ ) , C_R=1+\frac{6}{1-2\varepsilon},\qquad C_S=\frac{24}{\varepsilon}\left(1+\frac{\log2}{\zeta}\right), C R = 1 + 1 − 2 ε 6 , C S = ε 24 ( 1 + ζ log 2 ) , C M = C R [ C F + 12 C S ( L + 1 + ε C F 4 χ ) ] . C_M=C_R\left[C_F+12C_S\left(L+1+\frac{\varepsilon C_F}{4\chi}\right)\right]. C M = C R [ C F + 12 C S ( L + 1 + 4 χ ε C F ) ] . These constants depend only on the parameters in (14) .
Far captures and scale telescoping
Lemma 5.2 (Running minima and occupied bands). Along every vertex history, σ i ≤ t i ≤ 3 σ i \sigma_i\le t_i\le3\sigma_i σ i ≤ t i ≤ 3 σ i , the residuals Y i Y_i Y i are nested, and Λ i + 1 ≥ 2 Λ i \Lambda_{i+1}\ge2\Lambda_i Λ i + 1 ≥ 2 Λ i . For any subset J J J of capture indices, retain its last index j r j_r j r in each occupied dyadic band of σ i \sigma_i σ i , ordered chronologically. Then
∑ i ∈ J x i ≤ 12 ∑ r x j r , ∑ r log μ ( O j r ) μ ( I j r ) ≤ L log M 0 μ ( v ) . \sum_{i\in J}x_i\le12\sum_r x_{j_r},\qquad\sum_r\log\frac{\mu(O_{j_r})}{\mu(I_{j_r})}\le L\log\frac{M_0}{\mu(v)}. i ∈ J ∑ x i ≤ 12 r ∑ x j r , r ∑ log μ ( I j r ) μ ( O j r ) ≤ L log μ ( v ) M 0 .
Proof. If j < i j<i j < i realizes the running minimum, monotonicity of the root heights and (16) give t i ≤ h i ≤ h j + 1 ≤ 3 t j t_i\le h_i\le h_{j+1}\le3t_j t i ≤ h i ≤ h j + 1 ≤ 3 t j ; the case j = i j=i j = i is immediate. The next parent is the actual cell capturing the vertex, so its next pre-capture residual is contained in Y i Y_i Y i . Its threshold starts at 2 Λ i 2\Lambda_i 2 Λ i and can only increase.
Within one dyadic band, the scales t i t_i t i differ by less than a factor six. Backwards from the last selected index, the rates decrease by factors at least two. Summing this geometric series proves the first inequality. Successive occupied bands have strictly decreasing integer indices, so σ j r + L < 2 1 − L σ j r \sigma_{j_r+L}<2^{1-L}\sigma_{j_r} σ j r + L < 2 1 − L σ j r . The definition of L L L gives A 0 σ j r + L ≤ δ ε σ j r A_0\sigma_{j_r+L}\le\delta\varepsilon\sigma_{j_r} A 0 σ j r + L ≤ δ ε σ j r . Physical distances cannot decrease under induced restriction; hence O j r + L ⊆ I j r O_{j_r+L}\subseteq I_{j_r} O j r + L ⊆ I j r . For each residue class modulo L L L , the logarithmic mass ratios telescope to at most log ( M 0 / μ ( v ) ) \log(M_0/\mu(v)) log ( M 0 / μ ( v )) , since every inner ball contains v v v . Summing these bounds proves the second inequality.
Lemma 5.3 (Far captures). For every vertex,
∑ i : f a r x i ≤ C F log M 0 μ ( v ) . \sum_{i:\,\mathrm{far}}x_i\le C_F\log\frac{M_0}{\mu(v)}. i : far ∑ x i ≤ C F log μ ( v ) M 0 . This estimate includes far captures that are not flagged.
Proof. An eligible captured vertex satisfies R p / R v ≥ ( 1 − ε ) / ( 1 + ε ) ≥ 7 / 9 > 1 / χ R_p/R_v \ge(1-\varepsilon)/(1+\varepsilon) \ge7/9 > 1/\chi R p / R v ≥ ( 1 − ε ) / ( 1 + ε ) ≥ 7/9 > 1/ χ . Thus every far captured vertex is ineligible. By (15) , endpoint ineligibility and source eligibility imply
κ v > Λ R v 12 , κ p ≤ K p ≤ Λ R p 4 ≤ Λ R v 4 χ . \kappa_v > \frac{\Lambda R_v}{12}, \qquad\kappa_p \le K_p \le\frac{\Lambda R_p}{4} \le\frac{\Lambda R_v}{4\chi}. κ v > 12 Λ R v , κ p ≤ K p ≤ 4 Λ R p ≤ 4 χ Λ R v .
The fixed parent numerator cancels. Writing g χ = ( χ − 3 ) / ( 12 χ ) > 0 g_\chi= (\chi-3)/(12\chi) > 0 g χ = ( χ − 3 ) / ( 12 χ ) > 0 , we obtain
log m p m v > g χ Λ R v = g χ ε x i . (18) \log\frac{m_p}{m_v} > g_\chi\Lambda R_v = g_\chi\varepsilon x_i. \tag*{(18)} log m v m p > g χ Λ R v = g χ ε x i . (18)
The source core lies in O i O_i O i : its center is within 3 t i 3t_i 3 t i of v v v , its radius is at most δ ε t i / χ \delta\varepsilon t_i/\chi δ ε t i / χ , and t i ≤ 3 σ i t_i \le3\sigma_i t i ≤ 3 σ i . The parameter bounds give 3 ( 3 + δ ε / χ ) < 10 < A 0 3(3+\delta\varepsilon/\chi) < 10 < A_0 3 ( 3 + δ ε / χ ) < 10 < A 0 . The endpoint core contains I i I_i I i . Therefore
x i < 1 ε g χ log μ ( O i ) μ ( I i ) . (19) x_i < \frac{1}{\varepsilon g_\chi}\log\frac{\mu(O_i)}{\mu(I_i)}. \tag*{(19)} x i < ε g χ 1 log μ ( I i ) μ ( O i ) . (19)
Apply Lemma 5.2 first to all far indices and then to their representatives. The resulting constant is 12 L / ( ε g χ ) = C F 12L/(\varepsilon g_\chi)=C_F 12 L / ( ε g χ ) = C F .
Near captures and the seed argument
Lemma 5.4 (Seed reduction). Call a capture a seed if it is the vertex’s first capture, if the vertex is ineligible, or if the threshold increased inside its parent before the capture. Every seed satisfies
x i ≤ 8 ε K v ≤ 24 ε ( log 2 + B i ) . x_i \le\frac{8}{\varepsilon}K_v \le\frac{24}{\varepsilon}(\log2+B_i). x i ≤ ε 8 K v ≤ ε 24 ( log 2 + B i ) . If B i ≤ ζ B_i \le\zeta B i ≤ ζ , the seed is unflagged and its next capture, if any, follows an internal threshold increase. Consequently
∑ i : f l a g g e d x i ≤ C R ∑ i : s e e d B i > ζ x i , x i ≤ C S B i at each seed in the latter sum. \sum_{i:\,\mathrm{flagged}} x_i \le C_R \sum_{\substack{i:\,\mathrm{seed}\\ B_i>\zeta}} x_i,\qquad x_i \le C_S B_i \quad\text{at each seed in the latter sum.} i : flagged ∑ x i ≤ C R i : seed B i > ζ ∑ x i , x i ≤ C S B i at each seed in the latter sum.
Proof. At an ineligible seed, Λ R v < 4 K v \Lambda R_v < 4K_v Λ R v < 4 K v . Otherwise consider the last local doubling before the capture. Immediately before it, the vertex was still present and ineligible. Monotonicity of K v K_v K v gives Λ R v < 8 K v \Lambda R_v < 8K_v Λ R v < 8 K v at capture. At the initial root, if there was no increase, the initial minimum rate is at most the vertex’s current rate. This proves the first bound. The endpoint core contains I i I_i I i , so κ v ≤ log 2 + B i \kappa_v \le\log2+B_i κ v ≤ log 2 + B i ; now apply (15) .
Suppose B i ≤ ζ B_i \le\zeta B i ≤ ζ . Since 2 M i / m v ≤ 2 e B i ≤ 2 e ζ < 4 2M_i/m_v \le2e^{B_i} \le2e^\zeta< 4 2 M i / m v ≤ 2 e B i ≤ 2 e ζ < 4 , one has K v ≤ 2 K_v \le2 K v ≤ 2 and Λ R v ≤ 16 \Lambda R_v \le16 Λ R v ≤ 16 . The strict source-scale bound R p < 3 R v R_p < 3R_v R p < 3 R v and source eligibility yield K p ≤ Λ R p / 4 < 12 K_p \le\Lambda R_p/4 < 12 K p ≤ Λ R p /4 < 12 . Integrality gives K p ≤ 11 K_p \le11 K p ≤ 11 , and consequently
m p ≥ 2 M i 2 11 = 2 − 10 M i . m_p \ge\frac{2M_i}{2^{11}} = 2^{-10}M_i. m p ≥ 2 11 2 M i = 2 − 10 M i .
The complement of I i I_i I i in its parent has mass M i ( 1 − e − B i ) < 2 − 10 M i M_i(1-e^{-B_i}) < 2^{-10}M_i M i ( 1 − e − B i ) < 2 − 10 M i : for B i > 0 B_i>0 B i > 0 use 1 − e − B i < B i ≤ ζ ≤ 2 − 10 1-e^{-B_i}<B_i\le\zeta\le2^{-10} 1 − e − B i < B i ≤ ζ ≤ 2 − 10 , and for B i = 0 B_i=0 B i = 0 the complement has mass zero. Thus the source core intersects I i I_i I i . Also Λ R p ≥ 4 K p ≥ 4 \Lambda R_p \ge4K_p \ge4 Λ R p ≥ 4 K p ≥ 4 gives R p ≥ R v / 4 R_p \ge R_v/4 R p ≥ R v /4 . The intersection supplies
d Y i , ℓ ( p , v ) ≤ δ ( R p + R v ) < 4 δ R v ≤ 16 δ R p . d_{Y_{i,\ell}}(p,v) \le\delta(R_p+R_v) < 4\delta R_v \le16\delta R_p. d Y i , ℓ ( p , v ) ≤ δ ( R p + R v ) < 4 δ R v ≤ 16 δ R p .
Every vertex on a physical shortest p p p –v v v path has entry gauge less than 64 δ R p ≤ R p / 4 64\delta R_p \le R_p/4 64 δ R p ≤ R p /4 . The seed is therefore unflagged. The entire path is inside every allowed cell, and the source is on the new protected path, so t i + 1 ≤ 4 δ ε t i t_{i+1} \le4\delta\varepsilon t_i t i + 1 ≤ 4 δ ε t i . Every actual capture satisfies the floor Λ t X ( v ) > 4 / ( 3 ε ) \Lambda t_X(v)>4/(3\varepsilon) Λ t X ( v ) > 4/ ( 3 ε ) , because R p < 3 ε t X ( v ) R_p<3\varepsilon t_X(v) R p < 3 ε t X ( v ) and Λ R p ≥ 4 \Lambda R_p\ge4 Λ R p ≥ 4 . If the next capture used only the inherited threshold, then
x i + 1 ≤ 8 δ ε x i = 8 δ Λ R v ≤ 128 δ ≤ 1 2 < 4 3 ε , x_{i+1}\le8\delta\varepsilon x_i=8\delta\Lambda R_v\le128\delta\le\frac{1}{2}<\frac{4}{3\varepsilon}, x i + 1 ≤ 8 δ ε x i = 8 δ Λ R v ≤ 128 δ ≤ 2 1 < 3 ε 4 ,
a contradiction. Thus the next capture, if present, follows an internal threshold increase and is another seed.
Every maximal nonseed run follows a seed. Its first charge is at most six times that seed charge, by threshold inheritance and the arbitrary scale reset. Every subsequent predecessor is eligible, so later charges contract by 2 ε < 1 2\varepsilon<1 2 ε < 1 . The entire run costs at most 6 / ( 1 − 2 ε ) 6/(1-2\varepsilon) 6/ ( 1 − 2 ε ) times its preceding seed. A low-entropy seed is unflagged and cannot precede a nonseed run. Assigning the remaining flagged charges to high-entropy seeds proves the factor C R C_R C R . At these seeds, log 2 + B i ≤ ( 1 + log 2 / ζ ) B i \log2+B_i\le(1+\log2/\zeta)B_i log 2 + B i ≤ ( 1 + log 2/ ζ ) B i , proving the bound with C S C_S C S .
Proof of Theorem 5.1. Multiply Lemma 5.3 by μ ( v ) \mu(v) μ ( v ) and sum to obtain the far bound. For flagged charges apply Lemma 5.4 . Far high-entropy seeds are already paid by Q F Q_F Q F . For each vertex, group its remaining near high-entropy seeds into occupied dyadic bands and retain the last such seed in each band. Lemma 5.2 reduces their total charge to twelve times the sum over representatives. At each representative write
B i = log μ ( O i ) μ ( I i ) + log M i μ ( O i ) . B_i=\log\frac{\mu(O_i)}{\mu(I_i)}+\log\frac{M_i}{\mu(O_i)}. B i = log μ ( I i ) μ ( O i ) + log μ ( O i ) M i .
The first terms, summed with vertex masses, contribute at most L H 0 LH_0 L H 0 .
For the second terms fix an actual source cell X p X_p X p in a parent of mass M M M and partition it as
N p = { v ∈ X p : R v < χ R p } , F p = X p ∖ N p . N_p=\{v\in X_p:R_v<\chi R_p\},\qquad F_p=X_p\setminus N_p. N p = { v ∈ X p : R v < χ R p } , F p = X p ∖ N p .
Write N = μ ( N p ) N=\mu(N_p) N = μ ( N p ) , F = μ ( F p ) F=\mu(F_p) F = μ ( F p ) , and V = N + F V=N+F V = N + F . Every near representative in this cell belongs to N p N_p N p . For u , v ∈ N p u,v\in N_p u , v ∈ N p , (16) gives
t X ( p ) / 3 < t X ( u ) , t X ( v ) < χ t X ( p ) , d Y , ℓ ( u , v ) < 6 χ t X ( p ) . t_X(p)/3<t_X(u),t_X(v)<\chi t_X(p),\qquad d_{Y,\ell}(u,v)<6\chi t_X(p). t X ( p ) /3 < t X ( u ) , t X ( v ) < χ t X ( p ) , d Y , ℓ ( u , v ) < 6 χ t X ( p ) .
Since σ ( u ) ≥ t X ( u ) / 3 > t X ( p ) / 9 \sigma(u)\ge t_X(u)/3>t_X(p)/9 σ ( u ) ≥ t X ( u ) /3 > t X ( p ) /9 and A 0 ≥ 54 χ A_0\ge54\chi A 0 ≥ 54 χ , its outer ball contains all of N p N_p N p . Thus, for any representative subset of N p N_p N p ,
∑ u r e p r e s e n t e d i n N p μ ( u ) log M μ ( O ( u ) ) ≤ N log ( M / N ) ≤ V log ( M / V ) + F . \sum_{\substack{u\ \mathrm{represented}\\ \mathrm{in}\ N_p}}\mu(u)\log\frac{M}{\mu(O(u))}
\le N\log(M/N)
\le V\log(M/V)+F. u represented in N p ∑ μ ( u ) log μ ( O ( u )) M ≤ N log ( M / N ) ≤ V log ( M / V ) + F .
The last inequality follows because the derivative of z log ( M / z ) z\log(M/z) z log ( M / z ) is at least − 1 -1 − 1 on ( 0 , M ] (0,M] ( 0 , M ] . The terms V log ( M / V ) V\log(M/V) V log ( M / V ) sum to at most H 0 H_0 H 0 over the recursive partition: include the retained protected path as a terminal part of each parent partition and telescope the complete partition entropies. For v ∈ F p v\in F_p v ∈ F p ,
Λ t X ( v ) ≥ χ Λ R p ε ≥ 4 χ ε . \Lambda t_X(v)\ge\frac{\chi\Lambda R_p}{\varepsilon}\ge\frac{4\chi}{\varepsilon}. Λ t X ( v ) ≥ ε χ Λ R p ≥ ε 4 χ .
Therefore ∑ c e l l s F ≤ ε Q F / ( 4 χ ) \sum_{\mathrm{cells}}F\le\varepsilon Q_F/(4\chi) ∑ cells F ≤ ε Q F / ( 4 χ ) . Combining the local mass ratios, the global mass cancellation, the band factor twelve, and the seed reduction gives the stated C M C_M C M .
Credits for missed far opportunities
Theorem 5.5 (Far-opportunity credits). Fix an original edge e = u v e=uv e = uv . Follow its residual carving until first capture in each parent. Continue to a child only if both endpoints enter it, and stop at separation, a protected-path endpoint, or an earlier analytical stopping rule. At a reached source step set
τ = min { D p ( u ) , D p ( v ) } , ω = max { D p ( u ) , D p ( v ) } , I = 1 { max ( τ , R p / 2 ) < min ( ω , R p ) } . \tau= \min\{D_p(u),D_p(v)\}, \qquad\omega= \max\{D_p(u),D_p(v)\}, \qquad I = \mathbf{1}_{\{\max(\tau,R_p/2)<\min(\omega,R_p)\}}. τ = min { D p ( u ) , D p ( v )} , ω = max { D p ( u ) , D p ( v )} , I = 1 { m a x ( τ , R p /2 ) < m i n ( ω , R p )} . Label the step by an endpoint z z z attaining τ \tau τ , with a fixed tie rule. It is far-labelled z z z if I = 1 I=1 I = 1 and R p ≤ R z / χ R_p \le R_z/\chi R p ≤ R z / χ . Let θ j \theta_j θ j be any nonnegative pre-choice weights satisfying
∑ j r e a c h e d i n o n e p a r e n t θ j ≤ 1 3 \sum_{j\ \mathrm{reached\ in\ one\ parent}} \theta_j \le\frac{1}{3} j reached in one parent ∑ θ j ≤ 3 1 for every realized history. Then, separately for each fixed endpoint label,
∑ j f a r − l a b e l l e d z Λ j θ j ≤ C T log M 0 μ ( z ) , C T = 48 χ L ( χ − 3 ) ε . \sum_{j\ \mathrm{far-labelled}\ z} \Lambda_j\theta_j \le C_T\log\frac{M_0}{\mu(z)}, \qquad C_T=\frac{48\chi L}{(\chi-3)\varepsilon}. j far − labelled z ∑ Λ j θ j ≤ C T log μ ( z ) M 0 , C T = ( χ − 3 ) ε 48 χ L .
Proof. Boundary overlap implies that some allowed radius captures z z z . If z z z were eligible, (17) for that hypothetical capture would contradict R p ≤ R z / χ R_p \le R_z/\chi R p ≤ R z / χ . Thus z z z is ineligible, and the local gap (18) applies even if the actual radius misses both endpoints. Hypothetical capture also gives d Y , ℓ ( p , z ) ≤ 3 t X ( z ) d_{Y,\ell}(p,z) \le3t_X(z) d Y , ℓ ( p , z ) ≤ 3 t X ( z ) , so the source core lies in the local outer ball and the inner ball lies in the endpoint core. Use the running minimum over the entire surviving vertex history, including parents omitted below. The local certificate (19) therefore holds at every such step.
In each parent with far-labelled steps, retain the last reached one, with threshold Λ ‾ X z \overline{\Lambda}^{z}_X Λ X z and its pre-choice residual. Monotonicity of the threshold and the packing hypothesis imply
t X ( z ) ∑ j f a r − l a b e l l e d z Λ j θ j ≤ 1 3 Λ ‾ X z t X ( z ) . t_X(z)\sum_{j\ \mathrm{far-labelled}\ z}\Lambda_j\theta_j \le\frac{1}{3}\overline{\Lambda}^{z}_X t_X(z). t X ( z ) j far − labelled z ∑ Λ j θ j ≤ 3 1 Λ X z t X ( z ) .
These retained residuals are nested: the actual first-capturing cell lies in the residual at the retained step, and each subsequent residual lies inside that child. Their rates grow by at least two, because a child starts at twice the actual first-capture threshold, which is at least the preceding retained threshold. Omitting intervening parents preserves both properties.
Apply the occupied-band argument to these retained parents. Within one band the scales differ by less than six, and the rates form a backwards geometric series, giving factor twelve. At the last representative of each band the local mass certificate applies. After L L L occupied bands the later outer ball lies inside the earlier inner ball. The logarithms thus telescope to at most L log ( M 0 / μ ( z ) ) L\log(M_0/\mu(z)) L log ( M 0 / μ ( z )) . The product of constants is ( 1 / 3 ) ⋅ 12 L / ( ε g χ ) = C T (1/3)\cdot12L/(\varepsilon g\chi)=C_T ( 1/3 ) ⋅ 12 L / ( ε g χ ) = C T .
:::
Corollary 5.6 (Rational packed credits). Set θ p = m p 2 / ( 3 M 2 ) \theta_p=m_p^2/(3M^2) θ p = m p 2 / ( 3 M 2 ) at a source step. These weights satisfy Theorem 5.5 for every admissible radius sequence. If b e ≥ 0 b_e\ge0 b e ≥ 0 and μ ( v ) ≥ ∑ e ∋ v b e \mu(v)\ge\sum_{e\ni v}b_e μ ( v ) ≥ ∑ e ∋ v b e , then
∑ e b e ∑ active far steps of e t X ( z e ) Λ p θ p ≤ C T H 0 . \sum_e b_e\sum_{\text{active far steps of }e} t_X(z_e)\Lambda_p\theta_p \le C_T H_0. e ∑ b e active far steps of e ∑ t X ( z e ) Λ p θ p ≤ C T H 0 .
Proof. Guaranteed cores lie in their actual removed cells, which are disjoint in a parent. Hence
∑ p r e a c h e d θ p = 1 3 ∑ p ( m p M ) 2 ≤ 1 3 ( ∑ p m p M ) 2 ≤ 1 3 . \sum_{p\ \mathrm{reached}}\theta_p=\frac{1}{3}\sum_p\left(\frac{m_p}{M}\right)^2\le\frac{1}{3}\left(\sum_p\frac{m_p}{M}\right)^2\le\frac{1}{3}. p reached ∑ θ p = 3 1 p ∑ ( M m p ) 2 ≤ 3 1 ( p ∑ M m p ) 2 ≤ 3 1 .
Apply Theorem 5.5 separately to each fixed endpoint label of every edge. After multiplying by b e b_e b e and summing,
∑ e = u v b e ( log M 0 μ ( u ) + log M 0 μ ( v ) ) ≤ ∑ v μ ( v ) log M 0 μ ( v ) = H 0 . \sum_{e=uv} b_e\left(\log\frac{M_0}{\mu(u)}+\log\frac{M_0}{\mu(v)}\right)\leq\sum_v\mu(v)\log\frac{M_0}{\mu(v)}=H_0. e = uv ∑ b e ( log μ ( u ) M 0 + log μ ( v ) M 0 ) ≤ v ∑ μ ( v ) log μ ( v ) M 0 = H 0 .
Omitting stopped or separated edges only removes nonnegative terms.
A deterministic radius rule
We now construct a weighted low-stretch tree. Fix the parameters from (14) , and let
C g e o = 34 , W = 3 C g e o . C_{\mathrm{geo}}=34,\qquad W=3C_{\mathrm{geo}}. C geo = 34 , W = 3 C geo .
The constant C g e o C_{\mathrm{geo}} C geo is the protected-path guarantee in Lemma 4.3 ; W W W bounds the cost of separating an edge whose endpoint scales exceed its length. An auxiliary exponential density is used only to establish a local averaging inequality. The algorithm chooses a rational radius satisfying that inequality.
Fixed measure and active edges
For nonnegative demands a e a_e a e on the original edges, put b e = a e ℓ e b_e=a_e\ell_e b e = a e ℓ e and B = ∑ e b e B=\sum_e b_e B = ∑ e b e . If B = 0 B=0 B = 0 , return any spanning tree. Otherwise use the fixed measure
μ ( v ) = ∑ e ∋ v b e + B n , M 0 = 3 B , H 0 = ∑ v μ ( v ) log M 0 μ ( v ) ≤ 3 B log n . (20) \mu(v)=\sum_{e\ni v}b_e+\frac{B}{n},\qquad M_0=3B,\qquad H_0=\sum_v\mu(v)\log\frac{M_0}{\mu(v)}\leq3B\log n. \tag*{(20)} μ ( v ) = e ∋ v ∑ b e + n B , M 0 = 3 B , H 0 = v ∑ μ ( v ) log μ ( v ) M 0 ≤ 3 B log n . (20)
Start with a singleton protected path at an arbitrary root, and apply the source rule of § 5 . Complete each parent carving before processing its children in source order. Each child stores twice its own source-time threshold, regardless of subsequent increases in its parent.
Maintain a list of active original demand edges. At entry to a parent X X X , stop the analysis of an edge e = u v e=uv e = uv when
min { t X ( u ) , t X ( v ) } ≤ 16 χ ε ℓ e . (21) \min\{t_X(u),t_X(v)\}\leq\frac{16\chi}{\varepsilon}\ell_e. \tag*{(21)} min { t X ( u ) , t X ( v )} ≤ ε 16 χ ℓ e . (21)
Protected-path geometry and the connected-subtree invariant bound its final tree distance by K 0 ℓ e K_0\ell_e K 0 ℓ e , where
K 0 = C g e o ( 32 χ ε + 1 ) . K_0=C_{\mathrm{geo}}\left(\frac{32\chi}{\varepsilon}+1\right). K 0 = C geo ( ε 32 χ + 1 ) .
Pay this cost once. An active edge separated by a selected cell is also closed. An edge whose endpoints both enter a cell passes to the child; an edge whose endpoints both remain stays in the residual list. Stopping or closing an edge changes only this list, not the graph, measure, or source rule. Complete the structural construction even when the active list is empty.
For an active edge at a source step, either endpoint z z z gives the pre-choice bound
d T ( u , v ) ≤ C g e o ( 2 t X ( z ) + ℓ e ) ≤ W t X ( z ) . (22) d_T(u,v)\leq C_{\mathrm{geo}}(2t_X(z)+\ell_e)\leq Wt_X(z). \tag*{(22)} d T ( u , v ) ≤ C geo ( 2 t X ( z ) + ℓ e ) ≤ W t X ( z ) . (22)
Here t X ( z ) > ℓ e t_X(z)>\ell_e t X ( z ) > ℓ e follows from the stopping rule and (14) . These estimates hold in the final global tree because each recursively constructed tree is a connected subtree of it.
Implementation of the weighted construction. Algorithm 1 specifies the order of all graph operations. The original graph, physical lengths, demands, baseline weights, and vertex measure are shared by its recursive calls. Only the temporary lengths, distance tables, residual vertex sets, and active-edge lists belong to an individual instance. An active-edge list is bookkeeping for the radius objective: deleting an edge from this list never deletes it from an induced graph.
Algorithm 1. WeightedTree and its recursive procedure
Input: A connected loopless graph G = ( V , E ) G=(V,E) G = ( V , E ) , rational lengths ℓ e > 0 \ell_e>0 ℓ e > 0 , and rational demands a e ≥ 0 a_e\ge0 a e ≥ 0 . Use the parameters of (14) , with the concrete choice in § 11 , and W = 102 W=102 W = 102 .
Output: An original-edge spanning tree satisfying Theorem 6.5 .
Procedure WeightedTree( G , ℓ , a ) (G,\ell,a) ( G , ℓ , a ) .
Compute b e = a e ℓ e b_e=a_e\ell_e b e = a e ℓ e and B = ∑ e b e B=\sum_e b_e B = ∑ e b e . If B = 0 B=0 B = 0 , return any spanning tree, including the empty tree when ∣ V ∣ = 1 |V|=1 ∣ V ∣ = 1 .
Fix μ ( v ) = ∑ e ∋ v b e + B / ∣ V ∣ \mu(v)=\sum_{e\ni v}b_e+B/|V| μ ( v ) = ∑ e ∋ v b e + B /∣ V ∣ for the entire call. Choose the first vertex x x x in the fixed order and return Build( G , x , ( x ) , ⊥ , E ) (G,x,(x),\perp,E) ( G , x , ( x ) , ⊥ , E ) . The symbol ⊥ \perp ⊥ requests initialization of the root threshold.
Procedure Build( X , x , P , Λ i n , A ) (X,x,P,\Lambda_{\mathrm{in}},\mathcal{A}) ( X , x , P , Λ in , A ) . At entry, X = G [ V ( X ) ] X=G[V(X)] X = G [ V ( X )] is connected and P P P is a physical shortest path in X X X starting at x x x . Both endpoints of every edge in A \mathcal{A} A belong to X X X .
If V ( X ) = V ( P ) V(X)=V(P) V ( X ) = V ( P ) , stop all edges in A \mathcal{A} A and return E ( P ) E(P) E ( P ) ; their endpoint scales are zero.
Set w e = ℓ e / 2 w_e=\ell_e/2 w e = ℓ e /2 on E ( P ) E(P) E ( P ) and w e = ℓ e w_e=\ell_e w e = ℓ e otherwise. Compute current root distances h = d X , w ( x , ⋅ ) h=d_{X,w}(x,\cdot) h = d X , w ( x , ⋅ ) and a current shortest-path tree S S S containing P P P : fix the predecessors on P P P and choose a tight predecessor for each other vertex, with fixed ties. Compute t ( v ) = d X , ℓ ( v , P ) t(v)=d_{X,\ell}(v,P) t ( v ) = d X , ℓ ( v , P ) by a physical multi-source shortest-path computation, and set R v = ε t ( v ) R_v=\varepsilon t(v) R v = εt ( v ) and M = μ ( V ( X ) ) M=\mu(V(X)) M = μ ( V ( X )) . Keep h , S , t , R , M h,S,t,R,M h , S , t , R , M fixed throughout this parent.
Remove from A \mathcal{A} A each e = u v e=uv e = uv with min { t ( u ) , t ( v ) } ≤ 16 χ ℓ e / ε \min\{t(u),t(v)\}\le16\chi\ell_e/\varepsilon min { t ( u ) , t ( v )} ≤ 16 χ ℓ e / ε . Initialize Y = X Y=X Y = X , Λ = Λ i n \Lambda=\Lambda_{\mathrm{in}} Λ = Λ in , an empty ordered child list, and T = E ( P ) T=E(P) T = E ( P ) .
While V ( Y ) ∖ V ( P ) ≠ ∅ V(Y)\setminus V(P)\ne\varnothing V ( Y ) ∖ V ( P ) = ∅ :
(a) Compute all-pairs distances in ( Y , ℓ ) (Y,\ell) ( Y , ℓ ) and ( Y , w ) (Y,w) ( Y , w ) . For each candidate v ∈ V ( Y ) ∖ V ( P ) v\in V(Y)\setminus V(P) v ∈ V ( Y ) ∖ V ( P ) compute
m v = ∑ u : d Y , ℓ ( v , u ) ≤ δ R v μ ( u ) , K v = min { k ∈ Z ≥ 0 : 2 k m v ≥ 2 M } , λ v = 4 K v / R v . m_v=\sum_{u:\,d_{Y,\ell}(v,u)\le\delta R_v}\mu(u),\qquad K_v=\min\{k\in\mathbb{Z}_{\ge0}:2^k m_v\ge2M\},\qquad\lambda_v=4K_v/R_v. m v = u : d Y , ℓ ( v , u ) ≤ δ R v ∑ μ ( u ) , K v = min { k ∈ Z ≥ 0 : 2 k m v ≥ 2 M } , λ v = 4 K v / R v .
(b) If Λ = ⊥ \Lambda=\perp Λ =⊥ , set Λ = min v λ v \Lambda=\min_v\lambda_v Λ = min v λ v . If no v v v satisfies λ v ≤ Λ \lambda_v\le\Lambda λ v ≤ Λ , multiply Λ \Lambda Λ by the smallest power of two that makes one candidate eligible. Choose an eligible p p p maximizing h ( p ) − R p h(p)-R_p h ( p ) − R p , using the fixed tie rule.
(c) Call ChooseRadius at this state to obtain ( r , Z ) (r,Z) ( r , Z ) , where Z = { v ∈ Y : D p ( v ) < r } Z=\{v\in Y:D_p(v)<r\} Z = { v ∈ Y : D p ( v ) < r } . Put c = h ( p ) − r c=h(p)-r c = h ( p ) − r . On P S ( x , p ) P_S(x,p) P S ( x , p ) let a p a_p a p be the first vertex with h ( a p ) > c h(a_p)>c h ( a p ) > c , and let b p b_p b p be its predecessor.
(d) Append the child record
( G [ Z ] , a p , P S ( a p , p ) , 2 Λ , { e = u v ∈ A : u , v ∈ Z } ) . \left(G[Z],a_p,P_S(a_p,p),2\Lambda,\{e=uv\in\mathcal{A}:u,v\in Z\}\right). ( G [ Z ] , a p , P S ( a p , p ) , 2Λ , { e = uv ∈ A : u , v ∈ Z } ) .
Add the original edge a p b p a_pb_p a p b p to T T T . Close every active edge with exactly one endpoint in Z Z Z ; remove the child’s active edges from the parent list. Equivalently, replace A \mathcal{A} A by { e = u v ∈ A : u , v ∈ V ( Y ) ∖ Z } \{e=uv\in\mathcal{A}:u,v\in V(Y)\setminus Z\} { e = uv ∈ A : u , v ∈ V ( Y ) ∖ Z } . Replace Y Y Y by the induced graph Y [ V ( Y ) ∖ Z ] Y[V(Y)\setminus Z] Y [ V ( Y ) ∖ Z ] .
Process the stored child records in source order, calling Build on each record and adding its returned edges to T T T . Return T T T .
The threshold stored for a child is twice the value at its own source selection; later increases in the parent do not change that record. The loop continues even when A \mathcal{A} A is empty. Thus every vertex is included and every structural edge required by Lemma 4.2 is retained. Repeated threshold doublings can be batched because the residual and personal rates do not change during those doublings.
All shortest-path computations use exact rational arithmetic; the all-pairs implementation is the conservative one analyzed in Theorem 7.1 .
The one-state boundary identity
Fix the state before choosing a radius: the residual Y Y Y , source p p p , allowance R = R p R = R_p R = R p , parent mass M M M , core mass m = m p m = m_p m = m p , and threshold Λ \Lambda Λ . Source eligibility gives Λ R ≥ 4 log ( 2 M / m ) \Lambda R \ge4\log(2M/m) Λ R ≥ 4 log ( 2 M / m ) . For analysis only, consider
f ( r ) = Λ e − Λ ( r − R / 2 ) 1 − e − Λ R / 2 , R / 2 ≤ r ≤ R , β = ( e Λ R / 2 − 1 ) − 1 . f(r) = \frac{\Lambda e^{-\Lambda(r-R/2)}}{1-e^{-\Lambda R/2}}, \qquad R/2 \le r \le R, \qquad\beta= \left(e^{\Lambda R/2}-1\right)^{-1}. f ( r ) = 1 − e − Λ R /2 Λ e − Λ ( r − R /2 ) , R /2 ≤ r ≤ R , β = ( e Λ R /2 − 1 ) − 1 .
For an active edge e = u v e=uv e = uv , put
τ e = min { D p ( u ) , D p ( v ) } , ω e = max { D p ( u ) , D p ( v ) } , a e ∗ = max { τ e , R / 2 } , \tau_e = \min\{D_p(u),D_p(v)\}, \qquad\omega_e = \max\{D_p(u),D_p(v)\}, \qquad a_e^* = \max\{\tau_e,R/2\}, τ e = min { D p ( u ) , D p ( v )} , ω e = max { D p ( u ) , D p ( v )} , a e ∗ = max { τ e , R /2 } ,
w e = ( min { ω e , R } − a e ∗ ) + , I e = 1 { w e > 0 } , κ e = 1 − e − Λ w e . w_e = \left(\min\{\omega_e,R\}-a_e^*\right)_+, \qquad I_e = \mathbf{1}_{\{w_e>0\}}, \qquad\kappa_e = 1-e^{-\Lambda w_e}. w e = ( min { ω e , R } − a e ∗ ) + , I e = 1 { w e > 0 } , κ e = 1 − e − Λ w e .
Choose an endpoint z e z_e z e attaining τ e \tau_e τ e by a fixed tie rule. This label is determined before the radius and is captured whenever the step captures either endpoint.
Lemma 6.1 (Boundary identity). Conditional on the fixed state,
P c u t , e = κ e ( P c a p , e + β ) , κ e ≤ 4 ℓ e Λ I e , P_{\mathrm{cut},e}=\kappa_e(P_{\mathrm{cap},e}+\beta), \qquad\kappa_e\le4\ell_e\Lambda I_e, P cut , e = κ e ( P cap , e + β ) , κ e ≤ 4 ℓ e Λ I e , where P c a p , e P_{\mathrm{cap},e} P cap , e is the probability that at least one endpoint enters the selected cell. If I e = 1 I_e=1 I e = 1 and R − a e ∗ ≥ 1 / Λ R-a_e^*\ge1/\Lambda R − a e ∗ ≥ 1/Λ , then
P c u t , e ≤ ( 1 − e − 1 ) − 1 κ e P c a p , e < 2 κ e P c a p , e . P_{\mathrm{cut},e}\le(1-e^{-1})^{-1}\kappa_eP_{\mathrm{cap},e}<2\kappa_eP_{\mathrm{cap},e}. P cut , e ≤ ( 1 − e − 1 ) − 1 κ e P cap , e < 2 κ e P cap , e .
Proof. The gauge is 4-Lipschitz, so w e ≤ 4 ℓ e w_e\le4\ell_e w e ≤ 4 ℓ e . For I e = 1 I_e=1 I e = 1 , direct integration gives
P c a p , e = β ( e Λ ( R − a e ∗ ) − 1 ) , P c u t , e = κ e β e Λ ( R − a e ∗ ) . P_{\mathrm{cap},e}=\beta\left(e^{\Lambda(R-a_e^*)}-1\right), \qquad P_{\mathrm{cut},e}=\kappa_e\beta e^{\Lambda(R-a_e^*)}. P cap , e = β ( e Λ ( R − a e ∗ ) − 1 ) , P cut , e = κ e β e Λ ( R − a e ∗ ) .
The formulas include τ e ≤ R / 2 \tau_e\le R/2 τ e ≤ R /2 , when the capture probability is one. If I e = 0 I_e=0 I e = 0 , the cut probability and κ e \kappa_e κ e are both zero. The first assertion now follows from 1 − e − x ≤ x 1-e^{-x}\le x 1 − e − x ≤ x . In the stated nonterminal case, β ≤ P c a p , e / ( e − 1 ) \beta\le P_{\mathrm{cap},e}/(e-1) β ≤ P cap , e / ( e − 1 ) , proving the last inequality.
A boundary-capable edge is far at this step if R ≤ R z e / χ R\le R_{z_e}/\chi R ≤ R z e / χ , and near otherwise. A near edge is terminal if R − 1 / Λ < τ e < R R-1/\Lambda<\tau_e<R R − 1/Λ < τ e < R ; the other near boundary-capable edges are nonterminal. Since Λ R ≥ 4 \Lambda R\ge4 Λ R ≥ 4 , a near nonterminal edge satisfies R − a e ∗ ≥ 1 / Λ R-a_e^*\ge1/\Lambda R − a e ∗ ≥ 1/Λ . Let F \mathcal{F} F be the set of far edges and M \mathcal{M} M the union of far and near nonterminal edges. Both sets are fixed before choosing the radius.
Terminal windows and split entropy
For 0 ≤ V ≤ M 0\le V\le M 0 ≤ V ≤ M , define
H M ( V ) = V log ( M / V ) + ( M − V ) log ( M / ( M − V ) ) , H_M(V)=V\log(M/V)+(M-V)\log(M/(M-V)), H M ( V ) = V log ( M / V ) + ( M − V ) log ( M / ( M − V )) ,
with zero contributions at the endpoints.
Lemma 6.2 (Terminal credit). Let V ( r ) V(r) V ( r ) be the selected cell’s mass and D = R / ε D=R/\varepsilon D = R / ε . For any pre-choice subset Z \mathcal{Z} Z of terminal active edges,
D E ∑ e ∈ Z a e 1 { e c u t } ≤ 2 e c 0 ε E H M ( V ( r ) ) , c 0 = ( 4 − e ) / 3 > 0. D\,\mathbb{E}\sum_{e\in\mathcal{Z}}a_e\mathbf{1}_{\{e\ \mathrm{cut}\}}\le\frac{2e}{c_0\varepsilon}\mathbb{E}H_M(V(r)), \qquad c_0=(4-e)/3>0. D E e ∈ Z ∑ a e 1 { e cut } ≤ c 0 ε 2 e E H M ( V ( r )) , c 0 = ( 4 − e ) /3 > 0. :::`
Proof. Write x = Λ R x=\Lambda R x = Λ R and z = m / M z=m/M z = m / M . Eligibility gives
β ≤ z 2 4 − z 2 ≤ z 2 3 ≤ 1 3 . \beta\le\frac{z^2}{4-z^2}\le\frac{z^2}{3}\le\frac{1}{3}. β ≤ 4 − z 2 z 2 ≤ 3 z 2 ≤ 3 1 .
A terminal edge’s first entry threshold lies above R / 2 R/2 R /2 and within 1 / Λ 1/\Lambda 1/Λ of R R R . The boundary identity therefore gives
P c u t , e ≤ e κ e β ≤ 4 e ℓ e Λ β . P_{\mathrm{cut},e}\le e\kappa_e\beta\le4e\ell_e\Lambda\beta. P cut , e ≤ e κ e β ≤ 4 e ℓ e Λ β .
Both endpoints remain outside with probability at least 1 − β ( e − 1 ) ≥ c 0 1-\beta(e-1)\ge c_0 1 − β ( e − 1 ) ≥ c 0 . The endpoint budgets hold pointwise:
2 ∑ e ∈ Z b e 1 { both endpoints outside } ≤ μ ( Y ∖ X p ( r ) ) . 2\sum_{e\in\mathcal{Z}}b_e\mathbf{1}_{\{\text{both endpoints outside}\}}\le\mu(Y\setminus X_p(r)). 2 e ∈ Z ∑ b e 1 { both endpoints outside } ≤ μ ( Y ∖ X p ( r )) .
Taking expectations yields ∑ e ∈ Z b e ≤ E μ ( Y ∖ X p ( r ) ) / ( 2 c 0 ) \sum_{e\in\mathcal{Z}}b_e\le\mathbb{E}\mu(Y\setminus X_p(r))/(2c_0) ∑ e ∈ Z b e ≤ E μ ( Y ∖ X p ( r )) / ( 2 c 0 ) . The left side in the lemma is consequently at most
2 e c 0 ε x β E μ ( Y ∖ X p ( r ) ) . \frac{2e}{c_0\varepsilon}x\beta\,\mathbb{E}\mu(Y\setminus X_p(r)). c 0 ε 2 e x β E μ ( Y ∖ X p ( r )) .
The function x / ( e x / 2 − 1 ) x/(e^{x/2}-1) x / ( e x /2 − 1 ) decreases for x > 0 x>0 x > 0 . Using x ≥ 4 log ( 2 / z ) x\ge4\log(2/z) x ≥ 4 log ( 2/ z ) and max 0 < z ≤ 1 z log ( 2 / z ) = 2 / e \max_{0<z\le1}z\log(2/z)=2/e max 0 < z ≤ 1 z log ( 2/ z ) = 2/ e gives
x β ≤ 4 z 2 log ( 2 / z ) 4 − z 2 ≤ 4 3 z 2 log ( 2 / z ) ≤ 8 3 e z < z . x\beta\le\frac{4z^2\log(2/z)}{4-z^2}\le\frac{4}{3}z^2\log(2/z)\le\frac{8}{3e}z<z. x β ≤ 4 − z 2 4 z 2 log ( 2/ z ) ≤ 3 4 z 2 log ( 2/ z ) ≤ 3 e 8 z < z .
Every allowed cell contains its core, so V ≥ m V\ge m V ≥ m . Moreover μ ( Y ∖ X p ) ≤ M − V \mu(Y\setminus X_p)\le M-V μ ( Y ∖ X p ) ≤ M − V and − log ( 1 − y ) ≥ y -\log(1-y)\ge y − log ( 1 − y ) ≥ y for 0 ≤ y < 1 0\le y<1 0 ≤ y < 1 . Hence
m M μ ( Y ∖ X p ) ≤ V M ( M − V ) ≤ ( M − V ) log M M − V ≤ H M ( V ) . \frac{m}{M}\mu(Y\setminus X_p)\le\frac{V}{M}(M-V)\le(M-V)\log\frac{M}{M-V}\le H_M(V). M m μ ( Y ∖ X p ) ≤ M V ( M − V ) ≤ ( M − V ) log M − V M ≤ H M ( V ) .
Substitution proves the claim.
Lemma 6.3 (Rational entropy credits). For rational 0 < V < M 0<V<M 0 < V < M , put
s = min { V , M − V } , k = ⌈ log 2 ( M / s ) ⌉ , U M ( V ) = s ( 1 + k ) . s=\min\{V,M-V\},\qquad k=\lceil\log_2(M/s)\rceil,\qquad U_M(V)=s(1+k). s = min { V , M − V } , k = ⌈ log 2 ( M / s )⌉ , U M ( V ) = s ( 1 + k ) . Then
H M ( V ) ≤ U M ( V ) ≤ 3 log 2 H M ( V ) < 6 H M ( V ) . H_M(V)\le U_M(V)\le\frac{3}{\log2}H_M(V)<6H_M(V). H M ( V ) ≤ U M ( V ) ≤ log 2 3 H M ( V ) < 6 H M ( V ) . Over all actual source cells,
∑ X p U μ ( X ) ( μ ( X p ) ) ≤ 12 H 0 . \sum_{X_p}U_{\mu(X)}(\mu(X_p))\le12H_0. X p ∑ U μ ( X ) ( μ ( X p )) ≤ 12 H 0 .
Proof. By symmetry take V = s ≤ M / 2 V=s\le M/2 V = s ≤ M /2 . The inequality − log ( 1 − y ) ≤ y / ( 1 − y ) -\log(1-y)\le y/(1-y) − log ( 1 − y ) ≤ y / ( 1 − y ) implies
s log ( M / s ) ≤ H M ( s ) ≤ s ( log ( M / s ) + 1 ) ≤ s ( k + 1 ) . s\log(M/s)\le H_M(s)\le s(\log(M/s)+1)\le s(k+1). s log ( M / s ) ≤ H M ( s ) ≤ s ( log ( M / s ) + 1 ) ≤ s ( k + 1 ) .
Since k ≤ log ( M / s ) / log 2 + 1 k\le\log(M/s)/\log2+1 k ≤ log ( M / s ) / log 2 + 1 and log ( M / s ) ≥ log 2 \log(M/s)\ge\log2 log ( M / s ) ≥ log 2 , one also has U M ( s ) ≤ ( 3 / log 2 ) s log ( M / s ) U_M(s)\le(3/\log2)s\log(M/s) U M ( s ) ≤ ( 3/ log 2 ) s log ( M / s ) .
For any partition into masses v i v_i v i , ∑ i H M ( v i ) ≤ 2 ∑ i v i log ( M / v i ) \sum_i H_M(v_i)\le2\sum_i v_i\log(M/v_i) ∑ i H M ( v i ) ≤ 2 ∑ i v i log ( M / v i ) . To see this, normalize p i = v i / M p_i=v_i/M p i = v i / M and let f ( p ) = − p log p + ( 1 − p ) log ( 1 − p ) f(p)=-p\log p+(1-p)\log(1-p) f ( p ) = − p log p + ( 1 − p ) log ( 1 − p ) . On [ 0 , 1 / 2 ] [0,1/2] [ 0 , 1/2 ] the function is nonnegative and concave, with f ( 0 ) = 0 f(0)=0 f ( 0 ) = 0 , so f ( a ) + f ( b ) ≥ f ( a + b ) f(a)+f(b)\ge f(a+b) f ( a ) + f ( b ) ≥ f ( a + b ) when a + b ≤ 1 / 2 a+b\le1/2 a + b ≤ 1/2 . If every p i ≤ 1 / 2 p_i\le1/2 p i ≤ 1/2 , then ∑ i f ( p i ) ≥ 0 \sum_i f(p_i)\ge0 ∑ i f ( p i ) ≥ 0 . Otherwise the unique p 1 > 1 / 2 p_1>1/2 p 1 > 1/2 has complement q ′ = ∑ i > 1 p i < 1 / 2 q'=\sum_{i>1}p_i<1/2 q ′ = ∑ i > 1 p i < 1/2 , and ∑ i f ( p i ) ≥ f ( p 1 ) + f ( q ′ ) = 0 \sum_i f(p_i)\ge f(p_1)+f(q')=0 ∑ i f ( p i ) ≥ f ( p 1 ) + f ( q ′ ) = 0 . This is the required binary-entropy inequality.
Include each retained protected path as a terminal part of its parent’s partition. The complete partition entropies telescope to the entropy of the final parts, which is at most H 0 H_0 H 0 . The omitted part credits are nonnegative; multiplication by the factor six proves the total budget.
Selecting a radius
For a radius in the open support interval away from every vertex entry threshold, compute
C ( r ) = W ∑ e a c t i v e a e t X ( z e ) 1 { τ e < r < ω e } , C(r)=W\sum_{e\ \mathrm{active}}a_e t_X(z_e)\mathbf{1}_{\{\tau_e<r<\omega_e\}}, C ( r ) = W e active ∑ a e t X ( z e ) 1 { τ e < r < ω e } ,
Q ( r ) = 8 W Λ ∑ e ∈ M b e t X ( z e ) I e 1 { τ e < r } , Q(r)=8W\Lambda\sum_{e\in\mathcal{M}}b_e t_X(z_e)I_e\mathbf{1}_{\{\tau_e<r\}}, Q ( r ) = 8 W Λ e ∈ M ∑ b e t X ( z e ) I e 1 { τ e < r } ,
F = 4 W Λ m p 2 3 M 2 ∑ e ∈ F b e t X ( z e ) I e , F=4W\Lambda\frac{m_p^2}{3M^2}\sum_{e\in\mathcal{F}}b_e t_X(z_e)I_e, F = 4 W Λ 3 M 2 m p 2 e ∈ F ∑ b e t X ( z e ) I e ,
U ( r ) = U M ( μ { v ∈ Y : D p ( v ) < r } ) . U(r)=U_M\bigl(\mu\{v\in Y:D_p(v)<r\}\bigr). U ( r ) = U M ( μ { v ∈ Y : D p ( v ) < r } ) .
The cell is nonempty and excludes the nonempty protected path, so its mass lies strictly between zero and M M M .
Lemma 6.4 (Finite radius rule). Among the open intervals between consecutive distinct values of
{ R / 2 , R } ∪ { D p ( v ) : v ∈ Y , R / 2 < D p ( v ) < R } , \{R/2,R\}\cup\{D_p(v):v\in Y,\ R/2<D_p(v)<R\}, { R /2 , R } ∪ { D p ( v ) : v ∈ Y , R /2 < D p ( v ) < R } , there is one whose midpoint satisfies
C ( r ) ≤ Q ( r ) + 18 χ W ε U ( r ) + F . (23) C(r)\le Q(r)+\frac{18\chi W}{\varepsilon}U(r)+F. \tag*{(23)} C ( r ) ≤ Q ( r ) + ε 18 χ W U ( r ) + F . (23) Selecting the first such midpoint in increasing order uses only exact rational arithmetic and comparisons with powers of two.
Proof. Use the auxiliary exponential density for this state only. For far edges retain both terms of Lemma 6.1 ; for near nonterminal edges use its factor-two absorption. After multiplication by the pre-choice bound (22) , all main terms are bounded by E Q ( r ) \mathbb{E}Q(r) E Q ( r ) . The far remainder is at most F F F , since β ≤ m p 2 / ( 3 M 2 ) \beta\le m_p^2/(3M^2) β ≤ m p 2 / ( 3 M 2 ) .
For near terminal edges, t X ( z e ) < χ R / ε t_X(z_e)<\chi R/\varepsilon t X ( z e ) < χ R / ε . Their entire expected cut cost is therefore at most
2 e χ W c 0 ε E H M ( V ( r ) ) < 18 χ W ε E U ( r ) , \frac{2e\chi W}{c_0\varepsilon}\mathbb{E}H_M(V(r))<\frac{18\chi W}{\varepsilon}\mathbb{E}U(r), c 0 ε 2 eχ W E H M ( V ( r )) < ε 18 χ W E U ( r ) ,
by Lemma 6.2 , using e < 3 e<3 e < 3 and c 0 > 1 / 3 c_0>1/3 c 0 > 1/3 . Consequently E C ( r ) ≤ E Q ( r ) + ( 18 χ W / ε ) E U ( r ) + F \mathbb{E}C(r)\le\mathbb{E}Q(r)+(18\chi W/\varepsilon)\mathbb{E}U(r)+F E C ( r ) ≤ E Q ( r ) + ( 18 χ W / ε ) E U ( r ) + F . Each function is constant on every listed open interval, and every nonempty interval has positive auxiliary probability. Thus at least one interval satisfies (23) .
The portal is also constant on each interval: on the source path, D p ( v ) = h ( p ) − h ( v ) D_p(v)=h(p)-h(v) D p ( v ) = h ( p ) − h ( v ) , so all portal transitions occur at listed vertex thresholds. The midpoint realizes the same cell and portal. Neither an auxiliary probability nor a conditional expectation of future cost is evaluated by the algorithm.
An exact sweep for the radius. The following helper implements (23) without sampling. It receives the parent data and residual distance tables already computed by BUILD. In particular, M = μ ( V ( X ) ) M=\mu(V(X)) M = μ ( V ( X )) is the fixed parent mass, not the residual mass. All events at the same threshold are processed together; tested radii are strictly between consecutive thresholds.
Algorithm 2. ChooseRadius
Input: The current parent and residual data ( X , Y , P , S , h , t , R , μ ) (X,Y,P,S,h,t,R,\mu) ( X , Y , P , S , h , t , R , μ ) , threshold Λ \Lambda Λ , eligible source p p p , its core mass m p m_p m p , current distance table d Y , w d_{Y,w} d Y , w , and active residual edges A \mathcal{A} A . The fixed weights a e , b e a_e,b_e a e , b e are shared.
Output: A rational radius r ∈ ( R p / 2 , R p ) r \in(R_p/2,R_p) r ∈ ( R p /2 , R p ) satisfying (23) , and its cell Z Z Z .
Put R = R p R=R_p R = R p and M = μ ( V ( X ) ) M=\mu(V(X)) M = μ ( V ( X )) , and let Q p Q_p Q p be the root-to-p p p path of S S S . For every v ∈ Y v \in Y v ∈ Y compute
D p ( v ) = h ( p ) − 2 h ( v ) + min q ∈ Q p { h ( q ) + 2 d Y , w ( q , v ) } . D_p(v)=h(p)-2h(v)+\min_{q\in Q_p}\{h(q)+2d_{Y,w}(q,v)\}. D p ( v ) = h ( p ) − 2 h ( v ) + q ∈ Q p min { h ( q ) + 2 d Y , w ( q , v )} .
For each e = u v ∈ A e=uv\in\mathcal{A} e = uv ∈ A , compute τ e = min { D p ( u ) , D p ( v ) } \tau_e=\min\{D_p(u),D_p(v)\} τ e = min { D p ( u ) , D p ( v )} , ω e = max { D p ( u ) , D p ( v ) } \omega_e=\max\{D_p(u),D_p(v)\} ω e = max { D p ( u ) , D p ( v )} , and a fixed endpoint z e z_e z e attaining τ e \tau_e τ e . Set
I e = 1 { max { τ e , R / 2 } < min { ω e , R } } . I_e=\mathbf{1}\{\max\{\tau_e,R/2\}<\min\{\omega_e,R\}\}. I e = 1 { max { τ e , R /2 } < min { ω e , R }} .
Form the fixed sets
F = { e ∈ A : I e = 1 , R ≤ R z e / χ } , \mathcal{F}=\{e\in\mathcal{A}:I_e=1,\ R\le R_{z_e}/\chi\}, F = { e ∈ A : I e = 1 , R ≤ R z e / χ } ,
M = F ∪ { e ∈ A : I e = 1 , R > R z e / χ , τ e ≤ R − 1 / Λ } . \mathcal{M}=\mathcal{F}\cup\{e\in\mathcal{A}:I_e=1,\ R>R_{z_e}/\chi,\ \tau_e\le R-1/\Lambda\}. M = F ∪ { e ∈ A : I e = 1 , R > R z e / χ , τ e ≤ R − 1/Λ } .
Precompute
c e = W a e t ( z e ) , q e = 8 W Λ b e t ( z e ) 1 { e ∈ M } , F = 4 W Λ m p 2 3 M 2 ∑ e ∈ F b e t ( z e ) . c_e=Wa_e t(z_e),\qquad q_e=8W\Lambda b_e t(z_e)\mathbf{1}_{\{e\in\mathcal{M}\}},\qquad F=\frac{4W\Lambda m_p^2}{3M^2}\sum_{e\in\mathcal{F}}b_e t(z_e). c e = W a e t ( z e ) , q e = 8 W Λ b e t ( z e ) 1 { e ∈ M } , F = 3 M 2 4 W Λ m p 2 e ∈ F ∑ b e t ( z e ) .
Sort the distinct values in { R / 2 , R } ∪ { D p ( v ) : R / 2 < D p ( v ) < R } \{R/2,R\}\cup\{D_p(v):R/2<D_p(v)<R\} { R /2 , R } ∪ { D p ( v ) : R /2 < D p ( v ) < R } as β 0 < ⋯ < β s \beta_0<\cdots<\beta_s β 0 < ⋯ < β s . Record each vertex’s threshold group. Using these indices, bucket the vertex mass increments and the active-edge entry and exit events at interior thresholds in one pass; the update sums below are the resulting bucket totals. For the first midpoint r 0 = ( β 0 + β 1 ) / 2 r_0=(\beta_0+\beta_1)/2 r 0 = ( β 0 + β 1 ) /2 , initialize
Z = { v : D p ( v ) < r 0 } , V = μ ( Z ) , C = ∑ e : τ e < r 0 < ω e c e , Q = ∑ e : τ e < r 0 q e . Z=\{v:D_p(v)<r_0\},\qquad V=\mu(Z),\qquad C=\sum_{e:\tau_e<r_0<\omega_e}c_e,\qquad Q=\sum_{e:\tau_e<r_0}q_e. Z = { v : D p ( v ) < r 0 } , V = μ ( Z ) , C = e : τ e < r 0 < ω e ∑ c e , Q = e : τ e < r 0 ∑ q e .
For i = 0 , … , s − 1 i=0,\ldots,s-1 i = 0 , … , s − 1 , in increasing order, set r = ( β i + β i + 1 ) / 2 r=(\beta_i+\beta_{i+1})/2 r = ( β i + β i + 1 ) /2 and compute
u = min { V , M − V } , k = min { j ∈ Z ≥ 0 : 2 j u ≥ M } , U = u ( 1 + k ) . u=\min\{V,M-V\},\qquad k=\min\{j\in\mathbb{Z}_{\ge0}:2^j u\ge M\},\qquad U=u(1+k). u = min { V , M − V } , k = min { j ∈ Z ≥ 0 : 2 j u ≥ M } , U = u ( 1 + k ) .
If C ≤ Q + ( 18 χ W / ε ) U + F C\le Q+(18\chi W/\varepsilon)U+F C ≤ Q + ( 18 χ W / ε ) U + F , return ( r , Z ) (r,Z) ( r , Z ) . Otherwise, if another interval remains, process all events at d = β i + 1 d=\beta_{i+1} d = β i + 1 together:
Z ← Z ∪ { v : D p ( v ) = d } , V ← V + ∑ v : D p ( v ) = d μ ( v ) , Z\leftarrow Z\cup\{v:D_p(v)=d\},\qquad V\leftarrow V+\sum_{v:D_p(v)=d}\mu(v), Z ← Z ∪ { v : D p ( v ) = d } , V ← V + v : D p ( v ) = d ∑ μ ( v ) ,
C ← C + ∑ e : τ e = d c e − ∑ e : ω e = d c e , Q ← Q + ∑ e : τ e = d q e . C\leftarrow C+\sum_{e:\tau_e=d}c_e-\sum_{e:\omega_e=d}c_e,\qquad Q\leftarrow Q+\sum_{e:\tau_e=d}q_e. C ← C + e : τ e = d ∑ c e − e : ω e = d ∑ c e , Q ← Q + e : τ e = d ∑ q e .
Lemma 6.4 guarantees that the helper returns. The cell always has 0 < V < M 0<V<M 0 < V < M , so the integer defining U U U exists. An edge with equal entry thresholds contributes cancelling events to C C C , and has q e = 0 q_e=0 q e = 0 . Each vertex enters once and each active edge has at most two events. Consequently, once the gauges have been computed, sorting and sweeping use O ( ∣ E ( Y ) ∣ + ∣ Y ∣ log ∣ Y ∣ ) O(|E(Y)|+|Y|\log|Y|) O ( ∣ E ( Y ) ∣ + ∣ Y ∣ log ∣ Y ∣ ) arithmetic operations. Only the actual entry thresholds are sorted; no rational grid is enumerated. The portal chosen by Build is constant on the returned interval.
Theorem 6.5 (Deterministic weighted spanning-tree bound). For positive physical lengths and
nonnegative edge demands, the source rule and finite radius rule return a spanning tree satisfying
C a ( T ) ≤ K 0 B + K 1 H 0 ≤ C d e t B log n , C_a(T) \le K_0B + K_1H_0 \le C_{\mathrm{det}}B\log n, C a ( T ) ≤ K 0 B + K 1 H 0 ≤ C det B log n ,
where
K 1 = 8 W ( C M + C F ) + 4 W C T + 216 χ W ε , C d e t = K 0 log 2 + 3 K 1 . K_1 = 8W(C_M + C_F) + 4WC_T + \frac{216\chi W}{\varepsilon}, \qquad C_{\mathrm{det}} = \frac{K_0}{\log2} + 3K_1. K 1 = 8 W ( C M + C F ) + 4 W C T + ε 216 χ W , C det = log 2 K 0 + 3 K 1 .
The constant depends only on the fixed admissible parameters. For rational input, the algorithm runs in deterministic polynomial time in the standard Turing model.
Proof. Every edge is paid at its first small-scale parent or its first separation. These events exhaust its finite history, because the final retained pieces are protected paths, on which endpoint scales are zero. The first type contributes at most K 0 B K_0B K 0 B , and the second contributes at most ∑ C ( r ) \sum C(r) ∑ C ( r ) .
At a selected radius, every main term in Q ( r ) Q(r) Q ( r ) is assigned to an endpoint z e z_e z e actually captured by that source. A far term is covered by Q F Q_F Q F . For a near term, the stopping rule implies
R p > R z e / χ > 16 ℓ e . R_p > R_{z_e}/\chi> 16\ell_e. R p > R z e / χ > 16 ℓ e .
Boundary overlap gives D p ( z e ) = τ e > R p / 2 − 4 ℓ e > R p / 4 D_p(z_e) = \tau_e > R_p/2 - 4\ell_e > R_p/4 D p ( z e ) = τ e > R p /2 − 4 ℓ e > R p /4 . Thus the capture is flagged and covered by Q M Q_M Q M . At any capture of a vertex, all assigned incident baseline weights sum to at most its fixed μ \mu μ -mass. Hence
∑ Q ( r ) ≤ 8 W ( Q M + Q F ) ≤ 8 W ( C M + C F ) H 0 . \sum Q(r) \le8W(Q_M + Q_F) \le8W(C_M + C_F)H_0. ∑ Q ( r ) ≤ 8 W ( Q M + Q F ) ≤ 8 W ( C M + C F ) H 0 .
Corollary 5.6 gives ∑ F ≤ 4 W C T H 0 \sum F \le4WC_TH_0 ∑ F ≤ 4 W C T H 0 , including missed far opportunities, and Lemma 6.3 gives ∑ U ( r ) ≤ 12 H 0 \sum U(r) \le12H_0 ∑ U ( r ) ≤ 12 H 0 . Summing (23) proves the first bound. All three credit budgets hold for every admissible trajectory. Finally H 0 ≤ 3 B log n H_0 \le3B\log n H 0 ≤ 3 B log n and n ≥ 2 n \ge2 n ≥ 2 imply the stated value of C d e t C_{\mathrm{det}} C det . § 7 proves the running-time bound.
Running time
Theorem 7.1 (Polynomial construction with exact dyadic eligibility). Let the positive physical lengths and nonnegative edge demands be rational, and let L L L be their total binary input length, including the graph encoding. Use fixed rational parameters satisfying (14) , the fixed measure μ ( v ) = ∑ e ∋ v a e ℓ e + B / n \mu(v) = \sum_{e\ni v} a_e\ell_e + B/n μ ( v ) = ∑ e ∋ v a e ℓ e + B / n , and the exact rules
K v ( Y ) = ⌈ log 2 2 μ ( X ) m v ( Y ) ⌉ , λ v ( Y ) = 4 K v ( Y ) R v , R v = ε d X , ℓ ( v , P X ) . K_v(Y) = \left\lceil\log_2 \frac{2\mu(X)}{m_v(Y)} \right\rceil, \qquad\lambda_v(Y) = \frac{4K_v(Y)}{R_v}, \qquad R_v = \varepsilon d_{X,\ell}(v,P_X). K v ( Y ) = ⌈ log 2 m v ( Y ) 2 μ ( X ) ⌉ , λ v ( Y ) = R v 4 K v ( Y ) , R v = ε d X , ℓ ( v , P X ) . Here B = ∑ e a e ℓ e B = \sum_e a_e\ell_e B = ∑ e a e ℓ e ; if B = 0 B = 0 B = 0 , any spanning tree suffices. By Lemma 6.4, a suitable interval midpoint can be found using only exact rational arithmetic and integer comparisons. The whole weighted-tree construction can be implemented in
O ( n 4 + n m ) O(n^4 + nm) O ( n 4 + nm ) arithmetic operations on numbers with O ( L + n + log ( 2 m ) ) O(L + n + \log(2m)) O ( L + n + log ( 2 m )) bits. A conservative running-time bound is therefore
O ( ( n 4 + n m ) [ L + n + log ( 2 m ) ] 3 ) . O\left((n^4 + nm)[L + n + \log(2m)]^3\right). O ( ( n 4 + nm ) [ L + n + log ( 2 m ) ] 3 ) . The bound is not optimized. Numerical aspect ratios enter only through their binary encodings.
Proof. Let D D D be a common denominator of the original physical lengths, obtained by multiplying their denominators. Its bit length is O ( L ) O(L) O ( L ) . Every current length is either an original physical length or half that length: the construction never repeatedly halves an inherited edge. Thus current shortest-path distances and heights are on the grid 1 / ( 2 D ) 1/(2D) 1/ ( 2 D ) , and physical distances are on the grid 1 / D 1/D 1/ D . Every shortest path is simple and has at most n − 1 n-1 n − 1 edges, so their numerators have O ( L + log n ) O(L+\log n) O ( L + log n ) bits.
Write the fixed constants as ε = a / b \varepsilon=a/b ε = a / b and δ = c / d \delta=c/d δ = c / d with positive integers a , b , c , d a,b,c,d a , b , c , d . Then R v R_v R v is on the grid 1 / ( b D ) 1/(bD) 1/ ( b D ) , a physical core radius is on the grid 1 / ( b d D ) 1/(bdD) 1/ ( b d D ) , and vertex-entry gauges are on the grid 1 / ( 2 D ) 1/(2D) 1/ ( 2 D ) . The endpoints R p / 2 R_p/2 R p /2 , R p R_p R p and all entry gauges in the allowed interval are therefore on the grid 1 / ( 2 b D ) 1/(2bD) 1/ ( 2 b D ) . Their interval midpoints can be taken on the grid 1 / ( 4 b D ) 1/(4bD) 1/ ( 4 b D ) . The algorithm enumerates only the at most n + 1 n+1 n + 1 intervals defined by the actual gauges; it never enumerates the grid.
The fixed baseline weights a ℓ e ℓ a_\ell e_\ell a ℓ e ℓ and all fixed vertex masses have a common denominator of O ( L + log n ) O(L+\log n) O ( L + log n ) bits. Every subset mass has the same encoding bound. Since
μ ( v ) ≥ B / n , μ ( X ) ≤ 3 B , \mu(v) \ge B/n,\qquad\mu(X) \le3B, μ ( v ) ≥ B / n , μ ( X ) ≤ 3 B ,
every candidate satisfies
2 ≤ 2 μ ( X ) m v ( Y ) ≤ 6 n , 1 ≤ K v ( Y ) ≤ ⌈ log 2 ( 6 n ) ⌉ . 2 \le\frac{2\mu(X)}{m_v(Y)} \le6n,\qquad1 \le K_v(Y) \le\lceil\log_2(6n)\rceil. 2 ≤ m v ( Y ) 2 μ ( X ) ≤ 6 n , 1 ≤ K v ( Y ) ≤ ⌈ log 2 ( 6 n )⌉ .
The ceiling of the base-two logarithm of a positive rational u / v u/v u / v is found exactly: use the difference of the bit lengths of the positive integers u , v u,v u , v to obtain one of the two possible integer values, and resolve the choice by comparing u u u with the appropriate power-of-two multiple of v v v . This uses integer shifts and one comparison, and evaluates no logarithm.
Let ℓ min , ℓ max \ell_{\min},\ell_{\max} ℓ m i n , ℓ m a x denote the smallest and largest physical lengths. For all personal rates, in every recursive instance,
L ∗ : = ε ( n − 1 ) ℓ max 4 ≤ λ v ( Y ) ≤ 4 ⌈ log 2 ( 6 n ) ⌉ ε ℓ min = : U ∗ . L_* := \frac{\varepsilon(n-1)\ell_{\max}}{4} \le\lambda_v(Y) \le\frac{4\lceil\log_2(6n)\rceil}{\varepsilon\ell_{\min}}=:U_*. L ∗ := 4 ε ( n − 1 ) ℓ m a x ≤ λ v ( Y ) ≤ ε ℓ m i n 4 ⌈ log 2 ( 6 n )⌉ =: U ∗ .
The initial threshold is one of the root’s rational personal rates. All later thresholds are dyadic multiples of it. An internal doubling can occur only below U ∗ U_* U ∗ , and its resulting value is less than 2 U ∗ 2U_* 2 U ∗ . There are at most n − 1 n-1 n − 1 successive inheritances on a branch, so every threshold is at most 2 n + 1 U ∗ 2^{n+1}U_* 2 n + 1 U ∗ . Its exponent relative to the initial threshold is
O ( n + log ( U ∗ / L ∗ ) ) = O ( n + L + log n ) . O\left(n+\log(U_*/L_*)\right)=O(n+L+\log n). O ( n + log ( U ∗ / L ∗ ) ) = O ( n + L + log n ) .
Thus thresholds have O ( L + n + log n ) O(L+n+\log n) O ( L + n + log n ) -bit encodings. When a residual has no eligible source, all consecutive doublings can be performed at once: determine the smallest dyadic multiple of the current threshold at least the minimum personal rate. The same bit-length comparison computes this integer exponent exactly. This produces precisely the threshold that consecutive doublings would have produced.
Every selected source becomes part of its child’s protected path. That path is excluded from all subsequent source cells. Hence a vertex is selected as a source at most once; there are at most n − 1 n-1 n − 1 source selections and at most n n n recursive instances.
At each removal, recompute all-pairs physical distances and, if desired, all-pairs current distances by Floyd–Warshall. This uses O ( n 3 + m ) O(n^3+m) O ( n 3 + m ) arithmetic operations per residual and provides all core masses in another O ( n 2 ) O(n^2) O ( n 2 ) operations. All exact dyadic rates and the next source can then be determined in O ( n ) O(n) O ( n ) operations. One can instead compute the chosen source’s gauge by a multi-source shortest-path computation; the conservative all-pairs bound already suffices.
Sort the at most n n n entry gauges and sweep the open intervals. Maintain the cell mass, the cut-cost sum, and the capture-credit sum; each vertex enters once and each active edge has at most two entry events. The far credit is constant during this sweep. The rational binary-entropy majorant is computed from the current cell mass by the same exact power-of-two comparison. The sweep costs O ( m + n log n ) O(m+n\log n) O ( m + n log n ) arithmetic operations. Source-path portal changes are among these entry thresholds, so a selected midpoint determines both the cell and its portal.
All local cost and credit expressions use subset masses, the fixed baseline weights, physical distances, and the one current threshold. Their encoding lengths are O ( L + n + log ( 2 m ) ) O(L+n+\log(2m)) O ( L + n + log ( 2 m )) . For example, the far credit’s extra denominator is the square of the current parent mass; this contributes only O ( L + log n ) O(L+\log n) O ( L + log n ) bits. Parent-dependent credits are not accumulated by the algorithm, so their denominators do not multiply along a recursive chain. Summing the operation counts over at most n − 1 n-1 n − 1 removals proves the displayed arithmetic bound. Exact rational arithmetic with a conservative cubic time per arithmetic operation gives the stated running-time bound.
Composition with the rational master. Let L 0 L_0 L 0 be the encoding length of the fixed metric or capacity instance, and fix rational 0 < ϑ ≤ 1 0<\vartheta\le1 0 < ϑ ≤ 1 . The master in Theorem 8.1 , applied in dimension m m m , uses
N m a s t e r = O ( m log ( 2 m ) ϑ 2 ) N_{\mathrm{master}}=O\left(\frac{m\log(2m)}{\vartheta^2}\right) N master = O ( ϑ 2 m log ( 2 m ) )
oracle calls. Its per-coordinate query encoding is uniformly
Q = O ( log ( 2 m ) ϑ 2 [ 1 + bits ( ϑ ) ] + log ( 2 m ) ) . Q=O\left(\frac{\log(2m)}{\vartheta^2}[1+\operatorname{bits}(\vartheta)]+\log(2m)\right). Q = O ( ϑ 2 log ( 2 m ) [ 1 + bits ( ϑ )] + log ( 2 m ) ) .
For distances the oracle demands are p e / ℓ e p_e/\ell_e p e / ℓ e ; for capacities the auxiliary lengths are p e / c e p_e/c_e p e / c e and the demands are c e c_e c e . In either case the weighted-oracle input has total encoding length O ( L 0 + m Q ) O(L_0+mQ) O ( L 0 + m Q ) . In the capacity case every auxiliary length is strictly positive because every master query coordinate is strictly positive.
Consequently the total time spent in weighted-tree oracle calls is at most
O ( N m a s t e r ( n 4 + n m ) [ L 0 + m Q + n + log ( 2 m ) ] 3 ) . O\left(N_{\mathrm{master}}(n^4+nm)[L_0+mQ+n+\log(2m)]^3\right). O ( N master ( n 4 + nm ) [ L 0 + m Q + n + log ( 2 m ) ] 3 ) .
Add the polynomial master cost proved in Theorem 8.1 , and the optional exact support-reduction cost. Each returned tree’s fixed stretch or congestion vector can be computed in O ( m log n + n ) O(m\log n+n) O ( m log n + n ) arithmetic operations using tree depths, lowest common ancestors, and subtree sums; its entries have O ( L 0 + log ( 2 n m ) ) O(L_0+\log(2nm)) O ( L 0 + log ( 2 nm )) bits, independently of the changing oracle query. This yields an explicit deterministic polynomial-time construction for fixed ϑ \vartheta ϑ , and a polynomial dependence on 1 / ϑ 1/\vartheta 1/ ϑ when accuracy is part of the input. The displayed bound establishes polynomial running time; improving it to near-linear time would require additional algorithmic work.
Constructing an explicit distribution
The response-oracle approach of [10 Răc08, 3 , 21 , 4 ] converts a weighted tree oracle to an explicit distribution. This conversion can preserve its approximation factor up to any prescribed multiplicative loss. The number of oracle calls need not depend on the numerical size of the response bound, or on the aspect ratio of the graph data. The normalization of the stopping horizon is essential to this statement.
Theorem 8.1 (Rational adaptive-step master). Let d ≥ 1 d \ge1 d ≥ 1 , let A ⊆ Q ≥ 0 d \mathcal{A} \subseteq\mathbb{Q}_{\ge0}^{d} A ⊆ Q ≥ 0 d be a finite family of nonnegative rational vectors, and let ρ > 0 \rho> 0 ρ > 0 be rational. Suppose that for every strictly positive rational probability vector p p p an oracle returns a ∈ A a \in\mathcal{A} a ∈ A satisfying
p ⋅ a ≤ ρ . p \cdot a \le\rho. p ⋅ a ≤ ρ . For every rational 0 < ϑ ≤ 1 0 < \vartheta\le1 0 < ϑ ≤ 1 , a deterministic algorithm uses at most
N ≤ 1 + ⌈ 20 d ⌈ log 2 ( 2 d ) ⌉ ϑ 2 ⌉ N \le1 + \left\lceil\frac{20d\lceil\log_{2}(2d)\rceil}{\vartheta^{2}} \right\rceil N ≤ 1 + ⌈ ϑ 2 20 d ⌈ log 2 ( 2 d )⌉ ⌉ oracle calls and returns rational probabilities λ t \lambda_t λ t and oracle outputs a ( t ) a^{(t)} a ( t ) such that
∑ t λ t = 1 , ∑ t λ t a j ( t ) ≤ ( 1 + ϑ ) ρ ( 1 ≤ j ≤ d ) . \sum_t \lambda_t = 1,\qquad\sum_t \lambda_t a_j^{(t)} \le(1+\vartheta)\rho\quad(1 \le j \le d). t ∑ λ t = 1 , t ∑ λ t a j ( t ) ≤ ( 1 + ϑ ) ρ ( 1 ≤ j ≤ d ) . No exponential or logarithmic function is evaluated by the algorithm. Its oracle queries have a uniform polynomial encoding bound in d d d , 1 / ϑ 1/\vartheta 1/ ϑ , and the encoding of ϑ \vartheta ϑ , independent of the numerical magnitude and encoding length of ρ \rho ρ .
If the oracle runs in polynomial time in the fixed input and the query encoding, and its returned vectors have polynomial encoding length in these parameters, then the entire construction has polynomial running time in the fixed input, 1 / ϑ 1/\vartheta 1/ ϑ , and the encodings of ρ \rho ρ and ϑ \vartheta ϑ . In particular, numerical aspect ratios impose no pseudopolynomial dependence.
Proof. Put
k = ⌈ log 2 ( 2 d ) ⌉ , ξ = ϑ / 4 , b = 1 + ξ , H = 10 k ϑ 2 ρ . k = \lceil\log_{2}(2d)\rceil,\qquad\xi= \vartheta/4,\qquad b = 1+\xi,\qquad H = \frac{10k}{\vartheta^{2}\rho}. k = ⌈ log 2 ( 2 d )⌉ , ξ = ϑ /4 , b = 1 + ξ , H = ϑ 2 ρ 10 k .
The integer k k k is computed by comparisons with powers of two. All other quantities in this display are rational. Initialize x = 0 ∈ Q d x = 0 \in\mathbb{Q}^{d} x = 0 ∈ Q d and τ = 0 \tau= 0 τ = 0 . While τ < H \tau< H τ < H , form
w j = b ⌈ x j ⌉ , p j = w j ∑ i w i w_j = b^{\lceil x_j\rceil},\qquad p_j = \frac{w_j}{\sum_i w_i} w j = b ⌈ x j ⌉ , p j = ∑ i w i w j
and obtain a response a a a from the oracle. Write L = max j a j L = \max_j a_j L = max j a j and choose
Δ = { H − τ , L = 0 , min { L − 1 , H − τ } , L > 0. \Delta=
\begin{cases}
H-\tau, & L=0,\\
\min\{L^{-1},H-\tau\}, & L>0.
\end{cases} Δ = { H − τ , min { L − 1 , H − τ } , L = 0 , L > 0.
Record the pair ( Δ / H , a ) (\Delta/H,a) ( Δ/ H , a ) and update x ← x + Δ a x \leftarrow x+\Delta a x ← x + Δ a and τ ← τ + Δ \tau\leftarrow\tau+\Delta τ ← τ + Δ . Every coordinate increment lies in [ 0 , 1 ] [0,1] [ 0 , 1 ] .
For analysis only, define Φ ( x ) = ∑ j b x j \Phi(x)=\sum_j b^{x_j} Φ ( x ) = ∑ j b x j . Convexity of b z b^z b z on [ 0 , 1 ] [0,1] [ 0 , 1 ] gives b z ≤ 1 + ( b − 1 ) z = 1 + ξ z b^z \le1+(b-1)z=1+\xi z b z ≤ 1 + ( b − 1 ) z = 1 + ξ z . Moreover, b x j ≤ w j ≤ b b x j b^{x_j}\le w_j\le b b^{x_j} b x j ≤ w j ≤ b b x j , so
∑ j b x j a j ≤ ∑ j w j a j = ( ∑ j w j ) ( p ⋅ a ) ≤ b ρ Φ ( x ) . \sum_j b^{x_j}a_j \le\sum_j w_j a_j = \left(\sum_j w_j\right)(p\cdot a) \le b\rho\Phi(x). j ∑ b x j a j ≤ j ∑ w j a j = ( j ∑ w j ) ( p ⋅ a ) ≤ b ρ Φ ( x ) .
It follows that
Φ ( x + Δ a ) ≤ Φ ( x ) + ξ Δ ∑ j b x j a j ≤ Φ ( x ) ( 1 + ξ b ρ Δ ) ≤ Φ ( x ) exp ( ξ b ρ Δ ) . \Phi(x+\Delta a) \le\Phi(x)+\xi\Delta\sum_j b^{x_j}a_j \le\Phi(x)(1+\xi b\rho\Delta) \le\Phi(x)\exp(\xi b\rho\Delta). Φ ( x + Δ a ) ≤ Φ ( x ) + ξ Δ j ∑ b x j a j ≤ Φ ( x ) ( 1 + ξ b ρ Δ ) ≤ Φ ( x ) exp ( ξ b ρ Δ ) .
Induction therefore gives, at every reached state,
Φ ( x ) ≤ d exp ( ξ b ρ τ ) , x j ≤ ln d ln b + ξ b ρ τ ln b . \Phi(x) \le d\exp(\xi b\rho\tau), \qquad x_j \le\frac{\ln d}{\ln b}+\frac{\xi b\rho\tau}{\ln b}. Φ ( x ) ≤ d exp ( ξ b ρ τ ) , x j ≤ ln b ln d + ln b ξ b ρ τ .
Here are explicit elementary estimates for the two coefficients. Since 0 < ξ ≤ 1 / 4 0<\xi\le1/4 0 < ξ ≤ 1/4 ,
ln ( 1 + ξ ) ≥ ξ 1 + ξ ≥ 4 ξ 5 = ϑ 5 . \ln(1+\xi)\ge\frac{\xi}{1+\xi}\ge\frac{4\xi}{5}=\frac{\vartheta}{5}. ln ( 1 + ξ ) ≥ 1 + ξ ξ ≥ 5 4 ξ = 5 ϑ .
Also ln ( 1 + ξ ) ≥ ξ − ξ 2 / 2 \ln(1+\xi)\ge\xi-\xi^2/2 ln ( 1 + ξ ) ≥ ξ − ξ 2 /2 , and hence
ξ b ln b ≤ 1 + ξ 1 − ξ / 2 ≤ 1 + 2 ξ = 1 + ϑ 2 . \frac{\xi b}{\ln b}\le\frac{1+\xi}{1-\xi/2}\le1+2\xi=1+\frac{\vartheta}{2}. ln b ξ b ≤ 1 − ξ /2 1 + ξ ≤ 1 + 2 ξ = 1 + 2 ϑ .
The penultimate inequality follows by multiplying its positive denominator; the difference of the resulting right and left sides is ξ ( 1 / 2 − ξ ) ≥ 0 \xi(1/2-\xi)\ge0 ξ ( 1/2 − ξ ) ≥ 0 . Since ln d ≤ k \ln d\le k ln d ≤ k and τ ≤ H \tau\le H τ ≤ H , all coordinates satisfy
x j ≤ 5 k ϑ + ( 1 + ϑ 2 ) 10 k ϑ 2 = 10 k ϑ 2 + 10 k ϑ ≤ 20 k ϑ 2 = : Γ . x_j\le\frac{5k}{\vartheta}+\left(1+\frac{\vartheta}{2}\right)\frac{10k}{\vartheta^2}
=\frac{10k}{\vartheta^2}+\frac{10k}{\vartheta}\le\frac{20k}{\vartheta^2}=:\Gamma. x j ≤ ϑ 5 k + ( 1 + 2 ϑ ) ϑ 2 10 k = ϑ 2 10 k + ϑ 10 k ≤ ϑ 2 20 k =: Γ.
Every nonfinal iteration has L > 0 L>0 L > 0 and Δ = 1 / L \Delta=1/L Δ = 1/ L , and raises at least one coordinate of x x x by exactly one. Thus it raises ∑ j x j \sum_j x_j ∑ j x j by at least one. This sum never exceeds d Γ d\Gamma d Γ ; there are at most ⌊ d Γ ⌋ \lfloor d\Gamma\rfloor ⌊ d Γ ⌋ nonfinal iterations. The algorithm consequently terminates within the asserted call bound. At termination the recorded coefficients sum to one, and their mean vector is x / H x/H x / H . Finally,
x j H ≤ ln d H ln b + ξ b ln b ρ ≤ ϑ 2 + ( 1 + ϑ 2 ) ρ = ( 1 + ϑ ) ρ . \frac{x_j}{H}\le\frac{\ln d}{H\ln b}+\frac{\xi b}{\ln b}\rho\le\frac{\vartheta}{2}+\left(1+\frac{\vartheta}{2}\right)\rho=(1+\vartheta)\rho. H x j ≤ H ln b ln d + ln b ξ b ρ ≤ 2 ϑ + ( 1 + 2 ϑ ) ρ = ( 1 + ϑ ) ρ .
The remaining assertions follow from the running-time analysis below.
Uniform query encodings. Write b = P / Q b=P/Q b = P / Q in lowest terms and let B ϑ B_\vartheta B ϑ bound the bit lengths of P P P and Q Q Q . At every iteration set e j = ⌈ x j ⌉ e_j=\lceil x_j\rceil e j = ⌈ x j ⌉ and E = max j e j ≤ ⌈ Γ ⌉ E=\max_j e_j\le\lceil\Gamma\rceil E = max j e j ≤ ⌈ Γ ⌉ . The query can be computed using the positive integers
W j = P e j Q E − e j , p j = W j ∑ i W i . W_j=P^{e_j}Q^{E-e_j}, \qquad p_j=\frac{W_j}{\sum_i W_i}. W j = P e j Q E − e j , p j = ∑ i W i W j .
Each coordinate of p p p has bit length
O ( Γ B ϑ + log ( 2 d ) ) . O(\Gamma B_\vartheta+\log(2d)). O ( Γ B ϑ + log ( 2 d )) .
This bound contains neither a previous query size nor a previous oracle output size. In particular, the polynomial-time oracle hypothesis is not being iterated through a potentially growing sequence of polynomial bounds. Powers of P P P and Q Q Q through ⌈ Γ ⌉ \lceil\Gamma\rceil ⌈ Γ ⌉ can be precomputed once, after which each query uses O ( d ) O(d) O ( d ) integer arithmetic operations.
Accumulated rational encodings and an explicit master bound. Let B B B be a uniform bound on the numerator and denominator lengths of the returned vector entries, and let B H B_H B H be the encoding length of H H H . Every nonfinal step size is the reciprocal of one returned entry and has O ( B ) O(B) O ( B ) bits. The final step is H H H minus the preceding sum of step sizes. A common-denominator representation consequently bounds the encoding lengths of all accumulated times, load coordinates, and output probabilities by
O ( N B + B H + log ( 2 N ) ) . O(NB+B_H+\log(2N)). O ( N B + B H + log ( 2 N )) .
For example, before the final iteration each load coordinate is a sum of at most N N N products of two B B B -bit rationals. The last iteration adds one product involving the final step size; its denominator adds at most the same total bound. The normalization by H H H adds O ( B H ) O(B_H) O ( B H ) bits. The magnitudes of the accumulated load coordinates are bounded by Γ \Gamma Γ as proved above.
Thus, putting
S = O ( N B + B H + B ϑ + log ( 2 d N ) ) , S=O(NB+B_H+B\vartheta+\log(2dN)), S = O ( N B + B H + B ϑ + log ( 2 d N )) ,
the master needs O ( N d + Γ ) O(Nd+\Gamma) O ( N d + Γ ) arithmetic operations on numbers of at most S S S bits. With conservative elementary integer arithmetic, its additional running time is O ( ( N d + Γ ) S 3 ) O((Nd+\Gamma)S^3) O (( N d + Γ ) S 3 ) , besides the N N N oracle invocations and the computation of returned payoff vectors. This is an explicit polynomial bound; faster arithmetic can improve it. The dependence on ρ \rho ρ is through B H = O ( bits ( ρ ) + bits ( ϑ ) + log ( 2 k ) ) B_H=O(\operatorname{bits}(\rho)+\operatorname{bits}(\vartheta)+\log(2k)) B H = O ( bits ( ρ ) + bits ( ϑ ) + log ( 2 k )) , not its numerical value. The same uniform query bound shows that a polynomial-time oracle has a uniform polynomial response bound B B B even if this bound is not specified separately in advance.
Why the horizon may be smaller than one. The proof only requires H > 0 H>0 H > 0 ; there is no lower bound on its numerical value. Imposing H ≥ 1 H\ge1 H ≥ 1 would unnecessarily introduce a dependence on the magnitude of ρ \rho ρ into the maximum exponent and the call count. Choosing H H H proportional to 1 / ρ 1/\rho 1/ ρ cancels that dependence without changing the normalized output distribution.
Corollary 8.2 (Support at most d d d without increasing the mean). Given the list returned by Theorem 8.1, an elementary procedure using only exact rational linear algebra finds a distribution supported on at most d d d of its vectors whose mean is no larger in any coordinate. Stopping one step earlier preserves the entire mean vector and uses at most d + 1 d+1 d + 1 vectors. Both procedures have polynomial running time; no linear program or additional oracle call is needed.
Proof. First, while more than d + 1 d+1 d + 1 coefficients are positive, find an integer relation ∑ t α t ( 1 , a ( t ) ) = 0 \sum_t \alpha_t(1,a^{(t)})=0 ∑ t α t ( 1 , a ( t ) ) = 0 on any d + 2 d+2 d + 2 supported columns. It has both signs. With θ = min α t > 0 λ t / α t \theta=\min_{\alpha_t>0}\lambda_t/\alpha_t θ = min α t > 0 λ t / α t , replace λ t \lambda_t λ t by λ t − θ α t \lambda_t-\theta\alpha_t λ t − θ α t . This keeps all probabilities nonnegative, preserves their sum and their full mean, and removes a supported vector. The columns never change, so rational Gaussian elimination gives integer relations of uniformly polynomial bit length. If λ t = N t / D \lambda_t=N_t/D λ t = N t / D and i i i attains the minimum, the shared-denominator update
N t ′ = N t α i − N i α t , D ′ = D α i N'_t=N_t\alpha_i-N_i\alpha_t,\qquad D'=D\alpha_i N t ′ = N t α i − N i α t , D ′ = D α i
adds at most bits ( α i ) \operatorname{bits}(\alpha_i) bits ( α i ) to the denominator length. The nonnegative new numerators sum to D ′ D' D ′ , so their lengths are bounded by the same quantity. There are at most N N N elimination steps, proving the stated bit bound for exact mean preservation.
If at most d d d coefficients remain, we are done. Otherwise there are exactly d + 1 d+1 d + 1 positive coefficients. If their augmented columns ( 1 , a ( t ) ) (1,a^{(t)}) ( 1 , a ( t ) ) are dependent, apply the same elimination once more. If they are independent, solve the nonsingular rational system
∑ t h t = 0 , ∑ t h t a ( t ) = − 1 d . \sum_t h_t=0,\qquad\sum_t h_ta^{(t)}=-\mathbf{1}_d. t ∑ h t = 0 , t ∑ h t a ( t ) = − 1 d .
The vector h h h is nonzero and has both signs because its entries sum to zero. Put
t ∗ = min h t < 0 λ t − h t > 0 , λ t ′ = λ t + t ∗ h t . t_* = \min_{h_t<0} \frac{\lambda_t}{-h_t} > 0,\qquad\lambda'_t = \lambda_t + t_*h_t. t ∗ = h t < 0 min − h t λ t > 0 , λ t ′ = λ t + t ∗ h t .
All new probabilities are nonnegative and sum to one, at least one becomes zero, and the mean changes by − t ∗ 1 d -t_*\mathbf{1}_d − t ∗ 1 d . Thus every coordinate decreases and at most d d d vectors remain. Because the probabilities remain a convex combination of nonnegative vectors, this step cannot make any mean coordinate negative.
For completeness, the final system has a polynomial-bit rational solution by determinant bounds on its unchanged input columns. Clearing its denominators gives an integer vector α \alpha α and a positive integer q q q such that ∑ t α t = 0 \sum_t \alpha_t = 0 ∑ t α t = 0 and ∑ t α t a ( t ) = − q 1 d \sum_t \alpha_t a^{(t)} = -q\mathbf{1}_d ∑ t α t a ( t ) = − q 1 d . If i i i minimizes λ i / ( − α i ) \lambda_i/(-\alpha_i) λ i / ( − α i ) over negative α i \alpha_i α i , the shared-denominator update for this last move is
N t ′ = N t ( − α i ) + N i α t , D ′ = D ( − α i ) . N'_t = N_t(-\alpha_i) + N_i\alpha_t,\qquad D' = D(-\alpha_i). N t ′ = N t ( − α i ) + N i α t , D ′ = D ( − α i ) .
Its numerators are nonnegative and sum to D ′ D' D ′ . It therefore has the same polynomial bit-growth bound as an ordinary elimination.
Constructing the final probabilities. Algorithm 3 gives the adaptive mixture and support reduction. Each oracle response contains a tree and its payoff vector; the final loop changes only their probabilities.
Algorithm 3. Rational mixture and support reduction
Input: A dimension d ≥ 1 d \ge1 d ≥ 1 , rational parameters ρ > 0 \rho> 0 ρ > 0 and 0 < ϑ ≤ 1 0 < \vartheta\le1 0 < ϑ ≤ 1 , and an exact oracle which, for each strictly positive rational probability vector p p p , returns ( T , v ) (T,v) ( T , v ) with v ∈ Q ≥ 0 d v \in\mathbb{Q}_{\ge0}^d v ∈ Q ≥ 0 d and p ⋅ v ≤ ρ p \cdot v \le\rho p ⋅ v ≤ ρ .
Output: Trees T t T_t T t and positive rational probabilities q t q_t q t with support at most d d d , ∑ t q t = 1 \sum_t q_t = 1 ∑ t q t = 1 , and ∑ t q t v ( t ) ≤ ( 1 + ϑ ) ρ 1 \sum_t q_t v^{(t)} \le(1+\vartheta)\rho\mathbf{1} ∑ t q t v ( t ) ≤ ( 1 + ϑ ) ρ 1 coordinatewise.
Compute k = min { j ∈ Z ≥ 0 : 2 j ≥ 2 d } k = \min\{j \in\mathbb{Z}_{\ge0} : 2^j \ge2d\} k = min { j ∈ Z ≥ 0 : 2 j ≥ 2 d } , b = 1 + ϑ / 4 b = 1 + \vartheta/4 b = 1 + ϑ /4 , and H = 10 k / ( ϑ 2 ρ ) H = 10k/(\vartheta^2\rho) H = 10 k / ( ϑ 2 ρ ) . Write b = P b / Q b b = P_b/Q_b b = P b / Q b in lowest terms, set J = ⌈ 20 k / ϑ 2 ⌉ J = \lceil20k/\vartheta^2\rceil J = ⌈ 20 k / ϑ 2 ⌉ , and precompute the integer powers P b s , Q b s P_b^s,Q_b^s P b s , Q b s for 0 ≤ s ≤ J 0 \le s \le J 0 ≤ s ≤ J . Initialize x = 0 ∈ Q d x = 0 \in\mathbb{Q}^d x = 0 ∈ Q d , τ = 0 \tau= 0 τ = 0 , and an empty response list.
While τ < H \tau< H τ < H : compute
e j = ⌈ x j ⌉ , w j = P b e j Q b J − e j , p j = w j / ∑ i w i , e_j = \lceil x_j\rceil,\qquad w_j = P_b^{e_j}Q_b^{J-e_j},\qquad p_j = w_j/\sum_i w_i, e j = ⌈ x j ⌉ , w j = P b e j Q b J − e j , p j = w j / i ∑ w i ,
and obtain ( T , v ) (T,v) ( T , v ) from the oracle. Put L = max j v j L = \max_j v_j L = max j v j and set Δ = H − τ \Delta= H-\tau Δ = H − τ if L = 0 L=0 L = 0 , or Δ = min { 1 / L , H − τ } \Delta= \min\{1/L,H-\tau\} Δ = min { 1/ L , H − τ } otherwise. Append ( Δ / H , T , v ) (\Delta/H,T,v) ( Δ/ H , T , v ) to the list and update x ← x + Δ v x \gets x+\Delta v x ← x + Δ v and τ ← τ + Δ \tau\gets\tau+\Delta τ ← τ + Δ .
Coalesce identical pairs ( T , v ) (T,v) ( T , v ) , adding their probabilities. While more than d d d positive coefficients remain, let A A A have the supported augmented columns ( 1 , v ( t ) ) (1,v^{(t)}) ( 1 , v ( t ) ) .
(a) If the columns are dependent, find a nonzero integer vector α \alpha α with A α = 0 A\alpha=0 A α = 0 , supported on at most d + 2 d+2 d + 2 columns. Set s = min α t > 0 q t / α t s=\min_{\alpha_t>0}q_t/\alpha_t s = min α t > 0 q t / α t and replace q t q_t q t by q t − s α t q_t-s\alpha_t q t − s α t .
(b) Otherwise there are exactly d + 1 d+1 d + 1 columns. Solve A h = ( 0 , − 1 d ) Ah=(0,-\mathbf{1}_d) A h = ( 0 , − 1 d ) , set s = min h t < 0 q t / ( − h t ) s=\min_{h_t<0}q_t/(-h_t) s = min h t < 0 q t / ( − h t ) , and replace q t q_t q t by q t + s h t q_t+sh_t q t + s h t .
(c) Delete every zero coefficient and its response.
Return the remaining pairs ( q t , T t ) (q_t,T_t) ( q t , T t ) .
The integer weights in Step 2 give exactly the normalized prices p j = b ⌈ x j ⌉ / ∑ i b ⌈ x i ⌉ p_j=b^{\lceil x_j\rceil}/\sum_i b^{\lceil x_i\rceil} p j = b ⌈ x j ⌉ / ∑ i b ⌈ x i ⌉ . Theorem 8.1 gives e j ≤ J e_j \le J e j ≤ J ; the precomputed powers therefore give O ( d ) O(d) O ( d ) arithmetic operations per query, as in the uniform-encoding analysis. For either graph application take d = m d=m d = m and the response bound ρ = ρ n \rho=\rho_n ρ = ρ n from § 11 . For distances, the oracle calls WeightedTree( G , ℓ , ( p e / ℓ e ) e ) (G,\ell,(p_e/\ell_e)_e) ( G , ℓ , ( p e / ℓ e ) e ) and returns the stretch vector v = s T v=s_T v = s T of (25) . For capacities, it calls WeightedTree( G , ( p e / c e ) e , c ) (G,(p_e/c_e)_e,c) ( G , ( p e / c e ) e , c ) and returns the fundamental-cut congestion vector v = g T v=g_T v = g T of (26) . These payoffs are evaluated against the fixed original length or capacity vector, respectively. Each weighted-tree call constructs its own measure μ \mu μ , which remains fixed within that call. The output probabilities define q d i s t q^{\mathrm{dist}} q dist or q c a p q^{\mathrm{cap}} q cap ; in the latter case assign every supported tree edge its fundamental-cut capacity from (27) . The call bound and polynomial running time follow from Theorem 8.1 ; the last loop is precisely the rational support reduction of Corollary 8.2 .
Fundamental-cut capacities and routing
Theorem 6.5 supplies the weighted-tree hypothesis used below. We state the reduction for a generic response bound ρ \rho ρ . The distance–capacity equivalence is the family-preserving framework of Andersen–Feige [3 ] , following Räcke [22 ] ; we reproduce the identities and routing argument to specify all quantifiers.
The weighted hypothesis and the two payoff vectors
Let G = ( V , E ) G=(V,E) G = ( V , E ) be a finite connected undirected graph, with n ≥ 2 n\ge2 n ≥ 2 and m = ∣ E ∣ ≥ 1 m=|E|\ge1 m = ∣ E ∣ ≥ 1 . Loops may be deleted. Parallel edges cause no problem; a reduction to the simple-graph weighted hypothesis is given below. Write T ( G ) \mathcal{T}(G) T ( G ) for the finite collection of its actual spanning trees. The weighted hypothesis with parameter ρ ≥ 1 \rho\ge1 ρ ≥ 1 is the following assertion:
for all ℓ e > 0 , a e ≥ 0 , min T ∈ T ( G ) ∑ e = u v a e d T , ℓ ( u , v ) ≤ ρ ∑ e a e ℓ e . (24) \text{for all }\ell_e>0,\ a_e\ge0,\quad\min_{T\in\mathcal{T}(G)}\sum_{e=uv}a_e d_{T,\ell}(u,v)\le\rho\sum_e a_e\ell_e. \tag*{(24)} for all ℓ e > 0 , a e ≥ 0 , T ∈ T ( G ) min e = uv ∑ a e d T , ℓ ( u , v ) ≤ ρ e ∑ a e ℓ e . (24)
Theorem 6.5 gives a deterministic oracle for this hypothesis with ρ = ρ n \rho=\rho_n ρ = ρ n .
For a tree T T T and original edge e = u v e=uv e = uv , let P T ( e ) P_T(e) P T ( e ) be the unique u u u –v v v path in T T T . Set
M T ( e , f ) = 1 { f ∈ P T ( e ) } , e , f ∈ E . M_T(e,f)=\mathbf{1}_{\{f\in P_T(e)\}},\qquad e,f\in E. M T ( e , f ) = 1 { f ∈ P T ( e )} , e , f ∈ E .
For fixed original lengths ℓ > 0 \ell>0 ℓ > 0 the stretch payoff is
s T ( e ) = d T , ℓ ( e ) ℓ e = ∑ f M T ( e , f ) ℓ f ℓ e . (25) s_T(e)=\frac{d_{T,\ell}(e)}{\ell_e}=\frac{\sum_f M_T(e,f)\ell_f}{\ell_e}. \tag*{(25)} s T ( e ) = ℓ e d T , ℓ ( e ) = ℓ e ∑ f M T ( e , f ) ℓ f . (25)
For fixed original capacities c > 0 c>0 c > 0 the load and congestion payoffs are
L T ( f ) = ∑ e c e M T ( e , f ) , g T ( f ) = L T ( f ) c f . (26) L_T(f)=\sum_e c_eM_T(e,f),\qquad g_T(f)=\frac{L_T(f)}{c_f}. \tag*{(26)} L T ( f ) = e ∑ c e M T ( e , f ) , g T ( f ) = c f L T ( f ) . (26)
In particular L T ( f ) = 0 L_T(f)=0 L T ( f ) = 0 if f ∉ T f\notin T f ∈ / T . If f ∈ T f\in T f ∈ T and S f S_f S f is either component of T − f T-f T − f , then
L T ( f ) = c G ( δ G ( S f ) ) . (27) L_T(f)=c_G(\delta_G(S_f)). \tag*{(27)} L T ( f ) = c G ( δ G ( S f )) . (27)
We equip the spanning-tree topology with capacities c T ( f ) = L T ( f ) c_T(f)=L_T(f) c T ( f ) = L T ( f ) for f ∈ T f\in T f ∈ T . The tree uses original vertices and edges, while its capacities are the capacities of their fundamental cuts; in general c T ( f ) c_T(f) c T ( f ) differs from the original c f c_f c f . The probabilistic congestion parameter is max f E g T ( f ) \max_f\mathbb{E}g_T(f) max f E g T ( f ) , with the expectation taken before the maximum. A bound on this parameter need not bound E max f g T ( f ) \mathbb{E}\max_f g_T(f) E max f g T ( f ) for a single sampled tree. The routing guarantee below applies to the fixed fractional mixture of tree paths.
Lemma 9.1 (The weighted distance–capacity identity). For every strictly positive edge-price vector w w w and every tree T T T , the choices
ℓ f w = w f c f , a e = c e \ell_f^w=\frac{w_f}{c_f},\qquad a_e=c_e ℓ f w = c f w f , a e = c e satisfy
∑ e a e d T , ℓ w ( e ) = ∑ f w f g T ( f ) , ∑ e a e ℓ e w = ∑ e w e . (28) \sum_e a_e d_{T,\ell^w}(e)=\sum_f w_f g_T(f),\qquad\sum_e a_e\ell_e^w=\sum_e w_e. \tag*{(28)} e ∑ a e d T , ℓ w ( e ) = f ∑ w f g T ( f ) , e ∑ a e ℓ e w = e ∑ w e . (28) For fixed lengths ℓ \ell ℓ , the alternative choices a e = w e / ℓ e a_e=w_e/\ell_e a e = w e / ℓ e satisfy
∑ e a e d T , ℓ ( e ) = ∑ e w e s T ( e ) , ∑ e a e ℓ e = ∑ e w e . (29) \sum_e a_e d_{T,\ell}(e)=\sum_e w_e s_T(e),\qquad\sum_e a_e\ell_e=\sum_e w_e. \tag*{(29)} e ∑ a e d T , ℓ ( e ) = e ∑ w e s T ( e ) , e ∑ a e ℓ e = e ∑ w e . (29) Proof. Exchange the two finite sums:
∑ e c e ∑ f M T ( e , f ) w f c f = ∑ f w f c f ∑ e c e M T ( e , f ) . \sum_e c_e\sum_f M_T(e,f)\frac{w_f}{c_f}
=\sum_f\frac{w_f}{c_f}\sum_e c_eM_T(e,f). e ∑ c e f ∑ M T ( e , f ) c f w f = f ∑ c f w f e ∑ c e M T ( e , f ) . The second identity in (28) and both identities in (29) are direct substitutions. No tree transformation, Steiner vertex, or change of the admissible family occurs.
Exact existential bounds
Theorem 9.2 (Two distributions on the same family of spanning trees). Assume (24) . For every fixed ℓ > 0 \ell>0 ℓ > 0 there is a distribution q d i s t q^{\mathrm{dist}} q dist on T ( G ) \mathcal{T}(G) T ( G ) such that
E q d i s t d T , ℓ ( u , v ) ≤ ρ d G , ℓ ( u , v ) ( u , v ∈ V ) . \mathbb{E}_{q^{\mathrm{dist}}}d_{T,\ell}(u,v)\le\rho d_{G,\ell}(u,v)\qquad(u,v\in V). E q dist d T , ℓ ( u , v ) ≤ ρ d G , ℓ ( u , v ) ( u , v ∈ V ) . Every sampled tree also satisfies d T , ℓ ( u , v ) ≥ d G , ℓ ( u , v ) d_{T,\ell}(u,v)\ge d_{G,\ell}(u,v) d T , ℓ ( u , v ) ≥ d G , ℓ ( u , v ) . For every fixed c > 0 c>0 c > 0 there is a distribution q c a p q^{\mathrm{cap}} q cap on the same family such that
E q c a p L T ( f ) ≤ ρ c f ( f ∈ E ) . (30) \mathbb{E}_{q^{\mathrm{cap}}}L_T(f)\le\rho c_f\qquad(f\in E). \tag*{(30)} E q cap L T ( f ) ≤ ρ c f ( f ∈ E ) . (30) The two distributions are allowed to differ.
Proof. For the distance assertion, let p p p be any probability vector on E E E . Use a e = p e / ℓ e a_e=p_e/\ell_e a e = p e / ℓ e in (24) . This gives min T p ⋅ s T ≤ ρ \min_T p\cdot s_T\le\rho min T p ⋅ s T ≤ ρ . The finite minimax theorem produces a distribution satisfying E s T ( e ) ≤ ρ \mathbb{E}s_T(e)\le\rho E s T ( e ) ≤ ρ for every edge. Summing the tree triangle inequality along a fixed shortest graph path proves the all-pairs assertion. Each tree path is a graph path with the original lengths, which proves domination.
For capacities, first suppose p f > 0 p_f>0 p f > 0 for every edge and use Lemma 9.1 with w = p w=p w = p . Then min T p ⋅ g T ≤ ρ \min_T p\cdot g_T\le\rho min T p ⋅ g T ≤ ρ . To allow zero prices, for any η > 0 \eta>0 η > 0 use w f = p f + η / m w_f=p_f+\eta/m w f = p f + η / m . Nonnegativity gives
min T p ⋅ g T ≤ min T w ⋅ g T ≤ ρ ( 1 + η ) . \min_T p\cdot g_T\le\min_T w\cdot g_T\le\rho(1+\eta). T min p ⋅ g T ≤ T min w ⋅ g T ≤ ρ ( 1 + η ) . Let η ↓ 0 \eta\downarrow0 η ↓ 0 . The left side is fixed, so it is at most ρ \rho ρ . Finite minimax now gives E g T ( f ) ≤ ρ \mathbb{E}g_T(f)\le\rho E g T ( f ) ≤ ρ for all f f f , which is (30) .
All-cut domination and the full routing statement
Proposition 9.3 (Every tree dominates all cuts). For every spanning-tree T T T equipped with (27) , and every A ⊆ V A \subseteq V A ⊆ V ,
c G ( δ G ( A ) ) ≤ c T ( δ T ( A ) ) . (31) c_G(\delta_G(A)) \le c_T(\delta_T(A)). \tag*{(31)} c G ( δ G ( A )) ≤ c T ( δ T ( A )) . (31) Any distribution with E L T ( f ) ≤ α c f \mathbb{E}L_T(f) \le\alpha c_f E L T ( f ) ≤ α c f for every original edge also satisfies
E c T ( δ T ( A ) ) ≤ α c G ( δ G ( A ) ) ( A ⊆ V ) . (32) \mathbb{E}c_T(\delta_T(A)) \le\alpha c_G(\delta_G(A)) \qquad(A \subseteq V). \tag*{(32)} E c T ( δ T ( A )) ≤ α c G ( δ G ( A )) ( A ⊆ V ) . (32)
Proof. By interchanging finite sums,
c T ( δ T ( A ) ) = ∑ e c e ∣ P T ( e ) ∩ δ T ( A ) ∣ . c_T(\delta_T(A)) = \sum_e c_e \lvert P_T(e) \cap\delta_T(A) \rvert. c T ( δ T ( A )) = e ∑ c e ∣ P T ( e ) ∩ δ T ( A )∣ .
If e e e crosses A A A , its tree path crosses A A A at least once. This proves (31) ; paths for other edges may contribute additional nonnegative terms. For the upper bound,
c T ( δ T ( A ) ) = ∑ f ∈ δ G ( A ) L T ( f ) , c_T(\delta_T(A)) = \sum_{f \in\delta_G(A)} L_T(f), c T ( δ T ( A )) = f ∈ δ G ( A ) ∑ L T ( f ) ,
where non-tree edges contribute zero. Take expectations and apply the individual edge-load bounds. □
Corollary 9.4 (Cuts of prescribed cardinality). Given the capacity distribution of Theorem 1.1 , for every 1 ≤ s ≤ n − 1 1 \le s \le n-1 1 ≤ s ≤ n − 1 one can find in polynomial time a set S S S of exactly s s s vertices whose graph-cut capacity is at most α n \alpha_n α n times the minimum possible capacity of such a set. In particular, for even n n n this gives an O ( log n ) O(\log n) O ( log n ) approximation for minimum bisection with exact balance.
Proof. For each supported tree T T T , compute a minimum-capacity tree cut S T S_T S T with ∣ S T ∣ = s \lvert S_T\rvert= s ∣ S T ∣ = s . Return the candidate with the smallest original graph-cut capacity. If S ∗ S^* S ∗ is an optimal graph cut of cardinality s s s , the expectation bound implies that some supported tree satisfies
c T ( δ T ( S ∗ ) ) ≤ α n c G ( δ G ( S ∗ ) ) . c_T(\delta_T(S^*)) \le\alpha_n c_G(\delta_G(S^*)). c T ( δ T ( S ∗ )) ≤ α n c G ( δ G ( S ∗ )) .
For that tree, domination and the definition of S T S_T S T give
c G ( δ G ( S T ) ) ≤ c T ( δ T ( S T ) ) ≤ c T ( δ T ( S ∗ ) ) ≤ α n c G ( δ G ( S ∗ ) ) . c_G(\delta_G(S_T)) \le c_T(\delta_T(S_T)) \le c_T(\delta_T(S^*)) \le\alpha_n c_G(\delta_G(S^*)). c G ( δ G ( S T )) ≤ c T ( δ T ( S T )) ≤ c T ( δ T ( S ∗ )) ≤ α n c G ( δ G ( S ∗ )) .
It remains to describe the tree computation.
Root T T T arbitrarily. The state D P v ( σ , j ) \mathrm{DP}_v(\sigma,j) DP v ( σ , j ) is the minimum cut capacity inside the subtree at v v v , when exactly j j j of its vertices belong to the chosen set and v v v has membership σ ∈ { 0 , 1 } \sigma\in\{0,1\} σ ∈ { 0 , 1 } . Process the children w 1 , … , w d w_1,\ldots,w_d w 1 , … , w d of v v v in any fixed order. Before processing them, set D P v , 0 ( σ , j ) = 0 \mathrm{DP}_{v,0}(\sigma,j)=0 DP v , 0 ( σ , j ) = 0 if j = σ j=\sigma j = σ , and + ∞ +\infty + ∞ otherwise. Then use the knapsack recurrence
D P v , i ( σ , j ) = min τ ∈ { 0 , 1 } 0 ≤ h ≤ ∣ V ( T w i ) ∣ { D P v , i − 1 ( σ , j − h ) + D P w i ( τ , h ) + c T ( v w i ) 1 { σ ≠ τ } } , \mathrm{DP}_{v,i}(\sigma,j) =
\min_{\substack{\tau\in\{0,1\}\\0 \le h \le\lvert V(T_{w_i})\rvert}}
\left\{
\mathrm{DP}_{v,i-1}(\sigma,j-h)
+\mathrm{DP}_{w_i}(\tau,h)
+c_T(vw_i)\mathbf{1}_{\{\sigma\ne\tau\}}
\right\}, DP v , i ( σ , j ) = τ ∈ { 0 , 1 } 0 ≤ h ≤ ∣ V ( T w i )∣ min { DP v , i − 1 ( σ , j − h ) + DP w i ( τ , h ) + c T ( v w i ) 1 { σ = τ } } ,
with invalid counts assigned + ∞ +\infty + ∞ , and set D P v = D P v , d \mathrm{DP}_v=\mathrm{DP}_{v,d} DP v = DP v , d . The optimum at the root r r r is
min σ ∈ { 0 , 1 } D P r ( σ , s ) . \min_{\sigma\in\{0,1\}} \mathrm{DP}_r(\sigma,s). σ ∈ { 0 , 1 } min DP r ( σ , s ) .
Storing a minimizing choice in each state recovers S T S_T S T . The computation uses at most O ( n 3 ) O(n^3) O ( n 3 ) rational arithmetic operations per tree. Evaluating its graph cut takes O ( m ) O(m) O ( m ) more operations. The rational capacities have polynomial encodings, and the support has at most m m m trees, so the complete procedure runs in polynomial time. □
Räcke [22 ] established the logarithmic approximation for minimum bisection using hierarchical decomposition trees and a dynamic program that counts terminals at the leaves. The corollary recovers that approximation from spanning trees; here the cardinality state counts the original graph vertices directly.
Proposition 9.5 (Canonical capacities are the smallest cut dominator). Fix an actual spanning-tree topology T T T . A nonnegative capacity assignment b T b_T b T on its edges dominates every cut of G G G if and only if
b T ( f ) ≥ c G ( δ G ( S f ) ) ( f ∈ T ) . b_T(f) \ge c_G(\delta_G(S_f)) \qquad(f \in T). b T ( f ) ≥ c G ( δ G ( S f )) ( f ∈ T ) . Thus the fundamental-cut assignment c T c_T c T is the coordinatewise smallest capacity assignment on that topology with all-cut domination.
Proof. For necessity, the fundamental cut S f S_f S f crosses exactly one edge of T T T , namely f f f . Its domination inequality is precisely the displayed inequality. For sufficiency, the displayed inequalities give b T ( δ T ( A ) ) ≥ c T ( δ T ( A ) ) b_T(\delta_T(A)) \ge c_T(\delta_T(A)) b T ( δ T ( A )) ≥ c T ( δ T ( A )) for every cut A A A ; Proposition 9.3 supplies c T ( δ T ( A ) ) ≥ c G ( δ G ( A ) ) c_T(\delta_T(A)) \ge c_G(\delta_G(A)) c T ( δ T ( A )) ≥ c G ( δ G ( A )) . □ \square □
For precision, a demand vector has one nonnegative coordinate D u v D_{uv} D uv for each unordered pair of distinct vertices. Let F ( G , c ) \mathcal{F}(G,c) F ( G , c ) be the set of such vectors routable in G G G by a fractional multicommodity flow of congestion at most one. Let opt G , c ( D ) \operatorname{opt}_{G,c}(D) opt G , c ( D ) be its minimum possible congestion. For a cut S S S , write
D ( S ) = ∑ { u , v } : ∣ { u , v } ∩ S ∣ = 1 D u v . D(S)=\sum_{\{u,v\}:|\{u,v\}\cap S|=1}D_{uv}. D ( S ) = { u , v } : ∣ { u , v } ∩ S ∣ = 1 ∑ D uv .
Every flow routing D D D at congestion λ \lambda λ satisfies the necessary cut inequalities D ( S ) ≤ λ c G ( δ G ( S ) ) D(S)\le\lambda c_G(\delta_G(S)) D ( S ) ≤ λ c G ( δ G ( S )) .
Theorem 9.6 (Tree-flow mixtures and an oblivious routing operator). Suppose a distribution q q q satisfies E q L T ( f ) ≤ α c f \mathbb{E}_q L_T(f)\le\alpha c_f E q L T ( f ) ≤ α c f for all f ∈ E f\in E f ∈ E . Then the following assertions hold.
For every tree T T T and every demand vector D D D ,
opt T , c T ( D ) ≤ opt G , c ( D ) . \operatorname{opt}_{T,c_T}(D)\le\operatorname{opt}_{G,c}(D). opt T , c T ( D ) ≤ opt G , c ( D ) . In particular F ( G , c ) ⊆ F ( T , c T ) \mathcal{F}(G,c)\subseteq\mathcal{F}(T,c_T) F ( G , c ) ⊆ F ( T , c T ) .
For arbitrary demand vectors D T ∈ F ( T , c T ) D^T\in\mathcal{F}(T,c_T) D T ∈ F ( T , c T ) , one for each supported tree, the averaged vector ∑ T q T D T \sum_T q_TD^T ∑ T q T D T is routable in G G G at congestion at most α \alpha α . Consequently, with the sum denoting the weighted Minkowski sum,
F ( G , c ) ⊆ ∑ T q T F ( T , c T ) ⊆ α F ( G , c ) . (33) \mathcal{F}(G,c)\subseteq\sum_T q_T\mathcal{F}(T,c_T)\subseteq\alpha\mathcal{F}(G,c). \tag*{(33)} F ( G , c ) ⊆ T ∑ q T F ( T , c T ) ⊆ α F ( G , c ) . (33) Route each pair u , v u,v u , v by putting fraction q T q_T q T of its demand on P T ( u , v ) P_T(u,v) P T ( u , v ) . This defines a fixed fractional oblivious routing operator, independent of D D D , whose congestion for every D D D is at most α opt G , c ( D ) \alpha\operatorname{opt}_{G,c}(D) α opt G , c ( D ) .
Proof. Unique path routing in a tree puts load exactly D ( S f ) D(S_f) D ( S f ) on its edge f f f . If D D D can be routed in G G G at congestion λ \lambda λ , the cut inequality for S f S_f S f gives
D ( S f ) ≤ λ c G ( δ G ( S f ) ) = λ c T ( f ) . D(S_f)\le\lambda c_G(\delta_G(S_f))=\lambda c_T(f). D ( S f ) ≤ λ c G ( δ G ( S f )) = λ c T ( f ) .
Taking the minimum possible λ \lambda λ proves the first assertion. For the second assertion, route D T D^T D T on the original paths of T T T and scale that flow by q T q_T q T . Its load on an original edge f f f is at most q T L T ( f ) q_TL_T(f) q T L T ( f ) . Summing the flows gives the required averaged demand vector and edge load at most ∑ T q T L T ( f ) ≤ α c f \sum_T q_TL_T(f) \le\alpha c_f ∑ T q T L T ( f ) ≤ α c f . Taking D T = D D^T = D D T = D for all T T T also proves the first inclusion in (33) .
Finally, for a fixed D D D and λ = opt G , c ( D ) \lambda= \operatorname{opt}_{G,c}(D) λ = opt G , c ( D ) , the first assertion bounds every tree’s routing load by λ c T ( f ) \lambda c_T(f) λ c T ( f ) . The second argument then gives averaged edge load at most λ α c f \lambda\alpha c_f λ α c f . All pair paths and mixture coefficients were fixed before D D D was specified, as required for oblivious routing.
These are statements about one fixed distribution working for every cut and demand vector. The cut upper bound is an expectation for each fixed cut. It does not assert that one sampled tree simultaneously has distortion α \alpha α on all cuts. A demand feasible in a single supported tree need not be feasible in G G G at congestion α \alpha α ; the weighted mixture statement in Theorem 9.6 is the relevant converse.
Proof of Theorem 1.1 . Apply Theorem 6.5 with bound ρ n \rho_n ρ n . For the distance instance, a positive master query p p p uses a e = p e / ℓ e a_e = p_e/\ell_e a e = p e / ℓ e and returns the stretch vector s T s_T s T . For the capacity instance use ℓ e = p e / c e \ell_e = p_e/c_e ℓ e = p e / c e and a e = c e a_e = c_e a e = c e , which returns the congestion vector g T g_T g T by Lemma 9.1 . Both queries have baseline one. Theorem 8.1 therefore constructs an explicit rational distribution whose every coordinate is at most ( 1 + ϑ ) ρ n (1+\vartheta)\rho_n ( 1 + ϑ ) ρ n . Corollary 8.2 reduces its support to at most m m m without increasing any expected coordinate.
For distances, summing the tree triangle inequality along a fixed shortest graph path gives the all-pairs bound. Every tree path is an original graph path, so metric domination is deterministic. For capacities, Proposition 9.3 and Theorem 9.6 give the cut and routing conclusions with the same factor. All graph payoffs are sums of fixed input rationals divided by a fixed positive input rational. § 7 and the master’s bit analysis give the claimed running time, probabilities, and call bound. For positive real input, approximate the lengths or capacities by positive rationals and pass to a convergent subsequence of distributions on the finite tree family. There are only finitely many supports of size at most m m m , so the limit retains this support bound. All distance and load inequalities are continuous in the input.
Computing all fundamental-cut capacities. Root the returned tree and preprocess lowest common ancestors. For every graph edge u v uv uv , add c u v c_{uv} c uv at u u u and v v v , and subtract 2 c u v 2c_{uv} 2 c uv at their lowest common ancestor. Taking subtree sums makes the sum below each nonroot vertex equal to the capacity of its parent edge’s fundamental cut. This computes all tree capacities in O ( m log n + n ) O(m\log n+n) O ( m log n + n ) arithmetic operations.
Parallel edges and zero capacities. For parallel edges with common endpoints, retain one edge of minimum auxiliary length and replace their demands by their sum. The weighted tree cost is unchanged, while its baseline weakly decreases. Therefore a weighted theorem proved for simple graphs extends to multigraphs with the same constant, and the chosen representative is an original edge. For nonnegative capacities, first form the connected components of the positive-capacity subgraph. Run the capacity construction inside each such component and join those component trees by a fixed spanning tree of zero-capacity edges in the quotient. A joining edge has fundamental-cut capacity zero. A joint oracle call applies the weighted oracle separately to every positive component; their weighted costs add. The rational master on the coordinates of positive-capacity edges consequently produces a polynomial-support distribution for the whole graph, without taking an exponentially large product of component supports. The all-zero-capacity case is immediate. A one-vertex graph is also immediate.
Lower bounds
We first prove logarithmic lower bounds on diamonds, matching the three approximation guarantees. We then prove the support tradeoff and explain why the distance and capacity distributions may differ. None of these lower bounds uses the weighted upper theorem.
Diamond graphs force logarithmic stretch and cut distortion
The logarithmic order is already necessary on unweighted series-parallel graphs. Gupta, Newman, Rabinovich, and Sinclair [18 ] proved this for probabilistic approximation by arbitrary dominating tree metrics, using recursive diamond graphs. The argument below is a direct proof for spanning trees. Its uniform edge-average form also proves the corresponding lower bounds for fundamental-cut capacities.
Let D 0 D_0 D 0 consist of one edge with designated terminals. For h ≥ 1 h \ge1 h ≥ 1 , form D h D_h D h from four disjoint copies of D h − 1 D_{h-1} D h − 1 by putting two copies in series on each of two parallel branches. All edges of the resulting graph have length one and capacity one. Thus
m h = ∣ E ( D h ) ∣ = 4 h , n h = ∣ V ( D h ) ∣ = 2 ⋅ 4 h + 4 3 . m_h = |E(D_h)| = 4^h,\qquad n_h = |V(D_h)| = \frac{2\cdot4^h+4}{3}. m h = ∣ E ( D h ) ∣ = 4 h , n h = ∣ V ( D h ) ∣ = 3 2 ⋅ 4 h + 4 .
The canonical copies of D h D_h D h inside D k D_k D k will be called blocks of height h h h ; a block at depth i i i has height k − i k-i k − i . Only the two terminals of a block can have incident edges outside that block.
Proposition 10.1 (Uniform average stretch on diamonds). *For every k ≥ 1 k \ge1 k ≥ 1 and every spanning tree T ⊆ D k T \subseteq D_k T ⊆ D k ,
∑ e = u v ∈ E ( D k ) d T ( u , v ) ≥ k 4 k 6 − 4 k − 1 18 ≥ k 4 k 9 . \sum_{e=uv\in E(D_k)} d_T(u,v) \ge\frac{k4^k}{6}-\frac{4^k-1}{18}\ge\frac{k4^k}{9}. e = uv ∈ E ( D k ) ∑ d T ( u , v ) ≥ 6 k 4 k − 18 4 k − 1 ≥ 9 k 4 k .
Proof. We first record two features of a height-h h h block B B B with terminals a , b a,b a , b . Give the two terminals of D k D_k D k integer coordinates 0 0 0 and 2 k 2^k 2 k . Recursively, when the terminals of a height-h h h block have coordinates x x x and x + 2 h x+2^h x + 2 h , give both of its intermediate branch vertices coordinate x + 2 h − 1 x+2^{h-1} x + 2 h − 1 , and continue within the four children. This defines a single coordinate on every vertex of D k D_k D k . Every original edge changes the coordinate by exactly one. Consequently every path in the whole graph between a a a and b b b has length at least 2 h 2^h 2 h , even if it leaves B B B .
The block B B B also contains 2 h 2^h 2 h edge-disjoint paths between its terminals. This follows by induction: each pair of series children supplies 2 h − 1 2^{h-1} 2 h − 1 paths, by concatenating their edge-disjoint path families at the common terminal; the two parallel branches are edge-disjoint. Hence every cut inside B B B separating its terminals contains at least 2 h 2^h 2 h edges.
Fix T T T , and let T B T_B T B be the forest whose vertices are those of B B B and whose edges are E ( T ) ∩ E ( B ) E(T) \cap E(B) E ( T ) ∩ E ( B ) . Every component of T B T_B T B contains a terminal of B B B . Indeed, a component with neither terminal has no incident edge leaving B B B , and has no tree edge to another component inside B B B , so it would be disconnected from the rest of T T T . It follows that T B T_B T B has either one or two components. Call B B B closed in the first case and open in the second. An open block has one terminal in each of its two components.
Suppose a nonleaf block has c ∈ { 1 , 2 } c \in\{1,2\} c ∈ { 1 , 2 } components in its restricted forest and o o o open children. Each closed child contributes n h − 1 − 1 n_{h-1}-1 n h − 1 − 1 tree edges, and each open child contributes n h − 1 − 2 n_{h-1}-2 n h − 1 − 2 . Since ∣ V ( B ) ∣ = 4 n h − 1 − 4 |V(B)| = 4n_{h-1} - 4 ∣ V ( B ) ∣ = 4 n h − 1 − 4 , counting edges gives
∣ E ( T B ) ∣ = ( 4 − o ) ( n h − 1 − 1 ) + o ( n h − 1 − 2 ) = ∣ V ( B ) ∣ − o . |E(T_B)| = (4-o)(n_{h-1}-1) + o(n_{h-1}-2) = |V(B)| - o. ∣ E ( T B ) ∣ = ( 4 − o ) ( n h − 1 − 1 ) + o ( n h − 1 − 2 ) = ∣ V ( B ) ∣ − o .
The same forest has ∣ V ( B ) ∣ − c |V(B)|-c ∣ V ( B ) ∣ − c edges. Thus o = c o=c o = c : a closed parent has one open child, and an open parent has two (Figure 2 ).
Figure 2. The recursive diamond and the component count. Each side represents a height-( h − 1 ) (h-1) ( h − 1 ) child block of a height-h h h diamond. A solid side denotes a connected restriction of the spanning tree; a dashed side denotes a restriction with two components, one at each child terminal. A parent with c c c components has exactly c c c open children. The two drawings show the resulting forests on the four corner vertices.
Let O i O_i O i denote the number of open blocks at depth i i i . The root is closed, and there are 4 i 4^i 4 i blocks at depth i i i . The preceding count therefore gives
O 0 = 0 , O i + 1 = ( 4 i − O i ) + 2 O i = 4 i + O i , O i = 4 i − 1 3 . O_0 = 0,\qquad O_{i+1} = (4^i-O_i)+2O_i = 4^i+O_i,\qquad O_i = \frac{4^i-1}{3}. O 0 = 0 , O i + 1 = ( 4 i − O i ) + 2 O i = 4 i + O i , O i = 3 4 i − 1 .
For each open block B B B of height h h h , consider the cut between the two components of T B T_B T B and charge 2 h 2^h 2 h to every edge of B B B crossing that cut. There are at least 2 h 2^h 2 h such edges, so this block creates total charge at least 4 h 4^h 4 h .
If e = u v e=uv e = uv receives this charge, its endpoints lie in different components of T B T_B T B . The unique T T T -path from u u u to v v v must pass through both terminals of B B B : the only way out of either component is through its own terminal. The integer coordinate just constructed shows that the terminal-to-terminal part of this path has length at least 2 h 2^h 2 h . In particular, d T ( u , v ) ≥ 2 h d_T(u,v)\ge2^h d T ( u , v ) ≥ 2 h .
The blocks containing a fixed original edge form a single ancestor chain, with one block of each height. The charges received by that edge thus have distinct values among 1 , 2 , 4 , … , 2 k − 1 1,2,4,\ldots,2^{k-1} 1 , 2 , 4 , … , 2 k − 1 . Their sum is less than twice the largest received charge, and that largest charge is at most d T ( u , v ) d_T(u,v) d T ( u , v ) . (An edge receiving no charge contributes zero.) Writing Q Q Q for the total charge gives
2 ∑ e = u v ∈ E ( D k ) d T ( u , v ) ≥ Q ≥ ∑ i = 1 k O i 4 k − i = k 4 k 3 − 4 k − 1 9 . 2\sum_{e=uv\in E(D_k)} d_T(u,v) \ge Q \ge\sum_{i=1}^{k} O_i4^{k-i} = \frac{k4^k}{3} - \frac{4^k-1}{9}. 2 e = uv ∈ E ( D k ) ∑ d T ( u , v ) ≥ Q ≥ i = 1 ∑ k O i 4 k − i = 3 k 4 k − 9 4 k − 1 .
Dividing by two proves the first inequality. The second follows because their difference is ( ( k − 1 ) 4 k + 1 ) / 18 ≥ 0 ((k-1)4^k+1)/18 \ge0 (( k − 1 ) 4 k + 1 ) /18 ≥ 0 .
Corollary 10.2 (Three logarithmic lower bounds on the same diamonds). Let D D D be any distribution over spanning trees of D k D_k D k , with unit original edge lengths and capacities and with the fundamental-cut
capacities defined above. Then
max u v ∈ E ( D k ) E T ∼ D d T ( u , v ) d D k ( u , v ) ≥ k 9 , \max_{uv\in E(D_k)} \mathbb{E}_{T\sim\mathcal{D}}\frac{d_T(u,v)}{d_{D_k}(u,v)} \ge\frac{k}{9}, uv ∈ E ( D k ) max E T ∼ D d D k ( u , v ) d T ( u , v ) ≥ 9 k ,
max f ∈ E ( D k ) E T ∼ D g T ( f ) ≥ k 9 , \max_{f\in E(D_k)} \mathbb{E}_{T\sim\mathcal{D}}g_T(f) \ge\frac{k}{9}, f ∈ E ( D k ) max E T ∼ D g T ( f ) ≥ 9 k ,
max v ∈ V ( D k ) E T ∼ D c T ( δ T ( { v } ) ) c D k ( δ D k ( { v } ) ) ≥ k 9 . \max_{v\in V(D_k)}\frac{\mathbb{E}_{T\sim\mathcal{D}}c_T(\delta_T(\{v\}))}{c_{D_k}(\delta_{D_k}(\{v\}))} \ge\frac{k}{9}. v ∈ V ( D k ) max c D k ( δ D k ({ v })) E T ∼ D c T ( δ T ({ v })) ≥ 9 k .
In particular, expected stretch, maximum expected congestion, and maximum expected cut distortion each require order log n k \log n_k log n k . The cut lower bound already holds on singleton cuts.
Proof. For the metric assertion, d D k ( u , v ) = 1 d_{D_k}(u,v)=1 d D k ( u , v ) = 1 for every original edge. Average the proposition over T T T , and then average over the m k m_k m k original edges.
For capacities, an original edge u v uv uv crosses the fundamental cut of f ∈ T f\in T f ∈ T exactly when the T T T -path from u u u to v v v uses f f f . Double counting gives the pointwise identity
∑ f ∈ E ( D k ) L T ( f ) = ∑ u v ∈ E ( D k ) ∣ E ( P T ( u , v ) ) ∣ = ∑ u v ∈ E ( D k ) d T ( u , v ) . \sum_{f\in E(D_k)}L_T(f)=\sum_{uv\in E(D_k)}\lvert E(P_T(u,v))\rvert=\sum_{uv\in E(D_k)}d_T(u,v). f ∈ E ( D k ) ∑ L T ( f ) = uv ∈ E ( D k ) ∑ ∣ E ( P T ( u , v ))∣ = uv ∈ E ( D k ) ∑ d T ( u , v ) .
All c f = 1 c_f=1 c f = 1 . Averaging over T T T and then over f f f proves the congestion assertion.
Finally, sum the singleton tree cuts. Every selected tree edge is counted at its two endpoints, so
∑ v ∈ V ( D k ) c T ( δ T ( { v } ) ) = 2 ∑ f L T ( f ) , ∑ v ∈ V ( D k ) deg D k ( v ) = 2 m k . \sum_{v\in V(D_k)}c_T(\delta_T(\{v\}))=2\sum_f L_T(f),\qquad\sum_{v\in V(D_k)}\deg_{D_k}(v)=2m_k. v ∈ V ( D k ) ∑ c T ( δ T ({ v })) = 2 f ∑ L T ( f ) , v ∈ V ( D k ) ∑ deg D k ( v ) = 2 m k .
Therefore the degree-weighted average of the expected singleton-cut ratios is
∑ v deg D k ( v ) 2 m k E T ∼ D c T ( δ T ( { v } ) ) deg D k ( v ) = 1 m k E ∑ f L T ( f ) ≥ k 9 . \sum_v\frac{\deg_{D_k}(v)}{2m_k}\frac{\mathbb{E}_{T\sim\mathcal{D}}c_T(\delta_T(\{v\}))}{\deg_{D_k}(v)}=\frac{1}{m_k}\mathbb{E}\sum_fL_T(f)\ge\frac{k}{9}. v ∑ 2 m k deg D k ( v ) deg D k ( v ) E T ∼ D c T ( δ T ({ v })) = m k 1 E f ∑ L T ( f ) ≥ 9 k .
At least one ratio is at least this average. This argument does not require bounded degrees. Since n k = ( 2 ⋅ 4 k + 4 ) / 3 n_k=(2\cdot4^k+4)/3 n k = ( 2 ⋅ 4 k + 4 ) /3 , we have k = Θ ( log n k ) k=\Theta(\log n_k) k = Θ ( log n k ) .
The witness edge or singleton vertex may depend on the distribution. The last two assertions use the prescribed fundamental-cut capacities on spanning trees; they do not assert a lower bound for arbitrary capacities on arbitrary abstract trees.
A sharp support tradeoff for fundamental-cut congestion
Let A A A and B B B be disjoint sets of k ≥ 2 k\ge2 k ≥ 2 vertices. Put a star H A H_A H A on A A A and a star H B H_B H B on B B B , assign capacity M = k 2 M=k^2 M = k 2 to every star edge, and add all k 2 k^2 k 2 edges between A A A and B B B with capacity one. This is the graph G k G_k G k of Theorem 1.2 ; it has n = 2 k n=2k n = 2 k vertices and m = k 2 + 2 k − 2 m=k^2+2k-2 m = k 2 + 2 k − 2 edges. For a cross edge a b ab ab , its canonical tree is
T a b = H A ∪ H B ∪ { a b } . T_{ab}=H_A\cup H_B\cup\{ab\}. T ab = H A ∪ H B ∪ { ab } .
Figure 3 illustrates the construction.
Proof of Theorem 1.2 . Every atom has small probability. Every spanning tree contains a cross edge. Every cut separating the endpoints of a cross edge either splits a heavy star, and therefore crosses an edge of capacity M M M , or is the cut A ∣ B A \mid B A ∣ B , whose capacity is also M M M . Thus each tree T T T has a cross edge f ∈ T f \in T f ∈ T with g T ( f ) = L T ( f ) / c f ≥ M g_T(f) = L_T(f)/c_f \ge M g T ( f ) = L T ( f ) / c f ≥ M . If the tree has probability q T q_T q T , nonnegativity gives
E g T ′ ( f ) ≥ q T g T ( f ) ≥ q T M . \mathbb{E}g_{T'}(f) \ge q_T g_T(f) \ge q_T M. E g T ′ ( f ) ≥ q T g T ( f ) ≥ q T M .
A congestion bound ρ \rho ρ therefore requires q T ≤ ρ / M q_T \le\rho/M q T ≤ ρ / M for every tree, and hence ∣ supp q ∣ ≥ ⌈ M / ρ ⌉ |\operatorname{supp} q| \ge\lceil M/\rho\rceil ∣ supp q ∣ ≥ ⌈ M / ρ ⌉ .
Congestion approaching one. Choose a canonical tree T a b T_{ab} T ab by choosing its bridge. Its bridge has congestion M M M . Let f f f be a heavy spoke with leaf v ∈ A v \in A v ∈ A , and let r v r_v r v be the probability that the bridge’s A A A endpoint is v v v . If the endpoint is elsewhere, the fundamental-cut capacity of f f f is M + k M+k M + k . If the endpoint is v v v , it is M + k ( k − 1 ) M+k(k-1) M + k ( k − 1 ) . Consequently
E g T ( f ) = 1 + 1 k + k − 2 k r v . (34) \mathbb{E}g_T(f) = 1 + \frac{1}{k} + \frac{k-2}{k}r_v. \tag*{(34)} E g T ( f ) = 1 + k 1 + k k − 2 r v . (34)
The identical formula holds on side B B B with the column marginal. The uniform distribution over all k 2 k^2 k 2 bridges has every row and column marginal equal to 1 / k 1/k 1/ k . Every cross edge then has expected congestion one, and every heavy spoke has expected congestion
β k = 1 + 2 ( k − 1 ) k 2 . \beta_k = 1 + \frac{2(k-1)}{k^2}. β k = 1 + k 2 2 ( k − 1 ) .
For any capacitated graph and every spanning tree,
∑ f c f g T ( f ) = ∑ e c e ∣ E ( P T ( e ) ) ∣ ≥ ∑ e c e . \sum_f c_f g_T(f) = \sum_e c_e |E(P_T(e))| \ge\sum_e c_e. f ∑ c f g T ( f ) = e ∑ c e ∣ E ( P T ( e )) ∣ ≥ e ∑ c e .
Taking expectations shows that maximum expected congestion is always at least one. The optimum on G k G_k G k is therefore between 1 1 1 and β k \beta_k β k , and tends to one.
Attaining the support lower bound. Fix ρ ≥ γ k \rho\ge\gamma_k ρ ≥ γ k and put r = ⌈ k 2 / ρ ⌉ r = \lceil k^2/\rho\rceil r = ⌈ k 2 / ρ ⌉ . Choose r r r distinct cross edges whose degrees on either side are at most ⌈ r / k ⌉ \lceil r/k \rceil ⌈ r / k ⌉ . To construct them, decompose K k , k K_{k,k} K k , k into its k k k cyclic perfect matchings, take ⌊ r / k ⌋ \lfloor r/k \rfloor ⌊ r / k ⌋ whole matchings, and then take the required number of edges from the next matching. Use the corresponding r r r canonical trees uniformly.
Every selected cross edge has expected congestion k 2 / r ≤ ρ k^2/r \le\rho k 2 / r ≤ ρ . Each row and column marginal is at most ⌈ r / k ⌉ / r ≤ 1 / k + 1 / r \lceil r/k \rceil/r \le1/k + 1/r ⌈ r / k ⌉ / r ≤ 1/ k + 1/ r . By (34) , every spoke has expected congestion at most
1 + 2 ( k − 1 ) k 2 + k − 2 k r ≤ 1 + 2 ( k − 1 ) k 2 + ρ k 2 ≤ ρ . 1 + \frac{2(k-1)}{k^2} + \frac{k-2}{kr} \le1 + \frac{2(k-1)}{k^2} + \frac{\rho}{k^2} \le\rho. 1 + k 2 2 ( k − 1 ) + k r k − 2 ≤ 1 + k 2 2 ( k − 1 ) + k 2 ρ ≤ ρ .
The first inequality uses r ≥ k 2 / ρ r \ge k^2/\rho r ≥ k 2 / ρ and ( k − 2 ) / k ≤ 1 (k-2)/k \le1 ( k − 2 ) / k ≤ 1 ; the last is equivalent to ρ ≥ γ k \rho\ge\gamma_k ρ ≥ γ k . Together with the atom bound, this proves the exact support tradeoff.
Cut and metric approximation with one tree. Fix any canonical tree and equip it with fundamental-cut capacities. Each heavy spoke has tree capacity at most 2 M 2M 2 M , and the bridge has tree capacity M M M . For a cut S S S , let h h h be the number of heavy spokes crossing it. If h = 0 h=0 h = 0 , the cut is trivial or is A ∣ B A \mid B A ∣ B , and its tree and graph capacities are equal. If h ≥ 1 h \ge1 h ≥ 1 , then
c T ( δ T ( S ) ) ≤ 2 M h + M ≤ 3 M h ≤ 3 c G ( δ G ( S ) ) . c_T(\delta_T(S)) \le2Mh + M \le3Mh \le3c_G(\delta_G(S)). c T ( δ T ( S )) ≤ 2 M h + M ≤ 3 M h ≤ 3 c G ( δ G ( S )) .
The reverse inequality follows from fundamental-cut domination. Thus one canonical tree gives factor-three approximation of every cut.
For the metric assertion, assign unit length to every edge and choose the bridge joining the two star centers. The resulting spanning tree has diameter three. Distinct vertices are at graph distance at least one, so its distance distortion is at most three.
Corollary 10.3 (Nearly m m m trees for near-unit congestion). For the family G k G_k G k , congestion at most β k = 1 + 2 ( k − 1 ) / k 2 \beta_k = 1 + 2(k-1)/k^2 β k = 1 + 2 ( k − 1 ) / k 2 is attainable, and every distribution with congestion at most β k \beta_k β k has support at least
⌈ k 4 m ⌉ ≥ m − 2 n + 4 = ( 1 − o ( 1 ) ) m . \left\lceil\frac{k^4}{m} \right\rceil\ge m - 2n + 4 = (1-o(1))m. ⌈ m k 4 ⌉ ≥ m − 2 n + 4 = ( 1 − o ( 1 )) m . The same lower bound applies to an optimal-congestion distribution.
Proof. Here m = k 2 + 2 k − 2 m = k^2 + 2k - 2 m = k 2 + 2 k − 2 and β k = m / k 2 \beta_k = m/k^2 β k = m / k 2 . Apply Theorem 1.2 with ρ = β k \rho= \beta_k ρ = β k and use
k 4 m = m − 4 k + 4 + ( 2 k − 2 ) 2 m ≥ m − 2 n + 4. \frac{k^4}{m} = m - 4k + 4 + \frac{(2k-2)^2}{m} \ge m - 2n + 4. m k 4 = m − 4 k + 4 + m ( 2 k − 2 ) 2 ≥ m − 2 n + 4.
Since m = Θ ( n 2 ) m = \Theta(n^2) m = Θ ( n 2 ) , the bound is ( 1 − o ( 1 ) ) m (1-o(1))m ( 1 − o ( 1 )) m . The optimum congestion is at most β k \beta_k β k , so every optimal distribution also satisfies the premise.
The support lower bound uses the guarantee on each original edge. It therefore distinguishes that guarantee from expected cut approximation. It also separates the two orders of maximum and expectation: every tree of G k G_k G k has congestion at least k 2 k^2 k 2 , whereas one fixed mixture has maximum expected congestion 1 + O ( 1 / k ) 1 + O(1/k) 1 + O ( 1/ k ) .
Why a single distribution need not achieve both guarantees
The preceding results construct a distance distribution for a fixed length vector and a capacity distribution for a fixed capacity vector. One cannot require the same distribution to achieve both logarithmic guarantees for arbitrary prescribed data. The next general family is the obstruction in Andersen–Feige [3 Section 3.4] .
Proposition 10.4 (A simultaneous stretch–congestion obstruction). For every integer k ≥ 2 k \ge2 k ≥ 2 there is a graph on N = k ( k − 1 ) + 2 N = k(k-1) + 2 N = k ( k − 1 ) + 2 vertices, with unit lengths and unit capacities, such that every distribution over its spanning trees has either maximum expected stretch or maximum expected congestion at least ( k + 1 ) / 2 = Ω ( N ) (k+1)/2 = \Omega(\sqrt{N}) ( k + 1 ) /2 = Ω ( N ) .
Proof. Take k k k internally vertex-disjoint s s s –t t t paths, each consisting of k k k edges, and add the direct edge e 0 = s t e_0=st e 0 = s t . Write q = Pr [ e 0 ∈ T ] q=\Pr[e_0\in T] q = Pr [ e 0 ∈ T ] . If e 0 ∉ T e_0\notin T e 0 ∈ / T , the unique tree path between s s s and t t t is one complete long path and has length k k k ; if e 0 ∈ T e_0\in T e 0 ∈ T it has length one. Thus
E s T ( e 0 ) = k − ( k − 1 ) q . \mathbb{E}_{sT}(e_0)=k-(k-1)q. E s T ( e 0 ) = k − ( k − 1 ) q .
If e 0 ∈ T e_0\in T e 0 ∈ T , exactly one edge is missing from each long path: two missing edges on a path would disconnect an internal segment, and no missing edge would leave a cycle with e 0 e_0 e 0 . Each of those k k k missing edges has its tree path traverse e 0 e_0 e 0 , in addition to the path of e 0 e_0 e 0 itself. Hence L T ( e 0 ) = k + 1 L_T(e_0)=k+1 L T ( e 0 ) = k + 1 in this case. If e 0 ∉ T e_0\notin T e 0 ∈ / T , its load as a graph edge is zero. Therefore
E g T ( e 0 ) = ( k + 1 ) q . \mathbb{E}_{gT}(e_0)=(k+1)q. E g T ( e 0 ) = ( k + 1 ) q .
When q ≤ 1 / 2 q\le1/2 q ≤ 1/2 the expected stretch is at least ( k + 1 ) / 2 (k+1)/2 ( k + 1 ) /2 ; when q ≥ 1 / 2 q\ge1/2 q ≥ 1/2 the expected congestion is at least ( k + 1 ) / 2 (k+1)/2 ( k + 1 ) /2 . This proves the assertion for every k k k and every distribution.
Choice of parameters
Fix
ε = 1 100 , δ = 1 10000 , χ = 8 , A 0 = 1024 , ζ = 2 − 10 . \varepsilon=\frac{1}{100},\qquad\delta=\frac{1}{10000},\qquad\chi=8,\qquad A_0=1024,\qquad\zeta=2^{-10}. ε = 100 1 , δ = 10000 1 , χ = 8 , A 0 = 1024 , ζ = 2 − 10 .
We verify the sufficient conditions in (14) : ε ≤ 1 / 8 \varepsilon\le1/8 ε ≤ 1/8 , δ ≤ 1 / 256 \delta\le1/256 δ ≤ 1/256 , χ ≥ 4 \chi\ge4 χ ≥ 4 , A 0 = 1024 ≥ 54 χ = 432 A_0=1024\ge54\chi=432 A 0 = 1024 ≥ 54 χ = 432 , and ζ = 2 − 10 \zeta=2^{-10} ζ = 2 − 10 . These are the only restrictions on the five parameters used in the entropy and radius analysis.
For clarity, their roles can also be checked directly. First, ε ≤ 1 / 8 \varepsilon\le1/8 ε ≤ 1/8 implies
2 ε < 1 , 2 1 − 2 ε < 3 , 1 − ε 1 + ε > 1 χ . 2\varepsilon<1,\qquad\frac{2}{1-2\varepsilon}<3,\qquad\frac{1-\varepsilon}{1+\varepsilon}>\frac{1}{\chi}. 2 ε < 1 , 1 − 2 ε 2 < 3 , 1 + ε 1 − ε > χ 1 .
These inequalities justify geometric summation of nonseed runs, the uniform source-scale comparison, and the ineligibility of far captures. Second, δ ≤ 1 / 256 \delta\le1/256 δ ≤ 1/256 gives
4 δ < 1 2 , 64 δ ≤ 1 4 , 128 δ ≤ 1 2 < 4 3 ε . 4\delta<\frac{1}{2},\qquad64\delta\le\frac{1}{4},\qquad128\delta\le\frac{1}{2}<\frac{4}{3\varepsilon}. 4 δ < 2 1 , 64 δ ≤ 4 1 , 128 δ ≤ 2 1 < 3 ε 4 .
Thus source cores lie in every allowed cell, low-entropy seeds are unflagged, and their inherited successor charges contradict the capture floor. Third, χ > 3 \chi>3 χ > 3 makes the far mass gap g χ = ( χ − 3 ) / ( 12 χ ) g_\chi=(\chi-3)/(12\chi) g χ = ( χ − 3 ) / ( 12 χ ) positive. The condition A 0 ≥ 54 χ A_0\ge54\chi A 0 ≥ 54 χ contains each cell’s near set in the outer ball of any represented near vertex. It also contains far source cores, since 3 ( 3 + δ ε / χ ) < 10 < A 0 3(3+\delta\varepsilon/\chi)<10<A_0 3 ( 3 + δ ε / χ ) < 10 < A 0 . Finally ζ ≤ 2 − 10 < log 2 \zeta\le2^{-10}<\log2 ζ ≤ 2 − 10 < log 2 gives K v ≤ 2 K_v\le2 K v ≤ 2 at a low-entropy seed, while the complement of its inner ball has mass strictly below 2 − 10 M 2^{-10}M 2 − 10 M . Source eligibility and integrality give a source core of mass at least 2 − 10 M 2^{-10}M 2 − 10 M , so the two cores intersect.
The resulting entropy stride and constants are
L = 1 + ⌈ log 2 A 0 δ ε ⌉ = 31 , C F = 714240 < 10 6 , C T = 238080 < 3 ⋅ 10 5 . L=1+\left\lceil\log_2\frac{A_0}{\delta\varepsilon}\right\rceil=31,\qquad C_F=714240<10^6,\qquad C_T=238080<3\cdot10^5. L = 1 + ⌈ log 2 δ ε A 0 ⌉ = 31 , C F = 714240 < 1 0 6 , C T = 238080 < 3 ⋅ 1 0 5 .
Indeed 2 29 < A 0 / ( δ ε ) = 1024000000 < 2 30 2^{29}<A_0/(\delta\varepsilon)=1024000000<2^{30} 2 29 < A 0 / ( δ ε ) = 1024000000 < 2 30 . Using log 2 < 3 / 4 \log2<3/4 log 2 < 3/4 , we have
C R < 8 , C S ≤ 2400 ( 1 + 768 ) < 2 ⋅ 10 6 , L + 1 + ε C F 4 χ = 255.2 < 256. C_R<8,\qquad C_S\le2400(1+768)<2\cdot10^6,\qquad L+1+\frac{\varepsilon C_F}{4\chi}=255.2<256. C R < 8 , C S ≤ 2400 ( 1 + 768 ) < 2 ⋅ 1 0 6 , L + 1 + 4 χ ε C F = 255.2 < 256.
Consequently
C M < 8 ( 10 6 + 12 ⋅ 2 ⋅ 10 6 ⋅ 256 ) < 10 11 . C_M < 8(10^{6} + 12 \cdot2 \cdot10^{6} \cdot256) < 10^{11}. C M < 8 ( 1 0 6 + 12 ⋅ 2 ⋅ 1 0 6 ⋅ 256 ) < 1 0 11 .
The local constants are determined by geometry: C g e o = 34 C_{\mathrm{geo}} = 34 C geo = 34 and W = 3 C g e o = 102 W = 3C_{\mathrm{geo}} = 102 W = 3 C geo = 102 . Using log 2 > 1 / 2 \log2 > 1/2 log 2 > 1/2 and the preceding bounds in Theorem 6.5 gives
C d e t = 34 ( 32 χ / ε + 1 ) log 2 + 3 [ 8 W ( C M + C F ) + 4 W C T + 216 χ W ε ] < 3 ⋅ 10 14 < 10 15 . C_{\mathrm{det}} = \frac{34(32\chi/\varepsilon+ 1)}{\log2} + 3\left[8W(C_M + C_F) + 4WC_T + \frac{216\chi W}{\varepsilon}\right] < 3 \cdot10^{14} < 10^{15}. C det = log 2 34 ( 32 χ / ε + 1 ) + 3 [ 8 W ( C M + C F ) + 4 W C T + ε 216 χ W ] < 3 ⋅ 1 0 14 < 1 0 15 .
We may therefore take C 0 = 10 15 C_0 = 10^{15} C 0 = 1 0 15 in Theorem 1.1 . The following rational integer-valued response bound is valid:
ρ n = 10 15 ⌈ log 2 ( 2 n ) ⌉ , ∑ e = u v a e d T ( u , v ) ≤ ρ n ∑ e a e ℓ e . \rho_n = 10^{15}\left\lceil\log_2(2n)\right\rceil,\qquad\sum_{e=uv} a_e d_T(u,v) \leq\rho_n \sum_e a_e\ell_e. ρ n = 1 0 15 ⌈ log 2 ( 2 n ) ⌉ , e = uv ∑ a e d T ( u , v ) ≤ ρ n e ∑ a e ℓ e .
The large displayed constant is used only to make the exact rational master fully explicit; no numerical optimization is needed for the O ( log n ) O(\log n) O ( log n ) guarantee.
AI Use Statement. This work resulted from extensive interaction between the author and OpenAI Astra, with the author providing suggestions and critically checking Astra’s responses through iterative exchanges. OpenAI Astra was also used in writing and revising the manuscript through detailed prompts from the author, with Gemini Argon assisting in some parts. The author directed the process and evaluated the proposed ideas and arguments. A repository link for Lean files formalizing the main theorems are available on Github at https://github.com/hajiaghayi/Spanning-Tree-Embedding/blob/8dbf5dc3abc88ae1bc7074ff58482eab1aea73e8/SpanningTreeDecomposition1.zip .