A simple algorithm to generate the minimal separators and the maximal cliques of a chordal graph

A simple algorithm to generate the minimal separators and the maximal cliques of a chordal graph

Information Processing Letters 111 (2011) 508–511 Contents lists available at ScienceDirect Information Processing Letters www.elsevier.com/locate/i...

134KB Sizes 8 Downloads 143 Views

Information Processing Letters 111 (2011) 508–511

Contents lists available at ScienceDirect

Information Processing Letters www.elsevier.com/locate/ipl

A simple algorithm to generate the minimal separators and the maximal cliques of a chordal graph ✩ Anne Berry ∗ , Romain Pogorelcnik LIMOS UMR CNRS 6158, Ensemble Scientifique des Cézeaux, F-63 173 Aubière, France

a r t i c l e

i n f o

a b s t r a c t

Article history: Received 15 February 2010 Received in revised form 13 October 2010 Accepted 20 February 2011 Available online 3 March 2011 Communicated by W.-L. Hsu

We present a simple unified algorithmic process which uses either LexBFS or MCS on a chordal graph to generate the minimal separators and the maximal cliques in linear time in a single pass. © 2011 Elsevier B.V. All rights reserved.

Keywords: Graph algorithms LexBFS MCS Minimal separator Maximal clique Moplex ordering

1. Introduction Chordal graphs (which are the graphs with no chordless cycle on four or more vertices), were characterized by Fulkerson and Gross [7] as the graphs for which one can repeatedly find a simplicial vertex (a vertex whose neighborhood is a clique) and remove it from the graph, until no vertex is left; this process, called simplicial elimination, defines an ordering α on the vertices called a perfect elimination ordering (peo). (α (i ) denotes the vertex bearing number i, and α −1 (x) the number of vertex x.) At each step of the elimination process, a new transitory graph is defined. Rose [9] showed that for any given peo α of a chordal graph, any minimal separator of the graph is defined in the course of the simplicial elimination process as the transitory neighborhood of some vertex. Moreover, it is easy to see that in any maximal clique K , the vertex y of smallest ✩ Research supported by the French Agency for Research under the DEFIS program TODO, ANR-09-EMER-010. Corresponding author. E-mail addresses: [email protected] (A. Berry), [email protected] (R. Pogorelcnik).

*

0020-0190/$ – see front matter doi:10.1016/j.ipl.2011.02.013

© 2011

Elsevier B.V. All rights reserved.

number by α defines K as its transitory closed neighborhood. In the rest of the paper, we will call such vertices generators, and our goal will be to detect these generators efficiently. We will use graph search algorithms Lexicographic Breadth-First Search (LexBFS) [10] and Maximum Cardinality Search (MCS) [11], both tailored to generate a peo. The idea behind these algorithms is simple. Both algorithms number the vertices from n to 1 (n is the number of vertices of the graph). Each vertex x bears a label, which is the list of numbers of the neighbors of x with a higher number (LexBFS) or the cardinality of this list (MCS). At each step, a new vertex of maximum label is chosen to be numbered. If the label of this new vertex xi is not larger than the label of the previous vertex xi +1 , then we show that xi generates a minimal separator of the graph. The vertices of this minimal separator are the already numbered neighbors of xi . We also show that the previously numbered vertex, xi +1 , generates a maximal clique: xi +1 , together with its already numbered neighbors (excluding xi ), define a maximal clique of the graph. Thus the labels of algorithms LexBFS and MCS enable the user to detect the generators as soon as they are numbered. The minimal separators and maximal cliques can be

A. Berry, R. Pogorelcnik / Information Processing Letters 111 (2011) 508–511

found with a single pass, without requiring a preliminary pass of the algorithm to number all the vertices. These results have been proved for MCS in two separate papers. Blair and Peyton [5], while studying how MCS defines a clique tree of a chordal graph, showed how to generate the maximal cliques (which are the nodes of the clique tree), using the MCS labels. Kumar and Madhavan [8] showed how to generate the minimal separators of a chordal graph, given an MCS ordering. However, although a minimal separator and a maximal clique generator are actually detected at the same step of the algorithm, these results have not been unified. Moreover, Kumar and Madhavan [8] use a peo as input, but the minimal separators as well as the maximal cliques can be computed during the execution, which may be an important feature when handling very large graphs, since a global pre-numbering of the vertices is not necessary or even useful. LexBFS, as we will show, exhibits the same property as MCS regarding these generators. Our aim is to present a simple unified single-pass algorithm which generates the minimal separators and the maximal cliques of a chordal graph. 2. Preliminaries In the rest of the paper, we will consider all graphs to be connected. Symbol + denotes disjoint union. The neighborhood N G (x) of vertex x in graph G is N G (x) = { y = x | xy ∈ E } The  neighborhood N G ( X ) of a vertex set X ⊂ V is N G ( X ) = x∈ X N G (x) − X . A clique is a set of pairwise adjacent vertices. A clique X ⊆ V is maximal if ∀ y ∈ V − X , X + { y } fails to be a clique. A subset X of vertices is called a module if the vertices of X share the same external neighborhood: ∀x ∈ X , N (x) − X = N ( X ). A subset S of vertices of a connected graph G is called a separator if G ( V − S ) is not connected. A separator S is called an xy-separator if x and y lie in different connected components of G ( V − S ), a minimal xy-separator if S is an xy-separator and no proper subset of S is an xyseparator. A separator S is a minimal separator if there is some pair {x, y } such that S is a minimal xy-separator. Alternately, S is a minimal separator if and only if G ( V − S ) has at least two connected components C 1 and C 2 such that N (C 1 ) = N (C 2 ) = S; such components are called full components of S. 3. Main theorem Definition 3.1. Given a chordal graph G and a peo α of G, for any vertex x, Madj(x) is the set of neighbors of x with a number higher than that of x: Madj(x) = { y ∈ N (x) | α −1 ( y ) > α −1 (x)}. Definition 3.2. Given a chordal graph G and a peo we will call:

α of G,

• minimal separator generator any vertex xi with number i by α , such that Madj(xi ) is a minimal separator of G and label(xi )  label(xi +1 ), where xi +1 = α (i + 1);

509

• maximal clique generator a vertex y such that Madj( y )∪ { y } is a maximal clique of G. Theorem 3.3. Let α be a peo defined by either LexBFS or MCS, let xi be the vertex with number i, let xi +1 be the vertex with number i + 1. a) xi is a minimal separator generator if and only if label(xi )  label(xi +1 ). b) xi +1 is a maximal clique generator if and only if label(xi )  label(xi +1 ) or i + 1 = 1. We will now discuss the proof of Theorem 3.3. We will first present moplex elimination, which is a process that explains how both MCS and LexBFS scan the minimal separators and the maximal cliques of a chordal graph, then prove Theorem 3.3 for LexBFS. 3.1. Moplex elimination Definition 3.4. (See [1].) A moplex is a clique X such that X is a module and N ( X ) is a minimal separator. We extend this definition to a clique whose neighborhood is empty. A simplicial moplex is a moplex whose vertices are all simplicial. Property 3.5. (See [1].) Any chordal graph which is not a clique has at least two non-adjacent simplicial moplexes. From this was derived a variant of the characterization of Fulkerson and Gross for chordal graphs by simplicial elimination of vertices: Characterization 3.6. (See [1].) A graph is chordal if and only if one can repeatedly delete a simplicial moplex until the graph is a clique (which we will call the terminal moplex). We call this process moplex elimination. Note that moplex elimination on a chordal graph is a special case of simplicial elimination, since at each step a set of simplicial vertices is eliminated. Note also that for a connected graph G, the transitory elimination graph obtained at the end of each step remains connected. Moplex elimination defines an ordered partition ( X 1 , X 2 , . . . , X k ) of the vertices of the graph into the successive moplexes which are defined in the successive transitory elimination graphs. We will call this partition a moplex ordering. Theorem 3.7. 1 Let G be a chordal graph, let ( X 1 , X 2 , . . . , X k ) be a moplex ordering of G. At each step i < k of the elimination process finding moplex X i in transitory graph G i ,

• N G i ( X i ) is a minimal separator of G. • X i ∪ N G i ( X i ) is a maximal clique of G. The terminal moplex X k is a maximal clique. There are no other minimal separators or maximal cliques in G. 1

The results proved in this subsection were briefly presented in [2].

510

A. Berry, R. Pogorelcnik / Information Processing Letters 111 (2011) 508–511

To prove this, we will first recall the result from Rose [9]: Property 3.8. (See [9].) Let G be a chordal graph, let α be a peo of G, let S be a minimal separator of G. Then there is some vertex x such that Madj(x) = S, and in every full component C of S, there is some vertex y such that S ⊆ N ( y ) (such a vertex is called a confluence point of C ). Proof of Theorem 3.7. N G i ( X i ) is a minimal separator of G: in transitory graph G i , S i = N G i ( X i ) is by definition a minimal separator of G i , with at least two full components C 1 and C 2 , each containing a confluence point, which we will call x1 and x2 . Suppose S i is not a minimal separator of G. Then there must be a chordless path from x1 to x2 in G which contains no vertex of S i . Let y be the first vertex of this path to be eliminated, at step j < i; y must be simplicial in G j , but this is impossible, since y has two non-adjacent neighbors on the path. Thus every N G i ( X i ) (i < k) is a minimal separator of G, and by Property 3.8, all the minimal separators of G have thus been encountered. X i ∪ N G i ( X i ) = K i is a maximal clique of G: K i is a maximal clique of G i and thus a clique of G. Suppose it is not a maximal clique of G. Let x j be a vertex belonging to the moplex X j ( j < i ), of largest number which can be added to K i to make a larger clique in G. In G j , S j = N G j (x j ) is a minimal separator, which contains all the vertices of K i . But in G j , S j must have a full component C , disjoint from X j ; C must contain a confluence point z, which sees in G j all the vertices of S j , and thus all the vertices of K i . z must belong to a moplex of number larger than j, which contradicts the assumption that X j is of maximum number. The same argument can be given to prove that the terminal moplex is a maximal clique of G. 2 Note that for a given chordal graph G, there may be many different moplex orderings, but there is always the same number of moplexes, since this is the number of maximal cliques of the graph. With any moplex ordering ( X 1 , X 2 , . . . , X k ), we can associate a peo α by processing the moplexes from 1 to n and giving consecutive numbers to the vertices of a given moplex. Using α , we can define the minimal separators and maximal cliques as vertex neighborhoods: Property 3.9. Let G be a chordal graph, let ( X 1 , X 2 , . . . , X k ) be a moplex ordering of G, let α be a peo associated with this moplex ordering. Then in the course of a moplex elimination, at each step i that processes moplex X i ,

• the vertex x of moplex X i whose number is the smallest by α defines a maximal clique {x} ∪ Madj (x) of G; • the vertex y of moplex X i whose number is the largest by α defines a minimal separator Madj( y ) of G. 3.2. Generators of LexBFS LexBFS defines a moplex ordering and an associated peo: [1] showed that LexBFS always numbers as 1 a vertex belonging to a moplex (which we will call X 1 ). They

Fig. 1. Numbers and labels of a LexBFS execution on a chordal graph.

also proved that the vertices of X 1 receive consecutive numbers by LexBFS. These properties are true at each step of LexBFS in the transitory elimination graph. Therefore, LexBFS defines a moplex elimination ( X 1 , . . . , X k ), by numbering consecutively the vertices of X 1 , then numbering consecutively the vertices of X 2 , and so forth: Theorem 3.10. In a chordal graph, LexBFS defines a moplex ordering. Note that it is easy to deduce from [5] that MCS also defines a moplex ordering. Example 3.1. Fig. 1 shows the numbers and labels of a LexBFS execution on a chordal graph. Moplex

Minimal separator

Maximal clique

X1 = { g} X 2 = {a} X 3 = {d, h} X 4 = {c } X 5 = {b, e , f }

{b, c } {b, c , d} {b, c } {b, f }

{b, c , g } {a, b, c , d} {b, c , d, h} {b, c , f } {b, e , f }



Since the vertices of any transitory moplex X i are numbered consecutively, when running LexBFS (numbering the vertices from n to 1), as long as the labels increase, we are defining a moplex, X i . When the labels stop increasing (when numbering vertex y), then we have started a new moplex X i −1 which will contain y. Using Property 3.9, we can deduce that LexBFS and MCS generate the minimal separators and the maximal cliques. This completes the proof of Theorem 3.3. 4. Algorithm We can now derive from Theorem 3.3 a generalized algorithm to generate the minimal separators and the maximal cliques of a chordal graph. In the algorithm below, it is considered that either LexBFS or MCS is used. “Use i to increment the label of y” is translated as “add i to label of y” for LexBFS and as

A. Berry, R. Pogorelcnik / Information Processing Letters 111 (2011) 508–511

“add 1 to the label of y” for MCS. The labels are all considered initialized at the beginning, as ∅ for LexBFS and as 0 for MCS. G NUM ← G NUM + {xi } is shorthand for “V NUM ← V NUM + {xi }; G NUM ← G ( V NUM )”. In the same fashion, G ELIM ← G ELIM − {xi } is shorthand for V ELIM ← V ELIM − {xi }; G ELIM ← G ( V ELIM ); Algorithm Minseps-Maxcliques Input : A chordal graph G = ( V , E ). Output: Set S of minimal separators of G; Set S of minimal separator generators of G; Set K of maximal cliques of G; Set K of maximal clique generators of G; init: G NUM ← G (∅); G ELIM ← G; S ← ∅; S ← ∅; K ← ∅; K ← ∅; for i =n downto 1 do Choose a vertex xi of G ELIM of maximum label; G NUM ← G NUM + {xi }; if i = n and label(xi )  λ then //xi is a min. sep. generator and xi +1 is a max. clique //generator S ← S + {xi }; S ← S ∪ { N G NUM (xi )}; K ← K + {xi +1 }; K ← K + { N G NUM (xi +1 ) ∪ {xi +1 }};

λ ← label(xi ); foreach y ∈ N G ELIM (xi ) do

Use i to increment the label of y;

G ELIM ← G ELIM − {xi }; K ← K + {x1 }; K ← K + { N G (x1 ) ∪ x1 };

511

cliques of a chordal graph in a single pass of either LexBFS or MCS, without requiring the preliminary computation of a peo. Though both LexBFS and MCS yield an optimal lineartime process for this problem, it is important to note that they define a different set of peos of a chordal graph [4], and exhibit different local behaviors [3]. It may be important to use one or the other, depending on the intended application. We will end this paper by noting that two recently introduced search algorithms, LexDFS and MNS [6] fail to behave in the same fashion. We leave open the question as to whether other graph search algorithms are endowed with the properties described in this paper. Acknowledgements The authors thank the referees for their valuable suggestions. References [1] A. Berry, Separability generalizes Dirac’s theorem, Discrete Appl. Math. 84 (1–3) (1998) 43–53. [2] A. Berry, J.-P. Bordat, Moplex elimination orderings, Electron. Notes Discrete Math. 8 (2001) 6–9.

In Fig. 1, S = {{6, 7}, {5, 6}, {4, 5, 6}, {5, 6}}; S = {5, 4, 2, 1}; K = {{6, 7, 8}, {5, 6, 7}, {3, 4, 5, 6}, {2, 4, 5, 6}, {1, 5, 6}}; K = {6, 5, 3, 2, 1}. Note that Madj(h) is a minimal separator, but h = x3 is not a minimal separator generator because label(3) > label(4) with x4 = d. The complexity of the above algorithm is the same as for LexBFS or MCS, which is in O (n + m) time. A minimal separator may be generated several times, depending on the number of full components it defines:

[3] A. Berry, J.R.S. Blair, J.-P. Bordat, G. Simonet, Graph extremities defined by search algorithms, Algorithms 3 (2) (2010) 100–124, doi:10.3390/a3020100.

Property 4.1. Let α be a peo defined by LexBFS or MCS, let S be a minimal separator, let k be the number of full components of S. Then S has k − 1 generators by α .

[8] P. Sreenivasa Kumar, C.E. Veni Madhavan, Minimal vertex separators of chordal graphs, Discrete Appl. Math. 89 (1–3) (1998) 155–168.

[4] A. Berry, R. Krueger, G. Simonet, Maximal label search algorithms to compute perfect and minimal elimination orderings, SIAM J. Discrete Math. 23 (1) (2009) 428–446. [5] J.R.S. Blair, B.W. Peyton, An introduction to chordal graphs and clique trees, Graph Theory and Sparse Matrix Computation 84 (1–3) (1993) 1–29. [6] D.G. Corneil, R. Krueger, A unified view of graph searching, SIAM J. Discrete Math. 22 (4) (2008) 1259–1276. [7] D.R. Fulkerson, O.A. Gross, Incidence matrices and interval graphs, Pacific J. Math. 28 (1969) 565–570.

[9] D.J. Rose, Triangulated graphs and the elimination process, J. Math. Anal. Appl. 32 (3) (1970) 597–609.

5. Conclusion

[10] D.J. Rose, R.E. Tarjan, G.S. Lueker, Algorithmic aspects of vertex elimination on graphs, SIAM J. Comput. 5 (2) (1976) 266–283.

In this paper, we presented a simple and optimal process which generates the minimal separators and maximal

[11] R.E. Tarjan, M. Yannakakis, Simple linear-time algorithms to test chordality of graphs, test acyclicity of hypergraphs, and selectively reduce acyclic hypergraphs, SIAM J. Comput. 13 (3) (1984) 566–579.