Arranging Pairwise Disjoint Shapes to Partition Point Sets
Abstract
We consider the fine-grained complexity of covering a set of points in the Euclidean plane using a fixed set of geometric objects corresponding to rigid-body translations and, where permitted, rotations of a specified shape . Under the Exponential Time Hypothesis (ETH), and both with and without a pairwise disjointness constraint, we establish that no -time algorithm can exist for this problem in the following cases: (case 1) translatable unit disks; (case 2) translatable fixed-area axis-aligned squares; or (case 3) translatable and rotatable fixed-area equilateral triangles. Furthermore, by way of establishing the #P-completeness under parsimonious reductions of positive -in--SAT with a cubic planar -connected clause-variable incidence graph – pertinent to hardness reductions for counting tilings {(Moore & Robson; Discrete Comput. Geom. 26(4); 2001), (Pak & Yang; J. Comb. Theory. Ser. A 120(7); 2013)} – we establish in each case that there exists a quadratic time reduction from #SAT to counting the possible coverage-induced partitions of . Finally, we consider constraints on the density of the points in that make our coverage problems tractable. In particular, letting be any (not necessarily connected) subregion of a radius disk characterized by a semi-algebraic function, and letting be a set of points, we consider the density requirements that: (constraint 1) every points have a minimum bounding disk of radius greater than ; or (constraint 2) any points have a minimum bounding disk of radius at least . Here, when is part of the input, under both (constraint 1) and (constraint 2), and with and without a pairwise disjointness requirement, we show that finding a minimum cardinality set of translatable and/or rotatable instances of covering all points in is fixed-parameter tractable in the size of the semi-algebraic description of .
Keywords and phrases:
geometric covering, geometric packing, clustering, ply, bounded ply, planar geometry, frequency assignment problem, Exponential Time Hypothesis (ETH), Counting Exponential Time Hypothesis (#ETH)Funding:
Robert D. Barish: JSPS Kakenhi grant 25K21149.Copyright and License:
2012 ACM Subject Classification:
Mathematics of computing Discrete mathematics ; Theory of computation Problems, reductions and completenessEditor:
Pierre FraigniaudSeries and Publisher:
Leibniz International Proceedings in Informatics, Schloss Dagstuhl – Leibniz-Zentrum für Informatik
1 Introduction
One of the more elegant existence theorems in plane geometry is a circa 2008 result of Inaba [34, 35] that any arrangement of points in can be covered by pairwise interior-disjoint unit disks. This can be seen by noting that unit disks pack with an optimal density of in a hexagonal lattice arrangement [46], by using the union bound (i.e., Boole’s inequality) to determine the probability that points in the plane are covered by a randomly rotated and translated hexagonal lattice packing of unit disks is at least , and finally by observing that if (see, e.g., [39]). In the years since, this bound has been shown to hold for [2], and the upper bound for the number of points in that cannot always be covered by pairwise interior-disjoint unit disks, originally set at due to a triangular lattice construction of Winkler [49], has fallen from [21], to [39], and most recently to [2]. For many of these proofs, the authors either implicitly or explicitly assume the unit disks to be open.
Inaba’s result [34, 35] immediately raises the question of the computational complexity of deciding whether pairwise interior-disjoint (i.e., non-overlapping) or pairwise disjoint unit disks can cover some specified set of points in the plane. It also beckons us beyond recreational mathematics into much deeper waters, where we encounter the frequency assignment problem [1, 15, 20, 29] – a well-known applied combinatorial optimization challenge, in which the objective is to assign a scarce and expensive set of well-separated frequencies (often the subject of government-led auctions [13, 32]) to a set of radio stations in order to avoid interference between any two signals. Here, we can interpret the points we are trying to cover with unit disks in Inaba’s puzzle as the intended recipients of a signal emitted somewhere along the radio frequency spectrum (e.g., Hz for AM radio to Hz for modern Wi-Fi), and the disks used to cover the points as the broadcast ranges of radio frequency emitters with omnidirectional antennas. In this context, a coverage of the point set with pairwise disjoint unit disks becomes an explicit scheme for solving the special case of the frequency assignment problem where we have a finite number of points to cover and only one available broadcast frequency.
While we are unable to directly solve Inaba’s puzzle generalized to points, in this work we attempt to make progress towards its resolution by establishing the hardness of covering a specified set of points in with a fixed number of individually translatable and possibly rotatable geometric objects. More specifically, we concern ourselves with the problem of deciding the existence of and counting the following types of partitions:
Definition 1.
2D Geometric Cover Partition (2D-GCP). Let be a finite set of points in , and let be a finite collection of geometric objects embedded in the plane. A 2D Geometric Cover Partition (2D-GCP) induced by is a partitioning of into subsets (where if ) under the constraint that each subset is the set of all points covered by the same geometric object .
Definition 2.
Pairwise Disjoint 2D Geometric Cover Partition (PD-2D-GCP). An instance of a 2D-GCP where all geometric objects in are required to be pairwise disjoint (i.e., having no points in common even along their boundaries).
We clarify that, per Definition 1, each point in must belong to the closed set of points for exactly one geometric object in . Furthermore, when we refer to a 2D-GCP induced by a set of geometric shapes , we clarify that we are specifically and only referring to the partition of a point set into subsets for some .
We remark that the task of finding an inducing set for a 2D-GCP or PD-2D-GCP consisting of rigid-body isometries of a common geometric object is a natural generalization of the Minimum Size 1-Ply Set Cover of Unit Squares (MSP-SC-US) problem, very recently introduced and proven to be NP-hard by Govindarajan et al. [27], where the objective is to find a minimum set of axis-aligned unit squares covering a set of points under the constraint that the geometric system has ply-at-most-. In particular, we can determine that the NP-hardness proof for MSP-SC-US given by Govindarajan et al. [27] still holds under a pairwise disjointness requirement. Extensive precedent for 2D-GCPs and PD-2D-GCPs also comes from the broader literature on geometric packing and coverage problems [3, 4, 9, 6, 8, 11, 18, 19, 22, 23, 24, 26, 25, 28, 31, 48]. Although a full review of this literature is outside the scope of the current work due to space constraints, we briefly note an interesting connection to point set labeling and cartography problems [23, 25, 48], where, for instance, the objective may be to place non-overlapping geometric markers on a map that uniquely label individual features of interest. Finally, we remark that, in contrast to the search or decision problems for 2D-GCPs and PD-2D-GCPs, there appears to be little precedent for the problem of counting these types of partitions.
2 Paper overview
After providing some brief preliminaries and clarifications in Section 3, we lay a foundation for our fine-grained complexity analysis in Section 4 (see Theorem 1 and Corollary 2) by giving a quadratic time parsimonious reduction from #SAT to #CPP--Conn--in--SAT (see Section 3.1 concerning parsimonious reductions and Section 3.4 concerning #CPP--Conn--in--SAT). We briefly remark that this result is of independent interest in tiling theory, as it strengthens an earlier finding of Pak and Yang [40, 41] by establishing the #P-completeness under parsimonious reductions of counting tilings of finite simply connected regions using at most types of rectangles having integral edge lengths (see Corollary 3). This concomitantly establishes the completeness of the tiling problem for Ueda & Nagao’s Another Solution Problem (ASP) complexity class [42, 47, 50], implying that the problem of finding a second tiling given an initial one is as difficult as finding a second solution for SAT given an initial satisfying assignment. We briefly remark that our proof of the #P-completeness of #CPP--Conn--in--SAT under parsimonious reductions also answers an open question in the literature (see, e.g., [30]) concerning the hardness of -in--SAT instances with planar -connected clause-variable incidence graphs.
Subsequently, letting be a set of points embedded in , we consider the fine-grained complexity of finding and counting 2D-GCPs or PD-2D-GCPs in cases where the set of geometric objects used to cover corresponds to: (case 1) unit disks that can be individually translated as desired; (case 2) fixed-area axis-aligned squares that can be individually translated as desired; and (case 3) fixed-area equilateral triangles that can be individually translated and rotated as desired. In particular, via Theorem 2 and Corollary 3, under the ETH and #ETH, we rule out the existence of a -time algorithm in each of these cases for finding and counting 2D-GCPs or PD-2D-GCPs.
Finally, in Section 5 we investigate density constraints on that allow for tractable instances of the more general problem of covering with a minimum cardinality set of translatable and/or rotatable instances of a semi-algebraic subset of a disk. Specifically, we consider two types of density constraints for : (constraint 1) that any points have a minimum bounding disk of radius greater than , or (constraint 2) that any points have a minimum bounding disk of radius at least (where we earlier establish a bound in Theorem 2 for when (constraint 2) can lead to a polynomial-time algorithm under the ETH). Here, letting be a subregion of a radius disk (where we can allow for translations and/or rotations of ) admitting a description as a semi-algebraic set, and letting be the cardinality of , we show that this coverage problem – both with and without a pairwise disjointness constraint – is fixed-parameter tractable as a function of the size of the description of .
3 Preliminaries and clarifications
3.1 Complexity theoretic notions and terminology
Letting (where we usually have that ) be a pair of integer counting problems, and letting be a string specifying an instance of , a many-one counting reduction from to is a pair of polynomial time computable functions and such that . Here, a parsimonious reduction is a many-one counting reduction where is simply a projection onto its second argument. Recall that a problem is fixed-parameter tractable with respect to a parameter , if for some computable function , letting be a string specifying a problem instance, we have that the time complexity for the problem can be expressed as .
Unless otherwise stated, we assume the standard word RAM model, where geometric primitives are encoded as bounded-precision rationals (or exact algebraic extensions where necessitated, e.g., for equilateral triangles) and rotations are restricted to exact algebraic extensions of the rationals. In Section 5, we additionally allow symbolic algebraic computation over the reals for intermediate algebraic operations.
3.2 Graph theoretic notions and terminology
All graphs in this work should be assumed to be undirected and simple (meaning free of loops and multi-edges). While we refer the reader to Bondy & Murty [10] or Diestel [17] concerning basic graph theoretic terminology and notions, briefly recall that a geometric intersection graph for a set of geometric objects embedded in the plane constructs a vertex set having a bijection with , and connects a pair of vertices if and only if their associated geometric objects in have a point in common.
3.3 Plane geometry notions and terminology
Recall that a unit disk is a disk of radius embedded in the plane, that a unit line segment is a line segment of length embedded in the plane, and that a fixed-area axis-aligned square is a square embedded in the plane in such a manner that its vertex coordinates are given by the points for some and some specified . When we say that a geometric object covers a point in the plane, unless otherwise specified, we mean that belongs to the closed set of points defining the geometric object. In the context of this work, all motions of objects should be assumed to be rigid-body transformations which preserve the Euclidean distances between all pairs of points belonging to the object, and due to our focus on achiral shapes, we concern ourselves only with rigid transformations corresponding to translations and rotations. When we say that geometric objects are pairwise disjoint, we mean that no two objects have a point in common, regardless of whether any such common point falls internal to or along the boundary of one or both objects. Here, the ply [37] of a collection of geometric objects corresponds to the maximum number of objects sharing a common totally interior point.
3.4 Cubic Planar Positive -Connected -in- Satisfiability (CPP--Conn--in--SAT)
Letting be a Conjunctive Normal Form (CNF) formula (i.e., a conjunction of disjunctions) having exactly three literals per clause (where any literal can be either positive or negative), recall that the -in--SAT decision problem for input is the task of deciding whether admits a satisfying assignment in which each clause contains exactly one “True” literal and exactly two “False” literals. Recall also that the clause-variable incidence graph for is the bipartite graph , where (respectively, ) is in bijection with the variables (respectively, clauses) for , and where we add an edge between a vertex and if and only if the variable in corresponding to occurs as a (positive or negative) literal in the clause in corresponding to . We can now define the Cubic Planar Positive -Connected -in- Satisfiability (CPP--Conn--in--SAT) problem as a special case of -in--SAT in which the clause-variable incidence graph for any input formula is restricted to being a cubic (i.e., -regular) planar -connected graph, and where all literals in are positive. We also define #CPP--Conn--in--SAT as the #P problem of counting all satisfying assignments for instances of CPP--Conn--in--SAT. Finally, we define Cubic Planar Positive -in- Satisfiability (CPP--in--SAT) (respectively, #CPP--in--SAT) as a relaxation of CPP--Conn--in--SAT (respectively, #CPP--Conn--in--SAT) where we do not require the clause-variable incidence graph to be -connected.
3.5 Exponential Time Hypothesis (ETH) and Counting Exponential Time Hypothesis (#ETH)
The Exponential Time Hypothesis (ETH) of Impagliazzo & Paturi [33] is the conjecture that an variable instance of -SAT cannot be solved in time. In this context, using the Sparsification Lemma [33], we can assume that has at most clauses. The Counting Exponential Time Hypothesis (#ETH) of Dell et al. [16] is a relaxation of the ETH where we only assume that no -time algorithm exists for the -SAT problem of counting all solutions for a given -SAT formula . Here, the Sparsification Lemma can be adapted [16] to again allow us to assume that has at most clauses.
3.6 Orthogonal integer lattice embeddings
Letting be a graph of maximum vertex degree having vertex set and edge set , an orthogonal integer lattice embedding of is an embedding (or “drawing”) of on a grid graph . To construct , we map each vertex to a distinct vertex in , where we can then say that this lattice point in hosts or is hosting . Next, for every edge , we create a non-self-intersecting polyline composed of end-to-end adjacent horizontal and vertical unit line segments falling along the edges of and corresponding to a simple path between the lattice points hosting vertices (note that the polyline must not otherwise intersect a vertex in hosting a vertex in ). Here, the meeting of any horizontal and vertical unit line segments common to the same polyline is referred to as a bend in the embedding. As a clarification, polyline-polyline crossings are permitted, though we may refer to as a planar orthogonal integer lattice embedding if no such crossings occur. In this context, we define a cell in as a square centered on a point with vertex coordinates . Note that when we identify a cell at a coordinate with some set of points (e.g., as we do with various gadgets in the Theorem 2 proof argument), each point should be assumed to be translated to the coordinate . Finally, for some , we define a -fold enlargement of as the embedding generated by moving all lattice points hosting either vertices of an embedded graph or polyline segment endpoints to the coordinate , then reconnecting polylines by adding horizontal and vertical unit line segments without the introduction of additional bends.
For the purposes of visual intuition, in Fig. 1(b) we give an explicit example of an orthogonal integer lattice embedding of the -cube from Fig. 1(a) on a grid. Here, embedded vertices from Fig. 1(a) retain their original colorations in Fig. 1(b), (thick black) edges are used to indicate polylines, and a cell is indicated using a (thick dashed purple) box.
3.7 Algebraic and semi-algebraic sets
An algebraic set is a subset of an affine space that can be represented as the locus of zeroes (i.e., zero-set) of a finite collection of polynomial equations. A semi-algebraic set is a subset of an affine space represented as the union, intersection, and/or complement of finite collections of polynomial equations and inequalities. In this work, we will only consider semi-algebraic sets in .
4 Covering points in the plane with pairwise disjoint geometric objects
In this section, we first lay a foundation for our fine-grained analysis by giving a parsimonious reduction from #SAT to #CPP--Conn--in--SAT via Theorem 1 and Corollary 2. In the following Corollary 3, we briefly address an implication of Theorem 1 in tiling theory.
Theorem 1.
There exists an -time parsimonious reduction from any Planar -in--SAT formula with a -connected clause-variable incidence graph and clauses to a CPP--Conn--in--SAT formula.
Proof.
We proceed via a gadget-based reduction, where letting be an clause instance of Planar -in--SAT with a -connected clause-variable incidence graph , we will construct the clause-variable incidence graph for an clause instance of CPP--Conn--in--SAT. Specifically, the gadgets shown in Fig. 2(a–d) will be connected together via vertex identification operations between what we refer to as pole vertex pairs, corresponding to gadget vertices with labels of the form for some . Here, for some , we identify pole vertex pairs and by identifying the degree vertex with the degree vertex , and then identifying the degree vertex with the degree vertex . We clarify that a pole vertex pair will be referred to as unused if it has not yet participated in such a vertex identification operation. We now build via the following steps:
-
(step 1) For every degree variable vertex in , we construct an instance of the “-Output CPP--Conn--in--SAT variable gadget” shown in Fig. 2(a).
-
(step 2) For every clause vertex in , we construct an instance of the “CPP--Conn--in--SAT clause gadget” shown in Fig. 2(b).
-
(step 3) For every edge corresponding to a negated literal in , we construct an instance of the “CPP--Conn--in--SAT inverter gadget” shown in Fig. 2(c).
-
(step 4) For every degree or variable vertex in , we construct instances of the “-Output CPP--Conn--in--SAT variable gadget” shown in Fig. 2(a), then – as illustrated in Fig. 2(d) – we identify an unused pole vertex pair on each gadget with an unused pole vertex pair on the gadget , yielding a “-Output CPP--Conn--in--SAT variable gadget” with unused pole vertex pairs.
-
(step 5) For every edge between a variable vertex and a clause vertex in , where the variable corresponding to occurs as a positive literal in the clause corresponding to , we identify an unused pole vertex pair on (or if has degree or degree ) with an unused pole vertex pair on .
-
(step 6) For every edge between a variable vertex and a clause vertex in , where the variable corresponding to occurs as a negative literal in the clause corresponding to , we identify unused pole vertex pairs on the gadgets (or if has degree or degree ) and with unused pole vertex pairs on the gadget .
For the resulting graph , each pole vertex pair will correspond to a set of variable vertices in that, due to the properties of the gadgets in Fig. 2(a–d), will necessarily have the same “True” or “False” assignment in any witness for , and thus serve as a means of specifying inputs to gadgets that simulate clauses and variables in the input instance of Planar -in--SAT. Accordingly, to see that the -in--SAT formula corresponding to the clause-variable incidence graph represents a formula logically equivalent to , it suffices to make the following three observations concerning the -in--SAT formulas induced by the gadgets used to construct : (obs. 1) there are exactly two feasible satisfying assignments for each formula induced by a “-Output CPP--Conn--in--SAT variable gadget” or “-Output CPP--Conn--in--SAT variable gadget”, where one assignment sets all pole vertex pairs to “False” and the other sets all pole vertex pairs to “True” (simulating a variable vertex in ); (obs. 2) there are exactly three feasible satisfying assignments for each formula induced by a “CPP--Conn--in--SAT clause gadget”, where each assignment sets exactly one distinct pole vertex pair to “True” and all other pole vertex pairs to “False” (simulating a clause vertex in ); and (obs. 3) there are exactly two feasible satisfying assignments for formulas induced by the “CPP--Conn--in--SAT inverter gadget”, where each assignment sets exactly one of the two distinct pole vertex pairs to “True” (simulating an edge corresponding to a negated literal vertex in ).
It remains to observe that will necessarily be a cubic planar -connected graph corresponding to the clause-variable incidence graph of a CPP--Conn--in--SAT formula, and that (step 1) through (step 6) take time.
Corollary 2.
There exists an -time parsimonious reduction from any -SAT formula with clauses to a CPP--Conn--in--SAT formula with clauses.
Proof.
Observe first that we can parsimoniously reduce any clause instance of -SAT to an -clause instance of -SAT with a -connected clause-variable incidence graph in time by: (1) constructing two instances and of ; (2) for every variable in corresponding to a variable in , constructing the formula ; and (3) parsimoniously reducing the formula to a -SAT formula via standard methods. Next, recall that the reduction of Lichtenstein [36] can be used to parsimoniously reduce the clause instance of -SAT to an clause instance of Planar -SAT in time, and moreover, ensure that the clause-variable incidence graph of is -connected if the clause-variable incidence graph of is -connected. Thus, by Theorem 1, it now suffices to show the existence of an -time reduction from to an instance of Planar -in--SAT with a -connected clause-variable incidence graph. To proceed, let be a -in--SAT clause function which accepts a set of literals and returns “True” if and only if exactly one literal in the set is “True”, and observe that the -SAT clause is logically equivalent to the -in--SAT formula . Observe now that, after applying this reduction to an input Planar -SAT formula , the clause-variable incidence graph for the resulting Planar -in--SAT formula will be planar and -connected if the clause-variable incidence graph for is planar and -connected.
Corollary 3.
The following problems are -complete under parsimonious reductions: (P1) #CPP--Conn--in--SAT; (P2) counting tilings of finite simply connected regions using at most types of rectangles with integral edge lengths.
Proof.
Observe that (P1) and (P2) are trivially in #P. The claim for (P1) follows directly from Theorem 1 and Corollary 2. The claim for (P2) follows directly from a parsimonious reduction from #CPP--in--SAT to counting the stated type and number of rectangular tilings given in “Theorem 1.1” of Pak and Yang [40, 41]. We now proceed to establish Theorem 2 and Corollary 3.
Theorem 2.
Letting be the number of variables for an instance of #CPP--Conn--in--SAT, and letting be a set of points embedded in , there exists an -time parsimonious reduction from to counting either 2D-GCPs (or equivalently in the contexts we consider, PD-2D-GCPs) induced by a fixed number , where , geometric objects corresponding to:
-
(case 1) unit disks that can be individually translated;
-
(case 2) fixed-area axis-aligned squares that can be individually translated;
-
(case 3) fixed-area equilateral triangles that can be individually translated and rotated;
-
(case 4) disks of radius under the additional constraint that any points in have a minimum bounding disk of radius at least .
Proof.
We proceed by first treating (case 1) through (case 3), and then treating (case 4) by establishing its equivalence to (case 1) for the point sets we consider.
To begin, we consider an instance of #CPP--Conn--in--SAT with variables, clauses, and a cubic planar bipartite clause-variable incidence graph , where and are the partite sets of vertices in corresponding to clauses and variables, respectively. Using the algorithm of Biedl & Kant [7], we compute an orthogonal integer lattice embedding of having total area in time, and let be a -fold enlargement of (see Section 3.6 for an elaboration).
Next, we construct a set of points from the cells in , such that both 2D-GCPs and PD-2D-GCPs for in (case 1) through (case 3) will be in bijection with the satisfying assignments being counted by . Specifically, we begin by generating a set of points via the following steps:
-
(step 1) we identify every cell hosting a vertex with the “-in--SAT clause gadget” having point set {, , , , , , , , , , , , } – corresponding to the graph shown in Fig. 3(a.1) – rotated to ensure that points on the cell boundary fall along all outgoing polylines.
-
(step 2) We identify every cell hosting a vertex with the “variable gadget” having point set { , , , , , , , , , , , } – corresponding to the graph shown in Fig. 3(b.1) – rotated to ensure that points on the cell boundary fall along all outgoing polylines.
-
(step 3) For every cell not hosting a vertex and incident to two horizontal or two vertical polylines, we identify the cell with the “horizontal / vertical connector gadget” having point set {, , , , , , , , } – corresponding to the graph shown in Fig. 3(c.1) – rotated to ensure that points on the cell boundary fall along all outgoing polylines.
-
(step 4) For every cell not hosting a vertex and incident to a horizontal and a vertical polyline, we identify the cell with the “bend connector gadget” having point set {, , , , , , , } – corresponding to the graph shown in Fig. 3(d.1) – rotated to ensure that points on the cell boundary fall along all outgoing polylines.
-
(step 5) For every collection of cells identified with a “horizontal / vertical connector gadget” or “bend connector gadget” and corresponding to a common polyline in , to ensure that we replace at most of the point “horizontal / vertical connector gadgets” (which will necessarily exist since is a -fold enlargement of ) with the point “-shift horizontal / vertical connector gadget” having point set {, , , , , , , , , } – corresponding to the graph shown in Fig. 3(e.1) – rotated to ensure that points on the cell boundary fall along all outgoing polylines.
We now let be the set of points generated by moving every point to the coordinate , and make the following observations:
-
(obs. 1) No set of points in can be covered by a unit disk in (case 1), a single fixed-area axis-aligned square in (case 2), or a fixed-area equilateral triangle in (case 3).
-
(obs. 2) Letting be the number of point sets in corresponding to a -fold enlargement of point sets in arising from a “-in--SAT clause gadget” instance, there must be at least the same number of geometric objects in (case 1) through (case 3) covering distinct points each. Furthermore, if there are exactly instances of geometric objects in (case 1) through (case 3) covering distinct points each, then each of these geometric objects must contain some translation – due to a cell identification operation when constructing , and the subsequent -fold enlargement of when generating – of a point originally at coordinate in the “-in--SAT clause gadget”.
-
(obs. 3) If we require that all but geometric objects from (case 1) through (case 3) cover three distinct points each in , then the only possible manner (up to the uniqueness of the induced 2D-GCP or PD-2D-GCP) of covering the point sets in belonging to distinct instances (i.e., not falling along cell boundaries) of the “-in--SAT clause gadget”, “variable gadget”, “horizontal / vertical connector gadget”, “bend connector gadget”, or “-shift horizontal / vertical connector gadget” are shown in Fig. 3(a.2–a.4), Fig. 3(b.2–b.3), Fig. 3(c.2–c.4), Fig. 3(d.2–d.4), and Fig. 3(e.2–e.4), respectively. Here, unit disks for (case 1) are colored (blue), fixed-area axis-aligned squares are drawn as (black hollow squares), and fixed-area equilateral triangles correspond to (red hollow triangles). Briefly, under a constraint that all but geometric objects cover points each, it will otherwise be necessary in each case for some pair of geometric objects to cover the same point, where an example of this for the “variable gadget” is illustrated in Fig. 3(b.4) with the two (gray) disks sharing a common point in .
-
(obs. 4) Due to (step 5) of constructing , we can interpret the unit disk (case 1), fixed-area axis-aligned square (case 2), and fixed-area equilateral triangle (case 3) covers for the “-in--SAT clause gadget” in Fig. 3(a.2), Fig. 3(a.3), and Fig. 3(a.4) as corresponding to clauses with literals evaluating to {“True”, “False”, “False”}, {“False”, “True”, “False”}, and {“False”, “False”, “True”}, respectively. In this context, the output of “True” is propagated along polylines entering the cell for the gadget from the North in Fig. 3(a.2), West in Fig. 3(a.3), and East in Fig. 3(a.4). Correspondingly, we can interpret the unit disk (case 1), fixed-area axis-aligned square (case 2), and fixed-area equilateral triangle (case 3) covers for the “variable gadget” in Fig. 3(b.2) and Fig. 3(b.3) as variable assignments of “True” and “False”, respectively.
Putting (obs. 1) through (obs. 4) together, in (case 1) through (case 3) we can determine that each CPP--Conn--in--SAT satisfying assignment being counted by will correspond to the same number of geometric objects inducing a specific 2D-GCP or PD-2D-GCP, and furthermore, that geometric objects will be required to cover if no satisfying assignment exists. In (case 1) through (case 3) we moreover have a bijection between the satisfying assignments counted by and 2D-GCPs, as well as a bijection between the satisfying assignments counted by and PD-2D-GCPs for the same sets of geometric objects, as our reduction ensures these objects will be pairwise disjoint.
It remains to address (case 4). Here, we can determine that the minimum bounding disk for any points from the point sets in will have radius at least . Since , we can accordingly move each point to the coordinate to create a new point set where every points have a minimum bounding disk of radius at least . Finally, as we are proportionally increasing the radius of the unit disk from (case 1) to , we can observe that (case 4) is equivalent to (case 1) in the context of our construction.
Corollary 3.
Letting be a set of points embedded in , letting be the set of 2D-GCPs (or equivalently in the contexts we consider, PD-2D-GCPs) of induced by geometric objects in (case 1) through (case 4) of Theorem 2, and under the assumption that at most points can be covered by a geometric object, the following holds: (1) it is NP-hard and #P-hard under parsimonious reductions to decide whether and to determine , respectively; (2) no -time algorithm exists for deciding whether and for determining under the ETH and #ETH, respectively.
Proof.
First, given the polynomial time reduction from #CPP--Conn--in--SAT to counting 2D-GCPs (or equivalently in the contexts we consider, PD-2D-GCPs), in (case 1) through (case 4) of Theorem 2, the fact that CPP--in--SAT is NP-complete [38] establishes the NP-hardness of deciding whether , and Theorem 1 together with Corollary 2 establish that #CPP--Conn--in--SAT is #P-hard under parsimonious reductions. Next, establishing that no -time algorithm exists for deciding whether and determining under the ETH and #ETH, respectively, requires a bit more care. Naively, letting , it would appear that we require time to transitively parsimoniously reduce a variable and clause instance of -SAT to counting 2D-GCPs or to counting PD-2D-GCPs, as there is a quadratic overhead to reduce -SAT to #Planar--SAT via the reduction of Lichtenstein [36], and then another quadratic overhead in computing an orthogonal integer lattice embedding when reducing #Planar--SAT to counting 2D-GCPs or PD-2D-GCPs via the method given in the Theorem 2 proof argument. However, we can nevertheless achieve the required overhead by making use of the fact that the orthogonal integer lattice embedding algorithm of Biedl & Kant [7] can compute, in time, both area non-planar orthogonal integer lattice embeddings (i.e., embeddings with crossing polylines) and area embeddings of planar graphs with a specified outer face.
Specifically, we can: (step 1) follow exactly the Theorem 1 and Corollary 2 reductions to reduce an input -SAT formula to a -in--SAT (not necessarily planar) formula having a cubic clause-variable incidence graph with clauses containing only positive literals; (step 2) compute an orthogonal integer lattice embedding of using the algorithm of Biedl & Kant [7]; (step 3) use the Theorem 1 and Corollary 2 reduction to modify the -SAT to Planar--SAT crossing gadget of Lichtenstein [36] to correspond to a CPP--Conn--in--SAT formula, then compute an orthogonal integer lattice embedding of this gadget having the same outer face as in Lichtenstein’s original embedding [36]; and (step 4) for some constant factor enlargement of the clause-variable incidence graph for , use the gadget constructed in (step 3) to planarize the embedding from (step 2) while preserving its orthogonality. Here, (step 1) will take time, (step 2) will take time, (step 3) will take time, and finally, (step 4) will take at most time as there will be at most possible crossing polylines that require the placement of the gadget built in (step 3). Putting everything together, we can now rule out the existence of a -time algorithm for deciding whether and for determining under the ETH and #ETH, respectively.
5 Covering sparse sets of points with semi-algebraic subsets of a disk
In this section, we investigate density constraints on point sets in that allow for efficient algorithms to cover the points in – with or without a pairwise disjointness constraint – using a minimum cardinality set of translatable and/or rotatable instances of a semi-algebraic subregion of a disk (see Theorem 3). Due to the technical nature of the lemmas in this section (i.e., Lemma 4 through Lemma 6), as opposed to reading straight through, we recommend referring to them as they arise in the context of the proof argument for Theorem 3.
Lemma 4.
Letting be a set of unit disks embedded in , under the assumption that at most a fixed number of disks in can have center points in a minimum bounding disk of some fixed radius , the geometric intersection graph induced by can be computed in time.
Proof.
Let be the geometric intersection graph induced by , and observe that any requirement at most unit disk centers have a minimum bounding disk of some radius necessarily implies that . The current lemma now follows from the fact that we can construct a geometric intersection graph of unit disks having edges in time via a rounding and hash table approach [5].
Lemma 5.
Let be any finite set of points in , and let be a region admitting a description as a semi-algebraic set. In addition, let be the set of polynomial inequalities and/or equalities specifying and some fixed number of instances of , let be the number of quantified variables in , and let be the maximum degree of any polynomial in . With or without a pairwise disjointness constraint, there exists a constructive -time algorithm for finding individual translations and/or rotations of each instance of covering all points in or deciding that this is impossible.
Proof.
Let be a translation vector, let be a rotation angle, let be a rotation matrix with satisfying the constraint that , and let be a point in a subregion admitting a description as a semi-algebraic set. Observe that we can express any translation and/or rotation of as the function . Also observe that we can express the problem of translating and/or rotating to cover a set of points as a quantifier elimination problem for the first-order sentence over the reals “”. This now generalizes in a straightforward manner to instances where we have at most a fixed number of individually translatable and/or rotatable instances of , and where we can also choose to enforce a pairwise disjointness constraint by expressing in first-order logic that no point belongs to the closed set of points for any two translated and/or rotated instances of .
We now appeal to Collins’s (constructive) Cylindrical Algebraic Decomposition (CAD) algorithm for real quantifier elimination [12], which improves upon the original (also constructive) method of Tarski [44, 45]. Letting , , and be as stated, to establish the current lemma it remains to observe that Collins’s algorithm [12] has time complexity [14].
Lemma 6.
Letting be a connected geometric intersection graph of at most unit disks, the center points of the disks can always be covered by a pairwise interior-disjoint arrangement of unit disks such that at least one unit disk is tangent to the other two.
Proof.
Let be the set of center points for the unit disks inducing . The case where is trivial. For cases where , place a unit disk at the center of one edge to cover a pair of points, and if necessary, use a second unit disk tangent to the first to cover the third point, which must always be possible as is connected. In the case of distinct points, , assume without loss of generality that and induce a pair of adjacent vertices in , and use a single unit disk to cover both and . If also covers or , then we can proceed as in the cases where . Alternatively, if neither nor are covered by , then we can attempt to cover (respectively, ) with unit disks (respectively, ), such that and are tangent to and have no totally interior point in common.
Now assume that none of the aforementioned cases hold. To proceed, assume without loss of generality that (illustrated by the (red) unit disk in Fig. 4(a)) has its center point at the coordinate , fixed in position by at the coordinate and at the coordinate (illustrated by the diameter line between points and in Fig. 4(a)), and assume that is further from than . Now cover with a unit disk (illustrated by the (blue) disk in Fig. 4(a)), doing so in such a manner as to embed the center point of as far as possible from while ensuring and are tangent. Finally, embed a unit disk (illustrated by the (green) disk in Fig. 4(a)) such that it is tangent to both and , and has a center point as close as possible to . Here, recalling that we are assuming no three points are interior to a common unit disk, we can observe that will fall along the circumference of bordering the “hole” (i.e., the circular triangle) formed by the three mutually tangent disks, and must be an interior point in this hole. However, as illustrated in Fig. 4(a), letting be the tangency point of and and be the tangency point of and , we can observe that will form a triangle with side lengths of , , and at most . As this triangle fits in a unit disk and contains the points , this contradicts our earlier assumption that no three points can be covered by a single unit disk, yielding the current lemma.
Theorem 3.
For any finite set of points , letting be any semi-algebraic subset of a radius disk, under either of the following assumptions –
-
(constraint 1) any points have a minimum bounding disk of radius greater than ;
-
(constraint 2) any points have a minimum bounding disk of radius at least ;
– it is fixed-parameter tractable in the size of the semi-algebraic description of to find a minimum cardinality set (or determine that no such set exists) of translatable and rotatable instances of covering all points in both with and without a pairwise disjointness constraint.
Proof.
Let be an order geometric intersection graph of radius disks, each with a center at a distinct point in . Observe that, under both (constraint 1) and (constraint 2), Lemma 4 guarantees can be constructed in time. Also observe that, because is a subset of a radius disk, can only cover a collection of points in inducing a connected component in .
To first treat the case of (constraint 1), as a result of a straightforward triangle inequality argument, observe that this assumption ensures all connected components in have order at most . The assumption also ensures that a radius disk covering a pair of distinct points cannot have a point in common with another radius disk covering a point . Accordingly, no method of covering the points in one order-at-most- connected component can prohibit a method of covering the points in another order-at-most- connected component. Thus, we have reduced our task to computing an optimal cover of at most order-at-most- connected components using at most translated and/or rotated instances of , which, by Lemma 5, is fixed-parameter tractable in the size of the semi-algebraic description of . This yields the required fixed-parameter tractable algorithm in the case of (constraint 1).
Concerning now (constraint 2), first observe that, under this assumption, any set of points in inducing a connected path subgraph of will have a maximum bounding disk of radius , with no larger radius possible, if and only if these points are embedded at unit intervals along a line. Here, the if direction is trivial to verify, and the only if direction follows from the uniqueness of the equality case for the triangle inequality. As any other embedding of implies the existence of a minimum bounding disk of radius less than , this implies that any connected component in of order at least will correspond to points embedded at unit intervals along a line.
Next, observe that, under (constraint 2), Lemma 6 guarantees that any set of points in inducing an order-at-most- connected component in can always be covered using at most radius disks (i.e., the radius of the original disks from Lemma 6) having no totally interior points in common, and where at least one of these disks is tangent to the other two. As a consequence, no method of covering the points in one order-at-most- connected component can prohibit a method of covering the points in another order-at-most- connected component. To see this, let , , and be a set of three radius disks covering some set of points inducing an order-at-most- connected component in , and let be tangent to and . Observe that the existence of a radius disk covering a point , and having a totally interior point in common with , , or , will necessarily imply that the minimum bounding disk for has radius less than . For visual intuition we refer the reader to Fig. 4(b), where (colored (blue)), (colored (orange)), and (colored (green)) have collinear center points, maximizing the possible length of a line segment in their union. Nevertheless, when we embed (colored (red)) such that it has a totally interior point in common with , the diameter of a minimum bounding disk for the center points of through – illustrated as the line between points and in Fig. 4(b) – will be less than , implying a radius less than .
Thus, we have reduced our task to computing an optimal cover by of order-at-most- connected components, which Lemma 5 guarantees is fixed-parameter tractable in the size of the semi-algebraic description of , or of points embedded at unit intervals along a line. In the latter case, in combination with Lemma 5, it suffices to observe that a rotation of can cover at most points when embedded at the center of the unit interval between any two points, that a maximum matching for a path graph can be computed in time linear in its order, and that there is sufficient space for each point to be covered by a distinct instance of . This resolves (constraint 2), yielding the current theorem.
6 Concluding remarks
Recall that, on the one hand, (case 4) of Theorem 2 rules out the existence of a sub-exponential algorithm (assuming the ETH) for finding a minimum cardinality set of disks of radius covering a point set (with or without a pairwise disjointness constraint) under the additional constraint that any points in have a minimum bounding disk of radius at least . On the other hand, Theorem 3 tells us that this same problem can be solved in time if we use disks of radius to cover the points in . Accordingly we ask the following question: for a point set where any points have a minimum bounding disk of radius at least (or some other constant), and letting be some disk of radius , when can one efficiently find a minimum number of translated copies of covering all points in ? Furthermore, given that all of the results we report in this work hold with and without a ply constraint, we pose the question of which specific conditions make a pairwise disjointness constraint significant.
References
- [1] K. I. Aardal, S. P. M. van Hoesel, A. M. C. A. Koster, C. Mannino, and A. Sassano. Models and solution techniques for frequency assignment problems. Ann. Oper. Res., 153(1):79–129, 2007. doi:10.1007/s10479-007-0178-0.
- [2] G. Aloupis, R. A. Hearn, H. Iwasawa, and R. Uehara. Covering points with disjoint unit disks. Proc. 24th Canadian Conference on Computational Geometry (CCCG), pages 41–46, 2012.
- [3] P. Ashok, A. B. Roy, and S. Govindarajan. Local search strikes again: PTAS for variants of geometric covering and packing. J. Comb. Optim., 39:618–635, 2020. doi:10.1007/s10878-019-00432-y.
- [4] S. Bandyapadhyay, W. Lochet, S. Saurabh, and J. Xue. Minimum-membership geometric set cover, revisited. Proc. 39th Symposium on Computational Geometry (SoCG), pages 11:1–11:14, 2023. doi:10.4230/LIPIcs.SoCG.2023.11.
- [5] J. L. Bentley, D. F. Stanat, and E. H. Williams Jr. The complexity of finding fixed-radius near neighbors. Inf. Process. Lett., 6(6):209–212, 1977. doi:10.1016/0020-0190(77)90070-9.
- [6] T. Biedl, A. Biniaz, and A. Lubiw. Minimum ply covering of points with disks and squares. Comput. Geom., 94:101712:1–101712:11, 2021. doi:10.1016/j.comgeo.2020.101712.
- [7] T. Biedl and G. Kant. A better heuristic for orthogonal graph drawings. Comput. Geom., 9(3):159–180, 1998. doi:10.1016/S0925-7721(97)00026-6.
- [8] A. Biniaz and Z. Lin. Minimum ply covering of points with convex shapes. Proc. 32nd Canadian Conference on Computational Geometry (CCCG), pages 2–5, 2020.
- [9] A. Biniaz, P. Liu, A. Maheshwari, and M. Smid. Approximation algorithms for the unit disk cover problem in 2D and 3D. Comput. Geom., 60:8–18, 2017. doi:10.1016/j.comgeo.2016.04.002.
- [10] J. A. Bondy and U. S. R. Murty. Graph theory with applications (1st edition). Macmillan Press, New York, NY, 1976.
- [11] J. H. Chun, C. Kipp, and S. Roch. On exact covering with unit disks. Comput. Geom., 129:102193:1–102193:19, 2025. doi:10.1016/j.comgeo.2025.102193.
- [12] G. E. Collins. Quantifier elimination for real closed fields by cylindrical algebraic decompostion. Proc. Automata theory and formal languages: 2nd GI conference, pages 134–183, 1975. doi:10.1007/3-540-07407-4_17.
- [13] P. F. Cortese and M. Patrignani. Clustered planarity = flat clustered planarity. Proc. 26th International Symposium on Graph Drawing (GD), pages 23–38, 2018. doi:10.1007/978-3-030-04414-5_2.
- [14] J. H. Davenport and J. Heintz. Real quantifier elimination is doubly exponential. J. Symbolic Computation, 5(1–2):29–35, 1988. doi:10.1016/S0747-7171(88)80004-X.
- [15] D. de Werra and Y. Gay. Chromatic scheduling and frequency assignment. Discrete Appl. Math., 49(1–3):165–174, 1994. doi:10.1016/0166-218X(94)90207-0.
- [16] H. Dell, T. Husfeldt, D. Marx, N. Taslaman, and M. Wahlén. Exponential time complexity of the permanent and the Tutte polynomial. ACM Trans. Algorithms, 10(4):21:1–21:32, 2014. doi:10.1145/2635812.
- [17] R. Diestel. Graph theory (5th edition). Springer-Verlag, Heidelberg, 2017.
- [18] S. Durocher, J. M. Keil, and D. Mondal. Minimum ply covering of points with unit squares. Proc. 17th International Conference and Workshops on Algorithms and Computation (WALCOM), pages 23–35, 2023. doi:10.1007/978-3-031-27051-2_3.
- [19] S. Durocher, J. M. Keil, and D. Mondal. Approximation algorithms for minimum ply covering of points with unit squares and unit disks. Theoret. Comput. Sci., 1024:114906:1–114906:14, 2025. doi:10.1016/j.tcs.2024.114906.
- [20] A. Eisenblätter, M. Grötschel, and A. M. C. A. Koster. Frequency planning and ramifications of coloring. Discuss. Math. Graph Theory, 22(1):51–88, 2002.
- [21] V. Elser. Packing-constrained point coverings. Geombinatorics Quarterly, 21:73–80, 2011. arXiv:1101.3468.
- [22] T. Erlebach and E. J. van Leeuwen. Approximating geometric coverage problems. Proc. 19th Symposium on Discrete Algorithms (SODA), pages 1267–1276, 2008. doi:10.5555/1347082.1347220.
- [23] M. Formann and F. Wagner. A packing problem with applications to lettering of maps. Proc. 7th Symposium on Computational Geometry (SoCG), pages 281–288, 1991. doi:10.1145/109648.109680.
- [24] R. J. Fowler, M. S. Paterson, and S. L. Tanimoto. Optimal packing and covering in the plane are NP-complete. Inf. Process. Lett., 12(3):133–137, 1981. doi:10.1016/0020-0190(81)90111-3.
- [25] S. P. Gomes, G. M. Ribeiro, and L. A. N. Lorena. Dispersion for the point-feature cartographic label placement problem. Expert Syst. Appl., 40(15):5878–5883, 2013. doi:10.1016/j.eswa.2013.04.035.
- [26] T. F. Gonzalez. Covering a set of points in multidimensional space. Inf. Process. Lett., 40(4):181–188, 1991. doi:10.1016/0020-0190(91)90075-S.
- [27] S. Govindarajan, M. Patle, and S. Sarkar. Minimum membership geometric set cover in the continuous setting. Proc. 31st International Conference on Computing and Combinatorics (COCOON), pages 262–276, 2025. doi:10.1007/978-981-95-0215-8_20.
- [28] S. Govindarajan and S. Sarkar. Improved algorithms for minimum-membership geometric set cover. Proc. 10th Conference on Algorithms and Discrete Applied Mathematics (CALDAM), pages 103–116, 2024. doi:10.1007/978-3-031-52213-0_8.
- [29] W. K. Hale. Frequency assignment: theory and applications. Proc. IEEE, 68(12):1497–1514, 1980. doi:10.1109/PROC.1980.11899.
- [30] M. M. Hasan, D. Mondal, and M. S. Rahman. Positive planar satisfiability problems under -connectivity constraints. Theoret. Comput. Sci., 917:81–93, 2022. doi:10.1016/j.tcs.2022.03.013.
- [31] D. S. Hochbaum and W. Maass. Approximation schemes for covering and packing problems in image processing and VLSI. J. ACM, 32(1):130–136, 1985. doi:10.1145/2455.214106.
- [32] G. Illing and U. Klüh, editors. Spectrum auctions and competition in telecommunications. The MIT Press: Cambridge, MA, 2003. doi:10.7551/mitpress/6400.001.0001.
- [33] R. Impagliazzo and R. Paturi. On the complexity of k-SAT. J. Comput. Syst. Sci., 62(2):367–375, 2001. doi:10.1006/jcss.2000.1727.
- [34] N. Inaba. “10 points answer” (translated from Japanese), 2008. URL: http://inabapuzzle.com/hirameki/suuri_ans4.html.
- [35] N. Inaba. “10 points” (translated from Japanese), 2008. URL: http://inabapuzzle.com/hirameki/suuri_4.html.
- [36] D. Lichtenstein. Planar formulae and their uses. SIAM J. Comput., 11(2):329–343, 1982. doi:10.1137/0211025.
- [37] G. L. Miller, S. H. Teng, W. Thurston, and S. A. Vavasis. Separators for sphere-packings and nearest neighbor graphs. J. ACM, 44(1):1–29, 1997. doi:10.1145/256292.256294.
- [38] C. Moore and J. M. Robson. Hard tiling problems with simple tiles. Discrete Comput. Geom., 26(4):573–590, 2001. doi:10.1007/s00454-001-0047-6.
- [39] Y. Okayama, M. Kiyomi, and R. Uehara. On covering of any point configuration by disjoint unit disks. Proc. 23rd Canadian Conference on Computational Geometry (CCCG), pages 1–5, 2011.
- [40] I. Pak and J. Yang. Tiling simply connected regions with rectangles. J. Comb. Theory. Ser. A, 120(7):1804–1816, 2013. doi:10.1016/j.jcta.2013.06.008.
- [41] I. Pak and J. Yang. Tiling simply connected regions with rectangles, 2013. arXiv:1305.2796v1.
- [42] T. Seta. The complexities of puzzles, CROSS SUM and their another solution problems (ASP). Senior thesis, The University of Tokyo, Tokyo, Japan, 2002.
- [43] M. I. Shamos. Computational geometry. Ph.D. dissertation, Yale University, pages 1–236, 1978.
- [44] R. G. Swan. Tarski’s principle and the elimination of quantifiers, 2005. Unpublished; available at https://www.math.uchicago.edu/˜swan/expo/Tarski.pdf.
- [45] A. Tarski. A decision method for elementary algebra and geometry. Technical Report R-109, RAND Corporation, Santa Monica, CA, 1948. Available at https://www.rand.org/content/dam/rand/pubs/reports/2008/R109.pdf.
- [46] L. F. Tóth. Über die dichteste kugellagerung. Math. Z., 48:676–684, 1942. doi:10.1007/BF01180035.
- [47] N. Ueda and T. Nagao. NP-completeness results for NONOGRAM via parsimonious reductions. Technical Report TR96-0008, Department of Computer Science, Tokyo Institute of Technology, 1996.
- [48] M. van Kreveld, T. Strijk, and A. Wolff. Point set labeling with sliding labels. Proc. 14th Symposium on Computational Geometry (SoCG), pages 337–346, 1998. doi:10.1145/276884.276922.
- [49] P. Winkler. Puzzled: figures on a plane. Commun. ACM, 53(8):128, 2010. doi:10.1145/1810891.1810917.
- [50] T. Yato and T. Seta. Complexity and completeness of finding another solution and its application to puzzles. IEICE Trans. Fundamentals, E86-A(5):1052–1060, 2003. URL: http://search.ieice.org/bin/summary.php?id=e86-a_5_1052.
