arxiv: v1 [physics.soc-ph] 25 Dec 2018

Transcription

1 Optimal resource allocatio for etwork fuctioality Amikam Patro 1 Reuve Cohe 2 1 Departmet of Physics Bar-Ila Uiversity ramat-ga 5292 Israel 2 Departmet of Mathematics Bar-Ila Uiversity ramat-ga 5292 Israel Dated: December ) arxiv: v1 [physics.soc-ph] 25 Dec 218 The traditioal approach to etwork robustess is based o comparig etwork parameters before ad after a evet of odes removal such as the chage i etwork diameter the chage i giat compoet size ad the existece of giat compoet. Recetly a ew approach to etwork robustess was preseted where the etwork fuctioality durig its etire life spa durig the ode removal evet) is cosidered. This approach cosiders odes removal due to agig a state where odes ca be degeerated ad become ofuctioal whe their survival time duratio their lifetime is passed. Accordigly a problem that has to be solved is: i the etwork desig stage how to divide a budget of lifetime betwee the etwork s odes such that the etwork fuctioality durig all the stages of odes removal due to agig is maximized. To date the problem has bee solved oly partially ad umerically. I this paper we solve the problem aalytically ad derive a criterio for choosig the right set of odes o which the total lifetime budget should be divided. We also fid aalytically the best way of dividig the lifetime budget betwee the chose set odes such that the etwork robustess with cosideratio to its fuctioality i the etire life spa is maximized. Itroductio The robustess of a etwork which is the capability of a etwork to keep beig fuctioal eve whe some of its odes or liks) are removed is a issue that has bee researched ad studied widely [1 23]. Traditioally the measuremet of etwork robustess has bee performed i states where a removal of odes is a result of a attack o the etwork. There are two types of attacks i which the research has bee focused. The first type is radom attack where the choice of the etwork s odes for removal is performed radomly without preferece of aode to be removedoverayotherode i the etwork. The secod type is targeted attack where i cotrast to radom attack the attack is targeted agaist the odes with the cetral role of maitaiig the coectivity of the etwork. There have bee may studies about etwork robustess that are based o percolatio theory [1 1]. Accordig to this approach there are two states i which a etwork ca be foud the subcritical state where the etwork is fragmeted to may small compoets each of them cotais relatively small umber of odes ad the supercritical state where there exists a relative big compoet which cotais a fiite fractio of the etire etwork s odes i.e. scales as ON) amed the giat compoet ad aside from it possibly exist small compoets. The trasitio betwee the supercritical state to the subcritical state where the giat compoet is fragmeted ito small compoets occurs whe a critical combiatio of the etwork s odes is removed. The fractio of the etwork s odes that has ot to be removed i order to guaratee the existece of the giat compoet is amed the percolatio threshold of the etwork deoted by p c. I [1 3] the percolatio thresholds for two kids of radom etworks that have bee studied widely were measured. For Erdős-Réyi ER) etwork which was the first model of radom etwork to be studied [24 25] where the ode s degree follows a Poisso distributio Pk) = e λλk k! ad most of the odes degrees are close totheexpected degreeλ itwasfoud thatuderradom attackthepercolatiothresholdisp c = 1 λ. ForScale-Free SF) etwork a topology that was foud i may real etworks [26 29] where the ode s degree follows powerlaw distributio Pk) k γ such that beside most of the odes that have a very small degree there are some odes with a very high degree that are amed the hubs of the etwork it was foud that uder radom attack for γ > 3 p c has a fiite value but for γ 3 p c approaches. This meas that eve whe almost all of the odes are removed the etwork is still fuctioal ad a giat compoet exists i the etwork. For targeted attack it was foud that for ER etwork p c is close to its value uder radom attack sice most of the odes degrees are close to λ. I cotrast for SF etwork it was foud thatudertargetedattackp c ishigh siceremovigoly the very small group of the etwork s hubs that have a critical role of maitaiig the etwork coectivity causes the etwork to be fragmeted. The percolatio model for etwork robustess was also geeralized for iterdepedet etworks [11 12]. I additio to measurig etwork robustess by the percolatio threshold there have bee other parameters which were defied for measurig etwork robustess like chage i the etwork diameter which is the average short path legth betwee two radomly chose odes [ ] the relative size of the giat compoet which is the ratio betwee the giat compoet size after ad before the attack [15] chage i the betweeess cetrality which is the average of umber of shortest paths betwee pairs of the etwork s odes pass-

2 2 ig through radomly chose ode [16] coectivity loss which is the average of the ratio betwee the umber of odes coected to radomly chose ode after ad before the odes removal [17] ad f-robustess which is the umber of odes to be removed such that the giat compoet size is reduced by factor f f [1]) [18]. Efforts also have bee made to fid a topology for radom etwork o which the robustess is optimized both for radom ad targeted attack based o diameter chage [19] or percolatio threshold [6]. I additio to the studies about radom etworks robustess there are also some studies about robustess of real etworks like the North America power grid [17] ad optical commuicatio etworks [2]. A ew approach to etwork robustess has bee proposed i [21 22]. Accordig to it a measuremet of robustess should cosider the whole life spa of a etwork. This way is differet from the traditioal methods that measure robustess by comparig a fial ad iitial states after ad before a attack o the etwork took place respectively with o cosideratio to the itermediate state of the etwork fuctioality durig the attack. Accordig to the ew approach etwork robustess is measured by the itegral value of the fractio of odes i the giat compoet with respect to the umber of odes were removed from the etwork durig all the stages of the etwork attack. I [21 22] it was show by simulatios that choosig a small combiatio of pairs of the etwork s liks ad swappig the coectios of each pair for example the pair cosistig of a lik betwee odes i ad j ad a lik betwee odes k ad l is swapped ad becomes a lik betwee odes i ad k ad a lik betwee odes j ad l) ehaces the etwork robustess i.e. the itegral value metioed above very sigificatly. This approach was further elaborated i [23] where the effect of ode agig has bee added as a factor that should be cosidered i measurig etwork robustess. It was argued that each ode has a typical lifetime such thatfromthebegiigoftheetworklife aftertheode lifetime is passed the ode becomes degeerate ad i fact is removed from the etwork. Therefore ode agig is aother cause to ode removal i additio to the traditioal cause of a attack o the etwork. Accordigly a etwork survivability fuctio which is the itegral of the giat compoet size with respect to time durig the etire life spa of the etwork T has bee defied I = T Gt)dt 1) that is very similar i logic to the robustess itegral i [21 22] as was described above. For maximizig the etwork robustess accordig to it the followig problem has to be solved: i the etwork desig stage where the desiger allocates the odes with lifetime give a total budget of lifetime for all the etwork sodes what is the best way ofdividig the lifetime budget betwee them such that the survivability I i Eq. 1) is maximized? Ideed i [23] it was show umerically that whe a ode s lifetime is take from a distributio with expectatio that equals k α where k is the ode s degree the a critical value α c ca be foud depedet amog others o the etwork topology such that determiig the distributio expectatio to be k αc produces the optimized lifetime divisio with it the etwork survivability is maximized. However all the previous results are based o umerical solutios or o simulatios oly as metioed above ad to date there is still o aalytic solutio for a criterio or a algorithm accordig to it the etwork survivability would be maximized. I this paper we develop aalytically a solutio for the problem of etwork survivability maximizatio. We develop aalytically a criterio for choosig the right set of odes o which the total budget of lifetime should be divided ad also we fid aalytically the best way of dividig the lifetime betwee the odes of the chose set i such a way that the etwork survivability is maximized. Model Our model begis with a etwork with N odes ad a total budget of lifetime deoted by T. We choose a fiite descedig sequece of sets of the etwork odes A B C D E F.... Firstly we divide the etire lifetime budget T uiformly betwee set A odes. The we subtract a equal part of lifetime uits from each of the odes of set A B odes that belog to set A but ot belog to set B ad trasfer the etire lifetime was subtracted uiformly betwee set B odes. The we subtract a equal part of lifetime uits from each of the odes of set B C ad trasfer it uiformly betwee set C odes. We cotiue this lifetime trasfer i the same way from set C D to set D from set D E to set E etc. util we trasfer lifetime to the smallest set was chose. Fig. 1 is a basic illustratio of this model. Theory Presetatio of etwork survivability Forreadability we use Ψ = ABC... to ote the sets were chose where Ψ 1 = AΨ 2 = B etc. The umber of odes i set Ψ i ad set Ψ i Ψ i+1 would be oted Ψi ad Ψi Ψ i+1 respectively. The lifetime of set Ψ i would be oted t Ψi. The giat compoet size of set Ψ i whe set Ψ i is the oly set whose odes are fuctioal would be oted S Ψi. The lifetime subtracted from set Ψ i Ψ i+1 ad trasferred to set Ψ i+1 would be oted δt Ψi Ψ i+1. I geeral the etwork survivability with accordace

3 Due to the first subtractio of δt A B lifetime uits from each of the odes of set A B the lifetime of set A odes is reduced by δt A B ad becomest A = T δt A B ad the lifetime of set B odes icreases by δt A B lifetime uits over set A odes. Moreover due to the trasfer of the etire lifetime δt A B B that was subtracted from set A B each ode i set B receives i additio δt A B B B lifetime uits. Therefore the total lifetime of set B odes is t B = δt A A B B the summatio of δt A B ad δt A B A B B ). Therefore The etwork survivability after the first lifetime trasfer was performed is ) T I = S A δt A B = T S A +δt A B +S B δt A B SA S B B B 3 ). 4) FIG. 1: Basic model illustratio for choice of three sets A B C: a) Uiform divisio of the etire lifetime budget betwee the odes of set A. Node s lifetime is represeted by the height of the blue bars. b) Trasferrig lifetime from set A B to set B. Set B odes are those with the red bars above the blue bars. Red bars height represets the lifetime of set B odes after the lifetime trasfer to this set was performed. Blue bars height represets the remaied lifetime of set A. c) Trasferrig lifetime from set B C to set C. Set C odes are those with the gree bars above the red ad blue bars. Gree bars height represets the lifetime of set C odes after the lifetime trasfer to this set was performed. to its defiitio i Eq. 1) is I = S A t A +S B t B +S C t C +... = i S Ψi t Ψi. 2) I our model each stage of lifetime trasfer betwee sets chages the parameters of Eq. 2) ad thus chages the etwork survivability. Due to that i order to fid the fial etwork survivability after the completio of lifetime trasfers we calculate below the survivability i stages where i each stage we cosider aother lifetime trasfer i the sequece of trasfers. We begi firstly with the first stage of dividig the etire lifetime budget T uiformly betwee set A odes such that t A = T. Therefore the etwork survivability is I = T S A. 3) I similar way due to the secod trasfer of δt B C uits of lifetime from each of the odes of set B C to set C we get t B = δt A B B δt B C ad t C = δt B B C C ad the etwork survivability is S A SB I = T +δt A B S ) A + B +δt B C B SC C S B B ). 5) I geeral after the completio of a lifetime trasfer δt Ψi Ψ i+1 from each of the odes of set Ψ i Ψ i+1 to Ψi 1 set Ψ i+1 we get t Ψi = δt Ψi 1 Ψ i Ψi δt Ψi Ψ i+1 ad t Ψi+1 = δt Ψi Ψi Ψ i+1 Ψi+1. The fial etwork survivability after the last lifetime trasfer to the smallest set was chose is completed is I = T S Ψ1 Ψ1 + i SΨi+1 δt Ψi Ψi+1 Ψi S ) Ψ i. 6) Ψi+1 Ψi Figure 2 is a illustratio of the first three lifetime trasfers as described above for a etwork cotais 1 odes. For the sake of coveiece from ow o we develop the theory for oly four sets A B C D but the geeralizatio to ay umber of sets is simple as was show above. Sice δt A B is subtracted from the origial set A lifetime T ad δt B C is subtracted from the origial set B lifetime δt A A B B ad δt C D is subtracted from the origial set C lifetime δt B B C C therefore δt A B T δt B C δt A B B δt C D δt B C B C. 7)

4 was trasferred to set D. Based o Eq. 6) the etwork survivability is [ SA SB I = T +α A +α B α A SC C S B B S ) A B ) SD +α C α B α A S C D C 4 )]. 9) For each set we defie Z S as the fractio betwee its giat compoet size ad its odes umber. We also defieu A Z B Z A U B Z C Z B adu C Z D Z C. Accordigly we get a fial expressio for the etwork survivability I = T [Z A +α A U A +α B U B +α C U C ))]. 1) Maximizatio of etwork survivability Accordig to Eq. 1) maximum etwork survivability would be achieved by maximizig the followig max{α A U A +α B U B +α C U C ))}. 11) FIG. 2: Model illustratio for etwork with 1 odes. Nodes are siged by umbers form 1 to 1.a) First stage The etire lifetime budget T is divided uiformly betwee set A odes i this case all the etwork s odes). Lifetime of set A T is represetedby the height of blue bars attached to its odes. Liks betwee the odes are draw at the bottom of the figure. b) Secod stage subtractio of δt A B lifetime uits from each of the odes of set A B odes umber ad 1. Set A remais with T δt A B lifetime uits ad set B odes umber ad 9 gais automatically δt A B lifetime uits greater tha set A. Nodes i set B receive also δt A B A B B lifetime uits due to the lifetime trasfer from set A B. Total lifetime of set B odes is t B = δt A A B B. Lifetime of set B is represeted by the height of red bars attached to its odes. c) Third stage similar to b) but whe lifetime trasfer is performed from set B C odes umber 4 6 ad 8 ) to set C odes umber 2 5 ad 9 ). Lifetime of set C is represeted by the height of gree bars attached to its odes. We defie α A α B α C 1. With Eq. 7) we get δt A B = α A T δt B C B = α B δt A B = α B α A T 8) δt C D C = α C δt B C B = α C α B α A T. α A is the fractio of set A B origial lifetime T was trasferred to set B. α B is the fractio of set B C origial lifetime δt A A B B was trasferred to set C. α C is the fractio of set C D origial lifetime δt B B C C The maximizatio is implemeted by determiig the α s values from iside out accordig to the rule: if the expressio multiplied by α is greater tha the α = 1 ad if it is less or equal to the α =. We itroduce ow a specific example: firstly we have to maximize α C U C. Assume that U C > which meas Z D Z C = SD D SC C > ) the α C is determied to be 1 ad α C U C becomes U C. Accordigly The ext expressio to be maximized is α B U B +U C ). Assume that U B +U C which meas Z C Z B +Z D Z C = Z D Z B = SD D SB B ) the α B is determied to be ad α B U B +U C ) becomes. Accordigly the last expressio that has to be maximized is α A U A. Assume that U A > which meas Z B Z A = SB B SA > ) the α A is determied to be 1 ad α A U A becomes U A. Therefore accordig to Eqs. 1) ad 11) we get for the maximum of etwork survivability I max = T Z A +U A ) = T Z B = T S B B 12) that meas that maximum survivability is achieved by allocatig the etire lifetime T uiformly betwee the odes of set B. Note also that amog the four sets the ratio S is maximal for set B as was show Z D > Z C Z B Z D ad Z B > Z A ) ad this is the reaso for choosig this set to be allocated by the etire lifetime. The previous result ca be preseted i aother way. Accordig to Eqs. 2) ad 9) we get for the sets life-

5 5 times the followigs t A = 1 α A ) T t B = α A 1 α B ) T B t C = α A α B 1 α C ) T C Thus we ca defie a vector t D = α A α B α C T D. 13) 1 α A α A 1 α B ) α A α B 1 α C ) α A α B α C. The vector compoets represet the lifetime of sets ABC ad D respectively as a fractio of T. With the previous rule that maximum of the etwork survivability is achieved by determiig the α s values to be be or 1 we get that all the vectorcompoets equal except oe compoet that equals 1 the compoet that related to the set i which the ratio S is maximum. For the previous example with α A = 1 α B = ad α C = 1 the vector becomes 1) where oly the secod compoet related to set B equals 1. This presets as was cocluded above that i order to maximize the etwork survivability we have to allocate uiformly the etire lifetime to set B odes oly. The above example of four sets ca be geeralized to ay fial umber of sets performig the same aalysis of etwork survivability as was implemeted above. Let us summarize the previous aalysis with the geeral rule we proved: for each give group of sets of etwork s odes all are potetially to be allocated with lifetime the followigs has to be performed i order to maximize the etwork survivability: 1. Choice of the set i which the ratio betwee the giat compoet size ad the umber of odes is maximal. 2. Uiform divisio of the etire lifetime budget betwee the odes i the chose set. Simulatios ad Results Simulatios were implemeted o Erdős-Réyi ER) ad Scale-Free SF) radom etworks ad o real etworks. For each realizatio of each simulatio of each of the above etworks firstly a budget of lifetime is distributed o-uiformly betwee the various odes due to a distributio rule. Due to that the sets A B C etc. as they were defied above i the model sectio are idetified. The the etwork survivability simulatio begis ad a time clock is activated whereby the differet sets termiate their life accordig to their lifetimes t A t B t C etc. The time duratio betwee two successive sets failures is defied as oe stage of the whole etwork failure ad we ote the last stage of the etwork failure by m. The for each oe of the etwork failure stages represeted by i i = m) the product betwee the giat compoet size S Ψi ad the curret stage time duratio t Ψi is calculated. The etwork survivability is calculated to be I = m i=1 S Ψ i t Ψi. I the secod part of each simulatio realizatio our algorithm is activated. For each set of the etwork s odes the ratio S is calculated. We choose the set o which this ratio value is maximized ad divide the whole lifetime budget uiformly betwee this set s odes. The the etwork survivability simulatio begis agai ad a time clock is activated. Obviously the etire etwork collapses after the chose set s lifetime is passed. The etwork survivability would be the product of the chose set s giat compoet size by the chose set s lifetime. Figures 3a 1 ) ad 3a 2 ) preset simulatio results for ER etwork with N = 1 4 odes ad mea degree k = 2. The ouiform lifetime allocatig rule is each ode receives lifetime with proportio to its degree. A plot of the giat compoet size vs. the time durig the etwork failure is preseted i Fig. 3a 1 ) both for the ouiform lifetime divisio blue) ad for our theory with uiform lifetime divisio red). I the ouiform lifetime divisio the etwork survivability is spaed over four stages each of them is preseted by a blue rectagle. The giat compoet size for each stage is writte above the relevat rectagle ad the time duratio of each stage is preseted by braces below the t axis that are spaed over the stage s lifetime ad is writte below the braces. For example for the first stage A S A =.79 ad t u A =.49 the superscript u otes ouiform). Note that set A i the first stage icludes all the etwork s odes except odes with degree that obviously are ot allocated with lifetime) set B i the secod stage after odes with degree k = 1 were failed lifetime distributio with proportio to ode s degree causes the odes with degree k = 1 to receive miimum lifetime) icludes odes with degrees equal or greater tha 2 set C i the third stage icludes odes with degrees equal or greater tha 3 ad accordigly set D icludes odes with degrees equal or greater tha 4. The ext stages are ot preseted i the figure sice whe stage D is termiated ad odes with degree 4 are failed the etwork is fragmeted ad the giat compoet is ot exists aymore i the etwork. Hece the cotributio of these stages to the etwork survivability is. Fially the etwork survivability for ouifrom lifetime allocatio is I = S A t u A +S Bt u B +S Ct u C +S Dt u D. The uiform lifetime divisio accordig to our theory is preseted by the red rectagle. As was described above our algorithm calculates for each of the sets ABC ad D the ratio S of the giat compoet size ad the umber of odes i the set. Our algorithm fids that the

6 I t I t FIG. 3: [a)b)] Failure stages of ER etwork: a) Nouiform lifetime divisio proportioal to ode degree: Network size is N = 1 4 odes with average degree k = 2. a 1) Graph of giat compoet size vs. time. Four blue rectagles represet four failure stages with ouiform lifetime divisio. Red rectagle represets sigle failure stage with uiform lifetime divisio. Giat compoet size is writte o top of the relevat rectagle for each failure stage. Time duratio is writte ad oted by braces exteded across the relevat rectagle below the t axis for each failure stage. a 2) Graph of accumulated etwork survivability vs. time. Blue bars ad red bars represet ouiform ad uiform lifetime divisio respectively. The diagram eds at a time poit where the etwork collapses. For each divisio method the total etwork survivability is writte o top of its last graph poit. b) Nouiform lifetime divisio accordig to percetages of total lifetime: Network size is N = 1 4 odes with average degree k = 2.5. Percetages out of the total lifetime odes with degree k = 1 1% odes with degrees k = 23 15% ad odes with degrees equal or greater tha 4 75%. b 1) Graph of giat compoet size vs. time similar to a 1). b 2) graph of accumulated etwork survivability vs. time similar to b 1). Averages were take over 1 realizatios. For coveiece the averages were take o etworks were geerated with maximum degree k max = 1. The total lifetime budget T was take to be 1 4 similar to the etwork size i order to ormalize the etwork survivability to be betwee ad 1. Isets i a 1) ad b 1): bar graph of the ratio S betwee the giat compoet size ad the umber of odes for each of the sets. Maximum for this ratio is writte o top of the bar of the set o which this ratio is maximal. maximum for this ratio i our case is obtaied for set B. This result is show i the iset of the figure which is a bar graph with poits for the sets ABC ad D o the x axis ad a bar for each of the poits whose height represets the ratio S for the relevat set. We see that the highest bar belogs to set B with a value of.97. Accordigly the lifetime budget is fully divided uiformly betwee set B odes ad the whole etwork fails after oe stage oly failure of set B odes. The giat compoet size i this oly stageis writte abovethe red rectagle S B =.58 ad the time duratio of this stage is t u B = 1.68 the superscript u otes uiform) ad is preseted by braces below the t axis. Obviously the etwork survivability i this way is I = S B t u B. Figure 3a 2 ) is a graph of the etwork survivability I vs. the time durig the etwork failure both for the ouiform lifetime divisio blue) ad the uiform lifetime divisio red). For each time t o the x axis a bar is draw whose height represets the value of the expressio t St )dt which is the accumulated etwork survivability from the begiig of the etwork failure util time t. We ca see that although i the ouiform divisio the etwork survivability spas a loger time period tha for the uiform divisio the total survivability of the uiform divisio.97) is greater tha for the ouiform divisio.84). This result illustrates the idea of our theory that for each proposal of ouiform lifetime divisio betwee the etwork s odes there is other proposal of uiform lifetime divisio o which the etwork survivability is greater. Figures 3b 1 ) ad 3b 2 ) preset simulatio results for ER etwork with N = 1 4 odes ad mea degree

7 I t I t FIG. 4: a) Failure stages of SF etwork: Network size is N = odes with γ = 2.6. Network was divided to five sets odes with degree 1 odes with degree 2 odes with degrees 3 ad 4 odes with degrees betwee 5 ad 8 ad odes with degrees equal or greater tha 9. Each set was allocated by 2% of the total lifetime budget. a 1) Graph of giat compoet size vs. time. a 2) Graph of accumulated etwork survivability vs. time. Averages were take over 1 realizatios. b) Failure stages of Power-Grid real etwork: Network size is N = 4941 odes. Oly the first 85% of odes i the data set were cosidered. Percetages out of the total lifetime odes with degree k = 1 3% odes with degrees k = 2 25% odes with degrees k = 34 25% ad odes with degrees equal or greater tha 5 2%. b 1) Graph of giat compoet size vs. time. b 2) Graph of accumulated etwork survivability vs. time. Isets i a 1) ad b 1): The ratio S for each of the sets similar to the isets of Figs. 3a 1) ad 3b 1). k = 2.5. I this simulatio the ouiform lifetime divisio was implemeted accordig to percetages of the total lifetime such that 1% of the total lifetime was divided uiformly betwee the odes with degree 1 15% of the total lifetime was divided uiformly betwee the odes with degrees 2 ad 3 ad the remaied 75% of the lifetime was divided uiformly betwee the odes with degrees equal or greater tha 4. A iterestig poit is that accordig to this divisio although the total lifetime budget for odes with degree 1 1%) is less tha for odes with degrees 2 ad 3 15%) the lifetime of a sigle ode with degree 1 is greater tha the lifetime of a sigle ode with degree 2 or 3. This is because the etwork cotais sigificatly more odes with degrees 2 ad 3 tha odes with degree 1. As a result set A cotais all the etwork s odes except odes with degree set B cotais odes with degree 1 ad with degrees equal or greater tha 4 after the odes with degrees 2 ad 3 of set A B were failed ad set C cotais odes with degrees equal or greater tha 4. Figure 3b 2 ) validates agai our theory. We ca see that agaist the ouiform lifetime divisio proposal with etwork survivability value of.87 we propose a uiform lifetime divisio i which the etwork survivability is greater with value of.96. Figures 4a 1 ) ad 4a 2 ) preset simulatio results for SF etwork with N = odes ad γ = 2.6. The ouiform lifetime divisio was implemeted accordig to percetages of the total lifetime. We divide the etwork s odes to five sets odes with degree 1 odes with degree 2 odes with degrees 3 ad 4 odes with degrees betwee 5 ad 8 ad odes with degrees equal or greater tha 9. Each set was allocated by 2% of the total lifetime budget. I Fig. 4a 1 ) we ca see the iterestig result where our algorithm fids that maximum survivability would be achieved by uiformly lifetime divisio o set E with the high degrees that accordig to SF etwork properties cotais a very small part of the etwork s odes. This result ca be explaied due to the fact that the odes i set E have relatively high degrees ad the probability that all of them are coected together i oe compoet is very high. Therefore all or

8 8 at least most of set E odes are part of the giat compoet of this set ad the ratio SE E is very high ad would be greater tha this ratio value for the other sets. Aother iterestig poit is the very log time duratio of etwork survivability with uiform lifetime divisio o set E that is represeted by the very log legth of the redrectagleifig. 4a 1 ). Thisisduetothefactthato oe had set E receives 2% of the total lifetime equally to the other sets while o the other had it cotais a very small combiatio of odes relative to the other sets. Therefore each ode i this set receives a relatively high lifetime amout ad this set s lifetime is very high. Figures 4b 1 ) ad 4b 2 ) preset simulatio results for the Power-Grid PG) etwork [3] as a demostratio of our theory o real etworks too. Due to the ature of the odes belogig to this etwork which are electric elemets this etwork demostrates a classical example for odes removal due to agig. Sice PG etwork is fully coected the all its odes are cotaied i its giat compoet its ratio value S is 1 ad the solutio for maximum survivability is very trivial dividig the total lifetime budget uiformly betwee all the etwork s odes. I order to show otrivial solutios for this etwork too we choose to cosider oly part of the etwork s odes such that igorig the other odes causes the etwork to be ot fully coected. Therefore the results of Figs. 4b 1 ) ad 4b 2 ) are related to PG etwork where oly the first 85% of the odes i the data set were cosidered. The ouiform lifetime divisio was implemeted accordig to percetages of the total lifetime as follows: odes with degree 1 3% odes with degree 2 25% odes with degrees 3 ad 4 25% ad odes with degrees equal or greater tha 5-2%. The results i Fig. 4b 2 ) show agai that also i PG real etwork dividig the lifetime uiformly o set A odes the set that cotais all the etwork s odes except odes with degree gives a etwork survivability value of.88 that is greater tha the survivability value of.72 with the ouiform lifetime divisio. Summary I this work we developed a method for maximizig the robustess of a etwork i the dyamic approach where the robustess measuremet is performed by cosiderig the etwork fuctioality durig the etire odes removal evet ad accordig to a survivability fuctio Eq. 1). We proved aalytically that for a give group of sets of odes all of them are possibly for lifetime allocatio i order to maximize the etwork survivability we oly have to allocate lifetime to the set for which the ratio betwee its giat compoet size ad its umber of odes is maximal where the allocatio should be performed uiformly betwee the odes of the chose set. We hope that these fidigs could be useful at the stage of etwork desig as a tool for improvig etwork survivability. [1] R. Cohe K. Erez D. Be-Avraham ad S. Havli Physical review letters ). [2] D. S. Callaway M. E. Newma S. H. Strogatz ad D. J. Watts Physical review letters ). [3] R. Cohe K. Erez D. Be-Avraham ad S. Havli Physical review letters ). [4] R. Pastor-Satorras ad A. Vespigai Phys. Rev. E ). [5] M. E. J. Newma eprit arxiv:cod-mat/ ) cod-mat/2228. [6] G. Paul T. Taizawa S. Havli ad H. E. Staley The Europea Physical Joural B-Codesed Matter ad Complex Systems ). [7] L. K. Gallos R. Cohe P. Argyrakis A. Bude ad S. Havli Phys. Rev. Lett ). [8] R. Cohe ad S. Havli Complex Networks: Structure Robustess ad Fuctio Cambridge Uiversity Press 21). [9] F. Moroe ad H. A. Makse Nature ). [1] A. Braustei L. Dall Asta G. Semerjia ad L. Zdeborová Proceedigs of the Natioal Academy of Scieces ). [11] S. Buldyrev R. Parshai G. Paul H. Staley ad S. Havli Nature ). [12] S. V. Buldyrev N. W. Shere ad G. A. Cwilich Phys. Rev. E ). [13] R. Albert H. Jeog ad A.-L. Barabási Nature Lodo) ). [14] J.-H. Kim K.-I. Goh B. Kahg ad D. Kim Phys. Rev. Lett ). [15] A. E. Motter ad Y.-C. Lai Phys. Rev. E ). [16] P. Holme B. J. Kim C. N. Yoo ad S. K. Ha Phys. Rev. E ). [17] R. Albert I. Albert ad G. L. Nakarado Phys. Rev. E ). [18] P. Holme ad J. Zhao Phys. Rev. E ). [19] B. Shargel H. Sayama I. R. Epstei ad Y. Bar-Yam Physical review letters ). [2] L. Jahke J. W. Katelhardt R. Berkovits ad S. Havli Phys. Rev. Lett ). [21] C. M. Scheider A. A. Moreira J. S. Adrade S. Havli ad H. J. Herrma Proceedigs of the Natioal Academy of Scieces ). [22] H. J. Herrma C. M. Scheider A. A. Moreira J. S. A. Jr ad S. Havli Joural of Statistical Mechaics: Theory ad Experimet 211 P ). [23] Y. Li A. Patro S. Guo R. Kag D. Li S. Havli ad R. Cohe EPL Europhysics Letters) ). [24] P. Erdős ad A. Réyi Publ. Math. Debrece ). [25] P. Erdős ad A. Réyi Publ. Math. Ist. Hugar. Acad. Sci ). [26] S. Reder The Europea Physical Joural B - Codesed Matter ad Complex Systems ). [27] R. Albert H. Jeog ad A.-L. Barabási Nature 41

9 ). [28] A.-L. Barabási ad R. Albert Sciece ). [29] M. E. J. Newma Proceedigs of the Natioal Academy of Scieces ). [3] 9