Geometric Optimization Parameterized by Piercing Complexity
Abstract
Packing and Covering problems with geometric regions in the plane have been extensively studied and several notions of “complexity” of the regions involved have been developed and exploited to obtain good approximation algorithms. Examples of such complexity measures are VC-dimension, union complexity, shallow-cell complexity, fatness, etc. While these restrictions lead to constant-factor approximation algorithms in many cases, they typically do not lead to PTASs. In fact, several geometric Set Cover and Discrete Independent Set variants remain APX-hard even when these parameters are small, as demonstrated in earlier work by Chan and Grant (Exact algorithms and APX-hardness results for geometric packing and covering problems. Comput. Geom., 2014), and by Har-Peled and Quanrud (Approximation Algorithms for Polynomial-Expansion and Low-Density Graphs. SIAM J. Comput., 2017). A key feature of these hardness constructions is that many pairs of regions in the input pierce one another.
Motivated by this observation, we initiate a systematic study of geometric families parameterized by their piercing complexity. A connected region is said to pierce a connected region if has more than one connected component; we consider instances in which every region is pierced by at most a constant number of others. This framework smoothly interpolates between the classical non-piercing case-where local-search PTASs are known due to Raman and Ray (Constructing Planar Support for Non-Piercing Regions, Discret. Comput. Geom., 2020), and the fully general case, where APX-hardness persists.
Our main contribution is to show that bounded-piercing families admit efficient approximation schemes for fundamental geometric optimization problems. For regions in the plane with a constant piercing bound, we obtain PTASs for the (unweighted) Discrete Independent Set and Set Cover problems, and constant-factor approximation algorithms for their weighted variants. These results strictly generalize the known PTASs for non-piercing families and yield improved guarantees for several long-standing special cases, including Independent Set and Set Cover with axis-parallel rectangles under bounded piercing.
Overall, our work identifies piercing complexity as a robust and expressive topological parameter-distinct from geometric notions such as density or fatness-and demonstrates that bounding this parameter yields a broad family of geometric instances for which PTASs become achievable.
Keywords and phrases:
Geometric set cover, geometric discrete independent set, approximation algorithms, PTAS, parameterized complexity, piercing complexity, non-piercing regions, independent set, axis-parallel rectanglesCategory:
Track A: Algorithms, Complexity and GamesCopyright and License:
2012 ACM Subject Classification:
Theory of computation Computational geometryEditors:
Sayan Bhattacharya, Danupon Nanongkai, Michael Benedikt, and Gabriele PuppisSeries and Publisher:
Leibniz International Proceedings in Informatics, Schloss Dagstuhl – Leibniz-Zentrum für Informatik
1 Introduction
Given a set of elements and a collection of subsets of whose union covers , the Set Cover problem asks for a smallest sub-collection that also covers . In the Discrete Independent Set problem, the objective is to select a largest sub-collection such that each element of is contained in at most one set in . Both problems have been extensively studied, and tight bounds on their approximability are known in the general setting [15, 29, 25].
Improved approximation algorithms can often be obtained when these problems are restricted to geometric settings. Such improvements are typically enabled by bounding an appropriate notion of complexity of the underlying geometric objects. Well-studied complexity measures include dimension, VC-dimension, union complexity (and its generalization, shallow-cell complexity), and fatness, among others [32, 30, 20, 8]. In particular, small VC-dimension or shallow-cell complexity implies the existence of small -nets [32], which in turn yield approximation algorithms via LP-rounding.
A representative example is the Set Cover problem defined by unit disks and points in the plane, which admits a PTAS [5]. In contrast, when disks are replaced by geometrically similar regions – such as nearly congruent, nearly circular ellipses – the problem becomes APX-hard, despite remaining small with respect to the above complexity measures [7]. Chan and Grant [7] gave several further examples of simple geometric objects, including axis-parallel strips and wedges, for which both the Set Cover and Discrete Independent Set problems are APX-hard (see also [21]). A key structural feature underlying these hardness constructions is that many pairs of regions in the input pierce one another. Informally, a connected region is said to pierce another connected region if has two or more connected components.
In contrast, PTASs are known for a wide range of geometric packing and covering problems. In many cases, these algorithms are based on a simple local-search framework [9, 33, 3, 18, 5, 35]. A common requirement in these settings is that the regions are non-piercing111Examples of non-piercing regions include axis-parallel squares, disks, axis-parallel unit-height rectangles, and homothets of convex regions., meaning that no region pierces another (see Definition 2.1).
The fact that there exist PTASs for non-piercing regions, while known APX-hardness reductions involve many piercing pairs of regions motivates the central question of this paper: can one quantify the amount of piercing in an instance, and thereby interpolate between the tractable non-piercing regime and the intractable general case?
Quantifying the amount of piercing.
We quantify piercing using two parameters. For two regions and , we say that splits into pieces if has connected components. The fragmentation number is the maximum number of pieces into which any region can be split by another, and the piercing number is the maximum number of regions that pierce any single region. A collection of regions with fragmentation number at most and piercing number at most is called an -piercing set. The two parameters play different roles: bounds how many topologically “disruptive” interactions any one region has, while bounds the damage caused by one such interaction. This is a deliberately modest parameterization. For many simple geometric families, such as axis-parallel rectangles, the fragmentation number is already constant, so the piercing number is the parameter that interpolates between the non-piercing regime and highly piercing instances. We do not claim that these are the only possible parameters, but they are sufficient to recover local-search PTASs and they isolate a structural feature absent from standard measures such as VC-dimension or union complexity. This framework interpolates smoothly between the non-piercing case () and the fully general setting. Our main result, stated informally below, shows that bounded piercing suffices to recover PTASs for fundamental geometric optimization problems. Note that even if and are poly-logarithmic in , our results imply QPTASs for the two problems considered in this paper.
Theorem (Informal).
Let be a set of regions forming an -piercing set. For any , there is an algorithm with running time that yields a -approximation for the Set Cover problem defined by and a point set . Analogously, there is an algorithm with the same running time that yields a -approximation for the Discrete Independent Set problem.
Applications.
Many geometric optimization problems employ idealized shape models: sensor coverage as disks or wedges [39], map labels as rectangles [2], and facility service regions defined by -norm balls [19]. These models enable elegant algorithms but poorly model realistic scenarios. Non-piercing regions [35] partially address this by allowing highly irregular and non-convex geometric shapes, but the non-piercing requirement is often violated: obstacles may cause sensor regions to pierce one another, labels for intersecting roads can create piercing pairs, and competing service territories may fragment each other. The -piercing framework should therefore be viewed as a structural, parameterized relaxation of the non-piercing assumption. It applies to instances in which the piercing graph has bounded maximum degree and each piercing has bounded fragmentation; such restrictions can arise from placement rules, design constraints, or simply from the input being only mildly non-piercing. We use these scenarios as motivation rather than as a claim that the piercing number is universally small in practice.
Improved algorithms for special cases.
Our results imply improved approximation guarantees for several well-studied geometric problems. Both the Set Cover and Discrete Independent Set problems are APX-hard for axis-parallel rectangles [7]. For the continuous Independent Set problem, polynomial-time constant-factor approximations [31, 17] and a QPTAS [1] are known, while for the Set Cover problem, the best known approximation is via -nets [6, 14], which is tight under this framework [34]. We show that for both problems, if the piercing number of the rectangles is constant, then a PTAS exists – something not achievable by existing techniques.
A similar phenomenon arises for points and line segments. LP-based methods yield at best an -approximation [6, 14], and lower bounds on -nets rule out substantially better guarantees [4]. At the same time, the problem remains APX-hard even for axis-parallel segments [27]. Our results again yield PTASs for both the Set Cover and Discrete Independent Set problems when the piercing number is bounded.
Comparison with other parameters.
Parameterized algorithm design based on structural measures of the input is an active research direction [12]. Piercing complexity differs fundamentally from existing parameters in that it is purely topological, unlike geometric parameters such as fatness or density, or combinatorial parameters such as VC-dimension, shallow-cell complexity, or treewidth.
Treewidth, even when combined with bounded ply, can be useful for geometric approximation algorithms [28]. However, intersection graphs of -piercing regions can contain arbitrarily large grid minors even with bounded ply, implying unbounded treewidth. Thus, treewidth and piercing complexity are incomparable.
Har-Peled and Quanrud [22] introduced the notion of low-density intersection graphs, which leads to PTASs for several geometric problems. Their framework, however, relies on a geometric notion of object size, whereas -piercing is entirely topological, making the two notions incomparable.
Finally, while -piercing regions may have arbitrarily large union complexity, we show that they nevertheless have shallow-cell complexity linear in and (Theorem 3.3). Since shallow-cell complexity alone is insufficient to yield PTASs [22, 7], this demonstrates that piercing complexity captures a distinct structural property that enables efficient approximation schemes.
Technical contribution.
Our algorithms follow the local-search framework used in earlier work on geometric packing and covering problems. For non-piercing regions, prior work constructs planar supports, and the needed separators follow directly from planarity [35]. With bounded piercing, the corresponding exchange graphs need not be planar. The main technical step is to recover enough sparsity for local search: lens-bypassing operations simplify the arrangement while preserving the relevant coverage conditions and the -piercing property, and the resulting exchange graph is represented as a bounded-depth intersection graph of connected subgraphs of the planar dual arrangement graph. A Frederickson-style decomposition of this graph then supplies the local pieces used in the PTAS analysis. This is the point at which the proof extends, rather than merely reuses, the non-piercing machinery.
1.1 Related Work
For packing problems with fat objects222There are many definitions of fatness, but a definition that is sufficient for us is that region is fat if the ratio of the inscribing and circumscribing balls have radii that are at most a constant factor apart. in bounded dimension, eg. balls or cubes, Erlebach et al. [13] obtained a PTAS building on the work of Hochbaum and Maass [24] for unit disks. For regions that are not fat, Adamaszeck and Wiese [1] obtained a QPTAS (including weighted settings) for regions in the plane under the restriction that each region in the set is path connected. In the discrete setting, Chan and Grant [7] obtained hardness results for packing and covering problems with simple geometric regions in the plane such as axis-parallel strips.
For covering problems, approximation algorithms primarily build on -nets. Haussler and Welzl [23] showed -nets of size for set systems with bounded VC-dimension, leading to -approximations by Brönnimann and Goodrich [6] and Even et al.[14]. Clarkson and Varadarajan [11] improved -nets for geometric systems, culminating in Varadarajan’s [38] quasi-uniform sampling which was refined by Chan et al. [8] for systems with shallow-cell complexity , yielding -nets of size and an -approximation for weighted Set Cover, though with large constants even in simple cases. These algorithms also work in the weighted setting.
For unweighted problems, the unifying tool has been the local search framework, which iteratively improves a feasible solution with small changes. This yielded PTASs for independent set [9], hitting set [33], and Set Cover and dominating set [5]. Raman and Ray [35] showed that set systems from non-piercing regions admit planar supports, implying PTAS results for many packing and covering problems via local search, and also for demand/capacitated versions [36]. Raman and Singh [37] extended these results to higher-genus surfaces. The limitations of the local search approach were studied by Mustafa and Jartoux [26].
2 Preliminaries
In this section, we introduce notation and define terms that we use throughout the paper. We use the term “region” to refer to a closed and bounded set of points in the plane whose boundary is a simple Jordan curve. Thus, each region in our setting is a simply connected region333Raman and Ray [35] considered a slightly more general notion of regions, but for our purposes, we work with only simply connected regions. We denote the boundary of any region by .
Definition 2.1 (Non-piercing set [35]).
A set of regions in the plane is said to be non-piercing if for any , is a connected set.
The families of regions we consider in this paper may contain piercing pairs, where for two regions , the set consists of two or more connected components. We define below the fragmentation number which denotes the number of components in .
Definition 2.2 (Fragmentation number).
Let and be two regions. If their boundaries intersect, we define the fragmentation number of the ordered pair , denoted , as the number of connected components in (see Figure 1). If the boundaries of and do not intersect (which can happen if either one is contained in the other or they are disjoint), we define .
Definition 2.3 (Piercing).
We say that two regions and are piercing if and non-piercing otherwise.
Definition 2.4 (Dual arrangement graph).
Let be an arrangement of regions in the plane. The dual arrangement graph of is a graph in which there is a vertex corresponding to each cell in the arrangement of and two vertices are adjacent iff the boundaries of the corresponding cell touch in an arc of non-zero length. Observe that is a plane multi-graph.
Definition 2.5 (-piercing regions).
We say that a set of regions is -piercing if for any two regions , and ii) for any region , the number of regions that pierce is at most .
Definition 2.6 (Lens [5, 35]).
Let and be any two regions. A lens formed by and is a cell in the arrangement of the two regions that is contained in both and has two sides, one along the boundary of and the other along the boundary of (see Figure 2).
Definition 2.7 (Minimal lens [5, 35]).
Let be any set of regions. We say that lens defined by two of the regions is a minimal lens if does not contain any lens defined by any other pair of regions in .
Definition 2.8 (Expansion [5, 35]).
Let be any collection of regions. For any region (not necessarily in ), we define an “expansion of w.r.t. ” as the region where , denotes an open ball of radius around the origin, and denotes Minkowski sum. Here is chosen to be any appropriately small quantity so that does not contain any vertices of the arrangement of other than those lying in . When is clear from the context, we write instead of .
Definition 2.9 (Lens bypassing [5, 35]).
Let be any collection of regions. Let be a lens formed by the regions and in . We define the operation “bypassing of by ” as the modification of to .
The following definition is inspired by the decomposition theorem for planar graphs due to Frederickson [16].
Definition 2.10.
We say that a family of graphs admits a -decomposition for parameters if for any graph in with vertices, and any given parameter , we can find a set of size at most , and a partition of into sets s.t.
-
(i)
for all
-
(ii)
for all .
-
(iii)
Here, for any , denotes the set of all neighbors of the vertices in in .
Frederickson [16] showed that the family of planar graphs admits a -decomposition for some constants .
Definition 2.11 (Shallow cell complexity of a set of regions [8]).
The shallow cell complexity of a set of regions is a function s.t. for any , the number of combinatorially distinct cells of depth is . Two cells are considered combinatorially equivalent if they are contained in the same set of regions.
Shallow cell complexity is defined more generally in [8] for set systems. For our purposes, the above equivalent definition suffices.
3 Main Results
We now state the formal algorithmic consequences of bounded piercing complexity. Throughout this section, is an -piercing set of regions, where and are as in Definition 2.5. Our main contribution is a pair of local-search PTASs for the two unweighted problems considered in the paper. In the Discrete Independent Set problem, the objective is to select a largest subset of such that any point in is contained in at most one selected region. In the Set Cover problem, the goal is to select a smallest subset of regions such that each point of is in at least one selected region. If and are constants, the algorithms run in polynomial time; if and are polylogarithmic in , the same bounds give QPTASs.
The proofs use the same broad local-search paradigm as the known PTASs for non-piercing regions, but the structural step is different. Instead of obtaining a planar support, we use the lens-bypassing and decomposition lemmas in Section 4 to build exchange graphs that admit the decompositions required by the local-search analysis. The -local search algorithm referred to in the following two theorems is described in Section 5.
Theorem 3.1 (Discrete Independent Set).
The -local search algorithm with for an appropriately large and yields a -approximation for the Discrete Independent Set problem defined by a set of -piercing regions and a set of points . The algorithm runs in time .
Theorem 3.2 (Set Cover).
The -local search algorithm with for an appropriately large and yields a -approximation for the Set Cover problem defined by a set of -piercing regions and a set of points . The algorithm runs in time .
We also show that the shallow cell complexity of a set of regions with the above parameters and is from which the following approximation algorithms for the weighted setting follow via existing results of Chan and Har-Peled [9], and Chan et al. [8].
Theorem 3.3.
Let be an -piercing set of regions and be a set of points. For any point , let denote the subset of regions in containing . Then, the set system (often called the dual set system defined by and ) has shallow-cell complexity .
Note that the above result implies the existence of small -nets for the dual set system and approximation algorithms with corresponding small approximation factors for the set cover problem on the primal set system . See the results in [8, 38].
Theorem 3.4.
Let be a set system defined by a set of points in the plane and an -piercing set of regions , and let be a weight function. Then,
-
(i)
the weighted Set Cover problem admits an -approximation algorithm.
-
(ii)
the weighted Discrete Independent Set problem admits an -approximation algorithm.
4 -Piercing regions
In this section, we describe properties of an -piercing set of regions. Let be an -piercing set. We start with the notion of an inscribed region and use it to bound the number of components a region is split into by regions .
Definition 4.1 (inscribed region).
We say that a region is an inscribed region of another region if and touches at one or more (possibly infinite) points. The boundary consists of arcs that are along and arcs that are chords of . We refer to the chords as the “sides” of w.r.t. and denote their number by .
Lemma 4.2.
Let be any region with disjoint inscribed regions where . Then, the number of connected components in is .
Proof.
We prove the lemma by induction on . The base case is . Since each region of shares exactly one side of , the number of connected components of is . Thus, the lemma holds for . Now, suppose that and the lemma holds for smaller . Since the regions are disjoint, we can draw a chord of which does not intersect any of the regions in and splits into two regions and so that one of the regions lies in and remaining inscribed regions lie in . By the induction hypothesis, the number of connected components in is and the number of connected components in is . This implies that the number of connected components of is . We have a “” because the connected components of and that have on their boundaries correspond to the same connected region in . The lemma follows.
Let and be two regions whose boundaries intersect. Then, observe that every connected component of is an inscribed region of both and and furthermore .
Lemma 4.3.
Let and be two regions whose boundaries intersect. Let be the connected components in . Then, .
Proof.
Since for each , it suffices to prove that which follows directly from Lemma 4.2.
The above lemma implies that for any . If the boundaries of and don’t intersect, by definition we have and if they do then by Lemma 4.3, .
Lemma 4.4.
Let be a region whose boundary intersects the boundaries of a set of disjoint regions . Then, the number of connected components in is . In particular, if at most of the regions in pierce and for each , then the number of connected components in is at most .
Proof.
Let be set of connected components of which are all inscribed regions of and let . Since the ’s are disjoint, the ’s are disjoint sets. The number of connected components in is the same as the number of connected components in which by Lemma 4.2 is which is the same as by Lemma 4.3. Note that if does not pierce , then and therefore does not contribute to the sum above. Thus, if there are at most regions piercing , and for all , since the are disjoint, it follows that the number of connected components of is at most .
We now show that applying lens bypassing does not increase the fragmentation number of the arrangement.
Lemma 4.5 (bypassing a minimal lens does not increase fragmentation number).
Let be a set of regions and let be minimal lens (i.e., it does not contain another lens) in the arrangement. Let be the regions defining . If we now bypass by i.e., modify to , then for any region , .
Proof.
We first show that . To see this, note that the connected components of are the same as those in except that the lens we bypassed appears as a connected component in but not in . Since , does not contribute to the expression for in Lemma 4.3. Thus if the boundaries of and intersect, . Even if and don’t intersect, we have .
Now consider any region . The boundary of consists of an arc of and an arc of . Since is a minimal lens, the intersection of with consists of disjoint arcs each with one endpoint on and the other endpoint on . Since bypassing by effectively moves to 444more precisely a curve arbitrarily close to it, the intersection of with is combinatorially the same as the intersection of with . This implies that .
Definition 4.6 (Intersection graph of connected subgraphs and its depth).
Let be any graph and let be a set of connected subgraphs of . Then, the intersection graph defined by is a graph whose vertex set is and there is an edge corresponding to any pair iff and share a vertex of . We define the depth of the intersection graph as the maximum number of connected subgraphs in containing any specific vertex in .
The Lemma below follows from Frederickson’s decomposition [16].
Lemma 4.7.
Let and be fixed parameters. Let be the family of graphs such that each graph is the depth intersection graph of connected subgraphs in some planar graph s.t. the number of vertices in is at most times the number of vertices in . Then, admits a -decomposition where and .
Proof.
In order to show that admits a -decomposition, we need to show that for any given parameter and any graph with vertices, we can obtain and a partition of into satisfying the conditions stated in Definition 2.10.
We obtain the required by applying the result of Frederickson [16] which shows that the family of planar graphs admits a -decomposition for some constants . Let be the planar graph with s.t. is the depth intersection graph of connected subgraphs of . For any , let denote the connected subgraph in corresponding to the vertex in . Applying the result of Frederickson [16] with parameter , we obtain sets and a partition of into s.t. , , , , and .
Note that with . We next define as the set of all vertices in whose corresponding connected subgraph in contains one of the vertices in . Since any vertex in is contained in at most of the connected components, with . For any , we define as the set of vertices in whose corresponding connected subgraphs in only contain vertices from . Note that is a partition of since for any note that if does not contain any of the vertices in , then it can only contain vertices from exactly one group (this relies on the fact that for , ). Since every vertex in can be in at most connected subgraphs, . Next, note that with . To see this, note that for any s.t and , must contain a vertex of and each such vertex is contained in at most connected subgraphs. Finally, since for any and , and don’t have a vertex in common in – this follows from the fact that the vertices in are from and the vertices in are from and . The lemma follows.
We now define a two-colored version of an -piercing regions with additional properties. This will be useful in the analysis of the algorithms in Section 5.
Definition 4.8 (special -piercing regions).
A special -piercing set of regions is an -piercing set where each region is colored either red or blue, and any two regions of the same color are either disjoint or piercing, there are no monochromatic lenses i.e., lenses (not necessarily minimal) formed by pairs of regions of the same color and every bi-chromatic lens (i.e., a lens formed by a pair of regions of different color) is a minimal lens.
For any special -piercing set , let denote the bipartite intersection graph where there is an edge between any intersecting pair of regions of different colors in .
Lemma 4.9.
Let be an -piercing set of regions in which each region is colored either red or blue and let be a set of points s.t. and are in general position. Then, we can construct a special -piercing set in which corresponding to every , there is a region of the same color as so that the following condition holds for any . If is contained in at least one red and one blue region in , it is also contained in at least one red and one blue region in .
Proof.
We will do a sequence of minimal lens bypassing operations so that for any point the following coverage condition holds: if is contained in at least one red and one blue region before the bypassing, it is still contained in at least one red and one blue region after the bypassing. After such simplifications, we will show that the resulting arrangement is a special -piercing set. Note that in any lens bypassing any region is either not modified or modified to a region . Also, note that by Lemma 4.5, the set of regions remains -piercing at all times. First, note that upon bypassing any minimal lens formed by two regions of the same color, regardless of the region modified, the coverage condition holds. Therefore we can bypass all such monochromatic minimal lenses. Now consider a bi-chromatic minimal lens formed by a red and a blue region which is contained in a third region (red or blue). We can also bypass such lenses while maintaining the coverage condition by making sure that the region modified during the bypassing is the one whose color matches the third region. After we have bypassed all minimal lenses of the above two types, observe that we cannot have any monochromatic lenses (even non-minimal). To see this, consider a minimal monochromatic lens i.e., a monochromatic lens which does not contain another monochromatic lens. Such a lens is either a minimal lens or it contains a bi-chromatic lens. In the former case, the lens should have been bypassed, and in the latter case, the bi-chromatic lens is contained in a third region and should have been bypassed. This implies that any regions of the same color are either disjoint or are piercing. The lemma follows.
Definition 4.10 (simplified special -piercing set).
A simplified special -piercing set of regions is a special -piercing set s.t. any two regions in that pierce each other don’t form any lenses, and any two regions in that don’t pierce intersect in at most one lens.
Lemma 4.11 (obtaining a simplified special -piercing set).
Given any special -piercing set , we can obtain a simplified set in which for each region , there is a region of the same color as and the following condition holds. For any two regions of different colors intersect iff and intersect.
Proof.
Since is a special -piercing set, the arrangement of regions in does not contain any monochromatic lenses and all bi-chromatic lenses are minimal. We repeatedly bypass bi-chromatic minimal lenses formed by regions that also intersect outside the lens as this does not violate the conditions of the lemma. Thus if any lens is left it is the only cell in the intersection of two regions forming it. In particular, the two regions cannot be piercing and they cannot be intersecting in any other lens.
Lemma 4.12.
Let be a simplified special -set of regions. The number of cells in the arrangement of the regions in is where .
Proof.
It suffices to show that the number of intersections of the boundaries of pairs of regions in is . In order to do this, we partition the red regions into sets s.t. the regions in each set are disjoint. This is because regions of the same color are either piercing or disjoint and no region is pierced by more than other regions. Similarly, we partition the blue regions into sets s.t. the regions in each set are disjoint. Overall, we have sets of disjoint regions. Consider any two of these families and (which may or may not be of the same color).
We next argue that the number of intersections among the regions in is where . To see this, imagine replacing each region by the connected components in . The number of components any is split into is at most by Lemma 4.4. The overall number of regions now is at most which are interior disjoint but certain pairs touch along their boundaries.
Consider the touching graph defined by these regions in which there is a vertex corresponding to each region and two vertices are adjacent iff the boundaries of their corresponding regions overlap in an arc of non-zero length. Note that this graph is planar and hence, the number of edges in this graph is at most . Here it is important to note also that the touching graph does not have multi-edges due to the assumption that two piercing regions do not form any lenses and two non-piercing regions form at most one lens. Each edge in the touching graph corresponds to two intersection points among boundaries of regions in , showing that the number of intersections among the regions in is at most .
Since each set appears in pairs of the form and , the overall number of intersections among the boundaries of regions in is which in turn implies that the number of cells in the arrangement of regions in is also .
Lemma 4.13.
For any fixed constants , the family of bipartite intersection graphs admits a -decomposition where , and .
Proof.
Consider a special -set . By Lemma 4.11, we can assume without loss of generality, that is a simplified special -set. By Lemma 4.12, the number of cells in the arrangement of the regions in is Now consider the dual arrangement graph (Definition 2.4) of the regions in . Note that is planar and the cells contained in any region correspond to a connected subgraph of . Note also that any cell is contained in at most red regions and at most blue regions and therefore is in at most regions overall. In other words, each graph is the depth intersection graph of connected subgraphs in a planar graph having at most vertices where . By Lemma 4.7, admits a -decomposition where , and .
5 Discrete Independent Set and Set Cover
Let be a set of -piercing regions and let be a set of points in the plane.
In this section, we define the local search algorithms and show that they yield a PTAS for the Discrete Independent Set and Set Cover problems.
-Local Search
For any given parameter , the algorithm is the following. We always maintain a feasible solution which in the beginning is chosen arbitrarily. At any stage, we check if it is possible to replace some subset of the regions by another subset of regions so that and is a feasible solution with a better objective value than . If this is possible, we replace by and continue. If such a local improvement is not possible, we say that the current solution is -locally optimal in which case, we simply return the current solution . A straightforward implementation of the -local search algorithm has running time where is the total number of regions.
Theorem 3.1 (Discrete Independent Set). [Restated, see original statement.]
The -local search algorithm with for an appropriately large and yields a -approximation for the Discrete Independent Set problem defined by a set of -piercing regions and a set of points . The algorithm runs in time .
Proof.
Let Opt be an optimal solution and Local be the solution returned by the -local search algorithm. For the purpose of the analysis, it helps to assume that Opt and Local are disjoint. This can be done by removing the regions common to both and then comparing their sizes. Formally, let and let and . Let be the set of points not covered by the regions in . Note that and are disjoint and are the optimal and -locally optimal solutions respectively with respect to . For the following arguments, we will color the regions in red and the regions in blue.
By applying Lemma 4.9, to the set of regions and the set of points , we obtain a special -set where for each region , there is a corresponding shrunk region in . Let be a bipartite graph in which there is an edge for some and iff and (the shrunk regions corresponding to and ) intersect. Note that is isomorphic to the bipartite intersection graph of the regions in .
Note that every point is covered by at most one red region in and one blue region in and furthermore if is contained in a red region as well as a blue region , then there an edge in . Thus, the following local search condition holds: is a feasible solution for any where denotes the set of blue neighbors of the regions in in .
By Lemma 4.13, has a -decomposition where , and i.e., there is a set with and a partition of into sets, where , s.t. for all , for all , and . For , let and . We now show that for any , . For contradiction suppose that we have for some . By the local search condition: is a feasible solution. Since , is an improved feasible solution obtained from by adding elements and removing elements contradicting the assumption that is -locally optimal. Now, where from which it follows that for sufficiently large . This implies that . The theorem follows.
The proof of the following theorem is very similar to the proof of Theorem 3.1.
Theorem 3.2 (Set Cover). [Restated, see original statement.]
The -local search algorithm with for an appropriately large and yields a -approximation for the Set Cover problem defined by a set of -piercing regions and a set of points . The algorithm runs in time .
Proof.
The proof is very similar to the proof of Theorem 3.1 and we use the notation defined there. Since every point is covered by at least one red region and one blue region, the following local search condition holds: is a feasible solution for any . As in the proof of Theorem 3.1, we apply Lemma 4.13 and obtain the sets , and . We next show that for any , . For contradiction, suppose that this is not the case for some . By the local search condition, is a feasible solution. Since , is an improved feasible solution obtained by removing elements from and adding elements, contradicting the assumption that is -locally optimal. Now where . Assuming is sufficiently large, we obtain which implies that . The theorem follows.
Now, we show that the dual set system defined by a set of points and an -piercing set of regions has bounded shallow-cell complexity.
Theorem 3.3. [Restated, see original statement.]
Let be an -piercing set of regions and be a set of points. For any point , let denote the subset of regions in containing . Then, the set system (often called the dual set system defined by and ) has shallow-cell complexity .
Proof.
We first prove that if is a special -set then the number of combinatorially distinct bi-chromatic cells of depth (i.e., cells contained in exactly one red and one blue region) in the arrangement is . We then apply a standard argument due Clarkson and Shor [10] to prove the lemma. Let us first assume that is a special -set. Since we are only interested in combinatorially distinct bi-chromatic cells, we may assume by Lemma 4.11 that is a simplified special -set and by Lemma 4.12 the number of such cells is . Now, assume that is any -piercing set. Our goal is to bound the number of combinatorially distinct cells of depth in the arrangement . Since every cells has a vertex, it suffices to count the vertices of depth contained in a combinatorially distinct set of regions. Let denote this number. Let be set obtained by picking each region in independently with probability and coloring it red or blue uniformly at random. The probability that a vertex of depth appears as a bi-chromatic vertex (a vertex formed by intersection of the boundaries of a red and a blue region but not contained in any other region) in the arrangement is at least . Hence, in expectation, we get at least such vertices. By Lemma 4.9, the set can be modified to a special -piercing set with the same set of combinatorially distinct bi-chromatic cells. This means that number of bi-chromatic vertices in the arrangement is at most which is in expectation. Thus, implying that is .
6 Conclusion
We have introduced piercing complexity as a natural topological parameter for geometric optimization and shown that bounded piercing enables PTASs for fundamental packing and covering problems. Specifically, for -piercing regions, we obtained PTASs for unweighted Set Cover and Discrete Independent Set problems, along with -approximations for their weighted variants. These results generalize known algorithms for non-piercing regions and yield improved guarantees for important special cases, including axis-parallel rectangles and line segments under bounded piercing.
Our results build on the local-search framework developed for non-piercing regions [5, 35] but extend it to a significantly broader class of instances. Several important questions remain open. First, does a PTAS exist for the Hitting Set problem with -piercing regions? This appears non-trivial and likely requires deeper structural insights. Second, while known hardness constructions show that unbounded piercing can be an obstruction to PTASs, it remains unclear whether the fragmentation number is required to be bounded, or whether a weaker condition would suffice.
More broadly, piercing complexity captures a fundamental source of hardness that existing parameters – VC-dimension, shallow-cell complexity, fatness, and density – fail to measure. These parameters control combinatorial richness or geometric overlap but do not capture how intersections topologically fragment regions. By quantifying this fragmentation, piercing complexity directly measures a structural property that obstructs local-search techniques. In this sense, -piercing gives a broad structural extension of the non-piercing regime under which PTASs remain attainable, providing a unifying framework that explains both classical hardness constructions and previously isolated positive results. Finally, a finer structural understanding of piercing interactions may prove useful in resolving longstanding open problems, such as the Set Cover and Independent Set problems for axis-parallel rectangles, and the problems with line segments and points in the plane.
References
- [1] Anna Adamaszek and Andreas Wiese. A QPTAS for maximum weight independent set of polygons with polylogarithmically many vertices. In Proceedings of the Twenty-Fifth Annual ACM-SIAM Symposium on Discrete Algorithms, Proceedings of the Twenty-Ninth Annual ACM-SIAM Symposium on Discrete Algorithms, SODA, pages 645–656, 2014. doi:10.1137/1.9781611973402.49.
- [2] Pankaj K. Agarwal, Marc van Kreveld, and Subhash Suri. Label placement by maximum independent set in rectangles. Computational Geometry, 11(3):209–218, 1998. doi:10.1016/S0925-7721(98)00028-5.
- [3] Rom Aschner, Matthew J. Katz, Gila Morgenstern, and Yelena Yuditsky. Approximation schemes for covering and packing. In Subir Kumar Ghosh and Takeshi Tokuyama, editors, WALCOM: Algorithms and Computation, 7th International Workshop, WALCOM 2013, Kharagpur, India, February 14-16, 2013. Proceedings, volume 7748 of Lecture Notes in Computer Science, pages 89–100. Springer, 2013. doi:10.1007/978-3-642-36065-7_10.
- [4] Jozsef Balogh and Jozsef Solymosi. On the number of points in general position in the plane. Discrete Analysis, pages 1–20, 2018. doi:10.19086/da.4438.
- [5] Aniket Basu Roy, Sathish Govindarajan, Rajiv Raman, and Saurabh Ray. Packing and covering with non-piercing regions. Discrete & Computational Geometry, 2018.
- [6] Hervé Brönnimann and Michael T. Goodrich. Almost optimal set covers in finite VC-dimension. Discrete & Computational Geometry, 14(4):463–479, 1995. doi:10.1007/BF02570718.
- [7] Timothy M. Chan and Elyot Grant. Exact algorithms and APX-hardness results for geometric packing and covering problems. Computational Geometry Theory and Applications, 47(2):112–124, 2014. doi:10.1016/j.comgeo.2012.04.001.
- [8] Timothy M. Chan, Elyot Grant, Jochen Könemann, and Malcolm Sharpe. Weighted capacitated, priority, and geometric set cover via improved quasi-uniform sampling. In Proceedings of the Twenty-third Annual ACM-SIAM Symposium on Discrete Algorithms, Proceedings of the Twenty-Ninth Annual ACM-SIAM Symposium on Discrete Algorithms, SODA, pages 1576–1585, 2012. doi:10.1137/1.9781611973099.125.
- [9] Timothy M. Chan and Sariel Har-Peled. Approximation algorithms for maximum independent set of pseudo-disks. Discrete and Computational Geometry, 48(2):373–392, 2012. doi:10.1007/S00454-012-9417-5.
- [10] Kenneth L. Clarkson and Peter W. Shor. Applications of random sampling in computational geometry, II. Discrete & Computational Geometry, 4(5):387–421, 1989. doi:10.1007/BF02187740.
- [11] Kenneth L. Clarkson and Kasturi R. Varadarajan. Improved approximation algorithms for geometric set cover. Discrete & Computational Geometry, 37(1):43–58, 2007. doi:10.1007/s00454-006-1273-8.
- [12] Marek Cygan, Fedor V. Fomin, Łukasz Kowalik, Daniel Lokshtanov, Dániel Marx, Marcin Pilipczuk, Michał Pilipczuk, and Saket Saurabh. Parameterized algorithms, volume 5(4). Springer, 2015. doi:10.1007/978-3-319-21275-3.
- [13] Thomas Erlebach and Erik Jan van Leeuwen. PTAS for weighted set cover on unit squares. In Approximation, Randomization, and Combinatorial Optimization. Algorithms and Techniques, 13th International Workshop, APPROX 2010, and 14th International Workshop, RANDOM 2010, Barcelona, Spain, September 1-3, 2010. Proceedings, pages 166–177, 2010. doi:10.1007/978-3-642-15369-3_13.
- [14] Guy Even, Dror Rawitz, and Shimon Shahar. Hitting sets when the vc-dimension is small. Inf. Process. Lett., 95(2):358–362, 2005. doi:10.1016/j.ipl.2005.03.010.
- [15] Uriel Feige. A threshold of ln n for approximating set cover. Journal of the ACM (JACM), 45(4):634–652, 1998. doi:10.1145/285055.285059.
- [16] Greg N Frederickson. Fast algorithms for shortest paths in planar graphs, with applications. SIAM Journal on Computing, 16(6):1004–1022, 1987. doi:10.1137/0216064.
- [17] Waldo Gálvez, Arindam Khan, Mathieu Mari, Tobias Mömke, Madhusudhan Reddy Pittu, and Andreas Wiese. A 3-approximation algorithm for maximum independent set of rectangles. In Joseph (Seffi) Naor and Niv Buchbinder, editors, Proceedings of the 2022 ACM-SIAM Symposium on Discrete Algorithms, SODA, Virtual Conference / Alexandria, VA, USA, January 9 - 12, 2022, pages 894–905. SIAM, 2022. doi:10.1137/1.9781611977073.38.
- [18] Matt Gibson and Imran A. Pirwani. Algorithms for dominating set in disk graphs: Breaking the logn barrier. In Mark de Berg and Ulrich Meyer, editors, Algorithms – ESA 2010, pages 243–254, Berlin, Heidelberg, 2010. Springer Berlin Heidelberg.
- [19] Daniel Golovin, Anupam Gupta, Amit Kumar, and Kanat Tangwongsan. All-norms and all-l_p-norms approximation algorithms. In IARCS Annual Conference on Foundations of Software Technology and Theoretical Computer Science (2008), pages 199–210. Schloss Dagstuhl – Leibniz-Zentrum für Informatik, 2008. doi:10.4230/LIPIcs.FSTTCS.2008.1753.
- [20] Sariel Har-Peled. Geometric approximation algorithms. Number 173 in Mathematical Surveys and Monographs. American Mathematical Soc., 2011.
- [21] Sariel Har-Peled and Kent Quanrud. Approximation algorithms for polynomial-expansion and low-density graphs. In Algorithms – ESA 2015 - 23rd Annual European Symposium, Patras, Greece, September 14–16, 2015, Proceedings, pages 717–728, 2015. doi:10.1007/978-3-662-48350-3_60.
- [22] Sariel Har-Peled and Kent Quanrud. Approximation algorithms for polynomial-expansion and low-density graphs. SIAM J. Comput., 46(6):1712–1744, 2017. doi:10.1137/16M1079336.
- [23] David Haussler and Emo Welzl. epsilon-nets and simplex range queries. Discrete & Computational Geometry, 2:127–151, 1987. doi:10.1007/BF02187876.
- [24] Dorit S. Hochbaum and Wolfgang 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.
- [25] J Håstad. Clique is hard to approximate within . Acta Mathematica, 182(1):105–142, 1999.
- [26] Bruno Jartoux and Nabil H. Mustafa. Optimality of Geometric Local Search. In 34th International Symposium on Computational Geometry (SoCG 2018), Leibniz International Proceedings in Informatics (LIPIcs), pages 48:1–48:15, 2018. doi:10.4230/LIPIcs.SOCG.2018.48.
- [27] Katarzyna Kowalska and Michał Pilipczuk. Parameterized and approximation algorithms for coverings points with segments in the plane. In 41st International Symposium on Theoretical Aspects of Computer Science, 2024.
- [28] Daniel Lokshtanov, Fahad Panolan, Saket Saurabh, Jie Xue, and Meirav Zehavi. Subexponential parameterized algorithms on disk graphs. In Proceedings of the 2022 Annual ACM-SIAM Symposium on Discrete Algorithms (SODA), pages 2005–2031. SIAM, 2022.
- [29] Laszlo Lovász. On the ratio of optimal integral and fractional covers. Discrete Math, 13:383–390, 1975. doi:10.1016/0012-365X(75)90058-8.
- [30] Jiří Matoušek. Lectures on Discrete Geometry. Springer-Verlag New York, Inc., Secaucus, NJ, USA, 2002.
- [31] Joseph S. B. Mitchell. Approximating maximum independent set for rectangles in the plane. In 62nd IEEE Annual Symposium on Foundations of Computer Science, FOCS 2021, Denver, CO, USA, February 7-10, 2022, pages 339–350. IEEE, 2021. doi:10.1109/FOCS52979.2021.00042.
- [32] Nabil H. Mustafa. Sampling in combinatorial and geometric set systems, volume 265. American Mathematical Society, 2022.
- [33] Nabil H. Mustafa and Saurabh Ray. Improved results on geometric hitting set problems. Discrete & Computational Geometry, 44(4):883–895, 2010. doi:10.1007/S00454-010-9285-9.
- [34] János Pach and Gábor Tardos. Tight lower bounds for the size of epsilon-nets. In Proceedings of the twenty-seventh annual symposium on Computational geometry, pages 458–463, 2011. doi:10.1145/1998196.1998271.
- [35] Rajiv Raman and Saurabh Ray. Constructing planar support for non-piercing regions. Discrete and Computational Geometry, 64(3):1098–1122, 2020. doi:10.1007/S00454-020-00216-W.
- [36] Rajiv Raman and Saurabh Ray. On the geometric set multicover problem. Discret. Comput. Geom., 68(2):566–591, 2022. doi:10.1007/s00454-022-00402-y.
- [37] Rajiv Raman and Karamjeet Singh. Supports for outerplanar and bounded treewidth graphs, 2025. doi:10.48550/arXiv.2504.05039.
- [38] Kasturi Varadarajan. Weighted geometric set cover via quasi-uniform sampling. In Proceedings of the 42nd ACM Symposium on Theory of Computing, pages 641–648, 2010. doi:10.1145/1806689.1806777.
- [39] Bang Wang. Coverage problems in sensor networks: A survey. ACM Computing Surveys (CSUR), 43(4):1–53, 2011. doi:10.1145/1978802.1978811.
