Lagrangean relaxation

Transcription

1 Lagrangean relaxation Giovanni Righini Corso di Complementi di Ricerca Operativa

2 Joseph Louis de la Grange (Torino Paris 1813)

3 Relaxations Given a problem P, such as: minimize z P (x) s.t. x X P a problem R, such as: minimize z R (x) s.t. x X R is a relaxation of P if and only if these two conditions hold: X P X R z R (x) z P (x) x X P. As a consequence z R z P. The optimal value of the relaxation is not worse than the optimal value of the original problem.

4 Duality Given a problem P, such as: minimize z P (x) s.t. x X its dual problem D has the form: maximize w D (y) s.t. y Y and the following property holds (Weak Duality Theorem): x X, y Y w D (y) z P (x). In particular y Y w D (y) z P.

5 Dual bounds Combinatorial optimization problems are in general NP-hard and we usually resort to implicit enumeration to solve them to optimality (this is also useful for approximation purposes). Two main methods for implicit enumeration are: branch-and-bound, dynamic programming. In branch-and-bound algorithms we must associate a dual bound with every sub-problem. There are two main techniques to compute a dual bound to a given (sub-)problem: solving a relaxation to optimality (z R z P ), finding a feasible solution for its dual (sub-)problem ( y Y w D (y) z P ).

6 Consider this problem P: Lagrangean relaxation P)minimize z P (x) =c T x A Lagrangean relaxation of P is: s.t. A 1 x b 1 A 2 x b 2 x 0 integer LR)minimize z LR (x, λ) =c T x+λ T (b 1 A 1 x) s.t. A 2 x b 2 x 0 integer Violations of the relaxed constraints are penalized in the objective function by means of Lagrangean multipliers λ 0.

7 Is it a relaxation? P)min z P (x) =c T x s.t. A 1 x b 1 A 2 x b 2 x 0 integer LR)min z LR (x, λ) =c T x + λ T (b 1 A 1 x) s.t. A 2 x b 2 x 0 integer We indicate with Q the feasible region of the relaxed problem. X P Q because constraints A 1 x b 1 have been relaxed. x X P A 1 x b 1. So λ 0, λ T (b 1 A 1 x) 0 and hence z LR (x, λ) z P (x).

8 Variations Lagrangean relaxation can also be applied to equality constraints. In this case λ is unrestricted in sign. Lagrangean relaxation can also be applied to non-linear constraints.

9 Relaxations, relaxations, relaxations... Relaxations can be obtained in many different ways. Linear relaxation is: P)minimize z P (x) =c T x s.t. A 1 x b 1 A 2 x b 2 x 0 integer LP)minimize z LP (x) =c T x s.t. A 1 x b 1 A 2 x b 2 x 0

10 Relaxations, relaxations, relaxations... Relaxations can be obtained in many different ways. Lagrangean relaxation is: P)minimize z P (x) =c T x s.t. A 1 x b 1 A 2 x b 2 x 0 integer LR)minimize z LR (x, λ) =c T x + λ T (b 1 A 1 x) s.t. A 2 x b 2 x 0 integer

11 Relaxations, relaxations, relaxations... Relaxations can be obtained in many different ways. Surrogate relaxation is: P)minimize z P (x) =c T x s.t. A 1 x b 1 SR)minimize z SR (x) =c T x A 2 x b 2 x 0 integer s.t. λ T A 1 x λ T b 1 A 2 x b 2 x 0 integer

12 Lagrangean dual problem For each choice of the Lagrangean multipliers λ 0, we obtain a different instance, providing a different dual bound. Example (Set Covering). The Lagrangean dual problem asks for the best Lagrangean multipliers λ : LD) zld = max λ 0 {z LR (λ) = min {z LR(x, λ) = c T x + λ T (b 1 A 1 x)}}. x Q This gives us the best dual bound that can be obtained from the Lagrangean relaxation of constraints A 1 x b 1. But a duality gap may still exist between z LD and z P.

13 Complementary slackness If the optimal solution x for problem LR is also feasible for P, this does not guarantee that x is also optimal for P. This is true only if the penalty terms are null. Primal optimality requires both: feasibility: x X P, i.e. A 1 x b 1 complementarity: λ T (b 1 A 1 x ) = 0. When equality constraints are relaxed, the two conditions coincide.

14 Example (1) P) max z P (x) =7x 1 + 2x 2 s.t. x 1 + 2x 2 4 5x 1 + x x 1 2x 2 7 x 1 2 x 2 4 x Z 2 + x 2 conv(x) conv(q) X Q x 1

15 Example (2) LR) max z LR (x, λ) = =(7+λ)x 1 + (2 2λ)x 2 +4λ s.t. 5x 1 + x x 1 2x 2 7 x 1 2 x 2 4 x Z 2 + Q = {(2, 2),..., (4, 0)}. Q is the polyhedron defined by the linear constraints of LR. x 2 conv(x) conv(q) X Q x 1

16 Primal viewpoint The objective function z LR (x, λ) can be seen as... a function of x for λ fixed; a function of λ for x fixed. Primal viewpoint. For each given λ 0 we have z LR (λ) = max{z LR(x, λ) : x Q}. Since constraints and objectives are linear: z LR(λ) = max{z LR (x, λ) : x conv(q)}. This is a linear programming problem, where the orientation of the objective function depends on λ.

17 Primal viewpoint (example). For 0 λ 1 9 zlr (λ) = z((3, 4), λ) = 29 λ. x 2 Q For λ 1 9 zlr (λ) = z((4, 0), λ) = λ. Hence z LD = min λ 0 {z LR (λ)} = z LR ( 1 9 ) = z LR((3, 4), 1 9 ) = z LR ((4, 0), 1 9 ) = conv(x) conv(q) X x 1

18 Dual viewpoint For each given x i Q, z LR (x i, λ) is an affine function of λ. Hence z LR (λ) = max{z LR(x i, λ) : x i Q, λ 0} is piecewise linear and convex. z LR (λ) = z LR((3, 4), λ) for 0 λ 1 9 z LR (λ) = z LR((4, 0), λ) for λ /9 z LR(x i, ) z LR((4,0), ) z LR((3,1), ) z LR((3,2), ) z LR((3,3), ) Then z LR((3,4), ) zld = min{z LR (λ) : λ 0} = = z LR( 1 9 ) = z LR((2,2), ) z LR((2,3), ) 18 1/9 z LR((2,4), )

19 In our example when λ = 1 9 we have z LR (λ) = Consider the point x = 8 9 (3, 4) (4, 0), which is a convex combination of two integer points of Q. Point x is the intersection point between the segment joining the two integer points and the relaxed constraint = z LR((3, 4), 1 9 ) = z LR((4, 0), 1 9 ) = z LR(x, 1 9 ) = In general: = c T x (4 ( x 1 + 2x 2 )) = ct x. z LD = max{ct x : A 1 x b 1, x conv(q)}. There exists a convex combination x of points of Q, which satisfies the relaxed constraints and such that the value of the objective function in x is equal to zld. Finding this convex combination which satisfies the relaxed constraints and maximizes the original objective function is a linear problem and it is the dual of the Lagrangean dual.

20 How good is Lagrangean relaxation? Let define R = {x R n + : A 1x b 1 }. conv(x P ) = conv(q R) conv(q) R Hence: z P = max{ct x : x X P } = max{c T x : x conv(x P )} = max{c T x : x conv(q R)} max{c T x : x conv(q) R} = z LD. The duality gap z LD z P depends on the relationship between X P and Q (and then R) and from the coefficients c. z LD = z P c if and only if conv(q R) = conv(q) R. In our example they do not coincide: z LD = and z P = 28.

21 Lagrangean vs. linear relaxation When the relaxed problem is a discrete one, it is interesting to compare the dual bound zld that can be obtained from Lagrangean relaxation with the dual bound zlp which is (easily) obtained from linear relaxation. zld = z LP c if all extreme points of the polyhedron {x R n + : A 2 x b 2 } are integer, i.e. Q = conv(q). In our example this does not occur and we have z LP = > z LD = If we relaxed all constraints but the last two, we would obtain a feasible region Q with integer extreme points. In that case we would have z LP = z LD =

22 Lagrangean vs. linear relaxation Summing up: and then conv(x P ) conv(q) R {x R n + : Ax b} z P z LD z LP. In general the best dual bound given by Lagrangean relaxation is tighter than the dual bound given by linear relaxation. They are equivalent when the relaxed problem has the integrality property (all integral extreme points).

23 Special cases For some particular orientation of the objective function the bounds can be the same even if the general conditions for equivalence are not satisfied. With c 1 : z P = z LD = z LP. With c 2 : z P < z LD = z LP. With c 3 : z P < z LD < z LP. With c 4 : z P = z LD < z LP. x 2 c 4 X conv(x) Q c 3 c 2 conv(q) c 1 x 1

24 Approximation properties z LD z P ɛ λ 0, x X P : λ T (b 1 A 1 x) δ 1 z LR (λ) z LR(x, λ) δ 2 δ 1 + δ 2 ɛ. Therefore we can have z LD = z P only if there exist λ 0 and x X P such that these two conditions hold: Optimality of the solution of LR: z LR (λ) = z LR(x, λ) Complementary slackness: λ T (b 1 A 1 x) = 0. We can evaluate the distance from optimality if we know δ 1 : the error due to the choice of x in the LR δ 2 : the error due to the choice of λ.

25 Lagrangean relaxation and problem decomposition (1) A typical application of Lagrangean relaxation is to relax coupling constraints: this allows to decompose the problem into independent sub-problems. Example: Uncapacitated facility location problem. M N N min z P (x, y) = c ij x ij + f j y j i=1 j=1 j=1 N s.t. x ij = ( )1 j=1 x ij y j y B n i = 1,..., M i = 1,..., M j = 1,..., N

26 Lagrangean relaxation and problem decomposition (2) Relaxing the first set of constraints we have penalties of the form M i=1 λ i(1 N j=1 x ij): M N N N min z LR (x, y, λ) = (c ij λ i )x ij + f j y j + i=1 j=1 s.t. x ij y j y B n j=1 i=1 λ i i = 1,..., M j = 1,..., N This can be decomposed into N independent subproblems: M min z j (x, y, λ) = (c ij λ i )x ij + f j y j i=1 s.t. x ij y j x ij 0 y j {0, 1} i = 1,..., M i = 1,..., M that can be easily solved.

27 Which relaxation should we choose? There is a trade-off between three criteria: How many constraints do we relax? (n. of Lagrangean multipliers) How difficult the relaxed problem is? (computing time) Does the relaxed problem have the integrality property? (tightness of the dual bound)

28 An example Consider for instance the Generalized Assignment Problem. M N min z P (x) = c ij x ij i=1 j=1 N s.t. x ij = ( )1 j=1 M a ij x ij b j i=1 x B M N i = 1,..., M j = 1,..., N We can relax the red constraints, the blue constraints or both: we get three different relaxations.

29 Lagrangean relaxation 1 Penalty terms of constraints N j=1 x ij = ( )1 i = 1,..., M have the form: M N λ i (1 x ij ). Then we obtain: i=1 j=1 M N M min zlr(x, 1 λ) = (c ij λ i )x ij + i=1 j=1 i=1 M s.t. a ij x ij b j i=1 x B M N λ i j = 1,..., N This relaxation: requires M multipliers, implies solving N knapsack problem instances, does not have the integrality property.

30 Lagrangean relaxation 2 Penalty terms of constraints M i=1 a ijx ij b j j = 1,..., N have the form: N M µ j ( a ij x ij b j ). Then we obtain: min z 2 LR(x, µ) = This relaxation: s.t. j=1 i=1 j=1 i=1 M N N (c ij +a ij µ j )x ij µ j b j N x ij = ( )1 j=1 x B M N requires N multipliers, j=1 implies solving M trivial problem instances, does have the integrality property. i = 1,..., M

31 Lagrangean relaxation 3 Relaxing both sets of constraints, we obtain: M N M N min zlr(x, 3 λ, µ) = (c ij λ i +a ij µ j )x ij + λ i µ j b j i=1 j=1 i=1 j=1 s.t. x B M N This relaxation: requires M + N multipliers, implies solving M + N very trivial problem instances, does have the integrality property.

32 How do we solve the Lagrangean dual problem? The Lagrangean dual problem consists of minimizing a piecewise linear and convex function in a space with as many dimensions as the number of Lagrangean multipliers. Convexity: good news. There are no sub-optimal local minima. Piecewise linearity: bad news. First order partial derivatives (gradient vector) are not continuous everywhere. It is a linear programming problem with as many constraints as the number of integer points in Q. To solve this problem we can use one of several approximation methods: sub-gradient optimization, multiplier adjustment, dual ascent.

33 Sub-gradient optimization (1) It is a local search algorithm that iteratively moves from a current point λ (k) to a next point λ (k+1) with a step of length σ (k) in the direction opposite of the sub-gradient b Ax (λ (k) ). Begin Inizialize λ (0) ; k:=0; while not EndTest do Solve max z LR (x, λ (k) ) = c T x + λ (k)t (b Ax) s.t. x Q; Obtain x (λ (k) ); Choose σ (k) ; λ (k+1) := max{λ (k) σ (k) (b Ax (λ (k) )), 0}; k := k + 1; end while End The sub-gradient algorithm is not monotonic: it is necessary to save the best value encountered so far.

34 Sub-gradient optimization (2) Initialization is usually not important. The algorithm is not sensitive to it. The end test succeeds when no improvements are observed after a given number of iterations. The step can be chosen with some heuristic rule. Rule of thumb (a). If lim k k i=1 σ(i) = and lim k σ (k) = 0, then the sub-gradient algorithm converges to the optimal value of the Lagrangean dual. Rule of thumb (b). If σ (k) = σ (0) ρ k with ρ < 1, then it converges to optimality provided that µ (0) and ρ are large enough.

35 Beasley s rules of thumb: Sub-gradient optimization (3) Normalize the right hand sides of the relaxed constraints to 1; Choose 0 < π 2; Choose λ (0) arbitrarily; Define a scalar step σ = π z z LR (λ(0),x ) i max{g i, where,0} 2 z is the value of a feasible solution (possibly increased by a small percentage); x is the optimal solution of LR(λ (0) ); zlr(λ (0), x ) is the corresponding optimal value; i is the index spanning all relaxed constraints; G i = b i x ; n j=1 a ijxj Update λ i := max{0, λ i + σg i }; is the sub-gradient of constraint i evaluated in Halve π when the best value does not improve for a certain number of iterations; Stop after a prescribed number of iterations of when π becomes small enough.

36 Lagrangean heuristics The information provided by the optimal solution of the Lagrangean relaxation (corresponding to a dual bound), is used to build a primal feasible solution. Example (Set Covering). Lagrangean heuristics must be tailored to each specific problem and relaxation. This allows to compute primal bounds every time dual bounds are computed in a branch-and-bound algorithm.

37 Variable fixing After solving a Lagrangean relaxation, we have: a reduced cost for each variable, a gap between a primal and a dual bound. If the reduced cost is larger than the gap, the current value of the variable can be fixed. Example (Set covering).

38 Given a problem Lagrangean decomposition (1) min c T x s.t. A 1 x b 1 A 2 x b 2 x B n we can transform it into the equivalent problem min c T x + (c c ) T y s.t. A 1 x b 1 A 2 y b 2 x = y x, y B n Now we can relax the coupling constraint x = y and decompose.

39 Lagrangean decomposition (2) From the Lagrangean relaxation of the coupling constraints (equalities) we have: which decomposes into min c T x + (c c ) T y + λ T (x y) s.t. A 1 x b 1 A 2 y b 2 x, y B n min (c + λ) T x s.t. A 1 x b 1 x B n min (c c λ) T y s.t. A 2 y b 2 y B n