70 likes | 239 Views
6. Polyhedral Ties between Linear and Integer Programs. P = { xR + n : Ax b} S = P Z n , want conv (S) is a rational polyhedron P bounded S finite or conv (S) is polyhedron ( Weyl ) If P unbounded, may exist infinite number of points in S. Can we have Weyl type result?
E N D
6. Polyhedral Ties between Linear and Integer Programs • P = {xR+n: Ax b} S = P Zn, want conv(S) is a rational polyhedron P bounded S finite or conv(S) is polyhedron (Weyl) If P unbounded, may exist infinite number of points in S. Can we have Weyl type result? • Thm 6.1: P = {xR+n: Ax b} , S = P Zn, where (A, b) is an integer m(n+1) matrix, then the following statements are true. • finite {ql}lL of S and rays {rj}jJ of P such thatS = {xR+n: x = lL lql+ jJjrj, l=1, Z+|L|, Z+|J|} • If P is a cone, rays {vh}hH of P such that S = {xR+n: x = hH hvh, Z+|H|}.
Pf) I.P = {xR+n: x = kKkxk + jJjrj, k=1, k, j0, for kK and jJ} Assume {rj} for jJ are integer vectors. Let Q = {xZ+n: x = kK kxk + jJ jrj, k=1, k 0 for kK, 0j<1 for jJ} Q is a finite set, say Q = {qlZ+n: l L} and Q S. Observe that xi S xi Z+n and xi ={kKkixk+ jJ(ji - ji)rj}+{jJjirj}, k=1, k, j0, for kK and jJ. • xi = ql(i) + jJjirj, ji= j for all jJ. II. If P is a cone, ql S implies ql S Z+1. Take {vh : hH} = {ql : l L} {rj : j J}
Thm 6.2: P = {xR+n: Ax b}, S = P Zn conv(S) is a rational polyhedron. Pf) Convex combination of points {xiS, iI} can be written as x = iI ixi = iI i(ql(i) + jJ jirj) (iI i = 1) = lL({iI: l(i)=l}i)ql + jJ(iI iji)rj = lL lql + jJ jrj, Where l = {iI: l(i)=l}i 0 for lL, lLl = iIi = 1, and j =iI iji 0 for jJ. conv(S) = {xR+n: x = lL lql + jJ jrj, lLl= 1, l, j 0 for lL and jJ}, with ql, rj Z+n for lL and jJ. By Thm 4.13 (Weyl), conv(S) is a rational polyhedron. • Proof extends straightforwardly to mixed-integer sets with rational data. Hence, all of the following results apply to mixed-integer sets. • Extreme rays of P = {xR+n: Ax b} and conv(P Zn) coincide.
Thm 6.2 suggests that we can solve (IP) max{cx: xS, S = P Zn by solving the LP (CIP) max{cx: xconv(S)}. • Thm 6.3: Given S = PZn , P={xR+n: Ax b} and cRn, it follows that: • IP unbounded CIP unbounded. • If CIP has a bounded optimal value It has an optimal solution (namely, an extreme point of conv(S)) that is an optimal solution to IP. • IP has an optimal solution x0 x0 is an optimal solution to CIP. Pf) Let z0 and z* be the optimal values of IP and CIP, respectively. z* z0 since conv(S) S. a) z0 = z* = z* = extreme point x0conv(S) and ray rZ+n such that cr > 0. x0 + r conv(S) 0 x0 + r S Z+1 z0 = . b) CIP has an optimal solution extreme point optimal solution, say x0 x0 S z0 cx0 = z* z0 = z* (since z* z0) c) x0 optimal to IP CIP not unbounded (from a.) and x0 conv(S) . Hence, x0 is an optimal solution to CIP (from b.)
Cor 6.4: IP is either infeasible or unbounded or has an optimal solution. • We say that (, 0) is a valid inequality for a set S if x 0 for all xS. • Prop 6.5: If x 0 is valid for S, it is also valid for conv(S). Pf) Consider an xconv(S). Then x = jJ jxj, where xj S for jJ and jJj = 1 and j 0 for jJ. Hence x = jJj(xj) jJj0 = 0.
Prop 6.6: x 0 face of dimension k-1 of conv(S) k affinely independent points x1, … , xk S such that xi = 0, i = 1, … , k. Pf) affinely independent x1, … , xk conv(S) such that xi = 0, i =1, … ,k. Suppose x1S, x1 = jJ jx*j, x*jS, j>0 for all jJ, jJ j=1. x1= 0, x*j 0 for jJ x*j= 0 for all jJ x1, … , xkaffinely independent x*j* such that x*j*, x2, … , xkaffinely independent. Replace x1 with x*j* and repeat the procedure. • Relate IP to LP: z(d) = max{cx: Ax d, xZ+n} zLP(d) = max{cx: Ax d, xR+n}
Prop 6.7: • zLP(0) = z(0) • z(0) = 0 Q = {uR+n: uA c} • z(0) = Q = • If Q , then S = P Zn = or z(b) finite • If Q = , then S = or z(b) = Pf) Clearly 0 z(0) zLP(0) (primal feasible) If Q zLP(0) = 0 z(0) = zLP(0) = 0 If Q = ray rj of P with crj > 0. can take rj to be integer. z(0) = zLP(0) = (a, b, c) If Q P = or zLP(b) is finite, hence z(b) finite (d) If Q = P = or zLP(b) = zLP(0) = If S z(0) = z(b) = (from c) (e) • Cor 6.8: • P = S = • zLP(b) finite S = or z(b) finite • zLP(b) = S = or z(b) =