Cutting Planes for First Level RLT Relaxations of Mixed 0-1 Programs
|
|
- Samson York
- 6 years ago
- Views:
Transcription
1 Cutting Planes for First Level RLT Relaxations of Mixed 0-1 Programs 1 Cambridge, July Joint work with Franklin Djeumou Fomeni and Adam N. Letchford
2 Outline 1. Introduction 2. Literature Review 3. Cutting Planes 4. Separation 5. Computational Experiments 6. Concluding Remarks
3 Introduction The Reformulation-Linearization Technique (RLT), due to Adams and Sherali, is a technique for constructing hierarchies of LP relaxations for mixed-integer nonlinear programs. In practice, the first-level relaxation is usually the most useful. One can construct a first-level relaxation of any mixed 0-1 program with a linear or quadratic objective and any mixture of convex and quadratic constraints. We present several families of cutting planes, which can be used to strengthen such first-level relaxations.
4 Literature Review The seminal paper, Adams & Sherali (1986), was concerned with 0-1 quadratic programs of the form: { max x T Qx : Ax b, x {0,1} n}, where Q Q n n, A Q m n and b Q m. (Without loss of generality, we assume that Q is symmetric.) We will let V denote {1,...,n} and E denote {S V : S = 2}.
5 Literature Review (cont.) Let us introduce a new variable y ij for all {i,j} E, taking the value 1 if and only if both x i and x j do. The 0-1 QP can then be reformulated as the 0-1 LP: max s.t. i V Q ii x i + 2 e E Q e y e Ax b y ij x i ({i,j} E) y ij x j ({i,j} E) x i + x j y ij 1 ({i,j} E) (x,y) {0,1} V + E (In fact, this was already noted by Fortet, 1959.)
6 Literature Review (cont.) Replacing the condition (x,y) {0,1} V +E with the condition (x,y) [0,1] V +E yields an LP relaxation of the 0-1 QP. Adams and Sherali propose to strengthen the relaxation by adding 2mn additional inequalities as follows. Let α T x β be a constraint in the linear system Ax b, and let x k be one of the original variables. Then we have: (α T x)x k βx k = α i y ik (β α k )x k i V \{k} (α T x)(1 x k ) β(1 x k ) = α i (x i y ik ) β(1 x k ). i V \{k}
7 Literature Review (cont.) Sherali & Adams (1990) showed that, if the same approach is applied to a 0-1 LP, the LP relaxation of the reformulated problem dominates that of the original. They called the procedure the Reformulation-Linearization Technique. They also defined an entire hierarchy of relaxations, where the kth level RLT relaxation has variables representing products of up to k + 1 variables, and constraints that are derived by multiplying k + 1 of the original inequalities together.
8 Literature Review (cont.) There have been many developments since those early days: Lovász & Schrijver (1991) adapted the first-level relaxation to 0-1 problems with convex constraints, and also showed how to derive SDP relaxations that dominate the LP relaxation. In a series of papers, Adams & Sherali extended RLT to any bounded MINLP having a polynomial objective function and a mixture of convex and polynomial constraints. Padberg and others have studied the polytopes associated with various specific 0-1 QPs.
9 The New Cutting Planes The number of variables grows very rapidly with the level k. To keep the number of variables manageable, we concentrate on the first-level relaxation. We derive cutting planes using a two-phase approach: 1. Construct a valid cubic inequality; 2. Weaken it to make it valid for the first-level relaxation.
10 The New Cutting Planes: Phase 1 To construct valid cubic inequalities, we simply multiply triples of valid linear inequalities together. Suppose that the linear inequalities α i x β i, for i = 1,2,3, are valid for the original problem. We write them as β i α i x 0, and thereby derive the valid cubic inequality: (β 1 α 1 x)(β 2 α 2 x)(β 3 α 3 x) 0.
11 The New Cutting Planes: Phase 2 The key to phase 2 is the following lemma: Lemma 1 Suppose x i [0,1] for all i V, and y ij = x i x j for all e E. Then, for any triple {i,j,k} V, we have the following bounds on x i x j x k : x i x j x k min{y ij, y ik, y jk, 1 x i x j x k + y ij + y ik + y jk } x i x j x k max{0, y ij + y ik x i, y ij + y jk x j, y ik + y jk x k }.
12 The New Cutting Planes: Phase 2 (cont.) Due to the above lemma, we can convert a valid cubic inequality into a valid inequality for the first-level relaxation by replacing each term of the form x i x j x k with: any of 0, y ij + y ik x i, y ij + y jk x j or y ik + y jk x k (if it has a positive LHS coefficient). any of y ij, y ik, y jk or 1 x i x j x k + y ij + y ik + y jk (if the term has a negative LHS coefficient);
13 A special case: (s, t) inequalities Consider a single linear inequality α T x β and two binary variables x s and x t. We can form the following cubic inequality: (β α T x)x s x t 0. Since x 2 s x t = x s x 2 t = y st, we can re-write this as: α i x i x s x t (β α s α t )y st. i V \{s,t} To make further progress, we need more notation: for any S V, we let S + = {i S : α i > 0} and S = {i S : α i < 0}.
14 A special case: (s, t) inequalities (cont.) Then, for any partition of V into sets S,T,W and R, we can derive the (s, t) inequality: i S W α i y is + i T W α i y it α i x i α(w ) + α(s + W )x s i W + α(t + W )x t + ( β α({s,t} S + T + W R ) ) y st. Here, membership of i in S, T, W or R determines which of the weakenings x i x s x t min{y is, y it, y st, 1 x i x s x t + y st + y is + y it } x i x s x t max{0, y is + y st x s, y it + y st x t, y is + y it x i } we choose for the term x i x s x t.
15 More special cases: mixed and reverse (s, t) inequalities Following the same line of argument, we can derive: mixed (s, t) inequalities, by applying Lemma 1 to the cubic inequalities of form (β α T x)(1 x s )x t 0. reverse (s, t) inequalities, by applying Lemma 1 to the cubic inequalities of form (β α T x)(1 x s )(1 x t ) 0.
16 Remark I Using a simple disjunctive argument, it is possible to strengthen all three kinds of (s, t) inequalities. Basically, we compute upper bounds on the value that the LHS can take under the following four (mutually exclusive) conditions: x s = x t = y st = 0, x s = y st = 0 and x t = 1, x t = y st = 0 and x s = 1 x s = x t = y st = 1. We then use those bounds to tighten the RHS coefficients of x s, x t and y st.
17 Remark II The method can be applied to problems having: 1. A mixture of linear and quadratic constraints. 2. Non-linear constraints as long as they are convex.
18 Separation All three families of (s, t) inequalities are exponentially-large. So we need to think about separation. If the original problem has only linear constraints, then the separation problems for all three families can be solved in O(n 3 m) time.
19 Separation of standard (s, t) inequalities There are m ways to choose an inequality from Ax b. There are O(n 2 ) ways to select s and t. Then, for all i V +, we place i in S, T, W or R according to which of the following is largest: y is + y st x s y it + y st x t y is + y it x i 0. And, for all i V, we place i in S, T, W or R according to which of the following is largest: y is y it 1 x i x s x t + y is + y it + y st y st.
20 Separation (cont.) The separation problem for the mixed and reverse (s, t) inequalities can be tackled in a similar way. And the separation problems for the strengthened versions of all three families can be solved in O(n 3 m) time as well. For the cutting planes derived by multiplying together two linear inequalities and a variable or its complement, the separation time increases to O(n 3 m 2 ). For the ones derived by multiplying three linear inequalities, the separation time increases further to O(n 3 m 3 ).
21 Separation (cont.) Now consider the case in which there is a convex constraint x C, and we have a polynomial-time separation oracle for C (à la Lovász-Schrijver). We can show that the separation problems for all three kinds of strengthened (s, t) inequalities can be solved in polynomial time also in this case. On the other hand, we do not know how to solve the separation problem for the other inequalities in this case.
22 Computational Experiments We apply our scheme to the Quadratic Knapsack Problem: { max x T Qx : w T x c, x {0,1} n}, where Q Z n n +, w Zn + and c Z +. We randomly generate 5 instances for each different combination of size n {10,20,30,40,50} and density {25,50,75,100}. We separate all three families of strong (s, t) inequalities. We compare the resulting bounds against the benchmark upper bound given by the first-level RLT relaxation.
23 Computational Results Instances Gap (%) n RLTs RLTs + STs Improvement (%)
24 Computational Results Instances Gap (%) n RLTs RLTs + STs Improvement (%)
25 Concluding Remarks We have presented a general procedure for generating cutting planes for first-level RLT relaxations of mixed 0-1 programs. For the case in which all constraints are either linear or quadratic, we have polynomial-time separation routines for all families. For the case of convex constraints, we have polynomial-time separation routines only for the strengthened (s, t) inequalities. The computational results so far are quite encouraging.
Cutting Planes for RLT Relaxations of Mixed 0-1 Polynomial Programs
Cutting Planes for RLT Relaxations of Mixed 0-1 Polynomial Programs Franklin Djeumou Fomeni Konstantinos Kaparis Adam N. Letchford December 2013 Abstract The Reformulation-Linearization Technique (RLT),
More informationCutting Planes for RLT Relaxations of Mixed 0-1 Polynomial Programs
Cutting Planes for RLT Relaxations of Mixed 0-1 Polynomial Programs Franklin Djeumou Fomeni Konstantinos Kaparis Adam N. Letchford To Appear in Mathematical Programming Abstract The Reformulation-Linearization
More informationA new family of facet defining inequalities for the maximum edge-weighted clique problem
A new family of facet defining inequalities for the maximum edge-weighted clique problem Franklin Djeumou Fomeni June 2016 Abstract This paper considers a family of cutting planes, recently developed for
More informationA Note on Representations of Linear Inequalities in Non-Convex Mixed-Integer Quadratic Programs
A Note on Representations of Linear Inequalities in Non-Convex Mixed-Integer Quadratic Programs Adam N. Letchford Daniel J. Grainger To appear in Operations Research Letters Abstract In the literature
More informationNonconvex Quadratic Programming: Return of the Boolean Quadric Polytope
Nonconvex Quadratic Programming: Return of the Boolean Quadric Polytope Kurt M. Anstreicher Dept. of Management Sciences University of Iowa Seminar, Chinese University of Hong Kong, October 2009 We consider
More informationRLT-POS: Reformulation-Linearization Technique (RLT)-based Optimization Software for Polynomial Programming Problems
RLT-POS: Reformulation-Linearization Technique (RLT)-based Optimization Software for Polynomial Programming Problems E. Dalkiran 1 H.D. Sherali 2 1 The Department of Industrial and Systems Engineering,
More informationlinear programming and approximate constraint satisfaction
linear programming and approximate constraint satisfaction Siu On Chan MSR New England James R. Lee University of Washington Prasad Raghavendra UC Berkeley David Steurer Cornell results MAIN THEOREM: Any
More informationAn RLT Approach for Solving Binary-Constrained Mixed Linear Complementarity Problems
An RLT Approach for Solving Binary-Constrained Mixed Linear Complementarity Problems Miguel F. Anjos Professor and Canada Research Chair Director, Trottier Institute for Energy TAI 2015 Washington, DC,
More informationConvexification of Mixed-Integer Quadratically Constrained Quadratic Programs
Convexification of Mixed-Integer Quadratically Constrained Quadratic Programs Laura Galli 1 Adam N. Letchford 2 Lancaster, April 2011 1 DEIS, University of Bologna, Italy 2 Department of Management Science,
More informationSection Notes 9. IP: Cutting Planes. Applied Math 121. Week of April 12, 2010
Section Notes 9 IP: Cutting Planes Applied Math 121 Week of April 12, 2010 Goals for the week understand what a strong formulations is. be familiar with the cutting planes algorithm and the types of cuts
More informationLift-and-Project Inequalities
Lift-and-Project Inequalities Q. Louveaux Abstract The lift-and-project technique is a systematic way to generate valid inequalities for a mixed binary program. The technique is interesting both on the
More informationComparing Convex Relaxations for Quadratically Constrained Quadratic Programming
Comparing Convex Relaxations for Quadratically Constrained Quadratic Programming Kurt M. Anstreicher Dept. of Management Sciences University of Iowa European Workshop on MINLP, Marseille, April 2010 The
More informationQuadratic reformulation techniques for 0-1 quadratic programs
OSE SEMINAR 2014 Quadratic reformulation techniques for 0-1 quadratic programs Ray Pörn CENTER OF EXCELLENCE IN OPTIMIZATION AND SYSTEMS ENGINEERING ÅBO AKADEMI UNIVERSITY ÅBO NOVEMBER 14th 2014 2 Structure
More informationDecomposition Methods for Quadratic Zero-One Programming
Decomposition Methods for Quadratic Zero-One Programming Borzou Rostami This dissertation is submitted for the degree of Doctor of Philosophy in Information Technology Advisor: Prof. Federico Malucelli
More informationNotes on the decomposition result of Karlin et al. [2] for the hierarchy of Lasserre by M. Laurent, December 13, 2012
Notes on the decomposition result of Karlin et al. [2] for the hierarchy of Lasserre by M. Laurent, December 13, 2012 We present the decomposition result of Karlin et al. [2] for the hierarchy of Lasserre
More informationCSC Linear Programming and Combinatorial Optimization Lecture 12: The Lift and Project Method
CSC2411 - Linear Programming and Combinatorial Optimization Lecture 12: The Lift and Project Method Notes taken by Stefan Mathe April 28, 2007 Summary: Throughout the course, we have seen the importance
More informationHierarchies. 1. Lovasz-Schrijver (LS), LS+ 2. Sherali Adams 3. Lasserre 4. Mixed Hierarchy (recently used) Idea: P = conv(subset S of 0,1 n )
Hierarchies Today 1. Some more familiarity with Hierarchies 2. Examples of some basic upper and lower bounds 3. Survey of recent results (possible material for future talks) Hierarchies 1. Lovasz-Schrijver
More informationL R I R A P P O R T D E R E C H E R C H E KNAPSACK PROBLEM WITH PROBABILITY CONSTRAINTS GAIVORONSKI A A / LISSER A / LOPEZ R
R A P P O R T D E R E C H E R C H E L R I KNAPSACK PROBLEM WITH PROBABILITY CONSTRAINTS GAIVORONSKI A A / LISSER A / LOPEZ R Unité Mixte de Recherche 8623 CNRS-Université Paris Sud LRI 07/2008 Rapport
More informationLift-and-Project Techniques and SDP Hierarchies
MFO seminar on Semidefinite Programming May 30, 2010 Typical combinatorial optimization problem: max c T x s.t. Ax b, x {0, 1} n P := {x R n Ax b} P I := conv(k {0, 1} n ) LP relaxation Integral polytope
More informationSemidefinite Programming Basics and Applications
Semidefinite Programming Basics and Applications Ray Pörn, principal lecturer Åbo Akademi University Novia University of Applied Sciences Content What is semidefinite programming (SDP)? How to represent
More informationAM 121: Intro to Optimization! Models and Methods! Fall 2018!
AM 121: Intro to Optimization Models and Methods Fall 2018 Lecture 15: Cutting plane methods Yiling Chen SEAS Lesson Plan Cut generation and the separation problem Cutting plane methods Chvatal-Gomory
More informationCutting Planes in SCIP
Cutting Planes in SCIP Kati Wolter Zuse-Institute Berlin Department Optimization Berlin, 6th June 2007 Outline 1 Cutting Planes in SCIP 2 Cutting Planes for the 0-1 Knapsack Problem 2.1 Cover Cuts 2.2
More informationMINI-TUTORIAL ON SEMI-ALGEBRAIC PROOF SYSTEMS. Albert Atserias Universitat Politècnica de Catalunya Barcelona
MINI-TUTORIAL ON SEMI-ALGEBRAIC PROOF SYSTEMS Albert Atserias Universitat Politècnica de Catalunya Barcelona Part I CONVEX POLYTOPES Convex polytopes as linear inequalities Polytope: P = {x R n : Ax b}
More informationMoments and Positive Polynomials for Optimization II: LP- VERSUS SDP-relaxations
Moments and Positive Polynomials for Optimization II: LP- VERSUS SDP-relaxations LAAS-CNRS and Institute of Mathematics, Toulouse, France Tutorial, IMS, Singapore 2012 LP-relaxations LP- VERSUS SDP-relaxations
More informationOn Non-Convex Quadratic Programming with Box Constraints
On Non-Convex Quadratic Programming with Box Constraints Samuel Burer Adam N. Letchford July 2008 Abstract Non-Convex Quadratic Programming with Box Constraints is a fundamental N P-hard global optimisation
More informationCuts for mixed 0-1 conic programs
Cuts for mixed 0-1 conic programs G. Iyengar 1 M. T. Cezik 2 1 IEOR Department Columbia University, New York. 2 GERAD Université de Montréal, Montréal TU-Chemnitz Workshop on Integer Programming and Continuous
More informationDuality. Geoff Gordon & Ryan Tibshirani Optimization /
Duality Geoff Gordon & Ryan Tibshirani Optimization 10-725 / 36-725 1 Duality in linear programs Suppose we want to find lower bound on the optimal value in our convex problem, B min x C f(x) E.g., consider
More informationc 2000 Society for Industrial and Applied Mathematics
SIAM J. OPIM. Vol. 10, No. 3, pp. 750 778 c 2000 Society for Industrial and Applied Mathematics CONES OF MARICES AND SUCCESSIVE CONVEX RELAXAIONS OF NONCONVEX SES MASAKAZU KOJIMA AND LEVEN UNÇEL Abstract.
More informationCORC REPORT Approximate fixed-rank closures of covering problems
CORC REPORT 2003-01 Approximate fixed-rank closures of covering problems Daniel Bienstock and Mark Zuckerberg (Columbia University) January, 2003 version 2004-August-25 Abstract Consider a 0/1 integer
More informationSeparation Techniques for Constrained Nonlinear 0 1 Programming
Separation Techniques for Constrained Nonlinear 0 1 Programming Christoph Buchheim Computer Science Department, University of Cologne and DEIS, University of Bologna MIP 2008, Columbia University, New
More informationTravelling Salesman Problem
Travelling Salesman Problem Fabio Furini November 10th, 2014 Travelling Salesman Problem 1 Outline 1 Traveling Salesman Problem Separation Travelling Salesman Problem 2 (Asymmetric) Traveling Salesman
More informationCSCI 1951-G Optimization Methods in Finance Part 10: Conic Optimization
CSCI 1951-G Optimization Methods in Finance Part 10: Conic Optimization April 6, 2018 1 / 34 This material is covered in the textbook, Chapters 9 and 10. Some of the materials are taken from it. Some of
More informationOn the Polyhedral Structure of a Multi Item Production Planning Model with Setup Times
CORE DISCUSSION PAPER 2000/52 On the Polyhedral Structure of a Multi Item Production Planning Model with Setup Times Andrew J. Miller 1, George L. Nemhauser 2, and Martin W.P. Savelsbergh 2 November 2000
More informationOnline generation via offline selection - Low dimensional linear cuts from QP SDP relaxation -
Online generation via offline selection - Low dimensional linear cuts from QP SDP relaxation - Radu Baltean-Lugojan Ruth Misener Computational Optimisation Group Department of Computing Pierre Bonami Andrea
More informationA notion of Total Dual Integrality for Convex, Semidefinite and Extended Formulations
A notion of for Convex, Semidefinite and Extended Formulations Marcel de Carli Silva Levent Tunçel April 26, 2018 A vector in R n is integral if each of its components is an integer, A vector in R n is
More informationPolyhedral Approach to Integer Linear Programming. Tepper School of Business Carnegie Mellon University, Pittsburgh
Polyhedral Approach to Integer Linear Programming Gérard Cornuéjols Tepper School of Business Carnegie Mellon University, Pittsburgh 1 / 30 Brief history First Algorithms Polynomial Algorithms Solving
More informationLecture 2. Split Inequalities and Gomory Mixed Integer Cuts. Tepper School of Business Carnegie Mellon University, Pittsburgh
Lecture 2 Split Inequalities and Gomory Mixed Integer Cuts Gérard Cornuéjols Tepper School of Business Carnegie Mellon University, Pittsburgh Mixed Integer Cuts Gomory 1963 Consider a single constraint
More informationTwitter: @Owen134866 www.mathsfreeresourcelibrary.com Prior Knowledge Check 1) Simplify: a) 3x 2 5x 5 b) 5x3 y 2 15x 7 2) Factorise: a) x 2 2x 24 b) 3x 2 17x + 20 15x 2 y 3 3) Use long division to calculate:
More informationReduced rst-level representations via the reformulation-linearization technique: results, counterexamples, and computations
Discrete Applied Mathematics 101 (2000) 247 267 Reduced rst-level representations via the reformulation-linearization technique: results, counterexamples, and computations Hanif D. Sherali a;, Jonathan
More informationMoments and Positive Polynomials for Optimization II: LP- VERSUS SDP-relaxations
Moments and Positive Polynomials for Optimization II: LP- VERSUS SDP-relaxations LAAS-CNRS and Institute of Mathematics, Toulouse, France EECI Course: February 2016 LP-relaxations LP- VERSUS SDP-relaxations
More informationComparisons and enhancement strategies for linearizing mixed 0-1 quadratic programs
Discrete Optimization ( ) www.elsevier.com/locate/disopt Comparisons and enhancement strategies for linearizing mixed 0-1 quadratic programs Warren P. Adams a,, Richard J. Forrester b, Fred W. Glover c
More informationSuccessive Convex Relaxation Methods for Nonconvex Quadratic Optimization Problems
Successive Convex Relaxation Methods for Nonconvex Quadratic Optimization Problems Akiko TAKEDA Submitted in partial fulfillments of the requirement for the degree of DOCTOR OF SCIENCE Department of Mathematical
More informationarxiv: v15 [math.oc] 19 Oct 2016
LP formulations for mixed-integer polynomial optimization problems Daniel Bienstock and Gonzalo Muñoz, Columbia University, December 2014 arxiv:1501.00288v15 [math.oc] 19 Oct 2016 Abstract We present a
More informationCutting planes from extended LP formulations
Cutting planes from extended LP formulations Merve Bodur University of Wisconsin-Madison mbodur@wisc.edu Sanjeeb Dash IBM Research sanjeebd@us.ibm.com March 7, 2016 Oktay Günlük IBM Research gunluk@us.ibm.com
More informationConvex Optimization M2
Convex Optimization M2 Lecture 8 A. d Aspremont. Convex Optimization M2. 1/57 Applications A. d Aspremont. Convex Optimization M2. 2/57 Outline Geometrical problems Approximation problems Combinatorial
More informationPORTA, NEOS, and Knapsack Covers. Cover Inequalities. Prof. Jeff Linderoth. October 29, June 23, 2004 DIMACS Reconnect Conference on MIP Slide 1
PORTA, NEOS, and Knapsack Covers Cover Inequalities Prof. Jeff Linderoth October 29, 2003 June 23, 2004 DIMACS Reconnect Conference on MIP Slide 1 Today s Outline Knapsack Cover inequalities Facets Lifting
More informationSolving Box-Constrained Nonconvex Quadratic Programs
Noname manuscript No. (will be inserted by the editor) Solving Box-Constrained Nonconvex Quadratic Programs Pierre Bonami Oktay Günlük Jeff Linderoth June 13, 2016 Abstract We present effective computational
More informationMIP reformulations of some chance-constrained mathematical programs
MIP reformulations of some chance-constrained mathematical programs Ricardo Fukasawa Department of Combinatorics & Optimization University of Waterloo December 4th, 2012 FIELDS Industrial Optimization
More informationCORC Technical Report TR Cuts for mixed 0-1 conic programming
CORC Technical Report TR-2001-03 Cuts for mixed 0-1 conic programming M. T. Çezik G. Iyengar July 25, 2005 Abstract In this we paper we study techniques for generating valid convex constraints for mixed
More informationEasy and not so easy multifacility location problems... (In 20 minutes.)
Easy and not so easy multifacility location problems... (In 20 minutes.) MINLP 2014 Pittsburgh, June 2014 Justo Puerto Institute of Mathematics (IMUS) Universidad de Sevilla Outline 1 Introduction (In
More informationSequential pairing of mixed integer inequalities
Sequential pairing of mixed integer inequalities Yongpei Guan, Shabbir Ahmed, George L. Nemhauser School of Industrial & Systems Engineering, Georgia Institute of Technology, 765 Ferst Drive, Atlanta,
More informationCutting Plane Methods I
6.859/15.083 Integer Programming and Combinatorial Optimization Fall 2009 Cutting Planes Consider max{wx : Ax b, x integer}. Cutting Plane Methods I Establishing the optimality of a solution is equivalent
More informationA General Framework for Convex Relaxation of Polynomial Optimization Problems over Cones
Research Reports on Mathematical and Computing Sciences Series B : Operations Research Department of Mathematical and Computing Sciences Tokyo Institute of Technology 2-12-1 Oh-Okayama, Meguro-ku, Tokyo
More informationIntroduction to LP and SDP Hierarchies
Introduction to LP and SDP Hierarchies Madhur Tulsiani Princeton University Local Constraints in Approximation Algorithms Linear Programming (LP) or Semidefinite Programming (SDP) based approximation algorithms
More informationOn the knapsack closure of 0-1 Integer Linear Programs. Matteo Fischetti University of Padova, Italy
On the knapsack closure of 0-1 Integer Linear Programs Matteo Fischetti University of Padova, Italy matteo.fischetti@unipd.it Andrea Lodi University of Bologna, Italy alodi@deis.unibo.it Aussois, January
More informationOn Valid Inequalities for Quadratic Programming with Continuous Variables and Binary Indicators
On Valid Inequalities for Quadratic Programming with Continuous Variables and Binary Indicators Hongbo Dong and Jeff Linderoth Wisconsin Institutes for Discovery University of Wisconsin-Madison, USA, hdong6,linderoth}@wisc.edu
More information6-1 The Positivstellensatz P. Parrilo and S. Lall, ECC
6-1 The Positivstellensatz P. Parrilo and S. Lall, ECC 2003 2003.09.02.10 6. The Positivstellensatz Basic semialgebraic sets Semialgebraic sets Tarski-Seidenberg and quantifier elimination Feasibility
More informationE5295/5B5749 Convex optimization with engineering applications. Lecture 5. Convex programming and semidefinite programming
E5295/5B5749 Convex optimization with engineering applications Lecture 5 Convex programming and semidefinite programming A. Forsgren, KTH 1 Lecture 5 Convex optimization 2006/2007 Convex quadratic program
More informationThe Trust Region Subproblem with Non-Intersecting Linear Constraints
The Trust Region Subproblem with Non-Intersecting Linear Constraints Samuel Burer Boshi Yang February 21, 2013 Abstract This paper studies an extended trust region subproblem (etrs in which the trust region
More informationarxiv: v2 [math.oc] 3 Nov 2016
Compact Linearization for Binary Quadratic Problems subject to Assignment Constraints arxiv:1610.05375v2 [math.oc] 3 Nov 2016 Sven Mallach Institut für Informatik Universität zu Köln, 50923 Köln, Germany
More informationTwo-stage stochastic (and distributionally robust) p-order conic mixed integer programs: Tight second stage formulations
Two-stage stochastic (and distributionally robust p-order conic mixed integer programs: Tight second stage formulations Manish Bansal and Yingqiu Zhang Department of Industrial and Systems Engineering
More informationDeciding Emptiness of the Gomory-Chvátal Closure is NP-Complete, Even for a Rational Polyhedron Containing No Integer Point
Deciding Emptiness of the Gomory-Chvátal Closure is NP-Complete, Even for a Rational Polyhedron Containing No Integer Point Gérard Cornuéjols 1 and Yanjun Li 2 1 Tepper School of Business, Carnegie Mellon
More informationLecture: Examples of LP, SOCP and SDP
1/34 Lecture: Examples of LP, SOCP and SDP Zaiwen Wen Beijing International Center For Mathematical Research Peking University http://bicmr.pku.edu.cn/~wenzw/bigdata2018.html wenzw@pku.edu.cn Acknowledgement:
More informationIntroduction to Semidefinite Programming I: Basic properties a
Introduction to Semidefinite Programming I: Basic properties and variations on the Goemans-Williamson approximation algorithm for max-cut MFO seminar on Semidefinite Programming May 30, 2010 Semidefinite
More informationLagrange Duality. Daniel P. Palomar. Hong Kong University of Science and Technology (HKUST)
Lagrange Duality Daniel P. Palomar Hong Kong University of Science and Technology (HKUST) ELEC5470 - Convex Optimization Fall 2017-18, HKUST, Hong Kong Outline of Lecture Lagrangian Dual function Dual
More informationLMI Methods in Optimal and Robust Control
LMI Methods in Optimal and Robust Control Matthew M. Peet Arizona State University Lecture 02: Optimization (Convex and Otherwise) What is Optimization? An Optimization Problem has 3 parts. x F f(x) :
More informationNumerical Optimization. Review: Unconstrained Optimization
Numerical Optimization Finding the best feasible solution Edward P. Gatzke Department of Chemical Engineering University of South Carolina Ed Gatzke (USC CHE ) Numerical Optimization ECHE 589, Spring 2011
More informationCopositive Programming and Combinatorial Optimization
Copositive Programming and Combinatorial Optimization Franz Rendl http://www.math.uni-klu.ac.at Alpen-Adria-Universität Klagenfurt Austria joint work with I.M. Bomze (Wien) and F. Jarre (Düsseldorf) IMA
More informationHypergraph Matching by Linear and Semidefinite Programming. Yves Brise, ETH Zürich, Based on 2010 paper by Chan and Lau
Hypergraph Matching by Linear and Semidefinite Programming Yves Brise, ETH Zürich, 20110329 Based on 2010 paper by Chan and Lau Introduction Vertex set V : V = n Set of hyperedges E Hypergraph matching:
More informationLecture 18: Optimization Programming
Fall, 2016 Outline Unconstrained Optimization 1 Unconstrained Optimization 2 Equality-constrained Optimization Inequality-constrained Optimization Mixture-constrained Optimization 3 Quadratic Programming
More informationSherali-Adams Relaxations of Graph Isomorphism Polytopes. Peter N. Malkin* and Mohammed Omar + UC Davis
Sherali-Adams Relaxations of Graph Isomorphism Polytopes Peter N. Malkin* and Mohammed Omar + UC Davis University of Washington May 18th 2010 *Partly funded by an IBM OCR grant and the NSF. + Partly funded
More informationApplications of semidefinite programming in Algebraic Combinatorics
Applications of semidefinite programming in Algebraic Combinatorics Tohoku University The 23rd RAMP Symposium October 24, 2011 We often want to 1 Bound the value of a numerical parameter of certain combinatorial
More informationl p -Norm Constrained Quadratic Programming: Conic Approximation Methods
OUTLINE l p -Norm Constrained Quadratic Programming: Conic Approximation Methods Wenxun Xing Department of Mathematical Sciences Tsinghua University, Beijing Email: wxing@math.tsinghua.edu.cn OUTLINE OUTLINE
More informationOn the Lovász Theta Function and Some Variants
On the Lovász Theta Function and Some Variants Laura Galli Adam N. Letchford March 2017. To appear in Discrete Optimization. Abstract The Lovász theta function of a graph is a well-known upper bound on
More information5 Integer Linear Programming (ILP) E. Amaldi Foundations of Operations Research Politecnico di Milano 1
5 Integer Linear Programming (ILP) E. Amaldi Foundations of Operations Research Politecnico di Milano 1 Definition: An Integer Linear Programming problem is an optimization problem of the form (ILP) min
More informationTightening linearizations of non-linear binary optimization problems
Tightening linearizations of non-linear binary optimization problems Elisabeth Rodríguez Heck and Yves Crama QuantOM, HEC Management School, University of Liège Partially supported by Belspo - IAP Project
More informationConvex Optimization. (EE227A: UC Berkeley) Lecture 6. Suvrit Sra. (Conic optimization) 07 Feb, 2013
Convex Optimization (EE227A: UC Berkeley) Lecture 6 (Conic optimization) 07 Feb, 2013 Suvrit Sra Organizational Info Quiz coming up on 19th Feb. Project teams by 19th Feb Good if you can mix your research
More informationLecture #21. c T x Ax b. maximize subject to
COMPSCI 330: Design and Analysis of Algorithms 11/11/2014 Lecture #21 Lecturer: Debmalya Panigrahi Scribe: Samuel Haney 1 Overview In this lecture, we discuss linear programming. We first show that the
More informationLP formulations for mixed-integer polynomial optimization problems Daniel Bienstock and Gonzalo Muñoz, Columbia University, December 2014
LP formulations for mixed-integer polynomial optimization problems Daniel Bienstock and Gonzalo Muñoz, Columbia University, December 2014 Abstract We present a class of linear programming approximations
More informationPolynomial and Synthetic Division
Polynomial and Synthetic Division Polynomial Division Polynomial Division is very similar to long division. Example: 3x 3 5x 3x 10x 1 3 Polynomial Division 3x 1 x 3x 3 3 x 5x 3x x 6x 4 10x 10x 7 3 x 1
More informationRelations between Semidefinite, Copositive, Semi-infinite and Integer Programming
Relations between Semidefinite, Copositive, Semi-infinite and Integer Programming Author: Faizan Ahmed Supervisor: Dr. Georg Still Master Thesis University of Twente the Netherlands May 2010 Relations
More informationA finite branch-and-bound algorithm for nonconvex quadratic programming via semidefinite relaxations
Math. Program., Ser. A (2008) 3:259 282 DOI 0.007/s007-006-0080-6 FULL LENGTH PAPER A finite branch-and-bound algorithm for nonconvex quadratic programming via semidefinite relaxations Samuel Burer Dieter
More informationOperations Research Lecture 6: Integer Programming
Operations Research Lecture 6: Integer Programming Notes taken by Kaiquan Xu@Business School, Nanjing University May 12th 2016 1 Integer programming (IP) formulations The integer programming (IP) is the
More informationCopositive Programming and Combinatorial Optimization
Copositive Programming and Combinatorial Optimization Franz Rendl http://www.math.uni-klu.ac.at Alpen-Adria-Universität Klagenfurt Austria joint work with M. Bomze (Wien) and F. Jarre (Düsseldorf) and
More informationAn Introduction to Linear Matrix Inequalities. Raktim Bhattacharya Aerospace Engineering, Texas A&M University
An Introduction to Linear Matrix Inequalities Raktim Bhattacharya Aerospace Engineering, Texas A&M University Linear Matrix Inequalities What are they? Inequalities involving matrix variables Matrix variables
More informationThere are several approaches to solve UBQP, we will now briefly discuss some of them:
3 Related Work There are several approaches to solve UBQP, we will now briefly discuss some of them: Since some of them are actually algorithms for the Max Cut problem (MC), it is necessary to show their
More informationLecture 9: Dantzig-Wolfe Decomposition
Lecture 9: Dantzig-Wolfe Decomposition (3 units) Outline Dantzig-Wolfe decomposition Column generation algorithm Relation to Lagrangian dual Branch-and-price method Generated assignment problem and multi-commodity
More informationThe Retail Planning Problem Under Demand Uncertainty. UCLA Anderson School of Management
The Retail Planning Problem Under Demand Uncertainty George Georgiadis joint work with Kumar Rajaram UCLA Anderson School of Management Introduction Many retail store chains carry private label products.
More informationNetwork Flows. 6. Lagrangian Relaxation. Programming. Fall 2010 Instructor: Dr. Masoud Yaghini
In the name of God Network Flows 6. Lagrangian Relaxation 6.3 Lagrangian Relaxation and Integer Programming Fall 2010 Instructor: Dr. Masoud Yaghini Integer Programming Outline Branch-and-Bound Technique
More informationComputational Integer Programming. Lecture 2: Modeling and Formulation. Dr. Ted Ralphs
Computational Integer Programming Lecture 2: Modeling and Formulation Dr. Ted Ralphs Computational MILP Lecture 2 1 Reading for This Lecture N&W Sections I.1.1-I.1.6 Wolsey Chapter 1 CCZ Chapter 2 Computational
More informationarxiv:math/ v1 [math.co] 23 May 2000
Some Fundamental Properties of Successive Convex Relaxation Methods on LCP and Related Problems arxiv:math/0005229v1 [math.co] 23 May 2000 Masakazu Kojima Department of Mathematical and Computing Sciences
More informationminimize x x2 2 x 1x 2 x 1 subject to x 1 +2x 2 u 1 x 1 4x 2 u 2, 5x 1 +76x 2 1,
4 Duality 4.1 Numerical perturbation analysis example. Consider the quadratic program with variables x 1, x 2, and parameters u 1, u 2. minimize x 2 1 +2x2 2 x 1x 2 x 1 subject to x 1 +2x 2 u 1 x 1 4x
More informationThe Traveling Salesman Problem: Inequalities and Separation
The Traveling Salesman Problem: Inequalities and Separation Adam N. Letchford Department of Management Science, Lancaster University http://www.lancs.ac.uk/staff/letchfoa 1. The ILP formulation of the
More information13. Convex programming
CS/ISyE/ECE 524 Introduction to Optimization Spring 2015 16 13. Convex programming Convex sets and functions Convex programs Hierarchy of complexity Example: geometric programming Laurent Lessard (www.laurentlessard.com)
More informationExtended Linear Formulation for Binary Quadratic Problems
Noname manuscript No. (will be inserted by the editor) Extended Linear Formulation for Binary Quadratic Problems Fabio Furini Emiliano Traversi Received: date / Accepted: date Abstract We propose and test
More information3.7 Strong valid inequalities for structured ILP problems
3.7 Strong valid inequalities for structured ILP problems By studying the problem structure, we can derive strong valid inequalities yielding better approximations of conv(x ) and hence tighter bounds.
More information15.083J/6.859J Integer Optimization. Lecture 2: Efficient Algorithms and Computational Complexity
15.083J/6.859J Integer Optimization Lecture 2: Efficient Algorithms and Computational Complexity 1 Outline Efficient algorithms Slide 1 Complexity The classes P and N P The classes N P-complete and N P-hard
More informationAn exact reformulation algorithm for large nonconvex NLPs involving bilinear terms
An exact reformulation algorithm for large nonconvex NLPs involving bilinear terms Leo Liberti CNRS LIX, École Polytechnique, F-91128 Palaiseau, France. E-mail: liberti@lix.polytechnique.fr. Constantinos
More informationPolyhedral Results for A Class of Cardinality Constrained Submodular Minimization Problems
Polyhedral Results for A Class of Cardinality Constrained Submodular Minimization Problems Shabbir Ahmed and Jiajin Yu Georgia Institute of Technology A Motivating Problem [n]: Set of candidate investment
More informationDistributed Computation of Minimum Time Consensus for Multi-Agent Systems
Distributed Computation of Minimum Time Consensus for Multi-Agent Systems Ameer Mulla, Deepak Patil, Debraj Chakraborty IIT Bombay, India Conference on Decision and Control 2014 Los Angeles, California,
More information