Parallel Tree Contraction
In computer science, parallel tree contraction is a broadly applicable technique for the parallel solution of a large number of tree problems, and is used as an algorithm design technique for the design of a large number of parallel graph algorithms. Parallel tree contraction was introduced by Gary L. Miller and John H. Reif,[1] and has subsequently been modified to improve efficiency by X. He and Y. Yesha,[2] Hillel Gazit, Gary L. Miller and Shang-Hua Teng[3] and many others.[4]
Tree contraction has been used in designing many efficient parallel algorithms, including expression evaluation, finding lowest common ancestors, tree isomorphism, graph isomorphism, maximal subtree isomorphism, common subexpression elimination, computing the 3-connected components of a graph, and finding an explicit planar embedding of a planar graph[5]
Based on the research and work on parallel tree contraction, various algorithms have been proposed targeting to improve the efficiency or simplicity of this topic. This article hereby focuses on a particular solution, which is a variant of the algorithm by Miller and Reif, and its application.
Introduction
Over the past several decades there has been significant research on deriving new parallel algorithms for a variety of problems, with the goal of designing highly parallel (polylogarithmic depth), work-efficient (linear in the sequential running time) algorithms.[1] For some problems, tree turns out to be a nice solution. Addressing these problems, we can sometimes get more parallelism simply by representing our problem as a tree.
Considering a generic definition of a tree, there is a root vertex, and several child vertices attached to the root.[6] And the child vertices might have children themselves, and so on so forth. Eventually, the paths come down to leaves, which are defined to be the terminal of a tree. Then based on this generic tree, we can further come up with some special cases: (1) balanced binary tree; (2) linked list.[7] A balanced binary tree has exactly two branches for each vertex except for leaves. This gives a O(log n) bound on the depth of the tree.[8] A linked list is also a tree where every vertex has only one child. We can also achieve O(log n) depth using symmetry breaking.[9]
Given the general case of a tree, we would like to keep the bound at O(log n) no matter it is unbalanced or list-like or a mix of both. To address this problem, we make use of an algorithm called prefix sum by using the Euler tour technique.[10] With the Euler tour technique, a tree could be represented in a flat style, and thus prefix sum could be applied to an arbitrary tree in this format. In fact, prefix sum can be used on any set of values and binary operation which form a group: the binary operation must be associative, every value must have an inverse, and there exists an identity value.
With a bit of thought, we can find some exceptional cases where prefix sum becomes incapable or inefficient. Consider the example of multiplication when the set of values includes 0. Or there are some commonly desired operations are max() and min() which do not have inverses. The goal is to seek an algorithm which works on all trees, in expected O(n) work and O(log n) depth. In the following sections, a Rake/Compress algorithm will be proposed to fulfill this goal.[11]
Definitions
Before going into the algorithm itself, we first look at a few terminologies that will be used later.
- Rake[12] – Rake step joins every left leaf of binary nodes to the parent. By join, we mean that it undergoes a functional process which achieves the operation we want to make. An example of rake is given in Figure 1.
- Compress[12] – Compress step is actually a sequence of several events: (1) Find an independent set of unary nodes. (Independence here is defined such that no two are neighbors, meaning no parent to child relation) (2) Join each node in independent set with its child (Note that independent set is not unique). An example of compress is given in Figure 2.
And in order to solve actual problems using tree contraction, the algorithm has a structure:
Repeat until tree becomes a unary node { Rake; Compress; }
Analysis
For the moment, let us assume that all nodes have less than three children, namely binary. Generally speaking, as long as the degree is bounded, the bounds will hold.[13] But we will analyze the binary case for simplicity. In the two “degenerate” cases listed above, the rake is the best tool for dealing with balanced binary trees, and compress is the best for linked lists. However, arbitrary trees will have to require a combination of these operations. By this combination, we claim a theorem that
- Theorem: After O(log n) expected rake and compress steps, a tree is reduced to a single node.
Now rephrase the tree contraction algorithm as follows:
- Input: A binary tree rooted at r
- Output: A single node
- Operation: A sequence of contraction steps, each consisting of a rake operation and a compress operation (in any order). The rake operation removes all the leaf nodes in parallel. The compress operation finds an independent set of unary nodes and splice out the selected nodes.
To approach the theorem, we first take a look at a property of a binary tree. Given a binary tree T, we can partition the nodes of T into 3 groups: [math]\displaystyle{ T_0 }[/math] contains all leaf nodes, [math]\displaystyle{ T_1 }[/math] contains all nodes with 1 child, and [math]\displaystyle{ T_2 }[/math] contains all nodes with 2 children. It is easy to see that: [math]\displaystyle{ V(T) = T_0 \cup T_1 \cup T_2 }[/math]. Now we propose:
- Claim: [math]\displaystyle{ |T_0| = |T_2| + 1 }[/math]
This claim can be proved by strong induction on the number of nodes. It is easy to see that the base case of n=1 trivially holds. And we further assume the claim also holds for any tree with at most n nodes. Then given a tree with n+1 nodes rooted at r, there appears to be two cases:
- If r has only one subtree, consider the subtree of r. We know that the subtree has the same number of binary nodes and the same number of leaf nodes as the whole tree itself. This is true since the root is a unary node. And based the previous assumption, a unary node does not change either [math]\displaystyle{ T_0 }[/math] or [math]\displaystyle{ T_2 }[/math].
- If r has two subtrees, we define [math]\displaystyle{ T_0^L, T_2^L }[/math] to be the leaf nodes and binary nodes in the left subtree, respectively. Similarly, we define the same [math]\displaystyle{ T_0^R, T_2^R }[/math] for the right subtree. From previous, there is [math]\displaystyle{ |T_0^L| = |T_2^L| + 1 }[/math] and [math]\displaystyle{ |T_0^R| = |T_2^R| + 1 }[/math]. Also we know that T has [math]\displaystyle{ |T_0^L| + |T_0^R| }[/math] leaf nodes and [math]\displaystyle{ |T_2^L| + |T_2^R| + 1 }[/math] binary nodes. Thus, we can derive:
- [math]\displaystyle{ |T_0^L| + |T_0^R| = |T_2^L| + 1 + |T_2^R| + 1 = (|T_2^L| + |T_2^R| + 1) + 1 }[/math]
which proves the claim.
Following the claim, we then prove a lemma, which leads us to the theorem.
- Lemma: The number of nodes of after a contraction step is reduced by a constant factor in expectation.
Assume the number of nodes before the contraction to be m, and m' after the contraction. By definition, the rake operation deletes all [math]\displaystyle{ T_0 }[/math] and the compress operation deletes at least 1/4 of [math]\displaystyle{ T_1 }[/math] in expectation. All [math]\displaystyle{ T_2 }[/math] remains. Therefore, we can see:
- [math]\displaystyle{ E[m'] \leq |T_2| + \tfrac{3}{4}*|T_1| \leq \tfrac{3}{4} + \tfrac{3}{4}*|T_1| + \tfrac{3}{2}*|T_2| = \tfrac{3}{4}(1 + |T_1| + 2*|T_2|) = \tfrac{3}{4}(|T_0| + |T_1| + |T_2|) = \tfrac{3}{4}m }[/math]
Finally, based on this lemma, we can conclude that if the nodes are reduced by a constant factor in each iteration, after [math]\displaystyle{ O(\log n) }[/math], there will be only one node left.[14]
Applications
Expression Evaluation
To evaluate an expression given as a binary tree (this problem also known as binary expression tree),[15] consider that: An arithmetic expression is a tree where the leaves have values from some domain and each internal vertex has two children and a label from {+, x, %}. And further assume that these binary operations can be performed in constant time.
We now show the evaluation can be done with parallel tree contraction.[16]
- Step 1. Assign expressions to every node. The expression of a leaf is simply the value that it contains. Write L + R, L − R, or L × R for the operators, where L and R are the values of the expressions in the left and right subtrees, respectively.
- Step 2. When a left (right) child with 0 children is merged into an operator, replace L (R) with the value of the child.
- Step 3. When a node has 1 child, it has an expression that is a function of one variable. When a left (right) child with 1 child is merged into an operator, replace L (R) with the expression and change the variable in the expression to L (R) if appropriate.
In a node with 2 children, the operands in the expression are f(L) and g(R), where f and g are linear functions, and in a node with 1 child, the expression is h(x), where h is a linear function and x is either L or R. We prove this invariant by induction. At the beginning, the invariant is clearly satisfied. There are three types of merges that result in a not fully evaluated expression. (1) A 1-child node is merged into a 2-children node. (2) A leaf is merged into a 2-children node. (3) A 1-child node is merged into a 1-child node. All three types of merges do not change the invariant. Therefore, every merge simply evaluates or composes linear functions, which takes constant time [17]
References
- ↑ 1.0 1.1 Gary L. Miller and John H. Reif, Parallel Tree Contraction--Part I: Fundamentals., 1989
- ↑ X. He and Y. Yesha, "Binary tree algebraic computation and parallel algorithms for simple graphs.", Journal of Algorithms, 1988, pp 92-113
- ↑ Hillel Gazit, Gary L. Miller and Shang-Hua Teng, Optimal tree contraction in the EREW model, Springer, 1988
- ↑ Karl Abrahamson and et al., "A simple parallel tree contraction algorithm.", Journal of Algorithms, 1989, pp 287-302
- ↑ John H. Reif and Stephen R. Tate, Dynamic parallel tree contraction, Proceedings of the sixth annual ACM symposium on Parallel algorithms and architectures (ACM), 1994
- ↑ Thomas H. Cormen, Charles E. Leiserson, Ronald L. Rivest, and Clifford Stein. Introduction to Algorithms, Second Edition. MIT Press and McGraw-Hill, 2001. ISBN:0-262-03293-7 . Section 10.4: Representing rooted trees, pp. 214–217. Chapters 12–14 (Binary Search Trees, Red-Black Trees, Augmenting Data Structures), pp. 253–320.
- ↑ Donald Knuth. The Art of Computer Programming: Fundamental Algorithms, Third Edition. Addison-Wesley, 1997. ISBN:0-201-89683-4 . Section 2.3: Trees, pp. 308–423.
- ↑ "Chapter 6. Bounded height trees and tree-depth", Sparsity: Graphs, Structures, and Algorithms, Algorithms and Combinatorics, 28, Heidelberg: Springer, 2012, pp. 115–144, doi:10.1007/978-3-642-27875-4, ISBN 978-3-642-27874-7.
- ↑ Andrew Goldberg, Serge Plotkin, and Gregory Shannon, Parallel symmetry-breaking in sparse graphs, Proceedings of the nineteenth annual ACM symposium on Theory of computing (ACM), 1987
- ↑ Euler tour trees - in Lecture Notes in Advanced Data Structures. Prof. Erik Demaine; Scribe: Katherine Lai.
- ↑ Gary L. Miller and John H. Reif, Parallel tree contraction and its application, Defense Technical Information Center, 1985
- ↑ 12.0 12.1 Parallel Algorithms: Tree Operations, Guy Blelloch, Carnegie Mellon University, 2009
- ↑ MORIHATA, Akimasa, and Kiminori MATSUZAKI, A Parallel Tree Contraction Algorithm on Non-Binary Trees, MATHEMATICAL ENGINEERING TECHNICAL REPORTS, 2008
- ↑ Parallel Algorithms: Analyzing Parallel Tree Contraction, Guy Blelloch, 2007
- ↑ S Buss, Algorithms for boolean formula evaluation and for tree contraction, Arithmetic, Proof Theory, and Computational Complexity, 1993, pp. 96-115
- ↑ Bader, David A., Sukanya Sreshta, and Nina R. Weisse-Bernstein, Evaluating arithmetic expressions using tree contraction: A fast and scalable parallel implementation for symmetric multiprocessors (SMPs), High Performance Computing—HiPC 2002. Springer Berlin Heidelberg, 2002, pp. 63-75.
- ↑ Applications of Parallel Tree Contraction, Samuel Yeom, 2015
External links
- 6.851: Advanced Data Structures by Prof. Erik Demaine
de:Baum (Graphentheorie)