Kernelization

From HandWiki

In computer science, a kernelization is a technique for designing efficient algorithms that achieve their efficiency by a preprocessing stage in which inputs to the algorithm are replaced by a smaller input, called a "kernel". The result of solving the problem on the kernel should either be the same as on the original input, or it should be easy to transform the output on the kernel to the desired output for the original problem.

Kernelization is often achieved by applying a set of reduction rules that cut away parts of the instance that are easy to handle. In parameterized complexity theory, it is often possible to prove that a kernel with guaranteed bounds on the size of a kernel (as a function of some parameter associated to the problem) can be found in polynomial time. When this is possible, it results in a fixed-parameter tractable algorithm whose running time is the sum of the (polynomial time) kernelization step and the (non-polynomial but bounded by the parameter) time to solve the kernel. Indeed, every problem that can be solved by a fixed-parameter tractable algorithm can be solved by a kernelization algorithm of this type. This is also true for approximate kernelization.

Example: vertex cover

A standard example for a kernelization algorithm is the kernelization of the vertex cover problem by S. Buss.[1] In this problem, the input is an undirected graph [math]\displaystyle{ G }[/math] together with a number [math]\displaystyle{ k }[/math]. The output is a set of at most [math]\displaystyle{ k }[/math] vertices that includes an endpoint of every edge in the graph, if such a set exists, or a failure exception if no such set exists. This problem is NP-hard. However, the following reduction rules may be used to kernelize it:

  1. If [math]\displaystyle{ k\gt 0 }[/math] and [math]\displaystyle{ v }[/math] is a vertex of degree greater than [math]\displaystyle{ k }[/math], remove [math]\displaystyle{ v }[/math] from the graph and decrease [math]\displaystyle{ k }[/math] by one. Every vertex cover of size [math]\displaystyle{ k }[/math] must contain [math]\displaystyle{ v }[/math] since otherwise too many of its neighbors would have to be picked to cover the incident edges. Thus, an optimal vertex cover for the original graph may be formed from a cover of the reduced problem by adding [math]\displaystyle{ v }[/math] back to the cover.
  2. If [math]\displaystyle{ v }[/math] is an isolated vertex, remove it. An isolated vertex cannot cover any edges, so in this case [math]\displaystyle{ v }[/math] cannot be part of any minimal cover.
  3. If more than [math]\displaystyle{ k^2 }[/math] edges remain in the graph, and neither of the previous two rules can be applied, then the graph cannot contain a vertex cover of size [math]\displaystyle{ k }[/math]. For, after eliminating all vertices of degree greater than [math]\displaystyle{ k }[/math], each remaining vertex can only cover at most [math]\displaystyle{ k }[/math] edges and a set of [math]\displaystyle{ k }[/math] vertices could only cover at most [math]\displaystyle{ k^2 }[/math] edges. In this case, the instance may be replaced by an instance with two vertices, one edge, and [math]\displaystyle{ k=0 }[/math] , which also has no solution.

An algorithm that applies these rules repeatedly until no more reductions can be made necessarily terminates with a kernel that has at most [math]\displaystyle{ k^2 }[/math] edges and (because each edge has at most two endpoints and there are no isolated vertices) at most [math]\displaystyle{ 2k^2 }[/math] vertices. This kernelization may be implemented in linear time. Once the kernel has been constructed, the vertex cover problem may be solved by a brute force search algorithm that tests whether each subset of the kernel is a cover of the kernel. Thus, the vertex cover problem can be solved in time [math]\displaystyle{ O(2^{2k^2}+n+m) }[/math] for a graph with [math]\displaystyle{ n }[/math] vertices and [math]\displaystyle{ m }[/math] edges, allowing it to be solved efficiently when [math]\displaystyle{ k }[/math] is small even if [math]\displaystyle{ n }[/math] and [math]\displaystyle{ m }[/math] are both large.

Although this bound is fixed-parameter tractable, its dependence on the parameter is higher than might be desired. More complex kernelization procedures can improve this bound, by finding smaller kernels, at the expense of greater running time in the kernelization step. In the vertex cover example, kernelization algorithms are known that produce kernels with at most [math]\displaystyle{ 2k }[/math] vertices. One algorithm that achieves this improved bound exploits the half-integrality of the linear program relaxation of vertex cover due to Nemhauser and Trotter.[2] Another kernelization algorithm achieving that bound is based on what is known as the crown reduction rule and uses alternating path arguments.[3] The currently best known kernelization algorithm in terms of the number of vertices is due to (Lampis 2011) and achieves [math]\displaystyle{ 2k-c\log k }[/math] vertices for any fixed constant [math]\displaystyle{ c }[/math].

It is not possible, in this problem, to find a kernel of size [math]\displaystyle{ O(\log k) }[/math], unless P = NP, for such a kernel would lead to a polynomial-time algorithm for the NP-hard vertex cover problem. However, much stronger bounds on the kernel size can be proven in this case: unless coNP [math]\displaystyle{ \subseteq }[/math] NP/poly (believed unlikely by complexity theorists), for every [math]\displaystyle{ \epsilon\gt 0 }[/math] it is impossible in polynomial time to find kernels with [math]\displaystyle{ O(k^{2-\epsilon}) }[/math] edges.[4] It is unknown for vertex cover whether kernels with [math]\displaystyle{ (2-\epsilon)k }[/math] vertices for some [math]\displaystyle{ \epsilon\gt 0 }[/math] would have any unlikely complexity-theoretic consequences.

Definition

In the literature, there is no clear consensus on how kernelization should be formally defined and there are subtle differences in the uses of that expression.

Downey–Fellows notation

In the notation of (Downey Fellows), a parameterized problem is a subset [math]\displaystyle{ L\subseteq\Sigma^*\times\N }[/math] describing a decision problem.

A kernelization for a parameterized problem [math]\displaystyle{ L }[/math] is an algorithm that takes an instance [math]\displaystyle{ (x,k) }[/math] and maps it in time polynomial in [math]\displaystyle{ |x| }[/math] and [math]\displaystyle{ k }[/math] to an instance [math]\displaystyle{ (x',k') }[/math] such that

  • [math]\displaystyle{ (x,k) }[/math] is in [math]\displaystyle{ L }[/math] if and only if [math]\displaystyle{ (x',k') }[/math] is in [math]\displaystyle{ L }[/math],
  • the size of [math]\displaystyle{ x' }[/math] is bounded by a computable function [math]\displaystyle{ f }[/math] in [math]\displaystyle{ k }[/math], and
  • [math]\displaystyle{ k' }[/math] is bounded by a function in [math]\displaystyle{ k }[/math].

The output [math]\displaystyle{ (x',k') }[/math] of kernelization is called a kernel. In this general context, the size of the string [math]\displaystyle{ x' }[/math] just refers to its length. Some authors prefer to use the number of vertices or the number of edges as the size measure in the context of graph problems.

Flum–Grohe notation

In the notation of (Flum Grohe), a parameterized problem consists of a decision problem [math]\displaystyle{ L\subseteq\Sigma^* }[/math] and a function [math]\displaystyle{ \kappa:\Sigma^*\to\N }[/math], the parameterization. The parameter of an instance [math]\displaystyle{ x }[/math] is the number [math]\displaystyle{ \kappa(x) }[/math].

A kernelization for a parameterized problem [math]\displaystyle{ L }[/math] is an algorithm that takes an instance [math]\displaystyle{ x }[/math] with parameter [math]\displaystyle{ k }[/math] and maps it in polynomial time to an instance [math]\displaystyle{ y }[/math] such that

  • [math]\displaystyle{ x }[/math] is in [math]\displaystyle{ L }[/math] if and only if [math]\displaystyle{ y }[/math] is in [math]\displaystyle{ L }[/math] and
  • the size of [math]\displaystyle{ y }[/math] is bounded by a computable function [math]\displaystyle{ f }[/math] in [math]\displaystyle{ k }[/math].

Note that in this notation, the bound on the size of [math]\displaystyle{ y }[/math] implies that the parameter of [math]\displaystyle{ y }[/math] is also bounded by a function in [math]\displaystyle{ k }[/math].

The function [math]\displaystyle{ f }[/math] is often referred to as the size of the kernel. If [math]\displaystyle{ f=k^{O(1)} }[/math], it is said that [math]\displaystyle{ L }[/math] admits a polynomial kernel. Similarly, for [math]\displaystyle{ f={O(k)} }[/math], the problem admits linear kernel.

Kernelizability and fixed-parameter tractability are equivalent

A problem is fixed-parameter tractable if and only if it is kernelizable and decidable.

That a kernelizable and decidable problem is fixed-parameter tractable can be seen from the definition above: First the kernelization algorithm, which runs in time [math]\displaystyle{ O(|x|^c) }[/math] for some c, is invoked to generate a kernel of size [math]\displaystyle{ f(k) }[/math]. The kernel is then solved by the algorithm that proves that the problem is decidable. The total running time of this procedure is [math]\displaystyle{ g(f(k)) + O(|x|^c) }[/math], where [math]\displaystyle{ g(n) }[/math] is the running time for the algorithm used to solve the kernels. Since [math]\displaystyle{ g(f(k)) }[/math] is computable, e.g. by using the assumption that [math]\displaystyle{ f(k) }[/math] is computable and testing all possible inputs of length [math]\displaystyle{ f(k) }[/math], this implies that the problem is fixed-parameter tractable.

The other direction, that a fixed-parameter tractable problem is kernelizable and decidable is a bit more involved. Assume that the question is non-trivial, meaning that there is at least one instance that is in the language, called [math]\displaystyle{ I_{yes} }[/math], and at least one instance that is not in the language, called [math]\displaystyle{ I_{no} }[/math]; otherwise, replacing any instance by the empty string is a valid kernelization. Assume also that the problem is fixed-parameter tractable, i.e., it has an algorithm that runs in at most [math]\displaystyle{ f(k) \cdot |x|^c }[/math] steps on instances [math]\displaystyle{ (x, k) }[/math], for some constant [math]\displaystyle{ c }[/math] and some function [math]\displaystyle{ f(k) }[/math]. To kernelize an input, run this algorithm on the given input for at most [math]\displaystyle{ |x|^{c+1} }[/math] steps. If it terminates with an answer, use that answer to select either [math]\displaystyle{ I_{yes} }[/math] or [math]\displaystyle{ I_{no} }[/math] as the kernel. If, instead, it exceeds the [math]\displaystyle{ |x|^{c+1} }[/math] bound on the number of steps without terminating, then return [math]\displaystyle{ (x,k) }[/math] itself as the kernel. Because [math]\displaystyle{ (x,k) }[/math] is only returned as a kernel for inputs with [math]\displaystyle{ f(k)\cdot |x|^c \gt |x|^{c+1} }[/math], it follows that the size of the kernel produced in this way is at most [math]\displaystyle{ \max\{|I_{yes}|, |I_{no}|, f(k)\} }[/math]. This size bound is computable, by the assumption from fixed-parameter tractability that [math]\displaystyle{ f(k) }[/math] is computable.

More examples

  • Vertex cover parametrized by the size of the vertex cover: The vertex cover problem has kernels with at most [math]\displaystyle{ 2k }[/math] vertices and [math]\displaystyle{ O(k^2) }[/math] edges.[5] Furthermore, for any [math]\displaystyle{ \varepsilon \gt 0 }[/math], vertex cover does not have kernels with [math]\displaystyle{ O(k^{2-\varepsilon}) }[/math] edges unless [math]\displaystyle{ \text{coNP}\subseteq\text{NP/poly} }[/math].[4] The vertex cover problems in [math]\displaystyle{ d }[/math]-uniform hypergraphs has kernels with [math]\displaystyle{ O(k^{d}) }[/math] edges using the sunflower lemma, and it does not have kernels of size [math]\displaystyle{ O(k^{d-\varepsilon}) }[/math] unless [math]\displaystyle{ \text{coNP}\subseteq\text{NP/poly} }[/math].[4]
  • Feedback vertex set parametrized by the size of the feedback vertex set: The feedback vertex set problem has kernels with [math]\displaystyle{ 4k^2 }[/math] vertices and [math]\displaystyle{ O(k^2) }[/math] edges.[6] Furthermore, it does not have kernels with [math]\displaystyle{ O(k^{2-\varepsilon}) }[/math] edges unless [math]\displaystyle{ \text{coNP}\subseteq\text{NP/poly} }[/math].[4]
  • [math]\displaystyle{ k }[/math]-path: The [math]\displaystyle{ k }[/math]-path problem is to decide whether a given graph has a path of length at least [math]\displaystyle{ k }[/math]. This problem has kernels of size exponential in [math]\displaystyle{ k }[/math], and it does not have kernels of size polynomial in [math]\displaystyle{ k }[/math] unless [math]\displaystyle{ \text{coNP}\subseteq\text{NP/poly} }[/math].[7]
  • Bidimensional problems: Many parameterized versions of bidimensional problems have linear kernels on planar graphs, and more generally, on graphs excluding some fixed graph as a minor.[8]

Kernelization for structural parameterizations

While the parameter [math]\displaystyle{ k }[/math] in the examples above is chosen as the size of the desired solution, this is not necessary. It is also possible to choose a structural complexity measure of the input as the parameter value, leading to so-called structural parameterizations. This approach is fruitful for instances whose solution size is large, but for which some other complexity measure is bounded. For example, the feedback vertex number of an undirected graph [math]\displaystyle{ G }[/math] is defined as the minimum cardinality of a set of vertices whose removal makes [math]\displaystyle{ G }[/math] acyclic. The vertex cover problem parameterized by the feedback vertex number of the input graph has a polynomial kernelization:[9] There is a polynomial-time algorithm that, given a graph [math]\displaystyle{ G }[/math] whose feedback vertex number is [math]\displaystyle{ k }[/math], outputs a graph [math]\displaystyle{ G' }[/math] on [math]\displaystyle{ O(k^3) }[/math] vertices such that a minimum vertex cover in [math]\displaystyle{ G' }[/math] can be transformed into a minimum vertex cover for [math]\displaystyle{ G }[/math] in polynomial time. The kernelization algorithm therefore guarantees that instances with a small feedback vertex number [math]\displaystyle{ k }[/math] are reduced to small instances.

See also

Notes

  1. This unpublished observation is acknowledged in a paper of (Buss Goldsmith)
  2. (Flum Grohe)
  3. (Flum Grohe) give a kernel based on the crown reduction that has [math]\displaystyle{ 3k }[/math] vertices. The [math]\displaystyle{ 2k }[/math] vertex bound is a bit more involved and folklore.
  4. 4.0 4.1 4.2 4.3 (Dell van Melkebeek)
  5. (Chen Kanj)
  6. (Thomassé 2010)
  7. (Bodlaender Downey)
  8. (Fomin Lokshtanov)
  9. (Jansen Bodlaender)

References

Further reading