Abstract 1 Introduction 2 Paper overview 3 Preliminaries and clarifications

Arranging Pairwise Disjoint Shapes to Partition Point Sets

Robert D. Barish222Corresponding author. ORCID Division of Medical Data Informatics, Human Genome Center, Institute of Medical Science, University of Tokyo, Japan    Tetsuo Shibuya ORCID Division of Medical Data Informatics, Human Genome Center, Institute of Medical Science, University of Tokyo, Japan
Abstract

We consider the fine-grained complexity of covering a set of n 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 2o(n)-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 1-in-3-SAT with a cubic planar 3-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 12 disk characterized by a semi-algebraic function, and letting 𝒫 be a set of n points, we consider the density requirements that: (constraint 1) every 3 points have a minimum bounding disk of radius greater than 1; or (constraint 2) any 5 points have a minimum bounding disk of radius at least 2. 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.
Tetsuo Shibuya: JSPS Kakenhi grants {23K28035, 23K18501, 21H05052}.
Copyright and License:
[Uncaptioned image] © Robert D. Barish and Tetsuo Shibuya; licensed under Creative Commons License CC-BY 4.0
2012 ACM Subject Classification:
Mathematics of computing Discrete mathematics
; Theory of computation Problems, reductions and completeness
Editor:
Pierre Fraigniaud

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 n10 points in 2 can be covered by pairwise interior-disjoint unit disks. This can be seen by noting that unit disks pack with an optimal density of ρ=π23 in a hexagonal lattice arrangement [46], by using the union bound (i.e., Boole’s inequality) to determine the probability that n points in the plane are covered by a randomly rotated and translated hexagonal lattice packing of unit disks is at least (1n(1ρ)), and finally by observing that (1n(1ρ))>0 if n<6(63π)10.7411 (see, e.g., [39]). In the years since, this bound has been shown to hold for n12 [2], and the upper bound for the number of points in 2 that cannot always be covered by pairwise interior-disjoint unit disks, originally set at 60 due to a triangular lattice construction of Winkler [49], has fallen from 55 [21], to 53 [39], and most recently to 45 [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., 106 Hz for AM radio to 109 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 n10 points, in this work we attempt to make progress towards its resolution by establishing the hardness of covering a specified set of points 𝒫 in 2 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 2, and let S be a finite collection of geometric objects embedded in the plane. A 2D Geometric Cover Partition (2D-GCP) induced by S is a partitioning of 𝒫 into subsets {X1,,Xk} (where XiXj= if ij) under the constraint that each subset Xi𝒫 is the set of all points covered by the same geometric object sS.

Definition 2.

Pairwise Disjoint 2D Geometric Cover Partition (PD-2D-GCP). An instance of a 2D-GCP where all geometric objects in S 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 S. Furthermore, when we refer to a 2D-GCP induced by a set of geometric shapes S, we clarify that we are specifically and only referring to the partition of a point set 𝒫 into subsets {X1,,Xk} for some k.

We remark that the task of finding an inducing set S 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 (MS1P-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-1. In particular, we can determine that the NP-hardness proof for MS1P-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-3-Conn-1-in-3-SAT (see Section 3.1 concerning parsimonious reductions and Section 3.4 concerning #CPP-3-Conn-1-in-3-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 106 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-3-Conn-1-in-3-SAT under parsimonious reductions also answers an open question in the literature (see, e.g., [30]) concerning the hardness of 1-in-3-SAT instances with planar 3-connected clause-variable incidence graphs.

Subsequently, letting 𝒫 be a set of n points embedded in 2, we consider the fine-grained complexity of finding and counting 2D-GCPs or PD-2D-GCPs in cases where the set S 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 2o(n)-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 3 points have a minimum bounding disk of radius greater than 1, or (constraint 2) that any 5 points have a minimum bounding disk of radius at least 2 (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 12 disk (where we can allow for translations and/or rotations of Υ) admitting a description as a semi-algebraic set, and letting n 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 Σ={0,1}) be a pair of integer counting problems, and letting x be a string specifying an instance of ϕ, a many-one counting reduction from ϕ to ψ is a pair of polynomial time computable functions f:ΣΣ and h:Σ× such that ϕ(x)=h(x,ψ(f(x))). Here, a parsimonious reduction is a many-one counting reduction where h is simply a projection onto its second argument. Recall that a problem is fixed-parameter tractable with respect to a parameter k, if for some computable function f(k), letting x be a string specifying a problem instance, we have that the time complexity for the problem can be expressed as f(k)poly(|x|).

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 G for a set S of geometric objects embedded in the plane constructs a vertex set V having a bijection with S, and connects a pair of vertices va,vbV if and only if their associated geometric objects in S have a point in common.

3.3 Plane geometry notions and terminology

Recall that a unit disk is a disk of radius 1 embedded in the plane, that a unit line segment is a line segment of length 1 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 {(x,y),(x+k,y),(x+k,y+k),(x,y+k)} for some (x,y)2 and some specified k>0. When we say that a geometric object covers a point p in the plane, unless otherwise specified, we mean that p 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 1-in-3-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 Gϕ=(XY,E), where X (respectively, Y) is in bijection with the variables (respectively, clauses) for ϕ, and where we add an edge between a vertex vX and uY if and only if the variable in ϕ corresponding to v occurs as a (positive or negative) literal in the clause in ϕ corresponding to u. We can now define the Cubic Planar Positive 3-Connected 1-in-3 Satisfiability (CPP-3-Conn-1-in-3-SAT) problem as a special case of 1-in-3-SAT in which the clause-variable incidence graph for any input formula ϕ is restricted to being a cubic (i.e., 3-regular) planar 3-connected graph, and where all literals in ϕ are positive. We also define #CPP-3-Conn-1-in-3-SAT as the #P problem of counting all satisfying assignments for instances of CPP-3-Conn-1-in-3-SAT. Finally, we define Cubic Planar Positive 1-in-3 Satisfiability (CPP-1-in-3-SAT) (respectively, #CPP-1-in-3-SAT) as a relaxation of CPP-3-Conn-1-in-3-SAT (respectively, #CPP-3-Conn-1-in-3-SAT) where we do not require the clause-variable incidence graph to be 3-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 n variable instance of 3-SAT ϕ cannot be solved in 2o(n) time. In this context, using the Sparsification Lemma [33], we can assume that ϕ has at most 𝒪(n) clauses. The Counting Exponential Time Hypothesis (#ETH) of Dell et al. [16] is a relaxation of the ETH where we only assume that no 2o(n)-time algorithm exists for the #3-SAT problem of counting all solutions for a given 3-SAT formula ϕ. Here, the Sparsification Lemma can be adapted [16] to again allow us to assume that ϕ has at most 𝒪(n) clauses.

3.6 Orthogonal integer lattice embeddings

Letting H be a graph of maximum vertex degree 4 having vertex set VH and edge set EH, an orthogonal integer lattice embedding 𝒬 of H is an embedding (or “drawing”) of H on a 2 grid graph 𝒢. To construct 𝒬, we map each vertex viVH to a distinct vertex in 𝒢, where we can then say that this lattice point in 𝒢 hosts or is hosting viVH. Next, for every edge vavbEH, 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 va,vbVH (note that the polyline must not otherwise intersect a vertex in 𝒢 hosting a vertex in H). 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 1×1 square centered on a point (x,y)2 with vertex coordinates ((x12,y12),(x+12,y12),(x+12,y+12),(x12,y+12)). Note that when we identify a cell at a coordinate (x,y)2 with some set of points S (e.g., as we do with various gadgets in the Theorem 2 proof argument), each point piS should be assumed to be translated to the coordinate pi+(x,y). Finally, for some k>0, we define a k-fold enlargement of 𝒬 as the embedding generated by moving all lattice points (x,y) hosting either vertices of an embedded graph or polyline segment endpoints to the coordinate (kx,ky), 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 3-cube from Fig. 1(a) on a 6×6 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 1×1 cell is indicated using a (thick dashed purple) box.

Figure 1: Example orthogonal integer lattice embedding; (a) planar embedding of the 3-cube, where vertices are colored (blue) and (red) to indicate their partite set membership; (b) orthogonal integer lattice embedding of the 3-cube from (a) in a 6×6 integer lattice, where the embedded vertices retain their coloration from (a), (thick black) edges belong to polylines, and the (thick dashed purple) box indicates an example 1×1 cell in the embedding.

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

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-3-Conn-1-in-3-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 𝒪(m)-time parsimonious reduction from any Planar 1-in-3-SAT formula with a 2-connected clause-variable incidence graph and m>0 clauses to a CPP-3-Conn-1-in-3-SAT formula.

Proof.

We proceed via a gadget-based reduction, where letting ϕ be an m clause instance of Planar 1-in-3-SAT with a 2-connected clause-variable incidence graph Gϕ, we will construct the clause-variable incidence graph Hϕ for an 𝒪(m) clause instance ϕ of CPP-3-Conn-1-in-3-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 {p(r,1),p(r,2)} for some r>0. Here, for some r,s>0, we identify pole vertex pairs {p(r,1),p(r,2)} and {p(s,1),p(s,2)} by identifying the degree 2 vertex p(r,1) with the degree 1 vertex p(s,2), and then identifying the degree 2 vertex p(s,1) with the degree 1 vertex p(r,2). 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 Hϕ via the following 6 steps:

  • (step 1) For every degree k=3 variable vertex vi in Gϕ, we construct an instance αi of the “3-Output CPP-3-Conn-1-in-3-SAT variable gadget” shown in Fig. 2(a).

  • (step 2) For every clause vertex ui in Gϕ, we construct an instance βi of the “CPP-3-Conn-1-in-3-SAT clause gadget” shown in Fig. 2(b).

  • (step 3) For every edge ei corresponding to a negated literal in Gϕ, we construct an instance γi of the “CPP-3-Conn-1-in-3-SAT inverter gadget” shown in Fig. 2(c).

  • (step 4) For every degree k=2 or k4 variable vertex vi in Gϕ, we construct instances {q1,,qk} of the “3-Output CPP-3-Conn-1-in-3-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 qj with an unused pole vertex pair on the gadget q((jmodk)+1), yielding a “k-Output CPP-3-Conn-1-in-3-SAT variable gadget” δi with k unused pole vertex pairs.

  • (step 5) For every edge between a variable vertex vi and a clause vertex uj in Gϕ, where the variable corresponding to vi occurs as a positive literal in the clause corresponding to uj, we identify an unused pole vertex pair on αi (or δi if vi has degree 2 or degree 4) with an unused pole vertex pair on βj.

  • (step 6) For every edge ek between a variable vertex vi and a clause vertex uj in Gϕ, where the variable corresponding to vi occurs as a negative literal in the clause corresponding to uj, we identify unused pole vertex pairs on the gadgets αi (or δi if vi has degree 2 or degree 4) and βj with unused pole vertex pairs on the gadget γk.

For the resulting graph Hϕ, each pole vertex pair will correspond to a set of variable vertices in Hϕ 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 1-in-3-SAT. Accordingly, to see that the 1-in-3-SAT formula ϕ corresponding to the clause-variable incidence graph Hϕ represents a formula logically equivalent to ϕ, it suffices to make the following three observations concerning the 1-in-3-SAT formulas induced by the gadgets used to construct Hϕ: (obs. 1) there are exactly two feasible satisfying assignments for each formula induced by a “3-Output CPP-3-Conn-1-in-3-SAT variable gadget” or “k-Output CPP-3-Conn-1-in-3-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 Gϕ); (obs. 2) there are exactly three feasible satisfying assignments for each formula induced by a “CPP-3-Conn-1-in-3-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 Gϕ); and (obs. 3) there are exactly two feasible satisfying assignments for formulas induced by the “CPP-3-Conn-1-in-3-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 Gϕ).

It remains to observe that Hϕ will necessarily be a cubic planar 3-connected graph corresponding to the clause-variable incidence graph of a CPP-3-Conn-1-in-3-SAT formula, and that (step 1) through (step 6) take 𝒪(m) time.

Figure 2: Gadgets used in reducing Planar 1-in-3-SAT with a 2-connected clause-variable incidence graph to CPP-3-Conn-1-in-3-SAT, where (green) vertices correspond to 1-in-3-SAT clauses and (blue) vertices correspond to variables; (a)3-Output CPP-3-Conn-1-in-3-SAT variable gadget”; (b) “CPP-3-Conn-1-in-3-SAT clause gadget”; (c) “CPP-3-Conn-1-in-3-SAT inverter gadget”; (d) scheme for constructing a “k-Output CPP-3-Conn-1-in-3-SAT variable gadget” for arbitrary k>1, where (purple) edges connect vertices to be identified via a pole vertex pair identification operation, and where (black dashed) edges indicate a gadget cycle containing k total copies of the “3-Output CPP-3-Conn-1-in-3-SAT variable gadget”.
Corollary 2.

There exists an 𝒪(m2)-time parsimonious reduction from any 3-SAT formula with m>0 clauses to a CPP-3-Conn-1-in-3-SAT formula with 𝒪(m2) clauses.

Proof.

Observe first that we can parsimoniously reduce any m clause instance ψ of 3-SAT to an 𝒪(m)-clause instance ϕ of 3-SAT with a 2-connected clause-variable incidence graph in 𝒪(m) time by: (1) constructing two instances ψ1 and ψ2 of ψ; (2) for every variable xi in ψ1 corresponding to a variable yi in ψ2, constructing the formula Zi=(xi¬yi)(¬xiyi); and (3) parsimoniously reducing the formula ψ1ψ2i(Zi) to a 3-SAT formula ϕ via standard methods. Next, recall that the reduction of Lichtenstein [36] can be used to parsimoniously reduce the 𝒪(m) clause instance ϕ of 3-SAT to an m𝒪(m2) clause instance ϕ of Planar 3-SAT in 𝒪(m2) time, and moreover, ensure that the clause-variable incidence graph of ϕ is 2-connected if the clause-variable incidence graph of ϕ is 2-connected. Thus, by Theorem 1, it now suffices to show the existence of an 𝒪(m)-time reduction from ϕ to an instance ϕ′′ of Planar 1-in-3-SAT with a 2-connected clause-variable incidence graph. To proceed, let f be a 1-in-3-SAT clause function which accepts a set of 3 literals and returns “True” if and only if exactly one literal in the set is “True”, and observe that the 3-SAT clause (p1p2p3) is logically equivalent to the 1-in-3-SAT formula f(p1,x1,x2)f(p2,x2,x3)f(¬p3,x2,x4)f(x1,x3,x5). Observe now that, after applying this reduction to an input Planar 3-SAT formula ϕ, the clause-variable incidence graph for the resulting Planar 1-in-3-SAT formula ϕ′′ will be planar and 2-connected if the clause-variable incidence graph for ϕ is planar and 2-connected.

Corollary 3.

The following problems are #P-complete under parsimonious reductions: (P1) #CPP-3-Conn-1-in-3-SAT; (P2) counting tilings of finite simply connected regions using at most 106 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-1-in-3-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 n be the number of variables for an instance #ϕ of #CPP-3-Conn-1-in-3-SAT, and letting 𝒫 be a set of 𝒪(n2) points embedded in 2, there exists an 𝒪(n2)-time parsimonious reduction from #ϕ to counting either 2D-GCPs (or equivalently in the contexts we consider, PD-2D-GCPs) induced by a fixed number k, where k𝒪(n2), 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 65 under the additional constraint that any 5 points in 𝒫 have a minimum bounding disk of radius at least 2.

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-3-Conn-1-in-3-SAT with n variables, m𝒪(n) clauses, and a cubic planar bipartite clause-variable incidence graph Gϕ, where Vclause and Vvar are the partite sets of vertices in Gϕ corresponding to clauses and variables, respectively. Using the algorithm of Biedl & Kant [7], we compute an orthogonal integer lattice embedding 𝒬1 of Gϕ having total area 𝒪(n2) in 𝒪(n) time, and let 𝒬2 be a 3-fold enlargement of 𝒬1 (see Section 3.6 for an elaboration).

Next, we construct a set of points 𝒫2 from the cells in 𝒬2, 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 Ω2 via the following steps:

  • (step 1) we identify every cell hosting a vertex vVclause with the “1-in-3-SAT clause gadget” having point set {(1,0), (34,0), (59100,740), (720,740), (0,0), (0,34), (0,1), (740,720), (740,59100), (720,740), (59100,740), (34,0), (1,0)} – 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 vVvar with the “variable gadget” having point set { (1,0), (34,0), (916,0), (14,0), (0,14), (0,916), (0,34), (0,1), (14,0), (916,0), (34,0), (1,0)} – 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 vVclauseVvar and incident to two horizontal or two vertical polylines, we identify the cell with the “horizontal / vertical connector gadget” having point set {(1,0), (34,0), (12,0), (14,0), (0,0), (14,0), (12,0), (34,0), (1,0)} – 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 vVclauseVvar and incident to a horizontal and a vertical polyline, we identify the cell with the “bend connector gadget” having point set {(0,1), (0,34), (110,113200), (87400,151400), (151400,87400), (113200,116), (34,0), (1,0)} – 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 𝒬2, to ensure that (|ζ|2)0(mod3) we replace at most 2 of the 9 point “horizontal / vertical connector gadgets” (which will necessarily exist since 𝒬2 is a 3-fold enlargement of 𝒬1) with the 10 point “+1-shift horizontal / vertical connector gadget” having point set {(1,0), (34,0), (58,0), (38,0), (18,0), (18,0), (38,0), (58,0), (34,0), (1,0)} – 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 (x,y)Ω to the coordinate (4x,4y), and make the following observations:

  • (obs. 1) No set of 4 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 4-fold enlargement of point sets in Ω arising from a “1-in-3-SAT clause gadget” instance, there must be at least the same number φ of geometric objects in (case 1) through (case 3) covering 2 distinct points each. Furthermore, if there are exactly φ instances of geometric objects in (case 1) through (case 3) covering 2 distinct points each, then each of these geometric objects must contain some translation – due to a cell identification operation when constructing Ω, and the subsequent 4-fold enlargement of Ω when generating 𝒫 – of a point originally at coordinate (0,0) in the “1-in-3-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 “1-in-3-SAT clause gadget”, “variable gadget”, “horizontal / vertical connector gadget”, “bend connector gadget”, or “+1-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 3 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 “1-in-3-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-3-Conn-1-in-3-SAT satisfying assignment being counted by #ϕ will correspond to the same number k of geometric objects inducing a specific 2D-GCP or PD-2D-GCP, and furthermore, that k+1 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 5 points from the point sets in 𝒫 will have radius at least 53. Since (53)(65)=2, we can accordingly move each point (x,y)𝒫 to the coordinate (6x5,6y5) to create a new point set 𝒫 where every 5 points have a minimum bounding disk of radius at least 2. Finally, as we are proportionally increasing the radius of the unit disk from (case 1) to 65, we can observe that (case 4) is equivalent to (case 1) in the context of our construction.

Figure 3: Illustrations for the Theorem 2 parsimonious reduction from #CPP-3-Conn-1-in-3-SAT to counting 2D-GCPs (or equivalently in the contexts we consider, PD-2D-GCPs) induced by translatable unit disks (blue or gray disks), translatable fixed-area axis-aligned squares (black hollow squares), or translatable and rotatable fixed-area equilateral triangles (red hollow triangles); (a.1)–(a.4) illustrations for the “1-in-3-SAT clause gadget” point set; (b.1)–(b.4) illustrations for the “variable gadget” point set; (c.1)–(c.4) illustrations for the “horizontal / vertical connector gadget” point set; (d.1)–(d.4) illustrations for the “bend connector gadget” point set; (e.1)–(e.4) illustrations for the “+1-shift horizontal / vertical connector gadget” point set; see the proof argument for Theorem 2 for further details.
Corollary 3.

Letting 𝒫 be a set of n points embedded in 2, letting Z be the set of 2D-GCPs (or equivalently in the contexts we consider, PD-2D-GCPs) of 𝒫 induced by k𝒪(n) geometric objects in (case 1) through (case 4) of Theorem 2, and under the assumption that at most 3 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 |Z|1 and to determine |Z|, respectively; (2) no 2o(n)-time algorithm exists for deciding whether |Z|1 and for determining |Z| under the ETH and #ETH, respectively.

Proof.

First, given the polynomial time reduction from #CPP-3-Conn-1-in-3-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-1-in-3-SAT is NP-complete [38] establishes the NP-hardness of deciding whether |Z|1, and Theorem 1 together with Corollary 2 establish that #CPP-3-Conn-1-in-3-SAT is #P-hard under parsimonious reductions. Next, establishing that no 2o(n)-time algorithm exists for deciding whether |Z|1 and determining |Z| under the ETH and #ETH, respectively, requires a bit more care. Naively, letting n=q2, it would appear that we require 𝒪(q4) time to transitively parsimoniously reduce a q variable and 𝒪(q) clause instance of #3-SAT to counting 2D-GCPs or to counting PD-2D-GCPs, as there is a quadratic overhead to reduce #3-SAT to #Planar-3-SAT via the reduction of Lichtenstein [36], and then another quadratic overhead in computing an orthogonal integer lattice embedding when reducing #Planar-3-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 𝒪(q2) overhead by making use of the fact that the orthogonal integer lattice embedding algorithm of Biedl & Kant [7] can compute, in 𝒪(q) time, both 𝒪(q2) area non-planar orthogonal integer lattice embeddings (i.e., embeddings with crossing polylines) and 𝒪(q2) 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 3-SAT formula ϕ to a 1-in-3-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 3-SAT to Planar-3-SAT crossing gadget of Lichtenstein [36] to correspond to a CPP-3-Conn-1-in-3-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 𝒪(q) time, (step 2) will take 𝒪(q2) time, (step 3) will take 𝒪(1) time, and finally, (step 4) will take at most 𝒪(q2) time as there will be at most 𝒪(q2) 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 2o(n)-time algorithm for deciding whether |Z|1 and for determining |Z| 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 2 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 S be a set of n unit disks embedded in 2, under the assumption that at most a fixed number of disks in S can have center points in a minimum bounding disk of some fixed radius 1, the geometric intersection graph induced by S can be computed in 𝒪(n) time.

Proof.

Let G=(V,E) be the geometric intersection graph induced by S, and observe that any requirement at most k unit disk centers have a minimum bounding disk of some radius r1 necessarily implies that |E|𝒪(|V|). The current lemma now follows from the fact that we can construct a geometric intersection graph of unit disks having 𝒪(n) edges in 𝒪(n) time via a rounding and hash table approach [5].

Lemma 5.

Let 𝒫 be any finite set of points in 2, and let Υ2 be a region admitting a description as a semi-algebraic set. In addition, let S 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 S, and let d be the maximum degree of any polynomial in S. With or without a pairwise disjointness constraint, there exists a constructive 𝒪((|S|d)2ϑ)-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 (zx,zy) be a translation vector, let θ[0,2π] be a rotation angle, let R(θ)=[cos(θ)sin(θ)sin(θ)cos(θ)]=[αββα] be a rotation matrix with α,β satisfying the constraint that α2+β2=1, and let q be a point in a subregion Υ2 admitting a description as a semi-algebraic set. Observe that we can express any translation and/or rotation of Υ as the function f(zx,zy,α,β)(q)R(θ)q+(zx,zy). 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 “zx,zy,α,β(α2+β2=1p𝒫(pf(zx,zy,α,β)(Υ)))”. 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 S, ϑ, and d be as stated, to establish the current lemma it remains to observe that Collins’s algorithm [12] has time complexity 𝒪((|S|d)2ϑ) [14].

Lemma 6.

Letting G be a connected geometric intersection graph of at most 4 unit disks, the center points of the disks can always be covered by a pairwise interior-disjoint arrangement of 3 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 G. The case where |𝒫|=1 is trivial. For cases where 2|𝒫|3, 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 G is connected. In the case of 4 distinct points, pa,pb,pc,pd𝒫, assume without loss of generality that pa and pb induce a pair of adjacent vertices in G, and use a single unit disk Dred to cover both pa and pb. If Dred also covers pc or pd, then we can proceed as in the cases where |𝒫|3. Alternatively, if neither pc nor pd are covered by Dred, then we can attempt to cover pc (respectively, pd) with unit disks Dblue (respectively, Dgreen), such that Dblue and Dgreen are tangent to Dred 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 Dred (illustrated by the (red) unit disk in Fig. 4(a)) has its center point at the coordinate (0,0), fixed in position by pa at the coordinate (1,0) and pb at the coordinate (1,0) (illustrated by the diameter line between points p1=pa and p2=pb in Fig. 4(a)), and assume that pc is further from pa than pb. Now cover pc with a unit disk Dblue (illustrated by the (blue) disk in Fig. 4(a)), doing so in such a manner as to embed the center point of Dblue as far as possible from pd while ensuring Dred and Dblue are tangent. Finally, embed a unit disk Dgreen (illustrated by the (green) disk in Fig. 4(a)) such that it is tangent to both Dred and Dblue, and has a center point as close as possible to pd. Here, recalling that we are assuming no three points are interior to a common unit disk, we can observe that pc will fall along the circumference of Dblue bordering the “hole” (i.e., the circular triangle) formed by the three mutually tangent disks, and pd must be an interior point in this hole. However, as illustrated in Fig. 4(a), letting p3 be the tangency point of Dred and Dgreen and p4 be the tangency point of Dblue and Dgreen, we can observe that {p2,p3,p4} will form a triangle with side lengths of 1, 2, and at most 3. As this triangle fits in a unit disk and contains the points {pb,pc,pd}, this contradicts our earlier assumption that no three points can be covered by a single unit disk, yielding the current lemma.

Figure 4: Illustrations to accompany the (a) Lemma 6 and (b) Theorem 3 proof arguments.
Theorem 3.

For any finite set of n points 𝒫2, letting Υ be any semi-algebraic subset of a radius 12 disk, under either of the following assumptions –

  • (constraint 1) any 3 points have a minimum bounding disk of radius greater than 1;

  • (constraint 2) any 5 points have a minimum bounding disk of radius at least 2;

– 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 G=(V,E) be an order n geometric intersection graph of radius 12 disks, each with a center at a distinct point in 𝒫. Observe that, under both (constraint 1) and (constraint 2), Lemma 4 guarantees G can be constructed in 𝒪(n) time. Also observe that, because Υ is a subset of a radius 12 disk, Υ can only cover a collection of points in 𝒫 inducing a connected component in G.

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 G have order at most 2. The assumption also ensures that a radius 12 disk covering a pair of distinct points pa,pb𝒫 cannot have a point in common with another radius 12 disk covering a point pc𝒫{pa,pb}. Accordingly, no method of covering the points in one order-at-most-2 connected component can prohibit a method of covering the points in another order-at-most-2 connected component. Thus, we have reduced our task to computing an optimal cover of at most n order-at-most-2 connected components using at most 2 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 5 points in 𝒫 inducing a connected path subgraph H of G will have a maximum bounding disk of radius 2, 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 H implies the existence of a minimum bounding disk of radius less than 2, this implies that any connected component in G of order at least 5 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-4 connected component in G can always be covered using at most 3 radius 12 disks (i.e., 12 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-4 connected component can prohibit a method of covering the points in another order-at-most-4 connected component. To see this, let D1, D2, and D3 be a set of three radius 12 disks covering some set of points X𝒫 inducing an order-at-most-4 connected component in G, and let D1 be tangent to D2 and D3. Observe that the existence of a radius 12 disk D4 covering a point qX, and having a totally interior point in common with D1, D2, or D3, will necessarily imply that the minimum bounding disk for X{q} has radius less than 2. For visual intuition we refer the reader to Fig. 4(b), where D1 (colored (blue)), D2 (colored (orange)), and D3 (colored (green)) have collinear center points, maximizing the possible length of a line segment in their union. Nevertheless, when we embed D4 (colored (red)) such that it has a totally interior point in common with D3, the diameter of a minimum bounding disk for the center points of D1 through D4 – illustrated as the line between points p1 and p2 in Fig. 4(b) – will be less than 4, implying a radius less than 2.

Thus, we have reduced our task to computing an optimal cover by Υ of order-at-most-4 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 2 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 65 covering a point set 𝒫 (with or without a pairwise disjointness constraint) under the additional constraint that any 5 points in 𝒫 have a minimum bounding disk of radius at least 2. On the other hand, Theorem 3 tells us that this same problem can be solved in 𝒪(|𝒫|) time if we use disks of radius 12 to cover the points in 𝒫. Accordingly we ask the following question: for a point set 𝒫 where any k points have a minimum bounding disk of radius at least 2 (or some other constant), and letting Υ be some disk of radius 12<r<65, 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 3-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.