Mixed-integer vertex covers on bipartite graphs

Size: px
Start display at page:

Download "Mixed-integer vertex covers on bipartite graphs"

Transcription

1 Mxed-nteger vertex covers on bpartte graphs Mchele Confort, Bert Gerards, Gacomo Zambell November, 2006 Abstract Let A be the edge-node ncdence matrx of a bpartte graph G = (U, V ; E), I be a subset the nodes of G, and b be a vector such that 2b s ntegral. We consder the followng mxed-nteger set: X(G, b, I) = {x : Ax b, x 0, x nteger for all I}. We characterze conv(x(g, b, I)) n ts orgnal space. That s, we descrbe a matrx (A, b ) such that conv(x(g, b, I)) = {x : A x b }. Ths s accomplshed by computng the projecton onto the space of the x-varables of an extended formulaton, gven n [1], for conv(x(g, b, I)). We then gve a polynomal algorthm for the separaton problem for conv(x(g, b, I)), thus showng that the problem of optmzng a lnear functon over the set X(G, b, I) s polynomally solvable. 1 Introducton 1.1 The problem Gven a bpartte graph G = (U, V ; E), a vector b = (b e ) e E, wth the property that b s half-ntegral,.e. 2b e Z, e E, and a set I (U V ), we consder the problem of characterzng the convex hull of all nonnegatve x R V such that x + x j b for every E, x Z for every I. That s, gven the edge-node ncdence matrx A of a bpartte graph G, a partton (I, L) of ts column-set, and an half-ntegral vector b, we consder the followng mxed-nteger set: X(G, b, I) = {x : Ax b, x 0, x nteger for all I}. (1) Dpartmento d Matematca Pura e Applcata, Unverstá d Padova, Va Treste 63, Padova, Italy. Centrum voor Wskunde en Informatca, Kruslaan 413, 1098 SJ Amsterdam, The Netherlands & Technsche Unverstet Endhoven, Den Dolech 2, Endhoven 1

2 In ths paper we provde a formulaton for the polyhedron conv(x(g, b, I)). A formulaton for a polyhedron P (n ts orgnal space) s a descrpton of P as the ntersecton of a fnte number of half-spaces. So t conssts of a fnte set of nequaltes Cx d such that P = {x : Cx d}. In [1] a general technque was ntroduced to descrbe an extended formulaton for the set of solutons of a system Ax b, when A s a network matrx and some of the varables are restrcted to be nteger. A formulaton of P s extended whenever t defnes a polyhedron P n a hgher dmensonal space the ncludes the orgnal space, so that P s the projecton of ths polyhedral descrpton onto the orgnal space. In Secton 2 we derve the extended formulaton for conv(x(g, b, I)), whle n Secton 3 we descrbe a formulaton n the orgnal space by explctly computng the projecton of the polyhedron defned by the extended formulaton. Fnally, n Secton 4, we gve a polynomal tme algorthm to solve the separaton problem for conv(x(g, b, I)). 1.2 The man result Gven a bpartte graph G = (U, V ; E), a partton (I, L) of U V and an half-ntegral vector b, we say that a path P of G s an I-path f at least one endnode of P s n I, and no ntermedate node of P s n I. We say that P s odd f P has an odd number of edges e such that b e = 1 2 mod 1. In ths paper we show the followng: Theorem 1 The polyhedron conv(x(g, b, I)) s defned by the followng nequaltes: x + x j b E (2) 2x(V (P ) L) + x(v (P ) I) b(p ) P odd I-path (3) x 0 V (4) Esenbrand [2] conjectured that the nequaltes n (2)-(4) are suffcent to characterze conv(x(g, b, I)) when G s a path. So Theorem 1 shows that ths conjecture holds n a qute more general settng (and t certanly cannot be extended beyond that). Prelmnary results for the path case were obtaned by Skutella [9] and Esenbrand [2]. 1.3 Frst Chvátal closure The followng observaton allows us to descrbe X(G, b, I) n terms of a pure nteger set. Observaton 2 Let x be a vertex of conv(x(g, b, I)). Then 2 x s ntegral. 2

3 Proof: If not, let U and V be the sets of nodes n U and V, respectvely, such that 2 x s not nteger. Then, for ɛ small enough, the vectors x+ɛχ U ɛχ V and x ɛχ U + ɛχ V are both n conv(x(g, b, I)), where we denote by χ S the ncdence vector of S for any S U V. Let b = 2b, A be obtaned form A by multplyng by 2 the columns correspondng to nodes n I. By Observaton 2, the lnear transformaton x = x, I, x = 2x, L, maps X(G, b, I) nto {x : A x b, x 0, x ntegral} whch s a pure nteger set. Let P = v 1,... v n be an I-path. Notce that b(p ) = 1 2 mod 1 s equvalent to b (P ) odd. Then the nequalty b x (P ) (5) 2 V (P ) s a Gomory-Chvátal nequalty of {x : A x b, 0}. Indeed, assume v 1 I. If v n I, then (5) s obtaned from n (2x v 1 +x v 2 b 1 v 1 v 2 )+ 2 (x v +x v +1 b v v +1 )+ 1 2 (x v n 1 +2x v n b v n 1 v n ) =2 by roundng up the rght-hand-sde. If x n / I, then (5) s obtaned from n (2x v 1 + x v 2 b 1 v 1 v 2 ) + 2 (x v + x v +1 b v v +1 ) (x v n 0) =2 by roundng up the rght-hand-sde. Furthermore the nequaltes n (5) are equvalent to the nequaltes n (3). 1.4 The motvaton A (general) mxed-nteger set s a set of the form s {x Ax b, x nteger I} (6) where I s a subset of the columns of A and b s a vector that may contan fractonal components. In [1], t s shown that the problem of decdng f the above set s nonempty s NP-complete, even f b s an half-ntegral vector and A s a network matrx. (We refer the reader to [5] or [8] for defntons and results related to network matrces and, more generally, totally unmodular matrces.) However, t may be possble that, when A s the transpose of a network matrx, the assocated mxed-nteger programmng problem s polynomally solvable. Indeed, let MIX 2T U be a mxed-nteger set of the form (6) when 3

4 A s a network matrx. An extended formulaton of the polyhedron conv(mix 2T U ) was descrbed n [1]. The extended formulaton nvolves an addtonal varable for each possble fractonal parts taken by the varables at any vertex of conv(mix 2T U ). If ths number s polynomal n the sze of (A, b), t s shown n [1] that the formulaton s compact,.e. of polynomal sze n the sze of (A, b). Therefore the problem of optmzng a lnear functon over MIX 2T U can be effcently solved n ths case. However, t seems to be rather dffcult to compute the projecton n the orgnal x-space. It follows from Observaton 2 that f x s a vertex of conv(x(g, b, I)), then x x {0, 1 2 }. Therefore the extended formulaton for conv(x(g, b, I)) (whch wll be ntroduced n Secton 2) s compact. The man contrbuton of ths paper s the explct descrpton of the projecton of the polyhedron defned by ths extended formulaton n the orgnal x-space. The mxed-nteger set X(G, b, I) s related to some mxed-nteger sets that arse n the context of producton plannng (see [7]). The case when G s a star wth center node n L and leaves n I has been studed by Pochet and Wolsey n [6], where they gave an extended formulaton for the convex hull of feasble solutons whch s compact. Günlük and Pochet [3] projected ths formulaton onto the orgnal space, thus showng that the famly of mxng nequaltes gves the formulaton n the x-space. Mller and Wolsey [4] extended the results n [6] to general bpartte graphs, wth the restrcton that the partton (I, L) corresponds to the bpartton (U, V ) of the graph. Ther result shows that the mxng nequaltes assocated wth every sngle star of G havng center a node n L and leaf nodes all nodes n I gve a formulaton for ths case. 2 The extended formulaton We use here a modelng technque ntroduced by Pochet and Wolsey [6] and extensvely nvestgated n [1]. Observaton 2 allows to express each varable n L as x = µ δ, µ 0 nteger, 0 δ 1 nteger. (7) For now, we assume I =, that s, L = (U V ). Lemma 3 Let E, and suppose x, x j satsfy (7). If b = 1 2 mod 1, x, x j satsfy x + x j b f and only f µ + µ j b µ + δ + µ j + δ j b. If b = 0 mod 1, x, x j satsfy x + x j b f and only f µ + δ + µ j b µ + µ j + δ j b. (8) (9) 4

5 Proof. Assume x, x j satsfy (7). Then, f b = 1 2 mod 1, constrant x + x j b s satsfed f and only f µ + µ j b and δ + δ j = 1 whenever µ + µ j = b. If b = 0 mod 1, the constrant s satsfed f and only f µ + µ j b 1 and δ = δ j = 1 whenever µ + µ j = b 1. It s easy to see that these two condtons are modeled by the above constrants. Observaton 4 Gven E, the constrants (8) and (9) belong to the frst Chvátal closure of the polyhedron defned by µ δ + µ j δ j b µ, µ j 0 δ, δ j 1 δ, δ j 0 whenever b = 1 2 mod 1 and b = 0 mod 1, respectvely. By applyng the unmodular transformaton µ 0 = µ, µ 1 = µ + δ, the constrants x = µ δ, µ 0, 0 δ 1 become x 1 2 (µ0 + µ 1 ) = 0 (10) and constrants (8) and (9) become: µ µ 1 µ0 1 µ 0 + µ0 j b µ 1 + µ1 j b (11) (12) µ 1 + µ0 j b µ 0 + µ1 j b (13) Theorem 5 The projecton onto the space of the x varables of the polyhedron Q defned on the space of the varables (x, µ 0, µ 1 ) by the nequaltes (10), (11) for every U V, (12) for every E s.t. b = 1 2 mod 1 (13) for every E s.t. b = 0 mod 1 s the polyhedron conv(x(g, b, )). Proof: Snce the varable x s determned by (10) for all U V, we only need to show that the polyhedron defned by nequaltes (11) for every U V, (12) for every E s.t. b mod 1, and (13) for every 5 = 1 2

6 E s.t. b = 0 mod 1, s ntegral. Let A µ be the constrant matrx of the above system. Snce G s a bpartte graph, then the matrx Ā, obtaned by multplyng by 1 the columns of A µ relatve to the varables µ 0, µ1, V, has at most a 1 and at most a 1 n each row. Therefore Ā s the transpose of a network matrx, so A µ s totally unmodular (see [8]). Snce the left-hand-sdes of (11)-(13) are all nteger, the statement follows from the theorem of Hoffman and Kruskal. Observaton 6 Varable x s nteger valued f and only f δ = 0, U V. Therefore, for a gven I (U V ), the polyhedron conv(x(g, b, I)) s the projecton on the space of the x varables of the face Q I of Q defned by the equatons µ 1 1 µ0 = 0, I (whch correspond to δ = 0, I). 3 The formulaton n the orgnal space In ths secton we prove Theorem 1 by projectng onto the x-space the polyhedron Q I. Let p = µ0 µ1 2. The µ 0 = x + p and µ 1 = x p. The nequaltes defnng Q become: p + p j b x x j, b = 1 2 mod 1 p p j b x x j, b = 1 2 mod 1 p p j b x x j, b = 0 mod 1 p + p j b x x j, b = 0 mod 1 p 1 2 U V p 0 U V p x U V By Observaton 6, conv(x(g, B, I)) s the projecton onto the x-space of the polyhedron defned by the above nequaltes and by p = 0 for every I. Assocate multplers to the above constrants as follows: u ++ p + p j b x x j u p p j b x x j u + p p j b x x j u + p + p j b x x j u 1 2 p 1 2 u 0 p 0 u x p x (14) Any vald nequalty for conv(x(g, b, I)) has the form α u x β u, where α u x = + u )(x + x j ) + (u ++ b = 1 2 mod 1 6

7 b =0 mod 1 (u + + u + )(x + x j ) + u x x (15) β u = (u b + u ++ b ) + b = 1 2 mod 1 (u + + u + )b 1 2 u 1 2 (16) L b =0 mod 1 for some nonnegatve vector u = (u ++, u, u +, u +, u 1 2, u 0, ux ) such that up = 0, where P s the column-submatrx of the above system nvolvng columns correspondng to varables p, L (see e.g. Theorem 4.10 n [5]). For nstance the nequalty x + x j b where b = 1 2 mod 1 s obtaned by settng u ++ = u = 1 2, and all other entres of u to be 0. We are nterested n characterzng the nonnegatve vectors u such that up = 0 and α u x β u s facet defnng for conv(x(g, b, I)), and such that the nequalty α u x β u s not of the form x + x j b, for some E, or x 0, for some U V. From now on we wll assume, w.l.o.g., that the entres of u are nteger and relatvely prme. We defne an auxlary graph Γ u = (L {d}, F ), where d s a dummy node not n U V, and F u s defned as follows. For every edge E such that, j L, there are u ++ + u + u + parallel edges between and j n F, each edge correspondng to u + a varable among u ++, u, u +, u +. + For each node L, there are u u 0 + ux + j I : E (u++ + u + u + + u + ) parallel edges between d and n F, each edge correspondng to a varable among u 1 2, u 0, ux u +, for some j I., or u++, u, u +, We mpose a b-orentaton ω on Γ, that s, to each edge e F, and each endnode of e that belongs to L, we assocate the value ω(e, ) = tal f e corresponds to an nequalty of (14) where p has coeffcent 1, whle we assocate the value ω(e, ) = head f e corresponds to an nequalty of (14) where p has coeffcent +1. The dummy node d s nether a tal nor a head of any edge. Thus, each edge of Γ u can have one head and one tal, two heads, two tals, or, f d s one of the two endnodes, only one head and no tal or only one tal and no head. For each L, we denote wth δω n () the number of edges n F u of whch s a head, and wth δ out w () the number of edges n F of whch s a tal. ω () for every L. We say that Γ u s ω-euleran f δ n ω () = δ out 7

8 Observaton 7 Γ u s ω-euleran f and only f up = 0. We defne a closed ω-euleran walk n Γ u as a closed-walk n Γ u, v 0, e 0, v 1, e 1,..., v k, e k, v k+1, where v 0 = v k+1, wth the property that ω(e h 1, v h ) ω(e h, v h ) for every h such that v h s n L, h = 0,..., k, k + 1, where the ndces are taken modulo k. That s, f v h L, then v h s a head of e h 1 f and only f v h s a tal of e h. Observaton 8 Γ u s ω-euleran f and only f Γ u s the dsjont unon of closed ω-euleran walks. In partcular, every node n L {d} has even degree n Γ u. Observe that, f v 0, e 0,..., e k, v k+1 s a closed ω-euleran walk n Γ u, then both graphs Γ, Γ on L {d} wth edge-sets F = {e 1,..., e k } and F = F \ F, respectvely, are ω-euleran. Suppose F. Then there are nonnegatve nteger vectors u and u, both dfferent from zero, such that u P = 0, u P = 0, Γ = Γ u and Γ = Γ u, and u = u + u. By the fact that Γ and Γ are ω-euleran, and by the structure of the nequaltes n (14), the vectors (α u, β u ) and (α u, β u ) are both non-zero. Furthermore α u = α u + α u and β u = β u + β u, contradctng the fact that α u x β u s facet defnng and the entres of u are relatvely prme. Hence we have shown the followng. Observaton 9 Every closed ω-euleran walk of Γ u traverses all the edges n F. In partcular, there exsts a closed ω-euleran walk v 0, e 0,..., e k, v k+1 of Γ u such that F = {e h h = 1,..., k}. Suppose d has postve degree n Γ. Then we may assume, w.l.o.g., that v 0 = v k+1 = d. Suppose d = v h for some h = 1,..., k. Then v 0, e 0, v 1,..., e h 1 v h s a closed ω-euleran walk, contradctng the prevous observaton. Hence we have the followng. Observaton 10 Node d has degree 0 or 2 n Γ u. Next we show the followng. Lemma 11 Every node n L {d} has degree 0 or 2 n Γ u. Proof: We have already shown d has degree 0 or 2 n Γ u. If d has degree 2, we assume d = v 0 = v k+1, else v 0 s arbtrarly chosen. If there s a node n L wth degree at least 4, then there exsts dstnct ndces s, t {1,..., k} such that v s = v t. We choose s and t such that t s s postve and as small 8

9 as possble. Therefore C = v s, e s,..., e t 1, v t s a cycle of Γ u contanng only nodes n L. Snce G s a bpartte graph, C has even length, hence the edges n C can be parttoned nto two matchngs M 0, M 1 of cardnalty C /2. We wll denote wth HH, T T, HT the sets of edges of F wth, respectvely, two heads, two tals, one head and one tal. If v s s the head of exactly one among e s and e t 1, then C s a closed ω-euleran walk, contradctng Observaton 9. Hence v s s ether a head of both e s and e t 1 or a tal of both e s and e t 1. Ths shows that C T T = C HH ± 1. Therefore there s an odd number of edges e n C such that b e = 1 2 mod 1. By symmetry, we may assume e M 0 b e e M 1 b e Then the nequalty 2 x b e + 1 (17) 2 V (C) e C s vald for conv(x(g, b, I)), snce t s mpled by the vald nequaltes x + x j b, M 0, because 2 V (C) x = 2 (x +x j ) 2 b b e + b e = b e M 0 M 0 e M 0 e M 1 e C Case 1: Node v s s a tal of both e s and e t 1. Then C T T = C HH + 1, hence b e + b e + b e = b e (18) e C T T e C HH e C HT e C Let u be the vector obtaned from u as follows { u = u 1 u 0 v s = u 0 v s + 2 for every C all other components of u and u beng dentcal, where u among u ++, u, u +, u + correspondng to edge of C. s the varable Then one can easly see that Γ u s the graph obtaned from Γ u by removng the edges e s,..., e t, and addng two parallel edges v s d both wth tal n v s, hence Γ u s ω-euleran and u P = 0. By (18) β u = β u e C b e 1 2, whle by constructon α u x = α u x + 2 x. V (C) 9

10 Thus α u x β u can be obtaned by takng the sum of α u x β u and (17), contradctng the assumpton that α u x β u s facet defnng. Case 2: Node v s s a head of both e s and e t 1. Then C T T = C HH 1, hence e C T T b e + e C HH b e + e C HT Let u be the vector obtaned from u as follows b e = e C b e 1 2. (19) { u = u 1 u 1 2 vs = u 0 v s + 2 for every C all other components of u and u beng dentcal. Then one can easly see that Γ u s the graph obtaned from Γ u by removng the edges e s,..., e t, and addng two parallel edges v s d both wth head n v s, hence u P = 0. By (19) β u = β u e C b e , whle by constructon α u x = α u x + 2 x. V (C) Thus α u x β u can be obtaned by takng the sum of α u x β u and (17), contradctng the assumpton that α u x β u s facet defnng. We are now ready to gve the proof of the man theorem. Proof of Theorem 1. We show that all facet defnng nequaltes α u x β u, where u s nonnegatve, ntegral, and wth entres that are relatvely prme, that are not nequaltes n (2) or (4), are of the form (3). Frst we show the followng. E u > E u ++ In fact, we can wrte the nequalty α u x (u + u ++ )b + b = 1 2 mod V b =0 mod 1 u 1 2 (20) (u + + u + )b

11 as nonnegatve combnaton of nequaltes of the form (2) or (4), therefore we must have β u > (u + u ++ )b + (u + + u + )b. b = 1 2 mod 1 b =0 mod 1 Thus 0 < β u (u + u ++ )b + b = 1 2 mod 1 = 1 2 ( E whch proves (20). u E u ++ V u 1 2 ) b =0 mod 1 (u + + u + By Lemma (11) and Observaton (9), Γ u conssts of an nduced cycle C and solated nodes, where every node n V (C) L s a head of exactly one edge and a tal of exactly one edge. If d s an solated node, then each edge of C corresponds to a varable of the form u, and snce the total number of heads n C equals the number of tals, then E u = E u++ and V u 1 2 = 0, contradctng (20). Thus we may assume that C = v 0, e 0,..., e k, v k+1 where d = v 0 = v k+1. Clam: The followng are the only possble cases. 1. Edges dv 1, dv k of Γ u correspond to varables u x v 1 and u x v k, respectvely; 2. dv 1 corresponds to varable u or u + for some w I, and dv k corresponds to u x v k ; 3. dv 1 corresponds to varables u or u + for some w I, and dv k corresponds to varable u w v k or u + w v k for some w I. Proof of clam If v 1 s a head of e 0 and v k s a head of e k, then the number of edges among e 1,..., e k 1 wth two tals s one plus the number of edges wth two heads. Snce the former correspond to varables of type u for some E, and the latter correspond to to varables of type u ++ for some E, then by (20) dv 1 does not correspond to varable u 1 2 v1 )b or to a varable u ++ for any w I, and dv k does not correspond to varable u 1 2 vk or to a varable u ++ wv k for any w I, thus one of the above three cases holds. If v 1 s a tal of e 0 and v k s a head of e k, then the number of edges among e 1,..., e k 1 wth two tals s equal the number of edges wth two heads. By (20), dv 1 corresponds to varable u for some w I, and dv k corresponds to ether u x v k or to a varable u + w v k for some w I, thus case 2 or 3 holds. If v 1 s a tal of e 0 and v k s a tal of e k, then the number of edges among e 1,..., e k 1 wth two tals s equal one mnus the number of edges wth two 11

12 heads. By (20), dv 1 corresponds to varable u for some w I, and dv k corresponds to a varable u w v k for some w I, thus case 3 holds. Ths completes the proof of the clam. Case 1: Edges dv 1, dv k of Γ u correspond to varables u x v 1 and u x v k, respectvely. In ths case the path P = v 1, e 1,..., e k 1, v k of Γ u s also a path of G contanng only nodes n L, and P contans an odd number of edges e such that b e = 1 2 mod 1. The nequalty α ux β u s then 2x(V (P )) b(p ) The edges of P can be parttoned nto two matchngs M 0 and M 1, thus we may assume, w.l.o.g., e M 0 b e e M 1 b e Thus 2x(V (P )) 2 M 0 (x + x j ) 2 M 0 b e M 0 b e + e M 1 b e = b(p ) + 1 2, hence α u x β u s not facet defnng. Case 2: dv 1 corresponds to varable u or u + corresponds to u x v k. for some w I, and dv k In ths case, P = w, v 1, e 1,..., e k 1, v k s an odd I-path of G between w I and v k L. The nequalty α u x β u s 2x(V (P ) L)+x w b(p )+ 1 2, whch s one of the nequaltes n (3). Case 3: dv 1 corresponds to varables u or u + for some w I, and dv k corresponds to varable u w v k or u + w v k for some w I. If w w, then the path P = w, v 1, e 1,..., e k 1, v k, w s an odd I-path of G between w I and w I. The nequalty α u x β u s 2x(V (P ) L) + x w + x w b(p ) + 1 2, whch s one of the nequaltes n (3). If w = w, then we must have v 1 v k, snce otherwse v 1 would be ether the head or the tal of both edges of Γ u ncdent to v 1. Thus C = w, v 1,..., v k, w s a cycle of G. Snce G s a bpartte graph, C has even length, hence the edges n C can be parttoned nto two matchngs M 0, M 1 of cardnalty C /2. Snce C contans an odd number of edges e such that b w = 1 2 mod 1, then we may assume, w.l.o.g., e M 0 b e e M 1 b e The nequalty α u x β u s 2x(V (C )) b(c ) But 2x(V (C )) = 2 M 0 (x + x j ) 2 M 0 b e M 0 b e + e M 1 b e = b(c ) + 1 2, hence α u x β u s not facet defnng. 4 Separaton Theorem 5 and Observaton 6 mply that the problem of mnmzng a lnear functon over the set X(G, b, I) s solvable n polynomal tme, snce t reduces to solvng a lnear programmng problem over the set of feasble ponts for (10)-(13). In ths secton we gve a combnatoral polynomal tme algorthm for the separaton problem for the set conv(x(g, b, I)), thus gvng an alternatve 12

13 proof that the problem of optmzng a lnear functon over such polyhedron, and thus over X(G, b, I), s polynomal. Clearly, gven a nonnegatve vector x, we can check n polynomal-tme whether x satsfes (2) for every edge. Thus, by Theorem 1, we only need to descrbe a polynomal tme algorthm that, gven a nonnegatve vector x satsfyng (2), ether returns an nequalty of type (3) volated by x, or proves that none exsts. For every E, let s = x + x j b. Snce x satsfes (2), then s e s nonnegatve for every e E. Let P = v 1,... v n be an odd I-path. Clam 12 The vector x satsfes 2x (V (P ) L) + x (V (P ) I) b(p ) f and only f s (P ) + x ({v 1, v n } L) 1 2. Indeed, assume v 1 I. If v n I then n 1 n 1 s v v +1 = (x v + x v +1 b v v +1 ) =1 =1 gves the equalty s (P ) = 2x (V (P ) L) + x (V (P ) I) b(p ), hence 2x (V (P ) L) + x (V (P ) I) b(p ) f and only f s (P ) 1 2. If v n / I, then n 1 n 1 s v v +1 + x v n = (x v + x v +1 b v v +1 ) + x v n =1 =1 gves the equalty s (P )+x v n = 2x (V (P ) L)+x (V (P ) I) b(p ), hence 2x (V (P ) L) + x (V (P ) I) b(p ) f and only f s (P ) + x v n 1 2. Ths completes the proof of the Clam. Therefore, f we assgn length s e to every e E, we need to gve an algorthm that, for any two nodes r, t such that r I, ether determnes that the shortest odd I-path between r and t (f any) has length at least 1 2 x ({t} L), or returns an odd I-path P for whch 2x (V (P ) L) + x (V (P ) I) < b(p ) Observe that any walk W between r and t that contans an odd number of edges e such that b e = 1 2 mod 1 ether contans a sub-path P that s an odd I-path or t contans a cycle C that contans an odd number of edges e such that b e = 1 2 mod 1. In the former case, ether both endnodes of P are n I, or t s the only endnode of P n L. Hence, f s (W ) < 1 2 x ({t} L), then also s (P ) < 1 2 x ({t} L), hence 2x (V (P ) L)+x (V (P ) I) < b(p ) In the second case, snce G s bpartte, the edges of C can be parttoned nto two matchngs M 0 and M 1 such that b(m 0 ) b(m 1 ) Thus s (C) = C (x + x j b ) = 2x (V (C)) b(c) 2(x (V (C)) b(m 0 )) = 2 M 0 (x + x j b ) , hence s (W )

14 Thus we only need to fnd, for every par r, t V wth r I, the shortest walk W between r and t, w.r.t. the dstance s, among all such walks contanng an odd number of edges e such that b e = 1 2 mod 1. If, for a gven choce of r, t, s(w ) < 1 2 x ({t} L), then by the above argument we can fnd n polynomal tme a sub-path P of W such that P s an odd I-path and 2x (V (P ) L) + x (V (P ) I) < b(p ) + 1 2, otherwse we can conclude that x conv(x(g, b, I)). To conclude, we only need to show a polynomal tme algorthm that, gven an undrected graph Γ wth nonnegatve lengths on the edges l e, e E(Γ), a subset F E(Γ), and a par of nodes r, t V (Γ), determnes the walk W of mnmum length between r and t such that E(W ) F s odd, or determnes that no such walk exsts. Notce that ths problem easly reduces to the case where F = E(Γ), snce we can construct a graph Γ by subdvdng each edge uv E(Γ) \ F nto the path u, w, v, where w s a new node, and assgn lengths to uw and wv n such a way that the sum of such lengths equals l uv. Clearly, a walk W between two nodes r and t n Γ contans an odd number of edges n F f and only f the correspondng walk W n Γ has an odd number of edges. Furthermore, W and W have the same length. By the prevous argument, we are nterested n the problem of fndng a shortest walk wth an odd number of edges between a gven par of nodes. Ths problem can be solved n polynomal tme. Snce, as far as we know, ths fact s folklore, we brefly descrbe an algorthm. We construct a new graph Γ as follows. For every node v V (Γ), we have a par of nodes v, v n V (Γ ). For every edge uv E(Γ), we have two edges uv and u v n Γ, both wth length l uv. One can verfy that a walk W wth an odd number of edges between r and t exsts n Γ f and only f there exsts a walk of the same length between r and t n Γ. Hence we only need to fnd a shortest path between r and t n Γ, f any exsts, and output the correspondng walk n Γ. References [1] M. Confort, M. D Summa, F. Esenbrand, L.A. Wolsey, Network formulatons of mxed-nteger programs, In preparaton, [2] F. Esenbrand, Mxed Integer Programmng over TU systems, Manuscrpt, [3] O. Günlük and Y. Pochet, Mxng mxed nteger nequaltes, Mathematcal Programmng 90 (2001), [4] A. Mller and L.A. Wolsey, Tght formulatons for some smple MIPs and convex objectve IPs, Mathematcal Programmng B 98 (2003),

15 [5] G.L. Nemhauser, L.A. Wolsey, Integer and Combnatoral Optmzaton, Wley Interscence, New York, [6] Y. Pochet and L.A. Wolsey, Polyhedra for lot-szng wth Wagner-Whtn costs, Mathematcal Programmng 67 (1994), [7] Y. Pochet and L.A. Wolsey, Producton Plannng by Mxed Integer Programmng, Sprnger Seres n Operatons Research and Fnancal Engneerng, New York, [8] A. Schrver, Theory of Lnear and Integer Programmng, Wley, New York, [9] M. Skutella, Mxed Integer vertex cover on paths, Manuscrpt,

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

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

= 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

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

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

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

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

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

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

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

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

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

ISSN: ISO 9001:2008 Certified International Journal of Engineering and Innovative Technology (IJEIT) Volume 3, Issue 1, July 2013

ISSN: ISO 9001:2008 Certified International Journal of Engineering and Innovative Technology (IJEIT) Volume 3, Issue 1, July 2013 ISSN: 2277-375 Constructon of Trend Free Run Orders for Orthogonal rrays Usng Codes bstract: Sometmes when the expermental runs are carred out n a tme order sequence, the response can depend on the run

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

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

2.3 Nilpotent endomorphisms

2.3 Nilpotent endomorphisms s a block dagonal matrx, wth A Mat dm U (C) In fact, we can assume that B = B 1 B k, wth B an ordered bass of U, and that A = [f U ] B, where f U : U U s the restrcton of f to U 40 23 Nlpotent endomorphsms

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

FACTORIZATION IN KRULL MONOIDS WITH INFINITE CLASS GROUP

FACTORIZATION IN KRULL MONOIDS WITH INFINITE CLASS GROUP C O L L O Q U I U M M A T H E M A T I C U M VOL. 80 1999 NO. 1 FACTORIZATION IN KRULL MONOIDS WITH INFINITE CLASS GROUP BY FLORIAN K A I N R A T H (GRAZ) Abstract. Let H be a Krull monod wth nfnte class

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

Computing Correlated Equilibria in Multi-Player Games

Computing Correlated Equilibria in Multi-Player Games Computng Correlated Equlbra n Mult-Player Games Chrstos H. Papadmtrou Presented by Zhanxang Huang December 7th, 2005 1 The Author Dr. Chrstos H. Papadmtrou CS professor at UC Berkley (taught at Harvard,

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

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

Volume 18 Figure 1. Notation 1. Notation 2. Observation 1. Remark 1. Remark 2. Remark 3. Remark 4. Remark 5. Remark 6. Theorem A [2]. Theorem B [2].

Volume 18 Figure 1. Notation 1. Notation 2. Observation 1. Remark 1. Remark 2. Remark 3. Remark 4. Remark 5. Remark 6. Theorem A [2]. Theorem B [2]. Bulletn of Mathematcal Scences and Applcatons Submtted: 016-04-07 ISSN: 78-9634, Vol. 18, pp 1-10 Revsed: 016-09-08 do:10.1805/www.scpress.com/bmsa.18.1 Accepted: 016-10-13 017 ScPress Ltd., Swtzerland

More information

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal Inner Product Defnton 1 () A Eucldean space s a fnte-dmensonal vector space over the reals R, wth an nner product,. Defnton 2 (Inner Product) An nner product, on a real vector space X s a symmetrc, blnear,

More information

Min Cut, Fast Cut, Polynomial Identities

Min Cut, Fast Cut, Polynomial Identities Randomzed Algorthms, Summer 016 Mn Cut, Fast Cut, Polynomal Identtes Instructor: Thomas Kesselhem and Kurt Mehlhorn 1 Mn Cuts n Graphs Lecture (5 pages) Throughout ths secton, G = (V, E) s a mult-graph.

More information

Lecture 5 Decoding Binary BCH Codes

Lecture 5 Decoding Binary BCH Codes Lecture 5 Decodng Bnary BCH Codes In ths class, we wll ntroduce dfferent methods for decodng BCH codes 51 Decodng the [15, 7, 5] 2 -BCH Code Consder the [15, 7, 5] 2 -code C we ntroduced n the last lecture

More information

The Multiple Classical Linear Regression Model (CLRM): Specification and Assumptions. 1. Introduction

The Multiple Classical Linear Regression Model (CLRM): Specification and Assumptions. 1. Introduction ECONOMICS 5* -- NOTE (Summary) ECON 5* -- NOTE The Multple Classcal Lnear Regresson Model (CLRM): Specfcaton and Assumptons. Introducton CLRM stands for the Classcal Lnear Regresson Model. The CLRM s also

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

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

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

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

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

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

SL n (F ) Equals its Own Derived Group

SL n (F ) Equals its Own Derived Group Internatonal Journal of Algebra, Vol. 2, 2008, no. 12, 585-594 SL n (F ) Equals ts Own Derved Group Jorge Macel BMCC-The Cty Unversty of New York, CUNY 199 Chambers street, New York, NY 10007, USA macel@cms.nyu.edu

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

Anti-van der Waerden numbers of 3-term arithmetic progressions.

Anti-van der Waerden numbers of 3-term arithmetic progressions. Ant-van der Waerden numbers of 3-term arthmetc progressons. Zhanar Berkkyzy, Alex Schulte, and Mchael Young Aprl 24, 2016 Abstract The ant-van der Waerden number, denoted by aw([n], k), s the smallest

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

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

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

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

a b a In case b 0, a being divisible by b is the same as to say that

a b a In case b 0, a being divisible by b is the same as to say that Secton 6.2 Dvsblty among the ntegers An nteger a ε s dvsble by b ε f there s an nteger c ε such that a = bc. Note that s dvsble by any nteger b, snce = b. On the other hand, a s dvsble by only f a = :

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

Exercise Solutions to Real Analysis

Exercise Solutions to Real Analysis xercse Solutons to Real Analyss Note: References refer to H. L. Royden, Real Analyss xersze 1. Gven any set A any ɛ > 0, there s an open set O such that A O m O m A + ɛ. Soluton 1. If m A =, then there

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

= = = (a) Use the MATLAB command rref to solve the system. (b) Let A be the coefficient matrix and B be the right-hand side of the system.

= = = (a) Use the MATLAB command rref to solve the system. (b) Let A be the coefficient matrix and B be the right-hand side of the system. Chapter Matlab Exercses Chapter Matlab Exercses. Consder the lnear system of Example n Secton.. x x x y z y y z (a) Use the MATLAB command rref to solve the system. (b) Let A be the coeffcent matrx and

More information

Stanford University Graph Partitioning and Expanders Handout 3 Luca Trevisan May 8, 2013

Stanford University Graph Partitioning and Expanders Handout 3 Luca Trevisan May 8, 2013 Stanford Unversty Graph Parttonng and Expanders Handout 3 Luca Trevsan May 8, 03 Lecture 3 In whch we analyze the power method to approxmate egenvalues and egenvectors, and we descrbe some more algorthmc

More information

Perfect Competition and the Nash Bargaining Solution

Perfect Competition and the Nash Bargaining Solution Perfect Competton and the Nash Barganng Soluton Renhard John Department of Economcs Unversty of Bonn Adenauerallee 24-42 53113 Bonn, Germany emal: rohn@un-bonn.de May 2005 Abstract For a lnear exchange

More information

MMA and GCMMA two methods for nonlinear optimization

MMA and GCMMA two methods for nonlinear optimization MMA and GCMMA two methods for nonlnear optmzaton Krster Svanberg Optmzaton and Systems Theory, KTH, Stockholm, Sweden. krlle@math.kth.se Ths note descrbes the algorthms used n the author s 2007 mplementatons

More information

5 The Rational Canonical Form

5 The Rational Canonical Form 5 The Ratonal Canoncal Form Here p s a monc rreducble factor of the mnmum polynomal m T and s not necessarly of degree one Let F p denote the feld constructed earler n the course, consstng of all matrces

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

FINITELY-GENERATED MODULES OVER A PRINCIPAL IDEAL DOMAIN

FINITELY-GENERATED MODULES OVER A PRINCIPAL IDEAL DOMAIN FINITELY-GENERTED MODULES OVER PRINCIPL IDEL DOMIN EMMNUEL KOWLSKI Throughout ths note, s a prncpal deal doman. We recall the classfcaton theorem: Theorem 1. Let M be a fntely-generated -module. (1) There

More information

Simultaneous Optimization of Berth Allocation, Quay Crane Assignment and Quay Crane Scheduling Problems in Container Terminals

Simultaneous Optimization of Berth Allocation, Quay Crane Assignment and Quay Crane Scheduling Problems in Container Terminals Smultaneous Optmzaton of Berth Allocaton, Quay Crane Assgnment and Quay Crane Schedulng Problems n Contaner Termnals Necat Aras, Yavuz Türkoğulları, Z. Caner Taşkın, Kuban Altınel Abstract In ths work,

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

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

REGULAR POSITIVE TERNARY QUADRATIC FORMS. 1. Introduction

REGULAR POSITIVE TERNARY QUADRATIC FORMS. 1. Introduction REGULAR POSITIVE TERNARY QUADRATIC FORMS BYEONG-KWEON OH Abstract. A postve defnte quadratc form f s sad to be regular f t globally represents all ntegers that are represented by the genus of f. In 997

More information

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X Statstcs 1: Probablty Theory II 37 3 EPECTATION OF SEVERAL RANDOM VARIABLES As n Probablty Theory I, the nterest n most stuatons les not on the actual dstrbuton of a random vector, but rather on a number

More information

A combinatorial problem associated with nonograms

A combinatorial problem associated with nonograms A combnatoral problem assocated wth nonograms Jessca Benton Ron Snow Nolan Wallach March 21, 2005 1 Introducton. Ths work was motvated by a queston posed by the second named author to the frst named author

More information

COMPLEX NUMBERS AND QUADRATIC EQUATIONS

COMPLEX NUMBERS AND QUADRATIC EQUATIONS COMPLEX NUMBERS AND QUADRATIC EQUATIONS INTRODUCTION We know that x 0 for all x R e the square of a real number (whether postve, negatve or ero) s non-negatve Hence the equatons x, x, x + 7 0 etc are not

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

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

Homework Notes Week 7

Homework Notes Week 7 Homework Notes Week 7 Math 4 Sprng 4 #4 (a Complete the proof n example 5 that s an nner product (the Frobenus nner product on M n n (F In the example propertes (a and (d have already been verfed so we

More information

arxiv: v2 [cs.ds] 1 Feb 2017

arxiv: v2 [cs.ds] 1 Feb 2017 Polynomal-tme Algorthms for the Subset Feedback Vertex Set Problem on Interval Graphs and Permutaton Graphs Chars Papadopoulos Spyrdon Tzmas arxv:170104634v2 [csds] 1 Feb 2017 Abstract Gven a vertex-weghted

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

Math 261 Exercise sheet 2

Math 261 Exercise sheet 2 Math 261 Exercse sheet 2 http://staff.aub.edu.lb/~nm116/teachng/2017/math261/ndex.html Verson: September 25, 2017 Answers are due for Monday 25 September, 11AM. The use of calculators s allowed. Exercse

More information

Math 217 Fall 2013 Homework 2 Solutions

Math 217 Fall 2013 Homework 2 Solutions Math 17 Fall 013 Homework Solutons Due Thursday Sept. 6, 013 5pm Ths homework conssts of 6 problems of 5 ponts each. The total s 30. You need to fully justfy your answer prove that your functon ndeed has

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

Perron Vectors of an Irreducible Nonnegative Interval Matrix

Perron Vectors of an Irreducible Nonnegative Interval Matrix Perron Vectors of an Irreducble Nonnegatve Interval Matrx Jr Rohn August 4 2005 Abstract As s well known an rreducble nonnegatve matrx possesses a unquely determned Perron vector. As the man result of

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

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

Structure and Drive Paul A. Jensen Copyright July 20, 2003

Structure and Drive Paul A. Jensen Copyright July 20, 2003 Structure and Drve Paul A. Jensen Copyrght July 20, 2003 A system s made up of several operatons wth flow passng between them. The structure of the system descrbes the flow paths from nputs to outputs.

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

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

An Interactive Optimisation Tool for Allocation Problems

An Interactive Optimisation Tool for Allocation Problems An Interactve Optmsaton ool for Allocaton Problems Fredr Bonäs, Joam Westerlund and apo Westerlund Process Desgn Laboratory, Faculty of echnology, Åbo Aadem Unversty, uru 20500, Fnland hs paper presents

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

MATH 241B FUNCTIONAL ANALYSIS - NOTES EXAMPLES OF C ALGEBRAS

MATH 241B FUNCTIONAL ANALYSIS - NOTES EXAMPLES OF C ALGEBRAS MATH 241B FUNCTIONAL ANALYSIS - NOTES EXAMPLES OF C ALGEBRAS These are nformal notes whch cover some of the materal whch s not n the course book. The man purpose s to gve a number of nontrval examples

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

Appendix B. The Finite Difference Scheme

Appendix B. The Finite Difference Scheme 140 APPENDIXES Appendx B. The Fnte Dfference Scheme In ths appendx we present numercal technques whch are used to approxmate solutons of system 3.1 3.3. A comprehensve treatment of theoretcal and mplementaton

More information

Graph Reconstruction by Permutations

Graph Reconstruction by Permutations Graph Reconstructon by Permutatons Perre Ille and Wllam Kocay* Insttut de Mathémathques de Lumny CNRS UMR 6206 163 avenue de Lumny, Case 907 13288 Marselle Cedex 9, France e-mal: lle@ml.unv-mrs.fr Computer

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

Problem Do any of the following determine homomorphisms from GL n (C) to GL n (C)?

Problem Do any of the following determine homomorphisms from GL n (C) to GL n (C)? Homework 8 solutons. Problem 16.1. Whch of the followng defne homomomorphsms from C\{0} to C\{0}? Answer. a) f 1 : z z Yes, f 1 s a homomorphsm. We have that z s the complex conjugate of z. If z 1,z 2

More information

8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS

8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS SECTION 8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS 493 8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS All the vector spaces you have studed thus far n the text are real vector spaces because the scalars

More information

Formulas for the Determinant

Formulas for the Determinant page 224 224 CHAPTER 3 Determnants e t te t e 2t 38 A = e t 2te t e 2t e t te t 2e 2t 39 If 123 A = 345, 456 compute the matrx product A adj(a) What can you conclude about det(a)? For Problems 40 43, use

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

1 Matrix representations of canonical matrices

1 Matrix representations of canonical matrices 1 Matrx representatons of canoncal matrces 2-d rotaton around the orgn: ( ) cos θ sn θ R 0 = sn θ cos θ 3-d rotaton around the x-axs: R x = 1 0 0 0 cos θ sn θ 0 sn θ cos θ 3-d rotaton around the y-axs:

More information

Canonical transformations

Canonical transformations Canoncal transformatons November 23, 2014 Recall that we have defned a symplectc transformaton to be any lnear transformaton M A B leavng the symplectc form nvarant, Ω AB M A CM B DΩ CD Coordnate transformatons,

More information

LECTURE 9 CANONICAL CORRELATION ANALYSIS

LECTURE 9 CANONICAL CORRELATION ANALYSIS LECURE 9 CANONICAL CORRELAION ANALYSIS Introducton he concept of canoncal correlaton arses when we want to quantfy the assocatons between two sets of varables. For example, suppose that the frst set of

More information

Report on Image warping

Report on Image warping Report on Image warpng Xuan Ne, Dec. 20, 2004 Ths document summarzed the algorthms of our mage warpng soluton for further study, and there s a detaled descrpton about the mplementaton of these algorthms.

More information

Discrete Mathematics

Discrete Mathematics Dscrete Mathematcs 30 (00) 48 488 Contents lsts avalable at ScenceDrect Dscrete Mathematcs journal homepage: www.elsever.com/locate/dsc The number of C 3 -free vertces on 3-partte tournaments Ana Paulna

More information

European Journal of Combinatorics

European Journal of Combinatorics European Journal of Combnatorcs 0 (009) 480 489 Contents lsts avalable at ScenceDrect European Journal of Combnatorcs journal homepage: www.elsever.com/locate/ejc Tlngs n Lee metrc P. Horak 1 Unversty

More information

MATH 5707 HOMEWORK 4 SOLUTIONS 2. 2 i 2p i E(X i ) + E(Xi 2 ) ä i=1. i=1

MATH 5707 HOMEWORK 4 SOLUTIONS 2. 2 i 2p i E(X i ) + E(Xi 2 ) ä i=1. i=1 MATH 5707 HOMEWORK 4 SOLUTIONS CİHAN BAHRAN 1. Let v 1,..., v n R m, all lengths v are not larger than 1. Let p 1,..., p n [0, 1] be arbtrary and set w = p 1 v 1 + + p n v n. Then there exst ε 1,..., ε

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

Linear, affine, and convex sets and hulls In the sequel, unless otherwise specified, X will denote a real vector space.

Linear, affine, and convex sets and hulls In the sequel, unless otherwise specified, X will denote a real vector space. Lnear, affne, and convex sets and hulls In the sequel, unless otherwse specfed, X wll denote a real vector space. Lnes and segments. Gven two ponts x, y X, we defne xy = {x + t(y x) : t R} = {(1 t)x +

More information

GELFAND-TSETLIN BASIS FOR THE REPRESENTATIONS OF gl n

GELFAND-TSETLIN BASIS FOR THE REPRESENTATIONS OF gl n GELFAND-TSETLIN BASIS FOR THE REPRESENTATIONS OF gl n KANG LU FINITE DIMENSIONAL REPRESENTATIONS OF gl n Let e j,, j =,, n denote the standard bass of the general lnear Le algebra gl n over the feld of

More information

Support Vector Machines CS434

Support Vector Machines CS434 Support Vector Machnes CS434 Lnear Separators Many lnear separators exst that perfectly classfy all tranng examples Whch of the lnear separators s the best? + + + + + + + + + Intuton of Margn Consder ponts

More information

2 More examples with details

2 More examples with details Physcs 129b Lecture 3 Caltech, 01/15/19 2 More examples wth detals 2.3 The permutaton group n = 4 S 4 contans 4! = 24 elements. One s the dentty e. Sx of them are exchange of two objects (, j) ( to j and

More information

DIFFERENTIAL FORMS BRIAN OSSERMAN

DIFFERENTIAL FORMS BRIAN OSSERMAN DIFFERENTIAL FORMS BRIAN OSSERMAN Dfferentals are an mportant topc n algebrac geometry, allowng the use of some classcal geometrc arguments n the context of varetes over any feld. We wll use them to defne

More information

Salmon: Lectures on partial differential equations. Consider the general linear, second-order PDE in the form. ,x 2

Salmon: Lectures on partial differential equations. Consider the general linear, second-order PDE in the form. ,x 2 Salmon: Lectures on partal dfferental equatons 5. Classfcaton of second-order equatons There are general methods for classfyng hgher-order partal dfferental equatons. One s very general (applyng even to

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

Lecture 3. Ax x i a i. i i

Lecture 3. Ax x i a i. i i 18.409 The Behavor of Algorthms n Practce 2/14/2 Lecturer: Dan Spelman Lecture 3 Scrbe: Arvnd Sankar 1 Largest sngular value In order to bound the condton number, we need an upper bound on the largest

More information