On Partial Opimality by Auxiliary Submodular Problems

Size: px
Start display at page:

Download "On Partial Opimality by Auxiliary Submodular Problems"

Transcription

1 On Partal Opmalty by Auxlary Submodular Problem Alexander Shekhovtov Mc., engneer Václav Hlaváč Prof., head of the Center for Machne Percepton Center for Machne Percepton, Department of Cybernetc Faculty of Electrcal Engneerng, Czech Techncal Unverty n Prague Techncka 2, Prague 6, Czech Republc Abtract In th work, we prove everal relaton between three dfferent energy mnmzaton technque. A recently propoed method for determnng a provably optmal partal agnment of varable by Ivan Kovtun (IK), the lnear programmng relaxaton approach (LP) and the popular expanon move algorthm by Yur Boykov. We propoe a novel uffcent condton of optmal partal agnment, whch baed on LP relaxaton and called LP-autarky. We how that method of Kovtun, whch buld auxlary ubmodular problem, fulfll th uffcent condton. The followng lnk thu etablhed: LP relaxaton cannot be tghtened by IK. For non-ubmodular problem th a non-trval reult. In the cae of two label, LP relaxaton provde optmal partal agnment, known a pertency, whch, a we how, domnate IK. Relatng IK wth expanon move, we how that the et of fxed pont of expanon move wth any truncaton rule for the ntal problem and the problem retrcted by one-v-all method of IK would concde.e. expanon move cannot be mproved by th method. In the cae of two label, expanon move wth a partcular truncaton rule concde wth one-v-all method. Keyword: energy mnmzaton, partal optmalty, pertency, max-um, WCSP, MRF, autarky, LP-relaxaton, expanon move. 1 Introducton 1.1 Energy Mnmzaton In th work 1 we conder mnmzaton problem of the followng form: [ ] mn x L f 0 + V f (x ) + t E f t (x t ) = mn f(x). (1) x L 1 The work wa upported bu EU project FP7-ICT NIFT and FP7-ICT HUMAVIPS and the Czech project 1M0567 CAK. Here, V a fnte et and E V V. A concatenated vector of all varable x = (x V) called a labelng. Varable x take t value n a dcrete doman L, called label. Labelng x take value n L, the Cartean product of all doman L. In th paper all L wll have the ame number of label, but may have dfferent aocated orderng, etc. Notaton t denote the ordered par (, t) and x t denote the par of correpondng varable, (x, x t ). The objectve compoed of term f 0 R and functon f : L R and f t : L L t R. The problem (1) condered n everal feld. It 1

2 alo known a the labelng problem, the Weghted Contrant Satfacton (WCSP) and for the cae of two label ( L = 2, ) a the peudo-boolean 2 optmzaton [1]. Our termnology come from conderng probabltc model n the form of Gbb dtrbuton. There certan dfference between problem wth two label and more than two label, the later wll be referred to a mult-label problem. 1.2 Partal Optmalty Energy mnmzaton (1) an NP-hard problem n general. Technque whch allow u to fnd a part of the optmal labelng are of our central nteret here. The dea that t may be poble to fx a part of varable to take certan label uch that any optmal labelng wll provably have the ame partal agnment. More precely, we conder a ubet of varable A V and the agnment of label over th ubet y = (y A). The par (A, y) called a trong optmal partal agnment (trong pertency [3]), f for any mnmzer x t hold x A = y, where notaton x A the retrcton of x to A,.e. (x A). Lkewe, f there ext at leat one mnmzer x, for whch x A = y hold we ay that (A, y) a weak optmal partal agnment. Two or more trong optmal partal agnment can be combned together, becaue each of them preerve all optmal oluton. Th not true for weak agnment, even f they agn dfferent varable, they may not hare any globally optmal oluton n common. However, f we want to fnd a mnmzer of (1) (or at leat localze t a much a poble), a weak optmal partal agnment could be more helpful the bet one agn all varable. 1.3 Doman Contrant The dea of optmal partal agnment naturally extend to contranng a varable to a ubet of label K L. Let A V, let K L, A. Let K be the Cartean product of K, A. We ay that a par (A, K) a trong (rep. weak) optmal contrant f 2 Varable x {0, 1} are regarded a Boolean n th cae and peudo emphaze that a real-valued rather than Boolean functon of thee varable condered. x A K for all (rep. at leat one) mnmzer x. Th type of contrant called doman contrant. Obvouly, t nclude partal agnment a a pecal cae. 1.4 Autarke Some doman contrant follow from more pecfc properte called autarke. Th term occur n [3] for two-label problem and we conder t extenon [15] to mult-label problem. Let L = {0, 1... L} V, L N. Let x, y L. Defne component-we mnmum and maxmum of two labellng: (x y) = mn(x, y ), (2a) (x y) = max(x, y ). (2b) A par (x mn L, x max L) uch that x mn x max (component-we) called a weak autarky for problem (1), f x L f((x x mn ) x max ) f(x). (3) If addtonally for any x (x x mn ) x max trct nequalty f((x x mn ) x max ) < f(x) (4) hold, then the autarky called trong. The autarky provde doman contrant wth K = [x mn,..., x max ]. For any mnmzer x, we have that ˆx = (x x mn ) x max a mnmzer a well, and ˆx K. A trong autarky guarantee addtonally that x mut telf atfy x K. Indeed, f t wa not true then ˆx x and f(ˆx) < f(x ), whch a contradcton. Therefore a weak (rep. trong) autarky provde a weak (rep. trong) doman contrant. Determnng whether a gven par (x mn, x max ) a trong autarky an NP-hard decon problem [3]. Autarke can be combned together. A jon of two autarke (x 1, x 2 ), (y 1, y 2 ) the par (x 1 y 1, x 2 y 2 ). For trong autarke, the reult a trong autarky and th operaton commutatve, aocatve and dempotent, o that t defne a em-lattce. Proof. From defnton of autarke, we have f((((x x 1 ) x 2 ) y 1 ) y 2 ) f((x x 1 ) x 2 ) f(x) (5) 2

3 Note, that for x 1 x 2 we have (x x 1 ) x 2 = (x x 2 ) x 1. We can rewrte the labelng n the left hand de (LHS) a follow ((x x 1 ) x 2 ) y 1 ) y 2 = ((x x 2 ) (x 1 y 1 )) y 2. = (x (x 2 y 2 )) (x 1 y 1 ), (6) where doted equalty hold f y 2 x 1. Th atfed for trong autarke, becaue t would be a contradcton that all optmal labellng are below y 2 and above x 1. Thu there ext an autarky, whch provde the maxmal amount of doman contrant among trong autarke. It the jon of all trong autarke. It alo poble to jon non-contradctve weak autarke together, but let u leave t ade for now. We wll conder a pecal cae of autarke wth one-de contrant, of the form (x mn, L) or (0, x max ), where L and 0 repreent the labelng wth all component equal to L (rep. 0). For uch autarke nequalty (3) mplfe, becaue x 0 = x and x L = x. Method [10, 11] compute trong autarke of th form. By takng the jon of two trong autarke (x mn, L) and (0, x max ) we can obtan a trong autarky (x mn, x max ). However, the revere not true: f (x mn, x max ) a trong autarky, t doe not mply that (x mn, L) or (0, x max ) an autarky. And t the cae that other method (roof-dual [1] n the cae of two-label problem and t mult-label extenon [15]) can fnd an autarky of the form (x mn, x max ), whch not a jon of two one-de autarke. be t mnmzer. Then we have the followng properte: f(x x ) f(x), f(x x ) f(x). (9a) (9b) They ealy follow from ubmodularty, notng that f(x x ) f(x ) and f(x x ) f(x ). So, n fact, any par of optmal oluton (x 1, x 2 ) a weak autarky for th problem. Moreover, f we let x mn = arg mn x f(x), (10a) x max = arg mn x f(x), (10b) where arg mn the et of mnmzer, we ee that both x mn and x max are mnmzer of f and that (x mn, x max ) a trong autarky for f. In fact, t the jon of all trong autarke for f. Th trong autarky can be determned from a oluton of the correpondng maxflow problem. 2 Approach by Kovtun In th ecton, we revew technque [10, 11] for buldng autarke (and hence doman contran) by contructng auxlary problem. We take a omewhat dfferent perpectve on thee reult, however, our tatement and proof here are n a ene equvalent to one gven n [10, 11]. Theorem 1. Let f = g + h, let (x mn, x max ) be a trong autarky for g and a weak autarky for h. Then (x mn, x max ) a trong autarky for f. 1.5 Submodular Problem Functon f called ubmodular f x, y L f(x y) + f(x y) f(x) + f(y). (7) Proof. We have f((x x mn ) x max ) = g((x x mn ) x max ) + h((x x mn ) x max ) g(x) + h(x), (11) In the cae f defned by (1), t ubmodular ff (ee e.g. [18]) t E, x t, y t L t = L L t f t (x t ) + f t (y t ) f t (x t y t ) + f t (x t y t ). (8) Mnmzng a parwe ubmodular functon reduce to mncut problem [7], [13]. Let f be ubmodular and x and the nequalty trct f (x mn, x max ) trong for ether h or g. The dea of auxlary problem to contruct a ubmodular g, for whch, a we know, a trong autarky (x mn, L) can be found by choong x mn a the lowet mnmzer of g, gven by (10a). The trck to fnd 3

4 uch g that (x mn, L) at the ame tme an autarky for h = f g. The followng uffcent condton were propoed [11]: Statement 1. Let h atfy V, x L, ˆx K h (x ˆx ) h (ˆx ) t E, x t L t, ˆx t K t h t (x t ˆx t ) h t (x t ). (12a) (12b) Then for any x mn uch that x mn K, the par (x mn, L) a weak autarky for h. If addtonally V, ˆx K, x < ˆx h (x ˆx ) < h (ˆx ), then (x mn, L) a trong autarky. (12c) Proof. For any x L, ummng correpondng nequalte from (12a) and (12b), we obtan h (x x mn ) + t h t (x t x mn t ) h (x ) + t h t (x t ). If x x mn x, then V x < x mn mple trct nequalty. (13) and (12c) Two practcal method were propoed [11] to contruct g and (K V). We frt decrbe a more general approach. Algorthm 1: Sequental contructon of g, (K V), [11] 1. Start wth K =, V; 2. Fnd g uch that h = f g atfe (12) and g atfe ubmodularty contrant (8). 3. Fnd x mn = arg mn x g(x); 4. If x mn K for all V then top. 5. Set K K {x mn } V and go to tep 2. In tep 2 for each edge t E a ytem of lnear nequalte n g t ha to be olved. Whle [11] provde an explct oluton, for our conderaton t wll not be neceary. When the algorthm top, g ubmodular and (x mn, L) a trong autarky for g and a weak autarky for f g. By Theorem 1, t a trong autarky for f. It may top, however, wth x mn = 0 for all, o that effcently no contrant are derved. Beng a polynomal algorthm t cannot have a guarantee to mplfy the problem (1). A mpler non-teratve method propoed n [10] hown n Algorthm 2. It attempt to dentfy node where the label L better than any other label. The contructed auxlary problem g ha a property that t lowet mnmzer x mn = arg mn g(x) guar- Algorthm 2: One v all method, [10, 11] 1. For each choe uch orderng of L that 0 arg mn f (). L 2. Set g = f, V. 3. Set K = {0, L}. a t, = L, j = L, b t, = L, j L, 4. Set g t (, j) = c t, L, j = L, d t, L, j L, where a t, b t, c t, d t are uch that f t g t atfy (12b) and ubmodularty contrant. One of the oluton a follow: a t = f t (L, L), b t = mn j L f t(l, j), c t = mn f t(, L), L ( d t = mn b t + c t a t, mn L,j L [ f t (, j) + mn { b t f t (L, j), c t f t (, L) }]). (14) anteed to atfy x mn K V. (becaue all cot (g t (, j) < L, j < L) are equal and g (0) g () V, < L, ee proof n [10]). Therefore (x mn, L) a weak autarky for f g and Theorem (1) apple. Both method allow u to chooe varou orderng of et L. Strong doman contrant derved from varou orderng can be then combned. 4

5 3 LP-autarke In th ecton we ntroduce a pecal ubcla of autarke, whch preerve optmal oluton of the LPrelaxaton. Unlke wth general autarke, the memberhp to th ubcla polynomally verfable. We how that autarke contructed by algorthm 1, 2 belong to th ubcla. Th ha ueful mplcaton for LP relaxaton. 3.1 LP Relaxaton Let φ(x) be a vector wth component φ(x) 0 = 1, φ(x), = [[x =]] and φ(x) t, = [[x t =]], where [[ ]] 1 f the expreon nde true and 0 otherwe. Let f denote a vector wth component f 0, f, = f () and f t, = f t (). Wth repect to component of energy functon we wll be ung th ndex and parenthe notaton completely nterchangeably. Let, denote a calar product. Then we can wrte energy mnmzaton a mn f, φ(x). (15) x L It relaxaton to a lnear program wrtten a mn f, µ, (16) µ Λ where Λ the local polytope. It approxmate conv{φ(x) x L} from the outde, ee e.g. [18] for more detal. It gven by the lnear contrant µ 0 = 1, µ, 0, µ t, 0, L t µ t, = 1 t E, j L t µ t, = µ, L, t E, L µ t, = µ t,j j L t, t E. Vector µ Λ called a relaxed labelng. 3.2 LP-autarky (17) We now extend the noton of autarky to relaxed labellng. Defnton 1. A bnary operaton : Λ L Λ, defned a follow. Let y L and µ Λ. Then ν = µ y Λ contructed a: ν, = ν t, = µ,, < y, y µ,, = y, 0, > y ; µ t,, < y, j < y t, y µ t, j, = y, j < y t, j y t µ t,, < y, j = y t, µ t, j, = y, j = y t, y j y t 0, > y or j > y t. (18a) (18b) By contructon, the relaxed labelng ν ha non-zero weght only for label below y: ν, = 0 for > y and the ame for par t,. Let u check that ν Λ. Proof. Normalzaton contrant: ν, = µ, + µ, = µ, = 1. (19) <y y Margnalzaton contrant: µ t, + µ t, j, j < y t, <y y ν t, = µ t, + µ t, j, j = y t, <y y j y t j y t 0, j > y t, = ν t,j. (20) Operaton ν = µ y defned completely mlarly, havng ngleton component µ,, > y, (µ y), = µ,, = y, (21) y 0, < y. Defnton 2. We ay that a par (x mn, x max ) a weak LP-autarky for f, f µ Λ f, (µ x mn ) x max f, µ. (22) If addtonally for all µ uch that (µ x mn ) x max µ the trct nequalty hold then we ay that t a trong LP-autarky. 5

6 3.3 Properte of LP-autarke Statement 2. Any weak (rep. trong) LP-autarky a weak (rep. trong) autarky. Proof. By ubttutng µ = φ(x). Statement 3. Checkng whether (x mn, x max ) an LP-autarky for f can be olved n a polynomal tme. Proof. By contructon, (µ x mn ) x max a lnear map n µ, let u denote t Aµ. Inequalty (22) hold ff mn f, µ Aµ 0, (23) µ Λ whch a lnear program. To verfy whether A a trong LP-autarky we need to olve mn f, µ Aµ > 0 µ Λ,.t. µ, < V. x mn x max Statement 4. If f ubmodular, then µ Λ, y L µ, f + φ(y), f µ y, f + µ y, f. (24) (25) Proof. Scalar product n (25) are compoed of um of ngleton term and parwe term. We frt how that um of ngleton term are equal, expandng ngleton term n the rght hand de (RHS): [ ] (µ y), + (µ y), f() = µ,f () + µ, f (y )+ <y y µ,f () + µ, f (y ) = >y y µ,f () + ( ) µ, f (y ) = µ,f () + [[= y ]]f (y ). Now conder ubmodularty contrant: t E, L t, y t L t f t () + f t (y t ) f t ( y t ) + f t ( y t ). (26) (27) Multplyng th nequalty by µ t, and ummng over, we obtan on the LHS: µ t, f t () + f t (y t ) = µ t, f t () + [[=y t ]]f t (y t ) and on the RHS: [ µ t, ft ( y t ) + f t ( y t ) ] = [ ] (µ y)t, + (µ y) t, ft (), where the equalty verfed a follow: µ t, f t ( y t ) = µ t,f t() + µ t,f t(y, j)+ <y y j<y t j<y t µ t,f t(, y t ) + µ t,f t(y t) = <y y j y t j y t (µ y) t, f t,. (28) (29) (30) The term wth rewrtten mlarly. By ummng nequalte (28) (29) over t E and addng equalte (26) of the ngleton term, we get the reult. Statement 5. Let f be ubmodular and x arg mn f(x). Then µ Λ x µ x, f µ, f, (31a) µ x, f µ, f. (31b) Proof. Let u how (31a). For ubmodular problem LP-relaxaton (16) tght. Thu for any ν Λ there hold ν, f f(x ) = φ(x ), f. In partcular, for ν = µ y we have µ y, f φ(x ), f, whch when combned wth (25) mple the tatement. Statement 6. Let (x mn, L) be a trong LP-autarky for f, then: V, < x mn, µ arg mn µ, f µ Λ µ, = 0. (32) Proof. Let µ arg mn µ Λ µ, f and µ, > 0. Then µ x mn µ and f(µ x mn ) < f(µ ), whch contradct optmalty of µ. 6

7 3.4 Implcaton for Algorthm 1, 2 We have already een n tatement 5 that for a ubmodular problem g, takng y a a mnmzer (rep. the lowet mnmzer) of g gve a weak (rep. trong) LPautarky (y, L). Let u how that tatement 1 extend to LP-autarke too. Th would mply that autarke derved by algorthm 1, 2 are n fact LP-autarke for f = g + h. Statement 7. Let h atfy nequalte (12). Then for any y L uch that y K, the par (y, L) a weak LP-autarky for h. Proof. Let µ Λ. From nequalty (12a) we have ((µ y), µ, )h, 0. (33) Multplyng (12b) by µ t, and ummng over L t and over t E we obtan [ ] (µ y)t, µ t, ht, 0. (34) t Addng (33) and (34), we get: whch equvalent to (22). µ y µ, h 0, (35) We have hown that algorthm 1, 2 derve doman contrant n the form of trong LP-autarke. We know too that optmal oluton of LP-relaxaton wll obey doman contrant derved va trong LPautarke. Note, whle algorthm 1, 2 depend on the orderng of the label, oluton of the LP-relaxaton doe not. Hence, Corollary 1. Let (K L V) be a trong doman contrant derved by Algorthm 1, 2 w.r.t. any orderng of et L. Then the et of optmal oluton of LP relaxaton wth and wthout thee doman contrant would concde. We proved that LP relaxaton cannot be tghtened by algorthm 1, 2. It may only be mplfed by elmnatng all varable whch are guaranteed to be 0 n every optmal oluton. Th may be ueful n practcal method olvng LP relaxaton. For problem wth two label, the followng relaton alo hold. Let Λ = arg mn µ Λ f, µ. Let x mn = mn{ µ Λ µ, > 0}, x max = max{ µ Λ µ, > 0}, (36) then (x mn, x max ) a trong autarky for f. Th the roof-dual autarky [1]. Becaue for any other autarky derved va algorthm 1 and 2 tatement 6 hold, we conclude that roof-dual autarky domnate algorthm 1 and 2. 4 Expanon Move Expanon move algorthm [4] eek to mprove the current oluton x by conderng a move, whch for every V ether keep the current label x or change t to the label k. Algorthm 3: Expanon-Move [4] 1. Let x L, let k L. The move energy functon g(z) of bnary confguraton z {0, 1} V defned by: g 0 = f 0, g (0) = f (x ), g (1) = f (k), g t (1, 1) = f t (k, k), g t (1, 0) = f t (k, x t ), g t (0, 1) = f t (x, k), g t (0, 0) = f t (x, x t ). (37) 2. Let z arg mn g(z). z 3. If g(z ) < g(0), agn x { x, f z = 0, k, f z = 1. If the above procedure repeated for all label k L and no mprovement to x found then x ad to be a fxed pont of th method. In the cae f a metrc energy [4], the move energy g ubmodular for arbtrary x and tep 2 eay. Statement 8. Let f be metrc [4]. Let (x mn, L) be a trong autarky for f uch that x mn {0, L}, V. Then for any fxed pont x of the expanon-move algorthm there hold x x mn. (38) 7

8 Proof. Aume V uch that x < x mn. Then f(x x mn ) < f(x) and nce x {1, L}, t x x mn x, x mn = 1, = (39) L, x mn = L, whch a vald expanon move from x to label k = L, trctly mprovng the energy. In the cae when a move energy not ubmodular, t can be truncated to make t ubmodular whle tll preervng the property that the move doe not ncreae f(x) [12]. Let t = g t (1, 1) + g t (0, 0) g t (0, 1) g t (1, 0). Par t ubmodular ff t < 0. Defnton 3. Truncaton g of g dfferent from g only n non-ubmodular parwe component of g, whch are et a: g t,00 = g t,00 β t t, g t,01 = g t,01 + α t t, g t,10 = g t,10 + (1 α t β t ) t, g t,11 = g t,11, (40) where α t and β t are free parameter, atfyng α t 0, β t 0, α t + β t 1. It eay to verfy that g ubmodular, and g(z) g(0) g (z) g (0), (41) ayng that ncreae n g no more than ncreae n g when changng from 0 to z. Proof of (41). By contructon of g, for all t E uch that t > 0, enumeratng all z t, g t,00 g t,00 = 0, g t,01 g t,00 = g t,01 g t,00 + (α t + β t ) t, g t,10 g t,00 = g t,10 g t,00 + (1 α t ) t, g t,11 g t,00 = g t,11 g t,00 + β t t, (42) we ee that only potve value are added on RHS. It alo een that the added potve value do only ncreae wth β t. Th mean that the truncaton wth β t > 0 (let denote t g α,β ) never better than the truncaton wth β = 0 (let denote t g α ): z g(z) g(0) g α (z) g α (0) g α,β (z) g α,β (0). (43) Smlarly, the truncaton wth α = 0, β = 1 (g 0,1 ) not better than the truncaton g α,β : g α,β (z) g α,β (0) g 0,1 (z) g 0,1 (0). (44) Th verfed by examnng component: g 0,1 t (z t ) g 0,1 t (0) g α,β t (z t ) + g α,β t (0) = 0, z t = 00, t (1 (α + β)), z t = 01, t (1 (1 α)), z t = 10, t (1 β), z t = 11, 0. (45) If z an mprovng move for g 0,1 then t alo an mprovng move for any truncaton. We have the followng reult about Algorthm 2: Statement 9. Let (x mn, L) be a trong autarky for f obtaned by Algorthm 2. Let x be a fxed pont of the expanon-move algorthm wth any truncaton rule. Then x x mn. (46) Proof. We wll prove that the tatement hold for truncaton (α = 0, β = 1). We need to how that for a move from x to x x mn the truncated energy decreae at leat a much a doe auxlary problem bult by Alg. 2. Th can be verfed by npectng parwe component for the 4 cae z t = 00, 01, 10, Concluon We propoe a novel repreentaton of method [10, 11] a dervng doman contrant va LP-autarke. Th allow for comparon wth other method dervng doman contrant n the ame form [3, 15] and etablhng relaton wth common method of (approxmate) optmzaton. We alo beleve that label domnaton condton propoed by [5] can be nterpreted n the ame framework, allowng for the theoretcal comparon and or for the degn of combned method. Our reult open everal drecton for mprovement. A drect mprovement to Alg. 2 can be obtaned a follow. Alg. 2 contruct a mult-label auxlary problem, whch equvalent to a two-label problem (nce we know that there a mnmzer wth 8

9 x {0, L}, V). For two label problem, we alo know that the autarky contructed by roof-dual domnate the autarky by truncaton, o t wll be better to et [ d t = mn ft (, j) L,j L + mn { b t f t (L, j), c t f t (, L) }] (47) and olve for roof-dual ung reducton to maxflow [2]. Th would be a non-ubmodular auxlary problem. We can alo attempt to contruct auxlary problem wth mxed ubmodular and upermodular term a n [15] or degn an algorthm whch wll propoe an autarky n ome greedy way and then verfy t va olvng lnear program (24). Reference [1] E. Boro and P.L. Hammer. Peudo-boolean optmzaton. Dcrete Appled Mathematc, (123(1-3)): , [2] E. Boro, P. L. Hammer, and X. Sun. Network flow and mnmzaton of quadratc peudo-boolean functon. Techncal Report RRR , RUTCOR, May [3] E. Boro, P. L. Hammer, and G. Tavare. Preproceng of uncontraned quadratc bnary optmzaton. Techncal Report RRR , RUTCOR, Apr [4] Y. Boykov, O. Vekler, and R. Zabh. Fat approxmate energy mnmzaton va graph cut. IEEE Tranacton on Pattern Analy and Machne Intellgence, 23(11): , Nov [5] J. Demet, M. D. Maeyer, B. Haze, and I. Later. The dead-end elmnaton theorem and t ue n proten de-chan potonng. Nature, 356: , [6] P.L. Hammer, P. Hanen, and B. Smeone. Roof dualty, complementaton and pertency n quadratc 0-1 optmzaton. Math. Programmng, page , [7] Hroh Ihkawa. Exact optmzaton for Markov random feld wth convex pror. IEEE Tranacton on Pattern Analy and Machne Intellgence, 25(10): , [8] Vladmr Kolmogorov. Convergent tree-reweghted meage pang for energy mnmzaton. In Robert G. Cowell and Zoubn Ghahraman, edtor, AI and Stattc, page Socety for Artfcal Intellgence and Stattc, [9] Vladmr Kolmogorov and Martn Wanwrght. On the optmalty of tree-reweghted max-product meage pang. In To appear n 21t Conference on Uncertanty n Artfcal Intellgence (UAI), July [10] I. Kovtun. Partal optmal labelng earch for a NPhard ubcla of (max, +) problem. In DAGM- Sympoum, page , [11] I. Kovtun. Image egmentaton baed on uffcent condton of optmalty n NP-complete clae of tructural labellng problem. PhD the, IRTC ITS Natonal Academy of Scence, Ukrane, In Ukranan. [12] Carten Rother, Sanjv Kumar, Vladmr Kolmogorov, and Andrew Blake. Dgtal tapetry. In CVPR 05: Proceedng of the 2005 IEEE Computer Socety Conference on Computer Von and Pattern Recognton (CVPR 05) - Volume 1, page , Wahngton, DC, USA, IEEE Computer Socety. [13] Dmtr Schlenger and Bor Flach. Tranformng an arbtrary mnum problem nto a bnary one. Reearch Report TUD-FI06-01, Dreden Unverty of Technology, Aprl [14] M.I. Schlenger. Syntactc analy of two-dmenonal vual gnal n the preence of noe. Cybernetc and Sytem Analy, 12: , [15] A. Shekhovtov, V. Kolmogorov, P. Kohl, V. Hlavac, C. Rother, and P. Torr. LP-relaxaton of bnarzed energy mnmzaton.. Reearch Report CTU CMP , Czech Techncal Unverty, [16] Martn Wanwrght, Tomm Jaakkola, and Alan Wllky. Exact MAP etmate by (hyper)tree agreement. In S. Thrun S. Becker and K. Obermayer, edtor, Advance n Neural Informaton Proceng Sytem 15, page MIT Pre, [17] Tomáš Werner. A lnear programmng approach to max-um problem: A revew. Reearch Report CTU CMP , Center for Machne Percepton, Czech Techncal Unverty, Dec [18] Tomáš Werner. A lnear programmng approach to max-um problem: A revew. IEEE Tranacton on Pattern Analy and Machne Intellgence, 29(7): , July

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

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

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

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

On the U-WPF Acts over Monoids

On the U-WPF Acts over Monoids Journal of cence, Ilamc Republc of Iran 8(4): 33-38 (007) Unverty of Tehran, IN 06-04 http://jcence.ut.ac.r On the U-WPF ct over Monod. Golchn * and H. Mohammadzadeh Department of Mathematc, Unverty of

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

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

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

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

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

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

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

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

Solutions to exam in SF1811 Optimization, Jan 14, 2015

Solutions to exam in SF1811 Optimization, Jan 14, 2015 Solutons to exam n SF8 Optmzaton, Jan 4, 25 3 3 O------O -4 \ / \ / The network: \/ where all lnks go from left to rght. /\ / \ / \ 6 O------O -5 2 4.(a) Let x = ( x 3, x 4, x 23, x 24 ) T, where the varable

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

princeton univ. F 17 cos 521: Advanced Algorithm Design Lecture 7: LP Duality Lecturer: Matt Weinberg

princeton univ. F 17 cos 521: Advanced Algorithm Design Lecture 7: LP Duality Lecturer: Matt Weinberg prnceton unv. F 17 cos 521: Advanced Algorthm Desgn Lecture 7: LP Dualty Lecturer: Matt Wenberg Scrbe: LP Dualty s an extremely useful tool for analyzng structural propertes of lnear programs. Whle there

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

A METHOD TO REPRESENT THE SEMANTIC DESCRIPTION OF A WEB SERVICE BASED ON COMPLEXITY FUNCTIONS

A METHOD TO REPRESENT THE SEMANTIC DESCRIPTION OF A WEB SERVICE BASED ON COMPLEXITY FUNCTIONS UPB Sc Bull, Sere A, Vol 77, I, 5 ISSN 3-77 A METHOD TO REPRESENT THE SEMANTIC DESCRIPTION OF A WEB SERVICE BASED ON COMPLEXITY FUNCTIONS Andre-Hora MOGOS, Adna Magda FLOREA Semantc web ervce repreent

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

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

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

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

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

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

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

Approximating the Sum Operation for Marginal-MAP Inference

Approximating the Sum Operation for Marginal-MAP Inference Approxmatng the Sum Operaton or Margnal-MAP Inerence Qang Cheng, Feng Chen, Janwu Dong, Wenl Xu Alexander Ihler Tnghua Natonal Laboratory or Inormaton Scence and Technology Inormaton and Computer Scence

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

Variable Structure Control ~ Basics

Variable Structure Control ~ Basics Varable Structure Control ~ Bac Harry G. Kwatny Department of Mechancal Engneerng & Mechanc Drexel Unverty Outlne A prelmnary example VS ytem, ldng mode, reachng Bac of dcontnuou ytem Example: underea

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

m = 4 n = 9 W 1 N 1 x 1 R D 4 s x i

m = 4 n = 9 W 1 N 1 x 1 R D 4 s x i GREEDY WIRE-SIZING IS LINEAR TIME Chr C. N. Chu D. F. Wong cnchu@c.utexa.edu wong@c.utexa.edu Department of Computer Scence, Unverty of Texa at Autn, Autn, T 787. ABSTRACT In nterconnect optmzaton by wre-zng,

More information

a new crytoytem baed on the dea of Shmuley and roved t rovably ecure baed on ntractablty of factorng [Mc88] After that n 999 El Bham, Dan Boneh and Om

a new crytoytem baed on the dea of Shmuley and roved t rovably ecure baed on ntractablty of factorng [Mc88] After that n 999 El Bham, Dan Boneh and Om Weak Comote Dffe-Hellman not Weaker than Factorng Koohar Azman, azman@ceharfedu Javad Mohajer mohajer@harfedu Mahmoud Salmazadeh alma@harfedu Electronc Reearch Centre, Sharf Unverty of Technology Deartment

More information

Yong Joon Ryang. 1. Introduction Consider the multicommodity transportation problem with convex quadratic cost function. 1 2 (x x0 ) T Q(x x 0 )

Yong Joon Ryang. 1. Introduction Consider the multicommodity transportation problem with convex quadratic cost function. 1 2 (x x0 ) T Q(x x 0 ) Kangweon-Kyungk Math. Jour. 4 1996), No. 1, pp. 7 16 AN ITERATIVE ROW-ACTION METHOD FOR MULTICOMMODITY TRANSPORTATION PROBLEMS Yong Joon Ryang Abstract. The optmzaton problems wth quadratc constrants often

More information

Discrete Simultaneous Perturbation Stochastic Approximation on Loss Function with Noisy Measurements

Discrete Simultaneous Perturbation Stochastic Approximation on Loss Function with Noisy Measurements 0 Amercan Control Conference on O'Farrell Street San Francco CA USA June 9 - July 0 0 Dcrete Smultaneou Perturbaton Stochatc Approxmaton on Lo Functon wth Noy Meaurement Q Wang and Jame C Spall Abtract

More information

Solutions HW #2. minimize. Ax = b. Give the dual problem, and make the implicit equality constraints explicit. Solution.

Solutions HW #2. minimize. Ax = b. Give the dual problem, and make the implicit equality constraints explicit. Solution. Solutons HW #2 Dual of general LP. Fnd the dual functon of the LP mnmze subject to c T x Gx h Ax = b. Gve the dual problem, and make the mplct equalty constrants explct. Soluton. 1. The Lagrangan s L(x,

More information

On the Multicriteria Integer Network Flow Problem

On the Multicriteria Integer Network Flow Problem BULGARIAN ACADEMY OF SCIENCES CYBERNETICS AND INFORMATION TECHNOLOGIES Volume 5, No 2 Sofa 2005 On the Multcrtera Integer Network Flow Problem Vassl Vasslev, Marana Nkolova, Maryana Vassleva Insttute of

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

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

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0 MODULE 2 Topcs: Lnear ndependence, bass and dmenson We have seen that f n a set of vectors one vector s a lnear combnaton of the remanng vectors n the set then the span of the set s unchanged f that vector

More information

and decompose in cycles of length two

and decompose in cycles of length two Permutaton of Proceedng of the Natona Conference On Undergraduate Reearch (NCUR) 006 Domncan Unverty of Caforna San Rafae, Caforna Apr - 4, 007 that are gven by bnoma and decompoe n cyce of ength two Yeena

More information

Weak McCoy Ore Extensions

Weak McCoy Ore Extensions Internatonal Mathematcal Forum, Vol. 6, 2, no. 2, 75-86 Weak McCoy Ore Extenon R. Mohammad, A. Mouav and M. Zahr Department of Pure Mathematc, Faculty of Mathematcal Scence Tarbat Modare Unverty, P.O.

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

The Minimum Universal Cost Flow in an Infeasible Flow Network

The Minimum Universal Cost Flow in an Infeasible Flow Network Journal of Scences, Islamc Republc of Iran 17(2): 175-180 (2006) Unversty of Tehran, ISSN 1016-1104 http://jscencesutacr The Mnmum Unversal Cost Flow n an Infeasble Flow Network H Saleh Fathabad * M Bagheran

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

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

Lecture 10 Support Vector Machines. Oct

Lecture 10 Support Vector Machines. Oct Lecture 10 Support Vector Machnes Oct - 20-2008 Lnear Separators Whch of the lnear separators s optmal? Concept of Margn Recall that n Perceptron, we learned that the convergence rate of the Perceptron

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

A Novel Approach for Testing Stability of 1-D Recursive Digital Filters Based on Lagrange Multipliers

A Novel Approach for Testing Stability of 1-D Recursive Digital Filters Based on Lagrange Multipliers Amercan Journal of Appled Scence 5 (5: 49-495, 8 ISSN 546-939 8 Scence Publcaton A Novel Approach for Tetng Stablty of -D Recurve Dgtal Flter Baed on Lagrange ultpler KRSanth, NGangatharan and Ponnavakko

More information

ENTROPY BOUNDS USING ARITHMETIC- GEOMETRIC-HARMONIC MEAN INEQUALITY. Guru Nanak Dev University Amritsar, , INDIA

ENTROPY BOUNDS USING ARITHMETIC- GEOMETRIC-HARMONIC MEAN INEQUALITY. Guru Nanak Dev University Amritsar, , INDIA Internatonal Journal of Pure and Appled Mathematc Volume 89 No. 5 2013, 719-730 ISSN: 1311-8080 prnted veron; ISSN: 1314-3395 on-lne veron url: http://.jpam.eu do: http://dx.do.org/10.12732/jpam.v895.8

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

Resonant FCS Predictive Control of Power Converter in Stationary Reference Frame

Resonant FCS Predictive Control of Power Converter in Stationary Reference Frame Preprnt of the 9th World Congre The Internatonal Federaton of Automatc Control Cape Town, South Afrca. Augut -9, Reonant FCS Predctve Control of Power Converter n Statonary Reference Frame Lupng Wang K

More information

The Order Relation and Trace Inequalities for. Hermitian Operators

The Order Relation and Trace Inequalities for. Hermitian Operators Internatonal Mathematcal Forum, Vol 3, 08, no, 507-57 HIKARI Ltd, wwwm-hkarcom https://doorg/0988/mf088055 The Order Relaton and Trace Inequaltes for Hermtan Operators Y Huang School of Informaton Scence

More information

Separation Axioms of Fuzzy Bitopological Spaces

Separation Axioms of Fuzzy Bitopological Spaces IJCSNS Internatonal Journal of Computer Scence and Network Securty VOL3 No October 3 Separaton Axom of Fuzzy Btopologcal Space Hong Wang College of Scence Southwet Unverty of Scence and Technology Manyang

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

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

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

SELECTED SOLUTIONS, SECTION (Weak duality) Prove that the primal and dual values p and d defined by equations (4.3.2) and (4.3.3) satisfy p d.

SELECTED SOLUTIONS, SECTION (Weak duality) Prove that the primal and dual values p and d defined by equations (4.3.2) and (4.3.3) satisfy p d. SELECTED SOLUTIONS, SECTION 4.3 1. Weak dualty Prove that the prmal and dual values p and d defned by equatons 4.3. and 4.3.3 satsfy p d. We consder an optmzaton problem of the form The Lagrangan for ths

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

The Price of Anarchy in a Network Pricing Game

The Price of Anarchy in a Network Pricing Game The Prce of Anarchy n a Network Prcng Game John Muaccho and Shuang Wu Abtract We analyze a game theoretc model of competng network ervce provder that trategcally prce ther ervce n the preence of elatc

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

A Polyhedral Study of Quadratic Traveling Salesman Problems DISSERTATION

A Polyhedral Study of Quadratic Traveling Salesman Problems DISSERTATION A Polyhedral Study of Quadratc Travelng Saleman Problem DISSERTATION ubmtted to Department of Mathematc at Chemntz Unverty of Technology n accordance wth the requrement for the degree Dr. rer. nat. Dpl.-Math.

More information

Strong-Stability-Preserving, K-Step, 5- to 10-Stage, Hermite-Birkhoff Time-Discretizations of Order 12

Strong-Stability-Preserving, K-Step, 5- to 10-Stage, Hermite-Birkhoff Time-Discretizations of Order 12 Amercan Journal of Computatonal Mathematc 20 72-82 do:04236/acm202008 Publhed Onlne June 20 (http://wwwcrporg/ournal/acm) Strong-Stablty-Preervng K-Step 5- to 0-Stage Hermte-Brkhoff Tme-Dcretzaton of Order

More information

Seismic Reliability Analysis and Topology Optimization of Lifeline Networks

Seismic Reliability Analysis and Topology Optimization of Lifeline Networks The 4 th World Conference on Earthquake Engneerng October 2-7, 2008, Beng, Chna Semc Relablty Analy and Topology Optmzaton of Lfelne Network ABSTRACT: Je L and We Lu 2 Profeor, Dept. of Buldng Engneerng,

More information

ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION

ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION Advanced Mathematcal Models & Applcatons Vol.3, No.3, 2018, pp.215-222 ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EUATION

More information

728. Mechanical and electrical elements in reduction of vibrations

728. Mechanical and electrical elements in reduction of vibrations 78. Mechancal and electrcal element n reducton of vbraton Katarzyna BIAŁAS The Slean Unverty of Technology, Faculty of Mechancal Engneerng Inttute of Engneerng Procee Automaton and Integrated Manufacturng

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

1 Convex Optimization

1 Convex Optimization Convex Optmzaton We wll consder convex optmzaton problems. Namely, mnmzaton problems where the objectve s convex (we assume no constrants for now). Such problems often arse n machne learnng. For example,

More information

Method Of Fundamental Solutions For Modeling Electromagnetic Wave Scattering Problems

Method Of Fundamental Solutions For Modeling Electromagnetic Wave Scattering Problems Internatonal Workhop on MehFree Method 003 1 Method Of Fundamental Soluton For Modelng lectromagnetc Wave Scatterng Problem Der-Lang Young (1) and Jhh-We Ruan (1) Abtract: In th paper we attempt to contruct

More information

What Is Optimized in Tight Convex Relaxations for Multi-Label Problems?

What Is Optimized in Tight Convex Relaxations for Multi-Label Problems? What I Optmzed n Tght Convex Relaxaton for Mult-Label Problem? Chrtopher Zach Mcrooft Reearch Cambrdge, UK chzach@mcrooft.com Chrtan Häne ETH Zürch, Swtzerland chaene@nf.ethz.ch Marc Pollefey ETH Zürch,

More information

A CHARACTERIZATION OF ADDITIVE DERIVATIONS ON VON NEUMANN ALGEBRAS

A CHARACTERIZATION OF ADDITIVE DERIVATIONS ON VON NEUMANN ALGEBRAS Journal of Mathematcal Scences: Advances and Applcatons Volume 25, 2014, Pages 1-12 A CHARACTERIZATION OF ADDITIVE DERIVATIONS ON VON NEUMANN ALGEBRAS JIA JI, WEN ZHANG and XIAOFEI QI Department of Mathematcs

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

Computer Control Systems

Computer Control Systems Computer Control ytem In th chapter we preent the element and the bac concept of computercontrolled ytem. The dcretaton and choce of amplng frequency wll be frt examned, followed by a tudy of dcrete-tme

More information

Fixed points of IA-endomorphisms of a free metabelian Lie algebra

Fixed points of IA-endomorphisms of a free metabelian Lie algebra Proc. Indan Acad. Sc. (Math. Sc.) Vol. 121, No. 4, November 2011, pp. 405 416. c Indan Academy of Scences Fxed ponts of IA-endomorphsms of a free metabelan Le algebra NAIME EKICI 1 and DEMET PARLAK SÖNMEZ

More information

Research Article Runge-Kutta Type Methods for Directly Solving Special Fourth-Order Ordinary Differential Equations

Research Article Runge-Kutta Type Methods for Directly Solving Special Fourth-Order Ordinary Differential Equations Hndaw Publhng Corporaton Mathematcal Problem n Engneerng Volume 205, Artcle ID 893763, page http://dx.do.org/0.55/205/893763 Reearch Artcle Runge-Kutta Type Method for Drectly Solvng Specal Fourth-Order

More information

Outline and Reading. Dynamic Programming. Dynamic Programming revealed. Computing Fibonacci. The General Dynamic Programming Technique

Outline and Reading. Dynamic Programming. Dynamic Programming revealed. Computing Fibonacci. The General Dynamic Programming Technique Outlne and Readng Dynamc Programmng The General Technque ( 5.3.2) -1 Knapsac Problem ( 5.3.3) Matrx Chan-Product ( 5.3.1) Dynamc Programmng verson 1.4 1 Dynamc Programmng verson 1.4 2 Dynamc Programmng

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

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

Probabilistic & Unsupervised Learning

Probabilistic & Unsupervised Learning Probablstc & Unsupervsed Learnng Convex Algorthms n Approxmate Inference Yee Whye Teh ywteh@gatsby.ucl.ac.uk Gatsby Computatonal Neuroscence Unt Unversty College London Term 1, Autumn 2008 Convexty A convex

More information

T 11 ere we are aumng that the node label, 1, 2,, k all belong to or all belong to B, and fl( ; 1 ; 2 ;:::; k ) denote operad compoton. Th equvalence

T 11 ere we are aumng that the node label, 1, 2,, k all belong to or all belong to B, and fl( ; 1 ; 2 ;:::; k ) denote operad compoton. Th equvalence Tenor Product of 1 Operad Z. Fedorowcz Defnton. Let and B be operad actng on a topologcal pace X. We ay that the acton of and B are nterchangeable f for any k and any element ff 2B(k) the nduced map ff

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

The Second Anti-Mathima on Game Theory

The Second Anti-Mathima on Game Theory The Second Ant-Mathma on Game Theory Ath. Kehagas December 1 2006 1 Introducton In ths note we wll examne the noton of game equlbrum for three types of games 1. 2-player 2-acton zero-sum games 2. 2-player

More information

Modelli Clamfim Equazione del Calore Lezione ottobre 2014

Modelli Clamfim Equazione del Calore Lezione ottobre 2014 CLAMFIM Bologna Modell 1 @ Clamfm Equazone del Calore Lezone 17 15 ottobre 2014 professor Danele Rtell danele.rtell@unbo.t 1/24? Convoluton The convoluton of two functons g(t) and f(t) s the functon (g

More information

Vector Norms. Chapter 7 Iterative Techniques in Matrix Algebra. Cauchy-Bunyakovsky-Schwarz Inequality for Sums. Distances. Convergence.

Vector Norms. Chapter 7 Iterative Techniques in Matrix Algebra. Cauchy-Bunyakovsky-Schwarz Inequality for Sums. Distances. Convergence. Vector Norms Chapter 7 Iteratve Technques n Matrx Algebra Per-Olof Persson persson@berkeley.edu Department of Mathematcs Unversty of Calforna, Berkeley Math 128B Numercal Analyss Defnton A vector norm

More information

Deriving the X-Z Identity from Auxiliary Space Method

Deriving the X-Z Identity from Auxiliary Space Method Dervng the X-Z Identty from Auxlary Space Method Long Chen Department of Mathematcs, Unversty of Calforna at Irvne, Irvne, CA 92697 chenlong@math.uc.edu 1 Iteratve Methods In ths paper we dscuss teratve

More information

On the SO 2 Problem in Thermal Power Plants. 2.Two-steps chemical absorption modeling

On the SO 2 Problem in Thermal Power Plants. 2.Two-steps chemical absorption modeling Internatonal Journal of Engneerng Reearch ISSN:39-689)(onlne),347-53(prnt) Volume No4, Iue No, pp : 557-56 Oct 5 On the SO Problem n Thermal Power Plant Two-tep chemcal aborpton modelng hr Boyadjev, P

More information

Linear Feature Engineering 11

Linear Feature Engineering 11 Lnear Feature Engneerng 11 2 Least-Squares 2.1 Smple least-squares Consder the followng dataset. We have a bunch of nputs x and correspondng outputs y. The partcular values n ths dataset are x y 0.23 0.19

More information

Strong Efficient Domination in Graphs

Strong Efficient Domination in Graphs P P P IJISET - Internatonal Journal of Innovatve Scence, Engneerng & Technology, Vol Iue 4, June 04 wwwjetcom Strong Effcent Domnaton n Graph ISSN 348-7968 3 NMeenaP P, ASubramananP P, VSwamnathanP PDepartment

More information

An efficient algorithm for multivariate Maclaurin Newton transformation

An efficient algorithm for multivariate Maclaurin Newton transformation Annales UMCS Informatca AI VIII, 2 2008) 5 14 DOI: 10.2478/v10065-008-0020-6 An effcent algorthm for multvarate Maclaurn Newton transformaton Joanna Kapusta Insttute of Mathematcs and Computer Scence,

More information

10-801: Advanced Optimization and Randomized Methods Lecture 2: Convex functions (Jan 15, 2014)

10-801: Advanced Optimization and Randomized Methods Lecture 2: Convex functions (Jan 15, 2014) 0-80: Advanced Optmzaton and Randomzed Methods Lecture : Convex functons (Jan 5, 04) Lecturer: Suvrt Sra Addr: Carnege Mellon Unversty, Sprng 04 Scrbes: Avnava Dubey, Ahmed Hefny Dsclamer: These notes

More information

Outline. Bayesian Networks: Maximum Likelihood Estimation and Tree Structure Learning. Our Model and Data. Outline

Outline. Bayesian Networks: Maximum Likelihood Estimation and Tree Structure Learning. Our Model and Data. Outline Outlne Bayesan Networks: Maxmum Lkelhood Estmaton and Tree Structure Learnng Huzhen Yu janey.yu@cs.helsnk.f Dept. Computer Scence, Unv. of Helsnk Probablstc Models, Sprng, 200 Notces: I corrected a number

More information

The Study of Teaching-learning-based Optimization Algorithm

The Study of Teaching-learning-based Optimization Algorithm Advanced Scence and Technology Letters Vol. (AST 06), pp.05- http://dx.do.org/0.57/astl.06. The Study of Teachng-learnng-based Optmzaton Algorthm u Sun, Yan fu, Lele Kong, Haolang Q,, Helongang Insttute

More information

On Partial Optimality in Multi-label MRFs

On Partial Optimality in Multi-label MRFs On Partial Optimality in Multi-label MRFs P. Kohli 1 A. Shekhovtsov 2 C. Rother 1 V. Kolmogorov 3 P. Torr 4 1 Microsoft Research Cambridge 2 Czech Technical University in Prague 3 University College London

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

Predictors Using Partially Conditional 2 Stage Response Error Ed Stanek

Predictors Using Partially Conditional 2 Stage Response Error Ed Stanek Predctor ng Partally Condtonal Stage Repone Error Ed Stane TRODCTO We explore the predctor that wll relt n a mple random ample wth repone error when a dfferent model potlated The model we decrbe here cloely

More information

Extended Prigogine Theorem: Method for Universal Characterization of Complex System Evolution

Extended Prigogine Theorem: Method for Universal Characterization of Complex System Evolution Extended Prgogne Theorem: Method for Unveral Characterzaton of Complex Sytem Evoluton Sergey amenhchkov* Mocow State Unverty of M.V. Lomonoov, Phycal department, Rua, Mocow, Lennke Gory, 1/, 119991 Publhed

More information

VARIATION OF CONSTANT SUM CONSTRAINT FOR INTEGER MODEL WITH NON UNIFORM VARIABLES

VARIATION OF CONSTANT SUM CONSTRAINT FOR INTEGER MODEL WITH NON UNIFORM VARIABLES VARIATION OF CONSTANT SUM CONSTRAINT FOR INTEGER MODEL WITH NON UNIFORM VARIABLES BÂRZĂ, Slvu Faculty of Mathematcs-Informatcs Spru Haret Unversty barza_slvu@yahoo.com Abstract Ths paper wants to contnue

More information

A Preliminary Study on Material Utilization of Stiffened Cylindrical Shells

A Preliminary Study on Material Utilization of Stiffened Cylindrical Shells Reearch Journal of Appled Scence, Engneerng and echnology 6(5): 757-763, 03 ISSN: 040-7459; e-issn: 040-7467 Maxwell Scentfc Organzaton, 03 Submtted: December 8, 0 Accepted: February 08, 03 Publhed: Augut

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

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system Transfer Functons Convenent representaton of a lnear, dynamc model. A transfer functon (TF) relates one nput and one output: x t X s y t system Y s The followng termnology s used: x y nput output forcng

More information