Next Article in Journal
Convolution, Correlation, and Uncertainty Principles for the Quaternion Offset Linear Canonical Transform
Next Article in Special Issue
The Square of Some Generalized Hamming Graphs
Previous Article in Journal
Partial Slip Effects for Thermally Radiative Convective Nanofluid Flow
Previous Article in Special Issue
Discrete Integral and Discrete Derivative on Graphs and Switch Problem of Trees
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Efficient and Effective Directed Minimum Spanning Tree Queries

Cyberspace Institute of Advanced Technology, Guangzhou University, Guangzhou 511400, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Mathematics 2023, 11(9), 2200; https://doi.org/10.3390/math11092200
Submission received: 3 April 2023 / Revised: 1 May 2023 / Accepted: 2 May 2023 / Published: 6 May 2023
(This article belongs to the Special Issue Advances in Graph Theory: Algorithms and Applications)

Abstract

:
Computing directed Minimum Spanning Tree ( D M S T ) is a fundamental problem in graph theory. It is applied in a wide spectrum of fields from computer network and communication protocol design to revenue maximization in social networks and syntactic parsing in Natural Language Processing. State-of-the-art solutions are online algorithms that compute  D M S T  for a given graph and a root. For multi-query requirements, the online algorithm is inefficient. To overcome the drawbacks, in this paper, we propose an indexed approach that reuses the computation result to facilitate single and batch queries. We store all the potential edges of  D M S T  in a hierarchical tree in  O ( n )  space complexity. Furthermore, we answer the  D M S T  query of any root in  O ( n )  time complexity. Experimental results demonstrate that our approach can achieve a speedup of 2–3 orders of magnitude in query processing compared to the state-of-the-art while consuming O(n) index size.

1. Introduction

Finding a  D i r e c t e d M i n i m u m S p a n n i n g T r e e ( D M S T ) , which is also known as the Minimum Cost Arborescence problem, in a given directed graph is one of the fundamental problems in graph theory. For a directed graph G, given root r, the query aims to find a  D M S T  rooted at r that connects all the vertices.
D M S T  can be applied to many fields. In communication, it is implemented for minimum cost connectivity [1,2] and control [3]. In the database, it is utilized for reachability queries [4,5]. In natural language processing, it is a classic dependency parsing algorithm [6,7,8,9]. In visualization, it captures information [10,11] and facilitates genetic analysis [12,13,14]. In social networks, it evaluates influence [15] and information flows [16,17].
Chu-Liu and Edmonds [1,18] proposed an algorithm that answers the query in  O ( m n )  time complexity, where n is the number of vertices in the graph and m is the number of edges. Their algorithm is a two-phase algorithm that involves contraction and expansion. Tarjan [19] then proposed a faster implementation in  O ( m log n )  time for sparse graphs and  O ( n 2 )  for dense ones. Gabow et al. [20] improved their algorithms by implementing a Fibonacci heap. By utilizing the  O ( 1 ) u n i t e o p e r a t i o n  feature and a depth-first strategy, their algorithm is able to return a  D M S T  in  O ( m + n log n )  time.
Motivation. All of the above algorithms are online algorithms. For a series of queries on the same graph with different roots, they have to repeat the computing procedure for every single root in the graph. However, we find that for different roots, the corresponding  D M S T s always contain edges of minimum weight related to vertices. We build an index that chooses and stores the edge of the minimum weight for each vertex. Therefore, the computing time for the  D M S T  query can be saved by only referring to the edges in the index instead of searching the original graph.
Our Idea. In this paper, we propose an index-based approach for the  D M S T  problem rooted at any vertex on the given directed graph. We first keep track of all potential edges of any  D M S T  in the given graph and store them in a hierarchical tree. Then we find the  D M S T  of the given root by choosing edges in the tree index. The index computing time is  O ( m n ) , which is the same time cost as [1,18]. The query time for any given root is  O ( n ) . Furthermore, the space complexity is  O ( n ) .
The contribution of this paper can be summarized as follows:
  • We are the first to propose an efficient indexed approach for the  D M S T  problem. We can answer a single query at any root in the graph in  O ( n )  time. Furthermore, we can process batch queries even faster.
  • We prove the correctness of our algorithms. Furthermore, we prove that both single and batch algorithms take worst-case  O ( n )  space and time complexity.
  • We conduct our algorithms on different directed graph datasets to show the efficiency and effectiveness of our algorithms.

2. Related Work

D M S T  can be applied to a wide spectrum of fields. Since the directed maximum spanning tree can be found with the same algorithm by trivially negating the edge cost, we show its applications as follows:
Communication D M S T  of a communication network means the lowest cost way to propagate a message to all other nodes in G [1]. To address the connectivity issue in heterogeneous wireless networks, N. Li et al. [2] proposed a localized  D M S T  algorithm that preserves the network connectivity. The problem of containment control with the minimum number of leaders can be converted into a directed minimum spanning tree [3].
Database Management. To efficiently answer reachability queries, Jin et al. [5] created a path-graph and formed a path-tree cover for the directed graph by extracting a maximum directed spanning tree. Further, they introduced a novel tree-based index framework that utilizes the directed maximal weighted spanning tree algorithm and sampling techniques to maximally compress the generalized transitive closure for the labeled graphs [4].
Natural Language Processing. To generalize the projective parsing method to non-projective languages, McDonald et al. [7] formalized weighted dependency parsing as searching for maximum spanning trees in directed graphs. Furthermore, Smith et al. [8] found a directed maximum spanning tree for maximum  a p o s t e r i o r i  decoding. Moreover, in [9], the authors couched statistical sentence generation as a spanning tree problem and found the  D M S T  of a dependency tree with maximal probability. Yang et al. [6] used  D M S T  as a tool in the evaluation of the induced structure of their proposed structured summarization model.
Visualization. For the purpose of capturing 3D surface measurements with structured light, Brink et al. [10] considered all connections and adjacencies among stripe pixels in the form of a weighted directed graph and indexed the patterns by a maximum spanning tree. Mahshad et al. [11] applied  D M S T  in handwritten math formula recognition.  D M S T  is implemented in GrapeTree [12] to visualize genetic relations and helps genetic analysis in lineage tracing [13] and cancer evolution [14].
Social Networks. For social network hosts to achieve maximum revenue in viral marketing, Zehnder [15] extracted a  D M S T  to generate a most influential tree, which approximates a social network while preserving the most influential path. To contrast the spread of misinformation in online social networks, Marco et al. [16] modeled source identification as a maximal spanning branching problem. Furthermore, Peng et al. [17] extracted important information flows and the hierarchical structure of the networks with  D M S T .

3. Problem Statement

Let  G = ( V , E )  be a directed graph where  V ( G )  is the set of vertices and  E ( G )  is the set of arcs. Arc is a directed edge starting from u to v. We use V and A to denote  V ( G )  and  E ( G ) , and  n = | V |  and  m = | E |  to denote the number of vertices and arcs in the directed graph. We use edge to denote arc if the context is clear. We use  e = u , v  to denote the arc, where u is the tail of the arc, denoted as  t a i l ( e ) , and v is the head, denoted as  h e a d ( e ) . We use  i n - e d g e  to denote any arc incident on a vertex as the head and  o u t - e d g e  to denote any arc direct away from a vertex as the tail. Each arc  u , v  is associated with a positive cost  ϕ u , v . For each  v V , we use  i n ( v )  to denote all the in-edges of v, and we use  o u t ( v )  to denote all the out-edges of v. A path is a sequence of vertices  p = { v 1 , v 2 , , v e } , where  v i , v i + 1 E  and  v i , v j p v i v j . The summary of notations is in Table 1. If there is no path from a vertex to any other vertices in the graph, the  D M S T  rooted at the vertex will cost infinity, which is meaningless. Therefore, for simplicity, in the rest of this paper, we assume that G is strongly connected. As we can obtain the maximum directed spanning tree by simply negating each edge, we focus on answering the minimum one.
Problem Definition. Given a directed graph  G = ( V , E )  and a root vertex  r V , a  D i r e c t e d S p a n n i n g T r e e  is an acyclic subgraph T of G that has all the vertices of G. For each vertex  v T  and  v r :
(1) There is a path starting from r to v.
(2) v has only one in-edge.
D i r e c t e d M i n i m u m S p a n n i n g T r e e  is such a  D i r e c t e d S p a n n i n g T r e e  of the minimum total edge cost.
Example 1. 
In Figure 1, we show a directed graph G with 7 vertices and 15 edges. In Figure 1a, we number each edge of G. Furthermore, in Figure 1b, we show the cost of each edge. In Figure 1c, given root  v 1 , we show the directed minimum spanning tree  T v 1  of G rooted at  v 1 . Furthermore, its total cost is 29. For each vertex  v T v 1 v v 1  there is a path from  v 1  to v, and there is only one in-edge of each vertex in  T v 1  except  v 1 .

4. Existing Solution

In this section, we take a review of Chu-Liu and Edmonds’ algorithm (CLE) that runs in  O ( m n )  time, and our indexed approach is based on the observations. Given a directed graph  G = ( V , E )  and a root vertex  r V , CLE returns the Directed Minimum Spanning Tree  T  of r in a two-phase recursive manner. For each round of recursion,  C L E  chooses the minimum in-edges of each vertex except r, and checks if these in-edges are a cycle. If no cycle is found, then the edges chosen this round including the root vertex r are a  D M S T . Otherwise, the cycles are contracted to a new vertex and the algorithm goes to the next round.
Contraction Phase. For each round, the algorithm first selects minimum in-edges of each vertex  v V \ r , then finds cycles C of the selected edges. Each such cycle  c i C  found this round is contracted to a new vertex  v . All vertices  v  and  v V  and  v C  will be added to the new vertex set  V . The cost of edges  u , v v C  remains the same. The cost of edges  u , v u C v C  will be updated as  ϕ u , v m i n { ϕ ( i n ( v ) ) } . All edges will be added to a new edge set  E , and now we have  G = ( V , E ) . The algorithm finds and contracts cycles in a recursive manner until there are no cycles found, as the in-edge of root r is not considered. Then the algorithm starts to expand T of this round.
Expansion Phase. For the current round, the algorithm starts from root r and breaks cycles. For each  u , v T u c i v c i c i C , the algorithm recovers its original cost, deletes the in-edge  { w , v | w c i }  incident on v in cycle  c i  and adds it to T. Furthermore, the algorithm adds the edges in  c i  except  { w , v | w c i }  to T. The algorithm returns T to the previous round until the final T is found. Here is the framework of the Algorithm 1:
Algorithm 1: CLE  ( G , T , C , r )
Mathematics 11 02200 i001
Example 2. 
We show the contraction and expansion procedure of  C L E  algorithm at root  v 1  in Figure 2. We denote the edge j that is updated in the ith round as  A j i . In (a), for the first round, bold lines are the minimum in-edges of each vertex except the root  v 1 . In (b) and (c),  { v 2 , v 3 }  and  { v 4 , v 5 , v 6 }  are the cycles detected in the first round and contracted to a vertex. Then, the in-edge of each cycle is updated, shown by the dashed line. By now, G is contracted to  G . In (d), for the second round, in-edges of  G  have no cycle. Furthermore, the contraction phase ends. The expansion phase starts from (e). In (e), for graph  G , starting from out-edges of root  v 1 A 3 1  is the only out-edge that starts from the root.  A 3 1  recovers its original cost by adding back cost of  A 8  and breaks the cycle by deleting  A 8 . In (f),  A 13 1  recovers its cost by adding back cost of  A 4  and breaks the cycle by deleting  A 4 . Furthermore, the  D M S T  rooted at  v 1  is found.
By reviewing Chu-Liu and Edmonds’ algorithm, we have the following observations:
Observation 1. 
The cycles detected in each round have a hierarchical structure. CLE contracts and expands the cycles of the graph in a recursive manner, and generates a hierarchical structure naturally. The cycles contracted in the previous round may be a vertex of the cycles contracted in this round. Therefore, cycles of G and  G  in each round can form a hierarchical structure. We can build a hierarchical tree by trivially adapting the contraction phase, detailed in Section 5.
Observation 2. 
In the expansion phase, for each cycle to break, we only need to delete one edge with both ends in the cycle. Every vertex in T has only one in-edge. The cycle  c i C  will be contracted to a vertex  v  and there will be only one minimum in-edge of  v  incident on  v k c i . Furthermore v k  has an in-edge  v k 1 , v k  of  c i . We delete  v k 1 , v k  so that  v k  has only one in-edge.

5. Our Approach

In this section, we will first analyze the drawbacks of the  C L E  algorithm and show the hierarchical tree of the indexed edges. Then, we expand the  D M S T  with reference to the index. Finally, we elaborate on our approach to constructing the index and prove the correctness of our indexed algorithm.

5.1. Hierarchical Tree

Drawbacks of the  C L E  algorithm. We discuss the drawbacks of the  C L E  algorithm in terms of search space and result reuse.
Search space. For each query,  C L E  has to find the minimum in-edge of each vertex and retrieve cycles from the original graph. Therefore,  C L E  suffers from a large search space and spends  O ( m n )  time for the edges and cycles.
Result reuse. For a new query,  C L E  has to restart and recompute the minimum in-edges and detect cycles. The minimum in-edges computed last time are wasted.
Example 3. 
For example, in Figure 3, we show the first round contraction of  C L E  algorithm rooted at  v 7 . In (a), for the first round, we show the minimum in-edges of each vertex in bold lines except root  v 7 . In (b) and (c),  { v 2 , v 3 }  and  { v 4 , v 5 , v 6 }  are the cycles detected in the first round and contracted to a vertex. Then the in-edge of each cycle is updated, shown by the dashed line. Now we compare with the contraction phase in Figure 2 of root  v 1 . It is obvious that edges that form the cycles  { v 2 , v 3 }  and  { v 4 , v 5 , v 6 }  can be reused, though the  D M S T  for root  v 1  and  v 7  choose to delete different edges.
To solve these drawbacks, if we are able to identify the reusable edges and cycles and store them before expansion, we can choose to delete the useless edges for any given root and reduce the search space. Therefore, we propose an index that reuses the edges and cycles in the contraction phase with a trivial adaption and stores the potential edges on the hierarchical tree. By referring to the hierarchical tree, we can retrieve the  D M S T  for any given root instead of searching the entire graph. Therefore, we dramatically reduce the search space for finding the  D M S T . We make a trivial adaption by the contracting root to a vertex, and index the reusable edges and cycles on a hierarchical tree based on the following lemmas:
Lemma 1. 
For a given root r, contracting r to a vertex at any round still generates correct results.
Proof. 
For the final round h of  C L E , as we suppose the graph is strongly connected, we add back the in-edge of the root, and the graph will be contracted to a single vertex. When the expansion starts, the in-edge of the root is removed and the lemma is true. Suppose the lemma is true for round  2 = h ( h 2 ) . For round 1, the root r is contracted to a new vertex  r , since it is true for round 2; therefore, for root  r  the expansion generates the correct result. Now by deleting the in-edge of r, we still obtain the correct result. Therefore, the lemma still holds for any round in which the root is contracted.    □
Example 4. 
We show contraction with root in Figure 4. In (a), for the first round, we find cycle  c 1 = { v 4 , v 5 , v 6 } c 2 = { v 2 , v 3 } , and adapt in-edges of  c 1 , c 2 A 3 A 3 1 A 6 A 6 1 A 11 A 11 1 A 15 A 15 1 A 1 A 1 1 A 9 A 9 1 A 13 A 13 1 . In (b), for the second round, we find cycle  c 3 = { c 1 , v 1 } , and adapt in-edges of  c 3 A 6 1 A 6 2 A 11 1 A 11 2 A 15 1 A 15 2 . In (c), for the third round, we find cycle  c 4 = { c 3 , v 7 } , and adapt in-edges of  c 4 A 6 2 A 6 3 A 11 2 A 11 3 A 12 A 12 3 . Furthermore, next round we will find  c 5 = { c 4 , c 1 } c 5  is not shown since it is only a single vertex.
Lemma 2. 
For any two queries with  r 1  and  r 2 r 1 r 2 , the cycles detected in their contraction phase can be reused.
Proof. 
Since we proved in the previous lemma that contracting root at any round still generates the correct result, we follow the contraction of  r 1  and reuse the cycles detected in its contraction phase. We prove this lemma by contradiction. Suppose that the detected cycles of  r 1  can not be reused. Then for  r 2 , there should be a set of new edges that contracts  r 2  with less cost. However, this contradicts that for each round we find the minimum in-edges of each vertex. Therefore, the lemma holds.    □
For each cycle, it is related to a tree node in the hierarchical tree index H. We denote the correspondent tree node of  c i  as  T N i . Furthermore, we denote the tree node of the first cycle that contracts vertex v as  T N ( v ) ; see also Table 1. Suppose cycle  c i , c j  corresponds to  T N i , T N j . Furthermore,  T N j  is the parent tree node of  T N i . Furthermore, suppose  c i  is contracted to a vertex  v  and  v  is a member of  T N j . In the tree node  T N j v  has an minimum in-edge  u , v  and an out-edge  v , w . Then, edge  u , v  incidents on a member vertex of  T N i , and  v , w  incidents on a member vertex of  T N j . Therefore, we link the child tree node and its parent tree node. Furthermore, we denote the vertices and edges that link child tree nodes and their parents as  l i n k i n g v e r t i c e s  and  l i n k i n g e d g e s .
Example 5. 
In Figure 5, we show the hierarchical tree of Example 4. Furthermore, we mark the linking edge and linking vertex between the child tree node and its parent. In the first round, we detect cycles  c 1 = { v 4 , v 5 , v 6 }  and  c 2 = { v 2 , v 3 } . Then, we build tree node  T N 1  and  T N 2 . In the second round, we detect cycle  c 3 = { v 1 , c 1 } . Furthermore, we build tree node  T N 3 c 1  is a member of  T N 3 , the in-edge of  c 1  incident  c 1  on  v 5  and the out-edge of  c 1  is  A 2 . Therefore, we link  T N 1  with  T N 3  by linking edges  A 3 1  and  A 2  and linking vertices  v 5  and  v 1 . Repeat this procedure on each tree node and we have H in Figure 5.

5.2. Expansion

Now that we have all the potential edges of  D M S T  at an arbitrary root in graph G in H, we need to find all the edges of T for any given root. However, the vertices and cycles are in a hierarchical index structure. We have to design an order of expansion so that we can recover the  D M S T  of any given root and ensure correctness.
Suppose there is only one tree node in the hierarchical tree. For any given root r, we just need to delete the in-edge of r, break the cycle, and find the  D M S T  rooted at r. For more tree nodes in H, the cycles are organized hierarchically. We can expand r following Lemma 1. We first locate the tree node of r T N ( r ) , then break  T N ( r )  as we do with only one tree node. If  T N ( r )  is not the root node of the hierarchical tree, we find the parent tree nodes of  T N ( r )  along the linking edges and expand each parent tree node by regarding the linking vertex as the new root. If  T N ( r )  has child tree nodes, we find its child tree nodes along the linking edges and expand each child tree node by regarding the linking vertices as the new root.
Example 6. 
Given root  v 1 , we first locate its tree node  T N ( v 1 ) , delete  A 2  and break  c 1  and delete  A 8 . Then  c 3  in  T N 4  is the root in round 3, and we delete  A 15 2 . Then  c 4  in  T N 5  is the root in round 4, we delete  A 11 3 , break  c 2  and delete  A 4 . Now we have  T v 1  of total weight 29.
Here is how we expand with reference to the index in detail:
Delete the in-edge of the root. For any given root vertex r in H, we first locate the tree node  T N ( r )  it is in.  T N ( r )  is a cycle. Starting from r, we traverse along the out-edges of each vertex in the cycle and delete the in-edge of the root r in the cycle to break the cycle.
Locate new roots. In the traversing procedure in  T N ( r ) , when meeting a linking edge  e l  and linking vertex  v l , we treat  v l  as the new root in the child tree node  T N ( v l )  and repeat the procedure in  T N ( v l ) . Then we go up to the parent of  T N ( r )  along the linking edge, treat the linking vertex as the new root, and repeat the procedure.
In Algorithm 2, starting from root v and tree node  t n v  (line 3), we first break the cycle by removing the in-edge of root vertex v (line 5). We put all the linking edges in  T N ( v )  in queue Q (line 8) and add them to T (line 9). If  T N ( v )  is not the root node (line 10), we put its parent tree node  T N ( p v )  and its linking vertex  v  to queue Q (line 12). From the view of tree H, in each round we put all its neighbors, its child nodes, and parent node in the queue. It traverses the tree in a BFS manner.
Algorithm 2: Expansion  ( H , r )
Mathematics 11 02200 i002
For the  F i n d C y c l e  procedure, shown in Algorithm 3, we use all the minimum in-edges  M I  found this round and return the set of cycles. For each  v M I  (line 2), we start traversing backward along the minimum in-edge of v, and dye the vertices met with color i (lines 13–14). If we encounter any vertex with the same color of i, then a cycle is found (lines 4–5). We put all the vertices of color i in the cycle  c i  (lines 7–10). Then we start from the next vertex in  M I  until all the vertices in  M I  are visited.
Algorithm 3: FindCycles  ( M I )
Mathematics 11 02200 i003
Example 7. 
For the original graph in Figure 1, we have the hierarchical tree H in Figure 5. We show how we search for  T v 1  in the expansion phase. For root  v 1 , we add  A 3  to T then we enqueue  { T N 1 , v 5 }  and  { T N 4 , C 3 } . Next round, for  { T N 1 , v 5 }  we add  A 7  and  A 10  to T. Next round, for  { T N 4 , C 3 }  we add  A 14  to T then enqueue  { T N 5 , C 4 } . Next round, for  { T N 5 , C 4 }  we add  A 13  to T and enqueue  { T N 2 , v 3 } . Next round, for  { T N 2 , v 3 }  we add  A 5  to T. By now, we have obtained the  D M S T T v 1  at root  v 1 .
Theorem 1. 
Algorithm 2 correctly computes the  D M S T  at given root r.
Proof. 
Firstly, our algorithm traverses the hierarchical tree in a BFS manner. If BFS can traverse the entire tree, so can our algorithm. Therefore, our algorithm breaks all the cycles and returns a tree. Secondly, we prove this by contradiction. Suppose our algorithm returns  T  larger than T. This indicates that some edges of  T  are not the minimum in-edges of vertices in  T . It contradicts the fact that edges in H are the minimum in-edge related to each vertex found in each round in H. Therefore, our algorithm correctly computes the  D M S T  at root r.    □
Theorem 2. 
The query time of the expansion phase is  O ( n ) .
Proof. 
There are at most  n 1  cycles and at most  2 ( n 1 )  edges in H. We have to delete at most  n 1  edges and add at most  n 1  edges to obtain  D M S T . We traverse H in  O ( n )  time and add edges to T in  O ( n 1 )  time. Therefore, the query time of the expansion phase is  O ( n ) .    □

5.3. Hierarchical Tree Construction

The construction of the hierarchical tree is detailed as follows:
Contract the root. We do not specifically exclude the root vertex from the contraction phase. As we assume the graph is strongly connected, the graph will finally contract to a single vertex. The algorithm still generates the correct result as proved in Lemma 1.
Store potential tree edges. To reuse the edges and cycles in the contraction phase, we store all the edges that will be tree edges of any root. We store every edge of each cycle when we contract, and delete the edges that will not be in the tree when we query.
For each round of contraction, we select minimum in-edges of each vertex, find cycles, contract cycles into vertices and update corresponding in-edges in a recursive manner. The cycles are naturally hierarchical (Observation 1). We, therefore, build a hierarchical tree H with each tree node corresponding to a cycle. Then by linking child tree nodes and their parent tree nodes, we construct a tree H of cycles. Here, we build the hierarchical tree.
Find cycles. For round i, we first choose the minimum in-edge of each vertex, then we find cycles of this round.
Build the tree node. The cycles found in round 1 will be leaf nodes. For round i, we store all vertices in cycle  c i C  found this round into tree node  T N j . Furthermore, we contract the cycle to a new vertex.
Build Tree H. For round i, if  T N j  is not a leaf node, we link it with its child nodes by linking edges with linking vertices. Finally, we build the hierarchical tree of cycles H.
We introduce our contraction algorithm in Algorithm 4. For each round, we find minimum in-edges incident on v and find cycles. If no cycle is found this round, we have contracted the graph into a single vertex (lines 3–8). For every in-edge whose head is a vertex in cycles and a tail out of cycles found this round, we update their cost and add them to the new edge set  E  (lines 9–11). After that, we contract cycles into a vertex and add it to the new vertex set  V  (lines 12–13). We then add the vertices not in cycles and their minimum in-edge to the new graph  G , update the graph G, and put cycles to corresponding tree nodes (lines 14–18). Finally, the algorithm returns the hierarchical tree of the cycles.
Lemma 3. 
We have to contract at most  n 1  cycles in the directed graph G.
Proof. 
Lemma 3 is true for  | V | = 1 . Suppose it is true for  | V | = n 1 , and at most  n 2  cycles are contracted. When  | V | = n  and  V = { v 1 , . . . , v n 1 , v n } , we first pick any  n 1  vertices  V = { v 1 , . . . , v n 1 } V  can contract to a vertex  v , and at most  n 2  cycles are contracted. As the graph will finally contract to a single vertex, the contracted vertex  v  and the left vertex  v n  can contract to a cycle. Therefore, Lemma 3 holds.    □
Lemma 4. 
We have to store at most  2 ( n 1 )  potential tree edges in the hierarchical tree H.
Proof. 
According to Lemma 3, there will be at most  n 1  cycles in the directed graph. Each cycle has at least two edges. Therefore, we need to store at most  2 ( n 1 )  edges.    □
Lemma 5. 
We have to delete at most  n 1  edges from the hierarchical tree H to obtain  D M S T .
Proof. 
For the directed graph G with n vertices, the  D M S T  contains  n 1  edges. The  D M S T  contains n vertices each vertex has an in-edge except the root vertex, and the  D M S T  has  n 1  edges. According to Lemma 4, we store at most  2 ( n 1 )  edges in the hierarchical tree H, therefore we have to delete at most  n 1  edges from H.    □
Algorithm 4: Contraction  ( G )
Mathematics 11 02200 i004

6. Batch Query

In this section, we process a sequence of query vertices in a batch by utilizing the unaffected edges of different query vertices. Furthermore, we discuss the query scheduling to minimize the total query cost.

6.1. Batch Query Processing

For a sequence of query vertices, the two distinct query vertices may share many common edges. If we process each vertex independently, we have to break all the cycles and delete the edges for each query, which is costly.
Example 8. 
In Figure 5, for query vertex  v 4 T v 4  and next query vertex  v 5 T v 5 . The difference between  T v 4  and  T v 5  are only the in-edges of  v 4  and  v 5 . Actually, to obtain  T v 5 , we only need to add  A 7  and delete  A 8  from  T v 4 .
Observation 3. 
We derive this observation from Observation 2 and the expansion phase. For two distinct query vertices  q i  and  q i + 1 , given a cycle c they both have to break in their expansion phase, we identify the new roots  v i , v j c j i  for them q i  breaks the cycle c at vertex  v i , the in-edge of  v i  will be deleted q i + 1  breaks the cycle at  v j  and delete its in-edge. The edge difference of  q i  and  q i + 1  breaking cycle c are the in-edges of  v i  and  v j  in cycle c.
From the child tree node to the parent tree node, there is a linking vertex, and we treat it as the new root to break the parent tree node in the expansion phase if there is a query vertex in the child tree node. Therefore, different vertices in the parent tree node are related to different child tree nodes. Furthermore, vertices in the ancestor tree node are related to a subtree of child nodes.
Lemma 6. 
Given  q i T q i q i + 1 T q i + 1 , the difference between  T q i  and  T q i + 1  are in the subtree rooted at  q i  and  q i + 1 ’s Least Common Ancestors (LCA) in the hierarchical tree H.
Proof. 
We prove this by contradiction.  T N ( q i )  and  T N ( q i + 1 )  have linking vertices in their parent tree nodes. In their LCA tree node  T N l c a , their parent tree nodes have different linking vertices  v i  and  v i + 1  related to  q i  and  q i + 1  as the new roots. Suppose parent of  T N l c a  is  T N p , the linking vertex from  T N l c a  to  T N p  is  v l c a T N p . Suppose an arbitrary edge with a head vertex  v a v l c a  and a tail vertex  v b  in  T N p  is affected. Based on Observation 3, the in-edge of  v a  and in-edge of  v l c a  are affected. It indicates that there is one query vertex from a subtree of  T N p  related to  v a  and one query vertex from a subtree of  T N p  related to  v l c a , which contradicts that both  q i  and  q i + 1  are in the subtree related to  v l c a , and no query vertex is related to  v a .    □
From Lemma 6, we reduce the edges to be updated from the entire  D M S T  to the subtree of the  D M S T s of two distinct query vertices. Furthermore, we have to identify the cycles and edges to be updated. Based on Observation 3, only two edges are affected in each cycle related to the query vertex. Therefore, we only need to decide on the affected cycles and find the root vertex related to the query vertices when we break the cycles.
We identify the affected cycles by traversing along the linking edges and decide the new roots in each cycle by linking vertices. For two query vertices  q i , q i + 1 , we discuss the update from  q i + 1  to  q i  as the operations are symmetric. If  q i , q i + 1  are in the same tree node, we just update their in-edges. Otherwise, in their LCA tree node  T N l c a , we find the in-edge  A l c a  of the new root related to  q i  in  T N l c a , and identify the  h e a d ( A l c a )  as the new root in  T N ( h e a d ( A l c a ) ) . Then we repeat the above procedure with  q i  and  h e a d ( A l c a ) . Finally, we will identify all the affected edges and cycles and report the correct result.
We introduce our batch query algorithm in Algorithm 5. We find the LCA of two query vertices and their corresponding vertices in their LCA (line 3). Then we update the affected edges in their LCA cycle (lines 4–7). We locate the new roots in the next affected cycles (lines 8–9). Furthermore, we process the affected edges of cycles in a recursive way until all affected cycles are traversed (lines 10–13).
Algorithm 5: BatchExpansion  ( T , q i , q i + 1 )
Mathematics 11 02200 i005
Example 9. 
In Figure 5, for query vertex  v 7  and next query vertex  v 1 , the LCA of them is  T N 4 . We first update the edges from  v 1  to  v 7 . In  T N 4 v 7  and  C 3  are affected.  A 14  is added and  A 15  is deleted.  v 7  and  h e a d ( A 14 ) = v 7  are in the same cycle, and the procedure terminates. Then we update the edges from  v 7  to  v 1 . The LCA of  h e a d ( A 15 ) = v 6  and  v 1  is  T N 3 . In  T N 3 v 1  and  C 1  are affected.  A 3  is added and  A 2  is deleted.  v 1  and  h e a d ( A 2 ) = v 1  are in the same cycle, and the procedure terminates. The LCA of  h e a d ( A 15 ) = v 6  and  h e a d ( A 3 ) = v 5  is  T N 1 . In  T N 1 v 5  and  v 6  are affected.  A 10  is added and  A 8  is deleted. By now, we correctly update edges and obtain  T v 1  from  T v 7 .
Theorem 3. 
Algorithm 5 correctly update all the changed edges.
Proof. 
Based on Observation 3 and Lemma 6, we need to process all affected edges of the cycles contained in the sub-tree of the LCA. Algorithm 5 first finds the LCA of two query vertices and then traverses all the cycles in a recursive way. Therefore, all the cycles of the sub-tree and the affected edges are correctly updated. □
Theorem 4. 
Algorithm 5 runs in worst-case  O ( n )  time complexity.
Proof. 
We locate LCA in  O ( 1 )  time [21]. Suppose the LCA of  q i  and  q i + 1  is the root of H, and there are only two edges in each child cycle. Therefore, the worst-case time complexity is the same as re-running a single query by Algorithm 2. □

6.2. Query Scheduling

An optimal query scheduling can minimize the total cost of batch queries. However, obtaining an optimal order of query sequence is costly.
Instead, we adopt a simple heuristic method of query scheduling. We order the query sequence by its proximity. The closer the two vertices in the tree nodes of H, the lower the cost of querying the next root vertex. We traverse H post-order and label the tree node with the traversing order. Then we sort the query sequence by their corresponding traverse order. The total cost is  O ( n + k log k )  time complexity and  O ( n )  space complexity, where k is the size of the query sequence and n is the number of cycles. We evaluate the post-order query scheduling in the experiment.
Theorem 5. 
Post-order query scheduling is a two-approximation of optimal query scheduling.
Proof. 
Suppose the optimal sequence of query scheduling is  S = { q 1 , q 2 , , q k 1 , q k } , and they correspondent to tree nodes  T N ( q 1 ) T N ( q 2 ) , ⋯,  T N ( q k 1 )  and  T N ( q K ) . The optimal scheduling is a path on H that starts from  T N ( q 1 )  and ends at  T N ( q k ) . We denote the path as R. For a Steiner Tree that spans all the corresponding tree nodes in S, we denote the minimum one as  M S T . Both R and  M S T  are connected graphs while  M S T  is the minimum one. Therefore, we have  R M S T . For the post-order traverse on the  M S T , denoted as  P O , edges will be visited at most twice. So we have  M S T 1 2 P O . Therefore, we have  R 1 2 P O , and the post-order query scheduling is a two-approximation of optimal scheduling. □

7. Experiment

All algorithms were implemented in C++ and compiled with GNU GCC 4.4.7. All experiments were conducted on a machine with an Intel Xeon 2.8GHz CPU and 256 GB main memory running Linux (Red Hat Linux 4.4.7, 64bit).
In Table 2, we use the open source direct graph dataset from SNAP (https://snap.stanford.edu/data/, last time accessed on 19 February 2023) and KONECT (http://konect.cc/, last time accessed on 19 February 2023). We extracted the Strongly Connected Component from these directed datasets and removed all the self-loops. If the dataset was unweighted, we assigned random weights to the edges. We randomly generated 1000 queries and show the average cost as query time.
Exp-1. Index size and indexing time. We show the number of edges “EdgeNum” in the hierarchical tree of each dataset, the number of cycles “CycleNum”, and the preprocessing time “Pre” Figure 6. The number of edges and cycles in the hierarchical tree  T H  grows linearly with the size of the graph. For example, for dataset  U A , the number of edges is 2492, and the number of cycles is 1091, and for dataset SP, the number of edges in the tree is 2,210,189, and the number of cycles in the tree is 905,654.
Exp-2. Compare the single query time. In this experiment, we show the single query time of Chu-Liu and Edmonds’ as “CLE”. We also show the single query time of Gabow [20] implemented with a Fibonacci heap as “Gabow”. Furthermore, we show the average single query time of ours as “single” in Figure 7. The processing time of  C L E  increases dramatically as the size of datasets increases, owing to the  O ( m n )  time complexity. Though  G a b o w  shows good performance in relatively small graphs, it suffers from an increasing number of edges when the graph size grows. Meanwhile, our  s i n g l e  query time is linear to the growth of graph size because of  O ( n )  time complexity. For dataset  U A C L E ’s single query time is 0.1100 s.  G a b o w ’s single query time is 0.0070 s. Our single query time is 0.0016 s. For dataset  P G C L E ’s single query time is 3.8600 s.  G a b o w ’s single query time is 0.01100 s. Our single query time is  0.0142 . For dataset  S E C L E ’s single query time is 41.0300 s.  G a b o w ’s single query time is 0.0770 s. Our single query time is 0.04060s. For dataset  W T C L E ’s single query time is  706.92 G a b o w ’s single query time is 0.3640 s. Our single query time is 0.1607 s. For dataset  W B C L E ’s single query time is 3918.85 s.  G a b o w ’s single query time is 660.6230 s. Our single query time is 0.4851 s. For dataset  S P C L E ’s single query time is 251,089 s.  G a b o w ’s single query time is 1371.1320s. Our single query time is 2.0151 s. We can see that our single query time shows better performance than the online ones.
Exp-3. Compare single and batch query time. In this experiment, we compare the performance of our single query and batch query algorithms. The single query time grows linearly with the increase in the graph size. In the meantime, the batch query time is affected by the size and structure of the graph. As the batch query is related to the size of cycles of each dataset, more cycles in the graph indicate a greater time cost in query time. Though a larger graph size indicates more cycles, the number of cycles is related to the structure of the graph. Despite this, the batch query time is at least an order of magnitude faster than the single query.
Exp-4. Query scheduling of batch query. To evaluate the effect of query scheduling on the batch query. We conduct our experiment on different order schemes. We show the average number of updated edges of the random order as “Rand”, the average number of updated edges of scheduling by the node ID as “Node”, and the average number of updated edges of post-order traverse as “Post”. Furthermore, we construct a relatively worse case by selecting the next query vertex with less proximity to the previous one. We denote such a worse case as “Worse”. The results are shown in Figure 8. For all the datasets, random query scheduling updates a similar number of edges as the constructed relatively worse case. Furthermore, post-order scheduling performs slightly better than the sequence ordered by node ID. Both the post-order and node ID order are better than random order and the worse case. The good performance of scheduling by the node ID shows the close relationship between cycles during the construction of hierarchical tree H (Figure 9).

8. Conclusions

We propose an indexed approach to answer the  D i r e c t e d M i n i m u m S p a n n i n g T r e e  query of any root. We first pre-process the directed graph in  O ( m n )  time. In the procedure, we build a hierarchical tree that stores all the edges of potential  D M S T  with a space complexity of  O ( n ) . In the expansion phase, starting from the given root, we traverse all of its out-edges on H in a BFS manner to obtain the  D M S T . Then, we propose a batch expansion algorithm by utilizing the shared edges of two query vertices. The time complexity of both expansion algorithms is  O ( n ) .

Author Contributions

Methodology, D.O.; Writing—original draft, Z.W.; Writing—review & editing, Y.W., Q.L. and Z.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Guangzhou Research Foundation grant number SL2022A04J01445 and 202201020165.

Data Availability Statement

Not applicable.

Conflicts of Interest

The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Edmonds, J. Optimum branchings. J. Res. Natl. Bur. Stand. B 1967, 71, 233–240. [Google Scholar] [CrossRef]
  2. Li, N.; Hou, J.C. Topology control in heterogeneous wireless networks: Problems and solutions. In Proceedings of the IEEE INFOCOM 2004, Hong Kong, China, 7–11 March 2004; Volume 1. [Google Scholar]
  3. Gao, L.; Zhao, G.; Li, G.; Liu, Y.; Huang, J.; Deng, L. Containment control of directed networks with time-varying nonlinear multi-agents using minimum number of leaders. Phys. A Stat. Mech. Its Appl. 2019, 526, 120859. [Google Scholar] [CrossRef]
  4. Jin, R.; Hong, H.; Wang, H.; Ruan, N.; Xiang, Y. Computing label-constraint reachability in graph databases. In Proceedings of the 2010 ACM SIGMOD International Conference on Management of Data, Indianapolis, IN, USA, 6–10 June 2010; pp. 123–134. [Google Scholar]
  5. Jin, R.; Xiang, Y.; Ruan, N.; Wang, H. Efficiently answering reachability queries on very large directed graphs. In Proceedings of the 2008 ACM SIGMOD International Conference on Management of Data, Vancouver, DC, Canada, 9–12 June 2008; pp. 595–608. [Google Scholar]
  6. Liu, Y.; Titov, I.; Lapata, M. Single Document Summarization as Tree Induction. In Proceedings of the 2019 Conference of the North American Chapter of the Association for Computational Linguistics: Human Language Technologies, NAACL-HLT, Minneapolis, MN, USA, 2–7 June 2019; Association for Computational Linguistics: Toronto, ON, Canada, 2019; pp. 1745–1755. [Google Scholar]
  7. McDonald, R.; Pereira, F.; Ribarov, K.; Hajic, J. Non-projective dependency parsing using spanning tree algorithms. In Proceedings of the Human Language Technology Conference and Conference on Empirical Methods in Natural Language Processing, Vancouver, DC, USA, 6–8 October 2005; pp. 523–530. [Google Scholar]
  8. Smith, D.A.; Smith, N.A. Probabilistic models of nonprojective dependency trees. In Proceedings of the 2007 Joint Conference on Empirical Methods in Natural Language Processing and Computational Natural Language Learning (EMNLP-CoNLL), Prague, Czech Republic, 28–30 June 2007; pp. 132–140. [Google Scholar]
  9. Wan, S.; Dras, M.; Dale, R.; Paris, C. Improving grammaticality in statistical sentence generation: Introducing a dependency spanning tree algorithm with an argument satisfaction model. In Proceedings of the 12th Conference of the European Chapter of the ACL (EACL 2009), Athens, Greece, 30 March–3 April 2009; pp. 852–860. [Google Scholar]
  10. Brink, W.; Robinson, A.; Rodrigues, M.A. Indexing Uncoded Stripe Patterns in Structured Light Systems by Maximum Spanning Trees. In Proceedings of the BMVC, Leeds, UK, 1–4 September 2008; Citeseer: Princeton, NJ, USA, 2008; Volume 2018, pp. 1–10. [Google Scholar]
  11. Mahdavi, M.; Sun, L.; Zanibbi, R. Visual Parsing with Query-Driven Global Graph Attention (QD-GGA): Preliminary Results for Handwritten Math Formula Recognition. In Proceedings of the 2020 IEEE/CVF Conference on Computer Vision and Pattern Recognition, CVPR Workshops 2020, Seattle, WA, USA, 14–19 June 2020; pp. 2429–2438. [Google Scholar]
  12. Zhou, Z.; Alikhan, N.F.; Sergeant, M.J.; Luhmann, N.; Vaz, C.; Francisco, A.P.; Carriço, J.A.; Achtman, M. GrapeTree: Visualization of core genomic relationships among 100,000 bacterial pathogens. Genome Res. 2018, 28, 1395–1404. [Google Scholar] [CrossRef] [PubMed]
  13. Horns, F.; Vollmers, C.; Croote, D.; Mackey, S.F.; Swan, G.E.; Dekker, C.L.; Davis, M.M.; Quake, S.R. Lineage tracing of human B cells reveals the in vivo landscape of human antibody class switching. eLife 2016, 5, e16578. [Google Scholar] [CrossRef]
  14. Beerenwinkel, N.; Schwarz, R.F.; Gerstung, M.; Markowetz, F. Cancer evolution: Mathematical models and computational inference. Syst. Biol. 2015, 64, e1–e25. [Google Scholar] [CrossRef] [PubMed]
  15. Zehnder, B. Towards Revenue Maximization by VIRAL marketing: A Social Network Host’s Perspective. Master’s Thesis, ETH, Zürich, Switzerland, 2014. [Google Scholar]
  16. Amoruso, M.; Anello, D.; Auletta, V.; Cerulli, R.; Ferraioli, D.; Raiconi, A. Contrasting the Spread of Misinformation in Online Social Networks. J. Artif. Intell. Res. 2020, 69, 847–879. [Google Scholar] [CrossRef]
  17. Yue, P.; Cai, Q.; Yan, W.; Zhou, W. Information Flow Networks of Chinese Stock Market Sectors. IEEE Access 2020, 8, 13066–13077. [Google Scholar] [CrossRef]
  18. Chu, Y.J. On the shortest arborescence of a directed graph. Sci. Sin. 1965, 14, 1396–1400. [Google Scholar]
  19. Tarjan, R.E. Finding optimum branchings. Networks 1977, 7, 25–35. [Google Scholar] [CrossRef]
  20. Gabow, H.N.; Galil, Z.; Spencer, T.H.; Tarjan, R.E. Efficient algorithms for finding minimum spanning trees in undirected and directed graphs. Comb 1986, 6, 109–122. [Google Scholar] [CrossRef]
  21. Bender, M.A.; Farach-Colton, M. The LCA problem revisited. In Proceedings of the LATIN 2000: Theoretical Informatics: 4th Latin American Symposium, Punta del Este, Uruguay, 10–14 April 2000; Springer: Berlin/Heidelberg, Germany, 2000; pp. 88–94. [Google Scholar]
Figure 1. (a) The directed graph  G = ( V , E )  (b) the cost of each edge (c) the  D M S T  rooted at  v 1 .
Figure 1. (a) The directed graph  G = ( V , E )  (b) the cost of each edge (c) the  D M S T  rooted at  v 1 .
Mathematics 11 02200 g001
Figure 2. Contraction and expansion for  C L E  rooted at  v 1 . (a) minimum in-edges chosen in the 1st round contraction phase (b,c) detected cycles and updated edges in the 1st round (d) no cycle detected in the 2nd round, and contraction phase ends (e,f) edges chosen and cycles broken in the expansion phase.
Figure 2. Contraction and expansion for  C L E  rooted at  v 1 . (a) minimum in-edges chosen in the 1st round contraction phase (b,c) detected cycles and updated edges in the 1st round (d) no cycle detected in the 2nd round, and contraction phase ends (e,f) edges chosen and cycles broken in the expansion phase.
Mathematics 11 02200 g002
Figure 3. First round contraction for  C L E  rooted at  v 7 . (a) minimum in-edges chosen in the 1st round contraction phase (b,c) detected cycles and updated edges in the 1st round.
Figure 3. First round contraction for  C L E  rooted at  v 7 . (a) minimum in-edges chosen in the 1st round contraction phase (b,c) detected cycles and updated edges in the 1st round.
Mathematics 11 02200 g003
Figure 4. Contract root. (a) cycles detected and contracted in the 1st round contraction phase (b) cycles detected and contracted in the 2nd round (c) cycles detected and contracted in the 3rd round.
Figure 4. Contract root. (a) cycles detected and contracted in the 1st round contraction phase (b) cycles detected and contracted in the 2nd round (c) cycles detected and contracted in the 3rd round.
Mathematics 11 02200 g004
Figure 5. Hierarchical tree.
Figure 5. Hierarchical tree.
Mathematics 11 02200 g005
Figure 6. Indexsize and preprocessing time.
Figure 6. Indexsize and preprocessing time.
Mathematics 11 02200 g006
Figure 7. Compare single query time.
Figure 7. Compare single query time.
Mathematics 11 02200 g007
Figure 8. Compare single and batch query time.
Figure 8. Compare single and batch query time.
Mathematics 11 02200 g008
Figure 9. Average updated edges of scheduling.
Figure 9. Average updated edges of scheduling.
Mathematics 11 02200 g009
Table 1. The summary of notations.
Table 1. The summary of notations.
NotationDefinition
  G = ( V , E ) The directed graph G, vertex set V and arc set E
n, mnumber of vertices n, number of arcs m
  e = u , v arc e, from u to v, where u is tail, v is head
  ϕ u , v cost of arc  u , v
h e a d ( e ) t a i l ( e ) head and tail of the arc e
i n ( v ) o u t ( v ) all the in-edges whose head is v, all the out-edges whose tail is v
  T , T v any  D M S T  of graph G D M S T  rooted at v
Hhierarchical tree of G containing potential  D M S T  arcs
c i T N i T N ( v ) cycle, its corresponding tree node, and tree node containing v
Cset of cycles
Table 2. Dataset.
Table 2. Dataset.
NameAbbrvType#Vertices#Edges
US airportsUAInfrastructure Network140228,032
p2p-Gnutella30PGComputer Network13,37537,942
soc-Epinions1SEOnline Social Network32,223443,506
wiki-TalkWTCommunication Network111,8811,477,893
web-BerkStanWBHyperlink Network334,8564,523,219
soc-PokecSPOnline Social Network1,304,53629,183,654
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Wang, Z.; Ouyang, D.; Wang, Y.; Liang, Q.; Huang, Z. Efficient and Effective Directed Minimum Spanning Tree Queries. Mathematics 2023, 11, 2200. https://doi.org/10.3390/math11092200

AMA Style

Wang Z, Ouyang D, Wang Y, Liang Q, Huang Z. Efficient and Effective Directed Minimum Spanning Tree Queries. Mathematics. 2023; 11(9):2200. https://doi.org/10.3390/math11092200

Chicago/Turabian Style

Wang, Zhuoran, Dian Ouyang, Yikun Wang, Qi Liang, and Zhuo Huang. 2023. "Efficient and Effective Directed Minimum Spanning Tree Queries" Mathematics 11, no. 9: 2200. https://doi.org/10.3390/math11092200

APA Style

Wang, Z., Ouyang, D., Wang, Y., Liang, Q., & Huang, Z. (2023). Efficient and Effective Directed Minimum Spanning Tree Queries. Mathematics, 11(9), 2200. https://doi.org/10.3390/math11092200

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