Area-Optimal Simple Polygonalizations: The CG Challenge 2019
Abstract.
We give an overview of theoretical and practical aspects of finding a simple polygon of minimum (Min-Area) or maximum (Max-Area) possible area for a given set of points in the plane. Both problems are known to be -hard and were the subject of the 2019 Computational Geometry Challenge, which presented the quest of finding good solutions to more than 200 instances, ranging from all the way to .
1. Introduction
1.1. The Computational Geometry Challenge
The “CG:SHOP Challenge” (Computational Geometry: Solving Hard Optimization Problems) originated as a workshop at the 2019 Computational Geometry Week (CG Week) in Portland, Oregon in June, 2019. The goal was to conduct a computational challenge competition that focused attention on a specific hard geometric optimization problem, encouraging researchers to devise and implement solution methods that could be compared scientifically based on how well they performed on a database of instances. While much of computational geometry research has targeted theoretical research, often seeking provable approximation algorithms for -hard optimization problems, the goal of the CG Challenge was to set the metric of success based on computational results on a specific set of benchmark geometric instances. The 2019 CG Challenge focused on the problem of computing simple polygons whose vertices were a given set of points in the plane.
A tangible outcome is a new type of special issue, presented in this volume: a series of papers focusing on algorithm engineering methods for one difficult optimization problem. In this survey, we provide background and foundations of the underlying problem and give an overview of the results and contributions.
1.2. The 2019 Challenge Problem: Polygonizations of Optimal Area
The Traveling Salesman Problem (TSP) is one of the classic optimization problems: For a given set of locations and pairwise distances, find a shortest roundtrip that visits each position exactly once and returns to the start. In a geometric setting, in which locations correspond to a given set of points in the plane, and distances between points are induced by the Euclidean metric, it is a straighforward consequence of the triangle inequality that an optimal tour corresponds to a simple polygon with vertex set , such that has minimum total perimeter length.
This geometric motivation makes it natural to consider simple polygons with a given set of vertices that minimize another basic geometric measure: the enclosed area. This was considered in the past in the context of surface reconstruction, e.g., by O’Rourke (O’Rourke, 1980); as sketched in Section 2, there is also a close connection to point separation, which has gained importance in Artificial Intelligence. The same context has also raised interest in the corresponding maximization problem: For a given set of points in the plane, find a simple polygon with vertex set of maximum possible area. In the following, we refer to these two problems as Min-Area and Max-Area, respectively.
Problem: Min-Area and Max-Area
Given: A set of points in the plane.
Goal: A simple polygon with vertex set of minimum (Min-Area) or maximum (Max-Area) possible enclosed area.
There are a number of features that make these problems suitable for optimization challenges, based on notable similarities to and differences from the TSP. As shown by Fekete (Fekete, 1992; Fekete and Pulleyblank, 1993; Fekete, 2000), both Min-Area and Max-Area are -hard; however, while membership in for the Euclidean TSP is a long-standing, famous open problem (e.g., Problem #33 in The Open Problems Project (Demaine et al., 2001)), it is straightforward for area optimization, so issues of numerical stability and checking vailidity of solutions do not come into play. On the other hand, edges in a polygon with small area need not be short, as shown in Figure 1. This makes it challenging to restrict potential neighbors of a point in a good polygonization, increasing the difficulty of employing local search methods for efficient algorithms. This explains an apparent practical distinction to the TSP: While provably optimal solutions for TSP benchmark instances of considerable size have been known for a while, exact methods for area optimization appear more elusive.

1.3. Related Work
1.3.1. History and Background
The origins of the problem of finding area-optimal polygons can be traced back at least to the early days of Computational Geometry (O’Rourke (O’Rourke, 1980) in 1980). Resolving the complexity was posed by Suri in 1989 at CCCG, and restated in a different context by Mitchell (Mitchell, 1993) and Mitchell and Suri (Mitchell and Suri, 1995).
Studying the set of possible polygonizations is of great interest in various applications (Aichholzer et al., 2002; Zhu et al., 1996; Ackerman et al., 2009; Abellanas et al., 1993; Arkin et al., 2003). Auer and Held (Auer and Held, 1996) gave methods for generating random polygonizations. O’Rourke and Suri (O’Rourke et al., 2017) raised the question of estimating the number of polygonizations as a function of the number of points. At this point, this is known to lie between (Garcıa et al., 2000) and (Sharir et al., 2013).
1.3.2. Complexity
Fekete (Fekete, 2000, 1992) gave a proof of -completeness for both problems, based on a reduction from Hamiltonicity of Planar Cubic Directed Graphs, and generalized this proof to higher dimensions: For any fixed and , it is -hard to find the polyhedron with -dimensional faces of minimum total volume for given vertices in -dimensional space. He also showed that no polynomial-time approximation scheme (PTAS) exists for Min-Area and presented a -approximation algorithm for Max-Area (Fekete and Pulleyblank, 1993; Fekete, 1992). Moreover, he proved that the -hardness of the minimization problem also applies for higher dimensions. More specifically he showed that for given and it is -hard to find the minimal volume polyhedron with -dimensional faces for given vertices.
1.3.3. Heuristics
Recent work was mainly focused on finding new heuristics for both Min-Area and Max-Area. Taranilla et al. (Taranilla et al., 2011) proposed three different heuristics. Peethambaran et al. (Peethambaran et al., 2016, 2015) proposed randomized and greedy algorithms for Min-Area and -dimensional variants of both problems.
1.3.4. Other Challenges
The Open Problems Project (TOPP), maintained by Demaine, Mitchell and O’Rourke (Demaine et al., 2001), is a library of long-standing unsolved problems. On the more practical side, there have been different efforts, based on benchmark libraries, such as the TSPLIB (Reinelt, 1991). Since 1990, the DIMACS implementation challenges have addressed questions of determining realistic algorithm performance where worst-case analysis is overly pessimistic and probabilistic models are too unrealistic. Since 1994, the Graph Drawing (GD) community has held annual contests in conjunction with its annual symposium to monitor and challenge the current state of the graph-drawing technology and to stimulate new research directions for graph layout algorithm. More recently, a variety of implementation challenges have gained traction in the world of programming and optimization, but not yet in the field of Computational Geometry.
1.4. Outcomes
The contest generated a large number of contributions, both for Min-Area and Max-Area. In the aftermath, the top teams were invited to describe their methods in detailed papers, which form the substance of this special issue.
-
•
Julien Lepagnot, Laurent Moalic, Dominique Schmitt: Optimal area polygonization by triangulation and ray-tracing (Lepagnot et al., 2022)
-
•
Loïc Crombez, Guilherme D. da Fonseca, Yan Gerard: Greedy and Local Search Solutions to the Minimum and Maximum Are (Crombez et al., 2022)
-
•
Nir Goren, Efi Fogel, Dan Halperin: Area-optimal polygonization using simulated annealing (Goren et al., 2022)
-
•
Günther Eder, Martin Held, Steinpor Jasonarson, Philipp Mayer, Peter Palfrader: 2-Opt moves and flips for area-optimal polygonizations (Eder et al., 2022)
-
•
Natanael Ramos, Rai Caetan de Jesus, Pedro de Rezende, Cid de Souza, Fabio Luiz Usberti: Heuristics for area optimal polygonizations (Ramos et al., 2022)
In addition, there is one paper focusing on exact methods for computing provably optimal solutions.
-
•
Sándor P. Fekete, Andreas Haas, Phillip Keldenich, Michael Perk, Arne Schmidt: Computing area-optimal simple polygonization (Fekete et al., 2022)
2. Pick’s Theorem and Integrality
For a simple polygon with vertices, computing the Euclidean length of its perimeter involves evaluating a sum of square roots, for which membership in is a long-standing open problem (Demaine et al., 2001). This differs from computing the area of , which can be evaluated quite efficiently, e. g., see O’Rourke (O’Rourke, 1998). An elegant combinatorial answer is given by Pick’s theorem. (This also implies benign objective values, in particular for vertices whose coordinates are all even numbers, as chosen in the contest.)
Theorem 1 (Pick (Pick, 1899)).
Let be a simple polygon with integer vertices; let be the number of grid points contained in the interior of , and let be the number of grid points on the boundary of . Then

There are several elegant ways to prove Pick’s theorem, three of which can be found in (Funkenbusch, 1974; Coxeter, 1969; Gaskell et al., 1976). For a discussion of alternative approaches see the article by Niven and Zuckermann (Niven and Zuckerman, 1967). There are numerous generalizations to other than the orthogonal grid, e. g., by Ren and Reay (Ren and Reay, 1987); see Reeve (Reeve, 1957) for a generalization to higher dimensions.
Pick’s theorem also provides a bridge to issues of point separation: Maximizing the enclosed area amounts to finding a simple polygon that captures as many additional grid points as possible, while minimizing the area corresponds to excluding as many as possible. As the given points must lie on the boundary of , and only points within the convex hull of come into play, we get the following lower and upper bounds for the area.
Theorem 2.
Let be a set of points in the plane that all have integer coordinates. Let denote the number of points of the integer grid that are not contained in and strictly inside the convex hull, and let be the number of grid points not in that are on the boundary of the convex hull.
Then for any simple polygon on the vertex set , we have

3. Approximation
Using the area of the convex hull as an upper bound can also be used for approximating Max-Area.
Theorem 1 ((Fekete, 1992)).
Let be a set of points in the plane. We can determine a simple polygon on that has area larger than , where denotes the area of . This can be done in time .


Proof.
Let be a point on the convex hull of . In time , sort the points of by the slope of the lines , such that the neighbors of on the convex hull are the first and the last point, respectively. If there is a set of points for which the slope is the same, break the tie by ordering them in increasing distance from , except when those points have the smallest of all slopes, in which case we take them in order of decreasing distance from . Connecting the points in this order yields a simple polygon on .
If , we are done. Suppose this is not the case. Then the set has area at least . Now it is not hard to see that there is a simple polygon that contains , implying the claim. ∎


As shown in Figure 5, is a tight bound for this approach, even if all possible choices for are tested. Moreover, it was shown in (Fekete, 1992) that it is -complete to decide whether there is a simple polygon that contains strictly more than of the area of the convex hull;
At this time, no constant-factor approximation algorithm for Min-Area is known, hinting at a higher level of difficulty of the minimization problem. It may well be the case that no such approximation can be computed in polynomial time.
4. Contest and Outcomes
The 2019 Challenge was well received, with 28 teams from all over the world and a range of different scientific areas competing; participation was open to anyone. The contest itself was run through a dedicated server at TU Braunschweig, hosted at https://cgshop.ibr.cs.tu-bs.de/competition/cg-shop-2019/. It opened on February 28, 2019, and closed on May 31, 2019.
4.1. Instances
The contest started with a total of 247 benchmark instances, as follows. Each of these instances consisted of points in the plane with integer coordinates. For , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , there were six instances each. In addition, there is one instance of size .
The instances were of four different types; see Fekete et al. (Fekete et al., 2017) for a more detailed description of how to generate benchmark instances through illumination maps.
-
•
uniform: uniformly distributed at random from a square
-
•
edge: randomly generated according to the distribution of the rate of change (the “edges”) from various images
-
•
illumination: randomly generated according to the distribution of brightness of an image (such as an illumination map)
-
•
orthogonally collinear points: randomly generated on an integral grid, resulting in a large number of collinear points (similar to printed circuit boards and distorted blueprints). .
4.2. Evaluation
The comparison between different teams was based on an overall score. For each instance, this score is the ratio given by the achieved area divided by the area of the convex hull; thus, the score is a number between 0 and 1. The total score achieved by each team was the sum of all 247 individual instance scores. Feasibility of submitted solutions was checked at the time of upload; for instances without a feasible solution, a default score of 1 (for minimization) or 0 (for maximization) was used. For multiple submissions by the same team, only the best feasible solution submitted was used to compute the score. In case of ties, the tiebreaker was set to be the date/time a specific score was obtained. This turned out not to be necessary.
4.3. Results
In the end, the top 10 in the leaderboard looked as shown in Table 1; note that according to the scoring function, a lower score is better for Min-Area, while a higher score is better for Max-Area. The progress over time of each team’s score can be seen in Figure 6 (for Min-Area) and Figure 7 (for Max-Area).
Rk. | Team | Score (Min) | Score (Max) | # best | (unique) | # best | (unique) |
---|---|---|---|---|---|---|---|
Min | sols. | Max | sols. | ||||
1 | OMEGA/Mulhouse (FR) | 23.393 ( 9.47%) | 227.247 (92.00%) | 223 | (181) | 184 | (138) |
2 | lcrombez/Clermont (FR) | 25.751 (10.43%) | 226.691 (91.78%) | 66 | (23) | 109 | (4) |
3 | cgl@tau/Tel Aviv (IS) | 35.289 (14.29%) | 206.612 (83.65%) | 13 | (0) | 12 | (0) |
4 | CGA/Salzburg (AU) | 36.069 (14.60%) | 197.568 (79.99%) | 26 | (0) | 23 | (0) |
4 | UNICAMP/Campinas (BR) | 46.432 (18.80%) | 201.839 (81.72%) | 3 | (0) | 3 | (0) |
6 | mperk/Braunschweig (GE) | 68.431 (27.70%) | 191.483 (77.52%) | 24 | (0) | 23 | (0) |
6 | L’Aquila-Perugia (IT) | 57.373 (23.23%) | 179.752 (72.77%) | 19 | (0) | 18 | (0) |
8 | Stony Brook (US) | 85.179 (34.49%) | 162.031 (65.60%) | 1 | (0) | 1 | (0) |
9 | zhengdw (CA) | 89.437 (36.21%) | 154.723 (62.64%) | 0 | (0) | 1 | (0) |
10 | TGP/Eindhoven (NL) | 112.561 (45.57%) | 154.548 (62.57%) | 18 | (0) | 26 | (0) |
The top 5 finishers were invited for contributions to this special issue, as follows.
-
(1)
Team OMEGA/Mulhouse (France): Julien Lepagnot, Laurent Moalic, Dominique Schmitt (Lepagnot et al., 2022)
-
(2)
Team lcrombez/Clermont Auvergne(France): Loïc Crombez, Guilherme D. da Fonseca, Yan Gerard (Crombez et al., 2022)
-
(3)
Team cgl@tau/Tel Aviv (Israel): Nir Goren, Efi Fogel, Dan Halperin (Goren et al., 2022)
-
(4)
Team CGA/Salzburg (France): Günther Eder, Martin Held, Steinthor Jasonarson, Philipp Mayer, Peter Palfrader (Eder et al., 2022).
-
(5)
Team UNICAMP/Campinas (Brazil): Natanael Ramos, Rai Caetan de Jesus, Pedro de Rezende, Cid de Souza, Fabio Luiz Usberti (Ramos et al., 2022)
Figure 6 shows the development of total scores over time for Min-Area and all teams; it can be seen that OMEGA passed lcrombez shortly before the deadline, with cgl@tau just squeezing by CGA. A similar and even tighter outcome between OMEGA and lcrombez for Max-Area can be seen in Figure 7; for that problem, cgl@tau also placed third, but UNICAMP managed to beat out CGA for fourth place. Thus, the ranking was consistent for both Min-Area and Max-Area, except for Campinas doing better than Salzburg in Max-Area, so they both shared fourth place. All five teams engineered their solutions with the use of a variety of specific tools. Details of their methods and the engineering decisions they made are given in their respective papers.


Figure 8 shows the spread of results for all Min-Area instances; it can be seen that there is a strong deviation over all instance sizes, even for the small ones. This is surprising as one would expect simple heuristics to perform reasonably well for small instances; however, it appears that even for very small instances with points of Min-Area, more advanced ideas are necessary to achieve good results. (See the contribution by Fekete et al. (Fekete et al., 2022) for a detailed study of exact methods.) At around points, the mean score drops visibly; the best teams are able to obtain nearly the same score for all instances sizes above points. A similar overview for Max-Area is given in Figure 9. Here the deviation is very small for the small instances, showing that all teams where able to obtain reasonably good solutions for small instances. The deviation increases with the problem size until around points, after which it remains homogenous but smaller than for Min-Area. Like for Min-Area, the best teams obtained nearly equal scores for all instances sizes, while the mean score drops visibly after points.


5. Conclusions
The 2019 CG Challenge motivated a considerable number of teams to engage in intensive optimization studies. This success has not only led to practical progress on the problem of area optimization, but also turned the CG Challenge into a continuing feature of CGWeek, spawning considerable work through the 2020 Challenge problem (Minimum Convex Partition) and the 2021 problem (Coordinated Motion Planning). This promises to motivate further work on the involved problems, as well as other practical geometric optimization work. We are confident that this will further strengthen the bridges between optimization theory and practical algorithm engineering.
Acknowledgment
This work was supported by DFG project “Computational Geometry: Solving Hard Optimization Problems” (CG:SHOP), FE 407/21-1.
References
- (1)
- Abellanas et al. (1993) Manuel Abellanas, Jesús García, Gregorio Hernández Peñalver, Ferran Hurtado, Oriol Serra, and Jorge Urrutia. 1993. Updating polygonizations. In Computer Graphics Forum, Vol. 12. Wiley Online Library, 143–152. Issue 3.
- Ackerman et al. (2009) Eyal Ackerman, Oswin Aichholzer, and Balázs Keszegh. 2009. Improved upper bounds on the reflexivity of point sets. Computational Geometry 42, 3 (2009), 241–249.
- Aichholzer et al. (2002) Oswin Aichholzer, Franz Aurenhammer, and Hannes Krasser. 2002. Enumerating order types for small point sets with applications. Order 19, 3 (2002), 265–281.
- Arkin et al. (2003) Esther M. Arkin, Joseph S. B. Mitchell, Sándor P. Fekete, Ferran Hurtado, Marc Noy, Vera Sacristán, and Saurabh Sethia. 2003. On the reflexivity of point sets. In Discrete and Computational Geometry. Springer, 139–156.
- Auer and Held (1996) Thomas Auer and Martin Held. 1996. RPG-heuristics for the generation of random polygons. In Canadian Conference on Computational Geometry (CCCG). 38–44.
- Coxeter (1969) Harold Scott Macdonald Coxeter. 1969. Introduction to Geometry. Wiley, New York.
- Crombez et al. (2022) Loïc Crombez, Guilherme D. da Fonseca, and Yan Gerard. 2022. Greedy and Local Search Heuristics to Build Area-Optimal Polygons. Journal of Experimental Algorithmics (2022). This issue.
- Demaine et al. (2001) Erik D. Demaine, Joseph S. B. Mitchell, and Joseph O’Rourke. 2001. The Open Problems Project. , 573–582 pages. http://cs.smith.edu/~orourke/TOPP/.
- Eder et al. (2022) Günther Eder, Martin Held, Steinthor Jasonarson, Philipp Mayer, and Peter Palfrader. 2022. 2-Opt moves and flips for area-optimal polygonalizations. Journal of Experimental Algorithmics (2022). This issue.
- Fekete (1992) Sándor P. Fekete. 1992. Geometry and the Travelling Salesman Problem. Ph.D. Thesis. Department of Combinatorics and Optimization, University of Waterloo, Waterloo, ON, Canada.
- Fekete (2000) Sándor P. Fekete. 2000. On simple polygonizations with optimal area. Discrete & Computational Geometry 23, 1 (2000), 73–110.
- Fekete et al. (2017) Sándor P. Fekete, Andreas Haas, Michael Hemmer, Michael Hoffmann, Irina Kostitsyna, Dominik Krupke, Florian Maurer, Joseph S. B. Mitchell, Arne Schmidt, Christiane Schmidt, and Julian Troegel. 2017. Computing Nonsimple Polygons of Minimum Perimeter. Journal of Computational Geometry 8 (2017), 340–365. Issue 1.
- Fekete et al. (2022) Sándor P. Fekete, Andreas Haas, Phillip Keldenich, Michael Perk, and Arne Schmidt. 2022. Computing area-optimal simple polygonalization. Journal of Experimental Algorithmics (2022). This issue.
- Fekete and Pulleyblank (1993) Sándor P. Fekete and William R. Pulleyblank. 1993. Area optimization of simple polygons. In Symposium on Computational Geometry (SoCG). ACM, 173–182.
- Funkenbusch (1974) Walter W. Funkenbusch. 1974. From Euler’s formula to Pick’s formula using an edge theorem. The American Mathematical Monthly 81, 6 (1974), 647–648.
- Garcıa et al. (2000) Alfredo Garcıa, Marc Noy, and Javier Tejel. 2000. Lower bounds on the number of crossing-free subgraphs of KN. Computational Geometry 16, 4 (2000), 211–221.
- Gaskell et al. (1976) R. W. Gaskell, M. S. Klamkin, and P. Watson. 1976. Triangulations and Pick’s theorem. Mathematics Magazine 49, 1 (1976), 35–37.
- Goren et al. (2022) Nir Goren, Efi Fogel, and Dan Halperin. 2022. Area optimal polygonization using simulated annealing. Journal of Experimental Algorithmics (2022). This issue.
- Lepagnot et al. (2022) Julien Lepagnot, Laurent Moalic, and Dominique Schmitt. 2022. Optimal Area Polygonization by Triangulation and Visibility Search. Journal of Experimental Algorithmics (2022). This issue.
- Mitchell (1993) Joseph S. B. Mitchell. 1993. Approximation algorithms for geometric separation problems. Technical Report. SUNY Stony Brook.
- Mitchell and Suri (1995) Joseph S. B. Mitchell and Subhash Suri. 1995. Separation and approximation of polyhedral objects. Computational Geometry: Theory and Applications 5, 2 (1995), 95–114.
- Niven and Zuckerman (1967) Ivan Niven and Herbert S. Zuckerman. 1967. Lattice points and polygonal area. The American Mathematical Monthly 74, 10 (1967), 1195–1200.
- O’Rourke (1980) Joseph O’Rourke. 1980. Polyhedral object models from 3D points. Technical Report IFI-HH-M-77/80. Universität Hamburg.
- O’Rourke (1998) Joseph O’Rourke. 1998. Computational geometry in C. Cambridge university press.
- O’Rourke et al. (2017) Joseph O’Rourke, Subhash Suri, and Csaba D. Tóth. 2017. Polygons. In Handbook of discrete and computational geometry, Csaba D. Toth, Joseph O’Rourke, and Jacob E. Goodman (Eds.). CRC Press, Chapter 30, 787–810.
- Peethambaran et al. (2015) Jiju Peethambaran, Amal Dev Parakkat, and Ramanathan Muthuganapathy. 2015. A randomized approach to volume constrained polyhedronization problem. Journal of Computing and Information Science in Engineering 15, 1 (2015), 011009.
- Peethambaran et al. (2016) Jiju Peethambaran, Amal Dev Parakkat, and Ramanathan Muthuganapathy. 2016. An Empirical Study on Randomized Optimal Area Polygonization of Planar Point Sets. Journal of Experimental Algorithmics 21 (2016), 1–10.
- Pick (1899) Georg Pick. 1899. Geometrisches zur Zahlenlehre. Sitzungsberichte des Deutschen Naturwissenschaftlich-Medicinischen Vereines für Böhmen “Lotos” in Prag 19 (1899), 311–319.
- Ramos et al. (2022) Natanael Ramos, Rai Caetan de Jesus, Pedro de Rezende, Cid de Souza, and Fabio Luiz Usberti. 2022. Triangle-based heuristics for area optimal polygonizations. Journal of Experimental Algorithmics (2022). This issue.
- Reeve (1957) John E. Reeve. 1957. On the volume of lattice polyhedra. Proceedings of the London Mathematical Society 3, 1 (1957), 378–395.
- Reinelt (1991) Gerhard Reinelt. 1991. TSPlib – A Traveling Salesman Problem library. ORSA J. on Computing 3, 4 (1991), 376–384.
- Ren and Reay (1987) Ding Ren and John R. Reay. 1987. The boundary characteristic and Pick’s theorem in the Archimedean planar tilings. Journal of Combinatorial Theory, Series A 44, 1 (1987), 110–119.
- Sharir et al. (2013) Micha Sharir, Adam Sheffer, and Emo Welzl. 2013. Counting plane graphs: Perfect matchings, spanning cycles, and Kasteleyn’s technique. Journal of Combinatorial Theory, Series A 120, 4 (2013), 777–794.
- Taranilla et al. (2011) Maria Teresa Taranilla, Edilma Olinda Gagliardi, and Gregorio Hernández Peñalver. 2011. Approaching minimum area polygonization. In XVII Congreso Argentino de Ciencias de la Computacion. 31–40.
- Zhu et al. (1996) Chong Zhu, Gopalakrishnan Sundaram, Jack Snoeyink, and Joseph S. B. Mitchell. 1996. Generating random polygons with given vertices. Computational Geometry 6, 5 (1996), 277–290.