Dual bounds: can t get any better than...

Similar documents
to work with) can be solved by solving their LP relaxations with the Simplex method I Cutting plane algorithms, e.g., Gomory s fractional cutting

Discrete (and Continuous) Optimization WI4 131

and to estimate the quality of feasible solutions I A new way to derive dual bounds:

3.4 Relaxations and bounds

Integer Programming ISE 418. Lecture 8. Dr. Ted Ralphs

Integer Programming: Cutting Planes

Relaxations and Bounds. 6.1 Optimality and Relaxations. Suppose that we are given an IP. z = max c T x : x X,

Section Notes 8. Integer Programming II. Applied Math 121. Week of April 5, expand your knowledge of big M s and logical constraints.

Course Notes for MS4315 Operations Research 2

BBM402-Lecture 20: LP Duality

- Well-characterized problems, min-max relations, approximate certificates. - LP problems in the standard form, primal and dual linear programs

Integer Programming, Part 1

Planning and Optimization

LP Duality: outline. Duality theory for Linear Programming. alternatives. optimization I Idea: polyhedra

Discrete Optimization 2010 Lecture 7 Introduction to Integer Programming

Discrete Optimization 2010 Lecture 8 Lagrangian Relaxation / P, N P and co-n P

Introduction to Mathematical Programming IE406. Lecture 10. Dr. Ted Ralphs

3.10 Lagrangian relaxation

Maximum flow problem

Lecture 5 January 16, 2013

Lectures 6, 7 and part of 8

Cutting Plane Methods II

(P ) Minimize 4x 1 + 6x 2 + 5x 3 s.t. 2x 1 3x 3 3 3x 2 2x 3 6

AM 121: Intro to Optimization! Models and Methods! Fall 2018!

CS599: Convex and Combinatorial Optimization Fall 2013 Lecture 17: Combinatorial Problems as Linear Programs III. Instructor: Shaddin Dughmi

3. Linear Programming and Polyhedral Combinatorics

Discrete (and Continuous) Optimization Solutions of Exercises 2 WI4 131

Linear Programming Duality P&S Chapter 3 Last Revised Nov 1, 2004

Duality of LPs and Applications

CS675: Convex and Combinatorial Optimization Fall 2016 Combinatorial Problems as Linear and Convex Programs. Instructor: Shaddin Dughmi

Chapter 7 Network Flow Problems, I

MVE165/MMG631 Linear and integer optimization with applications Lecture 8 Discrete optimization: theory and algorithms

Modeling with Integer Programming

CS675: Convex and Combinatorial Optimization Fall 2014 Combinatorial Problems as Linear Programs. Instructor: Shaddin Dughmi

Linear Programming Duality

6. Linear Programming

where X is the feasible region, i.e., the set of the feasible solutions.

CO 250 Final Exam Guide

Linear and Integer Programming - ideas

Semidefinite and Second Order Cone Programming Seminar Fall 2001 Lecture 5

Linear Programming. Scheduling problems

Algorithmic Game Theory and Applications. Lecture 7: The LP Duality Theorem

MVE165/MMG630, Applied Optimization Lecture 6 Integer linear programming: models and applications; complexity. Ann-Brith Strömberg

The Primal-Dual Algorithm P&S Chapter 5 Last Revised October 30, 2006

Technische Universität München, Zentrum Mathematik Lehrstuhl für Angewandte Geometrie und Diskrete Mathematik. Combinatorial Optimization (MA 4502)

Discrete Optimization

4. Duality Duality 4.1 Duality of LPs and the duality theorem. min c T x x R n, c R n. s.t. ai Tx = b i i M a i R n

IE 5531: Engineering Optimization I

EXERCISES SHORTEST PATHS: APPLICATIONS, OPTIMIZATION, VARIATIONS, AND SOLVING THE CONSTRAINED SHORTEST PATH PROBLEM. 1 Applications and Modelling

Lecture 7: Lagrangian Relaxation and Duality Theory

3.7 Strong valid inequalities for structured ILP problems

15.081J/6.251J Introduction to Mathematical Programming. Lecture 24: Discrete Optimization

Advanced Algorithms 南京大学 尹一通

Optimization for Communications and Networks. Poompat Saengudomlert. Session 4 Duality and Lagrange Multipliers

Fall 2017 Qualifier Exam: OPTIMIZATION. September 18, 2017

Chapter 1 Linear Programming. Paragraph 5 Duality

Simplex method(s) for solving LPs in standard form

Integer Programming ISE 418. Lecture 12. Dr. Ted Ralphs

3.3 Easy ILP problems and totally unimodular matrices

Integer Programming ISE 418. Lecture 16. Dr. Ted Ralphs

3. Linear Programming and Polyhedral Combinatorics

Introduction to Mathematical Programming IE406. Lecture 21. Dr. Ted Ralphs

min 4x 1 5x 2 + 3x 3 s.t. x 1 + 2x 2 + x 3 = 10 x 1 x 2 6 x 1 + 3x 2 + x 3 14

Integer and Combinatorial Optimization: Introduction

Mathematical Programs Linear Program (LP)

14. Duality. ˆ Upper and lower bounds. ˆ General duality. ˆ Constraint qualifications. ˆ Counterexample. ˆ Complementary slackness.

Week 8. 1 LP is easy: the Ellipsoid Method

Week Cuts, Branch & Bound, and Lagrangean Relaxation

Introduction to integer programming II

Lecture 15 (Oct 6): LP Duality

In the original knapsack problem, the value of the contents of the knapsack is maximized subject to a single capacity constraint, for example weight.

Lecture #21. c T x Ax b. maximize subject to

Acyclic Semidefinite Approximations of Quadratically Constrained Quadratic Programs

ACO Comprehensive Exam March 17 and 18, Computability, Complexity and Algorithms

Travelling Salesman Problem

Algorithms and Theory of Computation. Lecture 13: Linear Programming (2)

CS 6820 Fall 2014 Lectures, October 3-20, 2014

SOLVING INTEGER LINEAR PROGRAMS. 1. Solving the LP relaxation. 2. How to deal with fractional solutions?

1 Overview. 2 Extreme Points. AM 221: Advanced Optimization Spring 2016

IP Duality. Menal Guzelsoy. Seminar Series, /21-07/28-08/04-08/11. Department of Industrial and Systems Engineering Lehigh University

New Integer Programming Formulations of the Generalized Travelling Salesman Problem

Compact vs. exponential-size LP relaxations

Workforce Scheduling. Outline DM87 SCHEDULING, TIMETABLING AND ROUTING. Outline. Workforce Scheduling. 1. Workforce Scheduling.

3.7 Cutting plane methods

EE364a Review Session 5

MATTHIAS GERDTS. COMBINATORIAL OPTIMISATION MSM 3M02b

1 Integer Decomposition Property

Convex Optimization and SVM

Lecture 8: Column Generation

3.8 Strong valid inequalities

Decision Procedures An Algorithmic Point of View

Computational Integer Programming Universidad de los Andes. Lecture 1. Dr. Ted Ralphs

Week 4. (1) 0 f ij u ij.

Advanced Linear Programming: The Exercises

MIT Algebraic techniques and semidefinite optimization February 14, Lecture 3

Mathematical Programming II Lecture 5 ORIE 6310 Spring 2014 February 7, 2014 Scribe: Weixuan Lin

A notion of Total Dual Integrality for Convex, Semidefinite and Extended Formulations

Cutting Plane Methods I

Lecture 5. x 1,x 2,x 3 0 (1)

The use of shadow price is an example of sensitivity analysis. Duality theory can be applied to do other kind of sensitivity analysis:

Transcription:

Bounds, relaxations and duality Given an optimization problem z max{c(x) x 2 }, how does one find z, or prove that a feasible solution x? is optimal or close to optimal? I Search for a lower and upper bound on the optimal objective value z apple z apple z. I If z z, then the optimal objective value is found. I If x? 2 satisfies z apple c(x? ) apple z, and z z apple, where > 0 is an optimality tolerance parameter, then x? is guaranteed to satisfy c(x? ) z. Primal bounds can you beat this? Lower(upper) bounds on the optimal objective function value of a max(min) problem Any feasible solution can be used to obtain a primal bound. IOE 51 Introduction to IP, Winter 2012 Bounds, relaxations, duality Page 37 c Marina A. Epelman Dual bounds can t get any better than... Upper(lower) bounds on the optimal objective function value of a max(min) problem are often obtained via relaxations of the problem Relaxation (IP) z max{c(x) x 2 R n } A problem (RP) z R max{f (x) x 2 R n } is a relaxation of (IP) if (i),and(ii) f (x) c(x) forallx 2. Proposition 2.1 If (RP) is a relaxation of (IP), then z R z. Proposition 2.3 (i) If a relaxation (RP) is infeasible, the original problem (IP) is infeasible. (ii) Let x? be an optimal solution of (RP). If x? 2 and f (x? )c(x? ), then x? is an optimal solution of (IP). IOE 51 Introduction to IP, Winter 2012 Bounds, relaxations, duality Page 3 c Marina A. Epelman

Linear Programming Relaxations LP realaxation For the integer program max{c x x 2 P \ Z n } with formulation P {x 2 R n Ax apple b}, thelinear programming relaxation is the linear program z LP max{c x x 2 P}. Clearly a relaxation P \ Z n P and objective function is unchanged. Proposition 2.2 Suppose, P 1, P 2 are two formulations for the IP max{c x x 2 Z n } with P 1 P 2.If max{c x x 2 P i }, i 1, 2, then z1 LP apple z2 LP for all c. z LP i If P is the ideal formulation, then z LP z?, and any extreme point optimal solution of the LP is an optimal solution of the IP. IOE 51 Introduction to IP, Winter 2012 Bounds, relaxations, duality Page 3 c Marina A. Epelman Combinatorial relaxations SP SP a SP tour is an assignment (or permutation) containing no subtours; hence z SP min c ij is a tour ; min c ij is an assignment ; zass IOE 51 Introduction to IP, Winter 2012 Bounds, relaxations, duality Page 40 c Marina A. Epelman

Combinatorial relaxations symmetric SP SSP aspwithc ij c ji ree and 1-tree A tree on a set of k nodes is an undirected graph that satisfies any two of the following three properties (i) It has k 1 edges (ii) It contains no cycles (iii) It is connected. A 1-tree on n nodes is a graph consisting of two edges adjacent to node 1, plus the edges of a tree on nodes {2,...,n}. Every tour is a 1-tree, so z SSP min min c ij is a tour ; c ij is a 1-tree ; z1-tree IOE 51 Introduction to IP, Winter 2012 Bounds, relaxations, duality Page 41 c Marina A. Epelman Lagrangian relaxation (IP) z max{c x Ax apple b, x 2 Z n }. Proposition 2.4 Let z(u) max{c x + u (b all u 0. Ax) x 2 }. henz apple z(u) for Above is called the Lagrangian relaxation of (IP). It is, indeed a relaxation larger feasible region, and for any x feasible for (IP), z(u) c x + u (b Ax) c x. IOE 51 Introduction to IP, Winter 2012 Bounds, relaxations, duality Page 42 c Marina A. Epelman

Duality a dual problem Duality he two problems (IP) z max{c x x 2 Z n },and (D) w min{w(u) u 2 U} form a (weak)-dual pair if c(x) apple w(u) forallx 2 and all u 2 U. Whenz w, theyform a strong-dual pair. Advantage any feasible solution of the dual provides a dual bound (whereas a relaxation needs to be solved to optimality ). Proposition 2.6 (i) If (D) is unbounded, the original problem (IP) is infeasible. (ii) If x? 2 and u? 2 U satisfy c(x? )w(u? ), the x? is optimal for (IP) and w? is optimal for (D). Proposition 2.5 he integer program z max{c x Ax apple b, x 2 Z n +} and the linear program w LP min{u b u A c, u 2 R m +} form a weak dual pair. IOE 51 Introduction to IP, Winter 2012 Bounds, relaxations, duality Page 43 c Marina A. Epelman A dual pair combinatorial perspective Given a graph G (V, E), a matching M E is a set of disjoint edges; a covering by nodes is a set R V of nodes such that every edge has at least one endpoint in R. Proposition 2.7 he problem of finding a maximum cardinality matching and the problem of finding a minimum cardinality covering by nodes form a weak-dual pair. Proof If M {(i 1, j 1 ),...,(i k, j k )} then the nodes i 1, j 1,...,i k, j k are distinct, and any covering by nodes R must contain at least one node from each pair. Hence R k M. IOE 51 Introduction to IP, Winter 2012 Bounds, relaxations, duality Page 44 c Marina A. Epelman

A dual pair LP perspective Node-edge incidence matrix Let n V and m E. Defineann by m node-edge incidence matrix A as follows a j,e 1ifnodej is an endpoint of edge e, and 0 otherwise. Maximum cardinality matching z max{e x Ax apple e, x 2 Z m +} Minimum cardinality cover by nodes w min{e y A y e, y 2 Z n +} Let z LP and w LP be the values of the corresponding LP relaxations. hen z apple z LP w LP apple w, establishing weak-duality. 2 Not a strong-dual pair! Consider A 4 0 1 1 1 0 1 1 1 0 3 5 IOE 51 Introduction to IP, Winter 2012 Bounds, relaxations, duality Page 45 c Marina A. Epelman