The level ancestor problem is the problem of preprocessing a given rooted tree T into a data structure that can determine the ancestor of a given node at a given depth from the root of the tree. Here depth of any node in a tree is the number of edges on the shortest path from the root of the tree to the node.
Given tree is represented as un-directed connected graph having n nodes and n-1 edges.
The idea to solve the above query is to use Jump Pointer Algorithm and pre-processes the tree in O( n log n ) time and answer level ancestor queries in O( logn ) time. In jump pointer, there is a pointer from node N to N’s j-th ancestor, for
j = 1, 2, 4, 8, …, and so on. We refer to these pointers as JumpN[i], where
Jumpu[i] = LA(N, depth(N) – 2i).
When the algorithm is asked to process a query, we repeatedly jump up the tree using these jump pointers. The number of jumps will be at most log n and therefore queries can be answered in O( logn ) time.
So we store 2ith ancestor of each node and also find the depth of each node from the root of the tree.
Now our task reduces to find the ( depth(N) – 2i )th ancestor of node N. Let’s denote X as ( depth(N) – 2i ) and let b bits of the X are set bits (1) denoted by s1, s2, s3, …sb.
X = 2(s1) + 2(s2) + … + 2(sb)
Now the problem is how to find 2j ancestors of each node and depth of each node from the root of the tree?
Initially, we know the 20th ancestor of each node is its parent. We can recursively compute 2j-th ancestor. We know 2j-th ancestor is 2j-1-th ancestor of 2j-1-th ancestor. To calculate the depth of each node we use the ancestor matrix. If we found the root node present in the array of the kth element at jth index then the depth of that node is simply 2j but if root node doesn’t present in the array of ancestors of the kth element than the depth of kth element is 2( index of last non zero ancestor at kth row ) + depth of ancestor present at last index of kth row.
Below is the algorithm to fill the ancestor matrix and depth of each node using dynamic programming. Here, we denote root node as R and initially assume the ancestor of root node as 0. We also initialize depth array with -1 means the depth of the current node is not set and we need to find its depth. If the depth of the current node is not equal to -1 means we have already computed its depth.
we know the first ancestor of each node so we take j>=1, For j>=1 ancstr[k][j] = 2jth ancestor of k = 2j-1th ancestor of (2j-1th ancestor of k) = ancstr[ancstr[i][j-1][j-1] if ancstr[k][j] == R && depth[k] == -1 depth[k] = 2j else if ancstr[k][j] == -1 && depth[k] == -1 depth[k] = 2(j-1) + depth[ ancstr[k][j-1] ]
Let’s understand this algorithm with below diagram.
In the given figure we need to compute 1st level ancestor of the node with value 8. First, we make ancestor matrix which stores 2ith ancestor of nodes. Now, 20 ancestor of node 8 is 10 and similarly 20 ancestor of node 10 is 9 and for node 9 it is 1 and for node 1 it is 5. Based on the above algorithm 1st level ancestor of node 8 is( depth(8)-1 )th ancestor of node 8.
We have pre computed depth of each node and depth of 8 is 5 so we finally need to find (5-1) = 4th ancestor of node 8 which is equal to 21th ancestor of [21 ancestor of node 8] and 21th ancestor of node 8 is 20th ancestor of [20th ancestor of node 8]. So, 20th ancestor of [20th ancestor of node 8] is node with value 9 and 21th ancestor of node 9 is node with value 5. Thus in this way we can compute all query in O(logn) time complexity.
Implementation:
C++
// CPP program to implement Level Ancestor Algorithm #include <bits/stdc++.h> using namespace std; int R = 0; // n -> it represent total number of nodes // len -> it is the maximum length of array to hold // ancestor of each node. In worst case, // the highest value of ancestor a node can have is n-1. // 2 ^ len <= n-1 // len = O(log2n) int getLen( int n) { return ( int )( log (n) / log (2)) + 1; } // ancstr represents 2D matrix to hold ancestor of node. // Here we pass reference of 2D matrix so that the change // made occur directly to the original matrix // depth[] stores depth of each node // len is same as defined above // n is total nodes in graph // R represent root node void setancestor(vector<vector< int > >& ancstr, vector< int >& depth, int * node, int len, int n) { // depth of root node is set to 0 depth[R] = 0; // if depth of a node is -1 it means its depth // is not set otherwise we have computed its depth for ( int j = 1; j <= len; j++) { for ( int i = 0; i < n; i++) { ancstr[node[i]][j] = ancstr[ancstr[node[i]][j - 1]][j - 1]; // if ancestor of current node is R its height is // previously not set, then its height is 2^j if (ancstr[node[i]][j] == R && depth[node[i]] == -1) { // set the depth of ith node depth[node[i]] = pow (2, j); } // if ancestor of current node is 0 means it // does not have root node at its 2nd power // on its path so its depth is 2^(index of // last non zero ancestor means j-1) + depth // of 2^(j-1) th ancestor else if (ancstr[node[i]][j] == 0 && node[i] != R && depth[node[i]] == -1) { depth[node[i]] = pow (2, j - 1) + depth[ancstr[node[i]][j - 1]]; } } } } // c -> it represent child // p -> it represent ancestor // i -> it represent node number // p=0 means the node is root node // R represent root node // here also we pass reference of 2D matrix and depth // vector so that the change made occur directly to // the original matrix and original vector void constructGraph(vector<vector< int > >& ancstr, int * node, vector< int >& depth, int * isNode, int c, int p, int i) { // enter the node in node array // it stores all the nodes in the graph node[i] = c; // to confirm that no child node have 2 ancestors if (isNode == 0) { isNode = 1; // make ancestor of x as y ancstr[0] = p; // ifits first ancestor is root than its depth is 1 if (R == p) { depth = 1; } } return ; } // this function will delete leaf node // x is node to be deleted void removeNode(vector<vector< int > >& ancstr, int * isNode, int len, int x) { if (isNode[x] == 0) cout << "node does not present in graph " << endl; else { isNode[x] = 0; // make all ancestor of node x as 0 for ( int j = 0; j <= len; j++) { ancstr[x][j] = 0; } } return ; } // x -> it represent new node to be inserted // p -> it represent ancestor of new node void addNode(vector<vector< int > >& ancstr, vector< int >& depth, int * isNode, int len, int x, int p) { if (isNode[x] == 1) { cout << " Node is already present in array " << endl; return ; } if (isNode[p] == 0) { cout << " ancestor not does not present in an array " << endl; return ; } isNode[x] = 1; ancstr[x][0] = p; // depth of new node is 1 + depth of its ancestor depth[x] = depth[p] + 1; int j = 0; // while we don't reach root node while (ancstr[x][j] != 0) { ancstr[x][j + 1] = ancstr[ancstr[x][j]][j]; j++; } // remaining array will fill with 0 after // we find root of tree while (j <= len) { ancstr[x][j] = 0; j++; } return ; } // LA function to find Lth level ancestor of node x void LA(vector<vector< int > >& ancstr, vector< int > depth, int * isNode, int x, int L) { int j = 0; int temp = x; // to check if node is present in graph or not if (isNode[x] == 0) { cout << "Node is not present in graph " << endl; return ; } // we change L as depth of node x - int k = depth[x] - L; // int q = k; // in this loop we decrease the value of k by k/2 and // increment j by 1 after each iteration, and check for set bit // if we get set bit then we update x with jth ancestor of x // as k becomes less than or equal to zero means we // reach to kth level ancestor while (k > 0) { // to check if last bit is 1 or not if (k & 1) { x = ancstr[x][j]; } // use of shift operator to make k = k/2 // after every iteration k = k >> 1; j++; } cout << L << "th level ancestor of node " << temp << " is = " << x << endl; return ; } int main() { // n represent number of nodes int n = 12; // initialization of ancestor matrix // suppose max range of node is up to 1000 // if there are 1000 nodes than also length // of ancestor matrix will not exceed 10 vector<vector< int > > ancestor(1000, vector< int >(10)); // this vector is used to store depth of each node. vector< int > depth(1000); // fill function is used to initialize depth with -1 fill(depth.begin(), depth.end(), -1); // node array is used to store all nodes int * node = new int [1000]; // isNode is an array to check whether a // node is present in graph or not int * isNode = new int [1000]; // memset function to initialize isNode array with 0 memset (isNode, 0, 1000 * sizeof ( int )); // function to calculate len // len -> it is the maximum length of array to // hold ancestor of each node. int len = getLen(n); // R stores root node R = 2; // construction of graph // here 0 represent that the node is root node constructGraph(ancestor, node, depth, isNode, 2, 0, 0); constructGraph(ancestor, node, depth, isNode, 5, 2, 1); constructGraph(ancestor, node, depth, isNode, 3, 5, 2); constructGraph(ancestor, node, depth, isNode, 4, 5, 3); constructGraph(ancestor, node, depth, isNode, 1, 5, 4); constructGraph(ancestor, node, depth, isNode, 7, 1, 5); constructGraph(ancestor, node, depth, isNode, 9, 1, 6); constructGraph(ancestor, node, depth, isNode, 10, 9, 7); constructGraph(ancestor, node, depth, isNode, 11, 10, 8); constructGraph(ancestor, node, depth, isNode, 6, 10, 9); constructGraph(ancestor, node, depth, isNode, 8, 10, 10); // function to pre compute ancestor matrix setancestor(ancestor, depth, node, len, n); // query to get 1st level ancestor of node 8 LA(ancestor, depth, isNode, 8, 1); // add node 12 and its ancestor is 8 addNode(ancestor, depth, isNode, len, 12, 8); // query to get 2nd level ancestor of node 12 LA(ancestor, depth, isNode, 12, 2); // delete node 12 removeNode(ancestor, isNode, len, 12); // query to get 5th level ancestor of node // 12 after deletion of node LA(ancestor, depth, isNode, 12, 1); return 0; } |
1th level ancestor of node 8 is = 5 2th level ancestor of node 12 is = 1 Node is not present in graph
Ready to dive in? Explore our Free Demo Content and join our DSA course, trusted by over 100,000 neveropen!