Matrix Multiplication II

Size: px
Start display at page:

Download "Matrix Multiplication II"

Transcription

1 Matrx Multplcaton II Yuval Flmu January 14, 201 Thee note tarted ther lfe a a lecture gven at the Toronto Student Semnar on February 9, The materal taken motly from the clac paper by Coppermth and Wnograd [CW]. Other ource are 15.7 of Algebrac Complexty Theory [ACT], Stother the [Sto], V. Wllam recent paper [Wl], and the paper by Cohn at al. [CKSU]. Starred ecton are the one we ddn t have tme to cover. We preent three dfferent algorthm, all taken from [CW], n rapd ucceon. All thee algorthm are baed on Straen groundbreakng laer method. Straen orgnal dea are decrbed n the appendx. 1 Algorthm 1 The tartng pont the dentty ɛ ɛ ( ( ) x [0] 0 y[1] z [1] + x [1] y [0] 0 z[1] + x [1] y [1] z [0] 0 + O(ɛ 4 ) = (x [0] 0 + ɛx[1] )(y [0] 0 + ɛy[1] )(z [0] 0 + ɛz[1] ) x [0] 0 + ɛ2 x [1] (1 qɛ)x [0] 0 y[0] 0 z[0] 0. ) ( y [0] 0 + ɛ2 y [1] ) ( z [0] 0 + ɛ2 On the left-hand de, we have a um of three matrx product 1, 1, q, q, 1, 1, 1, q, 1 harng ome varable. We can ndcate thee hared varable ung upercrpt: z [1] ɛ ( 1, 1, q 0,1,1 + q, 1, 1 1,0,1 + 1, q, 1 1,1,0 ) + O(ɛ 4 ). The notaton n, m, p,j,k ndcate a tenor of type n, m, p whoe x-varable have upercrpt, whoe y-varable have upercrpt j, and whoe z-varable have upercrpt k. The pont that dfferent upercrpt of x correpond to djont et of varable. The dentty a t tand n t very ueful, nce the three (rather trval) matrx product nvolved all hare varable. The dea of the laer method to compute the Nth tenor power of the dentty, and then omehow eparate the varable. The N th tenor power of the dentty ha (q + 2) N term on the rght-hand de. What about the left-hand de? It of the form ) + 1

2 ɛ N F + O(ɛ N+1 ), where F now the um of N dfferent matrx product tenor, each of volume q N (the volume of a tenor n, m, p nmp). The plan now to take a large ubet of the varable o that f we zero out all the other varable, F eparate nto a djont um. More pecfcally, we can thnk of F a a et of trple of ndce. We wll ngle out large et of ndce X, Y, Z uch that among the trple F X Y Z, no varable repeat. In other word, for each x X, there at mot one par (y, z) Y Z uch that (x, y, z) F. How many trple can we expect to obtan n th way? Let P be the ze of the projecton of F nto one of the coordnate (they are all the ame by ymmetry). We defntely cannot expect to get more than P trple. In our cae, P = 2 N. However, th bound a tad nave, a the followng argument from [CKSU] how. Clafy the trple accordng to ther ource dtrbuton, that, how many trple of each type 110, 101, 011 generated them. If there are ( a, b, c of thoe, repectvely (a + b + c = N), then the projecton on the coordnate have ze N ) ( c, N ) ( b, N ) a. Hence the number of trple of type a, b, c at mot the mnmum of thee, whch at mot ( ) N N/. There are O(N 2 ) type, hence there can be at mot O(N 2 ) ( N N/) 2 h(1/)n trple, where h(p) = p log 2 p (1 p) log 2 (1 p) the bnary entropy functon. What th argument teache u that t enough to conder one type of trple, n th cae N/, N/, N/. Techncally, th wll manfet telf n the followng way. For a ubet G of F, contruct a conflct graph by takng the et G a vertce, and connectng any two vertce whch hare a coordnate. The ubet G contng of all et of trple of type N/, N/, N/ ha the property that all vertce have the ame degree. Converely, the degree n the graph correpondng to F depend on the ource dtrbuton: a trple of type a, b, c ha roughly 2 N a + 2 N b + 2 N c neghbor. To ee th, conder for example the frt coordnate, whch cont of c zeroe and N c one. Each ndex equal to one could come from ether 110 or 101, whle each ndex equal to zero certanly came from 011. Surprngly, th method of reducng F to G by pecfyng a ource dtrbuton optmal: we can zero varable o that G eparate to a djont um of P 1 o(1) term, where P the projecton of G to each of the coordnate. We wll tate th hortly, but frt, a defnton. Defnton 1. Let T Z be a fnte collecton of trple, and let φ: T R + be arbtrary. We defne below the numercal quantty cap(t, φ), the capacty of T wth repect to φ. Let N 1. We extend φ to T N multplcatvely, and to ubet of T N addtvely. For {1, 2, }, let α : T N Z N denote the projecton to the th coordnate. Say that ubet X α (T N ) ( {1, 2, }) are good f α njectve on T N X 1 X 2 X for {1, 2, }. The N-capacty cap N (T, φ) the maxmum of φ(t N X 1 X 2 X ) over all good ubet X 1, X 2, X. The capacty cap(t, φ) defned a cap(t, φ) = lm up N N cap N (T, φ). For example, n our cae T = {(1, 1, 0), (1, 0, 1), (0, 1, 1)}, and we can take φ to be contant. In the next ecton, we wll ee an example n whch φ ha to be non-contant. We defned the capacty a a lmt uperor, but a tenor product contructon how that cap N1 +N 2 (T, φ) N1 (T, φ) N2 (T, φ), and th mple that the lmt actually ext. Lemma 1. Let T, φ be a n Defnton 1. Suppoe that furthermore, each trple n T um to a contant, that T cloed under permutaton, and that φ contant over permutaton. Then 2

3 cap(t, φ) atfe the followng nequalte: log 2 cap(t, φ) max H(π 1) + π(t) log π 2 φ(t), t T log 2 cap(t, φ) max H(π) max H(τ) + H(π 1 ) + π(t) log π τ : τ 1 =π 2 φ(t), 1 t T where π range over all permutaton-nvarant probablty dtrbuton over T, π 1 the dtrbuton of the frt coordnate under π, and H the entropy functon gven by H(π) = t T π(t) log 2 π(t). Note that the maxmum ext nce the range of π compact. In our preent cae, T permutaton-nvarant, and o the lemma how that cap(t, 1) = 2 H(1/,2/) = 2 h(1/). We are now n poton to apply Schönhage aymptotc um nequalty. Takng the N th tenor power of the dentty and zerong out varable approprately, we get on the left-hand ze a djont um of roughly 2 h(1/)n term of volume q N. Hence 2 h(1/)n q τn (q + 2) N. Rearrangng varable, we deduce q + 2 ω log q 2 h(1/). Pluggng q = 8, we get ω The contructon n Lemma 1 randomzed. In other word, the fnal algorthm n t explct. Th make no dfference n our model, nce t non-unform. The method of condtonal expectaton can determnze the lemma, n cae one really care about unformty, though fndng the algorthm, whle t take polynomal tme, mght take more tme than executng t. Tghtne of Lemma 1 When the range of T mall, τ = π. In partcular, th the cae (a can be checked ung lnear algebra) f all trple n T are non-negatve and um to a contant whch at mot 5. For a um of 6, we have the followng counterexample: X = S(0, 2, 4), S(1, 2, ) 2, Y = S(0,, ) 2, S(1, 1, 4) 2, S(2, 2, 2) 6. Here S(a, b, c) a horthand for all dtnct permutaton of (a, b, c). Both X and Y have the ame projecton nto each of the coordnate, but the conttuent are dfferent. We have H(X) = H( 1 18 [6], 1 9 [6] ).50, H(Y ) = H( 1 [6] 9, 1 ) Coppermth and Wnograd [CW], followed by Wllam [Wl], avod th ue by retrctng themelve to dtrbuton π for whch τ = π. That, ther dtrbuton are pecfed only by the margnal. A we remarked, n the former cae, th doen t matter. However, n the latter cae, t could matter. Indeed, the N th power of the algorthm contan all non-negatve trple (a, b, c) wth a + b + c = 2N. Th property follow from the fact that t hold for the bac algorthm.

4 1.1 Proof of Lemma 1 Before delvng nto the proof telf, we need ome auxlary reult. Lemma 2. Let π be a partton of n 2 nto k part. Then there a contant C k uch that ( n π) n C k 2 H(π/n)n nc k. Proof. Wthout lo of generalty, all part n π are non-empty. Accordng to Wkpeda, there ext contant K 1, K 2 uch that n! K 1 n n+1/2 e n K 2. Therefore ( n K 1 K2 k π) A K 2 K1 k, where A = nn+1/2 e n p π pp+1/2 e p n = p π p (n/p) p p π n =. p π p2h(π/n)n We can etmate Therefore K 1 K2 k n 1 k n (1 k)/2 A 2 H(π/n)n n. ( n π) 2 H(π/n)n K 2 K1 k n 1/2. Lemma. Let π be a fnte probablty dtrbuton. There ext a contant C uch that for all n, there a fnte probablty dtrbuton σ uch that nσ ntegral, σ(t) π(t) 1/n and H(π) H(σ) C/n. Proof. Let S be the upport of π. Conder σ L and σ H defned by σ L (t) = nπ(t) /n and σ H (t) = nπ(t) /n. Clearly σ L 1 whle σ H 1. It not hard to conclude that there a probablty dtrbuton σ uch that nσ ntegral and σ(t) π(t) 1/n. Now H(π)/ π(t) = log e(π(t)+1). In partcular, H(π) 1 = C contant. It follow that H(π) H(σ) C/n. We dvde the proof nto two part: upper bound and lower bound. 4

5 Upper bound Fx N 1, let X 1, X 2, X be good ubet of T N, and let F = T N X 1 X 2 X. Let D denote the et of all ource dtrbuton, that the et of all functon d: T N that um to N. Note D = poly(n). For each d D, let F d F cont of thoe element conformng to the dtrbuton d. Alo, let d denote the projecton dtrbuton on the th component. Snce the projecton are njectve on F, ung Lemma 2 we get ( ) F d mn N d mn 2H(d /N)N poly(n). Therefore φ(f ) = d D φ(t) d(t) F d t T poly(n) max d D t T φ(t) d(t) mn {1,2,} 2H(d /N)N. Now let be the et of all probablty dtrbuton on T. For any d D, d/n, and o Takng Nth root, we deduce that cap N (T, φ) poly(n) max π cap(t, φ) max π t T t T φ(t) π(t) φ(t) π(t)n mn {1,2,} 2H(p )N. mn {1,2,} 2H(π )N. Gven a dtrbuton π, let S(π) denote t ymmetrzaton. Note that S(π) = (π 1 +π 2 +π )/. Snce the entropy functon concave, H(S(π) 1 ) H(π 1) + H(π 2 ) + H(π ) On the other hand, nce φ permutaton-nvarant, φ(t) S(π)(t). t T φ(t) π(t) = t T mn(h(π 1 ), H(π 2 ), H(π )). Therefore the maxmum obtaned at ome ymmetrc dtrbuton. Lower bound Th the dffcult part of the proof. We wll ue Salem-Spencer et, whch are ubet of the nteger wthout three-term arthmetc progreon. Salem and Spencer [SS] howed contructvely that for all M, there ext uch ubet of {1,..., M} of ze M 1 o(1). Ther contructon wa mproved by Behrend [Beh], Moer [Mo] and Elkn [Elk]. Behrend and Elkn proof are randomzed, whle Moer proof contructve. For u, the orgnal contructon uffce. In fact, we wll employ Salem-Spencer et over the group Z M, for M odd. We can contruct uch et by takng Salem-Spencer ubet of {0,..., (M 1)/2}. Such et tll have ze M 1 o(1). Wthout lo of generalty, aume that all trple n T um to zero. Let π be any dtrbuton on T. Fx N 1, and let F be the ubet of T N correpondng to ource dtrbuton Nσ, where σ gven by Lemma. 5

6 Let M be an odd prme to be determned later. Let h: Z n M Z M be a random lnear functon, and chooe x 1, x 2, x Z M randomly under the contrant x 1 + x 2 + x = 0. Fnally, for x T N defne h (x) = h(α (x)) + x. Our contructon guarantee that f α (x) α (y) and j then (h (x), h (y), h j (x)) a unformly random element n Z p (nce M prme). Fnally, defne h 1 (x) = 2h 1 (x), h 2 (x) = 2h 2(x) and h (x) = h (x); the reaon behnd th werd defnton wll become apparent later on. Snce M odd, (h (x), h (y), h j (x)) alo a unformly random element. Let A be ome Salem-Spencer et for Z M of ze M 1 o(1). Defne et Y 1, Y 2, Y a follow: By contructon, for every x T N, Y = {x α (F ) : h (x) A}. h 1 (x) + h 2 (x) + h (x) = h(α 1 (x) + α 2 (x) + α (x)) = 0. Therefore h 1 (x) + h 2 (x) = 2h (x). Snce A a Salem-Spencer et, we deduce that h 1 (x) = h 2 (x) = h (x). Converely, f h 1 (x) = h 2 (x) then 2h 1 (x) = 2h (x) and o h 1 (x) = h (x). Th a property atfed by every element of T n Y 1 Y 2 Y. Th et therefore partton nto et (V a ) a A, where V a contan thoe member uch that h 1 (x) = h 2 (x) = h (x) = a. We contruct the et X out of the et Y eparately for each a. Gven a A, form a graph whoe vertex et V a, and two vertce are connected by an edge f they conflct, that, they hare one of the three coordnate. Denotng the et of edge by U a, there are at leat V a U a connected component. Chooe a repreentatve from each connected component. The et X contan all th coordnate of repreentatve. The graph wa contructed to guarantee that X 1, X 2, X good, and the number of trple a V a U a. Th contructon work for mple example, n whch F = T N α 1 (F ) α 2 (F ) α (F ); however, th not alway the cae. To refne th contructon, let V a = V a F, and let U a U a cont of thoe conflct touchng V a. When choong repreentatve, chooe a repreentatve from F poble. There are at leat V a U a connected component contanng trple from F (proof: frt add the edge n U a \U a), and o the reultng number of trple from F at leat a V a U a. Let u etmate now the ze of the et V a and U a. We need to defne a et G related to F : G = {(x 1, x 2, x ) : x α (F )}. We alo defne P = α 1 (F ). Each element n x F belong to V a wth probablty 1/M 2, nce h 1 (x), h 2 (x) unform over Z2 M, and h 1 (x) = h 2 (x) mple h 1 (x) = h 2 (x) = h (x). Next, we count the number of potental conflct nvolvng F. For each x α (F ), there are G /P trple havng x a ther th coordnate. Hence the number of conflct nvolvng F bounded by F G /P. Each conflct (x, y) belong to U a wth probablty 1/M, nce (aumng t a 1-conflct) h 1 (x), h 2 (x), h 2 (y) unform over Z M (nce α 1(x) = α 1 (y) and x y mple α 2 (x) α 2 (y)), and o h 1 (x) = h 2 (x) = h 2 (y) = a happen wth probablty 1/M. Hence the expected value of V a U a F F G M 2 P M. Choong M to be a prme of ze roughly 6 G /P, the reultng number of trple at leat M 1 o(1) F o(1) F = M 2M 2 2 G P. 6

7 Ung Lemma 2 and argument mlar to Lemma, we etmate P 2 H(σ 1)N 2 H(π 1)N. Smlarly, F 2 H(π)N and G 2 H(τ)N, where τ maxmze H(τ) under the contrant τ = π for {1, 2, }; the optmum permutaton-nvarant by the concavty of H. Note we etmate G by the larget contrbuton to t nce there are only polynomally many ource dtrbuton. Puttng everythng together, we get the tatement of the lemma. Remark An alternatve contructon, ung Straen orgnal approach and avodng Salem- Spencer et, appear n 15.7 of [ACT]. We dd take Straen connected component argument, whch replace a le elegant one due to Coppermth and Wnograd. Open problem Whch bound correct, the lower bound or the upper bound? 2 Algorthm 2 The dentty we tarted from can be lghtly tweaked to yeld even more: [ ] ( ) ɛ x [0] 0 y[1] z [1] + x [1] y [0] 0 z[1] + x [1] y [1] z [0] 0 + x [0] 0 y[0] 0 z[2] q+1 + x[0] 0 y[2] q+1 z[0] 0 + x[2] q+1 y[0] 0 z[0] 0 ɛ ( (x [0] 0 + ɛx[1] )(y [0] 0 + ɛy[1] )(z [0] 0 + ɛz[1] ) x [0] 0 + ɛ2 x [1] ) ( y [0] 0 + ɛ2 y [1] ) ( z [0] 0 + ɛ2 (1 qɛ)(x [0] 0 + ɛ x [2] q+1 )(y[0] 0 + ɛ y [2] q+1 )(z[0] 0 + ɛ z [2] q+1 ). z [1] ) + + O(ɛ 4 ) = The new dentty ha three more varable x [2] q+1, y[2] q+1, z[2] q+1 and three more term n the left-hand de 1, 1, 1 0,0,2 + 1, 1, 1 0,2,0 + 1, 1, 1 2,0,0. All th at the cot of no new term n the rght-hand de! We follow the analy of Algorthm 1, only now we have one more degree of freedom. Suppoe p gve weght α/ to each of the term of volume q, and weght β/ to each of the term of volume 1 (o α+β = 1). The projecton r gve probablte (α+2β)/, 2α/, β/ to 0, 1, 2, correpondngly. Followng our earler reaonng, Lemma 1 together wth the aymptotc um nequalty yeld ω log q α q H((α+2β)/,2α/,β/). Th tme the optmzaton more dffcult, but can be done numercally (or wth Lagrange multpler), gvng q = 6, α 0.952, β and ω Remark The tranton from Algorthm 1 to Algorthm 2 nvolve addng three new varable. In exactly the ame way, we could add even more varable, addng for example the term x [0] 0 y[0] 0 z[] q+2. However, th doe not ncreae the capacty. Indeed, uppoe the et X 1 contaned two member a, b dfferng only n z q+1 v. z q+2. Every trple n X 1 X 2 X n whch the frt coordnate a ha a doppelgänger n whch the frt coordnate b, o X 1, X 2, X n t good. 7

8 Algorthm Algorthm ue exactly the ame dentty a Algorthm 2, only quared. There an added complcaton, nce not all term are matrx product tenor. We wll ee how to handle th oon, but frt let ee how the dentty quared look lke. If we follow the recpe gven o far, the ndce wll be vector of length 2. Intead, we wll take the um of the two component (Coppermth and Wnograd call th couplng). Th enure that the x-, y- and z-ndce have a contant um..1 Squared dentty 1. The dentty quared ha (q + 2) 2 term on the rght-hand de, and 15 term on the left-hand de, whch break up a follow: (a) term mlar to 1, 1, 1 0,0,4, comng from (b) 6 term mlar to 1, 1, 2q 0,1,, comng from 1, 1, 1 0,0,2 1, 1, 1 0,0,2. 1, 1, q 0,1,1 1, 1, 1 0,0,2 1, 1, 1 0,0,2 1, 1, q 0,1,1. (c) term mlar to 1, 1, q ,2,2, comng from 1, 1, 1 0,2,0 1, 1, 1 0,0,2 1, 1, 1 0,0,2 1, 1, 1 0,2,0 1, 1, q 0,1,1 1, 1, q 0,1,1. (d) term mlar to T 1,1,2 4, comng from 1, q, 1 1,1,0 1, 1, 1 0,0,2 + 1, 1, 1 0,0,2 1, q, 1 1,1,0 + q, 1, 1 1,0,1 1, 1, q 0,1,1 + 1, 1, q 0,1,1 q, 1, 1 1,0,1. The fourth tenor, T 4, unfortunately not a matrx product tenor n telf, but t can be ued to produce matrx product tenor, n ome ene that we wll make prece. We wll have to generalze the aymptotc um nequalty from (n m p ) τ R ( n, m, p ) ω τ (1) to a more general form V τ (T ) R ( ) T ω τ. (2).2 The value of a tenor We frt preent the defnton of V τ, and then gve ome ntuton. The proof of (2) gven n the ubecton below. The operator ρ operate on tenor by rotatng the varable x, y, z to y, z, x. For example, f T = n, m, p = n m p x j y jk z k j=1 k=1 8

9 then ρ(t ), whch obtaned by replacng x y, y z, z x, ρ(t ) = n m p p n m x k y j z jk = x k y jz kj = p, n, m. j=1 k=1 k=1 j=1 Here x k = x k and z kj = z jk denote tranpoton. Next, we defne ymmetrzaton: S(T ) = T ρ(t ) ρ 2 (T ). So S( n, m, p ) = nmp, nmp, nmp. Symmetrzaton what allowed u to prove ω log nmp R( n, m, p ). The Nth tenor power of a tenor T T N = T T (N factor). We ay that a tenor T break up nto T = n, m, p f we can zero ome varable to get T. Th the ame proce we ued n Algorthm 1 and 2. We defne { V τ,n (T ) = max (n m p ) τ : T break up nto } n, m, p. For example, V τ,n ( n, m, p ) = (nmp) nτ. Fnally, V τ (T ) = up V τ,n (T ) 1/N. N So V τ ( n, m, p ) = (nmp) τ. The dea behnd th defnton that t mmedately mple that Indeed, uppoe that V τ,n (T ) 1/N R(T ). So for ome n, m, p, V τ (T ) = R(T ) ω τ. () V τ,n (T ) = (n m p ) τ, and S(T ) N reduce to T = n, m, p. The latter mple that R(T ) R(S(T ) N ) R(T ) N. The aymptotc um nequalty (1) tate that (n m p ) τ R(T ) ω τ. Snce V τ,n (T ) R(T ) N R(T ), the condton roughly hold, hence ω τ. 4 Calculatng the value It tme for u to take a cloer look at T 4. Wth the orgnal ndce, t look lke th: T 4 = 1, q, 1 10,10,02 + 1, q, 1 01,01,20 + q, 1, q 10,01,11 + q, 1, q 01,10,11. 9

10 Let rename the ndce to make thng clearer: T 4 = 1, q, 1,, + 1, q, 1 4,4,4 + q, 1, q,4,5 + q, 1, q 4,,5. The tenor T 4 not ymmetrc. Before ymmetrzng t, we take the Nth power, and only conder term wth (α/2)n, (α/2)n, (β/2)n, (β/2)n factor of each of the bac type, where α + β = 1. Each uch term ha volume q (α+2β)n. The projecton nto the x-ndex or the y-ndex ha N/2 coordnate equal to and N/2 coordnate equal to 4. The projecton nto the z-ndex ha (α/2)n coordnate equal to, (α/2)n coordnate equal to 4, and βn coordnate equal to 5. If we now ymmetrze, the volume of each bac term q (α+2β)n, and the projecton nto each coordnate ha ze ( ) N 2 ( N/2 α 2 N, α 2 N N, βn Ung Lemma 1, we get that the value at leat roughly the Nth root of ( ) N 2 ( N/2 Takng the Nth root, we get α 2 N, α 2 N N, βn ). ) q (α+2β)nτ 2 2N 2 H(α/2,α/2,β)N q (α+2β)nτ. V τ (T 4 ) 4 2 H(α/2,α/2,β) q (α+6β)τ = 4(α/2) α β β q (+β)τ. We can fnd the optmal value of β by ubttutng α = 1 β, dfferentatng the rght-hand de wth repect to β, and equatng to zero. We fnd that the optmal value of α, β are α = 2 q τ + 2, β = qτ q τ + 2. Subttutng th back, we get V τ (T 4 ) 4 1/ q τ (2 + q τ ) 1/. 4.1 Analyzng the algorthm Havng calculated the value of T 4, t tme to go back to the man plan. Prevouly we have optmzed the exact fracton to take of each of the bac term. However we have 15 of thee. In the end, we re only really ntereted n the ze of the projecton nto each coordnate (aumng the projecton ha the ame ze for all coordnate). Snce there are only 5 value of each coordnate, th reduce the number of degree of freedom to 4. In fact, nce n the dentty, the three ndce alway um to 4, the average of all component of the x-ndex hould be 4/. Th take away one more degree of freedom. Suppoe α + β + γ + δ = 1. Take term n whch each term of the frt knd appear (α/)n tme, each term of the econd type appear (β/6)n tme, each term of the thrd type appear (γ/)n tme, and each term of the fourth type appear (δ/)n tme. The normalzed htogram of each coordnate wll be 2α + β + γ : β + 2δ : 2γ + δ : β : α. 10

11 The aymptotc um nequalty equate (q + 2) 2 wth the entropy of th dtrbuton multpled by the value, whch V τ (T 1 ) α V τ (T 2 ) β V τ (T ) γ V τ (T 4 ) δ. We don t need to explctly compute the value of T 1, T 2, T nce they are matrx product tenor. If we optmze over q, α, β, γ, δ, we get that ω Proof* The generalzed form of the aymptotc um nequalty (2) mplctly ued n Coppermth and Wnograd orgnal paper [CW]. It tated explctly n Wllam work [Wl], but not proved there. The treatment below taken from Stother the [Sto], who explctly prove t (we lghtly modfy h notaton). Lemma 4. For every two tenor T 1, T 2, A a conequence, for every tenor T, Proof. Exerce. Lemma 5. For every tenor T, V τ,n (T 1 T 2 ) V τ,n (T 1 )V τ,n (T 2 ). V τ,cn (T ) V τ,n (T ) c. V τ (T ) = lm N V τ,n(t ) 1/N. Proof. We defned V τ (T ) a the upremum of V τ,n (T ) 1/N. Pck any ɛ > 0. For ome N, V τ,n (T ) 1/N V τ (T ) ɛ. Let M = cn + r, where c 1 and 0 r < N. Snce Lemma 4 how that c M = 1 N r MN > 1 N 1 M, V τ,m (T ) 1/M V τ,cn (T ) 1/M V τ,n (T ) c/m (V τ (T ) ɛ)v τ (T ) 1/M. If M large enough, V τ,m (T ) 1/M V τ (T ) 2ɛ. Corollary 6. For every two tenor T 1, T 2, V τ (T 1 T 2 ) V τ (T 1 )V τ (T 2 ). Lemma 7. For every two tenor T 1, T 2, V τ (T 1 T 2 ) V τ (T 1 ) + V τ (T 2 ). 11

12 Proof. Lemma 5 mple that for large N, V τ,n (T 1 ) V τ (T 1 ) N, V τ,n (T 2 ) V τ (T 2 ) N. The Nth tenor power S(T 1 T 2 ) N a djont um of varou term, ncludng S(T 1 ) N 1 S(T 2 ) N 2 for all N 1, N 2 atfyng N 1 + N 2 = N. Hence the defnton of value mple that If N 1 and N 2 are both large, then V τ,n (T 1 T 2 ) V τ,n1 (T 1 )V τ,n2 (T 2 ). V τ,n1 (T 1 )V τ,n2 (T 2 ) V τ (T 1 ) N 1 V τ (T 2 ) N 2. The rght-hand de maxmzed when N 1 /N 2 V τ (T 1 )/V τ (T 2 ), n whch cae V τ (T 1 ) N 1 V τ (T 2 ) N 2 (V τ (T 1 ) + V τ (T 2 )) N N + 1 When N large, both of N 1, N 2 are large, and o V 1/N τ,n (T 1 T 2 ) V τ (T 1 ) + V τ (T 2 ). Corollary 8. The generalzed aymptotc um nequalty (2) hold. 5 General formulaton* The etmate of ω obtaned ung Algorthm a two-tered proce. The frt tep to calculate everal value, and the econd tep to combne thee value and get an etmate of ω. Each tep requre combnng Lemma 1 wth the aymptotc um nequalty. We can phrae the method n the followng lemma. Lemma 9. Let T be a tenor whch a um of tenor {T : I}, where T ha ndce. Gven a probablty dtrbuton p on I, let r 1, r 2, r be the reultng projecton on the three coordnate. Then log 2 V τ (T ) H(r 1) + H(r 2 ) + H(r ) + E p [log 2 V τ (T )]. Proof ketch. Suppoe frt that the tenor T are matrx multplcaton tenor. Ung the method of proof of Lemma 1, we can demontrate the extence of ubet T N S(S) N of ze very roughly 1 N log 2 T N H(r 1) + H(r 2 ) + H(r ). Each ubet T N correpond to a djont um of matrx multplcaton tenor n, m, p uch that. (nmp) τ = V τ (T ) p()n. 12

13 The etmate of V τ (T ) now follow from the formula for V τ,n (T ). For the general cae, the reader can how that the defnton of V τ reman vald f we replace the orgnal defnton of V τ,n wth { V τ,n (T ) = max V τ (T ) : T break up nto } T. S In order to apply the lemma, we need to ue the bae cae tatng that V τ ( n, m, p ) = (nmp) τ. In order to ue the lemma to obtan a bound on ω, we ue the aymptotc um nequalty V ω (T ) R(T ). A an example, conder the tenor T 4 encountered n Algorthm. The ndex et and known value are gven by V τ (T,, ) = V τ (T 4,4,4 ) = q τ, V τ (T,4,5 ) = V τ (T 4,,5 ) = q 2τ. Convexty mple that n order to maxmze the entropy, we need the probablty p to depend only on the value (n th cae, only on the volume of the correpondng matrx multplcaton tenor), hence we gve each of the frt two trple probablty α/2, and each of the latter two trple probablty β/2 (where α + β = 1). We have H(r 1 ) = H(r 2 ) = H(1/2, 1/2) = 1, H(r ) = H(α/2, α/2, β). Regardng volume, we have E p [log 2 V τ (T )] = τ(α + 2β) log 2 q. The lemma now mple log 2 V τ (T 4 ) 2 + H(α/2, α/2, β) + τ(α + 2β) log 2 q. 6 Straen veron* Straen tart wth the dentty ɛ (x [1] y [0] 0 z[1] + x [0] 0 y[1] z [1] ) + O(ɛ 2 ) = (x [0] 0 + ɛx[1] )(y [0] 0 + ɛy[1] )z x [0] 0 y[0] 0 z. On the left, we have T = q, 1, 1 1,0,1 + 1, 1, q 0,1,1. Straen n effect calculate V τ (T ). Rae the dentty to the Nth power. The projecton have ze 2 N, 2 N, 1, o after ymmetrzng all the projecton have ze 4 N. All the ndvdual term have volume q N. Lemma 1 (whch Straen prove wth an explct contructon) how that V τ (T ) 4 1/ q τ. 1

14 Therefore, the aymptotc um nequalty how that ω τ where 4 1/ q τ = q + 1. For q = 5, th gve ω One dfference between Straen account and the later one the technque ued to eparate a tenor nto djont component. Coppermth and Wnograd do th by zerong out ome of the varable. Straen gve a dfferent ɛ-weght to each nput varable, and then obtan an approxmate equaton (an equaton of the type ued to defne border rank). H method obvate the need for Salem-Spencer et. Queton: The Coppermth-Wnograd contructon can be ealy emulated ung a Straen contructon. I the other drecton alo true? Reference [ACT] Peter Bürger, Mchael Clauen and M. Amn Shokrollah, Algebrac Complexty Theory, Sprnger, [Beh] Felx Behrend, On et of nteger whch contan no three term n arthmetc progreon, Proc. Nat. Acad. Sc. 2:1 2, [Elk] Mchael Elkn, An mproved contructon of progreon-free et, Irael J. Math. 184:9 128, [CKSU] Henry Cohn, Robert Klenberg, Baláz Szegedy and Chr Uman, Group-theoretc algorthm for matrx multplcaton, FOCS [CW] Dan Coppermth and Shmuel Wnograd, Matrx multplcaton va arthmetc progreon, J. Symb.Comput. 9(): , [Mo] Leo Moer, On non-averagng et of nteger, Canad. J. Math. 5:245 25, 195. [SS] Raphaël Salem and Donald Spencer, On et of nteger whch contan no three n arthmetc progreon, Proc. Nat. Acad. Sc. 28:561 56, [Sto] Andrew Jame Stother, On the complexty of matrx multplcaton, Ph.D. the (U. of Ednburgh), [Wl] Vrgna Valevka-Wllam, Breakng the Coppermth-Wnograd barrer, STOC

Additional File 1 - Detailed explanation of the expression level CPD

Additional File 1 - Detailed explanation of the expression level CPD Addtonal Fle - Detaled explanaton of the expreon level CPD A mentoned n the man text, the man CPD for the uterng model cont of two ndvdual factor: P( level gen P( level gen P ( level gen 2 (.).. CPD factor

More information

Improvements on Waring s Problem

Improvements on Waring s Problem Improvement on Warng Problem L An-Png Bejng, PR Chna apl@nacom Abtract By a new recurve algorthm for the auxlary equaton, n th paper, we wll gve ome mprovement for Warng problem Keyword: Warng Problem,

More information

Specification -- Assumptions of the Simple Classical Linear Regression Model (CLRM) 1. Introduction

Specification -- Assumptions of the Simple Classical Linear Regression Model (CLRM) 1. Introduction ECONOMICS 35* -- NOTE ECON 35* -- NOTE Specfcaton -- Aumpton of the Smple Clacal Lnear Regreon Model (CLRM). Introducton CLRM tand for the Clacal Lnear Regreon Model. The CLRM alo known a the tandard lnear

More information

Harmonic oscillator approximation

Harmonic oscillator approximation armonc ocllator approxmaton armonc ocllator approxmaton Euaton to be olved We are fndng a mnmum of the functon under the retrcton where W P, P,..., P, Q, Q,..., Q P, P,..., P, Q, Q,..., Q lnwgner functon

More information

Stanford University CS359G: Graph Partitioning and Expanders Handout 4 Luca Trevisan January 13, 2011

Stanford University CS359G: Graph Partitioning and Expanders Handout 4 Luca Trevisan January 13, 2011 Stanford Unversty CS359G: Graph Parttonng and Expanders Handout 4 Luca Trevsan January 3, 0 Lecture 4 In whch we prove the dffcult drecton of Cheeger s nequalty. As n the past lectures, consder an undrected

More information

Chapter 6 The Effect of the GPS Systematic Errors on Deformation Parameters

Chapter 6 The Effect of the GPS Systematic Errors on Deformation Parameters Chapter 6 The Effect of the GPS Sytematc Error on Deformaton Parameter 6.. General Beutler et al., (988) dd the frt comprehenve tudy on the GPS ytematc error. Baed on a geometrc approach and aumng a unform

More information

Scattering of two identical particles in the center-of. of-mass frame. (b)

Scattering of two identical particles in the center-of. of-mass frame. (b) Lecture # November 5 Scatterng of two dentcal partcle Relatvtc Quantum Mechanc: The Klen-Gordon equaton Interpretaton of the Klen-Gordon equaton The Drac equaton Drac repreentaton for the matrce α and

More information

U.C. Berkeley CS294: Spectral Methods and Expanders Handout 8 Luca Trevisan February 17, 2016

U.C. Berkeley CS294: Spectral Methods and Expanders Handout 8 Luca Trevisan February 17, 2016 U.C. Berkeley CS94: Spectral Methods and Expanders Handout 8 Luca Trevsan February 7, 06 Lecture 8: Spectral Algorthms Wrap-up In whch we talk about even more generalzatons of Cheeger s nequaltes, and

More information

Pythagorean triples. Leen Noordzij.

Pythagorean triples. Leen Noordzij. Pythagorean trple. Leen Noordz Dr.l.noordz@leennoordz.nl www.leennoordz.me Content A Roadmap for generatng Pythagorean Trple.... Pythagorean Trple.... 3 Dcuon Concluon.... 5 A Roadmap for generatng Pythagorean

More information

Improvements on Waring s Problem

Improvements on Waring s Problem Imrovement on Warng Problem L An-Png Bejng 85, PR Chna al@nacom Abtract By a new recurve algorthm for the auxlary equaton, n th aer, we wll gve ome mrovement for Warng roblem Keyword: Warng Problem, Hardy-Lttlewood

More information

Chapter 11. Supplemental Text Material. The method of steepest ascent can be derived as follows. Suppose that we have fit a firstorder

Chapter 11. Supplemental Text Material. The method of steepest ascent can be derived as follows. Suppose that we have fit a firstorder S-. The Method of Steepet cent Chapter. Supplemental Text Materal The method of teepet acent can be derved a follow. Suppoe that we have ft a frtorder model y = β + β x and we wh to ue th model to determne

More information

Problem Set 9 Solutions

Problem Set 9 Solutions Desgn and Analyss of Algorthms May 4, 2015 Massachusetts Insttute of Technology 6.046J/18.410J Profs. Erk Demane, Srn Devadas, and Nancy Lynch Problem Set 9 Solutons Problem Set 9 Solutons Ths problem

More information

Small signal analysis

Small signal analysis Small gnal analy. ntroducton Let u conder the crcut hown n Fg., where the nonlnear retor decrbed by the equaton g v havng graphcal repreentaton hown n Fg.. ( G (t G v(t v Fg. Fg. a D current ource wherea

More information

College of Computer & Information Science Fall 2009 Northeastern University 20 October 2009

College of Computer & Information Science Fall 2009 Northeastern University 20 October 2009 College of Computer & Informaton Scence Fall 2009 Northeastern Unversty 20 October 2009 CS7880: Algorthmc Power Tools Scrbe: Jan Wen and Laura Poplawsk Lecture Outlne: Prmal-dual schema Network Desgn:

More information

= z 20 z n. (k 20) + 4 z k = 4

= z 20 z n. (k 20) + 4 z k = 4 Problem Set #7 solutons 7.2.. (a Fnd the coeffcent of z k n (z + z 5 + z 6 + z 7 + 5, k 20. We use the known seres expanson ( n+l ( z l l z n below: (z + z 5 + z 6 + z 7 + 5 (z 5 ( + z + z 2 + z + 5 5

More information

Lecture 10 Support Vector Machines II

Lecture 10 Support Vector Machines II Lecture 10 Support Vector Machnes II 22 February 2016 Taylor B. Arnold Yale Statstcs STAT 365/665 1/28 Notes: Problem 3 s posted and due ths upcomng Frday There was an early bug n the fake-test data; fxed

More information

Maximizing the number of nonnegative subsets

Maximizing the number of nonnegative subsets Maxmzng the number of nonnegatve subsets Noga Alon Hao Huang December 1, 213 Abstract Gven a set of n real numbers, f the sum of elements of every subset of sze larger than k s negatve, what s the maxmum

More information

Assortment Optimization under MNL

Assortment Optimization under MNL Assortment Optmzaton under MNL Haotan Song Aprl 30, 2017 1 Introducton The assortment optmzaton problem ams to fnd the revenue-maxmzng assortment of products to offer when the prces of products are fxed.

More information

The multivariate Gaussian probability density function for random vector X (X 1,,X ) T. diagonal term of, denoted

The multivariate Gaussian probability density function for random vector X (X 1,,X ) T. diagonal term of, denoted Appendx Proof of heorem he multvarate Gauan probablty denty functon for random vector X (X,,X ) px exp / / x x mean and varance equal to the th dagonal term of, denoted he margnal dtrbuton of X Gauan wth

More information

U.C. Berkeley CS294: Beyond Worst-Case Analysis Luca Trevisan September 5, 2017

U.C. Berkeley CS294: Beyond Worst-Case Analysis Luca Trevisan September 5, 2017 U.C. Berkeley CS94: Beyond Worst-Case Analyss Handout 4s Luca Trevsan September 5, 07 Summary of Lecture 4 In whch we ntroduce semdefnte programmng and apply t to Max Cut. Semdefnte Programmng Recall that

More information

Calculation of time complexity (3%)

Calculation of time complexity (3%) Problem 1. (30%) Calculaton of tme complexty (3%) Gven n ctes, usng exhaust search to see every result takes O(n!). Calculaton of tme needed to solve the problem (2%) 40 ctes:40! dfferent tours 40 add

More information

This appendix presents the derivations and proofs omitted from the main text.

This appendix presents the derivations and proofs omitted from the main text. Onlne Appendx A Appendx: Omtted Dervaton and Proof Th appendx preent the dervaton and proof omtted from the man text A Omtted dervaton n Secton Mot of the analy provded n the man text Here, we formally

More information

Difference Equations

Difference Equations Dfference Equatons c Jan Vrbk 1 Bascs Suppose a sequence of numbers, say a 0,a 1,a,a 3,... s defned by a certan general relatonshp between, say, three consecutve values of the sequence, e.g. a + +3a +1

More information

Complete subgraphs in multipartite graphs

Complete subgraphs in multipartite graphs Complete subgraphs n multpartte graphs FLORIAN PFENDER Unverstät Rostock, Insttut für Mathematk D-18057 Rostock, Germany Floran.Pfender@un-rostock.de Abstract Turán s Theorem states that every graph G

More information

Two Approaches to Proving. Goldbach s Conjecture

Two Approaches to Proving. Goldbach s Conjecture Two Approache to Provng Goldbach Conecture By Bernard Farley Adved By Charle Parry May 3 rd 5 A Bref Introducton to Goldbach Conecture In 74 Goldbach made h mot famou contrbuton n mathematc wth the conecture

More information

More metrics on cartesian products

More metrics on cartesian products More metrcs on cartesan products If (X, d ) are metrc spaces for 1 n, then n Secton II4 of the lecture notes we defned three metrcs on X whose underlyng topologes are the product topology The purpose of

More information

Estimation of Finite Population Total under PPS Sampling in Presence of Extra Auxiliary Information

Estimation of Finite Population Total under PPS Sampling in Presence of Extra Auxiliary Information Internatonal Journal of Stattc and Analy. ISSN 2248-9959 Volume 6, Number 1 (2016), pp. 9-16 Reearch Inda Publcaton http://www.rpublcaton.com Etmaton of Fnte Populaton Total under PPS Samplng n Preence

More information

Polynomials. 1 More properties of polynomials

Polynomials. 1 More properties of polynomials Polynomals 1 More propertes of polynomals Recall that, for R a commutatve rng wth unty (as wth all rngs n ths course unless otherwse noted), we defne R[x] to be the set of expressons n =0 a x, where a

More information

Note on EM-training of IBM-model 1

Note on EM-training of IBM-model 1 Note on EM-tranng of IBM-model INF58 Language Technologcal Applcatons, Fall The sldes on ths subject (nf58 6.pdf) ncludng the example seem nsuffcent to gve a good grasp of what s gong on. Hence here are

More information

APPENDIX A Some Linear Algebra

APPENDIX A Some Linear Algebra APPENDIX A Some Lnear Algebra The collecton of m, n matrces A.1 Matrces a 1,1,..., a 1,n A = a m,1,..., a m,n wth real elements a,j s denoted by R m,n. If n = 1 then A s called a column vector. Smlarly,

More information

Generalized Linear Methods

Generalized Linear Methods Generalzed Lnear Methods 1 Introducton In the Ensemble Methods the general dea s that usng a combnaton of several weak learner one could make a better learner. More formally, assume that we have a set

More information

CHAPTER 9 LINEAR MOMENTUM, IMPULSE AND COLLISIONS

CHAPTER 9 LINEAR MOMENTUM, IMPULSE AND COLLISIONS CHAPTER 9 LINEAR MOMENTUM, IMPULSE AND COLLISIONS 103 Phy 1 9.1 Lnear Momentum The prncple o energy conervaton can be ued to olve problem that are harder to olve jut ung Newton law. It ued to decrbe moton

More information

Randić Energy and Randić Estrada Index of a Graph

Randić Energy and Randić Estrada Index of a Graph EUROPEAN JOURNAL OF PURE AND APPLIED MATHEMATICS Vol. 5, No., 202, 88-96 ISSN 307-5543 www.ejpam.com SPECIAL ISSUE FOR THE INTERNATIONAL CONFERENCE ON APPLIED ANALYSIS AND ALGEBRA 29 JUNE -02JULY 20, ISTANBUL

More information

Errors for Linear Systems

Errors for Linear Systems Errors for Lnear Systems When we solve a lnear system Ax b we often do not know A and b exactly, but have only approxmatons  and ˆb avalable. Then the best thng we can do s to solve ˆx ˆb exactly whch

More information

Lecture 12: Discrete Laplacian

Lecture 12: Discrete Laplacian Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly

More information

Root Locus Techniques

Root Locus Techniques Root Locu Technque ELEC 32 Cloed-Loop Control The control nput u t ynthezed baed on the a pror knowledge of the ytem plant, the reference nput r t, and the error gnal, e t The control ytem meaure the output,

More information

Example: (13320, 22140) =? Solution #1: The divisors of are 1, 2, 3, 4, 5, 6, 9, 10, 12, 15, 18, 20, 27, 30, 36, 41,

Example: (13320, 22140) =? Solution #1: The divisors of are 1, 2, 3, 4, 5, 6, 9, 10, 12, 15, 18, 20, 27, 30, 36, 41, The greatest common dvsor of two ntegers a and b (not both zero) s the largest nteger whch s a common factor of both a and b. We denote ths number by gcd(a, b), or smply (a, b) when there s no confuson

More information

Randomness and Computation

Randomness and Computation Randomness and Computaton or, Randomzed Algorthms Mary Cryan School of Informatcs Unversty of Ednburgh RC 208/9) Lecture 0 slde Balls n Bns m balls, n bns, and balls thrown unformly at random nto bns usually

More information

Edge Isoperimetric Inequalities

Edge Isoperimetric Inequalities November 7, 2005 Ross M. Rchardson Edge Isopermetrc Inequaltes 1 Four Questons Recall that n the last lecture we looked at the problem of sopermetrc nequaltes n the hypercube, Q n. Our noton of boundary

More information

Affine transformations and convexity

Affine transformations and convexity Affne transformatons and convexty The purpose of ths document s to prove some basc propertes of affne transformatons nvolvng convex sets. Here are a few onlne references for background nformaton: http://math.ucr.edu/

More information

Finding Dense Subgraphs in G(n, 1/2)

Finding Dense Subgraphs in G(n, 1/2) Fndng Dense Subgraphs n Gn, 1/ Atsh Das Sarma 1, Amt Deshpande, and Rav Kannan 1 Georga Insttute of Technology,atsh@cc.gatech.edu Mcrosoft Research-Bangalore,amtdesh,annan@mcrosoft.com Abstract. Fndng

More information

Team. Outline. Statistics and Art: Sampling, Response Error, Mixed Models, Missing Data, and Inference

Team. Outline. Statistics and Art: Sampling, Response Error, Mixed Models, Missing Data, and Inference Team Stattc and Art: Samplng, Repone Error, Mxed Model, Mng Data, and nference Ed Stanek Unverty of Maachuett- Amhert, USA 9/5/8 9/5/8 Outlne. Example: Doe-repone Model n Toxcology. ow to Predct Realzed

More information

THE CHINESE REMAINDER THEOREM. We should thank the Chinese for their wonderful remainder theorem. Glenn Stevens

THE CHINESE REMAINDER THEOREM. We should thank the Chinese for their wonderful remainder theorem. Glenn Stevens THE CHINESE REMAINDER THEOREM KEITH CONRAD We should thank the Chnese for ther wonderful remander theorem. Glenn Stevens 1. Introducton The Chnese remander theorem says we can unquely solve any par of

More information

Lecture 4: Universal Hash Functions/Streaming Cont d

Lecture 4: Universal Hash Functions/Streaming Cont d CSE 5: Desgn and Analyss of Algorthms I Sprng 06 Lecture 4: Unversal Hash Functons/Streamng Cont d Lecturer: Shayan Oves Gharan Aprl 6th Scrbe: Jacob Schreber Dsclamer: These notes have not been subjected

More information

Eigenvalues of Random Graphs

Eigenvalues of Random Graphs Spectral Graph Theory Lecture 2 Egenvalues of Random Graphs Danel A. Spelman November 4, 202 2. Introducton In ths lecture, we consder a random graph on n vertces n whch each edge s chosen to be n the

More information

Introduction to Algorithms

Introduction to Algorithms Introducton to Algorthms 6.046J/8.40J Lecture 7 Prof. Potr Indyk Data Structures Role of data structures: Encapsulate data Support certan operatons (e.g., INSERT, DELETE, SEARCH) Our focus: effcency of

More information

From Biot-Savart Law to Divergence of B (1)

From Biot-Savart Law to Divergence of B (1) From Bot-Savart Law to Dvergence of B (1) Let s prove that Bot-Savart gves us B (r ) = 0 for an arbtrary current densty. Frst take the dvergence of both sdes of Bot-Savart. The dervatve s wth respect to

More information

MULTIPLE REGRESSION ANALYSIS For the Case of Two Regressors

MULTIPLE REGRESSION ANALYSIS For the Case of Two Regressors MULTIPLE REGRESSION ANALYSIS For the Cae of Two Regreor In the followng note, leat-quare etmaton developed for multple regreon problem wth two eplanator varable, here called regreor (uch a n the Fat Food

More information

Section 8.3 Polar Form of Complex Numbers

Section 8.3 Polar Form of Complex Numbers 80 Chapter 8 Secton 8 Polar Form of Complex Numbers From prevous classes, you may have encountered magnary numbers the square roots of negatve numbers and, more generally, complex numbers whch are the

More information

NUMERICAL DIFFERENTIATION

NUMERICAL DIFFERENTIATION NUMERICAL DIFFERENTIATION 1 Introducton Dfferentaton s a method to compute the rate at whch a dependent output y changes wth respect to the change n the ndependent nput x. Ths rate of change s called the

More information

MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/15.070J Fall 2013 Lecture 12 10/21/2013. Martingale Concentration Inequalities and Applications

MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/15.070J Fall 2013 Lecture 12 10/21/2013. Martingale Concentration Inequalities and Applications MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.65/15.070J Fall 013 Lecture 1 10/1/013 Martngale Concentraton Inequaltes and Applcatons Content. 1. Exponental concentraton for martngales wth bounded ncrements.

More information

PHYS 705: Classical Mechanics. Calculus of Variations II

PHYS 705: Classical Mechanics. Calculus of Variations II 1 PHYS 705: Classcal Mechancs Calculus of Varatons II 2 Calculus of Varatons: Generalzaton (no constrant yet) Suppose now that F depends on several dependent varables : We need to fnd such that has a statonary

More information

1 The Mistake Bound Model

1 The Mistake Bound Model 5-850: Advanced Algorthms CMU, Sprng 07 Lecture #: Onlne Learnng and Multplcatve Weghts February 7, 07 Lecturer: Anupam Gupta Scrbe: Bryan Lee,Albert Gu, Eugene Cho he Mstake Bound Model Suppose there

More information

NP-Completeness : Proofs

NP-Completeness : Proofs NP-Completeness : Proofs Proof Methods A method to show a decson problem Π NP-complete s as follows. (1) Show Π NP. (2) Choose an NP-complete problem Π. (3) Show Π Π. A method to show an optmzaton problem

More information

Lecture 21: Numerical methods for pricing American type derivatives

Lecture 21: Numerical methods for pricing American type derivatives Lecture 21: Numercal methods for prcng Amercan type dervatves Xaoguang Wang STAT 598W Aprl 10th, 2014 (STAT 598W) Lecture 21 1 / 26 Outlne 1 Fnte Dfference Method Explct Method Penalty Method (STAT 598W)

More information

n ). This is tight for all admissible values of t, k and n. k t + + n t

n ). This is tight for all admissible values of t, k and n. k t + + n t MAXIMIZING THE NUMBER OF NONNEGATIVE SUBSETS NOGA ALON, HAROUT AYDINIAN, AND HAO HUANG Abstract. Gven a set of n real numbers, f the sum of elements of every subset of sze larger than k s negatve, what

More information

Introduction to Interfacial Segregation. Xiaozhe Zhang 10/02/2015

Introduction to Interfacial Segregation. Xiaozhe Zhang 10/02/2015 Introducton to Interfacal Segregaton Xaozhe Zhang 10/02/2015 Interfacal egregaton Segregaton n materal refer to the enrchment of a materal conttuent at a free urface or an nternal nterface of a materal.

More information

A Result on a Cyclic Polynomials

A Result on a Cyclic Polynomials Gen. Math. Note, Vol. 6, No., Feruary 05, pp. 59-65 ISSN 9-78 Copyrght ICSRS Pulcaton, 05.-cr.org Avalale free onlne at http:.geman.n A Reult on a Cyclc Polynomal S.A. Wahd Department of Mathematc & Stattc

More information

Advanced Quantum Mechanics

Advanced Quantum Mechanics Advanced Quantum Mechancs Rajdeep Sensarma! sensarma@theory.tfr.res.n ecture #9 QM of Relatvstc Partcles Recap of ast Class Scalar Felds and orentz nvarant actons Complex Scalar Feld and Charge conjugaton

More information

Appendix B. Criterion of Riemann-Stieltjes Integrability

Appendix B. Criterion of Riemann-Stieltjes Integrability Appendx B. Crteron of Remann-Steltes Integrablty Ths note s complementary to [R, Ch. 6] and [T, Sec. 3.5]. The man result of ths note s Theorem B.3, whch provdes the necessary and suffcent condtons for

More information

Random Walks on Digraphs

Random Walks on Digraphs Random Walks on Dgraphs J. J. P. Veerman October 23, 27 Introducton Let V = {, n} be a vertex set and S a non-negatve row-stochastc matrx (.e. rows sum to ). V and S defne a dgraph G = G(V, S) and a drected

More information

Economics 101. Lecture 4 - Equilibrium and Efficiency

Economics 101. Lecture 4 - Equilibrium and Efficiency Economcs 0 Lecture 4 - Equlbrum and Effcency Intro As dscussed n the prevous lecture, we wll now move from an envronment where we looed at consumers mang decsons n solaton to analyzng economes full of

More information

Spectral Graph Theory and its Applications September 16, Lecture 5

Spectral Graph Theory and its Applications September 16, Lecture 5 Spectral Graph Theory and ts Applcatons September 16, 2004 Lecturer: Danel A. Spelman Lecture 5 5.1 Introducton In ths lecture, we wll prove the followng theorem: Theorem 5.1.1. Let G be a planar graph

More information

Lecture 20: Lift and Project, SDP Duality. Today we will study the Lift and Project method. Then we will prove the SDP duality theorem.

Lecture 20: Lift and Project, SDP Duality. Today we will study the Lift and Project method. Then we will prove the SDP duality theorem. prnceton u. sp 02 cos 598B: algorthms and complexty Lecture 20: Lft and Project, SDP Dualty Lecturer: Sanjeev Arora Scrbe:Yury Makarychev Today we wll study the Lft and Project method. Then we wll prove

More information

Vapnik-Chervonenkis theory

Vapnik-Chervonenkis theory Vapnk-Chervonenks theory Rs Kondor June 13, 2008 For the purposes of ths lecture, we restrct ourselves to the bnary supervsed batch learnng settng. We assume that we have an nput space X, and an unknown

More information

Lecture Space-Bounded Derandomization

Lecture Space-Bounded Derandomization Notes on Complexty Theory Last updated: October, 2008 Jonathan Katz Lecture Space-Bounded Derandomzaton 1 Space-Bounded Derandomzaton We now dscuss derandomzaton of space-bounded algorthms. Here non-trval

More information

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix Lectures - Week 4 Matrx norms, Condtonng, Vector Spaces, Lnear Independence, Spannng sets and Bass, Null space and Range of a Matrx Matrx Norms Now we turn to assocatng a number to each matrx. We could

More information

C/CS/Phy191 Problem Set 3 Solutions Out: Oct 1, 2008., where ( 00. ), so the overall state of the system is ) ( ( ( ( 00 ± 11 ), Φ ± = 1

C/CS/Phy191 Problem Set 3 Solutions Out: Oct 1, 2008., where ( 00. ), so the overall state of the system is ) ( ( ( ( 00 ± 11 ), Φ ± = 1 C/CS/Phy9 Problem Set 3 Solutons Out: Oct, 8 Suppose you have two qubts n some arbtrary entangled state ψ You apply the teleportaton protocol to each of the qubts separately What s the resultng state obtaned

More information

Foundations of Arithmetic

Foundations of Arithmetic Foundatons of Arthmetc Notaton We shall denote the sum and product of numbers n the usual notaton as a 2 + a 2 + a 3 + + a = a, a 1 a 2 a 3 a = a The notaton a b means a dvdes b,.e. ac = b where c s an

More information

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons

More information

Limited Dependent Variables

Limited Dependent Variables Lmted Dependent Varables. What f the left-hand sde varable s not a contnuous thng spread from mnus nfnty to plus nfnty? That s, gven a model = f (, β, ε, where a. s bounded below at zero, such as wages

More information

Design of Recursive Digital Filters IIR

Design of Recursive Digital Filters IIR Degn of Recurve Dgtal Flter IIR The outut from a recurve dgtal flter deend on one or more revou outut value, a well a on nut t nvolve feedbac. A recurve flter ha an nfnte mule reone (IIR). The mulve reone

More information

Communication Complexity 16:198: February Lecture 4. x ij y ij

Communication Complexity 16:198: February Lecture 4. x ij y ij Communcaton Complexty 16:198:671 09 February 2010 Lecture 4 Lecturer: Troy Lee Scrbe: Rajat Mttal 1 Homework problem : Trbes We wll solve the thrd queston n the homework. The goal s to show that the nondetermnstc

More information

Linear Approximation with Regularization and Moving Least Squares

Linear Approximation with Regularization and Moving Least Squares Lnear Approxmaton wth Regularzaton and Movng Least Squares Igor Grešovn May 007 Revson 4.6 (Revson : March 004). 5 4 3 0.5 3 3.5 4 Contents: Lnear Fttng...4. Weghted Least Squares n Functon Approxmaton...

More information

A 2D Bounded Linear Program (H,c) 2D Linear Programming

A 2D Bounded Linear Program (H,c) 2D Linear Programming A 2D Bounded Lnear Program (H,c) h 3 v h 8 h 5 c h 4 h h 6 h 7 h 2 2D Lnear Programmng C s a polygonal regon, the ntersecton of n halfplanes. (H, c) s nfeasble, as C s empty. Feasble regon C s unbounded

More information

9 Characteristic classes

9 Characteristic classes THEODORE VORONOV DIFFERENTIAL GEOMETRY. Sprng 2009 [under constructon] 9 Characterstc classes 9.1 The frst Chern class of a lne bundle Consder a complex vector bundle E B of rank p. We shall construct

More information

Week 5: Neural Networks

Week 5: Neural Networks Week 5: Neural Networks Instructor: Sergey Levne Neural Networks Summary In the prevous lecture, we saw how we can construct neural networks by extendng logstc regresson. Neural networks consst of multple

More information

Solution Methods for Time-indexed MIP Models for Chemical Production Scheduling

Solution Methods for Time-indexed MIP Models for Chemical Production Scheduling Ian Davd Lockhart Bogle and Mchael Farweather (Edtor), Proceedng of the 22nd European Sympoum on Computer Aded Proce Engneerng, 17-2 June 212, London. 212 Elever B.V. All rght reerved. Soluton Method for

More information

P exp(tx) = 1 + t 2k M 2k. k N

P exp(tx) = 1 + t 2k M 2k. k N 1. Subgaussan tals Defnton. Say that a random varable X has a subgaussan dstrbuton wth scale factor σ< f P exp(tx) exp(σ 2 t 2 /2) for all real t. For example, f X s dstrbuted N(,σ 2 ) then t s subgaussan.

More information

REAL ANALYSIS I HOMEWORK 1

REAL ANALYSIS I HOMEWORK 1 REAL ANALYSIS I HOMEWORK CİHAN BAHRAN The questons are from Tao s text. Exercse 0.0.. If (x α ) α A s a collecton of numbers x α [0, + ] such that x α

More information

APPROXIMATE FUZZY REASONING BASED ON INTERPOLATION IN THE VAGUE ENVIRONMENT OF THE FUZZY RULEBASE AS A PRACTICAL ALTERNATIVE OF THE CLASSICAL CRI

APPROXIMATE FUZZY REASONING BASED ON INTERPOLATION IN THE VAGUE ENVIRONMENT OF THE FUZZY RULEBASE AS A PRACTICAL ALTERNATIVE OF THE CLASSICAL CRI Kovác, Sz., Kóczy, L.T.: Approxmate Fuzzy Reaonng Baed on Interpolaton n the Vague Envronment of the Fuzzy Rulebae a a Practcal Alternatve of the Clacal CRI, Proceedng of the 7 th Internatonal Fuzzy Sytem

More information

Feature Selection: Part 1

Feature Selection: Part 1 CSE 546: Machne Learnng Lecture 5 Feature Selecton: Part 1 Instructor: Sham Kakade 1 Regresson n the hgh dmensonal settng How do we learn when the number of features d s greater than the sample sze n?

More information

Module 9. Lecture 6. Duality in Assignment Problems

Module 9. Lecture 6. Duality in Assignment Problems Module 9 1 Lecture 6 Dualty n Assgnment Problems In ths lecture we attempt to answer few other mportant questons posed n earler lecture for (AP) and see how some of them can be explaned through the concept

More information

BOUNDEDNESS OF THE RIESZ TRANSFORM WITH MATRIX A 2 WEIGHTS

BOUNDEDNESS OF THE RIESZ TRANSFORM WITH MATRIX A 2 WEIGHTS BOUNDEDNESS OF THE IESZ TANSFOM WITH MATIX A WEIGHTS Introducton Let L = L ( n, be the functon space wth norm (ˆ f L = f(x C dx d < For a d d matrx valued functon W : wth W (x postve sem-defnte for all

More information

First Year Examination Department of Statistics, University of Florida

First Year Examination Department of Statistics, University of Florida Frst Year Examnaton Department of Statstcs, Unversty of Florda May 7, 010, 8:00 am - 1:00 noon Instructons: 1. You have four hours to answer questons n ths examnaton.. You must show your work to receve

More information

A new construction of 3-separable matrices via an improved decoding of Macula s construction

A new construction of 3-separable matrices via an improved decoding of Macula s construction Dscrete Optmzaton 5 008 700 704 Contents lsts avalable at ScenceDrect Dscrete Optmzaton journal homepage: wwwelsevercom/locate/dsopt A new constructon of 3-separable matrces va an mproved decodng of Macula

More information

Kernel Methods and SVMs Extension

Kernel Methods and SVMs Extension Kernel Methods and SVMs Extenson The purpose of ths document s to revew materal covered n Machne Learnng 1 Supervsed Learnng regardng support vector machnes (SVMs). Ths document also provdes a general

More information

arxiv: v1 [math.co] 1 Mar 2014

arxiv: v1 [math.co] 1 Mar 2014 Unon-ntersectng set systems Gyula O.H. Katona and Dánel T. Nagy March 4, 014 arxv:1403.0088v1 [math.co] 1 Mar 014 Abstract Three ntersecton theorems are proved. Frst, we determne the sze of the largest

More information

Lecture 2: Gram-Schmidt Vectors and the LLL Algorithm

Lecture 2: Gram-Schmidt Vectors and the LLL Algorithm NYU, Fall 2016 Lattces Mn Course Lecture 2: Gram-Schmdt Vectors and the LLL Algorthm Lecturer: Noah Stephens-Davdowtz 2.1 The Shortest Vector Problem In our last lecture, we consdered short solutons to

More information

Lecture 17: Lee-Sidford Barrier

Lecture 17: Lee-Sidford Barrier CSE 599: Interplay between Convex Optmzaton and Geometry Wnter 2018 Lecturer: Yn Tat Lee Lecture 17: Lee-Sdford Barrer Dsclamer: Please tell me any mstake you notced. In ths lecture, we talk about the

More information

Start Point and Trajectory Analysis for the Minimal Time System Design Algorithm

Start Point and Trajectory Analysis for the Minimal Time System Design Algorithm Start Pont and Trajectory Analy for the Mnmal Tme Sytem Degn Algorthm ALEXANDER ZEMLIAK, PEDRO MIRANDA Department of Phyc and Mathematc Puebla Autonomou Unverty Av San Claudo /n, Puebla, 757 MEXICO Abtract:

More information

Information Acquisition in Global Games of Regime Change (Online Appendix)

Information Acquisition in Global Games of Regime Change (Online Appendix) Informaton Acquton n Global Game of Regme Change (Onlne Appendx) Mchal Szkup and Iabel Trevno Augut 4, 05 Introducton Th appendx contan the proof of all the ntermedate reult that have been omtted from

More information

Notes on Frequency Estimation in Data Streams

Notes on Frequency Estimation in Data Streams Notes on Frequency Estmaton n Data Streams In (one of) the data streamng model(s), the data s a sequence of arrvals a 1, a 2,..., a m of the form a j = (, v) where s the dentty of the tem and belongs to

More information

Statistical Properties of the OLS Coefficient Estimators. 1. Introduction

Statistical Properties of the OLS Coefficient Estimators. 1. Introduction ECOOMICS 35* -- OTE 4 ECO 35* -- OTE 4 Stattcal Properte of the OLS Coeffcent Etmator Introducton We derved n ote the OLS (Ordnary Leat Square etmator ˆβ j (j, of the regreon coeffcent βj (j, n the mple

More information

Quick Visit to Bernoulli Land

Quick Visit to Bernoulli Land Although we have een the Bernoull equaton and een t derved before, th next note how t dervaton for an uncopreble & nvcd flow. The dervaton follow that of Kuethe &Chow ot cloely (I lke t better than Anderon).

More information

Lecture 10: May 6, 2013

Lecture 10: May 6, 2013 TTIC/CMSC 31150 Mathematcal Toolkt Sprng 013 Madhur Tulsan Lecture 10: May 6, 013 Scrbe: Wenje Luo In today s lecture, we manly talked about random walk on graphs and ntroduce the concept of graph expander,

More information

COS 521: Advanced Algorithms Game Theory and Linear Programming

COS 521: Advanced Algorithms Game Theory and Linear Programming COS 521: Advanced Algorthms Game Theory and Lnear Programmng Moses Charkar February 27, 2013 In these notes, we ntroduce some basc concepts n game theory and lnear programmng (LP). We show a connecton

More information

HMMT February 2016 February 20, 2016

HMMT February 2016 February 20, 2016 HMMT February 016 February 0, 016 Combnatorcs 1. For postve ntegers n, let S n be the set of ntegers x such that n dstnct lnes, no three concurrent, can dvde a plane nto x regons (for example, S = {3,

More information

Appendix for Causal Interaction in Factorial Experiments: Application to Conjoint Analysis

Appendix for Causal Interaction in Factorial Experiments: Application to Conjoint Analysis A Appendx for Causal Interacton n Factoral Experments: Applcaton to Conjont Analyss Mathematcal Appendx: Proofs of Theorems A. Lemmas Below, we descrbe all the lemmas, whch are used to prove the man theorems

More information

U.C. Berkeley CS294: Beyond Worst-Case Analysis Handout 6 Luca Trevisan September 12, 2017

U.C. Berkeley CS294: Beyond Worst-Case Analysis Handout 6 Luca Trevisan September 12, 2017 U.C. Berkeley CS94: Beyond Worst-Case Analyss Handout 6 Luca Trevsan September, 07 Scrbed by Theo McKenze Lecture 6 In whch we study the spectrum of random graphs. Overvew When attemptng to fnd n polynomal

More information