Source-linked AI summary

$c$-Packedness versus $\lambda$-Low-Density in Geometric Graphs: Theory and Practice

Gregor Diatzko, Félix Lasseux, Sabine Storandt

arXiv:2609.27231v1cs.DScs.CG

TL;DR

The paper asks whether c-packedness and λ-low-density support practical algorithms on large road networks, where parameter values and the usefulness of existing bounds were unclear. It develops parameter-computation, separator, and tree-decomposition methods; experiments find the parameters sufficiently small for parameterized bounds to matter, and yield an exact oracle with O(c) query time.

  • Problem

    The values of c and λ on large real-world networks were unknown, while some parameterized algorithms had impractical running times or separator bounds.

    Source: §1 · p. 2, §3 · p. 4
  • Method

    The paper develops scalable parameter-computation methods, improved balanced separators, and polynomial-time tree decompositions for c-packed graphs.

    Source: §1 · p. 2, §1.2 · p. 3, §4 · p. 7
  • Results

    Experiments find c and λ small enough for parameterized bounds to be meaningful, especially λ, and the exact oracle achieves O(c) query time.

    Source: §1 · p. 2, §6 · p. 17
  • Takeaways & Limitations

    The techniques produce small balanced separators and support concise exact distance oracles on large road networks while retaining theoretical guarantees.

    Source: §1 · p. 2, §1.2 · p. 3, §6 · p. 17

    The oracle's subsequent distance computations take O(c^3 n log n), and reducing this dependence on c remains open.

    Source: §6 · p. 18

Abstract

from arXiv · show

When designing algorithms for geometric graphs, exploiting structural parameters can lead to significantly improved bounds. Two prominent parameters in this context are $c$-packedness and $\lambda$-low density, both of which locally restrict graph complexity. Parameterized algorithms based on these parameters have been developed for computing well-separated pair decompositions, balanced separators, as well as distance oracles. Nevertheless the practical applicability of algorithms parameterized by $c$ or $\lambda$ remains unclear. While $c$-packed and $\lambda$-low-density graphs have been proposed as realistic models for road networks, the actual parameter values of large real-world instances have so far remained unknown, and existing theoretical guarantees are partially too loose for practical usage. In this paper we first devise scalable implementations for the approximate computation of $c$ and the exact computation of $\lambda$. Our experiments on road networks with millions of edges reveals a significant gap between the two parameters. On the theoretical side we prove that $c\in O(\lambda\sqrt n)$ which complements the known result that $\lambda\in O(c)$. Furthermore we present improved parameterized algorithms for balanced separator computation that reduce the separator size in theory and practice. We also show how to compute a tree decomposition with a width linear in the respective parameterized balanced separator size in polynomial time. This structural result yields a variety of new algorithmic consequences. Among them is an exact distance oracle with query time $O(c)$ for $c$-packed graphs after polynomial-time preprocessing, which improves upon the previous $O(c\log n)$ bound. Our experiments show that the proposed techniques efficiently produce small balanced separators and enable the construction of concise exact distance oracles on large road networks.

1 Introduction

The paper investigates whether geometric graph parameters support practical algorithms on large road networks, where their real-world values and practical applicability had been unclear. It measures these parameters and develops improved separator and tree-decomposition algorithms, including an exact distance oracle with O(c) query time.

1 Introduction

  • Road-network graph classes aim to explain small separators and support provable algorithms, but parameterized oracles can be impractical and real-world c and λ values were unknown.The paper defines practicality in terms of running time poly(p) · ˜O(n), noting that heuristic implementations may forfeit theoretical guarantees.Source: §1 · pp. 1–2

1.1 Related Work

  • Prior c and λ work offered parameterized algorithms, but reported experiments did not establish scalability to general graphs with millions of edges.Earlier c evaluations used small-to-moderate trajectory sets, while λ evaluations covered triangulations or smaller networks.Source: §1.1 · pp. 2–3

1.2 Contribution

  • The paper proves c ∈ O(λ√n), tightens parameter computation to exact λ and approximate c on large graphs, and finds λ consistently much smaller than c.The bound complements the known λ ∈ O(c) result.Source: §1.2 · p. 3
  • New techniques reduce c-packed separator sizes without increasing asymptotic running time, and the engineered algorithm consistently beats KaHIP on road networks.An improved PTAS for the smallest enclosing square problem is a key ingredient.Source: §1.2 · p. 3
  • A tree decomposition of width O(c) can be computed in O(c^2 n log n) time, enabling an exact distance oracle with O(c) query time after polynomial preprocessing.This improves on the previous O(c log n) query-time bound.Source: §1.2 · p. 3

2 Parameter Definitions and Relationships

The section defines low density and packedness through local geometric restrictions on edges. It relates the definitions to graph degree and outlines a scale-based argument for bounding edge length within a square.

2 Parameter Definitions and Relationships

  • A graph is λ-low-density when each radius-r square intersects at most λ edges longer than r; it is c-packed when clipped edge length in each such square is at most cr.Balls and axis-aligned squares give equivalent definitions up to a constant factor.Source: §2 · pp. 3–4
  • Both parameters upper-bound maximum degree, so graphs with constant λ or c have m ∈ O(n).The paper uses n for vertices and m for edges.Source: §2 · p. 4
  • The proof groups edge intersections with a unit-radius square by length scales, then accounts separately for edges shorter than the final scale.This scale decomposition covers all intersections and bounds their total length.Source: §2 · p. 4

3 Balanced Separators

The section targets balanced separators whose theoretical size bounds have hidden constants too large for practical use. It improves the construction’s balance and separator-size guarantees, including a square-based method with a 6c size bound for a 2/3-balanced separator.

3 Balanced Separators

  • An α-balanced separator leaves every connected component with at most αn nodes, and existing parameterized bounds have hidden constants too large for practical use.The section develops modified constructions to reduce separator sizes.Source: §3 · p. 4

3.1 From Ball Separators to Square Separators in Geometric Graphs

  • The c-packed separator algorithm also applies to λ-low-density graphs, achieving a size bound matching [32] without relying on randomization.The relationship between c and λ is used to transfer the bound.Source: §3.1 · p. 5
  • For the original iterative ball approach, reaching α = 1/2 requires about 475 iterations and yields a worst-case size bound of approximately 950c.Each iteration adds a separator, increasing the bound from 2c to 2kc.Source: §3.1 · p. 5
  • The refinements improve balance by adjusting the node counts inside and outside nested regions and using an approximation for the smallest enclosing ball.The ball approximation changes the radius relation and separator-size bound; switching to squares uses four half-side squares to cover a square and sets k = n/5.Source: §3.1 · p. 6
  • 6c bounds the size of a 2/3-balanced separator from the square-based method, over 60 times smaller than the original bound while retaining O(c^2n) running time.The bound uses ε = 0.05.Source: §3.1 · p. 7

4 Separator-Based Exact Distance Oracles

The section develops a polynomial-time route from balanced separators to tree decompositions and exact distance oracles for c-packed graphs. Its treewidth-based oracle achieves O(c) query time, independent of graph size.

4 Separator-Based Exact Distance Oracles

  • S-EDO stores separator-to-node distances and answers queries by checking bags along the separator tree, giving O(c log n) query time when separator sizes are O(c).Its space is O(cn log n), and its construction time is O(c^2n log^2 n).Source: §4 · p. 7

4.1 Tree Decompositions from Separators

  • A separator tree with bags expanded along paths to the root yields a valid tree decomposition, but its width grows to O(B log n).Reed’s algorithm avoids this logarithmic blow-up and returns width linear in B, though its balanced W-separator computation has exponential dependence on B.Source: §4.1 · pp. 7–8

4.2 Application to c-Packed Graphs

  • The adapted construction forms a balanced W-separator by uniting two separators, one based on all graph nodes and one based on W.The all-node separator is a 4/5-balanced square separator of size at most 3c; the resulting union has size at most 6c.Source: §4.2 · p. 9
  • Width 30c tree decompositions for c-packed graphs can be computed in O(c^2n log n) time.This enables tree-decomposition-based algorithms, including T-EDO, after polynomial-time preprocessing.Source: §4.2 · p. 9
  • O(c) query time is achieved by a c-packed exact distance oracle with O(c^3n log n) construction time and O(cn log n) space.Compared with S-EDO, T-EDO uses a factor of c more preprocessing time but improves query time by a factor of log n.Source: §4.2 · p. 9

5 Experimental Results

The experiments measure the density and packedness parameters on road networks, then evaluate engineered balanced separators and exact distance-oracle structures. They report a gap between c and λ and practical trade-offs in separator quality, runtime, and oracle query effort.

5.1 Efficient Computation of the Density

  • The λ algorithms count edges near points for each edge, using a quadtree to find nearby edges; the exact method computes an arrangement, while the 3-approximation samples points.Edges are processed by decreasing length or neighborhood size, and queries can run concurrently with limited synchronization.Source: §5.1 · pp. 10–11
  • The 3-approximation for λ is usually about 1.5 times worse than the optimum and is slower than the exact algorithm on the tested instances.The authors attribute the runtime difference to low instance densities; the multiple-core variant with limited synchronization is fastest in all cases.Source: §5.1 · p. 11, Figure 4 · p. 11

5.2 Efficient Computation of the Packedness

  • The 2-approximation sweeps outward from graph vertices, while the randomized 42-approximation uses a WSPD and translated-square quadtrees with adjusted constants for scalability.Kahan summation controls accumulated floating-point error in the 2-approximation; the adjusted 42-approximation’s observed factor remains below 5.Source: §5.2 · pp. 11–12
  • The 42-approximation is much faster than the 2-approximation, while a final sweep raises its packedness result to at least 75% of the 2-approximation’s on test instances.The 2-approximation usually finds a square twice as packed as the 42-approximation’s estimate; c is much larger than λ but much smaller than λ√n on tested networks.Source: §5.2 · p. 12, Figure 5 · p. 12

5.3 A Practical Separation Algorithm

  • The practical separator samples vertices, finds enclosing-ball pairs, and keeps the best separators; a precomputed spanning grid accelerates radius queries.A tested configuration used 100 sampled vertices and 50 separation runs, while a radius-ratio threshold of 1.2 bounds separator size by 6c.Source: §5.3 · pp. 13–14, Figure 6 · p. 13, Figure 8 · p. 15
  • On the small graph, the ball strategy yields better separators than KaHIP with fewer runs; on the large graph, it trades speed for quality against FAST or matches STRONG’s size faster.The comparison uses the same target balance, and the reported outcomes depend on the number of separation runs.Source: §5.3 · p. 15, Figure 10 · p. 16

5.4 Exact Distance Oracles

  • On the large dataset, both decompositions took around 4.5 hours to construct; the tree decomposition had 1.2 million bags and height 19.The separator decomposition had 1.7 million bags and height 24.Source: §5.4 · p. 15
  • T-EDO averages 20 distance comparisons per query, while S-EDO needs a bit less than four times more.The 75% quantile for T-EDO is 14, below S-EDO’s 25% quantile of 43.Source: §5.4 · p. 17, Figure 11 · p. 17

6 Conclusions and Future Work

The experiments find that λ and c remain moderate on large road networks, with λ consistently much smaller than c. The authors identify reducing the oracle’s dependence on c as an open problem.

6 Conclusions and Future Work

  • On large real-world road networks, λ and c are moderate and λ is consistently much smaller than c; an exact oracle with O(c) query time uses few distance computations on average.The maximum oracle time is also much smaller than the computed c-value.Source: §6 · p. 17
  • Har-Peled’s randomized O(c)-size separator algorithm may accelerate the implementation by avoiding max-flow, but its separator quality remains to be tested against this approach.Incorporating it into the tree-decomposition algorithm reduces that algorithm’s running time to O(cn log n).Source: §6 · p. 17
  • Oracle preprocessing remains dominated by distance computations taking O(c^3n log n), so reducing the dependence on c remains open.Source: §6 · p. 18

A Improved Approximation Algorithm for Smallest Enclosing Square

This appendix adapts a grid-based approximation method to find a smallest axis-aligned square enclosing k points. It returns a (1 + ε)-approximation in O(n + 1/ε) time when k ∈ Θ(n).

A Improved Approximation Algorithm for Smallest Enclosing Square

  • The square algorithm adapts the disk-based approach with minor modifications and improves its running-time dependence on ε.The authors state that the square version’s grid-cell point bound is 4k, compared with 5k for disks.Source: A Improved Approximation Algorithm for… · p. 22
  • The algorithm estimates the optimal radius, builds a grid, examines 3 × 3 cell clusters, and snaps cluster points to a finer grid before searching for the smallest qualifying square.For squares, at most 36k distinct points occur in a cluster before snapping; prefix sums support constant-time point counts during the search.Source: A Improved Approximation Algorithm for… · p. 22
  • For k ∈ Θ(n), the algorithm computes a (1 + ε)-approximation of the smallest axis-aligned k-enclosing square in O(n + 1/ε) time.Source: A Improved Approximation Algorithm for… · p. 22
Loading 2609.27231v1…