Running Time Analysis of the (1+1)-EA for OneMax and LeadingOnes under Bit-wise Noise

Size: px
Start display at page:

Download "Running Time Analysis of the (1+1)-EA for OneMax and LeadingOnes under Bit-wise Noise"

Transcription

1 Ruig Time Aalysis of the +-EA for OeMax ad LeadigOes uder Bit-wise Noise Chao Qia Uiversity of Sciece ad Techology of Chia Hefei 3007, Chia Wu Jiag Uiversity of Sciece ad Techology of Chia Hefei 3007, Chia ABSTRACT Previous ruig time aalyses of evolutioary algorithms EAs i oisy eviromets ofte studied the oe-bit oise model, which flips a radomly chose bit of a solutio before evaluatio. I this paper, we study a atural extesio of oe-bit oise, the bit-wise oise model, which idepedetly flips each bit of a solutio with some probability. We aalyze the ruig time of the +-EA solvig OeMax ad LeadigOes uder bit-wise oise for the first time, ad derive the rages of the oise level for polyomial ad superpolyomial ruig time bouds. The aalysis o LeadigOes uder bit-wise oise ca be easily trasferred to oebit oise, ad improves the previously kow results. CCS CONCEPTS Theory of computatio Theory of radomized search heuristics; KEYWORDS Noisy optimizatio, evolutioary algorithms, ruig time aalysis, computatioal complexity ACM Referece format: Chao Qia, Chao Bia, Wu Jiag, ad Ke Tag. 07. Ruig Time Aalysis of the +-EA for OeMax ad LeadigOes uder Bitwise Noise. I Proceedigs of GECCO 7, Berli, Germay, July 5-9, 07, 8 pages. DOI: INTRODUCTION Evolutioary algorithms EAs have bee widely applied to solve real-world optimizatio tasks, where the exact objective i.e., fitess evaluatio of cadidate solutios is ofte Permissio to make digital or hard copies of all or part of this work for persoal or classroom use is grated without fee provided that copies are ot made or distributed for profit or commercial advatage ad that copies bear this otice ad the full citatio o the first page. Copyrights for compoets of this work owed by others tha ACM must be hoored. Abstractig with credit is permitted. To copy otherwise, or republish, to post o servers or to redistribute to lists, requires prior specific permissio ad/or a fee. Request permissios from permissios@acm.org. GECCO 7, Berli, Germay 07 ACM /7/07... $5.00 DOI: Chao Bia Uiversity of Sciece ad Techology of Chia Hefei 3007, Chia biacht@mail.ustc.edu.c Ke Tag Uiversity of Sciece ad Techology of Chia Hefei 3007, Chia ketag@ustc.edu.c impossible, while we ca obtai oly a oisy oe [3, 3]. However, previous theoretical aalyses of EAs maily focused o oise-free optimizatio, where the fitess evaluatio is exact. For ruig time aalysis, a leadig theoretical aspect [, 4], oly a few pieces of work o oisy evolutioary optimizatio have bee reported. Droste [6] first aalyzed the +-EA o the OeMax problem i the presece of oe-bit oise ad showed the maximal oise probability p = log / allowig a polyomial ruig time, where is the problem size. Gieße ad Kötzig [] recetly studied the LeadigOes problem, ad proved that the expected ruig time is polyomial if p /6e ad expoetial if p = /. For iefficiet optimizatio of the +-EA uder high oise levels, some implicit mechaisms of EAs were proved to be robust to oise. I [], it was show that the µ+- EA with a small populatio of size Θlog ca solve Oe- Max i polyomial time eve if the probability of oe-bit oise reaches. The robustess of populatios to oise was also proved i the settig of o-elitist EAs [4, 7]. However, Friedrich et al. [8] showed the limitatio of populatios by provig that the µ+-ea eeds super-polyomial time for solvig OeMax uder additive Gaussia oise N 0, σ with σ 3. This difficulty ca be overcome by the compact geetic algorithm cga [8] ad a simple At Coloy Optimizatio ACO algorithm [9], both of which fid the optimal solutio i polyomial time with a high probability. The ability of explicit oise hadlig strategies was also theoretically studied. Qia et al. [9] proved that the threshold selectio strategy is robust to oise: the expected ruig time of the +-EA usig threshold selectio o OeMax uder oe-bit oise is always polyomial regardless of the oise level. For the +-EA solvig OeMax ad LeadigOes uder oe-bit or additive Gaussia oise, the resamplig strategy was show able to reduce the ruig time from expoetial to polyomial i high oise levels [8]. Akimoto et al. [] also proved that resamplig with a large sample size ca make optimizatio uder additive ubiased oise behave as optimizatio i a oise-free eviromet. The iterplay betwee resamplig ad implicit oise-hadlig mechaisms has bee statistically studied i [0]. The studies metioed above maily cosidered the oebit oise model, which flips a radom bit of a solutio before

2 GECCO 7, July 5-9, 07, Berli, Germay C. Qia et. al. Table : For the ruig time of the +-EA solvig OeMax ad LeadigOes uder prior oise models, the rages of oise parameters for a polyomial upper boud ad a super-polyomial lower boud are show below. +-EA bit-wise oise p, bit-wise oise, q oe-bit oise OeMax Olog /, ωlog / Olog /, ωlog / [] Olog /, ωlog / [6] LeadigOes Olog /, ωlog / Olog / 3, ωlog / /6e, = / []; Olog /, ωlog / evaluatio with probability p. However, the oise model, which ca chage several bits of a solutio simultaeously, may be more realistic ad eeds to be studied, as metioed i the first oisy theoretical work [6]. I this paper, we study the bit-wise oise model, which is characterized by a pair p, q of parameters. It happes with probability p, ad idepedetly flips each bit of a solutio with probability q before evaluatio. We aalyze the ruig time of the +-EA solvig OeMax ad LeadigOes uder bit-wise oise with two specific parameter settigs p, ad, q. The rages of p ad q for a polyomial upper boud ad a super-polyomial lower boud are derived, as show i the middle two colums of Table. For the +-EA o LeadigOes, we also trasfer the ruig time bouds from bit-wise oise p, to oe-bit oise by usig the same proof procedure. As show i the bottom right of Table, our results improve the previously kow oes []. The rest of this paper is orgaized as follows. Sectio itroduces some prelimiaries. The ruig time aalysis o OeMax ad LeadigOes is preseted i Sectios 3 ad 4, respectively. Sectio 5 cocludes the paper. PRELIMINARIES I this sectio, we first itroduce the optimizatio problems, evolutioary algorithms ad oise models studied i this paper, respectively, ad the preset the aalysis tools that we use throughout this paper.. OeMax ad LeadigOes I this paper, we use two well-kow pseudo-boolea fuctios OeMax ad LeadigOes. The OeMax problem as preseted i Defiitio. aims to maximize the umber of -bits of a solutio. The LeadigOes problem as preseted i Defiitio. aims to maximize the umber of cosecutive -bits coutig from the left of a solutio. Their optimal solutio is... briefly deoted as. It has bee show that the expected ruig time of the +-EA o OeMax ad LeadigOes is Θ log ad Θ, respectively [7]. Defiitio. OeMax. The OeMax Problem of size is to fid a bits biary strig x such that x = arg max x {0,} fx = i= xi. Defiitio. LeadigOes. The LeadigOes Problem of size is to fid a bits biary strig x such that x = arg max x {0,} fx = i. i= j= xj. Bit-wise Noise There are maily two kids of oise models: prior ad posterior [, 3]. The prior oise comes from the variatio o a solutio, while the posterior oise comes from the variatio o the fitess of a solutio. Previous theoretical aalyses ofte focused o a specific prior oise model, oe-bit oise. As preseted i Defiitio.3, it flips a radom bit of a solutio before evaluatio with probability p. However, i may realistic applicatios, oise ca chage several bits of a solutio simultaeously rather tha oly oe bit. We thus cosider the bit-wise oise model. As preseted i Defiitio.4, it happes with probability p, ad idepedetly flips each bit of a solutio with probability q before evaluatio. To the best of our kowledge, oly bit-wise oise with p = ad q [0, ] has bee recetly studied. Gieße ad Kötzig [] proved that for the +-EA o OeMax, the expected ruig time is polyomial if q = Olog / ad super-polyomial if q = ωlog /. I this paper, we will study two specific bit-wise oise models: p [0, ] q = ad p = q [0, ], which are briefly deoted as bit-wise oise p, ad bit-wise oise, q, respectively. Defiitio.3 Oe-bit Noise. Give a parameter p [0, ], let f x ad fx deote the oisy ad true fitess of a biary solutio x {0, }, respectively, the { f fx with probability p, x = fx with probability p, where x is geerated by flippig a uiformly radomly chose bit of x. Defiitio.4 Bit-wise Noise. Give parameters p, q [0, ], let f x ad fx deote the oisy ad true fitess of a biary solutio x {0, }, respectively, the { f fx with probability p, x = fx with probability p, where x is geerated by idepedetly flippig each bit of x with probability q..3 +-EA The +-EA as described i Algorithm is studied i this paper. For oisy optimizatio, oly a oisy fitess value f x istead of the exact oe fx ca be accessed, ad thus step 4 of Algorithm chages to be if f x f x. Note that the reevaluatio strategy is used as i [5, 6, ]. That is, besides evaluatig f x, f x will be reevaluated i each iteratio of the +-EA. The ruig time is usually defied as the umber of fitess evaluatios eeded to fid

3 Ruig Time Aalysis of the +-EA uder Bit-wise Noise GECCO 7, July 5-9, 07, Berli, Germay a optimal solutio w.r.t. the true fitess fuctio f for the first time [, 6, ]. ALGORITHM +-EA. Give a fuctio f over {0,} to be maximized, it cosists of the followig steps:. x := uiformly radomly selected from {0, }.. Repeat util the termiatio coditio is met 3. x :=flip each bit of x idepedetly with prob. /. 4. if fx fx 5. x := x..4 Aalysis Tools The process of the +-EA solvig OeMax or LeadigOes ca be directly modeled as a Markov chai {ξ t} + t=0. We oly eed to take the solutio space {0, } as the chai s state space i.e., ξ t X = {0, }, ad take the optimal solutio as the chai s optimal state i.e., X = { }. Give a Markov chai {ξ t} + t=0 ad ξˆt = x, we defie its first hittig time FHT as τ = mi{t ξˆt+t X, t 0}. The mathematical expectatio of τ, E[[τ ξˆt = x]] = + i=0 ip τ = i, is called the expected first hittig time EFHT startig from ξˆt = x. If ξ 0 is draw from a distributio π 0, E[[τ ξ 0 π 0]] = x X π0xe[[τ ξ0 = x]] is called the EFHT of the Markov chai over the iitial distributio π 0. Thus, the expected ruig time of the +-EA startig from ξ 0 π 0 is equal to + E[[τ ξ 0 π 0]], where the term correspods to evaluatig the iitial solutio, ad the factor correspods to evaluatig the offsprig solutio x ad reevaluatig the paret solutio x i each iteratio. Note that we cosider the expected ruig time of the +-EA startig from a uiform iitial distributio i this paper. I the followig, we give three drift theorems that will be used to derive the EFHT of Markov chais i the paper. LEMMA.5 ADDITIVE DRIFT []. Give a Markov chai {ξ t} + t=0 ad a distace fuctio V x, if for ay t 0 ad ay ξ t with V ξ t > 0, there exists a real umber c > 0 such that E[[V ξ t V ξ t+ ξ t]] c, the the EFHT satisfies that E[[τ ξ 0]] V ξ 0/c. LEMMA.6 SIMPLIFIED DRIFT [5, 6]. Let X t, t 0, be real-valued radom variables describig a stochastic process. Suppose there exists a iterval [a, b] R, two costats δ, ɛ > 0 ad, possibly depedig o l := b a, a fuctio rl satisfyig rl = ol/ logl such that for all t 0 the followig two coditios hold:. E[[X t X t+ a < X t < b]] ɛ,. P X t+ X t j X t > a rl + δ j for j N 0. The there is a costat c > 0 such that for T := mi{t 0 : X t a X 0 b} it holds P T cl/rl = Ωl/rl. LEMMA.7 SIMPLIFIED DRIFT WITH SELF-LOOPS [0]. Let X t, t 0, be real-valued radom variables describig a stochastic process. Suppose there exists a iterval [a, b] R, two costats δ, ɛ > 0 ad, possibly depedig o l := b a, a fuctio rl satisfyig rl = ol/ logl such that for all t 0 the followig two coditios hold:. a<i<b : E[[X t X t+ X t =i]] ɛ P X t+ i X t =i,. i > a, j N 0 : P X t+ X t j X t =i rl P Xt+ i Xt =i. +δ j The there is a costat c > 0 such that for T := mi{t 0 : X t a X 0 b} it holds P T cl/rl = Ωl/rl. 3 THE ONEMAX PROBLEM I this sectio, we aalyze the ruig time of the +- EA o OeMax uder bit-wise oise. Note that for bit-wise oise, q, it has bee proved that the maximal value of q allowig a polyomial ruig time is log /, as show i Theorem 3.. THEOREM 3.. [] For the +-EA o OeMax uder bitwise oise, q, the expected ruig time is polyomial if q = Olog / ad super-polyomial if q = ωlog /. For bit-wise oise p,, we prove i Theorems 3.4 ad 3.5 that the maximum value of p allowig a polyomial ruig time is log /. Istead of usig the origial drift theorems, we apply the upper ad lower bouds for the +-EA o stochastic OeMax i []. Let x k deote ay solutio with k umber of -bits, ad f x k deote its oisy objective value, which is a radom variable. Lemma 3. ituitively meas that if the probability of recogizig the true better solutio by oisy evaluatio is large, the ruig time ca be polyomially upper bouded. O the cotrary, Lemma 3.3 shows that if the probability of makig a right compariso is small, the ruig time ca be expoetially lower bouded. Both of them are proved by applyig stadard drift theorems, ad ca be used to simplify our aalysis. LEMMA 3.. [] Suppose there is a positive costat c /5 ad some < l / such that k < : P f x k < f x k+ l ; k < l : P f x k < f x k+ c k, the the +-EA optimizes f i expectatio i O log + Ol iteratios. LEMMA 3.3. [] Suppose there is some l /4 ad a costat c 6 such that l k < : P f x k < f x k+ c k, the the +-EA optimizes f i Ωl iteratios with a high probability. THEOREM 3.4. For the +-EA o OeMax uder bit-wise oise p,, the expected ruig time is polyomial if p = Olog /. PROOF. We prove it by usig Lemma 3.. For ay positive costat b, suppose that p b log /. We set the two parameters i Lemma 3. as c = mi{ b log, b} ad l =, ]. 5 c

4 GECCO 7, July 5-9, 07, Berli, Germay C. Qia et. al. For ay k <, f x k f x k+ implies that f x k k + or f x k+ k, either of which happes with probability at most p. By the uio boud, we get P f x k f x k+ p For ay k < l, we easily get b log = lc l. P f x k f x k+ lc < c k. By Lemma 3., we kow that the expected ruig time is O log + Ob log /c, i.e., polyomial. THEOREM 3.5. For the +-EA o OeMax uder bit-wise oise p,, the expected ruig time is super-polyomial if p = ωlog / ωlog / ad expoetial if p = Olog /. PROOF. We use Lemma 3.3 to prove it. Let c = 6. The case p = ωlog / ωlog / is first aalyzed. For ay positive costat b, let l = b log. For ay k l, we get P f x k f x k+ P f x k = k P f x k+ k. To make f x k = k, it is sufficiet that the oise does ot happe, i.e., P f x k = k p. To make f x k+ k, it is sufficiet to flip oe -bit ad keep other bits uchaged by oise, i.e., P f x k+ k p k+. Thus, P f x k f x k+ p p k + e Sice c k c l = ωlog /. = cb log, the coditio of Lemma 3.3 holds. Thus, the expected ruig time is Ωb log where b is ay costat, i.e., super-polyomial. For the case p = Olog /, let l =. We use aother lower boud p for P f x k = k, sice it is sufficiet that o bit flips by oise. Thus, we have P f x k f x k+ p p k + e = Ω. Sice c k c, the coditio of Lemma 3.3 holds. Thus, the expected ruig time is Ω, i.e., expoetial. 4 THE LEADINGONES PROBLEM I this sectio, we first aalyze the ruig time of the +- EA o LeadigOes uder bit-wise oise p, ad bit-wise oise, q, respectively. The, we trasfer the aalysis from bit-wise oise p, to oe-bit oise; the results are complemetary to the kow oes recetly derived i []. 4. Bit-wise Noise p, For bit-wise oise p,, we prove i the followig three theorems that the expected ruig time is polyomial if p = Olog / ad super-polyomial if p = ωlog /. Their proofs are accomplished by applyig additive drift aalysis, the simplified drift theorem with self-loops ad the simplified drift theorem, respectively. THEOREM 4.. For the +-EA o LeadigOes uder bitwise oise p,, the expected ruig time is polyomial if p = Olog /. PROOF. We use Lemma.5 to prove it. For ay positive costat b, suppose that p b log /. Let LOx deote the true umber of leadig -bits of a solutio x. We first costruct a distace fuctio V x as, for ay x with LOx = i, V x = + c + c i, where c = 4b log +. It is easy to verify that V x X = { } = 0 ad V x / X > 0. The, we ivestigate E[[V ξ t V ξ t+ ξ t = x]] for ay x with LOx < i.e., x / X. Assume that curretly LOx = i, where 0 i. We divide the drift ito two parts: positive E + ad egative E. That is, E[[V ξ t V ξ t+ ξ t = x]] = E + E. For the positive drift, we eed to cosider that the umber of leadig -bits is icreased. By mutatio, we have P LOx i + = i, sice it eeds to flip the i + -th bit which must be 0 of x ad keep the i leadig -bits uchaged. For ay x with LOx i +, f x < f x implies that f x i or f x i +. Note that, P f x i = p i, sice at least oe of the first i leadig -bits of x eeds to be flipped by oise; P f x i + = p i, 3 sice it eeds to flip the first 0-bit of x ad keep the leadig -bits uchaged by oise. By the uio boud, we get P f x f x = P f x < f x p i+ p i +, 4 where the last iequality is by p = Olog /. Furthermore, for ay x with V x i +, V x V x + c i+ + c i= c + c i. 5 By combiig Eqs., 4 ad 5, we have E + i c + c i c + c i, 6 where the last iequality is by i. e 3 For the egative drift, we eed to cosider that the umber of leadig -bits is decreased. By mutatio, we have P LOx i = i, 6 sice it eeds to flip at least oe leadig -bit of x. For ay x with LOx i where i, f x f x implies

5 Ruig Time Aalysis of the +-EA uder Bit-wise Noise GECCO 7, July 5-9, 07, Berli, Germay that f x i or f x i. Note that, P f x i p i, 7 sice for the first i bits of x, it eeds to flip the 0-bits whose umber is at least ad keep the -bits uchaged by oise; P f x i = p i, 8 sice at least oe leadig -bit of x eeds to be flipped by oise. By the uio boud, we get P f x f x p p i p i+.9 Furthermore, for ay x with LOx i, V x V x + c i. 0 By combiig Eqs. 6, 9 ad 0, we have E i p i + p + c i p e 3 + c i + c i. Thus, by subtractig E from E +, we have E[[V ξ t V ξ t+ ξ t = x]] + c i c 6 p 3 + c i 4b log + b log 6 3 6, where the secod iequality is by c = 4b log + ad p b log /. Note that V x + c e c = e 4b log + = e 4b. By Lemma.5, we get E[[τ ξ 0]] 6 e 4b = O 4b+, i.e., the expected ruig time is polyomial. THEOREM 4.. For the +-EA o LeadigOes uder bit-wise oise p,, if p = ωlog / o, the expected ruig time is super-polyomial. PROOF. We use Lemma.7 to prove it. Let X t = x 0 be the umber of 0-bits of the solutio x after t iteratios of the +-EA. Let c be ay positive costat. We cosider the iterval [0, c log ], i.e., the parameters a = 0 i.e., the global optimum ad b = c log i Lemma.7. The, we aalyze the drift E[[X t X t+ X t = i]] for i < c log. As i the proof of Theorem 4., we divide the drift ito two parts: positive E + ad egative E. That is, E[[X t X t+ X t = i]] = E + E. For the positive drift, we eed to cosider that the umber of 0-bits is decreased. For mutatio o x where x 0 = i, let X ad Y deote the umber of flipped 0-bits ad -bits, respectively. The, X Bi, ad Y B i,, where B, is the biomial distributio. Let P mutx, x deote the probability of geeratig x by mutatig x. To estimate a upper boud o E +, we assume that the offsprig solutio x with x 0 < i is always accepted. Thus, we have E + i P mutx, x i x 0 = k P X Y = k x : x 0 <i = i = i i k= k i j= j k= j=k k= P X = j P Y = j k k P X = j P Y = j k j= j P X = j = i. For the egative drift, we eed to cosider that the umber of 0-bits is icreased. We aalyze the i cases where oly oe -bit is flipped i.e., x 0 = i +, which happes with probability. Assume that LOx = k i. If the j-th where j k leadig -bit is flipped, the offsprig solutio x will be accepted i.e., f x f x if f x j ad f x j. Note that, P f x j = p+p j p j, where the equality is sice it eeds to keep the j leadig -bits of x uchaged, ad the last iequality is by p = o; P f x j = p j 3 = p j + j p e j pj 3, where the equality is sice at least oe of the first j leadig -bits of x eeds to be flipped by oise. Thus, we get P f x f x pj 6. 4 If oe of the i k o-leadig -bits is flipped, LOx = LOx = k. We ca use the same aalysis procedure as Eq. 4 i the proof of Theorem 4. to derive that P f x f x p k +, 5 where the secod iequality is by p = o. Combiig all the i cases, we get E k pj 6 + i k i+ i 6 j= pkk + + i k pk e 36 + i k. 6 By subtractig E from E +, we get E[[X t X t+ X t = i]] i pk 36 i k. 6 To ivestigate coditio of Lemma.7, we also eed to aalyze the probability P X t+ i X t = i. For X t+ i, it is ecessary that at least oe bit of x is flipped ad the offsprig x is accepted. We cosider two cases: at least oe of the k leadig -bits of x is flipped; the k leadig - bits of x are ot flipped ad at least oe of the last k bits is flipped. For case, the mutatio probability is k

6 GECCO 7, July 5-9, 07, Berli, Germay C. Qia et. al. ad the acceptace probability is at most p k+ by Eq. 9. For case, the mutatio probability is k k k ad the acceptace probability is at most. Thus, P X t+ i X t = i p + k. 7 Whe k < p, we have E[[X t X t+ X t = i]] i i k 6 k p/ 7c log k 6 4 p+ k 8, where the secod iequality is by k > p ad i < c log, the third iequality is by p = ωlog / ad the last is by k > p. Whe k p, we have E[[X t X t+ X t = i]] i pk 9 36 c log p 44 p 88 p + k, 576 where the secod iequality is by p = o ad i < c log, the third is by p = ωlog / ad the last is by k p. Combiig Eqs. 7, 8 ad 9, we get that coditio of Lemma.7 holds with ɛ =. 576 For coditio of Lemma.7, we eed to show P X t+ X t j X t = i rl P X +δ j t+ i X t = i for i. For P X t+ i X t = i, we aalyze the cases where oly oe bit is flipped. Usig the similar aalysis procedure as E, except that flippig ay bit rather tha oly -bit is cosidered here, we easily get P X t+ i X t = i pkk k 6. 0 For X t+ X t j, it is ecessary that at least j bits of x are flipped ad the offsprig solutio x is accepted. We cosider two cases: at least oe of the k leadig -bits is flipped; the k leadig -bits are ot flipped. For case, the mutatio j probability is at most k ad the acceptace probability is at most p k+ by Eq. 9. For case, the mutatio j probability is at most k k ad the acceptace j j probability is at most. Thus, we have P X t+ X t j X t = i k j p k + j + k k j j pkk+ 4 + k j pkk+ + k 44 j j By combiig Eq. 0 with Eq., we get that coditio of Lemma.7 holds with δ = ad rl = 44. Note that l = b a = c log. By Lemma.7, the expected ruig time is Ωc log, where c is ay positive costat. Thus, the expected ruig time is super-polyomial. THEOREM 4.3. For the +-EA o LeadigOes uder bitwise oise p,, the expected ruig time is expoetial if p = Ω. PROOF. We use Lemma.6 to prove it. Let X t = i be the umber of 0-bits of the solutio x after t iteratios of the +-EA. We cosider the iterval i [0, / ]. To aalyze the drift E[[X t X t+ X t = i]] = E + E, we use the same aalysis procedure as Theorem 4.. For the positive drift, we have E + i = o. For the egative drift, we re-aalyze Eqs. 4 ad 5. P f x j p j Thus, Eq. 4 becomes From Eqs. ad 3, we get that ad P f x j pj. 3 P f x f x p j 3 j. For Eq. 5, we eed to aalyze the acceptace probability for LOx = LOx = k. Sice it is sufficiet to keep the first k + bits of x ad x uchaged i oise, Eq. 5 becomes P f x f x p k+ p k+. 3 By applyig the above two iequalities to Eq. 6, we have k E p j j + + i k k =Ω, e 3 j= where the equality is by p = Ω. Thus, E + E = Ω. That is, coditio of Lemma.6 holds. Sice it is ecessary to flip at least j bits of x, we have P X t+ X t j X t j j j!, j which implies that coditio of Lemma.6 holds with δ = ad rl =. Note that l = /. Thus, by Lemma.6, the expected ruig time is expoetial. 4. Bit-wise Noise, q For bit-wise oise, q, we prove i the followig three theorems that the expected ruig time is polyomial if q = Olog / 3 ad super-polyomial if q = ωlog /. The proof idea is similar to that for bit-wise oise p,. The mai differece led by the chage of oise is the probability of acceptig the offsprig solutio, i.e., P f x f x. THEOREM 4.4. For the +-EA o LeadigOes uder bitwise oise, q, the expected ruig time is polyomial if q = Olog / 3. PROOF. The proof is very similar to that of Theorem 4.. The chage of oise oly affects the probability of acceptig the offsprig solutio i the aalysis. For ay positive costat b, suppose that q b log / 3. For the positive drift E +, we eed to re-aalyze P f x f x i.e., Eq. 4 i the proof of Theorem 4. for the paret x with LOx = i ad the offsprig x with LOx i +. By bit-wise oise, q, Eqs. ad 3 chage to P f x i = q i ; P f x i+= q i q. Thus, by the uio boud, Eq. 4 becomes P f x f x q i + q i q 4 = q i+ qi + /,

7 Ruig Time Aalysis of the +-EA uder Bit-wise Noise GECCO 7, July 5-9, 07, Berli, Germay where the last iequality is by q = Olog / 3. For the egative drift E, we eed to re-aalyze P f x f x i.e., Eq. 9 i the proof of Theorem 4. for the paret x with LOx = i where i ad the offsprig x with LOx i. By bit-wise oise, q, Eqs. 7 ad 8 chage to P f x i q q i, P f x i = q i. Thus, by the uio boud, Eq. 9 becomes P f x f x q q i + q i 5 = q i q i q q i+q, where the secod iequality is by q i i q ad q > 0 for q = Olog / 3. By applyig Eq. 4 ad Eq. 5 to E + ad E, respectively, Eq. chages to E[[V ξ t V ξ t+ ξ t = x]] + c i 4b log c i c 6 qi+ 3 b log That is, the coditio of Lemma.5 still holds with. Thus, 6 the expected ruig time is polyomial. THEOREM 4.5. For the +-EA o LeadigOes uder bit-wise oise, q, if q = ωlog / o/, the expected ruig time is super-polyomial. PROOF. We use the same aalysis process as Theorem 4.. The oly differece is the probability of acceptig the offsprig solutio due to the chage of oise. For the positive drift, we still have E + i, sice we assume that x is always accepted i the proof of Theorem 4.. For the egative drift, we eed to re-aalyze P f x f x for the paret solutio x with LOx = k ad the offsprig solutio x with LOx = j where j k +. For j k, to derive a lower boud o P f x f x, we cosider the j cases where f x = l ad f x l for 0 l j. Sice P f x = l = q l q ad P f x l = q l, Eq. 4 chages to j P f x f x q l q q l qj l=0 = qj + q q 6 j q j qj q qj, where the last iequality is by q j qj / sice q = o/. For j = k + i.e., LOx = LOx = k, we ca use the same aalysis as Eq. 4 to derive a lower boud /, sice the last iequality of Eq. 4 still holds with p = o/. Thus, Eq. 5 also holds here, i.e., P f x f x. 7 By applyig Eqs. 6 ad 7 to E, Eq. 6 chages to E qk + i k. 6 Thus, we have E[[X t X t+ X t = i]] = E + E i qk i k. 6 For the upper boud aalysis of P X t+ i X t = i i the proof of Theorem 4., we oly eed to replace the acceptace probability p k+ i the case of LOx < LOx with k + q i.e., Eq. 5. Thus, Eq. 7 chages to P X t+ i X t = i k + q + k q + k. To compare E[[X t X t+ X t = i]] with P X t+ i X t = i, we cosider two cases: k < q ad k q. By usig q = ωlog / ad applyig the same aalysis procedure as Eqs. 8 ad 9, we ca derive that coditio of Lemma.7 holds with ɛ =. 9 For the lower boud aalysis of P X t+ i X t = i, by applyig Eqs. 6 ad 7, Eq. 0 chages to qkk + P X t+ i X t = i + k 6. For the aalysis of X t+ X t j, by replacig the acceptace probability p k+ i the case of LOx < LOx with k + q, Eq. chages to P X t+ X t j X t = i qkk+ 4 + k j j qkk + + k j That is, coditio of Lemma.7 holds with δ =, rl = 48. Thus, the expected ruig time is super-polyomial. THEOREM 4.6. For the +-EA o LeadigOes uder bitwise oise, q, the expected ruig time is expoetial if q = Ω/. PROOF. We use Lemma.6 to prove it. Let X t = i be the umber of 0-bits of the solutio x after t iteratios of the +-EA. We cosider the iterval i [0, / ]. To aalyze the drift E[[X t X t+ X t = i]], we use the same aalysis procedure as the proof of Theorem 4.. We first cosider q = Ω/ o. We eed to aalyze the probability P f x f x, where the offsprig solutio x is geerated by flippig oly oe -bit of x. Let LOx = k. For the case where the j-th where j k leadig -bit is flipped, as the aalysis of Eq. 6, we get P f x f x If q j <, qj 4 qj q j qj q. ; otherwise, qj qj q qj. Thus, we have P f x f x mi{/4, qj/}. For the case that flips oe o-leadig -bit i.e., LOx = LOx = k, to derive a lower boud o P f x f x, we cosider f x = l ad f x l for 0 l k. Thus, k P f x f x q l q q l + q k+ q k qk l=0 + q k+ = + qk q,

8 GECCO 7, July 5-9, 07, Berli, Germay where the last iequality is by q = o. By applyig the above two iequalities to Eq. 6, we get k E { mi e 4, qj } + i k. j= If k, k j= mi{, qj } = Ω sice q = Ω/. If 4 k <, i k = Ω sice i. Thus, E = Ω. For q = Ω, we use the trivial lower boud q for the probability of acceptig the offsprig solutio x, sice it is sufficiet to flip the first leadig -bit of x by oise. The, E iq kq + i kq = = Ω. e e Thus, for q = Ω/, we have E[[X t X t+ X t = i]] = E + E i Ω = Ω. That is, coditio of Lemma.6 holds. Its coditio trivially holds with δ = ad rl =. Thus, the expected ruig time is expoetial. 4.3 Oe-bit Noise For the +-EA o LeadigOes uder oe-bit oise, it has bee kow that the ruig time is polyomial if p 6e ad expoetial if p = []. We exted this result by provig i Theorem 4.7 that the ruig time is polyomial if p = Olog / ad super-polyomial if p = ωlog /. The proof ca be accomplished as same as that of Theorems 4., 4. ad 4.3 for bit-wise oise p,. This is because although the probabilities P f x f x of acceptig the offsprig solutio are differet, their bouds used i the proofs for bit-wise oise p, still hold for oe-bit oise. THEOREM 4.7. For the +-EA o LeadigOes uder oe-bit oise, the expected ruig time is polyomial if p = Olog /, super-polyomial if p = ωlog / o ad expoetial if p = Ω. PROOF. We re-aalyze P f x f x for oe-bit oise, ad show that the bouds o P f x f x used i the proofs for bit-wise oise p, still hold for oe-bit oise. For the proof of Theorem 4., Eqs. ad 3 chage to P f x i = p i, P f x i + = p, ad thus Eq. 4 still holds; Eqs. 7 ad 8 chage to P f x i p, P f x i = p i, ad thus Eq. 9 still holds. For the proof of Theorem 4., Eqs. ad 3 chage to P f x j = p j, P f x j = p j, ad thus Eq. 4 still holds. For the proof of Theorem 4.3, Eq. still holds by the above two equalities; Eq. 3 still holds sice the probability of keepig the first k + bits of a solutio uchaged i oe-bit oise is p k+ p k+. 5 CONCLUSION C. Qia et. al. I this paper, we theoretically study the +-EA solvig OeMax ad LeadigOes uder bit-wise oise. We derive the rages of oise parameters for the ruig time beig polyomial ad super-polyomial, respectively. The previously kow rages for the +-EA solvig LeadigOes uder oe-bit oise are also improved. I the future, we shall improve the curretly derived bouds o LeadigOes, as they do ot cover the whole rage of oise parameters. ACKNOWLEDGMENTS We wat to thak the reviewers for their valuable commets. This work was supported by the NSFC , , the Fudametal Research Fuds for the Cetral Uiversities ad the CCF-Tecet Ope Research Fud. REFERENCES [] Y. Akimoto, S. Astete-Morales, ad O. Teytaud. 05. Aalysis of rutime of optimizatio algorithms for oisy fuctios over discrete codomais. Theoretical Computer Sciece , [] A. Auger ad B. Doerr. 0. Theory of Radomized Search Heuristics: Foudatios ad Recet Developmets. World Scietific, Sigapore. [3] L. Biachi, M. Dorigo, L. Gambardella, ad W. Gutjahr A survey o metaheuristics for stochastic combiatorial optimizatio. Natural Computig 8, 009, [4] D.-C. Dag ad P. K. Lehre. 05. Efficiet optimisatio of oisy fitess fuctios with populatio-based evolutioary algorithms. I Proceedigs of FOGA 5. Aberystwyth, UK, [5] B. Doerr, A. Hota, ad T. Kötzig. 0. Ats easily solve stochastic shortest path problems. I Proceedigs of GECCO. Philadelphia, PA, 7 4. [6] S. Droste Aalysis of the + EA for a oisy OeMax. I Proceedigs of GECCO 04. Seattle, WA, [7] S. Droste, T. Jase, ad I. Wegeer. 00. O the aalysis of the + evolutioary algorithm. Theoretical Computer Sciece 76, - 00, 5 8. [8] T. Friedrich, T. Kötzig, M. Krejca, ad A. Sutto. 06. The compact geetic algorithm is efficiet uder extreme Gaussia oise. IEEE Trasactios o Evolutioary Computatio 06, i press. [9] T. Friedrich, T. Kötzig, M. Krejca, ad A. Sutto. 06. Robustess of at coloy optimizatio to oise. Evolutioary Computatio 4, 06, [0] T. Friedrich, T. Kötzig, F. Quiza, ad A. Sutto. 07. Resamplig vs recombiatio: A statistical ru time estimatio. I Proceedigs of FOGA 7. Copehage, Demark. [] C. Gieße ad T. Kötzig. 06. Robustess of populatios i stochastic eviromets. Algorithmica 75, 3 06, [] J. He ad X. Yao. 00. Drift aalysis ad average time complexity of evolutioary algorithms. Artificial Itelligece 7, 00, [3] Y. Ji ad J. Brake Evolutioary optimizatio i ucertai eviromets-a survey. IEEE Trasactios o Evolutioary Computatio 9, 3 005, [4] F. Neuma ad C. Witt. 00. Bioispired Computatio i Combiatorial Optimizatio: Algorithms ad Their Computatioal Complexity. Spriger-Verlag, Berli, Germay. [5] P. Oliveto ad C. Witt. 0. Simplified drift aalysis for provig lower bouds i evolutioary computatio. Algorithmica 59, 3 0, [6] P. Oliveto ad C. Witt. 0. Erratum: Simplified drift aalysis for provig lower bouds i evolutioary computatio. CORR abs/ [7] A. Prugel-Beett, J. Rowe, ad J. Shapiro. 05. Ru-time aalysis of populatio-based evolutioary algorithm i oisy eviromets. I Proceedigs of FOGA 5. Aberystwyth, UK, [8] C. Qia, Y. Yu, K. Tag, Y. Ji, X. Yao, ad Z.-H. Zhou. 06. O the effectiveess of samplig for evolutioary optimizatio i oisy eviromets. Evolutioary Computatio 06, i press. [9] C. Qia, Y. Yu, ad Z.-H. Zhou. 05. Aalyzig evolutioary optimizatio i oisy eviromets. Evolutioary Computatio 05, i press. [0] J. E. Rowe ad D. Sudholt. 04. The choice of the offsprig populatio size i the,λ evolutioary algorithm. Theoretical Computer Sciece , 0 38.

arxiv: v1 [cs.ne] 2 Nov 2017

arxiv: v1 [cs.ne] 2 Nov 2017 Noame mauscript No. will be iserted by the editor) Ruig Time Aalysis of the +)-EA for OeMax ad LeadigOes uder Bit-wise Noise Chao Qia Chao Bia Wu Jiag Ke Tag Received: date / Accepted: date arxiv:7.00956v

More information

On the Effectiveness of Sampling for Evolutionary Optimization in Noisy Environments

On the Effectiveness of Sampling for Evolutionary Optimization in Noisy Environments O the Effectiveess of Samplig for Evolutioary Optimizatio i Noisy Eviromets Chao Qia,2 chaoqia@ustc.edu.c Yag Yu 2 yuy@ju.edu.c Ke Tag ketag@ustc.edu.c Yaochu Ji 3 yaochu.ji@surrey.ac.uk Xi Yao,4 x.yao@cs.bham.ac.uk

More information

arxiv: v1 [cs.ai] 20 Nov 2013

arxiv: v1 [cs.ai] 20 Nov 2013 Aalyzig Evolutioary Optimizatio i Noisy Eviromets Chao Qia, Yag Yu, Zhi-Hua Zhou Natioal Key Laboratory for Novel Software Techology Najig Uiversity, Najig 20023, Chia arxiv:3.4987v [cs.ai] 20 Nov 203

More information

General Lower Bounds for the Running Time of Evolutionary Algorithms

General Lower Bounds for the Running Time of Evolutionary Algorithms Geeral Lower Bouds for the Ruig Time of Evolutioary Algorithms Dirk Sudholt Iteratioal Computer Sciece Istitute, Berkeley, CA 94704, USA Abstract. We preset a ew method for provig lower bouds i evolutioary

More information

Analyzing Evolutionary Optimization in Noisy Environments

Analyzing Evolutionary Optimization in Noisy Environments Evolutioary Computatio /EVCO_a_0070-Qia Jauary 5, 206 6:7 Aalyzig Evolutioary Optimizatio i Noisy Eviromets Chao Qia qiac@lamda.ju.edu.c Natioal Key Laboratory for Novel Software Techology, Najig Uiversity,

More information

Selection Hyper-heuristics Can Provably be Helpful in Evolutionary Multi-objective Optimization

Selection Hyper-heuristics Can Provably be Helpful in Evolutionary Multi-objective Optimization Selectio Hyper-heuristics Ca Provably be Helpful i Evolutioary Multi-objective Optimizatio Chao Qia 1,2, Ke Tag 1, ad Zhi-Hua Zhou 2 1 UBRI, School of Computer Sciece ad Techology, Uiversity of Sciece

More information

A Rigorous View On Neutrality

A Rigorous View On Neutrality A Rigorous View O Neutrality Bejami Doerr Michael Gewuch Nils Hebbighaus Frak Neuma Algorithms ad Complexity Group Max-Plack-Istitut für Iformatik Saarbrücke, Germay Departmet of Computer Sciece Christia-Albrechts-Uiversity

More information

Reinforcement Learning Based Dynamic Selection of Auxiliary Objectives with Preserving of the Best Found Solution

Reinforcement Learning Based Dynamic Selection of Auxiliary Objectives with Preserving of the Best Found Solution Reiforcemet Learig Based Dyamic Selectio of Auxiliary Objectives with Preservig of the Best Foud Solutio arxiv:1704.07187v1 [cs.ne] 24 Apr 2017 Abstract Efficiecy of sigle-objective optimizatio ca be improved

More information

An Introduction to Randomized Algorithms

An Introduction to Randomized Algorithms A Itroductio to Radomized Algorithms The focus of this lecture is to study a radomized algorithm for quick sort, aalyze it usig probabilistic recurrece relatios, ad also provide more geeral tools for aalysis

More information

arxiv: v1 [cs.ne] 4 Sep 2017

arxiv: v1 [cs.ne] 4 Sep 2017 Theoretical Aalysis of Stochastic Search Algorithms Per Kristia Lehre School of Computer Sciece, Uiversity of Birmigham, Birmigham, UK Pietro S. Oliveto Departmet of Computer Sciece, Uiversity of Sheffield,

More information

Problem Set 2 Solutions

Problem Set 2 Solutions CS271 Radomess & Computatio, Sprig 2018 Problem Set 2 Solutios Poit totals are i the margi; the maximum total umber of poits was 52. 1. Probabilistic method for domiatig sets 6pts Pick a radom subset S

More information

On a Smarandache problem concerning the prime gaps

On a Smarandache problem concerning the prime gaps O a Smaradache problem cocerig the prime gaps Felice Russo Via A. Ifate 7 6705 Avezzao (Aq) Italy felice.russo@katamail.com Abstract I this paper, a problem posed i [] by Smaradache cocerig the prime gaps

More information

Upper and Lower Bounds on Unrestricted Black-Box Complexity of Jump n,l

Upper and Lower Bounds on Unrestricted Black-Box Complexity of Jump n,l Upper ad Lower Bouds o Urestricted Black-Box Complexity of Jump,l Maxim Buzdalov 1, Mikhail Kever 1, ad Bejami Doerr 1 ITMO Uiversity, 49 Kroverkskiy av., Sait-Petersburg, Russia, 197101 mbuzdalov@gmail.com,

More information

Drift analysis and average time complexity of evolutionary algorithms

Drift analysis and average time complexity of evolutionary algorithms Artificial Itelligece 127 (2001) 57 85 Drift aalysis ad average time complexity of evolutioary algorithms Ju He a,xiyao b, a Departmet of Computer Sciece, Norther Jiaotog Uiversity, Beijig 100044, PR Chia

More information

Application to Random Graphs

Application to Random Graphs A Applicatio to Radom Graphs Brachig processes have a umber of iterestig ad importat applicatios. We shall cosider oe of the most famous of them, the Erdős-Réyi radom graph theory. 1 Defiitio A.1. Let

More information

A collocation method for singular integral equations with cosecant kernel via Semi-trigonometric interpolation

A collocation method for singular integral equations with cosecant kernel via Semi-trigonometric interpolation Iteratioal Joural of Mathematics Research. ISSN 0976-5840 Volume 9 Number 1 (017) pp. 45-51 Iteratioal Research Publicatio House http://www.irphouse.com A collocatio method for sigular itegral equatios

More information

A Novel Genetic Algorithm using Helper Objectives for the 0-1 Knapsack Problem

A Novel Genetic Algorithm using Helper Objectives for the 0-1 Knapsack Problem A Novel Geetic Algorithm usig Helper Objectives for the 0-1 Kapsack Problem Ju He, Feidu He ad Hogbi Dog 1 arxiv:1404.0868v1 [cs.ne] 3 Apr 2014 Abstract The 0-1 kapsack problem is a well-kow combiatorial

More information

Lecture 11: Pseudorandom functions

Lecture 11: Pseudorandom functions COM S 6830 Cryptography Oct 1, 2009 Istructor: Rafael Pass 1 Recap Lecture 11: Pseudoradom fuctios Scribe: Stefao Ermo Defiitio 1 (Ge, Ec, Dec) is a sigle message secure ecryptio scheme if for all uppt

More information

Runtime Analysis of the (µ+1) EA on Simple Pseudo-Boolean Functions

Runtime Analysis of the (µ+1) EA on Simple Pseudo-Boolean Functions Rutime Aalysis of the (µ+1) EA o Simple Pseudo-Boolea Fuctios Carste Witt FB Iformatik, LS 2, Uiversität Dortmud, 44221 Dortmud, Germay carste.witt@cs.ui-dortmud.de November 14, 2005 Abstract Although

More information

Element sampling: Part 2

Element sampling: Part 2 Chapter 4 Elemet samplig: Part 2 4.1 Itroductio We ow cosider uequal probability samplig desigs which is very popular i practice. I the uequal probability samplig, we ca improve the efficiecy of the resultig

More information

Random Variables, Sampling and Estimation

Random Variables, Sampling and Estimation Chapter 1 Radom Variables, Samplig ad Estimatio 1.1 Itroductio This chapter will cover the most importat basic statistical theory you eed i order to uderstad the ecoometric material that will be comig

More information

Lecture 3: August 31

Lecture 3: August 31 36-705: Itermediate Statistics Fall 018 Lecturer: Siva Balakrisha Lecture 3: August 31 This lecture will be mostly a summary of other useful expoetial tail bouds We will ot prove ay of these i lecture,

More information

Problem Set 4 Due Oct, 12

Problem Set 4 Due Oct, 12 EE226: Radom Processes i Systems Lecturer: Jea C. Walrad Problem Set 4 Due Oct, 12 Fall 06 GSI: Assae Gueye This problem set essetially reviews detectio theory ad hypothesis testig ad some basic otios

More information

The Growth of Functions. Theoretical Supplement

The Growth of Functions. Theoretical Supplement The Growth of Fuctios Theoretical Supplemet The Triagle Iequality The triagle iequality is a algebraic tool that is ofte useful i maipulatig absolute values of fuctios. The triagle iequality says that

More information

Optimally Sparse SVMs

Optimally Sparse SVMs A. Proof of Lemma 3. We here prove a lower boud o the umber of support vectors to achieve geeralizatio bouds of the form which we cosider. Importatly, this result holds ot oly for liear classifiers, but

More information

Advanced Stochastic Processes.

Advanced Stochastic Processes. Advaced Stochastic Processes. David Gamarik LECTURE 2 Radom variables ad measurable fuctios. Strog Law of Large Numbers (SLLN). Scary stuff cotiued... Outlie of Lecture Radom variables ad measurable fuctios.

More information

Populations can be essential in tracking dynamic optima

Populations can be essential in tracking dynamic optima Populatios ca be essetial i trackig dyamic optima Duc-Cuog Dag 1, Thomas Jase 2, ad Per Kristia Lehre 1 arxiv:1607.03317v1 [cs.ne] 12 Jul 2016 1 School of Computer Sciece, Uiversity of Nottigham Jubilee

More information

Posted-Price, Sealed-Bid Auctions

Posted-Price, Sealed-Bid Auctions Posted-Price, Sealed-Bid Auctios Professors Greewald ad Oyakawa 207-02-08 We itroduce the posted-price, sealed-bid auctio. This auctio format itroduces the idea of approximatios. We describe how well this

More information

Let us give one more example of MLE. Example 3. The uniform distribution U[0, θ] on the interval [0, θ] has p.d.f.

Let us give one more example of MLE. Example 3. The uniform distribution U[0, θ] on the interval [0, θ] has p.d.f. Lecture 5 Let us give oe more example of MLE. Example 3. The uiform distributio U[0, ] o the iterval [0, ] has p.d.f. { 1 f(x =, 0 x, 0, otherwise The likelihood fuctio ϕ( = f(x i = 1 I(X 1,..., X [0,

More information

ACO Comprehensive Exam 9 October 2007 Student code A. 1. Graph Theory

ACO Comprehensive Exam 9 October 2007 Student code A. 1. Graph Theory 1. Graph Theory Prove that there exist o simple plaar triagulatio T ad two distict adjacet vertices x, y V (T ) such that x ad y are the oly vertices of T of odd degree. Do ot use the Four-Color Theorem.

More information

Lecture 2: Monte Carlo Simulation

Lecture 2: Monte Carlo Simulation STAT/Q SCI 43: Itroductio to Resamplig ethods Sprig 27 Istructor: Ye-Chi Che Lecture 2: ote Carlo Simulatio 2 ote Carlo Itegratio Assume we wat to evaluate the followig itegratio: e x3 dx What ca we do?

More information

Lecture 2 February 8, 2016

Lecture 2 February 8, 2016 MIT 6.854/8.45: Advaced Algorithms Sprig 206 Prof. Akur Moitra Lecture 2 February 8, 206 Scribe: Calvi Huag, Lih V. Nguye I this lecture, we aalyze the problem of schedulig equal size tasks arrivig olie

More information

Randomized Algorithms I, Spring 2018, Department of Computer Science, University of Helsinki Homework 1: Solutions (Discussed January 25, 2018)

Randomized Algorithms I, Spring 2018, Department of Computer Science, University of Helsinki Homework 1: Solutions (Discussed January 25, 2018) Radomized Algorithms I, Sprig 08, Departmet of Computer Sciece, Uiversity of Helsiki Homework : Solutios Discussed Jauary 5, 08). Exercise.: Cosider the followig balls-ad-bi game. We start with oe black

More information

Chapter 2 The Monte Carlo Method

Chapter 2 The Monte Carlo Method Chapter 2 The Mote Carlo Method The Mote Carlo Method stads for a broad class of computatioal algorithms that rely o radom sampligs. It is ofte used i physical ad mathematical problems ad is most useful

More information

Basics of Probability Theory (for Theory of Computation courses)

Basics of Probability Theory (for Theory of Computation courses) Basics of Probability Theory (for Theory of Computatio courses) Oded Goldreich Departmet of Computer Sciece Weizma Istitute of Sciece Rehovot, Israel. oded.goldreich@weizma.ac.il November 24, 2008 Preface.

More information

ON POINTWISE BINOMIAL APPROXIMATION

ON POINTWISE BINOMIAL APPROXIMATION Iteratioal Joural of Pure ad Applied Mathematics Volume 71 No. 1 2011, 57-66 ON POINTWISE BINOMIAL APPROXIMATION BY w-functions K. Teerapabolar 1, P. Wogkasem 2 Departmet of Mathematics Faculty of Sciece

More information

Monte Carlo Integration

Monte Carlo Integration Mote Carlo Itegratio I these otes we first review basic umerical itegratio methods (usig Riema approximatio ad the trapezoidal rule) ad their limitatios for evaluatig multidimesioal itegrals. Next we itroduce

More information

DETERMINATION OF MECHANICAL PROPERTIES OF A NON- UNIFORM BEAM USING THE MEASUREMENT OF THE EXCITED LONGITUDINAL ELASTIC VIBRATIONS.

DETERMINATION OF MECHANICAL PROPERTIES OF A NON- UNIFORM BEAM USING THE MEASUREMENT OF THE EXCITED LONGITUDINAL ELASTIC VIBRATIONS. ICSV4 Cairs Australia 9- July 7 DTRMINATION OF MCHANICAL PROPRTIS OF A NON- UNIFORM BAM USING TH MASURMNT OF TH XCITD LONGITUDINAL LASTIC VIBRATIONS Pavel Aokhi ad Vladimir Gordo Departmet of the mathematics

More information

arxiv: v1 [cs.ne] 26 Mar 2014

arxiv: v1 [cs.ne] 26 Mar 2014 How Crossover Speeds Up Buildig-Block Assembly i Geetic Algorithms arxiv:403.6600v [cs.ne] 26 Mar 204 Dirk Sudholt Departmet of Computer Sciece, Uiversity of Sheffield, Uited Kigdom Abstract We re-ivestigate

More information

Expectation and Variance of a random variable

Expectation and Variance of a random variable Chapter 11 Expectatio ad Variace of a radom variable The aim of this lecture is to defie ad itroduce mathematical Expectatio ad variace of a fuctio of discrete & cotiuous radom variables ad the distributio

More information

arxiv: v2 [cs.ne] 26 Nov 2014

arxiv: v2 [cs.ne] 26 Nov 2014 How Crossover Speeds Up Buildig-Block Assembly i Geetic Algorithms arxiv:403.6600v2 [cs.ne] 26 Nov 204 Dirk Sudholt Departmet of Computer Sciece, Uiversity of Sheffield, Uited Kigdom Abstract We re-ivestigate

More information

IP Reference guide for integer programming formulations.

IP Reference guide for integer programming formulations. IP Referece guide for iteger programmig formulatios. by James B. Orli for 15.053 ad 15.058 This documet is iteded as a compact (or relatively compact) guide to the formulatio of iteger programs. For more

More information

Rank Modulation with Multiplicity

Rank Modulation with Multiplicity Rak Modulatio with Multiplicity Axiao (Adrew) Jiag Computer Sciece ad Eg. Dept. Texas A&M Uiversity College Statio, TX 778 ajiag@cse.tamu.edu Abstract Rak modulatio is a scheme that uses the relative order

More information

Computational Intelligence Winter Term 2018/19

Computational Intelligence Winter Term 2018/19 Computatioal Itelligece Witer Term 28/9 Prof. Dr. Güter Rudolph Lehrstuhl für Algorithm Egieerig (LS ) Fakultät für Iformatik TU Dortmud Pla for Today Lecture Evolutioary Algorithms (EA) Optimizatio Basics

More information

Since X n /n P p, we know that X n (n. Xn (n X n ) Using the asymptotic result above to obtain an approximation for fixed n, we obtain

Since X n /n P p, we know that X n (n. Xn (n X n ) Using the asymptotic result above to obtain an approximation for fixed n, we obtain Assigmet 9 Exercise 5.5 Let X biomial, p, where p 0, 1 is ukow. Obtai cofidece itervals for p i two differet ways: a Sice X / p d N0, p1 p], the variace of the limitig distributio depeds oly o p. Use the

More information

Stat 421-SP2012 Interval Estimation Section

Stat 421-SP2012 Interval Estimation Section Stat 41-SP01 Iterval Estimatio Sectio 11.1-11. We ow uderstad (Chapter 10) how to fid poit estimators of a ukow parameter. o However, a poit estimate does ot provide ay iformatio about the ucertaity (possible

More information

MATH 320: Probability and Statistics 9. Estimation and Testing of Parameters. Readings: Pruim, Chapter 4

MATH 320: Probability and Statistics 9. Estimation and Testing of Parameters. Readings: Pruim, Chapter 4 MATH 30: Probability ad Statistics 9. Estimatio ad Testig of Parameters Estimatio ad Testig of Parameters We have bee dealig situatios i which we have full kowledge of the distributio of a radom variable.

More information

Vector Quantization: a Limiting Case of EM

Vector Quantization: a Limiting Case of EM . Itroductio & defiitios Assume that you are give a data set X = { x j }, j { 2,,, }, of d -dimesioal vectors. The vector quatizatio (VQ) problem requires that we fid a set of prototype vectors Z = { z

More information

FACULTY OF MATHEMATICAL STUDIES MATHEMATICS FOR PART I ENGINEERING. Lectures

FACULTY OF MATHEMATICAL STUDIES MATHEMATICS FOR PART I ENGINEERING. Lectures FACULTY OF MATHEMATICAL STUDIES MATHEMATICS FOR PART I ENGINEERING Lectures MODULE 5 STATISTICS II. Mea ad stadard error of sample data. Biomial distributio. Normal distributio 4. Samplig 5. Cofidece itervals

More information

Design and Analysis of Algorithms

Design and Analysis of Algorithms Desig ad Aalysis of Algorithms Probabilistic aalysis ad Radomized algorithms Referece: CLRS Chapter 5 Topics: Hirig problem Idicatio radom variables Radomized algorithms Huo Hogwei 1 The hirig problem

More information

Topic 9: Sampling Distributions of Estimators

Topic 9: Sampling Distributions of Estimators Topic 9: Samplig Distributios of Estimators Course 003, 2018 Page 0 Samplig distributios of estimators Sice our estimators are statistics (particular fuctios of radom variables), their distributio ca be

More information

Average Convergence Rate of Evolutionary Algorithms

Average Convergence Rate of Evolutionary Algorithms 1 Average Covergece Rate of Evolutioary Algorithms Ju He ad Guagmig Li Abstract arxiv:154.8117v3 [cs.ne] 2 Ju 215 I evolutioary optimizatio, it is importat to uderstad how fast evolutioary algorithms coverge

More information

Hypothesis Testing. Evaluation of Performance of Learned h. Issues. Trade-off Between Bias and Variance

Hypothesis Testing. Evaluation of Performance of Learned h. Issues. Trade-off Between Bias and Variance Hypothesis Testig Empirically evaluatig accuracy of hypotheses: importat activity i ML. Three questios: Give observed accuracy over a sample set, how well does this estimate apply over additioal samples?

More information

6.3 Testing Series With Positive Terms

6.3 Testing Series With Positive Terms 6.3. TESTING SERIES WITH POSITIVE TERMS 307 6.3 Testig Series With Positive Terms 6.3. Review of what is kow up to ow I theory, testig a series a i for covergece amouts to fidig the i= sequece of partial

More information

Lecture 7: October 18, 2017

Lecture 7: October 18, 2017 Iformatio ad Codig Theory Autum 207 Lecturer: Madhur Tulsiai Lecture 7: October 8, 207 Biary hypothesis testig I this lecture, we apply the tools developed i the past few lectures to uderstad the problem

More information

Lecture 14: Graph Entropy

Lecture 14: Graph Entropy 15-859: Iformatio Theory ad Applicatios i TCS Sprig 2013 Lecture 14: Graph Etropy March 19, 2013 Lecturer: Mahdi Cheraghchi Scribe: Euiwoog Lee 1 Recap Bergma s boud o the permaet Shearer s Lemma Number

More information

CS / MCS 401 Homework 3 grader solutions

CS / MCS 401 Homework 3 grader solutions CS / MCS 401 Homework 3 grader solutios assigmet due July 6, 016 writte by Jāis Lazovskis maximum poits: 33 Some questios from CLRS. Questios marked with a asterisk were ot graded. 1 Use the defiitio of

More information

w (1) ˆx w (1) x (1) /ρ and w (2) ˆx w (2) x (2) /ρ.

w (1) ˆx w (1) x (1) /ρ and w (2) ˆx w (2) x (2) /ρ. 2 5. Weighted umber of late jobs 5.1. Release dates ad due dates: maximimizig the weight of o-time jobs Oce we add release dates, miimizig the umber of late jobs becomes a sigificatly harder problem. For

More information

Econ 325/327 Notes on Sample Mean, Sample Proportion, Central Limit Theorem, Chi-square Distribution, Student s t distribution 1.

Econ 325/327 Notes on Sample Mean, Sample Proportion, Central Limit Theorem, Chi-square Distribution, Student s t distribution 1. Eco 325/327 Notes o Sample Mea, Sample Proportio, Cetral Limit Theorem, Chi-square Distributio, Studet s t distributio 1 Sample Mea By Hiro Kasahara We cosider a radom sample from a populatio. Defiitio

More information

Escaping Local Optima Using Crossover with Emergent Diversity

Escaping Local Optima Using Crossover with Emergent Diversity This article has bee accepted for publicatio i a future issue of this joural, but has ot bee fully edited Cotet may chage prior to fial publicatio Citatio iformatio: DOI 009/TEVC07740, IEEE Trasactios

More information

Topic 9: Sampling Distributions of Estimators

Topic 9: Sampling Distributions of Estimators Topic 9: Samplig Distributios of Estimators Course 003, 2018 Page 0 Samplig distributios of estimators Sice our estimators are statistics (particular fuctios of radom variables), their distributio ca be

More information

The coalescent coalescence theory

The coalescent coalescence theory The coalescet coalescece theory Peter Beerli September 1, 009 Historical ote Up to 198 most developmet i populatio geetics was prospective ad developed expectatios based o situatios of today. Most work

More information

Resampling Methods. X (1/2), i.e., Pr (X i m) = 1/2. We order the data: X (1) X (2) X (n). Define the sample median: ( n.

Resampling Methods. X (1/2), i.e., Pr (X i m) = 1/2. We order the data: X (1) X (2) X (n). Define the sample median: ( n. Jauary 1, 2019 Resamplig Methods Motivatio We have so may estimators with the property θ θ d N 0, σ 2 We ca also write θ a N θ, σ 2 /, where a meas approximately distributed as Oce we have a cosistet estimator

More information

Lecture 7: Properties of Random Samples

Lecture 7: Properties of Random Samples Lecture 7: Properties of Radom Samples 1 Cotiued From Last Class Theorem 1.1. Let X 1, X,...X be a radom sample from a populatio with mea µ ad variace σ

More information

The log-behavior of n p(n) and n p(n)/n

The log-behavior of n p(n) and n p(n)/n Ramauja J. 44 017, 81-99 The log-behavior of p ad p/ William Y.C. Che 1 ad Ke Y. Zheg 1 Ceter for Applied Mathematics Tiaji Uiversity Tiaji 0007, P. R. Chia Ceter for Combiatorics, LPMC Nakai Uivercity

More information

On forward improvement iteration for stopping problems

On forward improvement iteration for stopping problems O forward improvemet iteratio for stoppig problems Mathematical Istitute, Uiversity of Kiel, Ludewig-Mey-Str. 4, D-24098 Kiel, Germay irle@math.ui-iel.de Albrecht Irle Abstract. We cosider the optimal

More information

Last time, we talked about how Equation (1) can simulate Equation (2). We asserted that Equation (2) can also simulate Equation (1).

Last time, we talked about how Equation (1) can simulate Equation (2). We asserted that Equation (2) can also simulate Equation (1). 6896 Quatum Complexity Theory Sept 23, 2008 Lecturer: Scott Aaroso Lecture 6 Last Time: Quatum Error-Correctio Quatum Query Model Deutsch-Jozsa Algorithm (Computes x y i oe query) Today: Berstei-Vazirii

More information

1 Inferential Methods for Correlation and Regression Analysis

1 Inferential Methods for Correlation and Regression Analysis 1 Iferetial Methods for Correlatio ad Regressio Aalysis I the chapter o Correlatio ad Regressio Aalysis tools for describig bivariate cotiuous data were itroduced. The sample Pearso Correlatio Coefficiet

More information

Hashing and Amortization

Hashing and Amortization Lecture Hashig ad Amortizatio Supplemetal readig i CLRS: Chapter ; Chapter 7 itro; Sectio 7.. Arrays ad Hashig Arrays are very useful. The items i a array are statically addressed, so that isertig, deletig,

More information

ECE 901 Lecture 12: Complexity Regularization and the Squared Loss

ECE 901 Lecture 12: Complexity Regularization and the Squared Loss ECE 90 Lecture : Complexity Regularizatio ad the Squared Loss R. Nowak 5/7/009 I the previous lectures we made use of the Cheroff/Hoeffdig bouds for our aalysis of classifier errors. Hoeffdig s iequality

More information

Econ 325 Notes on Point Estimator and Confidence Interval 1 By Hiro Kasahara

Econ 325 Notes on Point Estimator and Confidence Interval 1 By Hiro Kasahara Poit Estimator Eco 325 Notes o Poit Estimator ad Cofidece Iterval 1 By Hiro Kasahara Parameter, Estimator, ad Estimate The ormal probability desity fuctio is fully characterized by two costats: populatio

More information

CS284A: Representations and Algorithms in Molecular Biology

CS284A: Representations and Algorithms in Molecular Biology CS284A: Represetatios ad Algorithms i Molecular Biology Scribe Notes o Lectures 3 & 4: Motif Discovery via Eumeratio & Motif Represetatio Usig Positio Weight Matrix Joshua Gervi Based o presetatios by

More information

Lecture 5: April 17, 2013

Lecture 5: April 17, 2013 TTIC/CMSC 350 Mathematical Toolkit Sprig 203 Madhur Tulsiai Lecture 5: April 7, 203 Scribe: Somaye Hashemifar Cheroff bouds recap We recall the Cheroff/Hoeffdig bouds we derived i the last lecture idepedet

More information

Topic 9: Sampling Distributions of Estimators

Topic 9: Sampling Distributions of Estimators Topic 9: Samplig Distributios of Estimators Course 003, 2016 Page 0 Samplig distributios of estimators Sice our estimators are statistics (particular fuctios of radom variables), their distributio ca be

More information

Similarity Solutions to Unsteady Pseudoplastic. Flow Near a Moving Wall

Similarity Solutions to Unsteady Pseudoplastic. Flow Near a Moving Wall Iteratioal Mathematical Forum, Vol. 9, 04, o. 3, 465-475 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/0.988/imf.04.48 Similarity Solutios to Usteady Pseudoplastic Flow Near a Movig Wall W. Robi Egieerig

More information

EECS564 Estimation, Filtering, and Detection Hwk 2 Solns. Winter p θ (z) = (2θz + 1 θ), 0 z 1

EECS564 Estimation, Filtering, and Detection Hwk 2 Solns. Winter p θ (z) = (2θz + 1 θ), 0 z 1 EECS564 Estimatio, Filterig, ad Detectio Hwk 2 Sols. Witer 25 4. Let Z be a sigle observatio havig desity fuctio where. p (z) = (2z + ), z (a) Assumig that is a oradom parameter, fid ad plot the maximum

More information

CS 330 Discussion - Probability

CS 330 Discussion - Probability CS 330 Discussio - Probability March 24 2017 1 Fudametals of Probability 11 Radom Variables ad Evets A radom variable X is oe whose value is o-determiistic For example, suppose we flip a coi ad set X =

More information

-ORDER CONVERGENCE FOR FINDING SIMPLE ROOT OF A POLYNOMIAL EQUATION

-ORDER CONVERGENCE FOR FINDING SIMPLE ROOT OF A POLYNOMIAL EQUATION NEW NEWTON-TYPE METHOD WITH k -ORDER CONVERGENCE FOR FINDING SIMPLE ROOT OF A POLYNOMIAL EQUATION R. Thukral Padé Research Cetre, 39 Deaswood Hill, Leeds West Yorkshire, LS7 JS, ENGLAND ABSTRACT The objective

More information

(b) What is the probability that a particle reaches the upper boundary n before the lower boundary m?

(b) What is the probability that a particle reaches the upper boundary n before the lower boundary m? MATH 529 The Boudary Problem The drukard s walk (or boudary problem) is oe of the most famous problems i the theory of radom walks. Oe versio of the problem is described as follows: Suppose a particle

More information

Lecture 2. The Lovász Local Lemma

Lecture 2. The Lovász Local Lemma Staford Uiversity Sprig 208 Math 233A: No-costructive methods i combiatorics Istructor: Ja Vodrák Lecture date: Jauary 0, 208 Origial scribe: Apoorva Khare Lecture 2. The Lovász Local Lemma 2. Itroductio

More information

On comparison of different approaches to the stability radius calculation. Olga Karelkina

On comparison of different approaches to the stability radius calculation. Olga Karelkina O compariso of differet approaches to the stability radius calculatio Olga Karelkia Uiversity of Turku 2011 Outlie Prelimiaries Problem statemet Exact method for calculatio stability radius proposed by

More information

Machine Learning Brett Bernstein

Machine Learning Brett Bernstein Machie Learig Brett Berstei Week 2 Lecture: Cocept Check Exercises Starred problems are optioal. Excess Risk Decompositio 1. Let X = Y = {1, 2,..., 10}, A = {1,..., 10, 11} ad suppose the data distributio

More information

Lecture 12: November 13, 2018

Lecture 12: November 13, 2018 Mathematical Toolkit Autum 2018 Lecturer: Madhur Tulsiai Lecture 12: November 13, 2018 1 Radomized polyomial idetity testig We will use our kowledge of coditioal probability to prove the followig lemma,

More information

Random assignment with integer costs

Random assignment with integer costs Radom assigmet with iteger costs Robert Parviaie Departmet of Mathematics, Uppsala Uiversity P.O. Box 480, SE-7506 Uppsala, Swede robert.parviaie@math.uu.se Jue 4, 200 Abstract The radom assigmet problem

More information

Recursive Algorithm for Generating Partitions of an Integer. 1 Preliminary

Recursive Algorithm for Generating Partitions of an Integer. 1 Preliminary Recursive Algorithm for Geeratig Partitios of a Iteger Sug-Hyuk Cha Computer Sciece Departmet, Pace Uiversity 1 Pace Plaza, New York, NY 10038 USA scha@pace.edu Abstract. This article first reviews the

More information

Sequences and Limits

Sequences and Limits Chapter Sequeces ad Limits Let { a } be a sequece of real or complex umbers A ecessary ad sufficiet coditio for the sequece to coverge is that for ay ɛ > 0 there exists a iteger N > 0 such that a p a q

More information

Berry-Esseen bounds for self-normalized martingales

Berry-Esseen bounds for self-normalized martingales Berry-Essee bouds for self-ormalized martigales Xiequa Fa a, Qi-Ma Shao b a Ceter for Applied Mathematics, Tiaji Uiversity, Tiaji 30007, Chia b Departmet of Statistics, The Chiese Uiversity of Hog Kog,

More information

The Rand and block distances of pairs of set partitions

The Rand and block distances of pairs of set partitions The Rad ad block distaces of pairs of set partitios Frak Ruskey 1 ad Jeifer Woodcock 1 Dept. of Computer Sciece, Uiversity of Victoria, CANADA Abstract. The Rad distaces of two set partitios is the umber

More information

A Note on Effi cient Conditional Simulation of Gaussian Distributions. April 2010

A Note on Effi cient Conditional Simulation of Gaussian Distributions. April 2010 A Note o Effi ciet Coditioal Simulatio of Gaussia Distributios A D D C S S, U B C, V, BC, C April 2010 A Cosider a multivariate Gaussia radom vector which ca be partitioed ito observed ad uobserved compoetswe

More information

Lecture 8: Solving the Heat, Laplace and Wave equations using finite difference methods

Lecture 8: Solving the Heat, Laplace and Wave equations using finite difference methods Itroductory lecture otes o Partial Differetial Equatios - c Athoy Peirce. Not to be copied, used, or revised without explicit writte permissio from the copyright ower. 1 Lecture 8: Solvig the Heat, Laplace

More information

ECE 901 Lecture 14: Maximum Likelihood Estimation and Complexity Regularization

ECE 901 Lecture 14: Maximum Likelihood Estimation and Complexity Regularization ECE 90 Lecture 4: Maximum Likelihood Estimatio ad Complexity Regularizatio R Nowak 5/7/009 Review : Maximum Likelihood Estimatio We have iid observatios draw from a ukow distributio Y i iid p θ, i,, where

More information

1 Statement of the Game

1 Statement of the Game ANALYSIS OF THE CHOW-ROBBINS GAME JON LU May 10, 2016 Abstract Flip a coi repeatedly ad stop wheever you wat. Your payoff is the proportio of heads ad you wish to maximize this payoff i expectatio. I this

More information

Recurrence Relations

Recurrence Relations Recurrece Relatios Aalysis of recursive algorithms, such as: it factorial (it ) { if (==0) retur ; else retur ( * factorial(-)); } Let t be the umber of multiplicatios eeded to calculate factorial(). The

More information

7.1 Convergence of sequences of random variables

7.1 Convergence of sequences of random variables Chapter 7 Limit theorems Throughout this sectio we will assume a probability space (Ω, F, P), i which is defied a ifiite sequece of radom variables (X ) ad a radom variable X. The fact that for every ifiite

More information

The picture in figure 1.1 helps us to see that the area represents the distance traveled. Figure 1: Area represents distance travelled

The picture in figure 1.1 helps us to see that the area represents the distance traveled. Figure 1: Area represents distance travelled 1 Lecture : Area Area ad distace traveled Approximatig area by rectagles Summatio The area uder a parabola 1.1 Area ad distace Suppose we have the followig iformatio about the velocity of a particle, how

More information

Random Walks on Discrete and Continuous Circles. by Jeffrey S. Rosenthal School of Mathematics, University of Minnesota, Minneapolis, MN, U.S.A.

Random Walks on Discrete and Continuous Circles. by Jeffrey S. Rosenthal School of Mathematics, University of Minnesota, Minneapolis, MN, U.S.A. Radom Walks o Discrete ad Cotiuous Circles by Jeffrey S. Rosethal School of Mathematics, Uiversity of Miesota, Mieapolis, MN, U.S.A. 55455 (Appeared i Joural of Applied Probability 30 (1993), 780 789.)

More information

Estimation of Backward Perturbation Bounds For Linear Least Squares Problem

Estimation of Backward Perturbation Bounds For Linear Least Squares Problem dvaced Sciece ad Techology Letters Vol.53 (ITS 4), pp.47-476 http://dx.doi.org/.457/astl.4.53.96 Estimatio of Bacward Perturbatio Bouds For Liear Least Squares Problem Xixiu Li School of Natural Scieces,

More information

Convergence of random variables. (telegram style notes) P.J.C. Spreij

Convergence of random variables. (telegram style notes) P.J.C. Spreij Covergece of radom variables (telegram style otes).j.c. Spreij this versio: September 6, 2005 Itroductio As we kow, radom variables are by defiitio measurable fuctios o some uderlyig measurable space

More information

Law of the sum of Bernoulli random variables

Law of the sum of Bernoulli random variables Law of the sum of Beroulli radom variables Nicolas Chevallier Uiversité de Haute Alsace, 4, rue des frères Lumière 68093 Mulhouse icolas.chevallier@uha.fr December 006 Abstract Let be the set of all possible

More information

Stochastic Simulation

Stochastic Simulation Stochastic Simulatio 1 Itroductio Readig Assigmet: Read Chapter 1 of text. We shall itroduce may of the key issues to be discussed i this course via a couple of model problems. Model Problem 1 (Jackso

More information