Next Article in Journal
Long Short-Term Memory Neural Network Applied to Train Dynamic Model and Speed Prediction
Next Article in Special Issue
A Compendium of Parameterized Problems at Higher Levels of the Polynomial Hierarchy
Previous Article in Journal
An Application of Manifold Learning in Global Shape Descriptors
Previous Article in Special Issue
Randomized Parameterized Algorithms for the Kidney Exchange Problem
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Practical Access to Dynamic Programming on Tree Decompositions †

1
Institute for Theoretical Computer Science, Universität zu Lübeck, 23562 Lübeck, Germany
2
Department of Computer Science, Kiel University, 24103 Kiel, Germany
*
Author to whom correspondence should be addressed.
This paper is an extended version of our paper published in ESA 2018.
Algorithms 2019, 12(8), 172; https://doi.org/10.3390/a12080172
Submission received: 18 July 2019 / Revised: 5 August 2019 / Accepted: 13 August 2019 / Published: 16 August 2019
(This article belongs to the Special Issue New Frontiers in Parameterized Complexity and Algorithms)

Abstract

:
Parameterized complexity theory has led to a wide range of algorithmic breakthroughs within the last few decades, but the practicability of these methods for real-world problems is still not well understood. We investigate the practicability of one of the fundamental approaches of this field: dynamic programming on tree decompositions. Indisputably, this is a key technique in parameterized algorithms and modern algorithm design. Despite the enormous impact of this approach in theory, it still has very little influence on practical implementations. The reasons for this phenomenon are manifold. One of them is the simple fact that such an implementation requires a long chain of non-trivial tasks (as computing the decomposition, preparing it, …). We provide an easy way to implement such dynamic programs that only requires the definition of the update rules. With this interface, dynamic programs for various problems, such as 3-coloring, can be implemented easily in about 100 lines of structured Java code. The theoretical foundation of the success of dynamic programming on tree decompositions is well understood due to Courcelle’s celebrated theorem, which states that every MSO-definable problem can be efficiently solved if a tree decomposition of small width is given. We seek to provide practical access to this theorem as well, by presenting a lightweight model checker for a small fragment of MSO 1 (that is, we do not consider “edge-set-based” problems). This fragment is powerful enough to describe many natural problems, and our model checker turns out to be very competitive against similar state-of-the-art tools.

1. Introduction

Parameterized algorithms aim to solve intractable problems on instances where the value of some parameter tied to the complexity of the instance is small. This line of research has seen enormous growth in the last few decades and produced a wide range of algorithms—see, for instance, [1]. More formally, a problem is fixed-parameter tractable (in fpt), if every instance I can be solved in time f ( κ ( I ) ) · poly ( | I | ) for a computable function f : N N , where κ ( I ) is the parameter of I. While the impact of parameterized complexity to the theory of algorithms and complexity cannot be overstated, its practical component is much less understood. Very recently, the investigation of the practicability of fixed-parameter tractable algorithms for real-world problems has started to become an important subfield (see e.g., [2,3]). We investigate the practicability of dynamic programming on tree decompositions—one of the most fundamental techniques of parameterized algorithms. A general result explaining the usefulness of tree decompositions was given by Courcelle [4], who showed that every property that can be expressed in monadic second-order logic (MSO) is fixed-parameter tractable if it is parameterized by treewidth of the input graph. By combining this result (known as Courcelle’s Theorem) with the f ( tw ( G ) ) · | G | algorithm of Bodlaender [5] to compute an optimal tree decomposition in fpt-time, a wide range of graph-theoretic problems is known to be solvable on these tree-like graphs. Unfortunately, both ingredients of this approach are very expensive in practice.
One of the major achievements concerning practical parameterized algorithms was the discovery of a practically fast algorithm for treewidth due to Tamaki [6]. Concerning Courcelle’s Theorem, there are currently two contenders concerning efficient implementations of it: D-Flat, an Answer Set Programming (ASP) solver for problems on tree decompositions [7]; and Sequoia, an MSO solver based on model checking games [8]. Both solvers allow to solve very general problems and the corresponding overhead might, thus, be large compared to a straightforward implementation of the dynamic programs for specific problems.

Our Contributions

In order to study the practicability of dynamic programs on tree decompositions, we expand our tree decomposition library Jdrasil with an easy to use interface for such programs: the user only needs to specify the update rules for the different kind of nodes within the tree decomposition. The remaining work—computing a suitable optimized tree decomposition and performing the actual run of the dynamic program—is done by Jdrasil. This allows users to implement a wide range of algorithms within very few lines of code and, thus, gives the opportunity to test the practicability of these algorithms quickly. This interface is presented in Section 3.
In order to balance the generality of MSO solvers and the speed of direct implementations, we introduce an MSO fragment (actually, a MSO 1 fragment), which avoids quantifier alternation, in Section 4. By concentrating on this fragment, we are able to build a model checker, called Jatatosk, that runs nearly as fast as direct implementations of the dynamic programs. To show the feasibility of our approach, we compare the running times of D-Flat, Sequoia, and Jatatosk for various problems. It turns out that Jatatosk is competitive against the other solvers and, furthermore, its behaviour is much more consistent (that is, it does not fluctuate greatly on similar instances). We conclude that concentrating on a small fragment of MSO gives rise to practically fast solvers, which are still able to solve a large class of problems on graphs of bounded treewidth.

2. Preliminaries

All graphs considered in this paper are undirected, that is, they consist of a set of vertices V and of a symmetric edge-relation E V × V . We assume the reader to be familiar with basic graph theoretic terminology—see, for instance, [9]. A tree decomposition of a graph G = ( V , E ) is a tuple ( T , ι ) consisting of a rooted tree T and a mapping ι from nodes of T to sets of vertices of G (which we call bags) such that (1) for all v V there is a non-empty, connected set { x V ( T ) v ι ( x ) } , and (2) for every edge { v , w } E there is a node y in T with { v , w } ι ( y ) . The width of a tree decomposition is the maximum size of one of its bags minus one, and the treewidth of G, denoted by tw ( G ) , is the minimum width any tree decomposition of G must have.
In order to describe dynamic programs over tree decompositions, it turns out be helpful to transform a tree decomposition into a more structured one. A nice tree decomposition is a triple ( T , ι , η ) where ( T , ι ) is a tree decomposition and η : V ( T ) { leaf , introduce , forget , join } is a labeling such that (1) nodes labeled “leaf” are exactly the leaves of T, and the bags of these nodes are empty; (2) nodes n labeled “introduce” or “forget” have exactly one child m such that there is exactly one vertex v V ( G ) with either v ι ( m ) and ι ( n ) = ι ( m ) { v } or v ι ( m ) and ι ( n ) = ι ( m ) \ { v } , respectively; (3) nodes n labeled “join” have exactly two children x , y with ι ( n ) = ι ( x ) = ι ( y ) . A very nice tree decomposition is a nice tree decomposition that also has exactly one node labeled “edge” for every e E ( G ) , which virtually introduces the edge e to the bag—that is, whenever we introduce a vertex, we assume it to be “isolated” in the bag until its incident edges are introduced. It is well known that any tree decomposition can efficiently be transformed into a very nice one without increasing its width (essentially, traverse through the tree and “pull apart” bags) [1]. Whenever we talk about tree decompositions in the rest of the paper, we actually mean very nice tree decompositions. Note that this increases the space needed to store the decomposition but only by a small factor of O ( tw ( G ) ) . Due to the exponential (in tw ( G ) ) space requirement typical for the dynamic programs, this overhead is usually negligible. However, we want to stress that all our interfaces also support “just” nice tree decompositions.
We assume the reader to be familiar with basic logic terminology and give just a brief overview over the syntax and semantic of monadic second-order logic (MSO)—see, for instance, [10] for a detailed introduction. A vocabulary (or signature) τ = ( R 1 a 1 , , R n a n ) is a set of relational symbols R i of arity a i 1 . A τ-structure is a set U—called universe—together with an interpretation R i U R a i of the relational symbols in U. Let x 1 , x 2 , be a sequence of first-order variables and X 1 , X 2 , be a sequence of second-order variables X i of arity ar ( X i ) . The atomic τ -formulas are x i = x j for two first-order variables and R ( x i 1 , , x i k ) , where R is either a relational symbol or a second-order variable of arity k. The set of τ -formulas is inductively defined by (1) the set of atomic τ -formulas; (2) Boolean connections ¬ ϕ , ( ϕ ψ ) , and ( ϕ ψ ) of τ -formulas ϕ and ψ ; (3) quantified formulas x ϕ and x ϕ for a first-order variable x and a τ -formula ϕ ; (4) quantified formulas X ϕ and X ϕ for a second-order variable X of arity 1 and a τ -formula ϕ . The set of free variables of a formula ϕ consists of the variables that appear in ϕ but are not bounded by a quantifier. We denote a formula ϕ with free variables x 1 , , x k , X 1 , , X as ϕ ( x 1 , , x k , X 1 , , X ) . Finally, we say a τ -structure S with an universe U is a model of an τ -formula ϕ ( x 1 , , x k , X 1 , , X ) if there are elements u 1 , , u k U and relations U 1 , , U with U i U ar ( X i ) with ϕ ( u 1 , , u k , U 1 , , U ) being true in S . We write S ϕ ( u 1 , , u k , U 1 , , U ) in this case.
Example 1.
Graphs can be modeled as { E 2 } -structures with universe U = V ( G ) and a symmetric interpretation of E. Properties such as “is 3-colorable” can then be described by formulas as:
ϕ ˜ 3 c o l = R G B ( x R ( x ) G ( x ) B ( x ) ) ( x y E ( x , y ) C { R , G , B } ¬ C ( x ) ¬ C ( y ) ) .
For instance, we have Algorithms 12 00172 i001 ϕ ˜ 3 c o l and Algorithms 12 00172 i002 ϕ ˜ 3 col . Here, and in the rest of the paper, we indicate with the tilde (as in ϕ ˜ ) that we will present a more refined version of ϕ later on.
The model checking problem asks, given a logical structure S and a formula ϕ , whether S ϕ holds. A model checker is a program that solves this problem. To be useful in practice, we also require that a model checker outputs an assignment of the free and existential bounded variables of ϕ in the case of S ϕ .

3. An Interface for Dynamic Programming on Tree Decompositions

It will be convenient to recall a classical viewpoint of dynamic programming on tree decompositions to illustrate why our interface is designed the way it is. We will do so by the guiding example of 3-coloring: Is it possible to color the vertices of a given graph with three colors such that adjacent vertices never share the same color? Intuitively, a dynamic program for 3-coloring will work bottom-up on a very nice tree decomposition and manages a set of possible colorings per node. Whenever a vertex is introduced, the program “guesses” a color for this vertex; if a vertex is forgotten, we have to remove it from the bag and identify configurations that become eventually equal; for join-bags, we just have to take the configurations that are present in both children; and, for edge bags, we have to reject colorings in which both endpoints of the introduced edge have the same color. To formalize this vague algorithmic description, we view it from the perspective of automata theory.

3.1. The Tree Automaton Perspective

Classically, dynamic programs on tree decompositions are described in terms of tree automata [10]. Recall that, in a very nice tree decomposition, the tree T is rooted and binary; we assume that the children of T are ordered. The mapping ι can then be seen as a function that maps the nodes of T to symbols from some alphabet Σ . A naïve approach to manage ι would yield a huge alphabet (depending on the size of the graph). We thus define the so-called tree-index, which is a map idx : V ( G ) { 0 , , tw ( G ) } such that no two vertices that appear in the same bag share a common tree-index. The existence of such an index follows directly from the property that every vertex is forgotten exactly once: We can simply traverse T from the root to the leaves and assign a free index to a vertex V when it is forgotten, and release the used index once we reach an introduce bag for v. The symbols of Σ then only contain the information for which tree-index there is a vertex in the bag. From a theoretician’s perspective, this means that | Σ | depends only on the treewidth; from a programmer’s perspective, the tree-index makes it much easier to manage data structures that are used by the dynamic program.
Definition 1 (Tree Automaton).
A nondeterministic bottom-up tree automaton is a tuple A = ( Q , Σ , Δ , F ) , where Q is the set of states of the automaton, F Q is the set of accepting states, Σ is an alphabet, and Δ ( Q { } ) × ( Q { } ) × Σ × Q is a transition relation in which Q is a special symbol to treat nodes with less than two children. The automaton is deterministic if, for every x , y Q { } and every σ Σ , there is exactly one q Q with ( x , y , σ , q ) Δ .
Definition 2 (Computation of a Tree Automaton).
The computation of a nondeterministic bottom-up tree automaton A = ( Q , Σ , Δ , F ) on a labeled tree ( T , ι ) with ι : V ( T ) Σ and root r V ( T ) is an assignment q : V ( T ) Q such that, for all n V ( T ) , we have (1) ( q ( x ) , q ( y ) , ι ( n ) , q ( n ) ) Δ if n has two children x, y; (2) ( q ( x ) , , ι ( n ) , q ( n ) ) Δ or ( , q ( x ) , ι ( n ) , q ( n ) ) Δ if n has one child x; (3) ( , , ι ( n ) , q ( n ) ) Δ if n is a leaf. The computation is accepting if q ( r ) F .

Simulating Tree Automata

A dynamic program for a decision problem can be formulated as a nondeterministic tree automaton that works on the decomposition—see the left side of Figure 1 for a detailed example. Observe that a nondeterministic tree automaton A will process a labeled tree ( T , ι ) with n nodes in time O ( n ) . When we simulate such an automaton deterministically, one might think that a running time of the form O ( | Q | · n ) is sufficient, as the automaton could be in any potential subset of the Q states at some node of the tree. However, there is a pitfall: for every node, we have to compute the set of potential states of the automaton depending on the sets of potential states of the children of that node, leading to a quadratic dependency on | Q | . This can be avoided for transitions of the form ( , , ι ( x ) , p ) , ( q , , ι ( x ) , p ) , and ( , q , ι ( x ) , p ) , as we can collect potential successors of every state of the child and compute the new set of states in linear time with respect to the cardinality of the set. However, transitions of the form ( q i , q j , ι ( x ) , p ) are difficult, as we now have to merge two sets of states. In detail, let x be a node with children y and z and let Q y and Q z be the set of potential states in which the automaton eventually is in at these nodes. To determine Q x , we have to check for every q i Q y and every q j Q z if there is a p Q such that ( q i , q j , ι ( x ) , p ) Δ . Note that the number of states | Q | can be quite large, as for tree decompositions with bags of size k the set Q is typically of cardinality 2 Ω ( k ) , and we will thus try to avoid this quadratic blow-up.
Observation 1.
A tree automaton can be simulated in time O ( | Q | 2 · n ) .
Unfortunately, the quadratic factor in the simulation cannot be avoided in general, as the automaton may very well contain a transition for all possible pairs of states. However, there are some special cases in which we can circumnavigate the increase in the running time.
Definition 3 (Symmetric Tree Automaton).
A symmetric nondeterministic bottom-up tree automaton is a nondeterministic bottom-up tree automaton A = ( Q , Σ , Δ , F ) in which all transitions ( l , r , σ , q ) Δ satisfy either l = , r = , or l = r .
Assume as before that we wish to compute the set of potential states for a node x with children y and z. Observe that, in a symmetric tree automaton, it is sufficient to consider the set Q y Q z and that the intersection of two sets can be computed in linear time if we take some care in the design of the underlying data structures.
Observation 2.
A symmetric tree automaton can be simulated in time O ( | Q | · n ) .
The right side of Figure 1 illustrates the deterministic simulation of a symmetric tree automaton. The massive time difference in the simulation of tree automata and symmetric tree automata significantly influenced the design of the algorithms in Section 4, in which we try to construct an automaton that is (1) “as symmetric as possible” and (2) allows for taking advantage of the “symmetric parts” even if the automaton is not completely symmetric.

3.2. The Interface

We introduce a simple Java-interface to our library Jdrasil, which originally was developed for the computation of tree decompositions only. The interface is built up from two classes: StateVectorFactory and StateVector. The only job of the factory is to generate StateVector objects for the leaves of the tree decomposition, or with the terms of the previous section: “to define the initial states of the tree automaton.” The StateVector class is meant to model a vector of potential states in which the nondeterministic tree automaton is at a specific node of the tree decomposition. Our interface does not define what a “state” is, or how a collection of states is managed. The only thing the interface requires a user to implement is the behaviour of the tree automaton when it reaches a node of the tree decomposition, i.e., given a StateVector for some node x of the tree decomposition and the information that the next node y reached by the automaton is of a certain type, the user has to compute the StateVector for y. To this end, the interface contains the methods shown in Listing 1.
Listing 1: The four methods of the interface describe the behaviour of the tree automaton. Here, “T” is a generic type for vertices. Each function obtains as a parameter the current bag and a tree-index “idx”. Other parameters correspond to bag-type specifics, e.g., the introduced or forgotten vertex v.
StateVector<T> introduce(Bag<T> b, T v, Map<T, Integer> idx);
StateVector<T> forget(Bag<T> b, T v, Map<T, Integer> idx);
StateVector<T> join(Bag<T> b, StateVector<T> o, Map<T, Integer> idx);
StateVector<T> edge(Bag<T> b, T v, T w, Map<T, Integer> idx);
		
This already rounds up the description of the interface, everything else is done by Jdrasil. In detail, given a graph and an implementation of the interface, Jdrasil will compute a tree decomposition (see [11] for the concrete algorithms used by Jdrasil), transform this decomposition into a very nice tree decomposition, potentially optimize the tree decomposition for the following dynamic program, and finally traverse through the tree decomposition and simulate the tree automaton described by the implementation of the interface. The result of this procedure is the StateVector object assigned to the root of the tree decomposition.

3.3. Example: 3-Coloring

Let us illustrate the usage of the interface with our running example of 3-coloring. A State of the automaton can be modeled as a simple integer array that stores a color (an integer) for every vertex in the bag. A StateVector stores a set of State objects, that is, essentially a set of integer arrays. Introducing a vertex v to a StateVector therefore means that three duplicates of each stored state have to be created, and, for every duplicate, a different color has to be assigned to v. Listing 2 illustrates how this operation could be realized in Java.
Listing 2: Exemplary implementation of the introduce method for 3-coloring. In the listing, the variable states is stored by the StateVector object and represents all currently possible states.
StateVector<T> introduce(Bag<T> b, T v, Map<T, Integer> idx) {
  Set<State> newStates = new HashSet<>();
  for (State state : states) { // ’states’ is the set of states
   for (int color = 1; color <= 3; color++) {
      State newState = new State(state); // copy the state
      newState.colors[idx.get(v)] = color;
      newStates.add(newState);
    }
  }
  states = newStates;
  return this;
}
		
The three other methods can be implemented in a very similar fashion: in the forget-method, we set the color of v to 0; in the edge-method, we remove states in which both endpoints of the edge have the same color; and, in the join-method, we compute the intersection of the state sets of both StateVector objects. Note that, when we forget a vertex v, multiple states may become identical, which is handled here by the implementation of the Java Set-class, which takes care of duplicates automatically.
A reference implementation of this 3-coloring solver is publicly available [12], and a detailed description of it can be found in the manual of Jdrasil [13]. Note that this implementation is only meant to illustrate the interface and that we did not make any effort to optimize it. Nevertheless, this very simple implementation (the part of the program that is responsible for the dynamic program only contains about 120 lines of structured Java-code) performs surprisingly well.

4. A Lightweight Model Checker for an MSO-Fragment

Experiments with the coloring solver of the previous section have shown a huge difference in the performance of general solvers as D-Flat and Sequoia against a concrete implementation of a tree automaton for a specific problem (see Section 5). This is not necessarily surprising, as a general solver needs to keep track of way more information. In fact, an MSO model checker on formula ϕ can probably (unless P = NP ) not run in time f ( | ϕ | + tw ) · poly ( n ) for any elementary function f [14]. On the other hand, it is not clear (in general) what the concrete running time of such a solver is for a concrete formula or problem (see e.g., [15] for a sophisticated analysis of some running times in Sequoia). We seek to close this gap between (slow) general solvers and (fast) concrete algorithms. Our approach is to concentrate only on a fragment of MSO, which is powerful enough to express many natural problems, but which is restricted enough to allow model checking in time that matches or is close to the running time of a concrete algorithm for the problem. As a bonus, we will be able to derive upper bounds on the running time of the model checker directly from the syntax of the input formula.
Based on the interface of Jdrasil, we have implemented a publicly available prototype called Jatatosk [16]. In Section 5, we describe various experiments on different problems on multiple sets of graphs. It turns out that Jatatosk is competitive against the state-of-the-art solvers D-Flat and Sequoia. Arguably, these two programs solve a more general problem and a direct comparison is not entirely fair. However, the experiments do reveal that it seems very promising to focus on smaller fragments of MSO (or perhaps any other description language) in the design of treewidth based solvers.

4.1. Description of the Used MSO-Fragment

We only consider vocabularies τ that contain the binary relation E 2 , and we only consider τ -structures with a symmetric interpretation of E 2 , that is, we only consider structures that contain an undirected graph (but may also contain further relations). The fragment of MSO that we consider is constituted by formulas of the form ϕ = X 1 X k i = 1 n ψ i , where the X j are second-order variables and the ψ i are first-order formulas of the form
ψ i { x y E ( x , y ) χ i , x y E ( x , y ) χ i , x y E ( x , y ) χ i , x y E ( x , y ) χ i , x χ i , x χ i } .
Here, the χ i are quantifier-free first-order formulas in canonical normal form. It is easy to see that this fragment is already powerful enough to encode many classical problems as 3-coloring ( ϕ ˜ 3 col from the introduction is part of the fragment), or vertex-cover (note that, in this form, the formula is not very useful—every graph is a model of the formula if we choose S = V ( G ) ; we will discuss how to handle optimization in Section 4.4):
ϕ ˜ vc = S x y E ( x , y ) S ( x ) S ( y ) .

4.2. A Syntactic Extension of the Fragment

Many interesting properties, such as connectivity, can easily be expressed in MSO, but not directly in the fragment that we study. Nevertheless, a lot of these properties can directly be checked by a model checker if it “knows” what kind of properties it actually checks. We present a syntactic extension of our MSO-fragment which captures such properties. The extension consists of three new second order quantifiers that can be used instead of X i .
The first extension is a k-partition quantifier, which quantifies over partitions of the universe:
partition X 1 , , X k X 1 X 2 X k ( x i = 1 k X i ( x ) ) ( x i = 1 k j i ¬ X i ( x ) ¬ X j ( x ) ) .
This quantifier has two advantages. First, formulas like ϕ ˜ 3 col can be simplified to
ϕ 3 col = partition R , G , B x y E ( x , y ) C { R , G , B } ¬ C ( x ) ¬ C ( y ) ,
and, second, the model checking problem for them can be solved more efficiently: the solver directly “knows” that a vertex must be added to exactly one of the sets.
We further introduce two quantifiers that work with respect to the symmetric relation E 2 (recall that we only consider structures that contain such a relation). The connected X quantifier guesses an X U that is connected with respect to E (in graph theoretic terms), that is, it quantifies over connected subgraphs. The forest F quantifier guesses a F U that is acyclic with respect to E (again in graph theoretic terms), that is, it quantifies over subgraphs that are forests. These quantifiers are quite powerful and allow, for instance, expressing that the graph induced by E 2 contains a triangle as minor:
ϕ triangle - minor = connected R connected G connected B · x ( ¬ R ( x ) ¬ G ( x ) ) ( ¬ G ( x ) ¬ B ( x ) ) ( ¬ B ( x ) ¬ R ( x ) ) x y E ( x , y ) R ( x ) G ( y ) x y E ( x , y ) G ( x ) B ( y ) x y E ( x , y ) B ( x ) R ( y ) .
We can also express problems that usually require more involved formulas in a very natural way. For instance, the feedback-vertex-set problem can be described by the following formula (again, optimization will be handled in Section 4.4):
ϕ ˜ fvs = S forest F x S ( x ) F ( x ) .

4.3. Description of the Model Checker

We describe our model checker in terms of a nondeterministic tree automaton that works on a tree decomposition of the graph induced by E 2 (note that, in contrast to other approaches in the literature (see, for instance, [10]), we do not work on the Gaifman graph). We define any state of the automaton as a bit-vector, and we stipulate that the initial state at every leaf is the zero-vector. For any quantifier or subformula, there will be some area in the bit-vector reserved for that quantifier or subformula and we describe how state transitions effect these bits. The “algorithmic idea” behind the implementation of these transitions is not new, and a reader familiar with folklore dynamic programs on tree decompositions (for instance for vertex-cover or steiner-tree) will probably recognize them. An overview over common techniques can be found in the standard textbooks [1,10].

4.3.1. The Partition Quantifier

We start with a detailed description of the k-partition quantifier partition X 1 , , X q (in fact, we do not implement an additional X quantifier, as we can always state X partition X , X ¯ ): Let k be the maximum bag-size of the tree decomposition. We reserve k · log 2 q bit in the state description (Here, and in the following, log ( x ) denotes the number of bits to store x, i.e., log ( x ) = log 2 ( x ) + 1 ), where each block of length log q indicates in which set X i the corresponding element of the bag is. On an introduce-bag (e.g., for v U ), the nondeterministic automaton guesses an index i { 1 , , q } and sets the log q bits that are associated with the tree-index of v to i. Equivalently, the corresponding bits are cleared when the automaton reaches a forget-bag. As the partition is independent of any edges, an edge-bag does not change any of the bits reserved for the k-partition quantifier. Finally, on join-bags, we may only join states that are identical on the bits describing the partition (as otherwise the vertices of the bag would be in different partitions)—meaning this transition is symmetric with respect to these bits (in terms of Section 3.1).

4.3.2. The Connected Quantifier

The next quantifier we describe is connected X , which has to overcome the difficulty that an introduced vertex may not be connected to the rest of the bag in the moment it got introduced, but may be connected to it when further vertices “arrive.” The solution to this dilemma is to manage a partition of the bag into k k connected components P 1 , , P k , for which we reserve k · log k bit in the state description. Whenever a vertex v is introduced, the automaton either guesses that it is not contained in X and clears the corresponding bits, or it guesses that v X and assigns some P i to v. Since v is isolated in the bag in the moment of its introduction (recall that we work on a very nice tree decomposition), it requires its own component and is therefore assigned to the smallest empty partition P i . When a vertex v is forgotten, there are four possible scenarios:
  • v X , then the corresponding bits are already cleared and nothing happens;
  • v X and v P i with | P i | > 1 , then v is just removed and the corresponding bits are cleared;
  • v X and v P i with | P i | = 1 and there are other vertices w in the bag with w X , then the automaton rejects the configuration, as v is the last vertex of P i and may not be connected to any other partition anymore;
  • v X is the last vertex of the bag that is contained in X, then the connected component is “done”, the corresponding bits are cleared and one additional bit is set to indicate that the connected component cannot be extended anymore.
When an edge { u , v } is introduced, components might need to be merged. Assume u , v X , u P i , and v P j with i < j (otherwise, an edge-bag does not change the state), then we essentially perform a classical union-operation from the well-known union-find data structure. Hence, we assign all vertices that are assigned to P j to P i . Finally, at a join-bag, we may join two states that agree locally on the vertices that are in X (that is, they have assigned the same vertices to some P i ); however, they do not have to agree in the way the different vertices are assigned to P i (in fact, there does not have to be an isomorphism between these assignments). Therefore, the transition at a join-bag has to connect the corresponding components analogous to the edge-bags—in terms of Section 3.1, this transition is not symmetric.
The description of the remaining quantifiers and subformulas is very similar and summarized in the following overview:
Algorithms 12 00172 i003
Algorithms 12 00172 i004

4.4. Extending the Model Checker to Optimization Problems

As the example formulas from the previous section already indicate, performing model checking alone will not suffice to express many natural problems. In fact, every graph is a model of the formula ϕ ˜ vc if S simply contains all vertices. It is therefore a natural extension to consider an optimization version of the model checking problem, which is usually formulated as follows [1,10]: Given a logical structure S , a formula ϕ ( X 1 , , X p ) of the MSO-fragment defined in the previous section with free unary second-order variables X 1 , , X p , and weight functions ω 1 , , ω p with ω i : U Z ; find S 1 , , S p with S i U such that i = 1 p s S i ω i ( s ) is minimized under S ϕ ( S 1 , , S p ) , or conclude that S is not a model for ϕ for any assignment of the free variables. We can now correctly express the (actually weighted) optimization version of vertex-cover as follows:
ϕ vc ( S ) = x y E ( x , y ) S ( x ) S ( y ) .
Similarly, we can describe the optimization version of dominating-set if we assume the input does not have isolated vertices (or is reflexive), and we can also fix the formula ϕ ˜ fvs :
ϕ ds ( S ) = x y E ( x , y ) S ( x ) S ( y ) , ϕ fvs ( S ) = forest F x S ( x ) F ( x ) .
We can also maximize the term i = 1 p s S i ω i ( s ) by multiplying all weights with 1 and, thus, express problems such as independent-set:
ϕ is ( S ) = x y E ( x , y ) ¬ S ( x ) ¬ S ( y ) .
The implementation of such an optimization is straightforward: there is a 2-partition quantifier for every free variable X i that partitions the universe into X i and X ¯ i . We assign a current value of i = 1 p s S i ω i ( s ) to every state of the automaton, which is adapted if elements are “added” to some of the free variables at introduce nodes. Note that, since we optimize an affine function, this does not increase the state space: even if multiple computational paths lead to the same state with different values at some node of the tree, it is well defined which of these values is the optimal one. Therefore, the cost of optimization only lies in the partition quantifier, that is, we pay with k bits in the state description of the automaton per free variable—independently of the weights (Of course, for each state, we still have to store the weight and we have to perform arithmetic operations on these weights).

4.5. Handling Symmetric and Non-Symmetric Joins

In Section 4.3, we have defined the states of our automaton with respect to a formula, Table 1 gives an overview of the number of bits we require for the different parts of the formula. Let bit ( ϕ , k ) be the number of bits that we have to reserve for a formula ϕ and a tree decomposition of maximum bag size k that is, the sum over the required bits of each part of the formula. By Observation 1, this implies that we can simulate the automaton (and hence solve the model checking problem) in time O * ( ( 2 bit ( ϕ , k ) ) 2 ) , or by Observation 2 in time O * ( 2 bit ( ϕ , k ) ) , if the automaton is symmetric (The notation O * supresses polynomial factors). Unfortunately, this is not always the case, in fact, only the quantifier partition X 1 , , X q , the bits needed to optimize over free variables, as well as the formulas that do not require any bits, yield a symmetric tree automaton. This means that the simulation is wasteful if we consider a mixed formula (for instance, one that contains a partition and a connected quantifier). To overcome this issue, we partition the bits of the state description into two parts: first, the “symmetric” bits of the quantifiers partition X 1 , , X q and the bits required for optimization, and in the “asymmetric” ones of all other elements of the formula. Let symmetric ( ϕ , k ) and asymmetric ( ϕ , k ) be defined analogously to bit ( ϕ , k ) . We implement the join of states as in the following lemma, allowing us to deduce the running time of the model checker for concrete formulas. Table 2 provides an overview for formulas presented here.
Lemma 1.
Let x be a join node of T with children y and z, and let Q y and Q z be sets of states in which the automaton may be at y and z. Then, the set Q x of states in which the automaton may be at node x can be computed in time O * ( 2 symmetric ( ϕ , k ) + 2 · asymmetric ( ϕ , k ) ) .
Proof. 
To compute Q x , we first split Q y into B 1 , , B q such that all elements in one B i share the same “symmetric bits.” This can be done in time | Q y | using bucket-sort. Note that we have q 2 symmetric ( ϕ , k ) and | B i | 2 asymmetric ( ϕ , k ) . With the same technique, we identify for every element v in Q z its corresponding partition B i . Finally, we compare v with the elements in B i to identify those for which there is a transition in the automaton. This yields an overall running time of the form | Q z | · max i = 1 q | B i | 2 bit ( ϕ , k ) · 2 asymmetric ( ϕ , k ) = 2 symmetric ( ϕ , k ) + 2 · asymmetric ( ϕ , k ) . □

5. Applications and Experiments

To show the feasibility of our approach, we have performed experiments for widely investigated graph problems: 3-coloring (Figure 2), vertex-cover (Figure 3), dominating-set (Figure 4), independent-set (Figure 5), and feedback-vertex-set (Figure 6). All experiments were performed on an Intel Core processor containing four cores of 3.2 GHz each and 8 Gigabyte RAM. Jdrasil was used with Java 8.1 and both Sequoia and D-Flat were compiled with gcc 7.2. All compilations were performed with the default optimization settings. The implementation of Jatatosk uses hashing to realize Lemma 1, which works well in practice. We use a data set assembled from different sources containing graphs with 18 to 956 vertices and treewidth 3 to 13. The first source is a collection of transit graphs from GTFS-transit feeds [17] that was also used for experiments in [18], the second source is real-world instances collected in [19], and the last one is that of the PACE challenge [2] with treewidth at most 11. In each of the experiments, the left picture always shows the difference of Jatatosk against D-Flat and Sequoia. A positive bar means that Jatatosk is faster by this amount in seconds, and a negative bar means that either D-Flat or Sequoia is faster by that amount. The bars are capped at 100 seconds. On every instance, Jatatosk was compared against the solver that was faster on this particular instance. The image also shows for every instance the treewidth of the input. The right image always shows a cactus plot that visualizes the number of instances that can be solved by each of the solvers in x seconds, that is, faster growing functions are better. For each experiment, there is a table showing the average, standard deviation, and median of the time (in seconds) each solver needed to solve the problem. The best values are highlighted.
The experiments reveal that Jatatosk is faster than its competitors on many instances. However, there are also formulas such as the one for the vertex cover problem on which one of the other solvers performs better on some of the instances. For an overall picture, the cactus plot in Figure 7 is the sum of all cactus plots from the experiments. It reveals that overall Jatatosk in fact outperforms its competitors. However, we stress once more that the comparison is not completely fair, as both Sequoia and D-Flat are powerful enough to model check the whole of MSO (and actually also MSO 2 ), while Jatatosk can only handle a fragment of MSO 1 .
As we can see, Jatatosk outperforms Sequoia and D-Flat if we consider the graph coloring problem. Its average time is more than a factor of 10 smaller than the average time of its competitors and Jatatosk solves about 30% of the instances more. This result is not surprising, as the fragment used by Jatatosk is directly tailored towards coloring and, thus, Jatatosk has a natural advantage.
The vertex-cover problem is solved best by Sequoia, which becomes apparent if we consider the difference plot. Furthermore, the average time used by Sequoia is better than the time used by Jatatosk. However, considering the cactus plot, the difference between Jatatosk and Sequoia with respect to solved instances is small, while D-Flat falls behind a bit. We assume that the similarity between Jatatosk and Sequoia is because both compile internally a similar algorithm.
In this experiment Jatatosk performs best with respect to all: the difference plot, the cactus plot, and the used average time. However, the difference to Sequoia is small and there are in fact a couple of instances that are solved faster by Sequoia. We are surprised by this result, as the worst-case running time of Jatatosk for dominating-set is O * ( 8 k ) and, thus, far from optimal. Furthermore, Sequoia promises in theory a better performance with a running time of the form O * ( 5 k ) [15].
This is the simplest test set that we consider within this paper, which is reflected by the fact that all solvers are able to solve almost the whole set. The difference between Jatatosk and Sequoia is minor: while Jatatosk has a slightly better peek-performance, there are more instances that are solved faster by Sequoia than the other way around.
This is the hardest test set that we consider within this paper, which is reflected by the fact that no solver is able to solve much more than 50% of the instances. Jatatosk outperforms its competitors, which is reflected in the difference plot, the cactus plot, and the used average time. We assume this is because Jatatosk uses the dedicated forest quantifier directly, while the other tools have to infer the algorithmic strategy from a more general formula.

6. Conclusions and Outlook

We investigated the practicability of dynamic programming on tree decompositions, which is arguably one of the cornerstones of parameterized complexity theory. We implemented a simple interface for such programs and used it to build a competitive graph coloring solver with just a few lines of code. We hope that this interface allows others to implement and explore various dynamic programs. The whole power of these algorithms is well captured by Courcelle’s Theorem, which states that there is an efficient version of such a program for every problem definable in monadic second-order logic. We took a step towards practice by implementing a “lightweight” version of a model checker for a small fragment of the logic. This fragment turns out to be powerful enough to express many natural problems such as 3-coloring, feedback-vertex-set, and more. Various experiments showed that the model checker is competitive against the state-of-the-art solvers D-Flat and Sequoia. It therefore seems promising, from a practical perspective, to study smaller fragments of MSO.

Author Contributions

Conceptualization, M.B. and S.B.; methodology, M.B. and S.B.; software, M.B.; validation, S.B.; formal analysis, M.B. and S.B.; investigation, S.B.; resources, S.B.; data curation, M.B. and S.B.; writing—original draft preparation, M.B. and S.B.; writing—review and editing, M.B. and S.B.; visualization, M.B.; supervision, M.B. and S.B.; project administration, M.B. and S.B.; funding acquisition, M.B. and S.B.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Cygan, M.; Fomin, F.V.; Kowalik, Ł.; Lokshtanov, D.; Marx, D.; Pilipczuk, M.; Pilipczuk, M.; Saurabh, S. Parameterized Algorithms; Springer: Berlin, Germany, 2015. [Google Scholar] [CrossRef]
  2. Dell, H.; Husfeldt, T.; Jansen, B.M.P.; Kaski, P.; Komusiewicz, C.; Rosamond, F.A. The first parameterized algorithms and computational experiments challenge. In Proceedings of the 11th International Symposium on Parameterized and Exact Computation (IPEC), Aarhus, Denmark, 24–26 August 2016; pp. 30:1–30:9. [Google Scholar] [CrossRef]
  3. Dell, H.; Komusiewicz, C.; Talmon, N.; Weller, M. The PACE 2017 parameterized algorithms and computational experiments challenge: The second iteration. In Proceedings of the 12th International Symposium on Parameterized and Exact Computation (IPEC), Vienna, Austria, 6–8 September 2017; pp. 30:1–30:12. [Google Scholar] [CrossRef]
  4. Courcelle, B. The monadic second-order logic of graphs. I. Recognizable sets of finite graphs. Inf. Comput. 1990, 85, 12–75. [Google Scholar] [CrossRef] [Green Version]
  5. Bodlaender, H.L. A linear-time algorithm for finding tree-decompositions of small treewidth. SIAM J. Comput. 1996, 25, 1305–1317. [Google Scholar] [CrossRef]
  6. Tamaki, H. Positive-Instance driven dynamic programming for treewidth. In Proceedings of the 25th Annual European Symposium on Algorithms (ESA), Vienna, Austria, 4–6 September 2017; pp. 68:1–68:13. [Google Scholar] [CrossRef]
  7. Abseher, M.; Bliem, B.; Charwat, G.; Dusberger, F.; Hecher, M.; Woltran, S. D-FLAT: Progress Report. Available online: https://www.dbai.tuwien.ac.at/research/report/dbai-tr-2014-86.pdf (accessed on 5 June 2019).
  8. Langer, A.J. Fast Algorithms for Decomposable Graphs. Ph.D. Thesis, The RWTH Aachen University, Aachen, Germany, 2013. [Google Scholar]
  9. Diestel, R. Graph Theory: Springer Graduate Text Gtm 173, 4th ed.; Springer: Berlin, Germany, 2012; Volume 173. [Google Scholar]
  10. Flum, J.; Grohe, M. Parameterized Complexity Theory; Springer: Berlin, Germany, 2006. [Google Scholar] [CrossRef]
  11. Bannach, M.; Berndt, S.; Ehlers, T. Jdrasil: A modular library for computing tree decompositions. In Proceedings of the 16th International Symposium on Experimental Algorithms (SEA), London, UK, 21–23 June 2017; pp. 28:1–28:21. [Google Scholar] [CrossRef]
  12. Bannach, M. Jdrasil for Graph Coloring. Available online: https://github.com/maxbannach/Jdrasil-for-GraphColoring (accessed on 23 January 2019). Commit: a5e52a8.
  13. Bannach, M.; Berndt, S.; Ehlers, T. Jdrasil. Available online: http://www.github.com/maxbannach/jdrasil (accessed on 5 June 2019). Commit: dfa1eee.
  14. Frick, M.; Grohe, M. The complexity of first-order and monadic second-order logic revisited. Ann. Pure Appl. Logic 2004, 130, 3–31. [Google Scholar] [CrossRef] [Green Version]
  15. Kneis, J.; Langer, A.; Rossmanith, P. Courcelle’s theorem—A game-theoretic approach. Discret. Optim. 2011, 8, 568–594. [Google Scholar] [CrossRef]
  16. Bannach, M.; Berndt, S. Jatatosk. Available online: https://github.com/maxbannach/Jatatosk/commit/45e306cfac5a273416870ec0bd9cd2c7f39a6932 (accessed on 8 April 2019).
  17. Fichte, J.K. gtfs2graphs—A Transit Feed to Graph Format Converter. Available online: https://github.com/daajoe/gtfs2graphs/commit/219944893f874b365de1ed87fc265fd5d19d5972 (accessed on 20 April 2018).
  18. Fichte, J.K.; Lodha, N.; Szeider, S. SAT-Based local improvement for finding tree decompositions of small width. In Proceedings of the International Conference on Theory and Applications of Satisfiability Testing (SAT), Melbourne, Australia, 28 August–1 September 2017; Springer: Cham, Switzerland, 2017; pp. 401–411. [Google Scholar]
  19. Abseher, M.; Dusberger, F.; Musliu, N.; Woltran, S. Improving the efficiency of dynamic programming on tree decompositions via machine learning. J. Artif. Intell. Res. 2015, 58, 275–282. [Google Scholar] [CrossRef]
Figure 1. The left picture shows a part of a tree decomposition of the grid graph Algorithms 12 00172 i005 with vertices { 0 , , 8 } (i.e., with rows { i , i + 1 , i + 2 } for i { 0 , 3 , 6 } and columns { i , i + 3 , i + 6 } for i { 0 , 1 , 2 } ). The index of a bag shows the type of the bag: a positive sign means “introduce”, a negative one “forget”, a pair represents an “edge”-bag, and the text is self explanatory. Solid lines represent real edges of the decomposition, while dashed lines illustrate a path (that is, some bags are skipped). On the left branch of the decomposition, a run of a nondeterministic tree automaton with tree-index 0 1 2 3 4 5 6 7 8 2 3 0 1 2 3 0 1 0 for 3-coloring is illustrated. To increase readability, states of the automaton are connected to the corresponding bags with gray lines, and, for some nodes, the states are omitted. In the right picture, the same automaton is simulated deterministically.
Figure 1. The left picture shows a part of a tree decomposition of the grid graph Algorithms 12 00172 i005 with vertices { 0 , , 8 } (i.e., with rows { i , i + 1 , i + 2 } for i { 0 , 3 , 6 } and columns { i , i + 3 , i + 6 } for i { 0 , 1 , 2 } ). The index of a bag shows the type of the bag: a positive sign means “introduce”, a negative one “forget”, a pair represents an “edge”-bag, and the text is self explanatory. Solid lines represent real edges of the decomposition, while dashed lines illustrate a path (that is, some bags are skipped). On the left branch of the decomposition, a run of a nondeterministic tree automaton with tree-index 0 1 2 3 4 5 6 7 8 2 3 0 1 2 3 0 1 0 for 3-coloring is illustrated. To increase readability, states of the automaton are connected to the corresponding bags with gray lines, and, for some nodes, the states are omitted. In the right picture, the same automaton is simulated deterministically.
Algorithms 12 00172 g001
Figure 2. Experiment for the 3-coloring problem.
Figure 2. Experiment for the 3-coloring problem.
Algorithms 12 00172 g002
Figure 3. Experiment for the vertex-cover problem.
Figure 3. Experiment for the vertex-cover problem.
Algorithms 12 00172 g003
Figure 4. Experiment for the dominating-set problem.
Figure 4. Experiment for the dominating-set problem.
Algorithms 12 00172 g004
Figure 5. Experiment for the independent-set problem.
Figure 5. Experiment for the independent-set problem.
Algorithms 12 00172 g005
Figure 6. Experiment for the feedback-vertex-set problem.
Figure 6. Experiment for the feedback-vertex-set problem.
Algorithms 12 00172 g006
Figure 7. The cactus plot is the sum of the plots from Figure 2, Figure 3, Figure 4, Figure 5 and Figure 6. Over all experiments, there are 300 instances in total that were tested. Jatatosk solves 88% of these, while Sequoia only manages to solve 76% and D-Flat 50% of the instances.
Figure 7. The cactus plot is the sum of the plots from Figure 2, Figure 3, Figure 4, Figure 5 and Figure 6. Over all experiments, there are 300 instances in total that were tested. Jatatosk solves 88% of these, while Sequoia only manages to solve 76% and D-Flat 50% of the instances.
Algorithms 12 00172 g007
Table 1. The table shows the precise number of bits we reserve in the description of a state of the tree automaton for different quantifiers and formulas. The values are with respect to a tree decomposition with maximum bag size k.
Table 1. The table shows the precise number of bits we reserve in the description of a state of the tree automaton for different quantifiers and formulas. The values are with respect to a tree decomposition with maximum bag size k.
Quantifier/FormulaNumber of Bits
free variables X 1 , , X q q · k
partition X 1 , , X q k · log q
connected X k · log k + 1
forest X k · log k
x y E ( x , y ) χ i 0
x y E ( x , y ) χ i k
x y E ( x , y ) χ i k + 1
x y E ( x , y ) χ i 1
x χ i 0
x χ i 1
Table 2. The table gives an overview of formulas ϕ used within this paper, together with the values symmetric ( ϕ , k ) and asymmetric ( ϕ , k ) , as well as the precise time our algorithm will require to model check an instance for that particular formula.
Table 2. The table gives an overview of formulas ϕ used within this paper, together with the values symmetric ( ϕ , k ) and asymmetric ( ϕ , k ) , as well as the precise time our algorithm will require to model check an instance for that particular formula.
ϕ symmetric ( ϕ , k ) asymmetric ( ϕ , k ) Time
ϕ 3 col k · log ( 3 ) 0 O * ( 3 k )
ϕ vc ( S ) k0 O * ( 2 k )
ϕ ds ( S ) kk O * ( 8 k )
ϕ triangle - minor 0 3 k · log ( k ) + 3 O * ( k 6 k )
ϕ fvs ( S ) k k · log ( k ) O * ( 2 k k 2 k )

Share and Cite

MDPI and ACS Style

Bannach, M.; Berndt, S. Practical Access to Dynamic Programming on Tree Decompositions. Algorithms 2019, 12, 172. https://doi.org/10.3390/a12080172

AMA Style

Bannach M, Berndt S. Practical Access to Dynamic Programming on Tree Decompositions. Algorithms. 2019; 12(8):172. https://doi.org/10.3390/a12080172

Chicago/Turabian Style

Bannach, Max, and Sebastian Berndt. 2019. "Practical Access to Dynamic Programming on Tree Decompositions" Algorithms 12, no. 8: 172. https://doi.org/10.3390/a12080172

APA Style

Bannach, M., & Berndt, S. (2019). Practical Access to Dynamic Programming on Tree Decompositions. Algorithms, 12(8), 172. https://doi.org/10.3390/a12080172

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop