Geometric separator

Last updated

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.

Contents

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.

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

HyperplaneGeometricSeparatorCounterExample.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).

HyperplaneGeometricSeparatorCounterExample2.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:
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 Nd-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.

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 . 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 large squares.

Next, consider all the "small squares" – the squares whose side-length is less than .

For every t, define: intersect(t) as the set of small squares intersected by the boundary of Rt. For every t1 and t2, if , then . Therefore, there is a gap of at least between the boundary of Rt1 and the boundary of Rt2. Therefore, intersect(t1) and intersect(t2) are disjoint. Therefore:

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

The separator we look for is the rectangle Rt, where . [7]

Application example

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

Generalizations

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

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 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 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

Notes

  1. 1 2 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. JSTOR   24492346.
  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. JSTOR   24492022.
  3. Pach, János; Tardos, Gábor (2001-10-28). "Separating convex sets by straight lines". 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 . MR   2060636.. See discussion following Theorem 2.2 and Fig. 1(a).
  5. 1 2 3 4 5 Smith, W. D.; Wormald, N. C. (1998). "Geometric separator theorems and applications". Proceedings 39th Annual Symposium on Foundations of Computer Science (Cat. No. 98CB36280). p. 232. doi:10.1109/sfcs.1998.743449. ISBN   978-0-8186-9172-0. S2CID   17962961.
  6. 1 2 Chan, T. M. (2003). "Polynomial-time approximation schemes for packing and piercing fat objects". Journal of Algorithms. 46 (2): 178–189. CiteSeerX   10.1.1.21.5344 . 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. 1 2 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. Miller, Gary L.; Teng, Shang-Hua; Thurston, William; Vavasis, Stephen A. (1997). "Separators for sphere-packings and nearest neighbor graphs". J. ACM. 44 (1): 1–29. doi: 10.1145/256292.256294 . S2CID   17331739..
  11. Kyncl, Jan. "Simultaneous geometric separator". MathOverflow. Retrieved 4 February 2014.

Related Research Articles

The aspect ratio of a geometric shape is the ratio of its sizes in different dimensions. For example, the aspect ratio of a rectangle is the ratio of its longer side to its shorter side—the ratio of width to height, when the rectangle is oriented as a "landscape".

<span class="mw-page-title-main">Circle</span> Simple curve of Euclidean geometry

A circle is a shape consisting of all points in a plane that are at a given distance from a given point, the centre. The distance between any point of the circle and the centre is called the radius.

<span class="mw-page-title-main">Convex set</span> In geometry, set whose intersection with every line is a single line segment

In geometry, a subset of a Euclidean space, or more generally an affine space over the reals, is convex if, given any two points in the subset, the subset contains the whole line segment that joins them. Equivalently, a convex set or a convex region is a subset that intersects every line into a single line segment . For example, a solid cube is a convex set, but anything that is hollow or has an indent, for example, a crescent shape, is not convex.

A histogram is a visual representation of the distribution of quantitative data. The term was first introduced by Karl Pearson. To construct a histogram, the first step is to "bin" the range of values— divide the entire range of values into a series of intervals—and then count how many values fall into each interval. The bins are usually specified as consecutive, non-overlapping intervals of a variable. The bins (intervals) are adjacent and are typically of equal size.

A perimeter is a closed path that encompasses, surrounds, or outlines either a two dimensional shape or a one-dimensional length. The perimeter of a circle or an ellipse is called its circumference.

<span class="mw-page-title-main">Rectangle</span> Quadrilateral with four right angles

In Euclidean plane geometry, a rectangle is a quadrilateral with four right angles. It can also be defined as: an equiangular quadrilateral, since equiangular means that all of its angles are equal ; or a parallelogram containing a right angle. A rectangle with four sides of equal length is a square. The term "oblong" is used to refer to a non-square rectangle. A rectangle with vertices ABCD would be denoted as  ABCD.

<span class="mw-page-title-main">Packing problems</span> Problems which attempt to find the most efficient way to pack objects into containers

Packing problems are a class of optimization problems in mathematics that involve attempting to pack objects together into containers. The goal is to either pack a single container as densely as possible or pack all objects using as few containers as possible. Many of these problems can be related to real-life packaging, storage and transportation issues. Each packing problem has a dual covering problem, which asks how many of the same objects are required to completely cover every region of the container, where objects are allowed to overlap.

<span class="mw-page-title-main">24-cell</span> Regular object in four dimensional geometry

In four-dimensional geometry, the 24-cell is the convex regular 4-polytope (four-dimensional analogue of a Platonic solid) with Schläfli symbol {3,4,3}. It is also called C24, or the icositetrachoron, octaplex (short for "octahedral complex"), icosatetrahedroid, octacube, hyper-diamond or polyoctahedron, being constructed of octahedral cells.

<span class="mw-page-title-main">Semicircle</span> Geometric shape

In mathematics, a semicircle is a one-dimensional locus of points that forms half of a circle. It is a circular arc that measures 180°. It has only one line of symmetry.

<span class="mw-page-title-main">Cupola (geometry)</span> Solid made by joining an n- and 2n-gon with triangles and squares

In geometry, a cupola is a solid formed by joining two polygons, one with twice as many edges as the other, by an alternating band of isosceles triangles and rectangles. If the triangles are equilateral and the rectangles are squares, while the base and its opposite face are regular polygons, the triangular, square, and pentagonal cupolae all count among the Johnson solids, and can be formed by taking sections of the cuboctahedron, rhombicuboctahedron, and rhombicosidodecahedron, respectively.

<span class="mw-page-title-main">16-cell</span> Four-dimensional analog of the octahedron

In geometry, the 16-cell is the regular convex 4-polytope (four-dimensional analogue of a Platonic solid) with Schläfli symbol {3,3,4}. It is one of the six regular convex 4-polytopes first described by the Swiss mathematician Ludwig Schläfli in the mid-19th century. It is also called C16, hexadecachoron, or hexdecahedroid [sic?].

<span class="mw-page-title-main">120-cell</span> Four-dimensional analog of the dodecahedron

In geometry, the 120-cell is the convex regular 4-polytope (four-dimensional analogue of a Platonic solid) with Schläfli symbol {5,3,3}. It is also called a C120, dodecaplex (short for "dodecahedral complex"), hyperdodecahedron, polydodecahedron, hecatonicosachoron, dodecacontachoron and hecatonicosahedroid.

<i>k</i>-d tree Multidimensional search tree for points in k dimensional space

In computer science, a k-d tree is a space-partitioning data structure for organizing points in a k-dimensional space. K-dimensional is that which concerns exactly k orthogonal axes or a space of any number of dimensions. k-d trees are a useful data structure for several applications, such as:

<span class="mw-page-title-main">Guillotine cutting</span> Process of producing small rectangular items of fixed dimensions

Guillotine cutting is the process of producing small rectangular items of fixed dimensions from a given large rectangular sheet, using only guillotine-cuts. A guillotine-cut is a straight bisecting line going from one edge of an existing rectangle to the opposite edge, similarly to a paper guillotine.

<span class="mw-page-title-main">Hyperplane separation theorem</span> On the existence of hyperplanes separating disjoint convex sets

In geometry, the hyperplane separation theorem is a theorem about disjoint convex sets in n-dimensional Euclidean space. There are several rather similar versions. In one version of the theorem, if both these sets are closed and at least one of them is compact, then there is a hyperplane in between them and even two parallel hyperplanes in between them separated by a gap. In another version, if both disjoint convex sets are open, then there is a hyperplane in between them, but not necessarily any gap. An axis which is orthogonal to a separating hyperplane is a separating axis, because the orthogonal projections of the convex bodies onto the axis are disjoint.

Discrepancy of hypergraphs is an area of discrepancy theory.

<span class="mw-page-title-main">Point reflection</span> Geometric symmetry operation

In geometry, a point reflection is a transformation of affine space in which every point is reflected across a specific fixed point. When dealing with crystal structures and in the physical sciences the terms inversion symmetry, inversion center or centrosymmetric are more commonly used.

In graph theory, the planar separator theorem is a form of isoperimetric inequality for planar graphs, that states that any planar graph can be split into smaller pieces by removing a small number of vertices. Specifically, the removal of vertices from an n-vertex graph can partition the graph into disjoint subgraphs each of which has at most vertices.

In geometry, specifically projective geometry, a blocking set is a set of points in a projective plane that every line intersects and that does not contain an entire line. The concept can be generalized in several ways. Instead of talking about points and lines, one could deal with n-dimensional subspaces and m-dimensional subspaces, or even more generally, objects of type 1 and objects of type 2 when some concept of intersection makes sense for these objects. A second way to generalize would be to move into more abstract settings than projective geometry. One can define a blocking set of a hypergraph as a set that meets all edges of the hypergraph.

In computational geometry, a maximum disjoint set (MDS) is a largest set of non-overlapping geometric shapes selected from a given set of candidate shapes.