UNIFORM RANDOM SAMPLING OF PLANAR GRAPHS IN LINEAR TIME ERIC´ FUSY Abstract. This article introduces new algorithms for the uniform random gen- eration of labelled planar graphs. Its principles rely on Boltzmann samplers, as recently developed by Duchon, Flajolet, Louchard, and Schaeffer. It combines the Boltzmann framework, a suitable use of rejection, a new combinatorial bi- jection found by Fusy, Poulalhon and Schaeffer, as well as a precise analytic description of the generating functions counting planar graphs, which was re- cently obtained by Gim´enez and Noy. This gives rise to an extremely efficient algorithm for the random generation of planar graphs. There is a preprocessing step of some fixed small cost; and the expected time complexity of generation is quadratic for exact-size uniform sampling and linear for approximate-size sampling. This greatly improves on the best previously known time complexity for exact-size uniform sampling of planar graphs with n vertices, which was a little over O(n7). This is the extended and revised journal version of a conference paper with the title “Quadratic exact-size and linear approximate-size random generation of planar graphs”, which appeared in the Proceedings of the International Con- ference on Analysis of Algorithms (AofA’05), 6-10 June 2005, Barcelona. 1. Introduction A graph is said to be planar if it can be embedded in the plane so that no two edges cross each other. In this article, we consider planar graphs that are labelled, i.e., the n vertices bear distinct labels in [1..n], and simple, i.e., with no loop nor multiple edges. Statistical properties of planar graphs have been intensively studied [6, 19, 20]. Very recently, Gim´enez and Noy [20] have solved exactly the difficult problem of the asymptotic enumeration of labelled planar graphs. They also provide exact analytic expressions for the asymptotic probability distribution of parameters such as the number of edges and arXiv:0705.1287v3 [math.CO] 18 Dec 2008 the number of connected components. However many other statistics on random planar graphs remain analytically and combinatorially intractable. Thus, it is an important issue to design efficient random samplers in order to observe the (asymptotic) behaviour of such parameters on random planar graphs. Moreover, random generation is useful to test the correctness and efficiency of algorithms on planar graphs, such as planarity testing, embedding algorithms, procedures for finding geometric cuts, and so on. Denise, Vasconcellos, and Welsh have proposed a first algorithm for the random gen- eration of planar graphs [8], by defining a Markov chain on the set n of labelled planar graphs with n vertices. At each step, two different vertices v and vG′ are chosen at ran- dom. If they are adjacent, the edge (v, v′) is deleted. If they are not adjacent and if the operation of adding (v, v′) does not break planarity, then the edge (v, v′) is added. By Key words and phrases. Planar graphs, Random generation, Boltzmann sampling. 1 2 ERIC´ FUSY Aux. mem. Preproc. time Time per generation Markov chains O(log n) O(log n) unknown {exact size} 5 ∗ 7 3 Recursive method O(n log n) O `n ´ O(n ) {exact size} 2 Boltzmann sampler O((log n)k) O((log n)k) O(n ) {exact size} O(n/ǫ) {approx. size} Figure 1. Complexities of the random samplers of planar graphs (O∗ stands for a big O taken up to logarithmic factors). symmetry of the transition matrix of the Markov chain, the probability distribution con- verges to the uniform distribution on n. This algorithm is very easy to describe but more difficult to implement, as there existsG no simple linear-time planarity testing algorithm. More importantly, the rate of convergence to the uniform distribution is unknown. A second approach for uniform random generation is the recursive method introduced by Nijenhuis and Wilf [25] and formalised by Flajolet, Van Cutsem and Zimmermann [15]. The recursive method is a general framework for the random generation of combinatorial classes admitting a recursive decomposition. For such classes, producing an object of the class uniformly at random boils down to producing the decomposition tree corresponding to its recursive decomposition. Then, the branching probabilities that produce the de- composition tree with suitable (uniform) probability are computed using the coefficients counting the objects involved in the decomposition. As a consequence, this method re- quires a preprocessing step where large tables of large coefficients are calculated using the recursive relations they satisfy. Bodirsky et al. have described in [5] the first polynomial-time random sampler for planar graphs. Their idea is to apply the recursive method of sampling to a well known combinatorial decomposition of planar graphs according to successive levels of connectiv- ity, which has been formalised by Tutte [33]. Precisely, the decomposition yields some recurrences satisfied by the coefficients counting planar graphs as well as subfamilies (con- nected, 2-connected, 3-connected), which in turn yield an explicit recursive way to generate planar graphs uniformly at random. As the recurrences are rather involved, the complex- ity of the preprocessing step is large. Precisely, in order to draw planar graphs with n vertices (and possibly also a fixed number m of edges), the random generator described in [5] requires a preprocessing time of order O n7(log n)2(log log n) and an auxiliary memory of size O(n5 log n). Once the tables have been computed, the complexity of each generation is O(n3). A more recent optimisation of the recursive method by Denise and Zimmermann [9] —based on controlled real arithmetics— should be applicable; it would improve the time complexity somewhat, but the storage complexity would still be large. In this article, we introduce a new random generator for labelled planar graphs, which relies on the same decomposition of planar graphs as the algorithm of Bodirsky et al. The main difference is that we translate this decomposition into a random generator using the framework of Boltzmann samplers, instead of the recursive method. Boltzmann samplers have been recently developed by Duchon, Flajolet, Louchard, and Schaeffer in [11] as a powerful framework for the random generation of decomposable combinatorial structures. The idea of Boltzmann sampling is to gain efficiency by relaxing the constraint of exact-size sampling. As we will see, the gain is particularly significant in the case of planar graphs, where the decomposition is more involved than for classical classes, such as trees. Given a combinatorial class, a Boltzmann sampler draws an object of size n with probability UNIFORM RANDOM SAMPLING OF PLANAR GRAPHS IN LINEAR TIME 3 proportional to xn (or proportional to xn/n! for labelled objects), where x is a certain real parameter that can be appropriately tuned. Accordingly, the probability distribution is spread over all the objects of the class, with the property that objects of the same size have the same probability of occurring. In particular, the probability distribution is uniform when restricted to a fixed size. Like the recursive method, Boltzmann samplers can be designed for any combinatorial class admitting a recursive decomposition, as there are explicit sampling rules associated with each classical construction (Sum, Product, Set, Substitution). The branching probabilities used to produce the decomposition tree of a random object are not based on the coefficients as in the recursive method, but on the values at x of the generating functions of the classes intervening in the decomposition. In this article, we translate the decomposition of planar graphs into Boltzmann samplers and obtain very efficient random generators that produce planar graphs with a fixed number of vertices or with fixed numbers of vertices and edges uniformly at random. Furthermore, our samplers have an approximate-size version where a small tolerance, say a few percents, is allowed for the size of the output. For practical purpose, approximate- size random sampling often suffices. The approximate-size samplers we propose are very efficient as they have linear time complexity. Theorem 1 (Samplers with respect to number of vertices). Let n N be a target size. ∈ An exact-size sampler An can be designed so as to generate labelled planar graphs with n vertices uniformly at random. For any tolerance ratio ǫ> 0, an approximate-size sampler An,ǫ can be designed so as to generate planar graphs with their number of vertices in [n(1 ǫ),n(1 + ǫ)], and following the uniform distribution for each size k [n(1 ǫ),n(1 + ǫ)]−. ∈ − Under a real-arithmetics complexity model, Algorithm An is of expected complexity 2 O(n ), and Algorithm An,ǫ is of expected complexity O(n/ǫ). Theorem 2 (Samplers with respect to the numbers of vertices and edges). Let n N be a target size and µ (1, 3) be a parameter describing the ratio edges-vertices. An exact-size∈ ∈ sampler An,µ can be designed so as to generate planar graphs with n vertices and µn edges uniformly at random. For any tolerance-ratio ǫ > 0, an approximate-size sampler⌊ ⌋ An,µ,ǫ can be designed so as to generate planar graphs with their number of vertices in [n(1 ǫ),n(1 + ǫ)] and their ratio edges/vertices in [µ(1 ǫ),µ(1 + ǫ)], and following the uniform− distribution for each fixed pair (number of vertices− , number of edges). Under a real-arithmetics complexity model, for a fixed µ (1, 3), Algorithm A is of ∈ n,µ expected complexity O (n5/2). For fixed constants µ (1, 3) and ǫ> 0, Algorithm A µ ∈ n,µ,ǫ is of expected complexity Oµ(n/ǫ) (the bounding constants depend on µ). The samplers are completely described in Section 6.1 and Section 6.2.