O the degree distributio of a growig etwork model Lida Farczadi Uiversity of Waterloo lidafarczadi@gmailcom Nicholas Wormald Moash Uiversity ickwormald@moashedu arxiv:14010933v1 [mathco] 5 Ja 2014 Abstract I this ote we make some specific observatios o the distributio of the degree of a give vertex i certai model of radomly growig etworks The rule for etwork growth is the followig Startig with a iitial graph of miimum degree at least k, ew vertices are added oe by oe Each ew vertex v first chooses a radom vertex w to joi to, where the probability of choosig w is proportioal to its degree The k edges are added from v to radomly chose eighbours of w 1 Itroductio I this ote we make some specific observatios o the distributio of the degree of a give vertex i certai model of radomly growig etworks Fix a iteger k 2 We start with a seed graph G 1 k cosistig of oe vertex v 1 with k loops For t 2 give G t 1 k we obtai G t k as follows: we add a ew vertex v t which coects first to a existig vertex v i chose by preferetial attachmet ad the to k 1 eighbours of v i chose uiformly at radom We let V t ad E t deote the vertex set ad edge set of G t k Note that V t = t ad E t = kt We deote by d t v i the degree of vertex v i i G t k We defie N tv i to be the set of eighbours of vertex v i i G t k We begi with a simple derivatio of the expected value of d t v i i Sectio 21, the describe some closely related existig results, ad apply them to get a more precise descriptio of the distributio of d t v i i Sectio 24 It is straightforward to modify our results for ay give iitial seed graph 2 Degree distributio 21 Expected degree of a give vertex Fix vertex v i We wat to study d v i, the degree of vertex i at the th step i the process For t i we have Takig expectatios of both sides gives Pv t+1 coects to v i d t v i = d tv i 2 E t + = d tv i 2kt = d tv i v j N tv i +d t v i d t v j k 2 E t d t v j k 1 2kt Pv t+1 coects to v i = E[d tv i ] Supported by a ARC Australia Laureate Fellowship Research supported partly by NSERC 1
We the have E[d t+1 v i ] = E[d t v i ]+Pv t+1 coects to v i = E[d t v i ]+ E[d tv i ] = 1+ 1 E[d t v i ] Sice each vertex has degree k whe it jois the graph we have E[d i v i ] = k We obtai for 1 i E[d v i ] = k 1 22 LCD model of Bollobás ad Riorda t=i 1+ 1 = kγiγ+1/2 ΓΓi + /2 = k /i1+o1/i The LCD model of Bollobás ad Riorda ca be described as follow: start with G 1 1 the graph with oe vertex ad oe loop; for t 2 give G t 1 1 obtai G t 1 by addig oe vertex v t ad oe edge coectig v t to a existig vertex v i chose radomly with probability give by Pv i = s { dt 1s 1 1 1 if 1 s t 1 if s = t The for a give parameter k > 1 obtai G k by first costructig Gk 1 o vertices v 1,v 2,,v k usig the process described above The idetify vertices v 1,,v k to form vertex v 1 of G k, vertices v k+1,,v 2k to form vertex v 2, ad so o We observe that both the k-eighbour model ad the Bollobás-Riorda model satisfy the followig coditio Pv t+1 coects to v i d t v i = d t v i t j=1 d tv j This is kows as the Barabási-Albert BA descriptio Hece the degree of a give vertex has the same distributio i both models I particular the followig result cocerig the degree sequece of Bollobás ad Riorda [BRS + 01], [BR03] applies to the k-eighbour model as well Theorem 1 Let N d be the umber of vertices of degree d i G k ad defie αk,d = 2kk +1 dd+1d+2 The for a fixed ǫ > 0 ad 0 d 1/15 the followig holds with high probability 1 ǫαk,d N d 1+ǫαk,d 23 Geeral preferetial attachmet models of Ostroumova et al We ca obtai results about the k eighbour model by observig that it belogs to a certai class of geeral preferetial attachmet models Specifically Ostroumova et al [ORS12] defie the PA-class by cosiderig all radom graph models Gk that fit the followig descriptio: G k is a graph with vertices ad k edges obtaied from the followig radom graph process: start at time 0 with a arbitrary seed graph G 0 k with 0 vertices ad k 0 edges; at time t obtai the graph G t k from Gt 1 k by addig a ew vertex ad k edges coectig this vertex to some k vertices of G t 1 k 2
The Gk belogs to the class PA-class if it satisfies the followig coditios for some costats A ad B: P d t+1 v i = d t v i G t d t v i k = 1 A B 1 +O dt v i 2 2 1 P d t+1 v i = d t v i +1 G t d t v i k = A +B 1 +O dt v i 2 2 2 P d t+1 v i = d t v i +j G t dt v i 2 k = O 2 2 j k 3 P d t+1 v t+1 = k +j G t 1 k = O 1 j k 4 The we ca observe that our model belogs to this PA-class with parameters A = 1/2 ad B = 0 The the followig two results from [ORS12] apply to our model Theorem 2 Let N d be the umber of vertices of degree d i G k ad θx be a arbitrary fuctio such that θx < X There exists a costat C > 0 such that for ay d k we have E[N d] = αk,d +θcd 4 where αk,d = 2kk+1 dd+1d+2 2kk +1d 3 Theorem 3 For ay δ > 0 there exists a fuctio ψ = o such that for ay k d 1 8 δ 4 lim P N d E[N d] ψ d 3 = 0 24 Ur models We ca obtai the distributio of d v i by usig a ur model Our ur cotais balls of two colours: white ad black White balls represet edge-eds icidet with vertex i ad black balls represet edge-eds ot icidet with vertex i Suppose that the ur iitially has a 0 white balls ad b 0 black balls where a 0 = k b 0 = 2i k t 0 = a 0 +b 0 = 2i At each step, oe ball is draw radomly from the ur If the draw ball is white, replace it ad put a additioal α white ad α black If it is black, replace ad put more black balls We ow itroduce some relevat results about ur models from Flajolet et al [FDP06] Cosider a triagular ur with replacemet matrix α α 0 LetH a0 a b 0 b betheumberofhistoriesoflegththatstarticofiguratioa0,b 0 adedicofiguratio a,b Note that a+b = a 0 +b 0 + The the geeratig fuctio of ur histories is defied as Hx,1,z :=,a H a0 a b 0 b x az! 3
ad is give by Hx,1,z = x a0 1 z b0/ 1 x α 1 1 z α/ a 0/α Lettig := 1 z 1/ we have from [FDP06, Equatio 74] Hx,1,z = x a0 b0 1 x α 1 α a 0/α The reader may otice that [FDP06] is ot etirely cosistet, but i that paper, the balls of the first type do always correspod to the first row of the replacemet matrix Now let A m be the umber of white balls i the ur after m trials The probability that this equals a 0 +xα for some 0 x m is give by [FDP06, Equatio 75] as x+ a 0 x PA m = a 0 +xα = α 1 x 1 i [m+b0 αi/ 1] m x i [m+t 0 / 1] m i=0 where [ ] m deotes fallig factorial We ext explai why the ur results apply to the radom etwork process At ay give step, let s deote the total degrees of the vertices If a vertex v has degree d the the probability it is chose as the first vertex is d/s The probability it is chose as the secod vertex via ay give oe of its d icidet edges is k 1/s Hece, the probability it receives a ew edge is dk/s This is similar to the derivatio i Sectio 21 So we may use a ur with kd white balls ad s black balls At each step, s icreases by 2k, whilst kd icreases by k if v receives a ew edge, ad by 0 otherwise Hece, at each step the umber of white balls is k times the degree of v, ad the parameters are α = k ad = 2k The iitial umber of white balls, a 0, is k times the iitial degree of the vertex, ad the iitial umber of black balls is 2k times the iitial umber of edges SupposetheiitialgraphisacopyofK j Theforayofthej iitialverticeswehaveα = k, a 0 = kj 1, t 0 = jj 1 ad b 0 = t 0 a 0 = j kj 1 Also = 2k For each of those j vertices v i which iitially have degree j 1, there are m = j trials So the probability that the correspodig ur process fiishes with kj 1+x balls is x [ ] j + j 12 ku x+j 2 x Pd v i = j 1+x = 1 u 2k 1 j x u [ j +R 1] j u=0 where R = jj 1/2k O the other had, if i > j the t 0 = jj 1+2i jk, a 0 = k 2, b 0 = t 0 a 0, m = i ad so x x+k 1 x [ j +R u+k/2 1] i Pd v i = k +x = 1 u x u [ j +R 1] i u=0 We ca also obtai the momets of A directly from the geeratig fuctio of ur histories I particular, as explaied i [FDP06], the first ad secod momets are give by E[A ] = Γ+1Γ t 0 Hx,1,z Γ t 0 + [z ] x x=1 = Γ+1Γ t 0 Γ t 0 + [z ]a 0 t0+α E [ A 2 ] Γ+1Γ t 0 = Γ t 0 + 2 Hx,1,z [z ] x 2 x=1 = Γ+1Γ t 0 Γ t 0 + [z ] a 0 a 0 +α t0+2α a 0 α+1 t0+α 4
Performig coefficiet extractio gives Ad by liearity: [z ]a 0 t0+α = a 0 [z ]1 z t 0 +α t 0+α = a 0 + t 0+α = a 1 0 = a 0 Γ + t0+α usig [t ]1+at r = Γ+1Γ t 0+α r a r +r 1 usig = x usig = y 1 Γx+1 Γy +1Γx y +1 [z ] a 0 a 0 +α t0+2α a 0 α+1 t0+α = a 0 a 0 +α[z ] t0+2α a 0 α+1[z ] t0+α [ = a 0 a 0 +α Γ + t0+2α Γ+1Γ t 0+2α Pluggig these coefficiets back i the equatios for the first ad secod momet we obtai E[A ] = a 0Γ t 0 Γ + t 0+α Γ t 0+α Γ t0 + E [ [ A 2 ] a0 +αγ t 0 Γ + t 0+2α = a0 Γ t 0+2α Γ t0 + + α+1γ t 0 Γ + t ] 0+α Γ t 0+α Γ t0 + α+1 Γ + t0+α Γ+1Γ t 0+α Applyig these results to our example where a 0 = k, b 0 = 2i k, α = 1, = 2 ad the umber of trials is i we obtai E[d i] = kγiγ+1/2 Γi+1/2Γ E [ d i 2] [ k +1ΓiΓ+1 = k Γi+1Γ kk +1 = +2E[d i] i + 2ΓiΓ+1/2 ] Γi+1/2Γ We ote that the value for the first momet matches the oe obtaied previously by solvig the simple recursio i Sectio 21 Refereces ] [BR03] B Bollobás ad OM Riorda Mathematical results o scale-free radom graphs Hadbook of graphs ad etworks: from the geome to the iteret, pages 1 34, 2003 [BRS + 01] BBollobás, ORiorda, JSpecer, G Tusády, et al The degreesequeceofascale-freeradom graph process Radom Structures & Algorithms, 183:279 290, 2001 [FDP06] P Flajolet, P Dumas, ad V Puyhaubert Some exactly solvable models of ur process theory Fourth Colloquium o Mathematics ad Computer Sciece Algorithms, Trees, Combiatorics ad Probabilities, 59 118, Discrete Math Theor Comput Sci Proc, AG, Assoc Discrete Math Theor Comput Sci, Nacy, 2006 [Gre11] EA Grechikov A estimate for the umber of edges betwee vertices of give degrees i radom graphs i the bollobás riorda model Moscow Joural of Combiatorics ad Number Theory, 12, 2011 5
[OG11] L Ostroumova ad E Grechikov The distributio of secod degrees i the bollob\ as riorda radom graph model Arxiv preprit arxiv:11085585, 2011 [ORS12] L Ostroumova, A Ryabcheko, ad E Samosvat Geeralized preferetial attachmet: tuable power-law degree distributio ad clusterig coefficiet Arxiv preprit arxiv:12053015, 2012 6