Linear Programming (Optimization)
Download
Report
Transcript Linear Programming (Optimization)
II.2
Strong Valid Inequalities and Facets for
Structured Integer Programs
1
1. Introduction
ο± Strength of the valid inequalities is important (want facet or high dimensional
face). Use the structure of the problem to obtain strong valid inequalities.
ο± Node packing problem: π = π₯ β π΅π : π₯π + π₯π β€ 1, β π, π β πΈ
dim(conv(π)) = π (π unit vector + 0 vector)
If πΆ β π is a clique,
πβπΆ π₯π
If πΆ is a maximal clique,
β€ 1 is a valid inequality.
πβπΆ π₯π
β€ 1 defines a facet of conv(π).
Pf) enough to show that there exist π affinely independent feasible points
satisfying the clique constraint at equality.
Let πΆ = {1, β¦ , π} be a maximal clique. For each node π β πΆ, β π(π) β πΆ
which is not connected to π since πΆ is a maximal clique.
Consider the characteristic vectors of the packings, {1}, β¦ , {π}, {π + 1, π(π +
1)}, β¦ , {π, π(π)}.
Integer Programming 2013
2
2
πΊ:
2
π»:
6
3
6
3
1
5
4
5
4
ο± π₯ 1 = 12(0, 1, 1, 1, 1, 1) is an extreme point of the initial polytope (clique
constraints and nonnegativities).
To cut off π₯ 1 , consider odd hole (chordless cycle, odd cycle with no crossing
edges). Let π» β π, H induces an odd hole, then
πβπ» π₯π
β€ ( π» β 1)/2 is a valid inequality. (π₯2 + π₯3 + π₯4 + π₯5 + π₯6 β€ 2)
Moreover, it gives a facet of the convex hull of node packings for the
subgraph with node set π» (five linearly independent vectors satisfying the
inequality at equality). But, it does not give a facet of conv(π) for πΊ.
Integer Programming 2013
3
ο± Consider Ξ±π₯1 + (π₯2 + π₯3 + π₯4 + π₯5 + π₯6 ) β€ 2.
Can we make ο‘ > 0 and large while the inequality remains valid?
When π₯1 = 0, it is valid for any ο‘ > 0.
When π₯1 = 1, πΌ β€ 2 - (π₯2 + π₯3 + π₯4 + π₯5 + π₯6 ). But π₯1 = 1 implies π₯2 =
β― = π₯6 = 0, so πΌ β€ 2.
Therefore, 2π₯1 + (π₯2 + π₯3 + π₯4 + π₯5 + π₯6 ) β€ 2 is a valid inequality.
Moreover, (1, 0, 0, 0, 0, 0) vector in addition to the five vectors we already
have are linearly independent and satisfy the constraint at equality βΉ facet.
ο± Lifting procedure: Systematic method to obtain facet (high dimensional face)
from a facet for the low dimensional restricted problem
Integer Programming 2013
4
ο± Prop 1.1: Suppose π β π΅π , π πΏ = πβ{π₯ β π΅π : π₯1 = πΏ}, πΏ = 0,1 , and
π
π=2 ππ π₯π
1
β€ π0 (1.7) is valid for π 0 . If π1 = β
, then π₯1 β€ 0 is valid for S. If
π β β
, then πΌ1 π₯1 + ππ=2 ππ π₯π β€ π0 (1.8) is valid for π for any πΌ1 β€ π0 β π,
where π = max{ ππ=2 ππ π₯π : π₯ β π1 }.
Moreover, if πΌ1 = π0 β π and (1.7) gives a face of dimension π of conv(π0),
then (1.8) gives a face of dimension π + 1 of conv(π). [If (1.7) facet of
conv(π0), (1.8) facet of conv(π).]
Pf) If π₯ β π 0 βΉ πΌ1 π₯1 +
If π₯ β π1 βΉ πΌ1 π₯1 +
π
π=2 ππ π₯π
π
π=2 ππ π₯π
β€ π0 valid.
= πΌ1 +
π
π=2 ππ π₯π
β€ πΌ1 + π β€ π0 valid.
(1.7) gives a π βdimensional face of conv(π0) βΉ β π₯ π β π 0 such that π₯1π =
0, π = 1, β¦ , π + 1. Also π₯ π satisfies πΌ1 π₯1 + ππ=2 ππ π₯π β€ π0 at equality.
Let π =
π
β
π=2 ππ π₯π ,
π₯ β β π1
πΌ1 = π0 β π βΉ π₯ β satisfies πΌ1 π₯1 +
π1
π
π=2 ππ π₯π
= π0 and π₯1β = 1 since π₯ β β
βΉ π₯ β cannot be written as affine combination of {π₯ 1 , β¦ , π₯ π+1 }, so the π + 2
vectors {π₯ β , π₯ 1 , β¦ , π₯ π+1 } are affinely independent.
οΏ
Integer Programming 2013
5
ο± Lifting is also possible from π1 to S.
ο± Prop 1.2: Suppose (1.7) is valid for π1 . If π 0 = β
, then π₯1 β₯1 is valid for π.
If π 0 β β
, then πΎ1 π₯1 + ππ=2 ππ π₯π β€ π0 + πΎ1 (1.9) is valid for π for any πΎ1 β₯
π β π0 , where π = max{ ππ=2 ππ π₯π : π₯ β π 0 }.
Moreover, if πΎ1 = π β π0 and (1.7) gives a face of dimension π of conv(π1),
then (1.9) gives a face of dimension π + 1 of conv(π). [If (1.7) facet of
conv(π0), (1.8) facet of conv(π).]
ο± When πΌ1 = π0 β π in Prop 1.1 or when πΎ1 = π β π0 in Prop 1.2, we say that
the lifting is maximum.
ο± Lifting is done sequentially one variable at a time and the strength of the lifted
inequality depends on the sequence of the variables lifted.
ο±
Integer Programming 2013
6
ο± Prop 1.3: Let π\π1 = {1, 2, β¦ , π‘} and suppose when Prop 1.1 is applied
sequentially using maximum lifting in the order (π1 , π2 , β¦ , ππβ1 , ππ , β¦ , ππ‘ ), the
inequality πβπ\π1 πΌπ π₯π + πβπ1 ππ π₯π β€ π0 is obtained.
β²
Then for any order (π1β² , π2β² , β¦ , ππβ1
, ππβ² , β¦ , ππ‘β² ) with ππβ² = ππ for π = 1, β¦ , π β 1, the
resulting inequality πβπ\π1 πΌπβ² π₯π + πβπ1 ππ π₯π β€ π0 obtained by maximum
lifting has πΌπβ²π β€ πΌππ .
Pf) ππ = ππ β² for some π > π. Then
πΌππ = π0 β max{
0 for π > π}
= π0 β max{
πβ1
π=1 πΌππ π₯ππ
πβ1
π=1 πΌππ π₯ππ
+
π : π₯ = 1 and π₯ =
π
π₯
βΆ
π₯
β
π
β©
{π₯
β
π΅
π
π
ππ
ππ
πβπ1
+
πβπ1 ππ π₯π
: π₯ β π β© {π₯ β π΅π : π₯ππ β² = 1 πππ π₯π β² =
+
πβπ1 ππ π₯π
+
π
0 for π β€ π β€ π‘, π β π }
β₯ π0 β max{
πβ1
π=1 πΌππ π₯ππ
π β1 β²
π=π πΌππβ² π₯ππβ²
: π₯ β π β© {π₯ β π΅π : π₯ππ β² =
1 πππ π₯π β² = 0 for π > π }
π
= πΌπβ²π
οΏ
Integer Programming 2013
7