7 min read
On this page

Line Arrangements and Duality

Line Arrangements

An arrangement of lines A(L)\mathcal{A}(\mathcal{L}) is the subdivision of the plane induced by a set of nn lines L\mathcal{L}. The arrangement consists of:

  • Vertices: intersection points of lines. At most (n2)\binom{n}{2} vertices
  • Edges: maximal connected portions of lines between vertices. At most n2n^2 edges
  • Faces: connected regions of R2L\mathbb{R}^2 \setminus \bigcup \mathcal{L}. At most (n2)+n+1\binom{n}{2} + n + 1 faces

By Euler's formula for planar subdivisions: VE+F=1V - E + F = 1 (connected) or 22 (if counting the unbounded face appropriately).

Complexity: Θ(n2)\Theta(n^2) for nn lines in general position (no two parallel, no three concurrent). This is tight.

Construction

Incremental algorithm: insert lines one at a time, updating the arrangement.

For each new line \ell:

  1. Find a face of the current arrangement that \ell enters
  2. Walk along \ell through the arrangement, splitting each traversed face into two
  3. Each face traversal takes O(1)O(1) amortized time

Total complexity: O(n2)O(n^2) time and space, which is optimal since the arrangement has Θ(n2)\Theta(n^2) features. The walk uses the zone of the new line.

DCEL Representation

Arrangements are typically stored as a doubly-connected edge list (DCEL):

  • Each edge is split into two half-edges with opposite orientations
  • Each half-edge stores: origin vertex, twin half-edge, incident face, next/prev half-edge in face boundary
  • Supports O(1)O(1) local navigation (adjacent faces, edges around a vertex)

Zone Theorem

The zone of a line \ell in an arrangement of nn lines is the set of faces whose closure intersects \ell.

Zone theorem: the total complexity (number of edges bounding faces) of the zone of any line in an arrangement of nn lines is O(n)O(n).

Proof sketch: charge each edge in the zone to one of the nn lines. A careful argument shows each line contributes O(1)O(1) edges to the zone, giving O(n)O(n) total.

Consequences:

  • Inserting a new line into an arrangement takes O(n)O(n) time (traverse its zone)
  • Incremental construction: i=1nO(i)=O(n2)\sum_{i=1}^n O(i) = O(n^2), matching the output size
  • The zone theorem generalizes to hyperplane arrangements in Rd\mathbb{R}^d: zone complexity is O(nd1)O(n^{d-1})

Point-Line Duality

A fundamental transform mapping points to lines and vice versa. Standard duality:

Point p=(a,b)Line p:y=axb\text{Point } p = (a, b) \quad \longleftrightarrow \quad \text{Line } p^* : y = ax - b

Line :y=cxdPoint =(c,d)\text{Line } \ell : y = cx - d \quad \longleftrightarrow \quad \text{Point } \ell^* = (c, d)

Key property (incidence preservation): point pp lies on line \ell iff point \ell^* lies on line pp^*. More specifically:

  • pp is above \ell iff \ell^* is above pp^*
  • pp is below \ell iff \ell^* is below pp^*

Order preservation: the duality preserves the vertical ordering relation.

Primal-Dual Correspondences

| Primal | Dual | |---|---| | Point | Line | | Line | Point | | Point above line | Dual line above dual point | | Line through two points | Point = intersection of two dual lines | | Convex hull vertex | Line bounding dual arrangement face | | Upper envelope of lines | Convex hull of dual points |

Applications of Duality

  1. Ham-sandwich cuts via duality: transformed to finding a point in an arrangement with specific level properties
  2. Minimum area triangle: three points from a set SS forming the minimum area triangle. Via duality, this relates to closest lines in the dual arrangement
  3. Slope selection: find the kk-th smallest slope among all (n2)\binom{n}{2} point pairs. Via duality, this becomes a kk-level problem. O(nlogn)O(n \log n) algorithm via parametric search or Cole's method
  4. Linear programming: the dual of finding the lowest point above nn lines is finding the intersection point of lines below nn dual points

Levels in Arrangements

The kk-level of an arrangement of lines is the set of points that have exactly kk lines strictly below them. Equivalently, it is a piecewise-linear xx-monotone curve.

  • Level 0 = lower envelope; level nn = upper envelope
  • Complexity of the kk-level: at most O(nk+1)O(n \sqrt{k+1}) (Dey, 1998). Best lower bound: n2Ω(logn)n \cdot 2^{\Omega(\sqrt{\log n})} (Toth). Exact complexity is a major open problem
  • The k\leq k-level (all points with k\leq k lines below) has complexity O(nk)O(nk)

Upper and Lower Envelopes

The upper envelope of nn lines: the pointwise maximum maxi(aix+bi)\max_i (a_i x + b_i).

  • Complexity: O(n)O(n) segments (each line contributes at most one segment)
  • Computation: O(nlogn)O(n \log n) via divide-and-conquer or incremental insertion
  • Equivalent to the convex hull of the dual point set (upper hull)

For curves (not lines): the upper envelope of nn curves that pairwise intersect at most ss times has complexity λs(n)\lambda_s(n), a slowly growing function related to Davenport-Schinzel sequences:

  • s=1s = 1: λ1(n)=n\lambda_1(n) = n (lines)
  • s=2s = 2: λ2(n)=2n1\lambda_2(n) = 2n - 1 (line segments)
  • s=3s = 3: λ3(n)=Θ(nα(n))\lambda_3(n) = \Theta(n \alpha(n)) where α\alpha is the inverse Ackermann function

Ham-Sandwich Theorem

Theorem (Stone-Tukey): given dd finite Borel measures in Rd\mathbb{R}^d, there exists a hyperplane that simultaneously bisects all dd measures.

In R2\mathbb{R}^2: given two point sets (red and blue), there exists a line that bisects both. Proof uses the intermediate value theorem or Borsuk-Ulam theorem.

Algorithm in R2\mathbb{R}^2: O(nlogn)O(n \log n) via duality. Transform to finding a point at a specific level in the dual arrangement.

Discrete version: for nn red and nn blue points in general position, a ham-sandwich cut exists. Not unique in general; the space of valid cuts has specific structure.

Helly's Theorem

Theorem (Helly, 1923): let C={C1,,Cn}\mathcal{C} = \{C_1, \ldots, C_n\} be a finite collection of convex sets in Rd\mathbb{R}^d with n>dn > d. If every d+1d+1 of them have a common point, then all of them have a common point.

i=1nCi    every (d+1)-subset has non-empty intersection\bigcap_{i=1}^n C_i \neq \emptyset \iff \text{every } (d+1)\text{-subset has non-empty intersection}

In R2\mathbb{R}^2: if every three convex sets have a common point, then they all do.

Applications:

  • Halfplane intersection: nn halfplanes have a common point iff every 3 do. Reduces feasibility of LP to checking O(n3)O(n^3) triples (not efficient, but theoretically fundamental)
  • Helly dimension: generalizes to abstract convexity spaces
  • Fractional Helly theorem (Barany, 1982): if at least α(nd+1)\alpha \binom{n}{d+1} of the (d+1)(d+1)-subsets have non-empty intersection (for α>0\alpha > 0), then at least βn\beta n of the sets have a common point, where β=1(1α)1/(d+1)\beta = 1 - (1-\alpha)^{1/(d+1)}

Radon's Theorem

Theorem (Radon, 1921): any d+2d + 2 points in Rd\mathbb{R}^d can be partitioned into two sets whose convex hulls intersect.

In R2\mathbb{R}^2: any 4 points can be partitioned into two pairs whose convex hulls (line segments) intersect. Equivalently, one point is inside the triangle of the other three, or the two diagonals of the quadrilateral cross.

Radon partition: the partition S=S1S2S = S_1 \cup S_2 with conv(S1)conv(S2)\text{conv}(S_1) \cap \text{conv}(S_2) \neq \emptyset. Can be computed by finding a non-trivial affine dependence among the points.

Radon's theorem implies Helly's theorem (classical proof by induction).

Hyperplane Arrangements

In Rd\mathbb{R}^d, an arrangement of nn hyperplanes has:

  • At most (nd)\binom{n}{d} vertices
  • Total complexity O(nd)O(n^d) (faces of all dimensions)
  • Can be constructed incrementally in O(nd)O(n^d) time (optimal)

Number of faces: the number of regions (full-dimensional faces) is at most k=0d(nk)\sum_{k=0}^d \binom{n}{k}, which is Θ(nd)\Theta(n^d) in general position. This is the Zaslavsky formula (computable via the characteristic polynomial of the intersection lattice).

Applications:

  • Configuration spaces: robot arm reachability
  • Machine learning: hyperplane arrangements partition feature space into regions with different linear functions (ReLU networks, decision trees)
  • Discrete geometry: many counting problems reduce to arrangement complexity

Arrangement of Segments

For nn line segments with kk intersection points:

  • Output complexity: O(n+k)O(n + k)
  • Construction: Bentley-Ottmann sweep gives O((n+k)logn)O((n + k) \log n) (see Chapter 07)
  • Bounding kk: at most (n2)\binom{n}{2}, but often much smaller in practice

Applications and Connections

  • Linear programming: the optimal LP solution is a vertex of the arrangement of constraint hyperplanes. Dual LP corresponds to envelope computation
  • Geometric optimization: many problems (smallest enclosing circle, widest corridor) reduce to computations in line arrangements
  • Machine learning: SVM margins correspond to distances between arrangement features; ReLU network decision boundaries are hyperplane arrangements
  • Discrete and combinatorial geometry: Szemeredi-Trotter theorem bounds point-line incidences to O(n2/3m2/3+n+m)O(n^{2/3}m^{2/3} + n + m), with deep consequences for arrangement complexity and additive combinatorics