arxiv: v1 [math.pr] 1 Jul 2017

Transcription

1 ON WEIGHTED DEPTHS IN RANDOM BINARY SEARCH TREES RAFIK AGUECH, ANIS AMRI, AND HENNING SULZBACH arxiv: v1 [math.pr] 1 Jul 2017 Abstract. Followig the model itroduced by Aguech, Lasmar ad Mahmoud [Probab. Egrg. Iform. Sci ], the weighted depth of a ode i a labelled rooted tree is the sum of all labels o the path coectig the ode to the root. We aalyze weighted depths of odes with give labels, the last iserted ode, odes ordered as visited by the depth first search process, the weighted path legth ad the weighted Wieer idex i a radom biary search tree. We establish three regimes of odes depedig o whether the secod order behaviour of their weighted depths follows from fluctuatios of the keys o the path, the depth of the odes, or both. Fially, we ivestigate a radom distributio fuctio o the uit iterval arisig as scalig limit for weighted depths of odes with at most oe child. 1. Itroductio The biary search tree is a importat data structure i computer sciece allowig for efficiet executio of database operatios such as isertio, deletio ad retrievig of data. Give a list of elemets x 1, x 2,..., x from a totally ordered set, it is the uique labelled rooted biary tree with odes costructed by successive isertio of all elemets satisfyig the followig property: for each ode i the tree with label or key, say y, all keys stored i its left right subtree are at most equal to strictly larger tha y. For a illustratio, see Figure Figure 1. biary search tree costructed from the list 4, 2, 6, 5, 7, 3, 1. Date: July 4, Mathematics Subject Classificatio. 60F05, 68P05, 68Q25. Key words ad phrases. aalysis of algorithm, data structures, biary search trees, cetral limit theorems, cotractio method, radom probability measures. 1

2 2 RAFIK AGUECH, ANIS AMRI, AND HENNING SULZBACH Properties of biary search trees are typically aalyzed uder the radom permutatio model where the data x 1,..., x are geerated by a uiformly chose permutatio of the first itegers. Amog the quatities studied i biary search trees, oe fids depths of ad distaces betwee odes related to the performace of search queries ad figer searches i the database, the total path legth measurig the cost of costructig the tree as well as the Wieer idex. Further, more complex parameters such as the height correspodig to worst case search times, the saturatio level ad the profile have bee studied thoroughly. We review literature relevat i the cotext of our work below. I this ote we complemet the wide literature o radom biary search trees by the aalysis of depths of odes, path legth ad Wieer idex i their weighted versios as itroduced by Aguech, Lasmar ad Mahmoud [1]. Here, the weighted depth of a ode is the sum of all keys stored o the path to the root. I [1], results about weighted depths of extremal paths have bee obtaied. Kuba ad Paholzer [19], [20] studied the problem i radom icreasig trees coverig the radom recursive tree ad the radom plae-orieted recursive tree. Weighted depths of odes ad the weighted height were also studied by Brouti ad Devroye [3] i a more geeral tree model, which relies o assigig weights to the edges of the tree. Further, the weighted path legth i this model was ivestigated by Rüschedorf ad Schopp [29]. Note that we deviate from the otatio itroduced i [1] ad [19] usig the term weighted depth for what is called weighted path legth there sice we also study a weighted versio of the total path legth of biary search trees. 2. Prelimiaries We itroduce some otatio. By the size of a fiite biary tree, we refer to its umber of odes. Upo embeddig a fiite rooted biary tree i the complete ifiite biary tree, a ode is called exteral if its graph distace to the biary tree is oe. Ay ode o level k 1 i a rooted biary tree is associated a vector v 1 v 2... v k {0, 1} k where v i = 0 if ad oly if the path from the root to the ode cotiues i the left subtree upo reachig level i 1. Let 1 ad 1 k. Uder the radom permutatio model short: permutatio model, let D k be the depth of the ode labelled k. By W k we deote the sum of all keys o the path from the root to the ode labelled k icludig the labels of both edpoits. For x = x 1 x 2... {0, 1}, let B x be the maximal depth amog odes of the form x 1... x k, k 0. We use X X to deote the weighted depth of the th iserted ode. Fially, we defie the height of the tree by H = sup{k N : D k > 0}. Throughout the paper, we deote by LX the distributio of a radom variable X. For real-valued X with fiite secod momet, we write σ X for its stadard deviatio. By N we deote a radom variable with the stadard ormal distributio,

3 ON WEIGHTED DEPTHS IN RANDOM BINARY SEARCH TREES 3 ad by µ the Dickma distributio o [0, characterized by its Fourier trasform, 1 e iλx e iλx 1 1 dµx = exp dx, λ R. x 0 The origis of the Dickma distributio go back to Dickma s [10] classical result o large prime divisors. Compare Hildebradt ad Teebaum [15] for a survey o the problem. I the probabilistic aalysis of algorithm, µ first arose i Hwag ad Tsai s [17] study of the complexity of Hoare s selectio algorithm. We refer to this work for a discussio of more details o the distributio, historical backgroud ad further refereces. Fially, we use the Ladau otatios little o, big O, little ω, big Ω ad big Θ as Depths ad height. We recall the followig fudametal property of radom biary search trees goig back to Devroye [6]: i probability ad with respect to all momets, we have 2 H log c, where c = is the larger of the two solutios to the trascedet equatio e = 2e x x. Next, by classical results due to Brow ad Shubert [4] ad Devroye [7], for ay x {0, 1}, i distributio, 3 B x log log N, X 2 log 2 log N. I [7, Theorem O1], the first covergece i the last display is formulated for x = 0 := The geeral case follows, sice, by symmetry, LB x = LB 0 for all x. The secod covergece was also claimed i a footote by Mahmoud ad Pittel [23]. Grübel [13] studied the process {B x : x {0, 1} }, the so-called silhouette, thereby obtaiig a fuctioal limit theorem for its itegrated versio. The asymptotic behaviour of depths of odes with give labels has bee aalyzed by Devroye ad Neiiger [9]: uiformly i 1 k ad as, 4 E [D k ] = logk k + O1, VarD k = logk k + O1. Moreover, for ay 1 k, which may deped o, i distributio 5 D k E [D k ] σ Dk N. Here, oe should also compare Grübel ad Stefaoski [14] for stroger results i the cotext of the correspodig Poisso approximatio. For a survey o depths ad distaces i biary search trees, we refer to Mahmoud s book [21]. Fially, the asymptotic behaviour of the weighted depths of the odes associated with the vectors 0 ad 1 := deoted by L ad R L ad R stad for left ad right were studied i [1]. I distributio, 6 L Y, R B 1 0, log

4 4 RAFIK AGUECH, ANIS AMRI, AND HENNING SULZBACH where Y has the Dickma distributio. The first covergece is closely related to the limit law i Theorem 3.1 i [17] Path legth ad Wieer idex. I a rooted tree, the path legth is defied as the sum over all depths of odes. Moreover, the Wieer idex is obtaied by summig all distaces of uordered pairs of vertices. For a radom biary search tree of size, we deote its path legth by P ad its Wieer idex by W. Deotig by γ the Euler-Mascheroi costat, we have 7 E [P ] = 2 log + 2γ 4 + o, VarP = goig back to Hoare [16] ad Kuth [18]. Further, by [25], 8 E [W ] = 2 2 log + 2γ o 2, VarW = 21 2π2 2 + o 2, π2 4 + o 4. 3 Cetral limit theorems for the path legth go back to Régier [27] ad Rösler [28], for the Wieer idex to Neiiger [25]. More precisely, by [25, Theorem 1.1], there exists a o-trivial radom variable Z o R 2 characterized by a stochastic fixed-poit equatio, such that, i distributio, W E[W ] 2, P E[P ] 9 Z The i.i.d. model. We also cosider biary search trees of size where the data are chose as the first values of a sequece of idepedet radom variables U 1, U 2,... each havig the uiform distributio o [0, 1]. Sice the vector raku 1,..., raku costitutes a uiformly chose permutatio, i distributio, both the permutatio model ad the i.i.d. model lead to the same ulabelled tree. We use the same otatio as i the permutatio model for quatities ot ivolvig the labels of odes, that is, X, h, H, P, W ad B x. Further, we defie the weighted path legth P as the sum of all weighted depths, ad the weighted Wieer idex W as the sum over all pairs of weighted distaces. Here, the weighted distace betwee two odes equals the sum of all labels o the path coectig them, labels of edpoits icluded. Notice that the weighted distace betwee a ode ad itself is equal to its label. Fially, aalogously to B x, we defie B x as the weighted depth of the ode of largest depth o the path x. We call {B x : x {0, 1} }, the weighted silhouette of the tree at time. 3. Mai results Our mai results are divided ito two groups: Theorems 1 ad 2 hold i the permutatio model while Theorems 3 ad 4 are formulated i the i.i.d. model Results i the permutatio model. We start with the expasios of the first two momets of the weighted depth W k. Uiformly i 1 k, as

5 ON WEIGHTED DEPTHS IN RANDOM BINARY SEARCH TREES 5, 10 E [W k ] = k logk k Ok + log, 11 VarW k = k 2 logk k Ok. It turs out that the asymptotic distributioal behavior of W k with respect to terms of secod order is etirely described by that of kd k if ad oly if k = ω/ log. Accordigly, i the remaider of this paper, we call odes with labels of order ω/ log large ad of order O/ log small. Theorem 1 Weighted depths of large odes. For k = ω/ log, 12 E [ W k kd k ] = oσ kdk. I particular, for 0 < α < 1 ad k/ α = olog 1/2, i distributio, Dk 2 log, W k 2α log 13 2 log α N, N. 2 log For the last iserted ode, i distributio, X 2 log X 14, 2 log 2 log N, ξ, where N ad ξ are idepedet ad ξ is uiformly distributed o [0, 1]. The asymptotic behavior of weighted depths of small odes is to be compared with the correspodig results i [19]. Here, aother phase trasitio occurs whe k = o/ log. Theorem 2 Weighted depths of small odes. Let k = O/ log. The, i distributio, Dk E [D k ], W k kd k 15 N, Y, σ Dk where N ad Y are idepedet ad Y has the Dickma distributio. Thus, if k log / β 0, i distributio, Dk 2 log, W k E [W k ] N, Y + 2βN 1. 2 log I particular, if k log / β = olog 1/2 with β > 0, the, i distributio, Dk 2 log, W k 2β log N, Y + 2βN. 2 log 3.2. Results i the i.i.d. model. Ay x {0, 1} correspods to a uique value x [0, 1] by x = i=0 x i2 i. This idetificatio becomes oe-to-oe upo allowig oly those x {0, 1} which cotai ifiitely may zeros ad x = 1. I the i.i.d. model, for ay x {0, 1}, k 1, the ode x 1... x k evetually appears i the sequece of biary search trees ad we write Ξ k x for its ultimate label. The followig theorem about the behavior of B x ivolves a radom cotiuous distributio fuctio arisig as the almost sure limit of Ξ k x, x [0, 1], as k. We believe that this process is of idepedet iterest ad state some of its

6 6 RAFIK AGUECH, ANIS AMRI, AND HENNING SULZBACH properties i Propositio 1 i Sectio 3.3. Figure 2 illustrate the scalig limit. The simulatios of Ξ 15 preseted i Theorem 3 Weighted silhouette. There exists a radom cotiuous ad strictly icreasig bijectio Ξx, x [0, 1], such that, almost surely, uiformly o the uit iterval, Ξ k x Ξx. For ay x [0, 1], i probability, 16 B x log Ξx. Also, for ay m 1, i probability 1 B x log Ξx 17 Further, i probability, 18 0 m dx 0. B x sup x [0,1] log c = with c as i 2. Fially, for ay x [0, 1], i distributio, B x log, B x 19 N, Ξx, log log where N ad Ξx are idepedet Figure 2. two simulatios of Ξ 15, the dotted lie beig the graph of the idetity fuctio. The ext theorem exteds the distributioal covergece result i Theorem 1.1 i [25], that is 9, by cetral limit theorems for the weighted path legth ad the weighted Wieer idex. Theorem 4 Weighted path legth ad Wieer idex. I the i.i.d. model, we have E [P ] = log + γ 3/2 + o, E [W ] = 2 log + γ 11/4 2 + o 2,

7 ON WEIGHTED DEPTHS IN RANDOM BINARY SEARCH TREES 7 ad VarP = 65 6π2 2 + o 2, 36 VarW = π2 4 + o The leadig costats i the expasios of the covariaces betwee P, W, P ad W are give i The leadig costat for CovP, W was already give i [25]. As, with covergece i distributio ad with respect to the first two momets i R 4, we have W E[W ] 2, W E[W ] 2, P E[P ], P E[P ] Z, where the limitig distributio LZ is the uique fixed-poit of the map T i 35. Coclusios. We have see that there exist three types of odes showig sigificatly differet behavior with respect to their weighted depths. By Theorem 1, for k = ω/ log, secod order fluctuatios of weighted depths are due to variatios of the depth of odes. I the secod regime, whe k = Θ/ log, variatios of weighted depths are determied by two idepedet cotributios, oe for the depths ad oe for the keys o the paths. Fially, whe k = o/ log oly fluctuatios of labels o paths ifluece secod order terms of weighted depths. The third regime ca be further subdivided with respect to the first order terms of W k ad kd k : for k = ω/ log, they coicide, for k = Θ/ log, they are of the same magitude, whereas, for k = o/ log, they are of differet scale. By Theorem 3, the weighted silhouette behaves cosiderably differet. Here, the lack of cocetratio aroud the mea leads to a iterestig radom distributio fuctio o the uit iterval as scalig limit Further results ad remarks. Model compariso. We decided to preset Theorems 3 ad 4 i the i.i.d. model rather tha i the permutatio model sice this allows for a stroger mode of covergece i 16, 17 ad a clearer presetatio of the proof of Theorem 4. I the i.i.d. model, deotig by W k the weighted depth of the ode of rak k amog the first iserted keys, Theorems 1 ad 2 remai valid upo replacig W k by W k. Similarly, Theorems 3 ad 4 hold i the permutatio model where weighted depths ad the weighted path legth are to be scaled dow by a factor ad the weighted Wieer idex by a factor 2. The covergeces i 16 ad 17 the oly hold i distributio. This ca be deduced most easily from the followig couplig of the two models: startig with the biary search tree i the i.i.d. model, also cosider the radom biary search tree i the permutatio model relyig o the permutatio raku 1,..., raku. The, for all 1 k, 20 W k W k H max 1 i U i raku i.

8 8 RAFIK AGUECH, ANIS AMRI, AND HENNING SULZBACH It is well-kow that the secod factor o right had side grows like 1/2, compare, e.g. Dosker s theorem for empirical distributio fuctios or the Dvoretzky-Kiefer- Wolfowitz iequality [11]. Combiig this, 20 ad 2 is sufficiet to trasfer all results i Sectio 3 betwee the two models. The depth first search process. I the permutatio model, let v 1,..., v +1 be the exteral odes as discovered by the depth first search process from left to right. By Dk ad W k, 1 k + 1, we deote depth ad weighted depth of the exteral ode v k. The, at the ed of Sectio 4.1, we show that, uiformly i 1 k, 21 E [ D k D k 2] = olog, E [ W k W k 2] = ovarw k. Thus, the results i Theorems 1 ad 2 also cover the secod order aalysis of the sequeces D k ad W k. Weighted distaces. I the permutatio model, let D k,l be the graph distace betwee the odes labelled 1 k l ad W k,l be the sum of all labels o the path from k to l, labels at the edpoits icluded. Asymptotic ormality for the sequece D k,l after rescalig uder the optimal coditio l k has bee obtaied i [9]. For uiformly chose odes, distributioal covergece results date back to Mahmoud ad Neiiger [22] ad Paholzer ad Prodiger [26]. Aalogously to Theorem 1, it is straightforward to prove cetral limit theorems joitly for weighted ad o-weighted distaces. We oly state the results. If l k = Ω ad k = ω/ log, the E [ W k,l kd k ld l ] = σ ldl. I particular, for 0 < s < t < 1 ad k/ s = olog 1/2, l/ t = olog 1/2, we have, i distributio, D k 2 log 2 log, D l 2 log 2 log, D k,l 4 log, W k,l 2s + t log 4 log 2 log N 1, N 2, N 1 + N 2 2, sn 1 + tn 2. Here, N 1, N 2 are idepedet radom variables both with the stadard ormal distributio. The limit process Ξ. The process Ξ i Theorem 3 is a radom distributio fuctio. I particular, it ca be regarded as a elemet i the set of càdlàg fuctios D[0, 1] cosistig of all f : [0, 1] R, such that, for all t [0, 1], ft = lim s t fs ad lim s t fs exists. The absolute value of f is defied by sup{ ft : t [0, 1]}. Edowed with Skorokhod s topology J 1, D[0, 1] becomes a Polish space. We refer to Chapter 3 i Billigsley s book [2] for detailed iformatio o this matter.

9 ON WEIGHTED DEPTHS IN RANDOM BINARY SEARCH TREES 9 Propositio 1 Properties of Ξ. The process Ξ is uique i distributio amog all càdlàg processes with fiite absolute secod momet satisfyig 22 LΞt t [0,1] = L 1 [0,1/2 tuξ2t + 1 [1/2,1 t 1 UΞ 2t 1 + U. t [0,1] Here, Ξ, Ξ, U are idepedet, U has the uiform distributio o [0, 1] ad Ξ is distributed like Ξ. We have i E [Ξt] = t for all t 0, 1; ii LΞt t [0,1] = L1 Ξ1 t t [0,1] ; iii Ξξ has the arcsie distributio with desity 1 π, x 0, 1, x1 x where Ξ, ξ are idepedet ad ξ has the uiform distributio o [0, 1]; iv for t 0, 1, LΞt has a smooth desity f t : 0, 1 0, ; v for t 0, 1/2, xf tx = f 2t x, x 0, 1, f t is strictly mootoically decreasig ad lim x 1 f t x = 0; vi with α i t := lim x 0 f i t x, i = 0, 1, t 0, 1/2 ad γ 0 = 1/4, γ 1 = 5/16, we have α i t α i t as t γ i. = 1 i for 0 < t γ i, α i t < for γ i < t < 1/2 ad Radom recursive trees. A radom recursive tree is costructed as follows: startig with the root labelled oe, i the kth step, k 2, a ode labelled k is iserted i the tree ad coected to a already existig ode chose uiformly at radom. Weighted depths i radom biary search trees differ substatially from those i radom recursive trees aalyzed i [19] where all odes show a asymptotic behaviour comparable to that of odes labelled k = o/ log i the biary search tree. The differece is highlighted by the weighted path legth. Beig of the same order as the path legth i biary search trees, it follows from results i [19] that the weighted path legth Q i a radom recursive tree of size is of order 2. The same is valid for its stadard deviatio. We cojecture that the sequece 2 Q coverges i distributio to a o-trivial limit; however, the recursive approach worked out i the proof of Theorem 4, which also applies to the aalysis of the path legth i radom recursive trees, seems ot to be fruitful i this cotext. Outlie. All results are proved i Sectio 4 startig with the proofs of Theorems 1 ad 2 as well as 21 i Sectio 4.1. Here, most argumets are based o represetatios of weighted depths as sums of bouded idepedet radom variables which go back to Devroye ad Neiiger [9]. Theorem 3 ad Propositio 1 are proved i Sectio 4.2. I this part, the costructio of the limitig process relies o suitable uiform L 1 -bouds o the icremets of the process Ξ k x x [0,1], k 1, while the properties of the limit laws formulated i Propositio 1 follow from the distributioal fixed-poit equatio 22. Fially, the proof of Theorem 4 relyig o the cotractio method is worked out i Sectio 4.3.

10 10 RAFIK AGUECH, ANIS AMRI, AND HENNING SULZBACH 4. Proofs 4.1. Weighted depths of labelled odes. I the permutatio model, let A j,k be the evet that the ode labelled k is i the subtree of the ode labelled j. The, D k = j=1 1 A j,k 1 ad W k = j=1 j1 A j,k. It is easy to see that A 1,k,..., A k 1,k ad A k+1,k,..., A,k are two families of idepedet evets; however, there exist subtle depedecies betwee the sets. Followig the approach i [9], let B j,k = A j,k 1 for j < k ad B j,k = A j,k+1 for j > k. For coveiece, let B k,k be a almost sure evet. The followig lemma summarizes results i [9] ad we refer to this paper for a proof. I this cotext, ote that Devroye [8] gives distributioal represetatios as sums of idepedet or m-idepedet idicator variables for quatities growig liearly i, such as the umber of leaves. Lemma 1. Let 1 k. The, the evets B j,k, j = 1,...,, are idepedet. For j k, we have j=1 P A j,k = 1 Bj,k \A j,k 1 k j + 1, P B j,k = 1 k j. From the lemma, it follows that E 2, ad E j=1 j1 Bj,k \A j,k 2k + log. The ideas i [9] ca also be used to aalyze secod mixed momets. Straightforward calculatios show the followig bouds: E 1 Bj,k 1 Bi,k \A i,k = O1, ad E i,j=1 ij1 Bj,k 1 Bi,k \A i,k = Ok 2 + klog 2. i,j=1 Here, both O-terms are uiform i 1 k. Defie D k = j=1 1 B j,k 1 ad W k = j=1 j1 B j,k. We make the followig observatio: O: The asymptotic statemets i 10, 11, Theorem 1 ad Theorem 2 are correct if ad oly if they are correct upo replacig D k by D k ad/or W k by W k. For i = 1, 2, 0 ad 1 k, set H i H i k. Usig Lemma 1, oe easily computes E [ Wk ] = kh 1 k, , := j=1 j i ad H i k, := Hi k 1 + Var W k = k 2 H 1 k, H2 2 k, k + 2kH1 k 1 2 H1 k 2 + k + 1. As H 1 = log+1+o1 ad H 2 = O1, both expasios 10 ad 11 follow from observatio O.

11 ON WEIGHTED DEPTHS IN RANDOM BINARY SEARCH TREES Weighted depths of large odes. We prove Theorem 1. First, 12 follows from 4 ad 23 E [ kd k W k ] k + k j P A j,k k +. j=1 For k = ω/ log, combiig 4, 5 ad 10, i distributio, Dk E [D k ], W k E [W k ] N, N. σ Dk σ Wk From here, statemet 13 follows from 4 ad 10. Cosiderig the last iserted ode with value Y, ote that, coditioally o Y = k, the correlatios betwee the evets A j,k, j < k ad A j,k, j > k vaish. More precisely, give Y = k, the family {1 Aj,k, j = 1,..., } is distributed like a family of idepedet Beroulli radom variables {V j,k : j = 1,..., } with P V j,k = 1 = k j 1 for j k ad P V k,k = 1 = 1. Thus, E [ Y X + 1 X ] 1 E k j 1 Aj,k Y = k = 1 k=1 k=1 j=1 E k j V j,k. By 3, we have X / log 2 i probability. Hece, i order to prove 14, it suffices to show that, i distributio, X 2 log, Y 24 N, ξ. 2 log For a sequece k satisfyig s k t for 0 < s < t < 1, let us coditio o the evet Y = k. The, by the cetral limit theorem for triagular arrays of row-wise idepedet uiformly bouded radom variables with divergig variace applied to V j,k, j = 1,...,, i distributio, j=1 X 2 log 2 log N. Hece, 24 follows from a applicatio of the theorem of domiated covergece otig that Y is uiformly distributed o {1,..., } Weighted depths of small odes. We prove Theorem 2. Let D> k = j=k+1 1 B ad W > j,k k = j=k+1 j1 B j,k. Sice k = O/ log, the same calculatio as i 23 shows that, 25 E [ W k W > k k D k D > k ] k 0,.

12 12 RAFIK AGUECH, ANIS AMRI, AND HENNING SULZBACH For λ, µ R, we have [ loge exp iλ D k > log / log + iµ W k > k D ] k > / = iλ λ log + log E exp i + µ j k B j,k j=k+1 log = iλ exp i λ log + µ j k 1 log + log 1 +. j k = iλ log + j=k+1 By a stadard Taylor expasio, the last display equals iλ exp i λ log + µ j k 1 log + + o1 j k j=k+1 exp 1 + iλ log j=k+1 = λ 2 / iλ log iµ j k 1 λ2 2 log j=0 j k λ2 2 log 1 exp iµ j o1 j o1 1 = λ 2 e iµx 1 /2 + dx + o1. 0 x Here, i the last step, we have used that the sum o the right had side is a Riema sum over the uit iterval whose mesh size 1 teds to zero. Thus, usig the otatio of the theorem, 1 ad Lévy s cotiuity theorem, i distributio, D> k log W k >, k D k > 26 N, Y. log I order to deduce 15 ote that, by Lemma 1, D k D > k ad D > k, W > k are idepedet while D k D > k E [ Dk D > k ] σ Dk D > k N, i distributio if ad oly if k usig the cetral limit theorem for sums of idepedet ad uiformly bouded radom variables. Sice D k E [ Dk ] = D k > E [D> k ] log log we deduce σ Dk σ Dk + D k D > k E [ Dk D > k ] σ Dk D > k Dk E [ Dk ] σ Dk, W k > k D k > N, Y, σ > Dk D σ Dk from 26 upo treatig the cases k = O1 ad k = ω1 separately. From here, the assertio 15 follows with the help of 25 ad observatio O. k,

13 ON WEIGHTED DEPTHS IN RANDOM BINARY SEARCH TREES Proof of 21. The mai observatio is that the kth exteral ode visited by the depth first search process is always cotaied i the subtree rooted at the ode labelled k. This ca be proved by iductio exploitig the decompositio of the tree at the root. Thus, deotig by H k the height of the subtree rooted at the ode labelled k, we have D k D k D k + H k, W k W k W k + M k H k. Here, M k stads for the largest label i the subtree rooted at the ode labelled k. Let T k be the size of the subtree rooted at k. The T k = 1 + T k < + T k > where T k < deotes the umber of elemets i the subtree rooted at k with values smaller tha k. By Lemma 1, for l k, we have P T k > l = P A k,k+l = 1/l + 1. Usig the same argumets for the quatity T k <, we deduce that, uiformly i 1 k, E [T k ] = Θlog, E [ T k 2] = Θ 1/2, E [ log T k 2] = O1. Thus, by a applicatio of 2, for some C 1 > 0, E [ D k D k 2] E [ H k 2] C 1 E [ log T k 2] = O1. By the same argumets, for some C 2 > 0, we have E [ W k W k 2] E [ M k H k 2] E [ k + T k 2 H k 2] From here, 21 follows from 10. C 2 k 2 + 2C 1 ke [ T k log T k 2] + C 1 E [ T k 2 log T k 2] = Ok 2 + log 2 1/ The weighted silhouette. We prove Theorem 3 ad Propositio 1. Proof of Theorem 3. We start with the uiform covergece of Ξ k. For all x [0, 1], Ξ k x Ξ k 1 x is distributed like the product of k + 1 idepedet radom variables, each of which havig the uiform distributio o [0, 1]. I particular, by the uio boud ad Markov s iequality, for ay m 1, P sup Ξ k x Ξ k 1 x t x [0,1] 2 k P k+1 i=1 U i t k 2 t m. m + 1 For k 1, let D k = {l2 k : l = 1,..., 2 k 1}. By costructio, for k 1, the map x Ξ k x is a right-cotiuous step fuctio. Further, it is cotiuous at x if ad

14 14 RAFIK AGUECH, ANIS AMRI, AND HENNING SULZBACH oly if x / D k. Next, for 0 < q < 1, [ ] E sup Ξ k x Ξ k 1 x x [0,1] = 0 q k + P sup Ξ k x Ξ k 1 x t x [0,1] k 2 t m dt q m + 1 k = q k m 1 m + 1 k q km 1. With m = 2 ad q = 2/3 the latter expressio is bouded by 2q k. By Markov s iequality, it follows that sup m sup x [0,1] Ξ m x Ξ x 0 i probability as. A applicatio of the triagle iequality shows that sup m,p sup x [0,1] Ξ m x Ξ p x 0 i probability as. By mootoicity, this covergece is almost sure. Thus, almost surely, Ξ k is uiformly Cauchy i the space of càdlàg fuctios edowed with the uiform topology. By completeess, Ξ k coverges to a limit deoted by Ξ with càdlàg paths. Moreover, Ξ is cotiuous at x / D where D = m 1 D m sice this is true for all Ξ k, k 1. For x D, let Φx be the key of the ode associated with x 1... x k 1 where k 1 is chose miimal with x D k. The, lim y x Ξx = Φx = Ξx. Thus, x Ξx is cotiuous. By the costructio of the tree, it is clear that Ξx < Ξy for ay x, y D with x < y. As D is dese i [0, 1], the process Ξ is strictly mootoically icreasig. Obviously, Ξ0 = 0 ad Ξ1 = 1; hece, Ξ is the distributio fuctio of a probability measure o [0, 1]. We tur to the covergece of B x. For ay fixed x [0, 1], display 3 implies that, as, i probability, B x/ log 1. Thus, 16 follows from the covergece Ξ k x Ξx. The covergece 16 is with respect to all momets sice B x H ad we have covergece of all momets i 2. By the theorem of domiated covergece, for ay m 1, agai usig 2, we have 1 [ B x E log Ξx m] dx 0. 0 This shows 17. To prove 18, ote that, for ay k 1, sup x [0,1] B x is larger tha the product of the height of the subtree rooted at the ode w k := o level k ad Ξ k 1 1. Let ε > 0. Fix k large eough such that P Ξ k 1 1 < 1 ε < ε. Coditioal o its size, the subtree rooted at w k is a radom biary search tree. Sice its size grows liearly i as, it follows from 2 that, for all sufficietly large, its height exceeds c ε log with probability at least 1 ε. For these values of, we have sup x [0,1] B x c 6ε log with probability at least 1 2ε. As ε was chose arbitrarily, this shows 18. For the joit covergece of B x ad B x for fixed x [0, 1], we abbreviate B := B x, B := B x, Ξ k := Ξ k x, Ξ = Ξx ad B = B log / log. Note that Ξ ad B are ot idepedet which causes the proof to be more techical. Deote by N k the time whe the ode associated with x 1... x k is iserted i the biary search tree. For ay ε > 0, we ca choose k, L 1 such that, for all dt

15 ON WEIGHTED DEPTHS IN RANDOM BINARY SEARCH TREES 15 sufficietly large, B P Ξ k Ξ ε + P N k L + P log Ξ ε ε. Further, there exists δ > 0 such that P Ξ k Ξ k 1 δ ε. The, for r, y R with P Ξ = y = 0, ad large eough, B P B r, log y 2ε + P B r, Ξ k y + 2ε, Ξ k Ξ k 1 δ, N k < L. Let x = x k+1 x k+2..., V 1, V 2,... be a idepedet copy of U 1, U 2,... ad Bi, p := 1 {Vi p}, 0, p [0, 1]. i=1 Give Ξ k, Ξ k Ξ k 1, N k, o N k <, B is distributed like B Bi N k, Ξ k Ξ k 1 x + k/ log where B x is distributed like B x ad idepedet from the remaiig quatities. We deduce B P B r, log y k 2ε + P + B log Bi L,δ x r, Ξ k y + 2ε, Ξ k Ξ k 1 δ, N k < L k 3ε + P + B log Bi L,δ x r P Ξ y + 2ε Usig the asymptotic ormality of B x after rescalig i 3, takig the limit superior as ad the lettig ε ted to zero, we obtai B lim sup P B r, log y P N r P Ξ y. The proof of the coverse directio establishig 19 is easier. It rus alog the same lies upo usig the trivial bouds Ξ k Ξ k 1 1 ad N k 0. Proof of Propositio 1. We start with the characterizatio of the distributio of the process. For a determiistic sequece of pairwise differet umbers u 1, u 2,... o the uit iterval, we defie ξ k x aalogously to Ξ k x i the ifiite biary search tree costructed from this sequece. Here, we abbreviate ξ k x = 0 if the ode x 1... x k is ot i the tree. Let m, m 1, be the subsequece defied by the elemets u m < u 1 ad u + m, m 1, be the subsequece defied by the elemets u + m > u 1. At least oe of these sequeces is ifiite. For m 1, let ym = u m /u 1 ad y m + = u + m u 1 /1 u 1. Next, defie ξ k ξ+ k, respectively aalogously to ξ k based o the sequece ym y m, + respectively. By costructio, for k 1, ξ k x = 1 [0,1/2 xu 1 ξ k 1 2x + 1 [1/2,1]x1 u 1 ξ + k 1 2x 1 + u 1. Applyig the costructio to the sequece U 1, U 2,... yields Ξ k x = 1 [0,1/2 xu 1 Ξ k 1 2x + 1 [1/2,1]x1 U 1 Ξ + k 1 2x 1 + U 1. Almost surely, the radom sequeces y m ad y + m are both ifiite ad Ξ k, Ξ+ k are idepedet copies of Ξ k. Further, both sequeces are idepedet of U 1.

16 16 RAFIK AGUECH, ANIS AMRI, AND HENNING SULZBACH Hece, lettig k i the last display, we obtai 22 o a almost sure level. The characterizatio of LΞ by 22 follows from a stadard cotractio argumet, the argumet o page 267 i [12] applies to our settig without ay modificatios. We move o to the statemets i vi o the margial distributios of the process. Here, we use otatio that was itroduced i the proof of Theorem 3. By cotiuity, it suffices to show i for x D. Let k 1. By symmetry, for 1 i 2 k 1, we have E [ Φi2 k ] = i2 k. Thus, the assertio follows for x D sice Φx = Ξx. The symmetry statemet ii is remiiscet of the fact that the uiform distributio o [0, 1] is symmetric aroud 1/2. More precisely, we apply the reflectio argumet from [1] which is at the core of the proof of the secod assertio i 6. Let U 1 = 1 U 1, U 2 = 1 U 2,... ad defie Ξ aalogously to Ξ i the biary search tree process relyig o the sequece U 1, U 2,... The, Ξ t + Ξ1 t = 1 for all t [0, 1] which proves ii. With Y = Ξξ, 22 yields LY = LUY + 1 A 1 U, where 1 A, U, Y are idepedet ad P A = 1/2. From [5], it follows that Y has the arcsie distributio, provig iii. We move o to the statemets about the distributio of Ξt. Let t 0, 1/2. Sice Ξ is strictly icreasig, we have Ξ2t 0, 1 almost surely. By 22, LΞt = LUΞ2t with coditios as i 22. Therefore, LΞt admits a desity. By symmetry, the same is true for t 1/2, 1. For t 0, 1/2, by coditioig o the value of U, oe fids the desity 27 [ ] 1[x,1] Ξ2t f t x = E, Ξ2t x 0, 1]. f t x is mootoically decreasig ad cotiuous o 0, 1] with f1 = 0. For t 1/2, 1, f t x = f 1 t 1 x, x 0, 1 is a desity of LΞt by ii. By 27, for t 0, 1/2, x 0, 1, 28 f t x = 1 x f 2t y dy, or xf y tx = f 2t x. Upo settig f 0 = f 1 = 0, the last idetity also holds for t = 0 ad t = 1/2 sice f 1/2 = 1 [0,1] is a desity of LΞ1/2. Thus, for ay t 0, 1, f t is smooth o 0, 1. Sice the uiform distributio takes values arbitrarily close to oe, it follows that, for all δ > 0, t 0, 1, we have P Ξt > 1 δ > 0. Hece, for all t 0, 1, the desity f t is strictly positive o 0, 1. Thus, for t 0, 1/2, f t is strictly mootoically decreasig. Summarizig, we have show iv ad v. For t 0, 1/4], the assertio α 0 t = i vi follows immediately from 28 sice α 0 2t > 0. Let 1/4 < t < 1/2. Assume α t <. The, f 21 2t1 <. By 28, it follows that f 1 2t1 is fiite ad hece f 2t0 is fiite. Thus, f 2t y/y is bouded i a eighhourhood of zero ad α 0 t <. For t > 3/8, we have t > 1/2, thus, α 0 t <. Iteratig this argumet leads to α 0 t < for all 1/3 < t < 1/2. I order to proceed further, ote that, for t > 1/4, there exists k N, such that, i probability, Ξt Z := U 1 U U 2 k l=1 U 2+l.

17 ON WEIGHTED DEPTHS IN RANDOM BINARY SEARCH TREES 17 Z admits a desity f Z give by Thus, f Z x = x rydy xrx, lim f Z x = 1 + x 0 rx = 1 x k x 2 P rydy <. l=1 U 2+l x v dv. 1 v It follows that α 0 t rxdx <. Sice Ξ is icreasig, the fuctio 0 t α 0 t is decreasig. Thus, by mootoicity ad cotiuity, it follows α 0 t as t 1/4. For t 1/4, α 0 t = follows immediately from 28 sice α 0 For 1/4 < t < 1/2, the remaiig statemets about α 1 t the results for α 0 t sice α 1 t = α t 2t <. are direct corollaries of. This fiishes the proof of vi. The curvature. We make a cocludig remark about the curvature of f t, t 0, 1/2. First, sice xf t x = f 2tx f tx, for 0 < t 1/4, the fuctio f t is covex. From 28 it is easy to deduce f 1/3 x = 21 x. Sice f 1/3 = f 1/2 = 0, it is plausible to cojecture that f t is covex for t 1/3 ad cocave for 1/3 t < 1/2. Cocavity at ratioal poits with small deomiator such as t = 3/8 or t = 5/12 ca be verified by had usig Weighted path legth ad Wieer idex. I order to obtai mea ad variace for the weighted path legth ad the weighted Wieer idex, we use the reflectio argumet from the proof of Propositio 1 ii. To this ed, let P ad W deote weighted path legth ad weighted Wieer idex i the biary search tree built from the sequece U1 = 1 U 1, U2 = 1 U 2,... The, P + P = P + ad W + W = W + 2 providig the claimed expasios for E [P ] ad E [W ] upo recallig 7 ad 8. For a fiite rooted labelled biary tree T, deote by pt its path legth, by pt its weighted path legth, by wt its Wieer idex ad by wt its weighted Wieer idex. Let T 1, T 2 be its left ad right subtree ad x the label of the root. The, deotig by T the size of T, for T 1, 29 pt = pt 1 + pt 2 + T 1, 30 wt = wt 1 + wt 2 + T 2 + 1pT 1 + T 1 + 1pT 2 + T + 2 T 1 T 2 1. The first statemet is obvious, the argumet for the secod ca be foud i [25]. For the weighted quatities, oe obtais 31 pt = pt 1 + pt 2 + T x, 32 wt = wt 1 + wt 2 + T 2 + 1pT 1 + T 1 + 1pT 2 + T + T 1 T 2 x. Agai, the first assertio is easy to see ad we oly justify the secod. The terms wt 1 ad wt 2 accout for weighted distaces withi the subtrees. The sum

18 18 RAFIK AGUECH, ANIS AMRI, AND HENNING SULZBACH of all weighted distaces betwee odes i the left subtree ad the root equals pt 1 + T 1 x. Replacig T 1 by T 2, we obtai the aalogous sum i the right subtree. The sum of all distaces betwee odes i differet subtrees equals T 1 pt 2 + T 2 pt 1 + T 1 T 2 x. Fially, we eed to add x for the weighted distace of the root to itself. Addig up the terms ad simplifyig leads to 32. For α, β > 0 let αt + β be the tree obtaied from T where each label y is replaced by αy + β. Obviously, pt = pαt +β with the aalogous idetity for the Wieer idex. For the weighted quatities, we have pαt + β = αpt + pt + T β, wαt + β = αwt + wt + T T + 1/2β. Let T be the biary search tree of size i the i.i.d. model. The, give I := raku 1, U := U 1, i distributio, the trees 1 U T 1 ad 1 1 U T 2 U 1 U are idepedet biary search trees of size I 1 ad I, costructed from idepedet sequeces of uiformly distributed radom variables o [0, 1]. Thus, combiig 29 34, for the vector Y = W, W, P, P T, we have Y d = + + U I U I 0 0 U U U I 1 U I U I U U I 3I + 2U/ I 1 I 2 I U 1, Y I 1 Y I where Y, Y, I, U are idepedet ad Y is distributed like Y. Here, d = idicates that left- ad righthad side are idetically distributed. We cosider the sequece Z 0 defied by Z := W E[W ] 2, W E[W ] 2, P E[P ], P T E[P ], 1,

19 ON WEIGHTED DEPTHS IN RANDOM BINARY SEARCH TREES 19 ad Z 0 = 0. Let α = E [W ], β = E [W ], γ = E [P ] ad δ = E [P ]. Further, let A 1 = A 2 = I 1 2 U 0 1 I 1 I 1 U 0 0 I I 1 I 1 I U 0, I I 2 1 U 1 I 2 U I 1 I 1 U I 1 I U 0 1 I 2 I 0 1 I I 1 U 1 I U, I ad C = C 1, C 2, C 3, C 4 T with C 1 = U 2 α I U 2 α I + U 2 β I + U + 1 I U I 2 γ I + U I C 2 = 1 2 β I β I I β, C 3 = U γ I U γ I + U δ I + C 4 = 1 δ I δ I δ. The, from the recurrece for Y, it follows γ I 1 2 δ I + U 2 + I 3I α, 1 I 1 1 δ I 1 + I 2 δ I 2 I U 1 γ, Z d = A 1 Z I 1 + A 2 Z I + C, 1, where Z, Z, I, U are idepedet ad Z is distributed like Z. We prove covergece of Z i distributio by a applicatio of the cotractio method. To this ed, ote that I / U almost surely by the strog law of large umbers. Thus, with covergece i L 2 ad almost surely, A 1 A 1 := A 2 A 2 := U 3 0 U 2 1 U 0 0 U 2 0 U1 U 0 0 U U, 1 U 3 U1 U 2 U1 U 2 U 2 1 U 0 1 U 2 0 U1 U U 2 U1 U U,

20 20 RAFIK AGUECH, ANIS AMRI, AND HENNING SULZBACH ad C C := U 2 log U + 1 U 2 log 1 U + U 14U 2 + 9U + 5/4 2U log U + 21 U log1 U + 6U1 U U 2 l U + 1 U 2 l1 U + U 2U l U + 21 U l 1 U + 1 For a quadratic matrix A, deote by A op its spectral radius. By calculatig the eigevalues of A 1 A T 1 ad A 2 A T 2, oe checks that A 1 op = U ad A 2 op = 1 U. Thus, E [ A 1 A T 1 op ] + E [ A2 A T 2 op ] E [ A1 2 op] + E [ A2 2 op] < 1. Moreover, we have P I {1,..., l} {} 0 for all fixed l. Thus, by Theorem 4.1 i [24], i distributio ad with covergece of the first two momets, we have Z W, W, P, P where LW, W, P, P is the uique fixed-poit of the map: 35 T : M 4 20 M 4 20, T µ = L A 1 Z + A 2 Z + C, with A 1, A 2, C defied above, where Z, Z, U are idepedet ad LZ = LZ = µ. Here, M 4 20 deotes the set of probability measures o R 4 with fiite absolute secod momet ad zero mea. Variaces ad covariaces ca be computed successively usig the fixed-poit equatio, e.g. i the followig order: E [ P 2], E [P W ], E [ W 2], E [P P], E [ P 2], E [PW ], E [P W], E [W W], E [PW], E [ W 2]. Additioally to the variaces give i the theorem, oe obtais π2 CovP, P 2, 6 10 π2 CovP, W 3, 3 10 π2 CovW, W 4, π2 CovP, W 3, 3 10 π2 CovP, W 3, π2 CovP, W Ackowledgemets The first author is grateful to the Kig Saud Uiversity, Deaship of Scietific Research, College of Sciece Research Ceter. The research of the third author was supported by a Feodor Lye Fellowship of the Alexader vo Humboldt- Foudatio. Refereces [1] Rafik Aguech, Nabil Lasmar, ad Hosam Mahmoud, Extremal weighted path legths i radom biary search trees, Probab. Egrg. Iform. Sci., 211, [2] Patrick Billigsley, Covergece of probability measures, Wiley Series i Probability ad Statistics: Probability ad Statistics. Joh Wiley & Sos, Ic., New York 1999 [3] Nicolas Brouti ad Luc Devroye, Large deviatios for the weighted height of a exteded class of trees, Algorithmica, 463-4, [4] Gerald G. Brow ad Bruo O. Shubert, O radom biary trees, Math. Oper. Res., 91, [5] Robert Che, Elsa Li ad Ala Zame, Aother arc sie law, Sakhyā Ser. A, 433,

21 ON WEIGHTED DEPTHS IN RANDOM BINARY SEARCH TREES 21 [6] Luc Devroye, A ote o the height of biary search trees, J. Assoc. Comput. Mach., 333, [7] Luc Devroye, Applicatios of the theory of records i the study of radom trees, Acta Iform., 261-2, [8] Luc Devroye, Limit laws for local couters i radom biary search trees, Radom Structures Algorithms, 23, [9] Luc Devroye ad Ralph Neiiger, Distaces ad figer search i radom biary search trees, SIAM J. Comput., 333, [10] Karl Dickma, O the frequecy of umbers cotaiig prime factors of a certai relative magitude, Arkiv för Mathematik, Astroomi och Fysik, 22A10, [11] Aryeh Dvoretzky, Jack Kiefer, ad Jacob Wolfowitz, Asymptotic miimax character of the sample distributio fuctio ad of the classical multiomial estimator, A. Math. Statist., 27, [12] Rudolf Grübel ad Uwe Rösler, Asymptotic distributio theory for Hoare s selectio algorithm, Adv. i Appl. Probab., 281, [13] Rudolf Grübel, O the silhouette of biary search trees, A. Appl. Probab., 195, [14] Rudolf Grübel ad Nikolče Stefaoski, Mixed Poisso approximatio of ode depth distributios i radom biary search trees, A. Appl. Probab., 151A, [15] Adolf Hildebrad ad Gérald Teebaum, Itegers without large prime factors, J. Théor. Nombres Bordeaux, 52, [16] C. A. R. Hoare, Quicksort, Comput. J., 5, [17] Hsie-Kuei Hwag ad Tsug-Hsi Tsai, Quickselect ad the Dickma fuctio, Combi. Probab. Comput., 114, [18] Doald. E. Kuth, The Art of Computer Programmig: Sortig ad Searchig, volume 3, Addiso-Wesley, Readig, MA 1973 [19] Markus Kuba ad Alois Paholzer, O weighted path legths ad distaces i icreasig trees, Probab. Egrg. Iform. Sci., 213, [20] Markus Kuba ad Alois Paholzer, O edge-weighted recursive trees ad iversios i radom permutatios, Discrete Math., 3084, [21] Hosam M. Mahmoud, Evolutio of radom search trees, Wiley-Itersciece Series i Discrete Mathematics ad Optimizatio. Joh Wiley & Sos, Ic., New York 1992 [22] Hosam M. Mahmoud ad Ralph Neiiger, Distributio of distaces i radom biary search trees, A. Appl. Probab., 131, [23] Hosam M. Mahmoud ad Boris Pittel, O the most probable shape of a search tree grow from a radom permutatio, SIAM J. Algebraic Discrete Methods, 51, [24] Ralph Neiiger, O a multivariate cotractio method for radom recursive structures with applicatios to Quicksort, Radom Structures Algorithms, 193-4, Aalysis of algorithms Kryica Morska, [25] Ralph Neiiger, The Wieer idex of radom trees, Combi. Probab. Comput., 11 6, [26] Alois Paholzer ad Helmut Prodiger, Spaig tree size i radom biary search trees. A. Appl. Probab., 142, [27] Mireille Régier, A limitig distributio for quicksort, RAIRO Iform. Théor. Appl., 233, [28] Uwe Rösler, A limit theorem for Quicksort, RAIRO Iform. Théor. Appl., 25, [29] Ludger Rüschedorf ad Eva-Maria Schopp, Note o the weighted iteral path legth of b-ary trees., Discrete Math. Theor. Comput. Sci., 91,

22 22 RAFIK AGUECH, ANIS AMRI, AND HENNING SULZBACH Rafik Aguech, Departmet of Statistic ad Operatio Research, College of Scieces, Kig Saud Uiversity, P.O. Box 2455, Riyadh 11451, Saudi Arabia address: Ais Amri, Uiversity of Moastir, Aveue Taher Hadded B.P 56, Moastir 5000, Tuisia address: Heig Sulzbach, McGill Uiversity, 3480 Uiversity Street, H3A 0E9 Motreal, QC, Caada address: Curret address: School of Mathematics, Uiversity of Birmigham, Birmigham B15 2TT, Great Britai