11email: {shinwooan,kyungjincho,eunjin.oh}@postech.ac.kr
Faster Algorithms for Cycle Hitting Problems on Disk Graphs††thanks: This work was supported by the National Research Foundation of Korea(NRF) grant funded by the Korea government(MSIT) (No.RS-2023-00209069)
Abstract
In this paper, we consider three hitting problems on a disk intersection graph: Triangle Hitting Set, Feedback Vertex Set, and Odd Cycle Transversal. Given a disk intersection graph , our goal is to compute a set of vertices hitting all triangles, all cycles, or all odd cycles, respectively. Our algorithms run in time , , and , respectively, where denotes the number of vertices of . These do not require a geometric representation of a disk graph. If a geometric representation of a disk graph is given as input, we can solve these problems more efficiently. In this way, we improve the algorithms for those three problem by Lokshtanov et al. [SODA 2022].
Keywords:
Disk graphs, feedback vertex set, triangle hitting set1 Introduction
In this paper, we present subexponential parameterized algorithms for the following three well-known parameterized problems on disk graphs: Triangle Hitting Set, Feedback Vertex Set, and Odd Cycle Transversal. Given a graph and an integer , these problems ask for finding a subset of of size that hits all triangles, cycles, and odd-cycles, respectively. On general graphs, the best-known algorithms for Triangle Hitting Set, Feedback Vertex Set, and Odd Cycle Transversal run in , , and time, respectively, where denotes the number of vertices of a graph [15, 16, 20]. All these problems are NP-complete, and moreover, no algorithm exists for these problems on general graphs unless ETH fails [10]. This motivates the study of these problems on special graph classes such as planar graphs and geometric intersection graphs.
Although these problems are NP-complete even for planar graphs, much faster algorithms exist for planar graphs. More specifically, they can be solved in time , , and , respectively, on planar graphs [11, 18]. Moreover, they are optimal unless ETH fails. The -time algorithms for Triangle Hitting Set and Feedback Vertex Set follow from the bidimensionality theory of Demaine at al [11]. A planar graph either has an grid graph as a minor, or its treewidth111The definition of the treewidth can be found in Section 2 is . This implies that for a YES-instance of Feedback Vertex Set, the treewidth of is . Then a standard dynamic programming approach gives a -time algorithm for Feedback Vertex Set. This approach can be generalized for -minor free graphs and bounded-genus graphs.
Recently, several subexponential-time algorithms for cycle hitting problems have been studied on geometric intersection graphs [2, 5, 6, 12, 13, 14]. Let be a set of geometric objects such as disks or polygons. The intersection graph of is the graph where a vertex of corresponds to an element of , and two vertices of are connected by an edge in if and only if their corresponding elements in intersect. In particular, the intersection graph of (unit) disks is called a (unit) disk graph, and the intersection graph of interior-disjoint polygons is called a map graph. Unlike planar graphs, map graphs and unit disk graphs can have large cliques. For unit disk graphs, Feedback Vertex Set can be solved in time [2], and Odd Cycle Transversal can be solved in time [4]. For map graphs, Feedback Vertex Set can be solved in time [12]. These are almost tight in the sense that no -time algorithm exists unless ETH fails [12, 14].
However, until very recently, little has been known for disk graphs, a broad class of graphs that generalizes planar graphs and unit disk graphs. Very recently, Lokshtanov et al. [17] presented the first subexponential-time algorithms for cycle hitting problems on disks graphs. The explicit running times of the algorithms are summarized in the first row of Table 1. On the other hand, the best-known lower bound on the computation time for these problems is assuming ETH. There is a huge gap between the upper and lower bounds.
1.0.1 Our results.
In this paper, we make progress on the study of subexponential-time parameterized algorithms on disk graphs by presenting faster algorithms for the cycle hitting problems. We say an algorithm is robust if this algorithm does not requires a geometric representation of a disk graph. We present -time, -time, and -time robust algorithms for Triangle Hitting Set, Feedback Vertex Set, and Odd Cycle Transversal, respectively. Furthermore, we present -time, -time, and -time algorithms for Triangle Hitting Set, Feedback Vertex Set, and Odd Cycle Transversal which are not robust. These results are summarized in the second and third rows of Table 1.
For Triangle Hitting Set, we devised a kernelization algorithm based on a crown decomposition by modifying the the algorithm in [1]. After a branching process, we obtain a set of instances such that has a set of triangles, which we call a core, for a value . Every triangle of shares at least two vertices with at least one triangle of a core. This allows us to remove several vertices from so that the number of vertices of is . Then we can obtain a tree decomposition of small treewidth, and then we apply dynamic programming on the tree decomposition.
For Feedback Vertex Set and Odd Cycle Transversal, we give an improved analysis of the algorithms in [17] by presenting improved bounds on one of the two main combinatorial results presented in [17] concerning the treewidth of disk graphs (Theorem 2.1). More specifically, Theorem 1 of [17] says that for any subset of such that for any two vertices and not in , the size of is , where is the ply of the disks represented by the vertices of . We improve an improved bound of . To obtain an improved bound, we classify the vertices of into three classes in a more sophisticated way, and then use the concept of the additively weighted higher-order Voronoi diagram. We believe the additively weighted higher-order Voronoi diagram will be useful in designing optimal algorithms for Triangle Hitting Set, Feedback Vertex Set, and Odd Cycle Transversal. Indeed, the Voronoi diagram (or Delaunay triangulation) is a main tool for obtaining an optimal algorithm for Feedback Vertex Set on unit disk graphs [2].
2 Preliminaries
For a graph , we let and be the sets of vertices and edges of , respectively. For a subset of , we use to denote the subgraph of induced by . Also, for a subset of , we simply denote by . For a vertex of , let be the set of the vertices of adjacent to . We call it the neighborhood of .
A triangle of is a cycle of consisting of three vertices. We denote a triangle consisting of three vertices , and by . We sometimes consider it as the set of vertices. For instance, the union of triangles is the set of all vertices of the triangles. A subset of is called a triangle hitting set of if has no triangle. Also, is called a feedback vertex set if has no cycle (and thus it is forest). Finally, is called a odd cycle transversal if has no odd cycle. In other words, a triangle hitting set, a feedback vertex set, and a odd cycle transversal hit all triangles, cycles and odd cycles, respectively. Notice that a feedback vertex set of is also a triangle hitting set.
2.0.1 Disk graphs.
Let be a disk graph defined by a set of disks. In this case, we say that is a geometric representation of . For a vertex of , we let denote the disk of represented by . The ply of is defined as the maximum number of disks of containing a common point. Note that the disk graph of a set of disks of ply has a clique of size at least . If it is clear from the context, we say that the ply of is . The arrangement of is the subdivision of the plane formed by the boundaries of the disks of that consists of vertices, edges and faces. The arrangement graph of , denoted by , is the plane graph where every face of the arrangement of contained in at least one disk is represented by a vertex, and vertices are adjacent if the faces they represent share an edge. Given a disk graph , it is NP-hard to compute its geometric representation [9]. We say an algorithm is robust if this algorithm does not requires a geometric representation of a disk graph.
2.0.2 Tree decomposition.
A tree decomposition of an undirected graph is defined as a pair , where is a tree and is a mapping from nodes of to subsets of (called bags) with the following conditions. Let be the set of bags of .
-
•
For , there is at least one bag in which contains .
-
•
For , there is at least one bag in which contains both and .
-
•
For , the nodes of containing in their bags are connected in .
The width of a tree decomposition is defined as the size of its largest bag minus one, and the treewidth of is the minimum width of a tree decomposition of . The treewidth of a disk graph is , where is the ply of [17].
2.0.3 Weighted Treewidth of a Disk Graph
One of the ingredients of our algorithms is to partition the vertices in each bag of a tree decomposition of into cliques. This technique was already used in [5] to design ETH-tight (non-parameterized) algorithms for various problems on geometric intersection graphs. To use this observation, they introduced the notion of clique-weighted treewidth.
We fix a geometric representation of . For a subset of , we say the partition of is a clique partition if all disks represented by the vertices of contain a common point in the plane for each index . We define the clique-weight of as the minimum of over all clique partitions of . Since a triangle hitting set contains all except for at most two vertices of , the number of intersections between and a triangle hitting set is , where is the clique-weight of .
A balanced separator of is a subset of vertices such that each connected component of has complexity for a constant . De berg et al. [5] showed that a geometric intersection graph of fat objects has a balanced separator of small clique-weight.
Lemma 1 ([5]).
A disk graph has a balanced separator and a clique partition of with weight . Moreover, if a geometric representation of is given, an and its clique partition can be computed in polynomial time.
For a tree decomposition of a graph, we define the clique-weighted width of as the maximum weight of a bag of the decomposition. Then the clique-weighted treewidth of a graph is defined as the minimum clique-weighted width of a tree decomposition of the graph. In the following, we simply use the term weight (and weighted width) instead of clique-weight (and clique-weighted width), if it is clear from the context. Then Lemma 1 implies the following lemma.
Lemma 2.
The weighted treewidth of a disk graph is . Given a geometric representation of , we can compute a tree decomposition of with weighted treewidth in polynomial time.
Proof. We construct a tree decomposition recursively for a disk graph with its geometric representation. Our construction satisfies that a subtree of rooted at a node together with the bags corresponding its nodes forms a tree decomposition of a subgraph of , say .
Our construction starts from the root node of with . For a current node of , we set if the number of vertices of is . Otherwise, we find a balanced separator of of size and its clique partition with weight in polynomial time by Lemma 1. For each connected component of , we recursively construct a tree decomposition of . Then we connect the root of to so that it becomes a child of . In addition to this, we add the vertices of to the bag of every node in the subtree rooted at including itself. There is no edge between two different components and . Thus, is a tree decomposition of . Furthermore, its weight treewidth is because is a balanced separator.
2.0.4 Bounded Ply.
If the ply of a geometric representation of is bounded, we can obtain a better bound on the weighted treewidth of .
Lemma 3.
The weighted treewidth of is at most times the treewidth of . Furthermore, the treewidth of is at most times the treewidth of .
Proof. Let be a tree decomposition of . We can construct a tree decomposition of as follows. Each vertex of corresponds to a face of the arrangement of . Moreover, the disks of containing a common face of the arrangement forms a clique in . For each node , let be the set of the vertices corresponding to disks of containing faces corresponding to the vertices of . It is not difficult to see that is a tree decomposition of . Then the weighted width (and treewidth) of is (and ) times the width of . This implies that the weighted treewidth (and treewidth) of is at most (and ) times the treewidth of .
Given a tree decomposition of with weighted width , Feedback Vertex Set can be solved in time using a rank-based approach as observed in [5].222De berg et al. [5] showed that Feedback Vertex Set and several problems can be solved in time for similarly sized objects. Although the other algorithms cannot be extended to an intersection graph of fat objects, the algorithm for Feedback Vertex Set can be extended to an intersection graph of fat object. We briefly show how to solve Feedback Vertex Set in Lemma 17. Similarly, Triangle Hitting Set and Odd Cycle Transversal can be solved in and time, respectively. They can be obtained by slightly modifying standard dynamic programming algorithms for those problems as observed in [5]. We briefly show how to do this in Section 5.
2.0.5 Higher-Order Voronoi Diagram.
To analyze the treewidth of of bounded ply, we use the concept of the higher-order Voronoi diagram. Given a set of weighted sites (points), its order- Voronoi diagram is defined as the subdivision of into maximal regions such that all points within a given region have the same nearest sites. Here, the distance between a site with weight and a point in the plane is defined as , where is the Euclidean distance between and .
Lemma 4 (Theorem 4 in [19]).
The complexity of the addictively weighted order- Voronoi diagram of point sites in the plane is .
3 Triangle Hitting Set
In this section, we present a robust -time algorithm for Triangle Hitting Set, which improves the -time algorithm of [17].
3.1 Two-Step Branching Process
We first apply a branching process as follows to obtain instances one of which is a YES-instance of Triangle Hitting Set where every clique of has size . For a clique of and a triangle hitting set , all except for at most two vertices of are contained in . In particular, if has a triangle hitting set of size , any clique of has size at most . For a clique of size at least , we branch on which vertices of the clique are not included in a triangle hitting set. After this, the solution size decrease by at least . We repeat this until every clique has size . In the resulting branching tree, every node has at most children since any clique of has size at most . And the branching tree has height . In this way, we can obtain instances of Triangle Hitting Set one of which is a YES-instance of Triangle Hitting Set where has a geometric representation of ply at most . Moreover, is an induced subgraph of , and . Using the EPTAS for computing a maximum clique in a disk graph [8], we can complete the branching step in time. Note that the algorithm in [8] does not require a geometric representation of a graph.

For each instance obtained from the first branching, we apply the second branching process to obtain instances one of which is a YES-instance having a core of size . A set of triangles of is called a core if for a triangle of , a triangle of shares at least two vertices with the triangle. See Fig. 1(a). During the branching process, we mark an edge to remember that one of its endpoints must be added to a triangle hitting set. The marking process has an invariant that no two marked edges share a common endpoint, and thus the number of marked edges is at most . The marks will be considered in the dynamic programming procedure. Initially, all edges are unmarked.
Let be a vertex such that has a matching of size at least , where be the set of neighbors of not incident to any marked edge. See Fig. 1(b). In this case, a triangle hitting set contains either or at least one endpoint of each edge in the matching. We branch on whether or not is added to a triangle hitting set. For the first case, we remove and its adjacent edges, and decrease by one. For the second case, we know that at least one endpoint of each edge in the matching must be contained in a triangle hitting set. However, we do not make a decision at this point. Instead, we simply mark all such edges.
Lemma 5.
The total number of instances from the two-step branching process is . Moreover, the branching process runs in time.
Proof. Recall that the first branching step produces instances. Let be an instance of Triangle Hitting Set we produce during the second branching process. Let be the number of instances produced by where has marked edges. By construction, we produce two instances: one has parameter , and one has at least marked edges. Thus we have
Since , we have . Therefore, the total number of instances obtained from the two-step branching process is . Moreover, since our algorithm runs in polynomial time per instance, the total running time is .
This branching process was already used in [17]. The following lemma is a key for our improvement over [17].
Lemma 6.
Let be a YES-instance obtained from the two-step branching. Then has a core of size , and we can compute one in polynomial time.
Proof. We construct a core of as follows. Let be the the union of a triangle hitting set of size at most and the set of all endpoints of the marked edges, which can be computed in polynomial time.333 We can find a hitting set of size at most as follows: Find a triangle, and add all its vertices to a triangle hitting set. Then remove all its vertices from the graph. Note that, the size of is at most . Then let be the set of triangles constructed as follows: for each edge of , we add an arbitrary triangle of formed by and to , if it exists. The number of triangles in is by Lemma 7.
At this moment, is not necessarily a core. See Fig. 1(c). Thus we add several triangles further to to compute a core of . By the branching process, for every vertex of , has a maximum matching of size at most . We add the triangles formed by and the edges of the maximum matching to for every vertex of . Note that we do not update during this phase, and thus triangles added to due to two different vertices might intersect. This algorithm clearly runs in polynomial time. Moreover, since the size of is , we add at most triangles to , and thus the size of is .
We claim that is a core of . Let be a triangle of not in . Since contains a triangle hitting set of size at most , it must contain at least one of and . If at least two of them, say and , are contained in , then there exists a triangle having edge in by construction. Thus, shares two vertices with such a triangle. The remaining case is that exactly one of them, say , is contained in . In this case, we have considered and a maximum matching of . The only reason why is not added to is that the edge is adjacent to another edge for some other triangle . Then is added to . Note that and share at least two vertices, and thus satisfies the condition for being a core.
Lemma 7.
For a subset of of size , has edges.
Proof. We use a charging scheme. For each edge of , we charge it to the smaller disk among the disks represented by and . Then each disk is charged by edges. This is because is intersected by disks larger than . To see this, recall that the ply of is at most due to the first branching process. Since the size of is , the number of edges of is .
The following observation will be used in the correctness proof of the kernelization step in Section 3.2. The observation holds because does not have any triangle not appearing in . {observation} Let be an induced subgraph of such that contains all vertices of a core of . Then is a core of .
For each instance obtained from the branching process, we apply the cleaning step that removes all vertices not hitting any triangle of from . Specifically, we remove a vertex of degree less than two. Also, we remove a vertex whose neighbors are independent. Whenever we remove a vertex, we also remove its adjacent edges.
3.2 Kernelization Using Crown Decomposition
Let be a YES-instance we obtained from the branching process. We show that if the number of vertices of contained in the triangles of is at least , we can construct a crown for the triangles of . Then using this, we can produce a YES-instance of Triangle Hitting Set in polynomial time where is a proper induced subgraph of . Thus by repeatedly applying this process (at most times) and then by removing all vertices not contained in any triangle of , we can obtain a YES-instance where is a disk graph of complexity .

We first define a crown decomposition of a graph for triangles. For illustration, see Fig. 2. A crown decomposition of a graph was initially introduced to construct a linear kernel for Vertex Cover. Later, Abu-Khzam [1] generalized this concept to hypergraphs. Using this, he showed that a triangle hitting set admits a quadratic kernel for a general graph. In our case, we will show that the size of a kernel is indeed due to the branching process. More specifically, it is due to the existence of a core of size .
Definition 1 ([1]).
A crown for the triangles of is a triple with a subset of , a subset of , and a matching of s.t.
-
•
no two vertices of are contained in the same triangle of ,
-
•
each edge of forms a triangle with some vertex of , and
-
•
every edge (vertex in ) of is matched under ,
where is the bipartite graph with vertex set such that and are connected by an edge if and only if and form a triangle.
If a crown for the triangles of exists, we can remove all vertices in , but instead, we mark all edges of . Then the resulting graph also has a triangle hitting set of size at most [1].
Lemma 8 ([1, Lemma 2]).
Let be a YES-instance of Triangle Hitting Set, and is a crown for the triangles of . Then the subgraph of obtained by removing all vertices of and by marking all edges of has a trianlge hitting set of size at most .
Proof. To make the paper self-contained, we give a proof of this lemma, but one can find this proof also in [1]. Let be a triangle hitting set of of size . We construct another hitting set from not containing any vertex of as follows. We initialize as . If there is an edge of and , then must be contained in , where the edge of is in . Then we replace with . We repeat this until at least one vertex of all edges of is contained in . Then .
We claim that is also a triangle hitting set. A triangle not intersecting is hit by by construction. Thus consider a triangle intersecting . Since no two vertices of are contained in the same triangle, exactly one of and , say , is contained in . Since every edge of is matched under , . Since contains either or , is hit by . Therefore, is a triangle hitting set of size at most .
Now we show that has a crown for its triangles if vertices are contained in the triangles of , where is a sufficiently large constant. Let be a core of of size , which exists due to Lemma 6 and Observation 3.1. Let be the set of vertices of not contained in any triangle of but contained in some triangle of . If vertices are contained in the triangles of , the size of is for a sufficiently large constant since . Then let be the set of all edges of which form triangles together with the vertices of . Note that for every edge of , its endpoints are contained in the same triangle of by the definition of the core. Thus the size of is at most . Therefore, if more than vertices are contained in the triangles of , .
Lemma 9.
If , there are two subsets and such that is a crown for the triangles of .
Proof. Let be the bipartite graph defined in Definition 1. Notice that is an independent set in since is bipartite. Let be a minimum vertex cover of . Also, let be a matching of size of such that every edge in a matching is incident to exactly one vertex of . Then let , and . Also, let be the set of edges of having their endpoints on . Since , is not fully contained in , and thus .
Then we show that is a crown for the triangles of . First, does not violate the condition in Definition 1 since . Then consider the condition for . Note that is contained in , and does not intersect . Thus for all triangles with , must be in , and then must be in . Therefore, is the set of all edges of that form triangles together with the vertices of , and thus satisfies the secons condition in Definition 1. Also, by construction, is a subgraph of such that all edges of are in , and thus are matched under .
By Lemma 8 and Lemma 9, we can obtain an instance equivalent to such that the union of the triangles has complexity for each instance obtained from Section 3.1 in polynomial time.
Theorem 3.1.
Given a disk graph with its geometric representation, we can find a triangle hitting set of of size in time, if it exists. Without a geometric representation, we can do this in time.
Proof. After the branching and kernelization processes, we have instances one of which is a YES-instance such that the size of the union of the triangles of is at most . For each instance , we remove all vertices not contained in any triangle of . Then the resulting graph has vertices. Then we compute a tree decomposition of of weighted treewidth using Lemma 2.
By applying dynamic programming on as described in Lemma 16, we can find a triangle hitting set of size in time if it exists. Therefore, the total running time is . By letting , we have -time.
If we are not given a geometric representation, we cannot use a tree decomposition of bounded weighted width. Instead, we can compute a tree decomposition of of treewidth . Then the total running time is . By letting , we have -time.
4 Feedback Vertex Set and Odd Cycle Transversal
In this section, we show that the algorithms in [17] for Feedback Vertex Set and Odd Cycle Transversal indeed take time and time algorithms respectively. It is shown in [17] that they take time and time, respectively, but we give a better analysis. We can obtain non-robust algorithms for these problems with better running times, but we omit the description of them.
We obtained better time bounds by classifying the vertices in a kernel used in [17] into two types, and by using the higher order additively weighted Voronoi diagrams. To make the paper self-contained, we present the algorithms in [17] here. The following lemma is a main observation of [17]. Indeed, the statement of the lemma given by [17] is stronger than this, but the following statement is sufficient for our purpose. We say a vertex is deep for a subset of if all neighbors of in are contained in . For a subset of , let be the graph obtained from by contracting each connected component of into a single vertex.
Lemma 10 (Theorem 1.2 of [17]).
Let be a disk graph that has a realization of ply . For a subset of , let be the union of and the set of all deep vertices for . If contains a core of , then the arrangement graph of has treewidth for a set , where is the treewidth of .
In Section 4.1, we show that the size of is . This is our main contribution in this section. Then we solve the two cycle hitting problems using dynamic programming on a tree decomposition of bounded treewidth.
4.0.1 Branching and Cleaning Process.
We first apply the two-step branching process as we did for Triangle Hitting Set. Note that if a vertex set is a feedback vertex set or odd cycle transversal, then is a triangle hitting set.
We apply the cleaning process for each instance we obtained from the earlier branching process. First, we remove all vertices of degree one. Then we keep vertices from each class of false twins as follows. A set of vertices is called a false twin if they are pairwise non-adjacent, and they have the same neighborhood. As observed in [17], for each class of false twins, every minimal feedback vertex set either contains the entire class except for at most one vertex, or none of the vertices in that class. Thus we keep only one vertex (an arbitrary one) from each class of false twins if is an instance of Feedback Vertex Set. Similarly, for each class of false twins, every odd cycle transversal contains the entire class except for at most two vertices, or none of the vertices in that class. Therefore, when we deal with Odd Cycle Transversal, we keep only two vertices from each class of false twins. We remember how many vertices each kept vertex represents and make use of this information in DP.
We have instances one of which is a YES-instance where has a core of size and a geometric representation of ply . Moreover, at most two vertices are in each class of false twins.
4.1 The Number of Deep Vertices
Let be an YES-instance obtained from the branching and cleaning process. Let be a subset of containing a core of . The disks of have ply at most two. In this section, we give an upper bound on the number of deep vertices for . For this, we classify the vertices of in two types: regular and irregular vertices. See Fig. 3(a).
Definition 2.
A vertex of is said to be irregular if belongs to one of the following types:
-
•
contains a vertex of ,
-
•
is contained in a face of ,
-
•
contains a disk represented by a vertex of , or
-
•
intersects three edges of incident to the same face of ,
where denotes the arrangement of the disks represented by . If does not belong to any of the types, we say is regular.

Lemma 11.
The number of deep and regular vertices is .
Proof. If is regular, the neighbors of in form two cliques. To see this, observe that the part of restricted to consists of parallel arcs. That is, the arcs can be sorted in a way that any two consecutive arcs come from the same face of . Also, any two arcs which are not consecutive in the sorted list are not incident to a common face of . In this case, there are two points and on the boundary of such that the line segment intersects all disks represented by intersecting . Then the disks of intersecting contains either or . Therefore, the neighbors of form at most two cliques.
Since each clique has size at most , a deep and regular vertex has at most neighbors in . Let be a deep and regular vertex having exactly neighbors. Consider the additively weighted order- Voronoi diagram of . Here, the distance between a disk of and a point in the plane is defined as , where is the center of and is the radius of . Then the center of is contained in a Voronoi region of the order- Voronoi diagram whose sites are exactly the neighbors of . Since no two deep vertices have the same neighborhood in , each Voronoi region contains at most one deep vertex. Therefore, the number of deep and regular vertices having exactly neighbors is linear in the complexity of for . By Lemma 4, its complexity is . Thus, the number of deep and regular vertices of is .
For irregular vertices, we make use of the fact that the ply of the disks represented by those vertices is at most two. It is not difficult to see that the number of irregular vertices of the first three types is . This is because the complexity of is . To analyze the number of irregular vertices of the fourth type, we use the following lemma.

Lemma 12.
Let be a connected region with edges in the plane, and be a set of interior-disjoint regions contained in such that each region intersects at least three edges of . Then the number of regions of is . Also, the total number of vertices of the regions of is .
Proof. We consider the following planar graph : each vertex of corresponds to an edge of . Also, for each region of with vertices in order, we add edges connecting the vertices corresponding to the edges of containing and for all indices . In addition to this, for each vertex of , we add an edge connecting the vertices of corresponding the two edges of incident to . See Fig. 4. There might be parallel edges in , but a pair of vertices has at most two parallel edges. Let be the set of faces bounded by a pair of parallel edges, and let be the set of the other faces of . Then each region of corresponds to a face of . To prove the first part of the lemma, it suffices to show that the size of is linear in the number of edges of .
By the Euler’s formula, , where and denote the set of vertices and edges of , respectively. Note that no vertex of has degree at most one. Since every face of is incident to at least three edges of , and every edge of is incident to exactly two faces of , . Thus we have . That is, . This implies the first part of the lemma.
Now consider the second part of the lemma. Since a pair of vertices has at most two parallel edges, . Therefore, . Since every edge of the regions of corresponds to an edge of , the total number of edges of the regions of is , and thus the total number of vertices of them is .
Lemma 13.
The number of irregular vertices of is .
Proof. Since the complexity of is , the number of vertices of belonging to the first three cases is . For the fourth case in Definition 2, let be a face of such that has at least three circular arcs which comes from the boundary of . In this case, we replace with . Since , all resulting regions have ply two. See Fig. 3(b). For a face having edges, there are such regions by Lemma 12. Therefore, the number of vertices belonging to the fourth case is linear in the complexity of , which is . Therefore, the lemma holds.
4.2 Feedback Vertex Set
In this section, we show how to compute a feedback vertex set of size in time, if it exists. Without a geometric representation of , we can do this in time. Let be an YES-instance of Feedback Vertex Set from the branching and cleaning processes. Then we have a core of of size . Let be the union of a feedback vertex set of size and the vertex set of all triangles of core of , and let be the union of and all deep vertices for . The size of is by Lemmas 11 and 13. Since contains a feedback vertex set, is a forest.
Therefore, the treewidth of is by Lemma 10. Then we compute a tree decomposition of treewidth , and then apply a dynamic programming on the tree decomposition to solve the problem in time. The total running time is by setting . The weighted treewidth of is by Lemmas 3 and 10. If we have a geometric representation of , we compute a tree decomposition of of weighted treewidth using Lemma 2. By applying dynamic programming on as described in Lemma 17, we can compute a feedback vertex set of size in time. Note that the dynamic programming on Feedback Vertex Set requires global connectivity information but we can track the connectivity by rank-based approach. See also [5]. Since we have instances from branching process, the total running time is by setting .
Theorem 4.1.
Given a disk graph , we can find a feedback vertex set of size in time, if it exists. Given a disk graph together with its geometric representation, we can do this in in time.
4.3 Odd Cycle Transversal
In this section, we present an -time randomised algorithm for Odd Cycle Transversal when a geometric representation of is given. In the case of Odd Cycle Transversal, we do not know how to solve the problem without using a geometric representation of . The algorithm in [17] also uses a geometric representation in this case.
Let be an YES-instance of Odd Cycle Transversal obtained from the branching and cleaning processes. We have a core of of size . Let be the union of the triangles of . Furthermore, let be the union of and the set of all deep vertices for of size by Lemmas 11 and 13. Let . Since does not contain a triangle, it is planar. Thus we use a contraction-decomposition theorem on of [3].
Lemma 14 (Lemma 1.1 of [3]).
Let be a planar graph. Then for any , there exist disjoint sets such that for every and every , treewidth of is . Moreover, these sets can be computed in polynomial time.
In particular, we set and compute sets of vertices of such that for any and any , has treewidth at most . Then there is an index such that for a fixed odd cycle transversal , has size . We iterate over every choice of , and every choice of of size at most . Due to the following lemma, the number of choices of we have to consider is .
Lemma 15 ([17]).
One can compute a candidate set of size for a solution of the Odd Cycle Transversal problem in polynomial time with success probability at least .
For each iteration, we contract and then apply a dynamic programming on a tree decomposition of . The number of iterations is , and has treewidth . By Lemma 3, the weighted treewidth and the treewidth of are and , respectively. We obtain the desired running times by setting (for a robust algorithm) and (for a non-robust algorithm). Then Lemma 18 implies the following theorem.
Theorem 4.2.
Given a disk graph , we can find a odd cycle transversal of size in time w.h.p. Given a disk graph together with its geometric representation, we can do this in in time w.h.p.
5 DP on a Tree Decomposition of Bounded Treewidth
In this section, we give dynamic programming algorithms for Triangle Hitting Set, Feedback Vertex Set, and Odd Cycle Transversalon bounded weighted treewidth graphs. Let be a tree decomposition of with weighted width . Recall that has a clique partition of weight . Let be the set of cliques in the clique partition of , and let be the set of disks represented by the vertices of . Also, we slightly abuse the notation so that denotes the set of disks represented by the vertices of a clique . For the set of disks, we again abuse the notation so that denotes the subgraph of induced by the set of vertices whose corresponding disk is contained in . A tree decomposition is nice if is a rooted tree and each node of belongs to one of the four categories:
-
•
Leaf node. is empty.
-
•
Introduce node. has exactly one child s.t for some
-
•
Forget node. has exactly one child s.t for some
-
•
Join node. has two children s.t .
It is not difficult to see that, given a tree decomposition with weighted width , we can compute a nice tree decomposition with weighted width in polynomial time. It is well known that we can convert a tree decomposition of width into a tree decomposition of width [10]. A slight modification of this conversion preserves the weighted width and the size of a tree asymptotically. Therefore, we may assume that is a nice tree decomposition of weighted width . Moreover, the number of nodes of is . Also, let be the union of for all nodes in the subtree of rooted at .
Lemma 16.
Given a tree decomposition of of weighted treewidth , we can solve Triangle Hitting Set in time.
Proof. For each node and each subset of , we define the subproblem of Triangle Hitting Set as
A solution of Triangle Hitting Set must contain all but two vertices from each clique. Therefore, since
for each node , all but subproblems are set to infinity. We say that a subproblem has a solution if it has a finite value. Moreover, we say is a partial solution of if and is triangle-free. We say that a partial solution is optimal for if its size is minimum among all partial solutions. We give formulas for each of the four cases, and then apply these formulas in a bottom-up manner on . Eventually, we compute , where is the root of .
Leaf node.
If is a leaf, we have only one subproblem .
Introduce node.
Let has one child with . We claim that the following formula holds:
Consider a possible optimal partial solution for the subproblem . Then is exactly , and is also a partial solution of . To see the is optimal for , suppose that there is an optimal partial solution for such that . Then is also a solution of because the disk in dose not intersect to a disk in . Therefore, is an optimal partial solution of , so the formula holds.
Forget node.
Let has one child with . We claim that the following formula holds:
For an optimal partial solution for , we set . Since , is a partial solution of . Conversely, suppose is minimum among all . Let be an optimal partial solution of . Then is also a partial solution of because and .
Join node.
Let has two children and . We claim that
Let is an optimal partial solution of . Then (and ) is a partial solution of (and ). Conversely, if there are optimal partial solutions of and of , we show that is a partial solution of . Suppose there are three disks that pairwise intersect. We may assume that . If , an induced subgraph of is a triangle. Then the partial solution of must hit one of the three disks and . Otherwise, both and intersect with , which implies that . Then, the partial solution of hits the triangle. For both cases, hits the set of three disks. Therefore, is a partial solution of .
Now we compute all finite of node in time. The total time complexity is because the number of nodes in tree is .
Lemma 17.
Given a tree decomposition of of weighted treewidth , we can solve Feedback Vertex Set in time.
Proof. (Sketch.) A dynamic programming for Feedback Vertex Set is more involved because we need to maintain global connectivity information. In particular, for each node , each subset of and each partition of , we define the subproblem of Feedback Vertex Set as
(1) | ||||
(2) | ||||
(3) |
We say that the subproblem has a solution if such a superset of exists. Then we say is a partial solution of . Similar to Lemma 16, we have choices of that has a solution.
Intuitively, if is a partial solution and are contained in the same partition class of , then two vertices of correspond to and are in the same connected component of . However, the number of choices of is because . To reduce the number of subproblems we have to consider, we use a rank-based approach [7]. The rank-based approach ensures that for the fixed , there are partitions among all possible choices of so that they cover all possibilities of the connected components of .
One small issue is that the rank-based approach is designed to get maximum connectivity whereas in Feedback Vertex Set we are aim to minimize the connectivity. As did in [7], we add a special universal vertex to the graph and increase the weighted width by 1. Now the task is to determine if there is a connected subgraph of that contains after deleting a partial solution from . In particular, the task is to compute maximum connectivity among graphs of vertices and edges. Now we use the rank-based approach and we compute partitions for each node and each subset . This completes the proof.
Lemma 18.
Given a tree decomposition of of weighted treewidth , we can solve Odd Cycle Transversal in time.
Proof. (Sketch.) For each node and each function , we define the subproblem of Odd Cycle Transversal as
Similar to Lemma 16, we say that is a partial solution of if satisfies the conditions of the definition above. Also, we say is optimal if the size of is minimum over all partial solutions of . Intuitively, as a solution of a subproblem , we remember not only the set of disks to be deleted as in Lemma 16 but also the bipartition of . The number of subproblems that have a solution is again because for each clique , all but two vertices should be mapped to zero. We focus on the join nodes because the formulas of other categories can be obtained in a straightforward way. Similar to Triangle Hitting Set, we claim that the following formula holds:
The disk in and the disk in do not intersect by definition. Therefore, if is an optimal partial solution of , the domain restriction (and ) of with respect to (and ) is a partial solution of (and ). Conversely, the domain union of two optimal partial solutions (and ) of (and ) makes a partial solution of because for all . This shows that the formula holds. Then Odd Cycle Transversal can be solved in time.
References
- [1] Abu-Khzam, F.N.: A kernelization algorithm for -hitting set. Journal of Computer and System Sciences 76(7), 524–531 (2010)
- [2] An, S., Oh, E.: Feedback Vertex Set on Geometric Intersection Graphs. In: Proceedings of the 32nd International Symposium on Algorithms and Computation (ISAAC 2021). pp. 47:1–47:12 (2021)
- [3] Bandyapadhyay, S., Lochet, W., Lokshtanov, D., Saurabh, S., Xue, J.: Subexponential parameterized algorithms for cut and cycle hitting problems on H-minor-free graphs. In: Proceedings of the 2022 Annual ACM-SIAM Symposium on Discrete Algorithms (SODA). pp. 2063–2084. SIAM (2022)
- [4] Bandyapadhyay, S., Lochet, W., Lokshtanov, D., Saurabh, S., Xue, J.: True contraction decomposition and almost eth-tight bipartization for unit-disk graphs. In: 38th International Symposium on Computational Geometry (SoCG 2022). vol. 224, pp. 11:1–11:16 (2022)
- [5] de Berg, M., Bodlaender, H.L., Kisfaludi-Bak, S., Marx, D., Zanden, T.C.v.d.: A framework for Exponential-Time-Hypothesis–tight algorithms and lower bounds in geometric intersection graphs. SIAM Journal on Computing 49(6), 1291–1331 (2020)
- [6] de Berg, M., Kisfaludi-Bak, S., Monemizadeh, M., Theocharous, L.: Clique-based separators for geometric intersection graphs. In: 32nd International Symposium on Algorithms and Computation (ISAAC 2021). pp. 22:1–22:15 (2021)
- [7] Bodlaender, H.L., Cygan, M., Kratsch, S., Nederlof, J.: Deterministic single exponential time algorithms for connectivity problems parameterized by treewidth. Inf. Comput. 243, 86–111 (2015)
- [8] Bonamy, M., Bonnet, E., Bousquet, N., Charbit, P., Giannopoulos, P., Kim, E.J., Rzążewski, P., Sikora, F., Thomassé, S.: EPTAS and subexponential algorithm for maximum clique on disk and unit ball graphs. Journal of the ACM 68(2) (2021)
- [9] Breu, H., Kirkpatrick, D.G.: Unit disk graph recognition is NP-hard. Computational Geometry 9(1-2), 3–24 (1998)
- [10] Cygan, M., Fomin, F.V., Kowalik, L., Lokshtanov, D., Marx, D., Pilipczuk, M., Pilipczuk, M., Saurabh, S.: Parameterized Algorithms. Springer Publishing Company, Incorporated (2015)
- [11] Demaine, E.D., Fomin, F.V., Hajiaghayi, M., Thilikos, D.M.: Subexponential parameterized algorithms on bounded-genus graphs and H-minor-free graphs. Journal of the ACM (JACM) 52(6), 866–893 (2005)
- [12] Fomin, F.V., Lokshtanov, D., Panolan, F., Saurabh, S., Zehavi, M.: Decomposition of map graphs with applications. In: Proceedings of the 46th International Colloquium on Automata, Languages, and Programming (ICALP 2019). pp. 60:1–60:15 (2019)
- [13] Fomin, F.V., Lokshtanov, D., Panolan, F., Saurabh, S., Zehavi, M.: Finding, hitting and packing cycles in subexponential time on unit disk graphs. Discrete & Computational Geometry 62(4), 879–911 (2019)
- [14] Fomin, F.V., Lokshtanov, D., Saurabh, S.: Bidimensionality and geometric graphs. In: Proceedings of the Twenty-Third Annual ACM-SIAM Symposium on Discrete Algorithms (SODA 2012). pp. 1563–1575 (2012)
- [15] Li, J., Nederlof, J.: Detecting feedback vertex sets of size in time. In: Proceedings of the Fourteenth Annual ACM-SIAM Symposium on Discrete Algorithms (SODA 2022). pp. 971–989 (2020)
- [16] Lokshtanov, D., Narayanaswamy, N., Raman, V., Ramanujan, M., Saurabh, S.: Faster parameterized algorithms using linear programming. ACM Transactions on Algorithms (TALG) 11(2), 1–31 (2014)
- [17] Lokshtanov, D., Panolan, F., Saurabh, S., Xue, J., Zehavi, M.: Subexponential parameterized algorithms on disk graphs (extended abstract). In: Proceedings of the 2022 Annual ACM-SIAM Symposium on Discrete Algorithms (SODA 2022). pp. 2005–2031
- [18] Lokshtanov, D., Saurabh, S., Wahlström, M.: Subexponential parameterized odd cycle transversal on planar graphs. In: IARCS Annual Conference on Foundations of Software Technology and Theoretical Computer Science (FSTTCS 2012). Schloss Dagstuhl-Leibniz-Zentrum fuer Informatik (2012)
- [19] Rosenberger, H.: Order- voronoi diagrams of sites with additive weights in the plane. Algorithmica 6(1), 490–521 (1991)
- [20] Wahlström, M.: Algorithms, measures and upper bounds for satisfiability and related problems. Ph.D. thesis, Department of Computer and Information Science, Linköpings universitet (2007)