EXTENDED FORMULATIONS FOR PACKING AND PARTITIONING ORBITOPES

Size: px
Start display at page:

Download "EXTENDED FORMULATIONS FOR PACKING AND PARTITIONING ORBITOPES"

Transcription

1 EXTENDED FORMULATIONS FOR PACKING AND PARTITIONING ORBITOPES YURI FAENZA AND VOLKER KAIBEL Abrac. We give compac exended formulaion for he packing and pariioning orbiope (wih repec o he full ymmeric group) decribed and analyzed in [6]. Thee polyope are he convex hull of all 0/1-marice wih lexicographically ored column and a mo, rep. exacly, one 1-enry per row. They are imporan objec for ymmery reducion in cerain ineger program. Uing he exended formulaion, we alo derive a raher imple proof of he fac [6] ha baically hifed-column inequaliie uffice in order o decribe hoe orbiope linearly. 1. Inroducion Exploiaion of ymmerie i crucial for many very difficul ineger programming model. Over he la few year, ignifican progre ha been achieved wih repec o general echnique for dealing wih ymmerie wihin branchand-cu algorihm. Very nice and effecive procedure have been devied, like iomorphim pruning [10, 12, 11, 13] and orbial branching [7, 8]. There ha alo been progre in underanding linear inequaliie o be added o cerain ineger program in order o remove ymmery. Toward hi end, orbiope have been inroduced in [6]. The packing orbiope O p,q and he pariioning orbiope O = p,q are he convex hull of all 0/1-marice of ize p q whoe column are in lexicographically decreaing order having a mo or exacly, repecively, one 1-enry per row. In [6], complee decripion wih linear inequaliie have been derived for hee polyope (ee Thm. 16 and 17 in [6]). Knowledge on orbiope urn ou o be quie ueful in pracical ymmery reducion for cerain ineger programming model. For inance, in a well-known formulaion of he graph pariioning problem (for graph having p node o be pariioned ino q par) he ymmery on he 0/1- variable x ij indicaing wheher node i i pu ino par j of he pariioning he ymmery ariing from permuing he par can be removed by requiring x O = p,q. We refer o [5] and [6] for a more deailed dicuion of he pracical ue of orbiope. Dae: June 2, Thi work ha been uppored by he European Union, FP6, MRTN-CT (ADONET). 1

2 2 FAENZA AND KAIBEL The opic of hi paper are exended formulaion for hee orbiope, i.e., (imple) linear decripion of higher dimenional polyope which can be projeced o O p,q and O = p,q. In fac, uch exended formulaion play imporan role in polyhedral combinaoric and ineger programming in general, becaue raher han olving a linear opimizaion problem over a polyhedron in he original pace, one may olve i over a (hopefully impler decribed) polyhedron of which he fir one i a linear projecion. For inance, he lif-and-projec approach [1] and oher general reformulaion cheme (e.g., he one due o Lováz and Schrijver [9] a well a Adam and Sherali [15]) are baed on exended formulaion. Recen example include work on mixed ineger programming for dual of nework marice [3, 4, 2]. More claical i he general heory on exended formulaion obained from (cerain) dynamic programming algorihm [14]. The reul of hi paper are much in he piri of he laer work. In order o give an overview on he conribuion of hi paper le u fir recall a few fac on orbiope. A no 0/1-marix wih a mo one 1-enry per row and lexicographically (decreaing) ored column ha a one above i main diagonal, we may aume wihou lo of generaliy ha O = p,q O p,q R Ip,q wih I p,q = {(i, j) [p] [q] : i j} (where [n] = {1, 2,..., n}). In fac, O = p,q i he face of O p,q defined by requiring ha all row-um inequaliie x(row i ) 1 for i [p] are aified wih equaliy, where row i = {(i, j) I p,q : j [q]}. The main reul of [6] i a complee decripion of O p,q and O = p,q by mean of linear equaion and inequaliie. Thi yem of conrain (he SCI-yem) coni, nex o nonnegaiviy conrain and row-um inequaliie (or row-um equaion for O = p,q), of he exponenially large cla of hifed column inequaliie (SCI ) ha will be defined a he end of Sec. 2. In [6] i i alo proved ha, up o a few excepion, hee exponenially many SCI define face of he orbiope. The proof given in [6] of he fac ha he SCI-yem compleely decribe hee orbiope i raher lenghy and omewha echnical. Exending over page 18 o 27, i hardly leave (no only) he reader wih a good idea of he reaon for he SCI-yem being ufficien o decribe he orbiope. In conra o hi, he conribuion of he preen work are he following: We provide a quie imple exended formulaion for O p,q (along wih a raher hor proof eablihing hi) and, moreover, we how by ome imple and naural (no echnical) argumen ha he SCI-yem decribe he projecion of he feaible region of ha exended formulaion o he original pace, hu providing a new proof howing ha he SCI-yem decribe O p,q. Thi laer proof i much horer han he original one, and i eem o provide much beer inigh ino he reaon for he SCI-yem o decribe he orbiope. Clearly, a O = p,q i a face of O p,q, he reul for he laer polyope immediaely yield correponding reul for he fir one. However, beide leading o ha impler proof, we believe ha our exended formulaion for O p,q i inereing ielf. I provide a decripion of a quie naural polyope (he orbiope O p,q) by

3 EXTENDED FORMULATIONS FOR PACKING AND PARTITIONING ORBITOPES 3 a yem of conrain in a pace whoe dimenion i roughly wice he original dimenion I p,q wih only linearly (in I p,q ) many nonzero coefficien, while every linear decripion of he orbiope in he original pace require exponenially many inequaliie. Thi may alo urn ou o be compuaionally aracive. The baic idea of our exended formulaion i o aign o each verex of O p,q a direced pah in a cerain acyclic digraph. The addiional variable in our exended formulaion are ued o uiably expre hee pah. The digraph we work wih i e up in Sec. 2, where we alo fix ome noaion and define SCI. In Sec. 3 we hen decribe he exended formulaion for O p,q and O = p,q (Thm. 6 and Cor. 8). The main work i done in Sec. 3.1, where he exended formulaion for O p,q wih addiional variable encoding he pah menioned above i inroduced and proved o define an inegral polyhedron (Thm. 4). From hi i i eay o conclude ha he formulaion indeed define a polyope ha projec down o O p,q (Thm. 6). Boh he exenion of uch reul o he pariioning cae O = p,q (Cor. 8 in Sec. 3.2), and he ranformaion of he yem in order o reduce he number of variable and nonzero coefficien (Thm. 10 in Sec. 3.3) are obained wihou much work. On he way, we alo derive linear (in I p,q ) ime algorihm for opimizing linear objecive funcion over O p,q and O = p,q (Cor. 7 and 9). In Sec. 4 we finally prove ha he projecion of he feaible region defined by he exended formulaion conain he polyope defined by he SCI-yem (Thm. 12), hu providing he new proof of he fac (Thm. 11) ha he laer polyope equal O p,q. We conclude wih a few remark and acknowledgemen in Sec The Seup Le u aume p q 1 hroughou he paper. We define a direced acyclic graph D p,q = (V p,q, A p,q ) wih node e V p,q = I p,q ([p] 0 {0}) {} {} (where [n] 0 = [n] {0} and mean dijoin union). Uing he noaion q(i) = min{i, q}, he e of arc of D p,q i where A p,q = A p,q Ap,q {(, (0, 0))} {((p, j), ) : j [q] 0 }, A p,q = {((i, j), (i + 1, j)) : i [p 1] 0, j [q(i)] 0 } i he e of verical arc ha are denoed by (i, j) = ((i, j), (i + 1, j)), and A p,q = {((i, j), (i + 1, j + 1)) : i [p 1] 0, j [q(i + 1) 1] 0 } i he e of diagonal arc ha are denoed by (i, j) = ((i, j), (i + 1, j + 1)). The crucial propery of D p,q i ha every verex of O p,q induce an --pah in D p,q a indicaed in Fig. 1. Noe ha differen verice may induce he ame pah.

4 4 FAENZA AND KAIBEL Figure 1. The digraph D 8,6 and a verex of O 8,6 along wih i --pah. For a ube W V p,q we ue he following noaion: ou(w ) = {(w, u) A p,q : w W, u W } ou (W ) = {(w, u) A p,q : w W, u W } in(w ) = {(u, w) A p,q : w W, u W } in (W ) = {(u, w) A p,q : w W, u W } in (W ) = {(u, w) Ap,q : w W, u W } For a direced pah Γ in D p,q, we denoe by V(Γ) V p,q he e of node on he pah Γ, S(Γ) V(Γ) he e of node on Γ no enered by Γ via diagonal arc, and T(Γ) V(Γ) he e of node on Γ lef by Γ via diagonal arc (ee Fig. 2). Noe ha S(Γ) alway conain he ar node of Γ, and T(Γ) alway exclude he end node of Γ. Remark 1. For every direced pah Γ in D p,q wih end node (i, j) I i,j, we have in (V(Γ)) = in (S(Γ)) and ou (V(Γ) \ {(i, j)}) = ou (T(Γ)). A ube S I p,q i a hifed column if and only if S = S(Γ) for ome (l, l)-(i 1, j 1)-pah Γ in D p,q wih i [p]\{1}, j [q]\{1}, and l [q]. The aociaed hifed-column inequaliy i x(bar i,j ) x(s), where bar i,j = {(i, l) I p,q : l j}

5 EXTENDED FORMULATIONS FOR PACKING AND PARTITIONING ORBITOPES 5 Figure 2. A pah Γ along wih he e S(Γ) (lef) and T(Γ) (righ). and, a uual, we wrie z(n) = e N z e for ome vecor z R M and a ube N M (ee Fig. 3). 3. Exended formulaion 3.1. The packing cae. Denoe by F p,q R Ap,q he e of all --flow (wihou any capaciy rericion) in D p,q wih flow value one. Clearly, F p,q i an inegral polyope. Since D p,q i acyclic, he verice of F p,q are he incidence vecor of he direced --pah (viewed a ube of arc) in D p,q. For a flow y F p,q and a node (i, j) I p,q, we denoe by y(i, j) = y(in(i, j)) = y(ou(i, j)) he amoun of flow paing node (i, j). For a ube W row i of node in he ame row, y(w ) = w W y(w) i he oal amoun of flow enering W (or, equivalenly, leaving W ). Lemma 2. For a direced (k, l)-(i, j)-pah Γ in D p,q wih (i, j) I p,q he following aemen hold for all y F p,q (ee Fig. 4): (1) If k = l 1 hen (2) If l = 0 hen y(in (S(Γ))) y(ou (T(Γ))) = y(bar i,j ). 1 y(ou (T(Γ))) = y(bar i,j ).

6 6 FAENZA AND KAIBEL Figure 3. Coefficien vecor of wo SCI wih he ame bar bar 8,5. Proof. For boh cae, le W = {(a, b) V p,q \{, } : a a for ome (a, b) V(Γ) bar i,j } be he e of all node (differen from and ) in or above V(Γ) bar i,j. For cae (1), we ar by oberving (1) y(in(w )) = y(ou(w )) (a y F p,q i an --flow). Becaue of (0, 0) W (hu being no adjacen o W ) we have in(w ) = in (V(Γ)), which according o Remark 1 equal in (S(Γ)), yielding (2) y(in(w )) = y(in (S(Γ))). Similarly, we have ou(w ) = ou (V(Γ) \ {(i, j)}) ou(bar i,j ), where Remark 1 give ou (V(Γ) \ {(i, j)}) = ou (T(Γ)). Thu, we obain (3) y(ou(w )) = y(ou (T(Γ))) + y(ou(bar i,j )). Equaion (1), (2), and (3) imply he aemen on cae (1). For cae (2), we exploi he fac ha he --flow y F p,q of value one aifie (4) 1 = y(ou(w {})) y(in(w {})). We have in(w {}) = and ou(w {}) = ou(w ) (due o (0, 0) W in hi cae). From ou(w ) = ou (V(Γ) \ {(i, j)}) ou(bar i,j )

7 EXTENDED FORMULATIONS FOR PACKING AND PARTITIONING ORBITOPES 7 Figure 4. Illuraion of par (1) (lef) and par (2) (righ) of Lemma 2 wih (i, j) = (7, 4). and ou (V(Γ) \ {(i, j)}) = ou (T(Γ)) (ee Remark 1), we hu derive he aemen on cae (2) from (4). For (i, j) I p,q we denoe by col i,j = {(k, j) : j k i} he upper par of he j-h column from (j, j) down o (i, j), including boh node. For he direced pah Γ wih V(Γ) = col i,j, we have S(Γ) = col i,j and T(Γ) =. Thu, par (1) of Lemma 2 implie he following. Remark 3. For all y F p,q and (i, j) I p,q, we have y(in (col i,j )) = y(bar i,j ). The cenral objec of udy of hi paper i he polyope P p,q = {(x, y) R Ip,q R Ap,q : y F p,q and (x, y) aifie (5) and (6) below} wih (5) y (i 1,j 1) x ij for all (i, j) I p,q and (6) x(bar i,j ) y(in (col i,j )) for all (i, j) I p,q (ee Fig. 5). Since F p,q i he e of all --flow of value one, all (x, y) P p,q

8 8 FAENZA AND KAIBEL Figure 5. Coefficien vecor of inequaliie (5) (lef) and (6) (righ). aify 0 y 1, and hu 0 x 1 due o (5) and (6). Furhermore, inequaliie (6) imply he row-um inequaliie x(row i ) 1 for all (x, y) P p,q. Theorem 4. The polyope P p,q i inegral. For he proof of hi heorem, we need he following reul. Lemma 5. Le x 1,..., x n 0 and y 1,..., y n R wih n n x l y l for all 1 j n. l=j l=j For all number α 1,..., α n R and 0 β 1 β 2 β n wih he inequaliy hold. α j β j n α j x j for all 1 j n n β j y j Proof of Lemma 5. We prove he claim by inducion on n. The cae n = 1 i rivial, hu le n 1. Ignoring index 1 and decreaing he remaining α j and β j

9 EXTENDED FORMULATIONS FOR PACKING AND PARTITIONING ORBITOPES 9 by β 1, from he inducion hypohei we obain n n (α j β 1 )x j (β j β 1 )y j. j=2 Due o α 1 β 1 and x 1 0 we have (α 1 β 1 )x 1 0, hu we deduce n n (α j β 1 )x j (β j β 1 )y j. Hence we have n α j x j j=2 n ( n β j y j β 1 x j n y j ), wih a nonpoiive righ-hand ide due o β 1 0 and n x j n y j. Proof of Theorem 4. In order o how ha P p,q i inegral, we how ha for an arbirary objecive funcion vecor c R Ip,q R Ap,q he opimizaion problem (7) max{ c, (x, y) : (x, y) P p,q } ha an opimal oluion wih 0/1-componen. We define wo vecor c (1), c (2) R Ip,q R Ap,q wih he following excepion: c (1) (i 1,j 1) which are zero in all componen = c (i,j) for all (i, j) I p,q c (2) (i 1,j) = max{0, c (i,1),..., c (i,j) } for all (i, j) I p,q We are going o eablih he following wo claim: (1) For each (x, y) P p,q we have c, (x, y) c + c (1) + c (2), (0, y). (2) For each --flow y F p,q {0, 1} Ap,q here i ome x {0, 1} Ip,q wih (x, y) P p,q and c + c (1) + c (2), (0, y) = c, (x, y). Wih hee wo claim, he exience of an inegral opimal oluion o (7) can be eablihed a follow: Le c R Ap,q be he y-par of c + c (1) + c (2). A F p,q i a 0/1-polyope, here i a 0/1-flow y F p,q {0, 1} Ap,q wih c, y = max{ c, y : y F p,q }. Due o c + c (1) + c (2), (0, y) = c, y for all y F p,q, claim (1) implie ha he opimal value of (7) i a mo c, y. On he oher hand, claim (2) enure ha here i ome x {0, 1} Ip,q wih (x, y ) P p,q and c, (x, y ) = c + c (1) + c (2), (0, y ) = c, y. Thu, (x, y ) i an inegral opimal oluion o (7).

10 10 FAENZA AND KAIBEL In order o prove claim (2), le y F p,q {0, 1} Ap,q, i.e., y i he incidence vecor of an --pah in D p,q. For he conrucion of ome x {0, 1} Ip,q a required we ar by iniializing x = 0. For each (i, j) V p,q wih y (i,j) = 1 e x i+1,j+1 = 1. For every (i, j) V p,q wih j 1 and y (i,j) = 1 chooe l [j] wih c i+1,l = max{c (i+1),1,..., c (i+1),j }, and e x i+1,l = 1 if c i+1,l 0. For claim (1) le (x, y) P p,q. Define x R Ip,q x ij = x ij y (i 1,j 1) for all (i, j) I p,q. A (x, y) aifie (5), x 0 hold. Furhermore, we have via (8) c (1), (0, y) = c, (x x, 0). Therefore, i uffice o how (9) c (2), (0, y) c, (x, 0), becaue (8) and (9) yield (10) c + c (1) + c (2), (0, y) = c, (0, y) + c (1), (0, y) + c (2), (0, y) c, (0, y) + c, (x x, 0) + c, (x, 0) = c, (x, y). In order o eablih (9), we prove for every i [p] q(i) c ij x ij q(i 1) c (2) (i 1,j) y (i 1,j), which by ummaion over i [p] yield (9). To ee (10) for ome i [p], oberve ha, for every j [q(i)], we have x (bar i,j ) = x(bar i,j ) y(in (bar i,j )) = x(bar i,j ) y(bar i,j ) + y(in (bar i,j )). Due o Remark 3, and a (x, y) aifie (6), hi implie x (bar i,j ) y(in (bar i,j )). Defining y (i,q(i)) = 0 in cae of q(i 1) < q(i), we hu have q(i) q(i) x il l=j for every j [q(i)]. Seing c (2) (i 1,q(i)) q(i 1) < q(i), we furhermore have l=j y (i 1,l) o he bigge componen of c in cae of 0 c (2) (i 1,1) c (2) (i 1,q(i))

11 EXTENDED FORMULATIONS FOR PACKING AND PARTITIONING ORBITOPES 11 and c ij c (2) ((i 1),j) for all j [q(i)]. Thu we can ue Lemma 5 (wih n = q(i), x j = x ij 0, y j = y (i 1,j), α j = c ij, and β j = c (2) (i 1,j) ) o deduce q(i) q(i) c ij x ij c (2) (i 1,j) y (i 1,j), which yield (10). From Theorem 4 one obain ha P p,q i an exended formulaion for O p,q. Theorem 6. The orbiope O p,q R I p,q i he orhogonal projecion of he polyope P p,q R Ip,q R Ap,q o he pace R Ip,q. Proof. Le x O p,q {0, 1} Ip,q be an arbirary verex of O p,q. The incidence vecor y {0, 1} Ap,q of he unique --pah uing all arc (i 1, j 1) wih (i, j) I p,q and x ij = 1 aifie (x, y) P p,q. Thu, O p,q i conained in he projecion of P p,q. To ee ha vice vera he projecion of P p,q i conained in O p,q, by Theorem 4 i uffice o oberve ha every 0/1-poin (x, y) P p,q i conained in O p,q. Clearly, for uch a poin x ha a mo one one-enry per row (ince he rowum inequaliie are implied by he fac (x, y) P p,q ). Furhermore, if he j-h column of x wa lexicographically larger han he (j 1)- column of x wih i being minimal uch ha x ij = 1 hold, hen one would find ha y(col i,j 1 ) = 0 hold (becaue of (5)), conradicing (6) for (i, j 1). From he proof of Theorem 4, we derive a combinaorial algorihm for he linear opimizaion problem (11) max{ d, x : x O p,q} wih d I p,q. Indeed, wih c = (d, 0) R Ip,q R Ap,q every opimal oluion (x, y ) o (12) max{ c, (x, y) : (x, y) P p,q } yield an opimal oluion x o (11). From he proof of Theorem 4 we know ha we can compue an opimal oluion (x, y ) o (12) by fir compuing he incidence vecor y {0, 1} Ap,q of a longe --pah in he digraph D p,q wih repec o arc lengh given by c (1) + c (2) (which can be done in linear ime ince D p,q i acyclic) and hen eing x {0, 1} Ip,q a decribed in he proof of claim (2) (in he proof of Theorem 4). Corollary 7. Linear opimizaion over O p,q can be olved in ime O(pq).

12 12 FAENZA AND KAIBEL 3.2. The pariioning cae. The previou reul can be eaily exended o he pariioning cae. Since he row-um inequaliie x(row i ) 1 are valid for P p,q, P = p,q = {(x, y) P p,q : x(row i ) = 1 for all i [p]} i a face of P p,q. Clearly, due o Theorem 6 hi face map o he face (ee Sec. 1) {x O p,q : x(row i ) = 1 for all i [p]} = O = p,q of O p,q via he orhogonal projecion ono he x-pace R Ip,q. Corollary 8. P = p,q i an exended formulaion for O = p,q. Suppoe we wan o olve (13) max{ d, x : x O = p,q} for ome d R Ip,q. A all poin x O = p,q aify he row-um equaion x(row i ) = 1 for all i [p], we may add, for each i, an arbirary conan o he objecive funcion coefficien of he variable belonging o row i wihou changing he opimal oluion o (13) (hough, of coure, changing he objecive funcion value of he oluion). Therefore, we may aume ha d ha only poiive componen. Bu hen max{ d, x : x O = p,q} = max{ d, x : x O p,q}, and all opimal oluion o he opimizaion problem over O p,q are poin in O = p,q. Thu we derive he following from Corollary 7. Corollary 9. Linear opimizaion over O = p,q can be olved in ime O(pq) Reducing he number of variable and nonzero elemen. Le u manipulae he defining yem of P p,q in order o decreae he number of variable and nonzero coefficien. Thi may be advanageou for pracical purpoe. I furhermore emphaize he impliciy of he exended formulaion. For he ake of readabiliy, we define bar i,j = and col i,j = whenever j > q(i). Since every y F p,q aifie y (i,j) = y(bar i,j ) y(bar i+1,j+1 ), we deduce from Remark 3 ha for all (x, y) P p,q y (i,j) = y(in (col i,j )) y(in (col i+1,j+1 )) hold for all verical arc (i, j) wih (i, j) I p,q (and i < p) a well a y ((p,j),) = y(in (col p,j )) y(in (col p,j+1 )) for all arc ((p, j), ) wih j [q] 0, where we defined y(in (col p,q+1 )) = 0. Similarly o he derivaion of Remark 3, one furhermore deduce ha every (x, y) P p,q aifie y (i,0) = 1 y(in (col i+1,1 )) for all i [p 1] 0. Finally, every (x, y) P p,q clearly aifie y (,(0,0)) = 1. Therefore, we can eliminae from he yem decribing P p,q all arc variable excep for he one correponding o diagonal arc.

13 EXTENDED FORMULATIONS FOR PACKING AND PARTITIONING ORBITOPES 13 We finally apply he linear ranformaion defined by o R Ip,q R A z ij = x(bar i,j ) and w ij = y(in (col i,j )) for all (i, j) I p,q, p,q, whoe invere i given by x ij = z i,j z i,j+1 for all (i, j) I p,q. (defining z i,q(i)+1 = 0, for all i [p]) and y (i,j) = w i+1,j+1 w i,j+1 for all i [p 1] 0, j [q(i + 1) 1] 0. Few calculaion are needed o check ha he previou ranformaion (bijecively) map P p,q ono he polyope P comp p,q R Ip,q R Ip,q defined by he following very compac e of conrain: (14) (15) (16) (17) (18) (19) w i+1,j+1 w i,j+1 0 for i [p 1] 0, j [q(i + 1) 1] 0 w i,j w i+1,j+1 0 for (i, j) I p,q, i < p w p,1 1 w i,j w i 1,j z ij + z i,j+1 0 for (i, j) I p,q z i,j w i,j 0 for (i, j) I p,q w i,q(i) 0 for i [p] Here, (14) repreen he nonnegaiviy conrain on he diagonal arc. Nonnegaiviy on he verical arc (i, j) wih i [p 1] 0 i refleced by (15) for j [q(i)] and by (16) (ogeher wih he nonnegaiviy of w, which i implied by (19) and (15)) for j = 0. Finally, equaion (5) and (6) ranlae o (17) and (18), repecively. Ignoring he nonnegaiviy conrain, yem (14) (19) ha le han 2pq variable and 4pq conrain, for a oal number of nonzero coefficien ha i maller han 10pq. Noe ha w 1,1 1 i a valid inequaliy for P comp p,q. The face of P comp p,q defined by w 1,1 = 1 i he image of he face P = p,q of P p,q. Thu, adding w 1,1 = 1 o he yem (14) (19) one arrive a anoher exended formulaion of O = p,q. We ummarize he reul of hi ubecion. Theorem 10. The polyope P comp p,q R Ip,q R Ip,q defined by (14) (19) i an exended formulaion of O p,q. The face of P comp p,q defined by w 1,1 = 1 i an exended formulaion of O = p,q. 4. The projecion Le Q p,q R Ip,q be he polyope defined by he nonnegaiviy conrain x 0, he row-um inequaliie x(row i ) 1 for all i [p] and all hifed-column inequaliie. By checking he verice (0/1-vecor) of O p,q i i eay o ee ha O p,q Q p,q hold. Thu, in order o prove Theorem 11. O p,q = Q p,q

14 14 FAENZA AND KAIBEL (which i Prop. 13 in [6]) i uffice (due o Theorem 6) o how he following: Theorem 12. For each x Q p,q here i ome y F p,q wih (x, y) P p,q. Proof. For x Q p,q conider he nework D p,q wih capaciy x ij on he diagonal arc (i 1, j 1) for each (i, j) I p,q and infinie capaciie on all oher arc. In hi nework, we conruc a feaible flow y F p,q of value one wih he propery (20) y (i 1,j 1) > 0 y (i 1,j 1) = x ij for all (i, j) I p,q. Phraed verbally, he flow y ue a verical arc only if he diagonal arc emanaing from i ail i auraed. Such a flow can eaily be conruced in he following way: ar by ending one uni of flow from o along he verical pah in column zero. A each ep, if he flow y conruced o far violae (20) for ome (i, j) I p,q, chooe uch a pair (i, j) wih minimal j, breaking ie by chooing i minimally a well. Wih ϑ = min{y (i 1,j 1), x ij y (i 1,j 1) } (i.e., he minimum of he flow on he verical arc and he reidual capaciy on he diagonal arc aring a (i 1, j 1)) reroue ϑ uni of he flow currenly ravelling on he verical arc (i 1, j 1) along he pah aring wih he diagonal arc (i 1, j 1) and hen uing he verical arc in column j. Noe ha hi affec only arc leaving node (k, l) wih k i and l j. Afer hi rerouing, (20) hold for (i, j). The minimaliy requiremen in he choice of (i, j) enure ha he flow on he wo arc leaving (i, j) i no changed again aferward. Thu, (20) will alway be aified for (i, j) in he fuure. Therefore, he procedure evenually end wih a flow a required. A (x, y) aifie (5) by conrucion, i uffice o how (6) in order o prove (x, y) P p,q. To hi end, le (i, j) I p,q. Due o Remark 3, we only need o prove (21) y(bar i,j ) x(bar i,j ). We conruc a direced (k, l)-(i, j)-pah Γ in he reidual nework wih repec o he flow y (conaining only hoe arc of D p,q ha are no auraed by y) wih k = l or l = 0 in he following way: Saring from he rivial (lengh zero) w- (i, j)-pah wih w = (i, j), in each ep we exend he pah a i curren ar node w by he diagonal arc enering w if hi arc i par of he reidual nework, and by he verical arc enering w oherwie. A he reidual nework conain all verical arc, we clearly can proceed hi way unil he ar node of he curren pah i ome node (k, l) wih k = l or wih l = 0. Since Γ i a pah in he reidual nework and due o (20), we have (22) y(ou (T(Γ))) = 0.

15 EXTENDED FORMULATIONS FOR PACKING AND PARTITIONING ORBITOPES 15 If l = 0, par (2) of Lemma 2 ogeher wih (22) yield y(bar i,j ) = 1, from which (21) follow ince x aifie he row-um inequaliy x(row i ) 1 and he nonnegaiviy conrain. If l 0, hen k = l 1. Thu, according o par (1) of Lemma 2 and due o (22), we have (23) y(in (S(Γ))) = y(bar i,j ). Since we preferred diagonal arc from he reidual nework in our backward conrucion of Γ, we find ha all arc from in (S(Γ)) are auraed by he flow y. Therefore, for he hifed column S = S(Γ), we have (uing (23)) (24) x(s) = y(in (S(Γ)) = y(bar i,j ). Le a A p,q be he arc in Γ enering (i, j), and denoe by Γ he pah ariing from Γ by removing a. If a i diagonal, hen he (l, l)-(i 1, j 1)-pah Γ aifie S(Γ ) = S, and he hifed-column inequaliy x(bar i,j ) x(s) (aified by x) eablihe (21) via (24). If a i verical, by conrucion of Γ, we have y (i 1,j 1) = x ij. Furhermore, he (l, l)-(i 1, j)-pah Γ aifie S(Γ ) = S \{(i, j}). Thu uing he hifed-column inequaliy x(bar i,j+1 ) x(s \ {(i, j)}) one obain from (24) he inequaliy y(bar i,j ) = x(s) = x(s \ {(i, j)}) + x ij x(bar i,j+1 ) + x ij = x(bar i,j ). Thu, (21) i eablihed alo in hi cae, which finally prove Theorem Remark In our view, he exended formulaion for he orbiope O p,q and O = p,q preened in hi paper once more demonrae he power ha lie in he concep of exended formulaion. No only do he exended formulaion provide a very compac way of decribing he orbiope, bu alo do hey allow o derive raher imple proof of he fac ha nonnegaiviy conrain, row-um inequaliie/equaion, and SCI uffice in order o linearly decribe O p,q and O = p,q. To u i eem ha hee proof beer reveal he reaon why SCI are neceary and (baically) ufficien in hee decripion. The conrucion of he flow in he proof of Thm. 12 i quie naural. The re of he proof (i.e., he backward conrucion of he pah Γ) one may alo have done wihou knowing he SCI in advance. Thu, knowing he exended formulaion, one poibly could alo have deeced SCI on he way rying o do hi proof. An inereing pracical queion i wheher he very pare and compac exended formulaion for orbiope lead o performance gain in branch-and-cu algorihm compared o verion ha dynamically add SCI via he linear ime eparaion algorihm (decribed in [6]). Acknowledgemen. We would like o hank Laura Sanià for ueful dicuion and Marc Pfech for valuable commen on an earlier verion of hi paper.

16 16 FAENZA AND KAIBEL Reference 1. Egon Bala, Sebaián Ceria, and Gérard Cornuéjol, A lif-and-projec cuing plane algorihm for mixed 0-1 program, Mah. Programming 58 (1993), no. 3, Ser. A, Michele Confori, Marco Di Summa, Friedrich Eienbrand, and Laurence A. Woley, Nework formulaion of mixed ineger program, CORE Dicuion Paper 2006/ Michele Confori, Marco Di Summa, and Laurence A. Woley, The mixing e wih flow, SIAM J. Dicree Mah. 21 (2007), no. 2, Michele Confori, Ber Gerard, and Giacomo Zambelli, Mixed-ineger verex cover on biparie graph, Proceeding of IPCO XII (Maeo Fichei and David Williamon, ed.), LNCS, vol. 4513, Springer-Verlag, 2007, pp Volker Kaibel, Mahia Peinhard, and Marc E. Pfech, Orbiopal fixing, Proceeding of IPCO XII (Maeo Fichei and David Williamon, ed.), LNCS, vol. 4513, Springer- Verlag, 2007, pp Volker Kaibel and Marc E. Pfech, Packing and pariioning orbiope, Mah. Programming, Ser. A 114 (2008), no. 1, Jeff Linderoh, Jame Orowki, Fabrizio Roi, and Sefano Smriglio, Orbial branching, Proceeding of IPCO XII (Maeo Fichei and David Williamon, ed.), LNCS, vol. 4513, Springer-Verlag, 2007, pp , Conrain orbial branching, Proceeding of IPCO XIII (Andrea Lodi and Giovanni Rinaldi, ed.), LNCS, Springer-Verlag, 2008, o appear. 9. L. Lováz and A. Schrijver, Cone of marice and e-funcion and 0-1 opimizaion, SIAM J. Opim. 1 (1991), no. 2, Franćoi Margo, Pruning by iomorphim in branch-and-cu, Mah. Program. 94 (2002), no. 1, , Exploiing orbi in ymmeric ILP, Mah. Program. 98 (2003), no. 1 3, , Small covering deign by branch-and-cu, Mah. Program. 94 (2003), no. 2 3, , Symmeric ILP: Coloring and mall ineger, Dicree Op. 4 (2007), no. 1, R. Kipp Marin, Ronald L. Rardin, and Brian A. Campbell, Polyhedral characerizaion of dicree dynamic programming, Oper. Re. 38 (1990), no. 1, Hanif D. Sherali and Warren P. Adam, A hierarchy of relaxaion and convex hull characerizaion for mixed-ineger zero-one programming problem, Dicree Appl. Mah. 52 (1994), no. 1, (Faenza) Diparimeno di Ingegneria dell Imprea, Univerià di Roma Tor Vergaa, Rome, Ialy addre, Faenza: faenza@dip.uniroma2.i (Kaibel) Oo-von-Guericke Univeriä Magdeburg, Fakulä für Mahemaik, Univeriäplaz 2, Magdeburg, Germany addre, Kaibel: kaibel@ovgu.de

Problem Set If all directed edges in a network have distinct capacities, then there is a unique maximum flow.

Problem Set If all directed edges in a network have distinct capacities, then there is a unique maximum flow. CSE 202: Deign and Analyi of Algorihm Winer 2013 Problem Se 3 Inrucor: Kamalika Chaudhuri Due on: Tue. Feb 26, 2013 Inrucion For your proof, you may ue any lower bound, algorihm or daa rucure from he ex

More information

Randomized Perfect Bipartite Matching

Randomized Perfect Bipartite Matching Inenive Algorihm Lecure 24 Randomized Perfec Biparie Maching Lecurer: Daniel A. Spielman April 9, 208 24. Inroducion We explain a randomized algorihm by Ahih Goel, Michael Kapralov and Sanjeev Khanna for

More information

Algorithmic Discrete Mathematics 6. Exercise Sheet

Algorithmic Discrete Mathematics 6. Exercise Sheet Algorihmic Dicree Mahemaic. Exercie Shee Deparmen of Mahemaic SS 0 PD Dr. Ulf Lorenz 7. and 8. Juni 0 Dipl.-Mah. David Meffer Verion of June, 0 Groupwork Exercie G (Heap-Sor) Ue Heap-Sor wih a min-heap

More information

18 Extensions of Maximum Flow

18 Extensions of Maximum Flow Who are you?" aid Lunkwill, riing angrily from hi ea. Wha do you wan?" I am Majikhie!" announced he older one. And I demand ha I am Vroomfondel!" houed he younger one. Majikhie urned on Vroomfondel. I

More information

Network Flows: Introduction & Maximum Flow

Network Flows: Introduction & Maximum Flow CSC 373 - lgorihm Deign, nalyi, and Complexiy Summer 2016 Lalla Mouaadid Nework Flow: Inroducion & Maximum Flow We now urn our aenion o anoher powerful algorihmic echnique: Local Search. In a local earch

More information

Introduction to Congestion Games

Introduction to Congestion Games Algorihmic Game Theory, Summer 2017 Inroducion o Congeion Game Lecure 1 (5 page) Inrucor: Thoma Keelheim In hi lecure, we ge o know congeion game, which will be our running example for many concep in game

More information

CSC 364S Notes University of Toronto, Spring, The networks we will consider are directed graphs, where each edge has associated with it

CSC 364S Notes University of Toronto, Spring, The networks we will consider are directed graphs, where each edge has associated with it CSC 36S Noe Univeriy of Torono, Spring, 2003 Flow Algorihm The nework we will conider are direced graph, where each edge ha aociaed wih i a nonnegaive capaciy. The inuiion i ha if edge (u; v) ha capaciy

More information

CS4445/9544 Analysis of Algorithms II Solution for Assignment 1

CS4445/9544 Analysis of Algorithms II Solution for Assignment 1 Conider he following flow nework CS444/944 Analyi of Algorihm II Soluion for Aignmen (0 mark) In he following nework a minimum cu ha capaciy 0 Eiher prove ha hi aemen i rue, or how ha i i fale Uing he

More information

Main Reference: Sections in CLRS.

Main Reference: Sections in CLRS. Maximum Flow Reied 09/09/200 Main Reference: Secion 26.-26. in CLRS. Inroducion Definiion Muli-Source Muli-Sink The Ford-Fulkeron Mehod Reidual Nework Augmening Pah The Max-Flow Min-Cu Theorem The Edmond-Karp

More information

Flow Networks. Ma/CS 6a. Class 14: Flow Exercises

Flow Networks. Ma/CS 6a. Class 14: Flow Exercises 0/0/206 Ma/CS 6a Cla 4: Flow Exercie Flow Nework A flow nework i a digraph G = V, E, ogeher wih a ource verex V, a ink verex V, and a capaciy funcion c: E N. Capaciy Source 7 a b c d e Sink 0/0/206 Flow

More information

1 Motivation and Basic Definitions

1 Motivation and Basic Definitions CSCE : Deign and Analyi of Algorihm Noe on Max Flow Fall 20 (Baed on he preenaion in Chaper 26 of Inroducion o Algorihm, 3rd Ed. by Cormen, Leieron, Rive and Sein.) Moivaion and Baic Definiion Conider

More information

Flow networks. Flow Networks. A flow on a network. Flow networks. The maximum-flow problem. Introduction to Algorithms, Lecture 22 December 5, 2001

Flow networks. Flow Networks. A flow on a network. Flow networks. The maximum-flow problem. Introduction to Algorithms, Lecture 22 December 5, 2001 CS 545 Flow Nework lon Efra Slide courey of Charle Leieron wih mall change by Carola Wenk Flow nework Definiion. flow nework i a direced graph G = (V, E) wih wo diinguihed verice: a ource and a ink. Each

More information

Reminder: Flow Networks

Reminder: Flow Networks 0/0/204 Ma/CS 6a Cla 4: Variou (Flow) Execie Reminder: Flow Nework A flow nework i a digraph G = V, E, ogeher wih a ource verex V, a ink verex V, and a capaciy funcion c: E N. Capaciy Source 7 a b c d

More information

4/12/12. Applications of the Maxflow Problem 7.5 Bipartite Matching. Bipartite Matching. Bipartite Matching. Bipartite matching: the flow network

4/12/12. Applications of the Maxflow Problem 7.5 Bipartite Matching. Bipartite Matching. Bipartite Matching. Bipartite matching: the flow network // Applicaion of he Maxflow Problem. Biparie Maching Biparie Maching Biparie maching. Inpu: undireced, biparie graph = (, E). M E i a maching if each node appear in a mo one edge in M. Max maching: find

More information

Network Flow. Data Structures and Algorithms Andrei Bulatov

Network Flow. Data Structures and Algorithms Andrei Bulatov Nework Flow Daa Srucure and Algorihm Andrei Bulao Algorihm Nework Flow 24-2 Flow Nework Think of a graph a yem of pipe We ue hi yem o pump waer from he ource o ink Eery pipe/edge ha limied capaciy Flow

More information

Admin MAX FLOW APPLICATIONS. Flow graph/networks. Flow constraints 4/30/13. CS lunch today Grading. in-flow = out-flow for every vertex (except s, t)

Admin MAX FLOW APPLICATIONS. Flow graph/networks. Flow constraints 4/30/13. CS lunch today Grading. in-flow = out-flow for every vertex (except s, t) /0/ dmin lunch oday rading MX LOW PPLIION 0, pring avid Kauchak low graph/nework low nework direced, weighed graph (V, ) poiive edge weigh indicaing he capaciy (generally, aume ineger) conain a ingle ource

More information

Graphs III - Network Flow

Graphs III - Network Flow Graph III - Nework Flow Flow nework eup graph G=(V,E) edge capaciy w(u,v) 0 - if edge doe no exi, hen w(u,v)=0 pecial verice: ource verex ; ink verex - no edge ino and no edge ou of Aume every verex v

More information

The Residual Graph. 11 Augmenting Path Algorithms. Augmenting Path Algorithm. Augmenting Path Algorithm

The Residual Graph. 11 Augmenting Path Algorithms. Augmenting Path Algorithm. Augmenting Path Algorithm Augmening Pah Algorihm Greedy-algorihm: ar wih f (e) = everywhere find an - pah wih f (e) < c(e) on every edge augmen flow along he pah repea a long a poible The Reidual Graph From he graph G = (V, E,

More information

The Residual Graph. 12 Augmenting Path Algorithms. Augmenting Path Algorithm. Augmenting Path Algorithm

The Residual Graph. 12 Augmenting Path Algorithms. Augmenting Path Algorithm. Augmenting Path Algorithm Augmening Pah Algorihm Greedy-algorihm: ar wih f (e) = everywhere find an - pah wih f (e) < c(e) on every edge augmen flow along he pah repea a long a poible The Reidual Graph From he graph G = (V, E,

More information

6.8 Laplace Transform: General Formulas

6.8 Laplace Transform: General Formulas 48 HAP. 6 Laplace Tranform 6.8 Laplace Tranform: General Formula Formula Name, ommen Sec. F() l{ f ()} e f () d f () l {F()} Definiion of Tranform Invere Tranform 6. l{af () bg()} al{f ()} bl{g()} Lineariy

More information

To become more mathematically correct, Circuit equations are Algebraic Differential equations. from KVL, KCL from the constitutive relationship

To become more mathematically correct, Circuit equations are Algebraic Differential equations. from KVL, KCL from the constitutive relationship Laplace Tranform (Lin & DeCarlo: Ch 3) ENSC30 Elecric Circui II The Laplace ranform i an inegral ranformaion. I ranform: f ( ) F( ) ime variable complex variable From Euler > Lagrange > Laplace. Hence,

More information

Algorithm Design and Analysis

Algorithm Design and Analysis Algorihm Deign and Analyi LECTURE 0 Nework Flow Applicaion Biparie maching Edge-dijoin pah Adam Smih 0//0 A. Smih; baed on lide by E. Demaine, C. Leieron, S. Rakhodnikova, K. Wayne La ime: Ford-Fulkeron

More information

, the. L and the L. x x. max. i n. It is easy to show that these two norms satisfy the following relation: x x n x = (17.3) max

, the. L and the L. x x. max. i n. It is easy to show that these two norms satisfy the following relation: x x n x = (17.3) max ecure 8 7. Sabiliy Analyi For an n dimenional vecor R n, he and he vecor norm are defined a: = T = i n i (7.) I i eay o how ha hee wo norm aify he following relaion: n (7.) If a vecor i ime-dependen, hen

More information

CS 473G Lecture 15: Max-Flow Algorithms and Applications Fall 2005

CS 473G Lecture 15: Max-Flow Algorithms and Applications Fall 2005 CS 473G Lecure 1: Max-Flow Algorihm and Applicaion Fall 200 1 Max-Flow Algorihm and Applicaion (November 1) 1.1 Recap Fix a direced graph G = (V, E) ha doe no conain boh an edge u v and i reveral v u,

More information

Introduction to SLE Lecture Notes

Introduction to SLE Lecture Notes Inroducion o SLE Lecure Noe May 13, 16 - The goal of hi ecion i o find a ufficien condiion of λ for he hull K o be generaed by a imple cure. I urn ou if λ 1 < 4 hen K i generaed by a imple curve. We will

More information

Matching. Slides designed by Kevin Wayne.

Matching. Slides designed by Kevin Wayne. Maching Maching. Inpu: undireced graph G = (V, E). M E i a maching if each node appear in a mo edge in M. Max maching: find a max cardinaliy maching. Slide deigned by Kevin Wayne. Biparie Maching Biparie

More information

FLAT CYCLOTOMIC POLYNOMIALS OF ORDER FOUR AND HIGHER

FLAT CYCLOTOMIC POLYNOMIALS OF ORDER FOUR AND HIGHER #A30 INTEGERS 10 (010), 357-363 FLAT CYCLOTOMIC POLYNOMIALS OF ORDER FOUR AND HIGHER Nahan Kaplan Deparmen of Mahemaic, Harvard Univeriy, Cambridge, MA nkaplan@mah.harvard.edu Received: 7/15/09, Revied:

More information

16 Max-Flow Algorithms and Applications

16 Max-Flow Algorithms and Applications Algorihm A proce canno be underood by opping i. Underanding mu move wih he flow of he proce, mu join i and flow wih i. The Fir Law of Mena, in Frank Herber Dune (196) There a difference beween knowing

More information

Algorithms and Data Structures 2011/12 Week 9 Solutions (Tues 15th - Fri 18th Nov)

Algorithms and Data Structures 2011/12 Week 9 Solutions (Tues 15th - Fri 18th Nov) Algorihm and Daa Srucure 2011/ Week Soluion (Tue 15h - Fri 18h No) 1. Queion: e are gien 11/16 / 15/20 8/13 0/ 1/ / 11/1 / / To queion: (a) Find a pair of ube X, Y V uch ha f(x, Y) = f(v X, Y). (b) Find

More information

Average Case Lower Bounds for Monotone Switching Networks

Average Case Lower Bounds for Monotone Switching Networks Average Cae Lower Bound for Monoone Swiching Nework Yuval Filmu, Toniann Piai, Rober Robere, Sephen Cook Deparmen of Compuer Science Univeriy of Torono Monoone Compuaion (Refreher) Monoone circui were

More information

They were originally developed for network problem [Dantzig, Ford, Fulkerson 1956]

They were originally developed for network problem [Dantzig, Ford, Fulkerson 1956] 6. Inroducion... 6. The primal-dual algorihmn... 6 6. Remark on he primal-dual algorihmn... 7 6. A primal-dual algorihmn for he hore pah problem... 8... 9 6.6 A primal-dual algorihmn for he weighed maching

More information

Laplace Transform. Inverse Laplace Transform. e st f(t)dt. (2)

Laplace Transform. Inverse Laplace Transform. e st f(t)dt. (2) Laplace Tranform Maoud Malek The Laplace ranform i an inegral ranform named in honor of mahemaician and aronomer Pierre-Simon Laplace, who ued he ranform in hi work on probabiliy heory. I i a powerful

More information

The multisubset sum problem for finite abelian groups

The multisubset sum problem for finite abelian groups Alo available a hp://amc-journal.eu ISSN 1855-3966 (prined edn.), ISSN 1855-3974 (elecronic edn.) ARS MATHEMATICA CONTEMPORANEA 8 (2015) 417 423 The muliube um problem for finie abelian group Amela Muraović-Ribić

More information

CSE 521: Design & Analysis of Algorithms I

CSE 521: Design & Analysis of Algorithms I CSE 52: Deign & Analyi of Algorihm I Nework Flow Paul Beame Biparie Maching Given: A biparie graph G=(V,E) M E i a maching in G iff no wo edge in M hare a verex Goal: Find a maching M in G of maximum poible

More information

EECE 301 Signals & Systems Prof. Mark Fowler

EECE 301 Signals & Systems Prof. Mark Fowler EECE 31 Signal & Syem Prof. Mark Fowler Noe Se #27 C-T Syem: Laplace Tranform Power Tool for yem analyi Reading Aignmen: Secion 6.1 6.3 of Kamen and Heck 1/18 Coure Flow Diagram The arrow here how concepual

More information

Greedy. I Divide and Conquer. I Dynamic Programming. I Network Flows. Network Flow. I Previous topics: design techniques

Greedy. I Divide and Conquer. I Dynamic Programming. I Network Flows. Network Flow. I Previous topics: design techniques Algorihm Deign Technique CS : Nework Flow Dan Sheldon April, reedy Divide and Conquer Dynamic Programming Nework Flow Comparion Nework Flow Previou opic: deign echnique reedy Divide and Conquer Dynamic

More information

Today s topics. CSE 421 Algorithms. Problem Reduction Examples. Problem Reduction. Undirected Network Flow. Bipartite Matching. Problem Reductions

Today s topics. CSE 421 Algorithms. Problem Reduction Examples. Problem Reduction. Undirected Network Flow. Bipartite Matching. Problem Reductions Today opic CSE Algorihm Richard Anderon Lecure Nework Flow Applicaion Prolem Reducion Undireced Flow o Flow Biparie Maching Dijoin Pah Prolem Circulaion Loweround conrain on flow Survey deign Prolem Reducion

More information

18.03SC Unit 3 Practice Exam and Solutions

18.03SC Unit 3 Practice Exam and Solutions Sudy Guide on Sep, Dela, Convoluion, Laplace You can hink of he ep funcion u() a any nice mooh funcion which i for < a and for > a, where a i a poiive number which i much maller han any ime cale we care

More information

Lecture 26. Lucas and Stokey: Optimal Monetary and Fiscal Policy in an Economy without Capital (JME 1983) t t

Lecture 26. Lucas and Stokey: Optimal Monetary and Fiscal Policy in an Economy without Capital (JME 1983) t t Lecure 6. Luca and Sokey: Opimal Moneary and Fical Policy in an Economy wihou Capial (JME 983. A argued in Kydland and Preco (JPE 977, Opimal governmen policy i likely o be ime inconien. Fiher (JEDC 98

More information

NECESSARY AND SUFFICIENT CONDITIONS FOR LATENT SEPARABILITY

NECESSARY AND SUFFICIENT CONDITIONS FOR LATENT SEPARABILITY NECESSARY AND SUFFICIENT CONDITIONS FOR LATENT SEPARABILITY Ian Crawford THE INSTITUTE FOR FISCAL STUDIES DEPARTMENT OF ECONOMICS, UCL cemmap working paper CWP02/04 Neceary and Sufficien Condiion for Laen

More information

Mon Apr 2: Laplace transform and initial value problems like we studied in Chapter 5

Mon Apr 2: Laplace transform and initial value problems like we studied in Chapter 5 Mah 225-4 Week 2 April 2-6 coninue.-.3; alo cover par of.4-.5, EP 7.6 Mon Apr 2:.-.3 Laplace ranform and iniial value problem like we udied in Chaper 5 Announcemen: Warm-up Exercie: Recall, The Laplace

More information

Network Flows UPCOPENCOURSEWARE number 34414

Network Flows UPCOPENCOURSEWARE number 34414 Nework Flow UPCOPENCOURSEWARE number Topic : F.-Javier Heredia Thi work i licened under he Creaive Common Aribuion- NonCommercial-NoDeriv. Unpored Licene. To view a copy of hi licene, vii hp://creaivecommon.org/licene/by-nc-nd/./

More information

Discussion Session 2 Constant Acceleration/Relative Motion Week 03

Discussion Session 2 Constant Acceleration/Relative Motion Week 03 PHYS 100 Dicuion Seion Conan Acceleraion/Relaive Moion Week 03 The Plan Today you will work wih your group explore he idea of reference frame (i.e. relaive moion) and moion wih conan acceleraion. You ll

More information

We just finished the Erdős-Stone Theorem, and ex(n, F ) (1 1/(χ(F ) 1)) ( n

We just finished the Erdős-Stone Theorem, and ex(n, F ) (1 1/(χ(F ) 1)) ( n Lecure 3 - Kövari-Sós-Turán Theorem Jacques Versraëe jacques@ucsd.edu We jus finished he Erdős-Sone Theorem, and ex(n, F ) ( /(χ(f ) )) ( n 2). So we have asympoics when χ(f ) 3 bu no when χ(f ) = 2 i.e.

More information

Geometric Path Problems with Violations

Geometric Path Problems with Violations Click here o view linked Reference 1 1 1 0 1 0 1 0 1 0 1 Geomeric Pah Problem wih Violaion Anil Mahehwari 1, Subha C. Nandy, Drimi Paanayak, Saanka Roy and Michiel Smid 1 1 School of Compuer Science, Carleon

More information

EE Control Systems LECTURE 2

EE Control Systems LECTURE 2 Copyrigh F.L. Lewi 999 All righ reerved EE 434 - Conrol Syem LECTURE REVIEW OF LAPLACE TRANSFORM LAPLACE TRANSFORM The Laplace ranform i very ueful in analyi and deign for yem ha are linear and ime-invarian

More information

Optimal State-Feedback Control Under Sparsity and Delay Constraints

Optimal State-Feedback Control Under Sparsity and Delay Constraints Opimal Sae-Feedback Conrol Under Spariy and Delay Conrain Andrew Lamperki Lauren Leard 2 3 rd IFAC Workhop on Diribued Eimaion and Conrol in Neworked Syem NecSy pp. 24 29, 22 Abrac Thi paper preen he oluion

More information

Selfish Routing. Tim Roughgarden Cornell University. Includes joint work with Éva Tardos

Selfish Routing. Tim Roughgarden Cornell University. Includes joint work with Éva Tardos Selfih Rouing Tim Roughgarden Cornell Univeriy Include join work wih Éva Tardo 1 Which roue would you chooe? Example: one uni of raffic (e.g., car) wan o go from o delay = 1 hour (no congeion effec) long

More information

Exponential Sawtooth

Exponential Sawtooth ECPE 36 HOMEWORK 3: PROPERTIES OF THE FOURIER TRANSFORM SOLUTION. Exponenial Sawooh: The eaie way o do hi problem i o look a he Fourier ranform of a ingle exponenial funcion, () = exp( )u(). From he able

More information

Notes for Lecture 17-18

Notes for Lecture 17-18 U.C. Berkeley CS278: Compuaional Complexiy Handou N7-8 Professor Luca Trevisan April 3-8, 2008 Noes for Lecure 7-8 In hese wo lecures we prove he firs half of he PCP Theorem, he Amplificaion Lemma, up

More information

Mathematische Annalen

Mathematische Annalen Mah. Ann. 39, 33 339 (997) Mahemaiche Annalen c Springer-Verlag 997 Inegraion by par in loop pace Elon P. Hu Deparmen of Mahemaic, Norhweern Univeriy, Evanon, IL 628, USA (e-mail: elon@@mah.nwu.edu) Received:

More information

Soviet Rail Network, 1955

Soviet Rail Network, 1955 Sovie Rail Nework, 1 Reference: On he hiory of he ranporaion and maximum flow problem. Alexander Schrijver in Mah Programming, 1: 3,. Maximum Flow and Minimum Cu Max flow and min cu. Two very rich algorihmic

More information

arxiv: v1 [cs.cg] 21 Mar 2013

arxiv: v1 [cs.cg] 21 Mar 2013 On he rech facor of he Thea-4 graph Lui Barba Proenji Boe Jean-Lou De Carufel André van Renen Sander Verdoncho arxiv:1303.5473v1 [c.cg] 21 Mar 2013 Abrac In hi paper we how ha he θ-graph wih 4 cone ha

More information

2. VECTORS. R Vectors are denoted by bold-face characters such as R, V, etc. The magnitude of a vector, such as R, is denoted as R, R, V

2. VECTORS. R Vectors are denoted by bold-face characters such as R, V, etc. The magnitude of a vector, such as R, is denoted as R, R, V ME 352 VETS 2. VETS Vecor algebra form he mahemaical foundaion for kinemaic and dnamic. Geomer of moion i a he hear of boh he kinemaic and dnamic of mechanical em. Vecor anali i he imehonored ool for decribing

More information

Please Complete Course Survey. CMPSCI 311: Introduction to Algorithms. Approximation Algorithms. Coping With NP-Completeness. Greedy Vertex Cover

Please Complete Course Survey. CMPSCI 311: Introduction to Algorithms. Approximation Algorithms. Coping With NP-Completeness. Greedy Vertex Cover Pleae Complee Coure Survey CMPSCI : Inroducion o Algorihm Dealing wih NP-Compleene Dan Sheldon hp: //owl.oi.uma.edu/parner/coureevalsurvey/uma/ Univeriy of Maachue Slide Adaped from Kevin Wayne La Compiled:

More information

Max Flow, Min Cut COS 521. Kevin Wayne Fall Soviet Rail Network, Cuts. Minimum Cut Problem. Flow network.

Max Flow, Min Cut COS 521. Kevin Wayne Fall Soviet Rail Network, Cuts. Minimum Cut Problem. Flow network. Sovie Rail Nework, Max Flow, Min u OS Kevin Wayne Fall Reference: On he hiory of he ranporaion and maximum flow problem. lexander Schrijver in Mah Programming, :,. Minimum u Problem u Flow nework.! Digraph

More information

FIXED POINTS AND STABILITY IN NEUTRAL DIFFERENTIAL EQUATIONS WITH VARIABLE DELAYS

FIXED POINTS AND STABILITY IN NEUTRAL DIFFERENTIAL EQUATIONS WITH VARIABLE DELAYS PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 136, Number 3, March 28, Page 99 918 S 2-9939(7)989-2 Aricle elecronically publihed on November 3, 27 FIXED POINTS AND STABILITY IN NEUTRAL DIFFERENTIAL

More information

CONTROL SYSTEMS. Chapter 10 : State Space Response

CONTROL SYSTEMS. Chapter 10 : State Space Response CONTROL SYSTEMS Chaper : Sae Space Repone GATE Objecive & Numerical Type Soluion Queion 5 [GATE EE 99 IIT-Bombay : Mark] Conider a econd order yem whoe ae pace repreenaion i of he form A Bu. If () (),

More information

On the Exponential Operator Functions on Time Scales

On the Exponential Operator Functions on Time Scales dvance in Dynamical Syem pplicaion ISSN 973-5321, Volume 7, Number 1, pp. 57 8 (212) hp://campu.m.edu/ada On he Exponenial Operaor Funcion on Time Scale laa E. Hamza Cairo Univeriy Deparmen of Mahemaic

More information

u(t) Figure 1. Open loop control system

u(t) Figure 1. Open loop control system Open loop conrol v cloed loop feedbac conrol The nex wo figure preen he rucure of open loop and feedbac conrol yem Figure how an open loop conrol yem whoe funcion i o caue he oupu y o follow he reference

More information

Linear Algebra Primer

Linear Algebra Primer Linear Algebra rimer And a video dicuion of linear algebra from EE263 i here (lecure 3 and 4): hp://ee.anford.edu/coure/ee263 lide from Sanford CS3 Ouline Vecor and marice Baic Mari Operaion Deerminan,

More information

MAXIMUM FLOW. introduction Ford-Fulkerson algorithm maxflow-mincut theorem

MAXIMUM FLOW. introduction Ford-Fulkerson algorithm maxflow-mincut theorem MAXIMUM FLOW inroducion Ford-Fulkeron algorihm maxflow-mincu heorem Mincu problem Inpu. An edge-weighed digraph, ource verex, and arge verex. each edge ha a poiive capaciy capaciy 9 10 4 15 15 10 5 8 10

More information

CHAPTER 7: SECOND-ORDER CIRCUITS

CHAPTER 7: SECOND-ORDER CIRCUITS EEE5: CI RCUI T THEORY CHAPTER 7: SECOND-ORDER CIRCUITS 7. Inroducion Thi chaper conider circui wih wo orage elemen. Known a econd-order circui becaue heir repone are decribed by differenial equaion ha

More information

CSE 421 Introduction to Algorithms Winter The Network Flow Problem

CSE 421 Introduction to Algorithms Winter The Network Flow Problem CSE 42 Inroducion o Algorihm Winer 202 The Nework Flow Problem 2 The Nework Flow Problem 5 a 4 3 x 3 7 6 b 4 y 4 7 6 c 5 z How much uff can flow from o? 3 Sovie Rail Nework, 955 Reference: On he hiory

More information

Flow networks, flow, maximum flow. Some definitions. Edmonton. Saskatoon Winnipeg. Vancouver Regina. Calgary. 12/12 a.

Flow networks, flow, maximum flow. Some definitions. Edmonton. Saskatoon Winnipeg. Vancouver Regina. Calgary. 12/12 a. Flow nework, flow, maximum flow Can inerpre direced graph a flow nework. Maerial coure hrough ome yem from ome ource o ome ink. Source produce maerial a ome eady rae, ink conume a ame rae. Example: waer

More information

Chapter 7: Inverse-Response Systems

Chapter 7: Inverse-Response Systems Chaper 7: Invere-Repone Syem Normal Syem Invere-Repone Syem Baic Sar ou in he wrong direcion End up in he original eady-ae gain value Two or more yem wih differen magniude and cale in parallel Main yem

More information

Soviet Rail Network, 1955

Soviet Rail Network, 1955 7.1 Nework Flow Sovie Rail Nework, 19 Reerence: On he hiory o he ranporaion and maximum low problem. lexander Schrijver in Mah Programming, 91: 3, 00. (See Exernal Link ) Maximum Flow and Minimum Cu Max

More information

Maximum Flow and Minimum Cut

Maximum Flow and Minimum Cut // Sovie Rail Nework, Maximum Flow and Minimum Cu Max flow and min cu. Two very rich algorihmic problem. Cornerone problem in combinaorial opimizaion. Beauiful mahemaical dualiy. Nework Flow Flow nework.

More information

INDEPENDENT SETS IN GRAPHS WITH GIVEN MINIMUM DEGREE

INDEPENDENT SETS IN GRAPHS WITH GIVEN MINIMUM DEGREE INDEPENDENT SETS IN GRAPHS WITH GIVEN MINIMUM DEGREE JAMES ALEXANDER, JONATHAN CUTLER, AND TIM MINK Absrac The enumeraion of independen ses in graphs wih various resricions has been a opic of much ineres

More information

Let. x y. denote a bivariate time series with zero mean.

Let. x y. denote a bivariate time series with zero mean. Linear Filer Le x y : T denoe a bivariae ime erie wih zero mean. Suppoe ha he ime erie {y : T} i conruced a follow: y a x The ime erie {y : T} i aid o be conruced from {x : T} by mean of a Linear Filer.

More information

Timed Circuits. Asynchronous Circuit Design. Timing Relationships. A Simple Example. Timed States. Timing Sequences. ({r 6 },t6 = 1.

Timed Circuits. Asynchronous Circuit Design. Timing Relationships. A Simple Example. Timed States. Timing Sequences. ({r 6 },t6 = 1. Timed Circuis Asynchronous Circui Design Chris J. Myers Lecure 7: Timed Circuis Chaper 7 Previous mehods only use limied knowledge of delays. Very robus sysems, bu exremely conservaive. Large funcional

More information

Chapter 6. Laplace Transforms

Chapter 6. Laplace Transforms Chaper 6. Laplace Tranform Kreyzig by YHLee;45; 6- An ODE i reduced o an algebraic problem by operaional calculu. The equaion i olved by algebraic manipulaion. The reul i ranformed back for he oluion of

More information

Physics 240: Worksheet 16 Name

Physics 240: Worksheet 16 Name Phyic 4: Workhee 16 Nae Non-unifor circular oion Each of hee proble involve non-unifor circular oion wih a conan α. (1) Obain each of he equaion of oion for non-unifor circular oion under a conan acceleraion,

More information

Price of Stability and Introduction to Mechanism Design

Price of Stability and Introduction to Mechanism Design Algorihmic Game Theory Summer 2017, Week 5 ETH Zürich Price of Sabiliy and Inroducion o Mechanim Deign Paolo Penna Thi i he lecure where we ar deigning yem which involve elfih player. Roughly peaking,

More information

23 Maximum Flows and Minimum Cuts

23 Maximum Flows and Minimum Cuts A proce canno be underood by opping i. Underanding mu move wih he flow of he proce, mu join i and flow wih i. The Fir Law of Mena, in Frank Herber Dune (196) Conrary o expecaion, flow uually happen no

More information

! Abstraction for material flowing through the edges. ! G = (V, E) = directed graph, no parallel edges.

! Abstraction for material flowing through the edges. ! G = (V, E) = directed graph, no parallel edges. Sovie Rail Nework, haper Nework Flow Slide by Kevin Wayne. opyrigh Pearon-ddion Weley. ll righ reerved. Reference: On he hiory of he ranporaion and maximum flow problem. lexander Schrijver in Mah Programming,

More information

5.2 GRAPHICAL VELOCITY ANALYSIS Polygon Method

5.2 GRAPHICAL VELOCITY ANALYSIS Polygon Method ME 352 GRHICL VELCITY NLYSIS 52 GRHICL VELCITY NLYSIS olygon Mehod Velociy analyi form he hear of kinemaic and dynamic of mechanical yem Velociy analyi i uually performed following a poiion analyi; ie,

More information

t )? How would you have tried to solve this problem in Chapter 3?

t )? How would you have tried to solve this problem in Chapter 3? Exercie 9) Ue Laplace ranform o wrie down he oluion o 2 x x = F in x = x x = v. wha phenomena do oluion o hi DE illurae (even hough we're forcing wih in co )? How would you have ried o olve hi problem

More information

Maximum Flow in Planar Graphs

Maximum Flow in Planar Graphs Maximum Flow in Planar Graph Planar Graph and i Dual Dualiy i defined for direced planar graph a well Minimum - cu in undireced planar graph An - cu (undireced graph) An - cu The dual o he cu Cu/Cycle

More information

T L. t=1. Proof of Lemma 1. Using the marginal cost accounting in Equation(4) and standard arguments. t )+Π RB. t )+K 1(Q RB

T L. t=1. Proof of Lemma 1. Using the marginal cost accounting in Equation(4) and standard arguments. t )+Π RB. t )+K 1(Q RB Elecronic Companion EC.1. Proofs of Technical Lemmas and Theorems LEMMA 1. Le C(RB) be he oal cos incurred by he RB policy. Then we have, T L E[C(RB)] 3 E[Z RB ]. (EC.1) Proof of Lemma 1. Using he marginal

More information

Algorithm Design and Analysis

Algorithm Design and Analysis Algorihm Deign and Analyi LECTURES 17 Nework Flow Dualiy of Max Flow and Min Cu Algorihm: Ford-Fulkeron Capaciy Scaling Sofya Rakhodnikova S. Rakhodnikova; baed on lide by E. Demaine, C. Leieron, A. Smih,

More information

23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes

23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes Represening Periodic Funcions by Fourier Series 3. Inroducion In his Secion we show how a periodic funcion can be expressed as a series of sines and cosines. We begin by obaining some sandard inegrals

More information

Syntactic Complexity of Suffix-Free Languages. Marek Szykuła

Syntactic Complexity of Suffix-Free Languages. Marek Szykuła Inroducion Upper Bound on Synacic Complexiy of Suffix-Free Language Univeriy of Wrocław, Poland Join work wih Januz Brzozowki Univeriy of Waerloo, Canada DCFS, 25.06.2015 Abrac Inroducion Sae and ynacic

More information

Additional Methods for Solving DSGE Models

Additional Methods for Solving DSGE Models Addiional Mehod for Solving DSGE Model Karel Meren, Cornell Univeriy Reference King, R. G., Ploer, C. I. & Rebelo, S. T. (1988), Producion, growh and buine cycle: I. he baic neoclaical model, Journal of

More information

Notes on cointegration of real interest rates and real exchange rates. ρ (2)

Notes on cointegration of real interest rates and real exchange rates. ρ (2) Noe on coinegraion of real inere rae and real exchange rae Charle ngel, Univeriy of Wiconin Le me ar wih he obervaion ha while he lieraure (mo prominenly Meee and Rogoff (988) and dion and Paul (993))

More information

Selfish Routing and the Price of Anarchy. Tim Roughgarden Cornell University

Selfish Routing and the Price of Anarchy. Tim Roughgarden Cornell University Selfih Rouing and he Price of Anarchy Tim Roughgarden Cornell Univeriy 1 Algorihm for Self-Inereed Agen Our focu: problem in which muliple agen (people, compuer, ec.) inerac Moivaion: he Inerne decenralized

More information

CHAPTER. Forced Equations and Systems { } ( ) ( ) 8.1 The Laplace Transform and Its Inverse. Transforms from the Definition.

CHAPTER. Forced Equations and Systems { } ( ) ( ) 8.1 The Laplace Transform and Its Inverse. Transforms from the Definition. CHAPTER 8 Forced Equaion and Syem 8 The aplace Tranform and I Invere Tranform from he Definiion 5 5 = b b {} 5 = 5e d = lim5 e = ( ) b {} = e d = lim e + e d b = (inegraion by par) = = = = b b ( ) ( )

More information

Bernoulli numbers. Francesco Chiatti, Matteo Pintonello. December 5, 2016

Bernoulli numbers. Francesco Chiatti, Matteo Pintonello. December 5, 2016 UNIVERSITÁ DEGLI STUDI DI PADOVA, DIPARTIMENTO DI MATEMATICA TULLIO LEVI-CIVITA Bernoulli numbers Francesco Chiai, Maeo Pinonello December 5, 206 During las lessons we have proved he Las Ferma Theorem

More information

Research Article An Upper Bound on the Critical Value β Involved in the Blasius Problem

Research Article An Upper Bound on the Critical Value β Involved in the Blasius Problem Hindawi Publihing Corporaion Journal of Inequaliie and Applicaion Volume 2010, Aricle ID 960365, 6 page doi:10.1155/2010/960365 Reearch Aricle An Upper Bound on he Criical Value Involved in he Blaiu Problem

More information

Generalized Orlicz Spaces and Wasserstein Distances for Convex-Concave Scale Functions

Generalized Orlicz Spaces and Wasserstein Distances for Convex-Concave Scale Functions Generalized Orlicz Space and Waerein Diance for Convex-Concave Scale Funcion Karl-Theodor Surm Abrac Given a ricly increaing, coninuou funcion ϑ : R + R +, baed on he co funcional ϑ (d(x, y dq(x, y, we

More information

Fractional Ornstein-Uhlenbeck Bridge

Fractional Ornstein-Uhlenbeck Bridge WDS'1 Proceeding of Conribued Paper, Par I, 21 26, 21. ISBN 978-8-7378-139-2 MATFYZPRESS Fracional Ornein-Uhlenbeck Bridge J. Janák Charle Univeriy, Faculy of Mahemaic and Phyic, Prague, Czech Republic.

More information

arxiv: v2 [cs.cg] 30 Apr 2015

arxiv: v2 [cs.cg] 30 Apr 2015 Diffue Reflecion Diameer in Simple Polygon arxiv:1302.2271v2 [c.cg] 30 Apr 2015 Gill Bareque Deparmen of Compuer Science Technion Haifa, Irael bareque@c.echnion.ac.il Eli Fox-Epein Deparmen of Compuer

More information

ARTIFICIAL INTELLIGENCE. Markov decision processes

ARTIFICIAL INTELLIGENCE. Markov decision processes INFOB2KI 2017-2018 Urech Univeriy The Neherland ARTIFICIAL INTELLIGENCE Markov deciion procee Lecurer: Silja Renooij Thee lide are par of he INFOB2KI Coure Noe available from www.c.uu.nl/doc/vakken/b2ki/chema.hml

More information

The Primal Simplex Approach to the QoS Routing Problem

The Primal Simplex Approach to the QoS Routing Problem The Primal Simplex Approach o he QoS Rouing Problem Ying Xiao, Krihnaiyan Thulairaman Univeriy of Oklahoma, Norman, OK {ying_xiao, hulai}@ou.edu Guoliang Xue Arizona Sae Univeriy, Tempe, AZ xue@au.edu

More information

Inventory Analysis and Management. Multi-Period Stochastic Models: Optimality of (s, S) Policy for K-Convex Objective Functions

Inventory Analysis and Management. Multi-Period Stochastic Models: Optimality of (s, S) Policy for K-Convex Objective Functions Muli-Period Sochasic Models: Opimali of (s, S) Polic for -Convex Objecive Funcions Consider a seing similar o he N-sage newsvendor problem excep ha now here is a fixed re-ordering cos (> 0) for each (re-)order.

More information

11!Hí MATHEMATICS : ERDŐS AND ULAM PROC. N. A. S. of decomposiion, properly speaking) conradics he possibiliy of defining a counably addiive real-valu

11!Hí MATHEMATICS : ERDŐS AND ULAM PROC. N. A. S. of decomposiion, properly speaking) conradics he possibiliy of defining a counably addiive real-valu ON EQUATIONS WITH SETS AS UNKNOWNS BY PAUL ERDŐS AND S. ULAM DEPARTMENT OF MATHEMATICS, UNIVERSITY OF COLORADO, BOULDER Communicaed May 27, 1968 We shall presen here a number of resuls in se heory concerning

More information

Piecewise-Defined Functions and Periodic Functions

Piecewise-Defined Functions and Periodic Functions 28 Piecewie-Defined Funcion and Periodic Funcion A he ar of our udy of he Laplace ranform, i wa claimed ha he Laplace ranform i paricularly ueful when dealing wih nonhomogeneou equaion in which he forcing

More information

THE MATRIX-TREE THEOREM

THE MATRIX-TREE THEOREM THE MATRIX-TREE THEOREM 1 The Marix-Tree Theorem. The Marix-Tree Theorem is a formula for he number of spanning rees of a graph in erms of he deerminan of a cerain marix. We begin wih he necessary graph-heoreical

More information

Today: Max Flow Proofs

Today: Max Flow Proofs Today: Max Flow Proof COSC 58, Algorihm March 4, 04 Many of hee lide are adaped from everal online ource Reading Aignmen Today cla: Chaper 6 Reading aignmen for nex cla: Chaper 7 (Amorized analyi) In-Cla

More information

On the N P-hardness of GRacSim drawing

On the N P-hardness of GRacSim drawing Journal of Graph Algorihm and Applicaion hp://jgaa.info/ vol. 0, no. 0, pp. 0 0 (0) DOI: 10.7155/jgaa.00456 On he N P-hardne of GRacSim drawing and k-sefe Problem L. Grilli 1 1 Deparmen of Engineering,

More information