Prepared with an evaluation version of PCT E X. Visit

Size: px
Start display at page:

Download "Prepared with an evaluation version of PCT E X. Visit"

Transcription

1 cspringer (Original version in the Proceedings of FPSAC') A Generalized Cover Tie for Rando Walks on Graphs Cyril Banderier and Robert P. Dobrow 2 Algoriths Project. INRIA (Rocquencourt), France. Cyril.Banderier@inria.fr, 2 Clarkson University (Potsda, NY), USA. dobrowb@clarkson.edu, Abstract. Given a rando walk on a graph, the cover tie is the rst tie (nuber of steps) that every vertex has been hit (covered) by the walk. Dene the arking tie for the walk as follows. When the walk reaches vertex v i, a coin is ipped and with probability p i the vertex is arked (or colored). We study the tie that every vertex is arked. (When all the p i's are equal to, this gives the usual cover tie proble.) General forulas are given for the arking tie of a graph. Connections are ade with the generalized coupon collector's proble. Asyptotics for sall p i's are given. Techniques used include cobinatorics of rando walks, theory of deterinants, analysis and probabilistic considerations. Introduction The following proble was subitted during a supper at the eeting Analysis of Algoriths in June 999, at Barcelona: Conjecture (Supper Conjecture) Iagine guests around a table, soeone has the water carafe and decides to pour soe water in his glass with probability p. Then he gives the water carafe randoly to his right or left neighbor. This one does the sae, and so on. Call T (p) the nuber of carafe oves before everyone has got soe water. What can one say about the average tie E(T (p))? In particular, is it true that pe(t (p))! H when p!? (H is the -th haronic nuber). One will show that this conjecture is indeed true! In fact, we will tackle the proble for a slightly ore general proble: think about a dinner where soe people are ore or less not in speaking ters with soe others and so they give the carafe preferentially to their friends! The proble we consider was otivated by the gae Trivial Pursuit. Players ove pieces around a gae board answering questions on various topics (e.g., history, sports, etc.). On certain positions in the gae board, if a player answers a question correctly he gets a colored piece. And in our siplied version of the gae, when a player gets all the colored pieces he wins. Assue that for each topic there is soe probability of answering a question correctly. How long will a typical Trivial Pursuit gae take to play? Prepared with an evaluation version of PCT E. Visit

2 The reader should be able to see the connection with the following odel. Consider the usual (discrete) rando walk on a connected directed graph G with vertex set V = fv ; :::;v g. When the rando walk reaches vertex v i, that vertex is arked with probability p i, i = ;:::;. We are interested in the arking tie T (fp ; :::;p g), the rst tie that all the vertices have been arked. When all the p i 's are equal to, this is the usual cover tie proble. When the graph is the coplete graph and the rando walk oves to any vertex uniforly at rando, the proble reduces to the classical coupon collector's proble. (E.g., How any packs of Pokeon cards should you buy in order to collect the full set of Pokeon characters?) It is well know P that the expected k tie to cover a graph with vertices is H ; where H k := i= i is the kth Haronic nuber. The generalized coupon collector's proble asks for the tie to cover the coplete graph when the transition probabilities (weights) for the underlying rando walk are not constant (see the survey [2] for a lot of applications). There are classical results established either by cobinatorics (inclusion/exclusion principles [9], shue product [5]) or probability (artingales [8], tails estiates [7]). We also give soe references for the coupon collector proble for arbitrary graphs [4]; soe articles deal with asyptotic considerations [3]. Our proble (which involves two levels of randoness!) sees quite new and of course allows us to rederive/iprove previous results (in the peculiar case p=). When all the p i 's are equal (say, to p), a nice but nally predictable result is that pe(t (p)) is a rational nuber (as usual for Markovian process with rational transition probabilities). Of course, in principle all the questions concerning the rando walk can be addressed by solving appropriate linear equations using the transition atrix for the rando walk Markov chain but such an approach quickly becoes infeasible for even sall values of. On the coplete graph, in the case of equal arking probabilities (p i = p for all i), and of probabilities i for the coupons, a natural conjecture is that E[T (p)] (=p)h when p!. The reason is because for sall p, one would expect the tie to ark a particular vertex to be Geoetric with paraeter p (and thus ean =p). A quick (but erroneous) probabilistic reasoning consists in seeing the arking probability p as a change of tie-scale and thus one gets (for any graph) E(T (p)) = E(T ())=p. This is false (even asyptocally)! However, we will show this intuition to be true for regular graphs and prove a stronger result: For an arbitrary graph G, one has pe(t (p))! K as p!, where K is the expected cover tie for the generalized coupon collector's proble where the set of weights is the stationary distribution for the rando walk. This result will be proven in Section 4. Notation (Graph) Through all the paper, G is a directed connected graph with vertices v ; :::;v and with a transition atrix A (that is, a ij is the probability of a transition fro v i to v j ). The stationary distribution of this graph is noted ; :::;. The probability to ark the vertex v i whenever it is visited is In the USA, kids have a great craze for these cards. They trade the and stare at the a lot... :-) Prepared with an evaluation version of PCT E. Visit

3 Linear graph with reflexion at buttos (=8) Cyrclic graph (=8) Linear graph with loops at buttos (=8) Coplete graph (=4) Fig.. Soe of the oriented connected graphs considered in our exaples. noted p i (i =;:::;), and one sets q i := p i. When all the p i 's are equal, one siply notes p the probability of arking a vertex. e G is a coplete graph with transition atrix e A (related in soe sense to A, as explained later). Notation 2 (Rando variables) T (resp. e T ) is the \waiting tie" rando variable that represents the rst tie when all the vertices of the graph G (resp. eg) have been arked. The rando variable i (resp. e i ) gives the rst tie the i-th vertex of the graph G (resp. e G) is arked. Notation 3 (Operators) One notes [z n ]S the coecient of z n in a given series S. For any subset of f;:::;g, jj stands for the nuber of eleents in and one notes the substitution u i for i 2 and u i stays unchanged for i 62, the substitution u i q i for i 2 and u i for i 62. The sae notation will be used with respect to any other set of foral variables (for exaple with t ; :::;t instead of u ; :::;u ). We will also denote the identity atrix as Id and U as the diagonal atrix with diagonal eleents u ; :::;u. Exaples. To illustrate the above notations: [z 2 ] (+4z 2 + z 3 ) =4.For =4 and = f; 2; 4g, one has (3u + u u 2 + u u 4 ) = u 3 3 and (3u + u u 2 + u u 4 ) =3q + q q 2 + +q 4. In the sequel, all the graphs considered have vertices ( > ) and are directed and irreducible: there is a sequence of connected edges linking any pair of vertices. First, we establish a cobinatorial forula for the arking tie in section 2. One will explain in section 3 a probabilistic intuition which allows us to siplify the proble and thus to prove the Supper Conjecture in Section 4, another proof is given in Section 5. The last sections deal with peculiar cases. 2 Marking Tie on any Graph We rst give a generating function-based forula (based on nite dierences) for E(T (p)), the expected arking tie. For the cover tie on the coplete graph, another approach can be found in [5]. We refer the reader to Figure 4 for an exaple of our approach on two graphs of size 3. Prepared with an evaluation version of PCT E. Visit

4 Theore (General forula for average arking tie) The average tie for arking all the vertices of an arbitrary graph G is given by E(T (fp ; :::;p g)) = ( ) jj+ u (Id AU) i ; 6=; i= where A, U, the p i 's, and the 's are dened as in Notation and 3. Proof. For the transition atrix A, the entry A n i;j of A n, gives the probability of oving fro vertex v i to v j in n steps. With the atrix U dened as in Notation 3, the coecient of the onoial u k uk in (AU) n i;j gives the probability that the rando walk oves fro v i to v j in n steps such that vertex v t is visited k t ties, t =;:::;. Thus the probability generating function for the walks on the graph (beginning in v ), where z encodes the length of the walk and the u i 's encode the nuber of ties the walk visits v i is F (z;u ; :::;u ) = u z k (AU) k ;i = u (Id zau) ;i : i= k= i= Recall that p i (respectively, q i := p i ) is the probability to ark (respectively, not to ark) the vertex v i. Taking into account the fact that the walk has visited all the vertices arking the at least once leads to the substitution u n i qi n. This justi- es the introduction of the dierence operator i f(u i ) := f() f(q i ): So F + := 2 : :: F gives the probability generating function of the walks that arked all the vertices and one has in fact F + (z) = f;:::;g ( ) jj F (z;u ; :::;u ); i.e., F + is the su of F evaluated the set (dened in Notation 3). Therefore the probability generating function for T (fp ; :::;p g) is ( z)f + (z), so E(T (fp ; :::;p jz=( z)f + (z): A change of variable z = t and a local developent in t = ( z)u = (Id AU) jz= Id zau ij Note that Id A is never invertible (whereas Id za is always invertible), so one has to deal apart with the substitution ;. ut Note that setting the p i 's to gives a forula for the coupon collector proble. 3 The Probabilistic Intuition We restrict here the discussion to the case when all the p i 's are sall (equivalently, you can think p sall with p := ax p i ). One argues that the rando walk behaves, for sall values of p, like a walk on the coplete graph (with loops). We obtain a new Markov chain that is \equivalent" to the original Markov Chain, in the sense that they both have the sae (liiting) stationary distributions. Prepared with an evaluation version of PCT E. Visit

5 Let f i g i= be the stationary distribution for the rando walk on G. Note that i gives the probability that the walk is in vertex v i in stationarity (i.e., after a \long tie"). When p is sall, the tie to ark the graph is high. And the proportion of tie n that the walk is in vertex v i will be i n + O( p n) for large n, with probability close to. (This can be shown, for instance, by a central liit theore for Markov chains, or fro results on large deviations.) Thus, for \long enough" walks over G (when p is sall, all walks are \long enough"), the average length of tie to ark all the vertices should be very nearly the average length of tie to ark the vertices of the coplete graph (with loops), but with transition probabilities corresponding to the stationary distribution for the original graph. Let G e denote the coplete graph on vertices with loops. Dene its transition atrix A e with the transition probabilities ea i;j := j (for i; j =;:::;). Consider a rando walk process for the general arking proble that begins in the vertex v i with probability i. For such a process de- ne k e to be the rst tie that vertex ev k is arked. Then the k e are geoetric rando variables with paraeter k p k and P ( k e = n) = k p k ( k p k ) n, but there are not independent (as i ~ 6= j ~ for i 6= j). Observe that T = ax( ; :::; ) and T e = ax( e ; :::; e ). Consider now anew process on the graph. Instead of one \rando walker," consider particles at each of the vertices all oving siultaneously to neighboring vertices according to the sae transition echanis, but independently of each other. If we let Y i (and Yi e ) denote the rst tie that vertex v i is arked then observe that the Y i has the sae distribution as i but the Y i 's are independent and the i 's are not. Dene now Z = ax(y ; :::;Y ) and Z e = ax( Y e ; :::;e Y ). Our intuition is as follows: As the Y e i s and the Y i 's will behave siilarly for sall p, one should have that E( Z) e E(Z): As siultaneous arkings (for the Y process) occur with probability O(p 2 ), one has that E(Z) E(T ). Also E( Z) e E( T e ) and thus E(T ) E( T e ). Thus the study of the expected arking tie on the graph G e should answer our question as to the expected arking tie on the graph G. We ake the above rigorous in the next section. 4 Algebraic and Cobinatorial Proof Theore 2 (Closed for forula with stationary distribution) Let be a graph with transition probabilities ea i;j := j. The expected tie E(e T ) for arking the whole graph G e is =p ties the tie needed for visiting all the vertices E(e T )= p E(coupon collector on G); e and an inclusion-exclusion forula holds E(e T )= p 6=; where the 's are dened as in Notation 3. ( ) jj ( :::+ ) ; e G Prepared with an evaluation version of PCT E. Visit

6 Proof. Consider the i s, the q i's, and the p i 's as foral variables, then Prob(e T n) =fp > g :::fp> g( p + q + :::+ p + q ) n ; where fp >Sg (with S 2 IR[[p ; :::;p ; q ; :::;q ]]) stands for the su of onoials in p of positive exponent in S. Thus, we have Prob(e T n) = f;:::;g ( ) jj ( ( p + q + :::+ p + q )) n ; where the substitutions are dened as in Notation 3. Multiplying by z n and suing for n, when all the p i 's are equal to p, with the substitutions now taken to act on the i 's, leads to Prob( e T n)z n = f;:::;g ( ) jj z(q + p ( :::+ )) : Multiplying by z and dierentiating in z = gives the expected tie E(e T )= 6=; ( ) jj p( + ( :::+ )) : ut graph G graph G e /2 /2 3/8 3/8 /4 transition atrix A = /3 /3 /3 ea = 3/8 3/8 /4 /4 /4 /2 3/8 3/8 /4 stationary distribution (3/8, 3/8, /4) (3/8, 3/8, /4) u! u! U = u 2 U = u 2 u 3 u 3 generating atrix of the walks (Id zau) (Id z AU) e expected covering tie 95=2 ' 7:9 29=5 ' 5:8 expected arking tie (closed for) expected arking tie (asyptotics)! p + 87p p 3 75p +6p 2 + 9p 3 29=5p E(T (p)) = 29=5p + O() E(e T (p)) = 29=5p Fig. 2. Exaple of our approach on a graph of size 3. Most of the results in the literature are about the (expected) covering tie of the e G-colun (and other higher oents). In our situation (the G-colun), the lack of independence is the ain diculty. Our paper shows that in order to get the last entries in the G-colun, one can proceed as follow: Consider the rst 2 entries of the G-colun, then deal with eg-colun whose last entries gives the wanted result (for G).! Prepared with an evaluation version of PCT E. Visit

7 Proposition The arking tie T can be approxiated by e T : E(e T ) E(T )=O(); (p! ): Proof. In order to prove that 6=;( ) jj+ i= u (Id AU) e i ( ) jj+ u (Id AU) i = O() 6=; i= where the substitutions are dened in Notation 3, it is sucient to check that [p ] Id e AU In fact, one has ( za) = ea ij [p ] (Id AU) ij = O(): () z + O(): This is established by writing A = P + R, which decoposes A as an eigenprojection on the eigensubspace related to its eigenvector, plus the projection R on the suppleentary subspace (as it is well known P n = P = e A and PR = RP = ), thus P z n (P + R) n = ea z + (Id zr) where (Id zr) is actually regular in z = (as, by Perron{ Frobenius theory, is an eigenvalue of ultiplicity ), so ( zr) = O(). So this sets the case of the equality () when = f;:::;g (siply consider z = p). Theq other cases appear as a perturbation of two noninvertible atrices (naely Id A and Id e A), which gives (Id A + B) = ea () + O() = (Id e A + e B) where B is A with its i-th colun set to whenever i 2 (so A B = P AU), and where () is the perturbation of the eigenvalue, so () = i2 i. See [6] for the analicity of the perturbations of projections, eigenvalues, etc. ut For people with probabilistic anities, we give below another approach which also proves that E(T )=E( e T )+O() as p!. 5 Analytical and Probabilistic Proof The following proposition (identical to Proposition ) is the key point Proposition 2 (Equivalence of expected ties) The arking tie for the graph G and G e have the sae rst order asyptotics, naely Proof. The ain idea is the following E(T )=E(e T )+O(): T (T without its tails) (e T without its tails ) e T ; where eans here that the expectations have the sae rst order asyptotics as p!. Note that P (T n) = p c ; =(n ;:::;n) Prepared with an evaluation version of PCT E. Visit

8 Pwhere the su is over all nonnegative -tuples = (n ; :::;n ) such that k n k = n, p is the probability that the walk takes n steps and visits each vertex v k n k ties, k =;:::;, and c Q is the probability that all vertices have been arked by such a tour, that is c = interval" I as I := pj ln pj ( ( p)nj j= ). Dene the \central : ; j ln pj p Further dene the \ultidiensional box" B as B := Y i= [n i p n ln n; n i + p n ln n]: For sall p and for large n, by a classic result in large deviation theory, one has n j p n ln n<n j < n j + p n ln n (that is, 2 B) with probability exp( c ln 2 n) (with c >, see []). As Prob(T = n) is the probability that a success occurs exactly at the n-th step, one has E(T )= n( Prob(T n)): This su can be split as follows (with a lot of abusive but natural notations!) E(T )= (2) n<i n2i;2b n2i;62b n>i;2b n>i;62b where, for exaple, n < I eans for n before the \central interval" I. The rst su is bounded by the length of the interval of suation, which is o(=p) as p!, the sus for 62 B are bounded by exp( c ln 2 n) and so is the reaining su for n>i. One now focuses on the su over I and B. By a liit theore on Markov Chains [], p follows a ultidiensional Gaussian law g(), thus E(T )= n2i 2B g()c! + o(p ): (3) The sae schee can be applied to G, e with the sae central interval I and box B since the two rando walks have the sae stationary distribution. Thus E(T )=E(e T )+o(p )=E(e T )+O(), since E( T e ) and E(T ) are rational fractions in p (see Theore and 2). ut The Theore below is the ain result of the paper and answers precisely to the Supper Conjecture. The following \integral forula" is well known (cf. [5], where it is established for the coplete graph by eans of shue products and Laplace transfor), but we give an alternate derivation here, valid for any graph. Prepared with an evaluation version of PCT E. Visit

9 Theore 3 (First order asyptotics. Integral for) For any graph G with a stationary distribution ( ; :::; ), the expected arking tie is Z E(T )= K p + O(); where K = ( Y j= Proof. With the sae notations as above, the starting point is c = c (n;:::;n) = where p := ln( p). Thus E(T )= n2i Y 2B ( ( p) nj ) = g() Y j= Y j= ( exp( j x))) dx: ( e p nj ); exp( p n j j ) A + O(): In the box B, n j = n j + j p n ln n (where jj j < ), so E(T )= n2i c = Y ( exp( p n j p j p n ln n)); and thus 2B g() Y j= Now, ajoring P g() by exp( c ln 2 n) yields E(T )= n2i ( exp( p n j p j p n ln n)) A + O(): exp( p n j p p j n ln n) A + O() + o(p ): j= The sus being for p j ln pj n p j ln pj, this leads to exp( p (O(n))( exp( p o(n))) ( exp( p o(n)) and copleting the tails gives p o(n) p p j ln pjo(p j ln pj) =o(p ); E(T )= n2i = n( exp( p n j ) j= Y j= A + O() + o(p ): ( exp( p n j ))) + O() + o(p ): (4) Indeed the P rst introduced part is p =j ln pj = o(p ) and the last introduced part is n>p j ln Q pj 2 exp( p in( i )) n = o(p ). Set f(x) := ( j= exp( p x j ). As f (x) has a fast enough decay to near, the Euler-Maclaurin forula gives: Prepared with an evaluation version of PCT E. Visit

10 n= f(n) Z f(x)dx = f() + f() 2 Applying this to the forula (4) leads to Z E(T )= p ( Y j= + Z Since p = p + o(p), one has E(T )=K=p + O(). (x bxc =2)f (x)dx = O(): ( exp( j x))) dx + O() when p! : The \integral forula" allows us to copute E(T ) to any precision in linear tie (as all the integrated functions have a nice behavior), whereas the forulae of Theore and 2 are ipracticable because they coprise 2 suands. 6 What about Balanced or Regular Graphs? The ost coon odel for rando walk on a graph (the unweighted case) is to suppose that at a particular vertex the walk oves with probability proportional to the degree of that vertex. It is well known that the stationary distribution for such a walk is related to the degree of the vertices. Call a graph \balanced" if the outdegree of each vertex is equal to its indegree and if each outgoing edge is equally likely. For balanced graphs, there is a siple relation between the stationary distribution (the left eigenvector associated to eigenvalue of the transition atrix of the graph) and the degrees of edges: Proposition 3 (Stationary distribution for balanced graphs) For balanced graphs, one has i = P Ni Ni (N i is the nuber of incoing edges of vertex v i ). Proof. Let N ij be the nuber of edges fro v i to v j, and note N j the nuber of incoing edges to v j and N i the nuber of outgoing edges fro v i. The left eigenvector (for eigenvalue ) satises C A = ( ; 2 ; :::; ); N =N : :: N =N ( ; 2 ; :::; )... P N =N : ::N =N hence i Nij i = j. This equation is indeed satised by Ni i := P Ni Ni when N i = N i. ut Theore 2 rewrites in these cases: P Corollary (Balanced graphs) For balanced graphs with N := N i edges, one has E(T )= N ( ) jj p (N + N 2 + :::+ N ) + O(); 6=; where the substitutions operate on the fn j g. Proof. Direct consequence of Theore 2 and Proposition 3. ut ut Prepared with an evaluation version of PCT E. Visit

11 There are several classes of graphs for which the forula can be siplied. The ore interesting one, the class of regular graphs (graphs whose all vertices have the sae nuber of incoing and outgoing edges), is the object of the following corollary. Corollary 2 (Regular graphs) When all the i 's are equal (in particular if G is a regular graph and outgoing edges are chosen uniforly at rando), one has E(T )= H p + O(). Proof. ( ) i i p( + i ) = ( ) i+ p i i = H : p E(e T )= i= The last equality follows as iteration of forward nite dierence operators and Euler's R transfor. Q Equivalently, one could use Theore 3. Here the integral becoes ( ( exp( x=))) dx, which siplies to H j=. A third proof of this forula is probabilistic and considers the axiu of i.i.d. geoetric rando variables. Thus when i = = (for i =;:::;), one has Prob(e Z = n) = Prob( i e Y = n) Prob( j e < n) j6=i i i= = p=( p=) n ( ( p=) n ) : Since E(e T )= Pn nprob( e T = n), one is interested by the coecient of the lowest ter in the Laurent developent at p =of n nz( z )n ( ( z )n ) = k= ( ) k k z ( ( z=) k+ ) 2 : As ( ( z=) k+ ) 2 = O(z 2 ) as z!, the valuation of the developent in Laurent series of z ( z=) k+ 2 is. Thus the coecient of the lowest ter of F (z) is the su of the residues of the rational fractions, then k Res(z =)= ( ) k k= (k +) 2 = H : 7 Exaples The forula for the expected tie in the coupon collector proble reduces to for the cyclic graph C, to for the line graph L (with reection) and to ( )H + for the coplete graph without loops K. For L, there is an equivalence between the coupon collector proble and the rando walks (with jups +; ) fro to of height, the last ones being linked to continued fraction theory and hence to a quotient of Chebyshev polynoials. We apply below the forula of Theore for all nonisoorphic unoriented connected graphs without ultiplicity having at ost 4 vertices. We give for each graph, the stationary distribution, the average tie for the classical coupon Prepared with an evaluation version of PCT E. Visit

12 collector proble and, in the last colun, the coecient of p (the leading coecient) in the asyptotics of the average tie for arking the whole graph: size graph stationary expected expected distributions cover tie arking tie = K () =2 K 2 (/2, /2) 2 3 =3 K 3 (/3, /3, /3) 4 3H 3 =/2 =3 L 3 (/4, /2, /4) 5 9/3 =4 L 4 (/6, /3, /3, /6) 99/ =4 K 4 (/4, /4, /4, /4) 3/2 4H 4 =25/3 =4 C 4 (/4, /4, /4, /4) 7 4H 4 =25/3 =4 Y 4 (/5, 2/5, /5, /5) / =4 T 4 (/8, 3/8, /4, /4) 63/5 369/35 =4 Q 4 (3/, /5, /5, 3/) 69/ 62/7 The graphs K, L and C are dened as above whereas T 4 is the triangle with a tail, Y 4 is the Y shape, Q 4 is the square with one diagonal. We recall that, as seen in previous sections, for the coplete graph with loops, the coupon collector tie is H and the arking tie is H p. This is nally the behavior of all regular graphs, those for which the Supper Conjecture holds. Acknowledgeent. The rst of the two authors would like to thank Philippe Flajolet for a lot of helpful rearks and does not resist quoting one of his pseudo-reforulations of our proble \You are a PhD student and you attept to prove a theore; your chance of success at each attept is p (p very close to, of course). Let T be the tie it takes you to get the theore right..." References. Aldous (David J.) and Fill (Jaes). { Reversible Markov Chains and Rando Walks on Graphs. { Book in preparation. Available at 2. Boneh (Arnon) and Hofri (Micha). { The coupon-collector proble revisited. Co. Statist. Stochastic Models, vol. 3, n, 997, pp. 39{ Csorg}o (Sandor). { A rate of convergence for coupon collectors. Acta Sci. Math. (Szeged), vol. 57, n-4, 993, pp. 337{ Feige (Uriel). { Collecting coupons on trees, and the cover tie of rando walks. Coput. Coplexity, vol. 6, n4, 996=97, pp. 34{ Flajolet (Philippe), Gardy (Daniele), and Thionier (Loys). { Birthday paradox, coupon collectors, caching algoriths and self-organizing search. Discrete Applied Matheatics, vol. 39, n3, 992, pp. 27{ Kato (Tosio). { Perturbation theory for linear operators. { Springer-Verlag, Nath (Harindar B.). { Waiting tie in the coupon-collector's proble. Austral. J. Statist., vol. 5, 973, pp. 32{ Sen (Pranab Kuar). { Invariance principles for the coupon collector's proble: a artingale approach. Ann. Statist., vol. 7, n2, 979, pp. 372{ von Schelling (Herann). { Coupon collecting for unequal probabilities. Aer. Math. Monthly, vol. 6, 954, pp. 36{3. Prepared with an evaluation version of PCT E. Visit

13.2 Fully Polynomial Randomized Approximation Scheme for Permanent of Random 0-1 Matrices

13.2 Fully Polynomial Randomized Approximation Scheme for Permanent of Random 0-1 Matrices CS71 Randoness & Coputation Spring 018 Instructor: Alistair Sinclair Lecture 13: February 7 Disclaier: These notes have not been subjected to the usual scrutiny accorded to foral publications. They ay

More information

Solutions of some selected problems of Homework 4

Solutions of some selected problems of Homework 4 Solutions of soe selected probles of Hoework 4 Sangchul Lee May 7, 2018 Proble 1 Let there be light A professor has two light bulbs in his garage. When both are burned out, they are replaced, and the next

More information

. The univariate situation. It is well-known for a long tie that denoinators of Pade approxiants can be considered as orthogonal polynoials with respe

. The univariate situation. It is well-known for a long tie that denoinators of Pade approxiants can be considered as orthogonal polynoials with respe PROPERTIES OF MULTIVARIATE HOMOGENEOUS ORTHOGONAL POLYNOMIALS Brahi Benouahane y Annie Cuyt? Keywords Abstract It is well-known that the denoinators of Pade approxiants can be considered as orthogonal

More information

Block designs and statistics

Block designs and statistics Bloc designs and statistics Notes for Math 447 May 3, 2011 The ain paraeters of a bloc design are nuber of varieties v, bloc size, nuber of blocs b. A design is built on a set of v eleents. Each eleent

More information

CSE525: Randomized Algorithms and Probabilistic Analysis May 16, Lecture 13

CSE525: Randomized Algorithms and Probabilistic Analysis May 16, Lecture 13 CSE55: Randoied Algoriths and obabilistic Analysis May 6, Lecture Lecturer: Anna Karlin Scribe: Noah Siegel, Jonathan Shi Rando walks and Markov chains This lecture discusses Markov chains, which capture

More information

A note on the multiplication of sparse matrices

A note on the multiplication of sparse matrices Cent. Eur. J. Cop. Sci. 41) 2014 1-11 DOI: 10.2478/s13537-014-0201-x Central European Journal of Coputer Science A note on the ultiplication of sparse atrices Research Article Keivan Borna 12, Sohrab Aboozarkhani

More information

Chaotic Coupled Map Lattices

Chaotic Coupled Map Lattices Chaotic Coupled Map Lattices Author: Dustin Keys Advisors: Dr. Robert Indik, Dr. Kevin Lin 1 Introduction When a syste of chaotic aps is coupled in a way that allows the to share inforation about each

More information

The Weierstrass Approximation Theorem

The Weierstrass Approximation Theorem 36 The Weierstrass Approxiation Theore Recall that the fundaental idea underlying the construction of the real nubers is approxiation by the sipler rational nubers. Firstly, nubers are often deterined

More information

Feature Extraction Techniques

Feature Extraction Techniques Feature Extraction Techniques Unsupervised Learning II Feature Extraction Unsupervised ethods can also be used to find features which can be useful for categorization. There are unsupervised ethods that

More information

Handout 7. and Pr [M(x) = χ L (x) M(x) =? ] = 1.

Handout 7. and Pr [M(x) = χ L (x) M(x) =? ] = 1. Notes on Coplexity Theory Last updated: October, 2005 Jonathan Katz Handout 7 1 More on Randoized Coplexity Classes Reinder: so far we have seen RP,coRP, and BPP. We introduce two ore tie-bounded randoized

More information

Quantum algorithms (CO 781, Winter 2008) Prof. Andrew Childs, University of Waterloo LECTURE 15: Unstructured search and spatial search

Quantum algorithms (CO 781, Winter 2008) Prof. Andrew Childs, University of Waterloo LECTURE 15: Unstructured search and spatial search Quantu algoriths (CO 781, Winter 2008) Prof Andrew Childs, University of Waterloo LECTURE 15: Unstructured search and spatial search ow we begin to discuss applications of quantu walks to search algoriths

More information

The Euler-Maclaurin Formula and Sums of Powers

The Euler-Maclaurin Formula and Sums of Powers DRAFT VOL 79, NO 1, FEBRUARY 26 1 The Euler-Maclaurin Forula and Sus of Powers Michael Z Spivey University of Puget Sound Tacoa, WA 98416 spivey@upsedu Matheaticians have long been intrigued by the su

More information

Bipartite subgraphs and the smallest eigenvalue

Bipartite subgraphs and the smallest eigenvalue Bipartite subgraphs and the sallest eigenvalue Noga Alon Benny Sudaov Abstract Two results dealing with the relation between the sallest eigenvalue of a graph and its bipartite subgraphs are obtained.

More information

}, (n 0) be a finite irreducible, discrete time MC. Let S = {1, 2,, m} be its state space. Let P = [p ij. ] be the transition matrix of the MC.

}, (n 0) be a finite irreducible, discrete time MC. Let S = {1, 2,, m} be its state space. Let P = [p ij. ] be the transition matrix of the MC. Abstract Questions are posed regarding the influence that the colun sus of the transition probabilities of a stochastic atrix (with row sus all one) have on the stationary distribution, the ean first passage

More information

arxiv: v1 [math.pr] 17 May 2009

arxiv: v1 [math.pr] 17 May 2009 A strong law of large nubers for artingale arrays Yves F. Atchadé arxiv:0905.2761v1 [ath.pr] 17 May 2009 March 2009 Abstract: We prove a artingale triangular array generalization of the Chow-Birnbau- Marshall

More information

Lecture 9 November 23, 2015

Lecture 9 November 23, 2015 CSC244: Discrepancy Theory in Coputer Science Fall 25 Aleksandar Nikolov Lecture 9 Noveber 23, 25 Scribe: Nick Spooner Properties of γ 2 Recall that γ 2 (A) is defined for A R n as follows: γ 2 (A) = in{r(u)

More information

Model Fitting. CURM Background Material, Fall 2014 Dr. Doreen De Leon

Model Fitting. CURM Background Material, Fall 2014 Dr. Doreen De Leon Model Fitting CURM Background Material, Fall 014 Dr. Doreen De Leon 1 Introduction Given a set of data points, we often want to fit a selected odel or type to the data (e.g., we suspect an exponential

More information

arxiv: v1 [math.nt] 14 Sep 2014

arxiv: v1 [math.nt] 14 Sep 2014 ROTATION REMAINDERS P. JAMESON GRABER, WASHINGTON AND LEE UNIVERSITY 08 arxiv:1409.411v1 [ath.nt] 14 Sep 014 Abstract. We study properties of an array of nubers, called the triangle, in which each row

More information

Polygonal Designs: Existence and Construction

Polygonal Designs: Existence and Construction Polygonal Designs: Existence and Construction John Hegean Departent of Matheatics, Stanford University, Stanford, CA 9405 Jeff Langford Departent of Matheatics, Drake University, Des Moines, IA 5011 G

More information

A1. Find all ordered pairs (a, b) of positive integers for which 1 a + 1 b = 3

A1. Find all ordered pairs (a, b) of positive integers for which 1 a + 1 b = 3 A. Find all ordered pairs a, b) of positive integers for which a + b = 3 08. Answer. The six ordered pairs are 009, 08), 08, 009), 009 337, 674) = 35043, 674), 009 346, 673) = 3584, 673), 674, 009 337)

More information

Sequence Analysis, WS 14/15, D. Huson & R. Neher (this part by D. Huson) February 5,

Sequence Analysis, WS 14/15, D. Huson & R. Neher (this part by D. Huson) February 5, Sequence Analysis, WS 14/15, D. Huson & R. Neher (this part by D. Huson) February 5, 2015 31 11 Motif Finding Sources for this section: Rouchka, 1997, A Brief Overview of Gibbs Sapling. J. Buhler, M. Topa:

More information

A Note on Scheduling Tall/Small Multiprocessor Tasks with Unit Processing Time to Minimize Maximum Tardiness

A Note on Scheduling Tall/Small Multiprocessor Tasks with Unit Processing Time to Minimize Maximum Tardiness A Note on Scheduling Tall/Sall Multiprocessor Tasks with Unit Processing Tie to Miniize Maxiu Tardiness Philippe Baptiste and Baruch Schieber IBM T.J. Watson Research Center P.O. Box 218, Yorktown Heights,

More information

4 = (0.02) 3 13, = 0.25 because = 25. Simi-

4 = (0.02) 3 13, = 0.25 because = 25. Simi- Theore. Let b and be integers greater than. If = (. a a 2 a i ) b,then for any t N, in base (b + t), the fraction has the digital representation = (. a a 2 a i ) b+t, where a i = a i + tk i with k i =

More information

Topic 5a Introduction to Curve Fitting & Linear Regression

Topic 5a Introduction to Curve Fitting & Linear Regression /7/08 Course Instructor Dr. Rayond C. Rup Oice: A 337 Phone: (95) 747 6958 E ail: rcrup@utep.edu opic 5a Introduction to Curve Fitting & Linear Regression EE 4386/530 Coputational ethods in EE Outline

More information

The concavity and convexity of the Boros Moll sequences

The concavity and convexity of the Boros Moll sequences The concavity and convexity of the Boros Moll sequences Ernest X.W. Xia Departent of Matheatics Jiangsu University Zhenjiang, Jiangsu 1013, P.R. China ernestxwxia@163.co Subitted: Oct 1, 013; Accepted:

More information

arxiv: v1 [cs.ds] 3 Feb 2014

arxiv: v1 [cs.ds] 3 Feb 2014 arxiv:40.043v [cs.ds] 3 Feb 04 A Bound on the Expected Optiality of Rando Feasible Solutions to Cobinatorial Optiization Probles Evan A. Sultani The Johns Hopins University APL evan@sultani.co http://www.sultani.co/

More information

New upper bound for the B-spline basis condition number II. K. Scherer. Institut fur Angewandte Mathematik, Universitat Bonn, Bonn, Germany.

New upper bound for the B-spline basis condition number II. K. Scherer. Institut fur Angewandte Mathematik, Universitat Bonn, Bonn, Germany. New upper bound for the B-spline basis condition nuber II. A proof of de Boor's 2 -conjecture K. Scherer Institut fur Angewandte Matheati, Universitat Bonn, 535 Bonn, Gerany and A. Yu. Shadrin Coputing

More information

Using EM To Estimate A Probablity Density With A Mixture Of Gaussians

Using EM To Estimate A Probablity Density With A Mixture Of Gaussians Using EM To Estiate A Probablity Density With A Mixture Of Gaussians Aaron A. D Souza adsouza@usc.edu Introduction The proble we are trying to address in this note is siple. Given a set of data points

More information

Study on Markov Alternative Renewal Reward. Process for VLSI Cell Partitioning

Study on Markov Alternative Renewal Reward. Process for VLSI Cell Partitioning Int. Journal of Math. Analysis, Vol. 7, 2013, no. 40, 1949-1960 HIKARI Ltd, www.-hikari.co http://dx.doi.org/10.12988/ia.2013.36142 Study on Markov Alternative Renewal Reward Process for VLSI Cell Partitioning

More information

Poly-Bernoulli Numbers and Eulerian Numbers

Poly-Bernoulli Numbers and Eulerian Numbers 1 2 3 47 6 23 11 Journal of Integer Sequences, Vol. 21 (2018, Article 18.6.1 Poly-Bernoulli Nubers and Eulerian Nubers Beáta Bényi Faculty of Water Sciences National University of Public Service H-1441

More information

Completeness of Bethe s states for generalized XXZ model, II.

Completeness of Bethe s states for generalized XXZ model, II. Copleteness of Bethe s states for generalized XXZ odel, II Anatol N Kirillov and Nadejda A Liskova Steklov Matheatical Institute, Fontanka 27, StPetersburg, 90, Russia Abstract For any rational nuber 2

More information

M ath. Res. Lett. 15 (2008), no. 2, c International Press 2008 SUM-PRODUCT ESTIMATES VIA DIRECTED EXPANDERS. Van H. Vu. 1.

M ath. Res. Lett. 15 (2008), no. 2, c International Press 2008 SUM-PRODUCT ESTIMATES VIA DIRECTED EXPANDERS. Van H. Vu. 1. M ath. Res. Lett. 15 (2008), no. 2, 375 388 c International Press 2008 SUM-PRODUCT ESTIMATES VIA DIRECTED EXPANDERS Van H. Vu Abstract. Let F q be a finite field of order q and P be a polynoial in F q[x

More information

A new type of lower bound for the largest eigenvalue of a symmetric matrix

A new type of lower bound for the largest eigenvalue of a symmetric matrix Linear Algebra and its Applications 47 7 9 9 www.elsevier.co/locate/laa A new type of lower bound for the largest eigenvalue of a syetric atrix Piet Van Mieghe Delft University of Technology, P.O. Box

More information

a a a a a a a m a b a b

a a a a a a a m a b a b Algebra / Trig Final Exa Study Guide (Fall Seester) Moncada/Dunphy Inforation About the Final Exa The final exa is cuulative, covering Appendix A (A.1-A.5) and Chapter 1. All probles will be ultiple choice

More information

Analysis of Polynomial & Rational Functions ( summary )

Analysis of Polynomial & Rational Functions ( summary ) Analysis of Polynoial & Rational Functions ( suary ) The standard for of a polynoial function is ( ) where each of the nubers are called the coefficients. The polynoial of is said to have degree n, where

More information

This model assumes that the probability of a gap has size i is proportional to 1/i. i.e., i log m e. j=1. E[gap size] = i P r(i) = N f t.

This model assumes that the probability of a gap has size i is proportional to 1/i. i.e., i log m e. j=1. E[gap size] = i P r(i) = N f t. CS 493: Algoriths for Massive Data Sets Feb 2, 2002 Local Models, Bloo Filter Scribe: Qin Lv Local Models In global odels, every inverted file entry is copressed with the sae odel. This work wells when

More information

ORIGAMI CONSTRUCTIONS OF RINGS OF INTEGERS OF IMAGINARY QUADRATIC FIELDS

ORIGAMI CONSTRUCTIONS OF RINGS OF INTEGERS OF IMAGINARY QUADRATIC FIELDS #A34 INTEGERS 17 (017) ORIGAMI CONSTRUCTIONS OF RINGS OF INTEGERS OF IMAGINARY QUADRATIC FIELDS Jürgen Kritschgau Departent of Matheatics, Iowa State University, Aes, Iowa jkritsch@iastateedu Adriana Salerno

More information

The degree of a typical vertex in generalized random intersection graph models

The degree of a typical vertex in generalized random intersection graph models Discrete Matheatics 306 006 15 165 www.elsevier.co/locate/disc The degree of a typical vertex in generalized rando intersection graph odels Jerzy Jaworski a, Michał Karoński a, Dudley Stark b a Departent

More information

Curious Bounds for Floor Function Sums

Curious Bounds for Floor Function Sums 1 47 6 11 Journal of Integer Sequences, Vol. 1 (018), Article 18.1.8 Curious Bounds for Floor Function Sus Thotsaporn Thanatipanonda and Elaine Wong 1 Science Division Mahidol University International

More information

Combining Classifiers

Combining Classifiers Cobining Classifiers Generic ethods of generating and cobining ultiple classifiers Bagging Boosting References: Duda, Hart & Stork, pg 475-480. Hastie, Tibsharini, Friedan, pg 246-256 and Chapter 10. http://www.boosting.org/

More information

Physics 139B Solutions to Homework Set 3 Fall 2009

Physics 139B Solutions to Homework Set 3 Fall 2009 Physics 139B Solutions to Hoework Set 3 Fall 009 1. Consider a particle of ass attached to a rigid assless rod of fixed length R whose other end is fixed at the origin. The rod is free to rotate about

More information

In this chapter, we consider several graph-theoretic and probabilistic models

In this chapter, we consider several graph-theoretic and probabilistic models THREE ONE GRAPH-THEORETIC AND STATISTICAL MODELS 3.1 INTRODUCTION In this chapter, we consider several graph-theoretic and probabilistic odels for a social network, which we do under different assuptions

More information

Chapter 6 1-D Continuous Groups

Chapter 6 1-D Continuous Groups Chapter 6 1-D Continuous Groups Continuous groups consist of group eleents labelled by one or ore continuous variables, say a 1, a 2,, a r, where each variable has a well- defined range. This chapter explores:

More information

Probability Distributions

Probability Distributions Probability Distributions In Chapter, we ephasized the central role played by probability theory in the solution of pattern recognition probles. We turn now to an exploration of soe particular exaples

More information

Least Squares Fitting of Data

Least Squares Fitting of Data Least Squares Fitting of Data David Eberly, Geoetric Tools, Redond WA 98052 https://www.geoetrictools.co/ This work is licensed under the Creative Coons Attribution 4.0 International License. To view a

More information

arxiv:math/ v1 [math.co] 22 Jul 2005

arxiv:math/ v1 [math.co] 22 Jul 2005 Distances between the winning nubers in Lottery Konstantinos Drakakis arxiv:ath/0507469v1 [ath.co] 22 Jul 2005 16 March 2005 Abstract We prove an interesting fact about Lottery: the winning 6 nubers (out

More information

A Bernstein-Markov Theorem for Normed Spaces

A Bernstein-Markov Theorem for Normed Spaces A Bernstein-Markov Theore for Nored Spaces Lawrence A. Harris Departent of Matheatics, University of Kentucky Lexington, Kentucky 40506-0027 Abstract Let X and Y be real nored linear spaces and let φ :

More information

Physics 215 Winter The Density Matrix

Physics 215 Winter The Density Matrix Physics 215 Winter 2018 The Density Matrix The quantu space of states is a Hilbert space H. Any state vector ψ H is a pure state. Since any linear cobination of eleents of H are also an eleent of H, it

More information

Lecture 21 Nov 18, 2015

Lecture 21 Nov 18, 2015 CS 388R: Randoized Algoriths Fall 05 Prof. Eric Price Lecture Nov 8, 05 Scribe: Chad Voegele, Arun Sai Overview In the last class, we defined the ters cut sparsifier and spectral sparsifier and introduced

More information

arxiv: v2 [math.co] 8 Mar 2018

arxiv: v2 [math.co] 8 Mar 2018 Restricted lonesu atrices arxiv:1711.10178v2 [ath.co] 8 Mar 2018 Beáta Bényi Faculty of Water Sciences, National University of Public Service, Budapest beata.benyi@gail.co March 9, 2018 Keywords: enueration,

More information

Acyclic Colorings of Directed Graphs

Acyclic Colorings of Directed Graphs Acyclic Colorings of Directed Graphs Noah Golowich Septeber 9, 014 arxiv:1409.7535v1 [ath.co] 6 Sep 014 Abstract The acyclic chroatic nuber of a directed graph D, denoted χ A (D), is the iniu positive

More information

A method to determine relative stroke detection efficiencies from multiplicity distributions

A method to determine relative stroke detection efficiencies from multiplicity distributions A ethod to deterine relative stroke detection eiciencies ro ultiplicity distributions Schulz W. and Cuins K. 2. Austrian Lightning Detection and Inoration Syste (ALDIS), Kahlenberger Str.2A, 90 Vienna,

More information

Left-to-right maxima in words and multiset permutations

Left-to-right maxima in words and multiset permutations Left-to-right axia in words and ultiset perutations Ay N. Myers Saint Joseph s University Philadelphia, PA 19131 Herbert S. Wilf University of Pennsylvania Philadelphia, PA 19104

More information

Hermite s Rule Surpasses Simpson s: in Mathematics Curricula Simpson s Rule. Should be Replaced by Hermite s

Hermite s Rule Surpasses Simpson s: in Mathematics Curricula Simpson s Rule. Should be Replaced by Hermite s International Matheatical Foru, 4, 9, no. 34, 663-686 Herite s Rule Surpasses Sipson s: in Matheatics Curricula Sipson s Rule Should be Replaced by Herite s Vito Lapret University of Lublana Faculty of

More information

lecture 36: Linear Multistep Mehods: Zero Stability

lecture 36: Linear Multistep Mehods: Zero Stability 95 lecture 36: Linear Multistep Mehods: Zero Stability 5.6 Linear ultistep ethods: zero stability Does consistency iply convergence for linear ultistep ethods? This is always the case for one-step ethods,

More information

Probability and Stochastic Processes: A Friendly Introduction for Electrical and Computer Engineers Roy D. Yates and David J.

Probability and Stochastic Processes: A Friendly Introduction for Electrical and Computer Engineers Roy D. Yates and David J. Probability and Stochastic Processes: A Friendly Introduction for Electrical and oputer Engineers Roy D. Yates and David J. Goodan Proble Solutions : Yates and Goodan,1..3 1.3.1 1.4.6 1.4.7 1.4.8 1..6

More information

A Simple Regression Problem

A Simple Regression Problem A Siple Regression Proble R. M. Castro March 23, 2 In this brief note a siple regression proble will be introduced, illustrating clearly the bias-variance tradeoff. Let Y i f(x i ) + W i, i,..., n, where

More information

NORMAL MATRIX POLYNOMIALS WITH NONSINGULAR LEADING COEFFICIENTS

NORMAL MATRIX POLYNOMIALS WITH NONSINGULAR LEADING COEFFICIENTS NORMAL MATRIX POLYNOMIALS WITH NONSINGULAR LEADING COEFFICIENTS NIKOLAOS PAPATHANASIOU AND PANAYIOTIS PSARRAKOS Abstract. In this paper, we introduce the notions of weakly noral and noral atrix polynoials,

More information

Lecture 21 Principle of Inclusion and Exclusion

Lecture 21 Principle of Inclusion and Exclusion Lecture 21 Principle of Inclusion and Exclusion Holden Lee and Yoni Miller 5/6/11 1 Introduction and first exaples We start off with an exaple Exaple 11: At Sunnydale High School there are 28 students

More information

Supplementary Information for Design of Bending Multi-Layer Electroactive Polymer Actuators

Supplementary Information for Design of Bending Multi-Layer Electroactive Polymer Actuators Suppleentary Inforation for Design of Bending Multi-Layer Electroactive Polyer Actuators Bavani Balakrisnan, Alek Nacev, and Elisabeth Sela University of Maryland, College Park, Maryland 074 1 Analytical

More information

Note on generating all subsets of a finite set with disjoint unions

Note on generating all subsets of a finite set with disjoint unions Note on generating all subsets of a finite set with disjoint unions David Ellis e-ail: dce27@ca.ac.uk Subitted: Dec 2, 2008; Accepted: May 12, 2009; Published: May 20, 2009 Matheatics Subject Classification:

More information

ADVANCES ON THE BESSIS- MOUSSA-VILLANI TRACE CONJECTURE

ADVANCES ON THE BESSIS- MOUSSA-VILLANI TRACE CONJECTURE ADVANCES ON THE BESSIS- MOUSSA-VILLANI TRACE CONJECTURE CHRISTOPHER J. HILLAR Abstract. A long-standing conjecture asserts that the polynoial p(t = Tr(A + tb ] has nonnegative coefficients whenever is

More information

arxiv: v1 [cs.ds] 29 Jan 2012

arxiv: v1 [cs.ds] 29 Jan 2012 A parallel approxiation algorith for ixed packing covering seidefinite progras arxiv:1201.6090v1 [cs.ds] 29 Jan 2012 Rahul Jain National U. Singapore January 28, 2012 Abstract Penghui Yao National U. Singapore

More information

Bootstrapping Dependent Data

Bootstrapping Dependent Data Bootstrapping Dependent Data One of the key issues confronting bootstrap resapling approxiations is how to deal with dependent data. Consider a sequence fx t g n t= of dependent rando variables. Clearly

More information

Stochastic vertex models and symmetric functions

Stochastic vertex models and symmetric functions Stochastic vertex odels and syetric functions Alexey Bufetov MIT 6Noveber,2017 Applications of the algebra of syetric functions to probability: Schur easures (Okounkov), Schur processes (Okounkov-Reshetikhin),

More information

Lectures 8 & 9: The Z-transform.

Lectures 8 & 9: The Z-transform. Lectures 8 & 9: The Z-transfor. 1. Definitions. The Z-transfor is defined as a function series (a series in which each ter is a function of one or ore variables: Z[] where is a C valued function f : N

More information

Supplementary Material for Fast and Provable Algorithms for Spectrally Sparse Signal Reconstruction via Low-Rank Hankel Matrix Completion

Supplementary Material for Fast and Provable Algorithms for Spectrally Sparse Signal Reconstruction via Low-Rank Hankel Matrix Completion Suppleentary Material for Fast and Provable Algoriths for Spectrally Sparse Signal Reconstruction via Low-Ran Hanel Matrix Copletion Jian-Feng Cai Tianing Wang Ke Wei March 1, 017 Abstract We establish

More information

Research Article On the Isolated Vertices and Connectivity in Random Intersection Graphs

Research Article On the Isolated Vertices and Connectivity in Random Intersection Graphs International Cobinatorics Volue 2011, Article ID 872703, 9 pages doi:10.1155/2011/872703 Research Article On the Isolated Vertices and Connectivity in Rando Intersection Graphs Yilun Shang Institute for

More information

An Extension to the Tactical Planning Model for a Job Shop: Continuous-Time Control

An Extension to the Tactical Planning Model for a Job Shop: Continuous-Time Control An Extension to the Tactical Planning Model for a Job Shop: Continuous-Tie Control Chee Chong. Teo, Rohit Bhatnagar, and Stephen C. Graves Singapore-MIT Alliance, Nanyang Technological Univ., and Massachusetts

More information

Soft Computing Techniques Help Assign Weights to Different Factors in Vulnerability Analysis

Soft Computing Techniques Help Assign Weights to Different Factors in Vulnerability Analysis Soft Coputing Techniques Help Assign Weights to Different Factors in Vulnerability Analysis Beverly Rivera 1,2, Irbis Gallegos 1, and Vladik Kreinovich 2 1 Regional Cyber and Energy Security Center RCES

More information

Closed-form evaluations of Fibonacci Lucas reciprocal sums with three factors

Closed-form evaluations of Fibonacci Lucas reciprocal sums with three factors Notes on Nuber Theory Discrete Matheatics Print ISSN 30-32 Online ISSN 2367-827 Vol. 23 207 No. 2 04 6 Closed-for evaluations of Fibonacci Lucas reciprocal sus with three factors Robert Frontczak Lesbank

More information

Generalized eigenfunctions and a Borel Theorem on the Sierpinski Gasket.

Generalized eigenfunctions and a Borel Theorem on the Sierpinski Gasket. Generalized eigenfunctions and a Borel Theore on the Sierpinski Gasket. Kasso A. Okoudjou, Luke G. Rogers, and Robert S. Strichartz May 26, 2006 1 Introduction There is a well developed theory (see [5,

More information

A Simplified Analytical Approach for Efficiency Evaluation of the Weaving Machines with Automatic Filling Repair

A Simplified Analytical Approach for Efficiency Evaluation of the Weaving Machines with Automatic Filling Repair Proceedings of the 6th SEAS International Conference on Siulation, Modelling and Optiization, Lisbon, Portugal, Septeber -4, 006 0 A Siplified Analytical Approach for Efficiency Evaluation of the eaving

More information

Midterm 1 Sample Solution

Midterm 1 Sample Solution Midter 1 Saple Solution NOTE: Throughout the exa a siple graph is an undirected, unweighted graph with no ultiple edges (i.e., no exact repeats of the sae edge) and no self-loops (i.e., no edges fro a

More information

Lecture October 23. Scribes: Ruixin Qiang and Alana Shine

Lecture October 23. Scribes: Ruixin Qiang and Alana Shine CSCI699: Topics in Learning and Gae Theory Lecture October 23 Lecturer: Ilias Scribes: Ruixin Qiang and Alana Shine Today s topic is auction with saples. 1 Introduction to auctions Definition 1. In a single

More information

1 Generalization bounds based on Rademacher complexity

1 Generalization bounds based on Rademacher complexity COS 5: Theoretical Machine Learning Lecturer: Rob Schapire Lecture #0 Scribe: Suqi Liu March 07, 08 Last tie we started proving this very general result about how quickly the epirical average converges

More information

Algebraic Montgomery-Yang problem: the log del Pezzo surface case

Algebraic Montgomery-Yang problem: the log del Pezzo surface case c 2014 The Matheatical Society of Japan J. Math. Soc. Japan Vol. 66, No. 4 (2014) pp. 1073 1089 doi: 10.2969/jsj/06641073 Algebraic Montgoery-Yang proble: the log del Pezzo surface case By DongSeon Hwang

More information

Fast Montgomery-like Square Root Computation over GF(2 m ) for All Trinomials

Fast Montgomery-like Square Root Computation over GF(2 m ) for All Trinomials Fast Montgoery-like Square Root Coputation over GF( ) for All Trinoials Yin Li a, Yu Zhang a, a Departent of Coputer Science and Technology, Xinyang Noral University, Henan, P.R.China Abstract This letter

More information

THE POLYNOMIAL REPRESENTATION OF THE TYPE A n 1 RATIONAL CHEREDNIK ALGEBRA IN CHARACTERISTIC p n

THE POLYNOMIAL REPRESENTATION OF THE TYPE A n 1 RATIONAL CHEREDNIK ALGEBRA IN CHARACTERISTIC p n THE POLYNOMIAL REPRESENTATION OF THE TYPE A n RATIONAL CHEREDNIK ALGEBRA IN CHARACTERISTIC p n SHEELA DEVADAS AND YI SUN Abstract. We study the polynoial representation of the rational Cherednik algebra

More information

Uniform Approximation and Bernstein Polynomials with Coefficients in the Unit Interval

Uniform Approximation and Bernstein Polynomials with Coefficients in the Unit Interval Unifor Approxiation and Bernstein Polynoials with Coefficients in the Unit Interval Weiang Qian and Marc D. Riedel Electrical and Coputer Engineering, University of Minnesota 200 Union St. S.E. Minneapolis,

More information

The proofs of Theorem 1-3 are along the lines of Wied and Galeano (2013).

The proofs of Theorem 1-3 are along the lines of Wied and Galeano (2013). A Appendix: Proofs The proofs of Theore 1-3 are along the lines of Wied and Galeano (2013) Proof of Theore 1 Let D[d 1, d 2 ] be the space of càdlàg functions on the interval [d 1, d 2 ] equipped with

More information

Ocean 420 Physical Processes in the Ocean Project 1: Hydrostatic Balance, Advection and Diffusion Answers

Ocean 420 Physical Processes in the Ocean Project 1: Hydrostatic Balance, Advection and Diffusion Answers Ocean 40 Physical Processes in the Ocean Project 1: Hydrostatic Balance, Advection and Diffusion Answers 1. Hydrostatic Balance a) Set all of the levels on one of the coluns to the lowest possible density.

More information

What are primes in graphs and how many of them have a given length?

What are primes in graphs and how many of them have a given length? What are pries in graphs an how any of the have a given length? Aurey Terras Math. Club Oct. 30, 2008 A graph is a bunch of vertices connecte by eges. The siplest exaple for the talk is the tetraheron

More information

Multi-Dimensional Hegselmann-Krause Dynamics

Multi-Dimensional Hegselmann-Krause Dynamics Multi-Diensional Hegselann-Krause Dynaics A. Nedić Industrial and Enterprise Systes Engineering Dept. University of Illinois Urbana, IL 680 angelia@illinois.edu B. Touri Coordinated Science Laboratory

More information

G G G G G. Spec k G. G Spec k G G. G G m G. G Spec k. Spec k

G G G G G. Spec k G. G Spec k G G. G G m G. G Spec k. Spec k 12 VICTORIA HOSKINS 3. Algebraic group actions and quotients In this section we consider group actions on algebraic varieties and also describe what type of quotients we would like to have for such group

More information

CS Lecture 13. More Maximum Likelihood

CS Lecture 13. More Maximum Likelihood CS 6347 Lecture 13 More Maxiu Likelihood Recap Last tie: Introduction to axiu likelihood estiation MLE for Bayesian networks Optial CPTs correspond to epirical counts Today: MLE for CRFs 2 Maxiu Likelihood

More information

Generalized AOR Method for Solving System of Linear Equations. Davod Khojasteh Salkuyeh. Department of Mathematics, University of Mohaghegh Ardabili,

Generalized AOR Method for Solving System of Linear Equations. Davod Khojasteh Salkuyeh. Department of Mathematics, University of Mohaghegh Ardabili, Australian Journal of Basic and Applied Sciences, 5(3): 35-358, 20 ISSN 99-878 Generalized AOR Method for Solving Syste of Linear Equations Davod Khojasteh Salkuyeh Departent of Matheatics, University

More information

Linear recurrences and asymptotic behavior of exponential sums of symmetric boolean functions

Linear recurrences and asymptotic behavior of exponential sums of symmetric boolean functions Linear recurrences and asyptotic behavior of exponential sus of syetric boolean functions Francis N. Castro Departent of Matheatics University of Puerto Rico, San Juan, PR 00931 francis.castro@upr.edu

More information

Characterization of the Line Complexity of Cellular Automata Generated by Polynomial Transition Rules. Bertrand Stone

Characterization of the Line Complexity of Cellular Automata Generated by Polynomial Transition Rules. Bertrand Stone Characterization of the Line Coplexity of Cellular Autoata Generated by Polynoial Transition Rules Bertrand Stone Abstract Cellular autoata are discrete dynaical systes which consist of changing patterns

More information

Elastic Force: A Force Balance: Elastic & Gravitational Force: Force Example: Determining Spring Constant. Some Other Forces

Elastic Force: A Force Balance: Elastic & Gravitational Force: Force Example: Determining Spring Constant. Some Other Forces Energy Balance, Units & Proble Solving: Mechanical Energy Balance ABET Course Outcoes: 1. solve and docuent the solution of probles involving eleents or configurations not previously encountered (e) (e.g.

More information

Birthday Paradox Calculations and Approximation

Birthday Paradox Calculations and Approximation Birthday Paradox Calculations and Approxiation Joshua E. Hill InfoGard Laboratories -March- v. Birthday Proble In the birthday proble, we have a group of n randoly selected people. If we assue that birthdays

More information

COS 424: Interacting with Data. Written Exercises

COS 424: Interacting with Data. Written Exercises COS 424: Interacting with Data Hoework #4 Spring 2007 Regression Due: Wednesday, April 18 Written Exercises See the course website for iportant inforation about collaboration and late policies, as well

More information

MODULAR HYPERBOLAS AND THE CONGRUENCE ax 1 x 2 x k + bx k+1 x k+2 x 2k c (mod m)

MODULAR HYPERBOLAS AND THE CONGRUENCE ax 1 x 2 x k + bx k+1 x k+2 x 2k c (mod m) #A37 INTEGERS 8 (208) MODULAR HYPERBOLAS AND THE CONGRUENCE ax x 2 x k + bx k+ x k+2 x 2k c (od ) Anwar Ayyad Departent of Matheatics, Al Azhar University, Gaza Strip, Palestine anwarayyad@yahoo.co Todd

More information

THE SUPER CATALAN NUMBERS S(m, m + s) FOR s 3 AND SOME INTEGER FACTORIAL RATIOS. 1. Introduction. = (2n)!

THE SUPER CATALAN NUMBERS S(m, m + s) FOR s 3 AND SOME INTEGER FACTORIAL RATIOS. 1. Introduction. = (2n)! THE SUPER CATALAN NUMBERS S(, + s FOR s 3 AND SOME INTEGER FACTORIAL RATIOS XIN CHEN AND JANE WANG Abstract. We give a cobinatorial interpretation for the super Catalan nuber S(, + s for s 3 using lattice

More information

Fixed-to-Variable Length Distribution Matching

Fixed-to-Variable Length Distribution Matching Fixed-to-Variable Length Distribution Matching Rana Ali Ajad and Georg Böcherer Institute for Counications Engineering Technische Universität München, Gerany Eail: raa2463@gail.co,georg.boecherer@tu.de

More information

Extension of CSRSM for the Parametric Study of the Face Stability of Pressurized Tunnels

Extension of CSRSM for the Parametric Study of the Face Stability of Pressurized Tunnels Extension of CSRSM for the Paraetric Study of the Face Stability of Pressurized Tunnels Guilhe Mollon 1, Daniel Dias 2, and Abdul-Haid Soubra 3, M.ASCE 1 LGCIE, INSA Lyon, Université de Lyon, Doaine scientifique

More information

ON THE TWO-LEVEL PRECONDITIONING IN LEAST SQUARES METHOD

ON THE TWO-LEVEL PRECONDITIONING IN LEAST SQUARES METHOD PROCEEDINGS OF THE YEREVAN STATE UNIVERSITY Physical and Matheatical Sciences 04,, p. 7 5 ON THE TWO-LEVEL PRECONDITIONING IN LEAST SQUARES METHOD M a t h e a t i c s Yu. A. HAKOPIAN, R. Z. HOVHANNISYAN

More information

The Distribution of the Covariance Matrix for a Subset of Elliptical Distributions with Extension to Two Kurtosis Parameters

The Distribution of the Covariance Matrix for a Subset of Elliptical Distributions with Extension to Two Kurtosis Parameters journal of ultivariate analysis 58, 96106 (1996) article no. 0041 The Distribution of the Covariance Matrix for a Subset of Elliptical Distributions with Extension to Two Kurtosis Paraeters H. S. Steyn

More information

E0 370 Statistical Learning Theory Lecture 6 (Aug 30, 2011) Margin Analysis

E0 370 Statistical Learning Theory Lecture 6 (Aug 30, 2011) Margin Analysis E0 370 tatistical Learning Theory Lecture 6 (Aug 30, 20) Margin Analysis Lecturer: hivani Agarwal cribe: Narasihan R Introduction In the last few lectures we have seen how to obtain high confidence bounds

More information

Pseudo-marginal Metropolis-Hastings: a simple explanation and (partial) review of theory

Pseudo-marginal Metropolis-Hastings: a simple explanation and (partial) review of theory Pseudo-arginal Metropolis-Hastings: a siple explanation and (partial) review of theory Chris Sherlock Motivation Iagine a stochastic process V which arises fro soe distribution with density p(v θ ). Iagine

More information