Difference between revisions of "Segment tree"

Revision as of 09:22, 10 July 2010

The segment tree is a highly versatile data structure, based upon the divide-and-conquer paradigm, which can be thought of as a tree of intervals of an underlying array, constructed so that queries on ranges of the array as well as modifications to the array's elements may be efficiently performed.

Structure

This segment tree.

Suppose that we use the function defined above to evaluate $f(1,N)$ , where $N$ is the number of elements in the array. When $N$ is large, this recursive call has two "children", one of which is the recursive call $f(1,\lfloor\frac{N+1}{2}\rfloor)$ , and the other one of which is $f(\lfloor\frac{N+1}{2}\rfloor+1,N)$ . Each of these children will then have two children of its own, and so on, down until the base case is reached. If we represent these recursive calls with a tree structure, the call $f(1,N)$ would be the root, it would have two children, each child would have two more children, and so on; the base cases would be the leaves of the tree. We are now ready to specify the structure of the segment tree:

it is a binary tree which represents some underlying array;
each node is associated with some interval of the array and contains the value(s) of one or more functions of the elements in that interval;
the root node is associated with the entire array (i.e. the interval $[1,N]$ );
each leaf is associated with an individual element;
each non-leaf node has two children whose associated intervals are disjoint, and the union of the intervals associated with the two children is the interval associated with the parent;
each child's interval has approximately half the size of the parent's interval;
the data stored in each non-leaf node is not only a function of the elements in its associated interval but also a function of the data stored in its children.

That is, the structure of the segment tree is exactly that of the recursive call tree of the function $f(1,N)$ defined above (up to some minor details, such as whether the middle element of an interval of odd size belongs to the left child interval or the right).
Hence, for example, the root node of the array $[9,2,6,3,1,5,0,7]$ would contain the number 0: the minimum of the entire array. Its left child would contain the minimum of $[9,2,6,3]$ , that is, 2, and the right would contain the minimum of $[1,5,0,7]$ , that is, 0. Evidently, the datum stored in the root node, 0, is the minimum of the data stored at its immediate children, 0 and 2. Similarly, the datum 2 stored in the left child is the minimum of the data stored at its immediate children, 2 and 3 (which are the minima in $[9,2]$ and $[6,3]$ , respectively). Each individual element in the array has a leaf node, which merely contains that element itself. (The minimum of a range containing a single element is, of course, the element itself.)

Operations

There are three fundamental operations on a segment tree.

Construction

To construct a segment tree is to initialize it so that it represents some array; we can query and modify this array later but first we need to construct a valid segment tree. We can construct a segment tree either top-down or bottom-up. Top-down construction is recursive; we attempt to fill in the root node, which leads to a recursive call on each of the two children, and so on; the base cases are the leaf nodes, which may be immediately filled in with the corresponding values from the array. After the two recursive calls made on behalf of a non-leaf node return, that node's value is set to the minimum of the values stored at the children. Bottom-up construction is left as an exercise for the reader; the difference in speed is likely to be negligible.

Update

0 is changed to 8 in the segment tree above. Yellow nodes are those on the path from the root to the modified leaf; they are recursively traversed and their data are modified as the update occurs.

To update a segment tree is to modify the value of one element in the underlying array. To do so, we first modify the corresponding leaf node. The other leaf nodes are not affected, since each leaf node is associated only with an individual element. The modified node's parent is affected, since its associated interval contains the modified element, and so is the grandparent, and so on up to the root, but no other nodes are affected. To execute a top-down update, we begin by requesting an update of the root; this leads to a recursive call of one of the two children --- the one whose associated interval contains the modified element. (The other child and its subtree are unaffected.) The same is done recursively for the child; the base case is the leaf node associated with the element to update. After a recursive call has completed, the value in a non-leaf node is re-evaluated as the minimum of the values at its two children. Updating a segment tree bottom-up is possible too and is again left as an exercise for the reader.

Query

Only the two nodes marked in yellow must be accessed to find the minimum of the elements corresponding to the leaves marked in grey.

To query a segment tree is to use it to determine a function of a range in the underlying array (in this case, the minimum element of that range). The execution of a query is more complex than the execution of an update and will be illustrated by example. Suppose we wish to know the minimum element between the first and sixth, inclusive. We shall represent this query as $f(1,6)$ . Each node in the segment tree contains the minimum in some interval: for example, the root node contains $f(1,8)$ , its left child $f(1,4)$ , its right $f(5,8)$ , and so on, with each leaf containing $f(x,x)$ for some $x$ . There is no node that contains $f(1,6)$ , but we notice that $f(1,6) = \min(f(1,4),f(5,6))$ , and that there are nodes in the segment tree containing those two values (shown in yellow). (This expression for $f(1,6)$ is not the one given by the definition of $f$ , but it is fairly clear that $f(x,y) = \min(f(x,z),f(z+1,y))$ where $x \le z < y$ , regardless of the actual value of $z$ .)
When querying a segment tree, therefore, we select a subset of the nodes with the property that the union of their sets of descendants is exactly the interval whose minimum we seek. To do so, we start at the root and recurse over nodes whose corresponding intervals have at least one element in common with the query interval. Hence, in our example of $f(1,6)$ , we notice that both the left and right subtrees contain descendants that contain elements in the query interval; hence, we recurse on both. The left child of the root serves as a base case, since its interval is contained entirely within the query interval; hence it is chosen (marked in yellow). At the right child of the root, we notice that its left child has descendants in the query interval, but its right does not; hence we recurse on the former and not the latter. The former is now another base case, and it too is chosen and marked in yellow. The recursion has now terminated, and the desired minimum is the minimum of the nodes that have been chosen.

Analysis

Some important performance bounds for segment trees are as follows:

Space

A node whose depth is $d$ is associated with an interval of size no greater than $\left\lceil{\frac{N}{2^d}}\right\rceil$ . This is easily proven by induction on $d$ . We therefore see that all nodes at depth $\lceil\lg N \rceil$ are responsible for no more than one element: they are leaves. Hence, a segment tree is a fully balanced binary tree; it has the minimum possible height for the number of nodes it contains. Because of this, we may store the tree as a breadth-first traversal of its nodes, where the left child is (by convention) explored before the right. The segment tree above for $[9,2,6,3,1,5,8,7]$ , then, would be stored as $[1,2,1,2,3,1,7,9,2,6,3,1,5,8,7]$ . Suppose the root is stored at index 1. Then, if a non-leaf node has index $i$ , its left child will be stored at $2i$ and its right at $2i+1$ . Note that there will be some wasted space on the lowest level of the tree; this however is generally acceptable. A binary tree with a depth of $\lceil\lg N\rceil$ may have up to $2^{\lceil\lg N\rceil + 1}-1$ elements, so some mathematical analysis shows that a segment tree on an array of $N$ elements requires an array of no more than $4N$ elements for its own storage in this manner. Asymptotically, a segment tree uses $\mathcal{O}(N)$ space.