Heilbronn triangle problem

Short description: On point sets with no small-area triangles

Unsolved problem in mathematics: What is the asymptotic growth rate of the area of the smallest triangle determined by three out of [math]\displaystyle{ n }[/math] points in a square, when the points are chosen to maximize this area? (more unsolved problems in mathematics)

Six points in the unit square, with the smallest triangles (red) having area 1/8, the optimal area for this number of points. Other larger triangles are colored blue. These points are an affine transformation of a regular hexagon, but for larger numbers of points the optimal solution does not form a convex polygon.

In discrete geometry and discrepancy theory, the Heilbronn triangle problem is a problem of placing points in the plane, avoiding triangles of small area. It is named after Hans Heilbronn, who conjectured that, no matter how points are placed in a given area, the smallest triangle area will be at most inversely proportional to the square of the number of points. His conjecture was proven false, but the asymptotic growth rate of the minimum triangle area remains unknown.

Definition

The Heilbronn triangle problem concerns the placement of [math]\displaystyle{ n }[/math] points within a shape in the plane, such as the unit square or the unit disk, for a given number [math]\displaystyle{ n }[/math]. Each triple of points form the three vertices of a triangle, and among these triangles, the problem concerns the smallest triangle, as measured by area. Different placements of points will have different smallest triangles, and the problem asks: how should [math]\displaystyle{ n }[/math] points be placed to maximize the area of the smallest triangle?

More formally, the shape may be assumed to be a compact set [math]\displaystyle{ D }[/math] in the plane, meaning that it stays within a bounded distance from the origin and that points are allowed to be placed on its boundary. In most work on this problem, [math]\displaystyle{ D }[/math] is additionally a convex set of nonzero area. When three of the placed points lie on a line, they are considered as forming a degenerate triangle whose area is defined to be zero, so placements that maximize the smallest triangle will not have collinear triples of points. The assumption that the shape is compact implies that there exists an optimal placement of [math]\displaystyle{ n }[/math] points, rather than only a sequence of placements approaching optimality. The number [math]\displaystyle{ \Delta_D(n) }[/math] may be defined as the area of the smallest triangle in this optimal placement. An example is shown in the figure, with six points in a unit square. These six points form [math]\displaystyle{ \tbinom63=20 }[/math] different triangles, four of which are shaded in the figure. Six of these 20 triangles, with two of the shaded shapes, have area 1/8; the remaining 14 triangles have larger areas. This is the optimal placement of six points in a unit square: all other placements form at least one triangle with area 1/8 or smaller. Therefore, [math]\displaystyle{ \Delta_D(6)=\tfrac18 }[/math].

Although researchers have studied the value of [math]\displaystyle{ \Delta_D(n) }[/math] for specific shapes and specific small numbers of points, Heilbronn was concerned instead about its asymptotic behavior: if the shape [math]\displaystyle{ D }[/math] is held fixed, but [math]\displaystyle{ n }[/math] varies, how does the area of the smallest triangle vary with [math]\displaystyle{ n }[/math]? That is, Heilbronn's question concerns the growth rate of [math]\displaystyle{ \Delta_D(n) }[/math], as a function of [math]\displaystyle{ n }[/math]. For any two shapes [math]\displaystyle{ D }[/math] and [math]\displaystyle{ D' }[/math], the numbers [math]\displaystyle{ \Delta_D(n) }[/math] and [math]\displaystyle{ \Delta_{D'}(n) }[/math] differ only by a constant factor, as any placement of [math]\displaystyle{ n }[/math] points within [math]\displaystyle{ D }[/math] can be scaled by an affine transformation to fit within [math]\displaystyle{ D' }[/math], changing the minimum triangle area only by a constant. Therefore, in bounds on the growth rate of [math]\displaystyle{ \Delta_D(n) }[/math] that omit the constant of proportionality of that growth, the choice of [math]\displaystyle{ D }[/math] is irrelevant and the subscript may be omitted.

Heilbronn's conjecture and its disproof

Heilbronn conjectured prior to 1951 that the minimum triangle area always shrinks rapidly as a function of [math]\displaystyle{ n }[/math]—more specifically, inversely proportional to the square of [math]\displaystyle{ n }[/math]. In terms of big O notation, this can be expressed as the bound [math]\displaystyle{ \Delta(n)=O\left(\frac{1}{n^2}\right). }[/math]

Solutions to the no-three-in-line problem, large sets of grid points with no three collinear points, can be scaled into a unit square with minimum triangle area [math]\displaystyle{ \Omega(1/n^2) }[/math].

In the other direction, Paul Erdős found examples of point sets with minimum triangle area proportional to [math]\displaystyle{ 1/n^2 }[/math], demonstrating that, if true, Heilbronn's conjectured bound could not be strengthened. Erdős formulated the no-three-in-line problem, on large sets of grid points with no three in a line, to describe these examples. As Erdős observed, when [math]\displaystyle{ n }[/math] is a prime number, the set of [math]\displaystyle{ n }[/math] points [math]\displaystyle{ (i,i^2\bmod n) }[/math] on an [math]\displaystyle{ n\times n }[/math] integer grid (for [math]\displaystyle{ 0\le i\lt n }[/math]) have no three collinear points, and therefore by Pick's formula each of the triangles they form has area at least [math]\displaystyle{ \tfrac12 }[/math]. When these grid points are scaled to fit within a unit square, their smallest triangle area is proportional to [math]\displaystyle{ 1/n^2 }[/math], matching Heilbronn's conjectured upper bound. If [math]\displaystyle{ n }[/math] is not prime, then a similar construction using a prime number close to [math]\displaystyle{ n }[/math] achieves the same asymptotic lower bound.

(Komlós Pintz) eventually disproved Heilbronn's conjecture, by using the probabilistic method to find sets of points whose smallest triangle area is larger than the ones found by Erdős. Their construction involves the following steps:

• Randomly place [math]\displaystyle{ n^{1+\varepsilon} }[/math] points in the unit square, for some [math]\displaystyle{ \varepsilon\gt 0 }[/math].
• Remove all pairs of points that are unexpectedly close together.
• Prove that there are few remaining low-area triangles and therefore only a sublinear number of cycles formed by two, three, or four low-area triangles. Remove all points belonging to these cycles.
• Apply a triangle removal lemma for 3-uniform hypergraphs of high girth to show that, with high probability, the remaining points include a subset of [math]\displaystyle{ n }[/math] points that do not form any small-area triangles.

The area resulting from their construction grows asymptotically as [math]\displaystyle{ \Delta(n)=\Omega\left(\frac{\log n}{n^2}\right). }[/math] The proof can be derandomized, leading to a polynomial-time algorithm for constructing placements with this triangle area.

Upper bounds

Every set of [math]\displaystyle{ n }[/math] points in the unit square forms a triangle of area at most inversely proportional to [math]\displaystyle{ n }[/math]. One way to see this is to triangulate the convex hull of the given point set [math]\displaystyle{ S }[/math], and choose the smallest of the triangles in the triangulation. Another is to sort the points by their [math]\displaystyle{ x }[/math]-coordinates, and to choose the three consecutive points in this ordering whose [math]\displaystyle{ x }[/math]-coordinates are the closest together. In the first paper published on the Heilbronn triangle problem, in 1951, Klaus Roth proved a stronger upper bound on [math]\displaystyle{ \Delta(n) }[/math], of the form [math]\displaystyle{ \Delta(n)=O\left(\frac{1}{n\sqrt{\log\log n}}\right). }[/math] The best bound known to date is of the form [math]\displaystyle{ \Delta(n)\leq\frac{\exp{\left(c\sqrt{\log n}\right)}}{n^{8/7}}, }[/math] for some constant [math]\displaystyle{ c }[/math], proven by (Komlós Pintz).

A new upper bound equal to [math]\displaystyle{ n^{-\frac{8}{7}-\frac{1}{2000}} }[/math] was proven by (Cohen Pohoata).

Specific shapes and numbers

(Goldberg 1972) has investigated the optimal arrangements of [math]\displaystyle{ n }[/math] points in a square, for [math]\displaystyle{ n }[/math] up to 16. Goldberg's constructions for up to six points lie on the boundary of the square, and are placed to form an affine transformation of the vertices of a regular polygon. For larger values of [math]\displaystyle{ n }[/math], (Comellas Yebra) improved Goldberg's bounds, and for these values the solutions include points interior to the square. These constructions have been proven optimal for up to seven points. The proof used a computer search to subdivide the configuration space of possible arrangements of the points into 226 different subproblems, and used nonlinear programming techniques to show that in 225 of those cases, the best arrangement was not as good as the known bound. In the remaining case, including the eventual optimal solution, its optimality was proven using symbolic computation techniques.

The following are the best known solutions for 7–12 points in a unit square, found through simulated annealing; the arrangement for seven points is known to be optimal.

• 

  7 points in a square, all 8 minimal triangles shaded ([math]\displaystyle{ A\approx 0.0839 }[/math])

• 

  8 points in a square, 5 of 12 minimal triangles shaded ([math]\displaystyle{ A\approx 0.0724 }[/math])

• 

  9 points in a square, 6 of 11 minimal triangles shaded ([math]\displaystyle{ A\approx 0.0549 }[/math])

• 

  10 points in a square, 3 of 16 minimal triangles shaded ([math]\displaystyle{ A\approx 0.0465 }[/math])

• 

  11 points in a square, 8 of 28 minimal triangles shaded ([math]\displaystyle{ A\approx 0.0370 }[/math])

• 

  12 points in a square, 3 of 20 minimal triangles shaded ([math]\displaystyle{ A\approx 0.0326 }[/math])

Instead of looking for optimal placements for a given shape, one may look for an optimal shape for a given number of points. Among convex shapes [math]\displaystyle{ D }[/math] with area one, the regular hexagon is the one that maximizes [math]\displaystyle{ \Delta_D(6) }[/math]; for this shape, [math]\displaystyle{ \Delta_D(6)=\tfrac16 }[/math], with six points optimally placed at the hexagon vertices. The convex shapes of unit area that maximize [math]\displaystyle{ \Delta_D(7) }[/math] have [math]\displaystyle{ \Delta_D(7)=\tfrac19 }[/math].

Variations

There have been many variations of this problem including the case of a uniformly random set of points, for which arguments based on either Kolmogorov complexity or Poisson approximation show that the expected value of the minimum area is inversely proportional to the cube of the number of points. Variations involving the volume of higher-dimensional simplices have also been studied.

Rather than considering simplices, another higher-dimensional version adds another parameter [math]\displaystyle{ k }[/math], and asks for placements of [math]\displaystyle{ n }[/math] points in the unit hypercube that maximize the minimum volume of the convex hull of any subset of [math]\displaystyle{ k }[/math] points. For [math]\displaystyle{ k=d+1 }[/math] these subsets form simplices but for larger values of [math]\displaystyle{ k }[/math], relative to [math]\displaystyle{ d }[/math], they can form more complicated shapes. When [math]\displaystyle{ k }[/math] is sufficiently large relative to [math]\displaystyle{ \log n }[/math], randomly placed point sets have minimum [math]\displaystyle{ k }[/math]-point convex hull volume [math]\displaystyle{ \Omega(k/n) }[/math]. No better bound is possible; any placement has [math]\displaystyle{ k }[/math] points with volume [math]\displaystyle{ O(k/n) }[/math], obtained by choosing some [math]\displaystyle{ k }[/math] consecutive points in coordinate order. This result has applications in range searching data structures.

See also

• Danzer set, a set of points that avoids empty triangles of large area

External links

• Weisstein, Eric W.. "Heilbronn Triangle Problem". http://mathworld.wolfram.com/HeilbronnTriangleProblem.html. 
• Erich's Packing Center, by Erich Friedman, including the best known solutions to the Heilbronn problem for small values of [math]\displaystyle{ n }[/math] for squares, circles, equilateral triangles, and convex regions of variable shape but fixed area