Geometric separator

From HandWiki
Revision as of 16:52, 6 February 2024 by TextAI2 (talk | contribs) (change)
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)

A geometric separator is a line (or another shape) that partitions a collection of geometric shapes into two subsets, such that proportion of shapes in each subset is bounded, and the number of shapes that do not belong to any subset (i.e. the shapes intersected by the separator itself) is small. When a geometric separator exists, it can be used for building divide-and-conquer algorithms for solving various problems in computational geometry.

Separators that are lines

General question

In 1979, Helge Tverberg[1] raised the following question. For two positive integers k, l, what is the smallest number n(k,l) such that, for any family of pairwise-disjoint convex objects in the plane, there exists a straight line that has at least k objects on one side and at least l on the other side?

The following results are known.

  • Obviously, n(1,1)=1.
  • Hope and Katchalski[2] proved that n(k,1) ≤ 12(k-1) for all k ≥ 2.
  • Villanger[1] proved that n(2,2) = ∞: he showed an infinite family of pairwise-disjoint segments such that no straight line has two segments in each side. Pach and Tardos[3] showed a simpler construction using only unit segments, and another construction using only discs (or squares).

Separators for axes-parallel rectangles

Given a set of N=4k disjoint axis-parallel rectangles in the plane, there is a line, either horizontal or vertical, such that at least N/4 rectangles lie entirely to each side of it (thus at most N/2 rectangles are intersected by the separator line).

Proof

Define W as the most western vertical line with at least N/4 rectangles entirely to its west. There are two cases:

  • If there are at least N/4 rectangles entirely to the east of W, then W is a vertical separator.
  • Otherwise, by moving W slightly to the west, we get a vertical line that intersects more than N/2 rectangles. Find a point on this line that has at least N/4 rectangles above and N/4 rectangles below it, and draw a horizontal separator through it.

Optimality

HyperplaneׂGeometricSeparatorCounterExample.svg

The number of intersected shapes, guaranteed by the above theorem, is O(N). This upper bound is asymptotically tight even when the shapes are squares, as illustrated in the figure to the right. This is in sharp contrast to the upper bound of O(N) intersected shapes, which is guaranteed when the separator is a closed shape (see previous section).

HyperplaneׂGeometricSeparatorCounterExample2.svg

Moreover, when the shapes are arbitrary rectangles, there are cases in which no line that separates more than a single rectangle can cross less than N/4 rectangles, as illustrated in the figure to the right.[4]

Generalizations

The above theorem can be generalized from disjoint rectangles to k-thick rectangles. Additionally, by induction on d, it is possible to generalize the above theorem to d dimensions and get the following theorem:[5]

Given N axis-parallel d-boxes whose interiors are k-thick, there exists an axis-parallel hyperplane such that at least:
[math]\displaystyle{ \lfloor(N+1-k)/(2d)\rfloor }[/math]
of the d-box interiors lie to each side of the hyperplane.

For the special case when k = N − 1 (i.e. each point is contained in at most N − 1 boxes), the following theorem holds:[5]

Given N axis-parallel d-boxes whose interiors are (N − 1)-thick, there exists an axis-parallel hyperplane that separates two of them.

The objects need not be boxes, and the separators need not be axis-parallel:

Let C be a collection of possible orientations of hyperplanes (i.e. C = {horizontal,vertical}). Given N d-objects, such that every two disjoint object are separated by a hyperplane with an orientation from C, whose interiors are k-thick, there exists a hyperplane with an orientation from C such that at least: (N + 1 − k)/O(C) of the d-objects interiors lie entirely to each side of the hyperplane.

Algorithmic versions

It is possible to find the hyperplanes guaranteed by the above theorems in O(Nd) steps. Also, if the 2d lists of the lower and upper endpoints of the intervals defining the boxes's ith coordinates are pre-sorted, then the best such hyperplane (according to a wide variety of optimality measures) may be found in O(Nd) steps.

Separators that are closed shapes

A simple case in which a separator is guaranteed to exist is the following:[5][6]

Given a set of n disjoint axis-parallel squares in the plane, there is a rectangle R such that, at most 2n/3 of the squares are inside R, at most 2n/3 of the squares are outside R, and at most O(sqrt(n)) of the squares are not inside and not outside R (i.e. intersect the boundary of R).

Thus, R is a geometric separator that separates the n squares into two subset ("inside R" and "outside R"), with a relatively small "loss" (the squares intersected by R are considered "lost" because they do not belong to any of the two subsets).

Proof

Define a 2-fat rectangle as an axis-parallel rectangle with an aspect ratio of at most 2.

Let R0 be a minimal-area 2-fat rectangle that contains the centers of at least n/3 squares. Thus every 2-fat rectangle smaller than R0 contains fewer than n/3 squares.

For every t in [0,1), let Rt be a 2-fat rectangle with the same center as R0, inflated by 1 + t.

  • Rt contains R0, so it contains the centers of at least n/3 squares.
  • Rt is less than twice as large as R0, so it can be covered by two 2-fat rectangles that are smaller than R0. Each of these 2-fat rectangles contains the centers of less than n/3 squares. Therefore Rt contains the centers of less than 2n/3 squares.

Now it remains to show that there is a t for which Rt intersects at most O(sqrt(n)) squares.

First, consider all the "large squares" – the squares whose side-length is at least [math]\displaystyle{ \operatorname{width}(R_0) / 2 \sqrt{n} }[/math]. For every t, the perimeter of Rt is at most 2·perimeter(R0) which is at most 6·width(R0), so it can intersect at most [math]\displaystyle{ 12 \sqrt{n} }[/math] large squares.

Next, consider all the "small squares" – the squares whose side-length is less than [math]\displaystyle{ \operatorname{width}(R_0) / 2 \sqrt{n} }[/math].

For every t, define: intersect(t) as the set of small squares intersected by the boundary of Rt. For every t1 and t2, if [math]\displaystyle{ |t_1-t_2| \geq 1/\sqrt{n} }[/math], then [math]\displaystyle{ |\operatorname{width}(R_{t_1})-\operatorname{width}(R_{t_2})| \geq \operatorname{width}(R_0)/\sqrt{n} }[/math]. Therefore, there is a gap of at least [math]\displaystyle{ \operatorname{width}(R_0) / 2 \sqrt{n} }[/math] between the boundary of Rt1 and the boundary of Rt2. Therefore, intersect(t1) and intersect(t2) are disjoint. Therefore:

[math]\displaystyle{ \sum_{j=0}^{\sqrt{n}-1}{|\operatorname{intersect}(j/\sqrt{n})|} \leq n }[/math]

Therefore, by the pigeonhole principle there is a certain j0 for which:

[math]\displaystyle{ |\operatorname{intersect}(j_0/\sqrt{n})| \leq \sqrt{n} }[/math]

The separator we look for is the rectangle Rt, where [math]\displaystyle{ t = j_0/\sqrt{n} }[/math].[7]

Application example

Using this separator theorem, we can solve certain problems in computational geometry in the following way:

  • Separate the input set of squares to two disjoint subsets;
  • Solve the problem on each subset separately;
  • Combine the solutions to the two sub-problems and get an approximate solution to the original problem.

Generalizations

The above theorem can be generalized in many different ways, with possibly different constants. For example:

  • Instead of squares, the input collection can contain arbitrary fat objects, such as: circles, rectangles with a bounded aspect ratio, etc.
  • Instead of two-dimensional shapes in a plane, the input collection can contain objects of any dimension, and they can be situated in a d-dimensional torus.
  • Instead of requiring that the shapes in the input collection be disjoint, we can put a weaker requirement, that the collection is:[5]
    • k-thick, i.e., each point is covered by at most k different shapes.
    • l-k-thick, i.e., each point is covered by at most k different shapes with a size ratio (size of largest shape divided by size of smallest shape) at most l.
    • k-overloaded, i.e., for any subcollection of shapes, the sum of their individual measures is at most k times the measure of their union.
  • Instead of a rectangle separator, the separator can be any shape that can be covered by smaller copies of itself.
  • Instead of bounding the number of shapes in each side of the separator, it is possible to bound any measure which satisfies certain axioms.[6]

Optimality

The ratio of 1:2, in the square separator theorem above, is the best that can be guaranteed: there are collections of shapes that cannot be separated in a better ratio using a separator that crosses only O(sqrt(n)) shapes. Here is one such collection (from theorem 34 of [5]):

Consider an equilateral triangle. At each of its 3 vertices, put N/3 shapes arranged in an exponential spiral, such that the diameter increases by a constant factor every turn of the spiral, and each shape touches its neighbours in the spiral ordering. For example, start with a 1-by-Φ rectangle, where Φ is the golden ratio. Add an adjacent Φ-by-Φ square and get another golden rectangle. Add an adjacent (1+Φ)-by-(1+Φ) square and get a larger golden rectangle, and so on.

Now, in order to separate more than 1/3 of the shapes, the separator must separate O(N) shapes from two different vertices. But to do this, the separator must intersect O(N) shapes.

Separators that are width-bounded strips between parallel hyperplanes

[8]

Let Q be a set of n points in the plane such that the minimal distance between points is d. Let a>0 be a constant.
There is a pair of parallel lines of distance a, such that at most 2n/3 points lie to each side of the strip, and at most [math]\displaystyle{ 1.3 {a \over d} \sqrt{n} }[/math] points lie inside the strip.
Equivalently: there is a line such that at most 2n/3 points lie to each side of it and at most [math]\displaystyle{ 1.3 {a \over d} \sqrt{n} }[/math] points lie at a distance of less than a/2 from it.

Proof sketch

Define the centerpoint of Q as a point o such that every line through it has at most 2n/3 points of Q in each side of it. The existence of a centerpoint can be proved using Helly's theorem.

For a given point p and constant a>0, define Pr(a,p,o) as the probability that a random line through o lies at a distance of less than a from p. The idea is to bound this probability and thus bound the expected number of points at a distance less than a from a random line through o. Then, by the pigeonhole principle, at least one line through o is the desired separator.

Applications

Bounded-width separators can be used for approximately solving the protein folding problem.[9] It can also be used for an exact sub-exponential algorithm to find a maximum independent set, as well as several related covering problems, in geometric graphs.[8]

Geometric separators and planar graph separators

The planar separator theorem may be proven by using the circle packing theorem to represent a planar graph as the contact graph of a system of disks in the plane, and then by finding a circle that forms a geometric separator for those disks.[10]

See also

  • Ham sandwich theorem: given n measurable objects in n-dimensional space, it is possible to divide all of them in half (with respect to their measure, i.e. volume) with a single (n − 1)-dimensional hyperplane.
  • Guillotine separation: the problem of separating convex objects in the plane using guillotine cuts.
  • Other Separation theorems.
  • Simultaneous separator: a separator that simultaneously separates the shapes in several collections, while simultaneously intersecting a small number of shapes in each collection, may not always exist.[11]

Notes

  1. 1.0 1.1 Tverberg, Helge (1979). "A separation property of plane convex sets". Mathematica Scandinavica 45 (2): 255–260. doi:10.7146/math.scand.a-11840. ISSN 0025-5521. 
  2. Hope, Rafael; Katchalski, Meir (1990). "Separating plane convex sets". Mathematica Scandinavica 66 (1): 44–46. doi:10.7146/math.scand.a-12291. ISSN 0025-5521. 
  3. Pach, János; Tardos, Gábor (2001-10-28). "Separating convex sets by straight lines" (in en). Discrete Mathematics 241 (1–3): 427–433. doi:10.1016/S0012-365X(01)00128-5. ISSN 0012-365X. 
  4. Dumitrescu, Adrian; Mitchell, Joseph S. B.; Sharir, Micha (2004). "Binary space partitions for axis-parallel segments, rectangles, and hyperrectangles". Discrete & Computational Geometry 31 (2): 207–227. doi:10.1007/s00454-003-0729-3. . See discussion following Theorem 2.2 and Fig. 1(a).
  5. 5.0 5.1 5.2 5.3 5.4 Smith, W. D.; Wormald, N. C. (1998). "Geometric separator theorems and applications". Proceedings 39th Annual Symposium on Foundations of Computer Science (Cat. No.98CB36280). pp. 232. doi:10.1109/sfcs.1998.743449. ISBN 978-0-8186-9172-0. http://users.monash.edu.au/~nwormald/papers/geomsep.ps.gz. 
  6. 6.0 6.1 Chan, T. M. (2003). "Polynomial-time approximation schemes for packing and piercing fat objects". Journal of Algorithms 46 (2): 178–189. doi:10.1016/s0196-6774(02)00294-8. 
  7. This proof is based on the more general proof of Chan (2003), but with the better constants of Smith&Wormald (1998).
  8. 8.0 8.1 Fu, B. (2011). "Theory and application of width bounded geometric separators". Journal of Computer and System Sciences 77 (2): 379–392. doi:10.1016/j.jcss.2010.05.003. 
  9. Fu, B.; Wang, W. (2007). "Geometric Separators and Their Applications to Protein Folding in the HP-Model". SIAM Journal on Computing 37 (4): 1014. doi:10.1137/s0097539704440727. 
  10. "Separators for sphere-packings and nearest neighbor graphs". J. ACM 44 (1): 1–29. 1997. doi:10.1145/256292.256294. .
  11. Kyncl, Jan. "Simultaneous geometric separator". MathOverflow. https://mathoverflow.net/q/155540.