arxiv: v1 [cs.dm] 7 Apr 2011

Transcription

1 Explosive Percolatio i Erdős-Réyi-Like Radom Graph Processes Kostatios Paagiotou Reto Spöhel Max Plack Istitute for Iformatics Saarbrücke, Germay {kpaagio rspoehel}@mpi-if.mpg.de arxiv: v1 [cs.dm] 7 Apr 2011 Agelika Steger Heig Thomas Istitute of Theoretical omputer Sciece ETH Zurich, 8092 Zurich, Switzerlad {steger hthomas}@if.ethz.ch Abstract. The evolutio of the largest compoet has bee studied itesely i a variety of radom graph processes, startig i 1960 with the Erdős-Réyi process. It is well kow that this process udergoes a phase trasitio at /2 edges whe, asymptotically almost surely, a liear-sized compoet appears. Moreover, this phase trasitio is cotiuous, i.e., i the limit the fuctio f(c) deotig the fractio of vertices i the largest compoet i the process after c edge isertios is cotiuous. A variatio of the Erdős-Réyi process are the so-called Achlioptas processes i which i every step a radom pair of edges is draw, ad a fixed edge-selectio rule selects oe of them to be icluded i the graph while the other is put back. Recetly, Achlioptas, D Souza ad Specer [1] gave strog umerical evidece that a variety of edge-selectio rules exhibit a discotiuous phase trasitio. However, Riorda ad Warke [11] very recetly showed that all Achlioptas processes have a cotiuous phase trasitio. I this work we prove discotiuous phase trasitios for a class of Erdős-Réyi-like processes i which i every step we coect two vertices, oe chose radomly from all vertices, ad oe chose radomly from a restricted set of vertices. 1. Itroductio I their semial paper [8] from 1960 Erdős ad Réyi aalyze the size of the largest compoet i the promiet radom graph model G,m, a graph draw uiformly at radom from all graphs o vertices with m edges. For ay graph G let L 1 (G) deote the size of its largest compoet. We say that a evet occurs asymptotically almost surely (a.a.s.) if it occurs with probability 1 o(1) as teds to ifiity. Theorem 1 ([8]). For ay costat c > 0 the followig holds. If c < 0.5, the a.a.s. L 1 (G, c ) = O(log ). If c > 0.5, the a.a.s. L 1 (G, c ) = Ω(), ad all other compoets have O(log ) vertices. This result ca be see from a radom graph process perspective. Startig with the empty graph o the vertex set [] := {1, 2,..., }, we add a sigle edge chose uiformly at radom from all o-edges i every step. It is ot hard to see that the graph after isertig m edges is distributed as The author was supported by a fellowship of the Swiss Natioal Sciece Foudatio The author was supported by the Swiss Natioal Sciece Foudatio, grat

2 2 2E E+07 ER MP MS 2E E+07 ER HR (0.25) HR (0.5) HR (0.9) 1E+07 1E+07 5E+06 5E+06 0E+00 0E+00 5E+06 1E E+07 2E+07 0E+00 0E+00 5E+06 1E E+07 2E+07 Figure 1. Evolutio of the largest compoet over the first edges of Erdős- Réyi process, mi-product ad mi-sum rule ad the half-restricted process with parameters 0.25, 0.5 ad 0.9 o 20 millio vertices G,m. I this cotext Theorem 1 states that, asymptotically, a liear-sized compoet (a so called giat ) appears aroud the time whe we have iserted /2 edges. We say that the Erdős-Réyi process has a phase trasitio at /2. This phase trasitio has bee studied i great detail (for a survey see hapter 3 i [9]). It is kow that L 1 (G,c ) a.a.s. satisfies L 1 (G,c ) = (1 + o(1))f(c) for some cotiuous fuctio f(c) with f(c) = 0 for every c < 0.5 ad lim c f(c) = 0. Thus, the phase trasitio i the Erdős-Réyi process is cotiuous 1 (see Figure 1). A variat of the Erdős-Réyi process which gaied much attetio over the last decade, mostly cocerig the questio whether oe ca delay or accelerate the appearace of the giat compoet [2, 12], is give by the so-called Achlioptas processes. I this variat oe starts with the empty graph o the vertex set [], ad i every step gets preseted a pair of edges chose uiformly at radom from all pairs of o-edges i the curret graph. A fixed edge-selectio rule the selects exactly oe of them to be iserted ito the graph while the other is put back ito the pool of o-edges. Recetly, Achlioptas, D Souza ad Specer [1] provided strog umerical evidece that the miproduct (select the edge that miimizes the product of the compoet sizes of the edpoits) ad mi-sum rule (select the edge that miimizes the respective sum) exhibit discotiuous trasitios, i cotrast to a variety of closely related edge-selectio rules, i particular the oes aalyzed i [12]. The authors cojectured that a.a.s. the umber of edge isertios betwee the appearace of a compoet of size ad oe of size /2 is at most 2 2/3 = o(), that is, at the phase trasitio a costat fractio of the vertices is accumulated ito a sigle giat compoet withi a subliear umber of steps (see Figure 1). This pheomeo is also called explosive percolatio ad of great iterest for may physicists. Thus, a series of papers has bee devoted to uderstadig this phase trasitio (see e.g. [3, 4, 6, 10]), most of the argumets ot beig rigorous but supported by computer simulatios. outerig the umerical evidece it was claimed i [5] that the trasitio is actually cotiuous. Recetly, Riorda ad Warke i [11] ideed cofirmed this claim with a rigorous proof. I fact, their argumet shows cotiuous phase trasitios for a eve larger class of processes. A radom graph process (G(T )) T 0 is called a l-vertex rule if G(0) is empty, ad G(T ) is obtaied from G(T 1) by drawig a set V T of l vertices uiformly at radom ad addig oly edges withi V T, 1 I the literature such a phase trasitio is also called secod order while discotiuous phase trasitios are called first order.

3 3 where at least oe edge has to be added if all vertices i V T are i differet compoets. Observe that every Achlioptas process is a 4-vertex rule. 2. Our Results I this work we prove that a variat of the Erdős-Réyi process which we call half-restricted process exhibits a discotiuous phase trasitio. Ituitively speakig, a discotiuous trasitio ca oly occur if oe avoids coectig two compoets that are already large. The idea to avoid this i a Achlioptas process is to preset two edges ad choose the oe that coects smaller compoets, which essetially is what the mi-product ad mi-sum rule do. However, we pursue a differet approach here. We coect two vertices i every step, but we restrict oe of them to be withi the smaller compoets. For every 0 < β 1 ad every labeled graph G with := V (G) we defie the restricted vertex set R β (G) to be the β vertices i smallest compoets. Precisely, let v 1, v 2,..., v be the vertices of G sorted ascedig i the sizes of the compoets they are cotaied i, where vertices with the same compoet size are sorted lexicographically. The R β (G) := {v 1, v 2,..., v β }. The half-restricted process has a parameter 0 < β 1 ad starts with the empty graph R β (0) = R,β (0) o the vertex set []. I every step T 1 we draw oe urestricted vertex v 1 [] uiformly at radom ad, idepedetly, oe restricted vertex v 2 R β (R β (T 1)) uiformly at radom. We obtai R β (T ) by isertig a edge betwee v 1 ad v 2 if the edge is ot already preset (i which case we do othig). Note that for β < 1 the half-restricted process is ot a l-vertex rule. For a half-restricted process let α T be the radom variable that deotes the maximum size of all compoets that the restricted vertex ca be draw from i step T. learly, α T is icreasig i T. For every positive iteger k we deote the radom variable for the last step whe α T is still below k by T k := max{t : α T < k}. Our mai result is that for ay parameter d < 1 the half-restricted process exhibits explosive percolatio. Thus, eve though Figure 1 suggests that the phase trasitios of the mi-product (or mi-sum) rule ad the half-restricted process behave similarly, their mathematical structure is fudametally differet. Theorem 2 (Mai Result). Let 0 < β < 1. For every K = K() with l() 1.02 K, every = () with 1 l(k), ad every ε > 0 we have a.a.s. that (i) L 1 (R β (T )) K, ad (ii) L 1 (R β (T + / l())) (1 ε)(1 β). Note that for K = l() 2 ad = l l l the theorem shows that, a.a.s., the umber of steps from the first appearace of a compoet of size l() 2 to the appearace of a compoet of size ((1 β)/2) is O(/ l l l l ) = o(). 3. Proof of the Mai Result I our proof we eed a rather techical lemma itroduced i the followig. Let N be a positive iteger, ad for every 0 i < N let X i Geom( N i N ) be a geometrically distributed radom variable ad set X(a, b) := b i=a X i for every 0 a < b < N. All subsequet statemets ad argumets are about sums of geometrically distributed radom variables, but we will use that they have the followig combiatorial iterpretatio to a coupo collector sceario with N coupos. (I a coupo collector sceario, we have a umber of differet coupos ad repeatedly

4 4 draw oe uiformly at radom with replacemet. We are iterested i how ofte we have to draw util we have see every coupo.) Observe that X i is distributed like the umber of coupos we eed to draw while holdig exactly i differet coupos (waitig for the (i + 1)st), ad thus X(a, b) ca be see as the umber of coupos we draw while holdig betwee a ad b differet coupos or, equivaletly, to obtai b a + 1 coupos of a fixed subset of N a coupos. Note that E[X(a, b)] = b i=a N a N N i = N i=n b 1 i = N(H N a H N b 1 ), (1) where H = i=1 1/i deotes the -th harmoic umber for every 1. The followig lemma cocers the situatio whe a = N ω(1) ad b = N 2. I terms of coupo collector, we wat to collect all but 1 coupo of a fixed subset of ω(1) coupos. We show a expoetial boud o the probability that the umber of coupos we eed to draw for this is much lower tha its expectatio. Lemma 3. Let k = k(n) = ω(1) ad s = s(n) N l(k). The for N large eough Proof. First ote that by (1) we have Pr[X(N k, N 2) s] e k0.99. E[X(N k, N 2)] = N(H k H 1 ) k=ω(1) = (1 + o(1))n l(k). Thus, s E[X(N k, N 2)]. osider N differet coupos ad let {c 1, c 2,..., c k } be a fixed subset of them. We ow draw s coupos with replacemet ad call the resultig set F. For every 1 j k let Y j be a idicator variable for the evet that coupo c j is ot draw withi these s trials, i.e., c j / F. By the commets precedig this lemma, the probability of the evet X(N k, N 2) s equals the probability that at most 1 of these coupos is ot draw withi the s trials. Hece, for Y := k j=1 Y j we observe that Pr[X(N k, N 2) s] = Pr[Y 1]. (2) Usig the idetity 1 x e 2x which holds for all 0 x 1/2 we have for every 1 j k that ( E[Y j ] = Pr[Y j = 1] = 1 1 ) s e 2s N, N ad thus, for N large eough, E[Y ] ke 2s N = k 1 s N l(k) s N l(k) 8k (3) Oe ca check that the radom variables Y 1, Y 2,..., Y k are egatively associated (see e.g. hapter 3 i [7]). Hece, we ca apply heroff bouds to Y ad obtai for N large eough that [ Pr[Y 1] (3) ( Pr Y 1 1 ) ] (3) E[Y ]/8 E[Y ] e e k0.99, 2 which together with (2) fiishes the proof. We ow tur to the proof of Theorem 2. Proof (of Theorem 2). We fix β < 1, K = K() with l() 1.02 K, = () with 1 l(k), ad ε > 0. To simplify otatio we write R(T ) istead of R β (T ). We first address (i). We eed to show that at the step whe the restricted vertex ca for the first time be i a compoet of size, there is a.a.s. o compoet of size larger tha K. The mai idea is that a large compoet eeds to be draw by the urestricted vertex so ofte that this is ulikely to happe withi T steps.

5 5 Note that up to step T two compoets of size at least ca ever be merged by a edge sice the restricted vertex i every step is draw from vertices i compoets of size less tha. Hece, we ca easily keep track of these compoets. We call a compoet A i R(T ) a chuk if it has size at least. Let A 1, A 2,... deote all chuks i order of appearace durig the process. By the pigeohole priciple, there ca be at most / chuks. For every 1 i / we deote by E i the evet that chuk A i has size larger tha K i R(T ). We will show that Pr[E i ] 1/ for every 1 i. By applyig the uio boud this implies / E i Pr[R(T ) cotais a comp. of size > K] = Pr i=1 1 1 = o(1). (4) It remais to boud Pr[E i ] for every 1 i /. Let 1 i / be fixed for the remaider of the proof. Sice for every T T the restricted vertex is draw from vertices i compoets of size smaller tha, the chuk A i ca grow by at most i every step. Hece, the chuk has size at most (j + 1) before a vertex from the chuk is draw for the jth time. Sice moreover, a chuk has size at most 2 whe it appears, a vertex from the chuk eeds to be draw i at least K/ 1 steps after its appearace for E i to happe. Let X 2, X 4,..., X K/ deote the umber of steps betwee steps i which we draw a vertex from A i. That is, X 2 is the umber of steps from the appearace of the chuk util a vertex from A i is draw for the first time, X 3 is the time from that step util a vertex from A i is draw for the secod time, ad so o. Furthermore, let X := K/ j=2 X j. The, Pr[E i ] Pr[X T ], (5) ad it thus suffices to show that Pr[X T ] 1/. learly, oly the urestricted vertex i every step ca be i A i, ad thus, X j is stochastically domiated by a geometrically distributed radom variable with parameter j/. Hece, let Y j be a radom variable with Y j Geom(j/) ad set Y := K/ j=2 Y j. The Pr[X T ] Pr[Y T ]. (6) This is exactly the setup of Lemma 3 with N = /, k = K/ ad s = T. ocerig the prerequisites of the lemma, we have k = ω(1) sice l(k) K. Furthermore, it is ot hard to show that T 4 (see e.g. Lemma 3 ad the remark followig the proof i [11]) ad thus N l(k) = ( K ) l l(k) = (1 o(1)) l(k) l(k) 4 T = s. Hece, we ca apply Lemma 3, which gives us for large eough that Sice K l() 1.02 ad l(k) we have Pr[Y T ] e (K/)0.99. Pr[Y T ] e l() = 1. Usig (6), (5) ad (4) this settles the proof of (i), ad it remais to prove (ii), i.e., we have to show that R(T + / l()) cotais a compoet of size (1 ε)(1 β) with high probability. We set a := /(2 l()) ad split the proof ito two parts. I the first part (the first a additioal steps after T ) we collect a suitable amout of vertices i compoets of size at least, ad i the secod part (the remaiig a steps) we actually build a giat compoet o these vertices. osider the steps T + 1 to T + a i the graph process. Let U(T ) deote set of vertices i compoets of size at least i R(T ). We show a lower boud o the size of U(T + a) that holds with high probability. Note that by defiitio of T we have U(T + 1) (1 β). If U(T + a) (1 β/2), we have a boud that is by far good eough for our purposes, hece we restrict ourselves to the case U(T + i) (1 β/2) for every 1 i a.

6 6 For every 1 i a let X i := U(T + i) U(T + i 1) deote the umber of vertices added to the compoets of size at least i the ith additioal step. Furthermore, let X := a i=1 X i. We lower boud the probability that X i cotributes at least 1 vertex. learly, this happes if the urestricted vertex is draw from compoets of size at least, which by defiitio of T happes with probability at least 1 β, ad if the restricted vertex is draw from compoets of size smaller tha, which happes with probability at least 1/2 sice we assume U(T + i) (1 β/2). Hece, for every 1 i a we have Pr[X i 1] 1 β 2 Hece, E[X] a(1 β)/2 = (1 β)/(4 l()). I particular, X stochastically domiates a sum of a idepedet Beroulli radom variables with parameter (1 β)/2, such that we ca use heroff bouds to obtai [ Pr X < 1 ] (1 β) e Θ( l() ) = o(1). 2 4 l() Hece, we ca i the followig coditio o the evet that. U(T + a) (1 β) + 1 β 8 l(). (7) We ow look at the secod half of additioal steps, i.e., steps T + a + 1 to T + 2a, ad show that a.a.s. i these steps a sufficietly large compoet is created withi U := U(T + a). Note that U is a fixed set of vertices which does ot chage from step to step. Assume that R(T + 2a) has o compoet of size (1 ε)(1 β). We ow show that a.a.s. this assumptio leads to a cotradictio. We call a step successful if it coects two compoets i U. Sice every compoet i U has size at least we have that (1 β) successful steps will coect all compoets i U such that U forms oe giat compoet of size at least (1 β). We ow compute the probability to have a successful step if we do ot have a compoet of size (1 ε)(1 β). For a successful step, the restricted vertex eeds to be i U which happes with probability at least U (1 β) β (7) 1 β 8β l(), ad the urestricted vertex eeds to be draw from a differet compoet i U, which happes with probability at least U (1 ε)(1 β) (7) ε(1 β). Hece, the probability to have a successful step is at least p := ε(1 β)2 8β l(). Thus, the umber T of successful steps stochastically domiates a sum of a idepedet Beroulli distributed radom variables with parameter p, ad we obtai that E[T ] ap = ε(1 β)2 2 l() 8β l() = ε(1 β)2 16β (8) l() 2, (9)

7 7 ad by heroff bouds [ (1 β) Pr T ] [ ( = Pr T 1 (1 [ ( (9) Pr T 1 (1 l() 2 /=o(1) E[T ]/8 e (9) = e Ω( l() 2 ) = o(1). 16β ) ) l()2 ε(1 β) 2 (1 β)ε 16β ) ) E[T ]] 16β l()2 (1 β)ε Pr [T 1 ] 2 E[T ] ] l() 2 Hece, we a.a.s. have at least (1 β)/ successful steps, which together with the cosideratios aroud (8) cotradicts the assumptio to have o compoet of size (1 ε)(1 β). Refereces [1] D. Achlioptas, R. M. D Souza, ad J. Specer. Explosive percolatio i radom etworks. Sciece, 323(5920): , [2] T. Bohma ad A. Frieze. Avoidig a giat compoet. Radom Structures Algorithms, 19(1):75 85, [3] Y. S. ho, B. Kahg, ad D. Kim. luster aggregatio model for discotiuous percolatio trasitios. Phys. Rev. E, 81(3):030103, Mar [4] Y. S. ho, J. S. Kim, J. Park, B. Kahg, ad D. Kim. Percolatio trasitios i scale-free etworks uder the Achlioptas process. Phys. Rev. Lett., 103(13):135702, Sep [5] R. A. da osta, S. N. Dorogovtsev, A. V. Goltsev, ad J. F. F. Medes. Explosive percolatio trasitio is actually cotiuous. Phys. Rev. Lett., 105(25):255701, Dec [6] R. M. D Souza ad M. Mitzemacher. Local cluster aggregatio models of explosive percolatio. Phys. Rev. Lett., 104(19):195702, May [7] D. Dubhashi, A. Pacoesi, ad ambridge Uiversity Press. ocetratio of measure for the aalysis of radomized algorithms. ambridge Uiversity Press, [8] P. Erdős ad A. Réyi. O the evolutio of radom graphs. Magyar Tud. Akad. Mat. Kutató It. Közl., 5:17 61, [9] S. Jaso, T. Luczak, ad A. Ruciński. Radom graphs. Wiley-Itersciece, New York, [10] F. Radicchi ad S. Fortuato. Explosive percolatio: A umerical aalysis. Phys. Rev. E, 81(3):036110, Mar [11] O. Riorda ad L. Warke. Achlioptas process phase trasitios are cotiuous. arxiv: v1, [12] J. Specer ad N. Wormald. Birth cotrol for giats. ombiatorica, 27(5): , 2007.