Mixed Integer Programming (MIP) for Causal Inference and Beyond
|
|
- Osborn Campbell
- 5 years ago
- Views:
Transcription
1 Mixed Integer Programming (MIP) for Causal Inference and Beyond Juan Pablo Vielma Massachusetts Institute of Technology Columbia Business School New York, NY, October, 2016.
2 Traveling Salesman Problem (TSP): Visit Cities Fast
3 Traveling Salesman Problem (TSP): Visit Cities Fast Paradoxes, Contradictions, and the Limits of Science Many research results define boundaries of what cannot be known, predicted, or described. Classifying these limitations shows us the structure of science and reason. Noson S. Yanofsky S A computer would have to check all these possible routes to find the shortest one.
4 MIP = Avoid Enumeration Number of tours for 49 cities Fastest supercomputer flops Assuming one floating point operation per tour: > years times the age of the universe! How long does it take on an iphone? Less than a second! 4 iterations of cutting plane method! = 48!/ Dantzig, Fulkerson and Johnson 1954 did it by hand! For more info see tutorial in ConcordeTSP app Cutting planes are the key for effectively solving (even NPhard) MIP problems in practice.
5 Using IP to visit Germany 45 cities(1832) 120 cities(1977) 15,112 cities(2004)
6 50+ Years of MIP = Significant Solver Speedups Algorithmic Improvements (Machine Independent): CPLE v1.2 (1991) v11 (2007): 29,000x speedup Gurobi v1 (2009) v6.5 (2015): 48.7x speedup Commercial, but free for academic use (Reasonably) effective free / open source solvers: GLPK, CBC and SCIP (free only for non-commercial) Easy to use, fast and versatile modeling languages Julia based JuMP modelling language Linear MIP solvers very mature and effective: Convex nonlinear MIP getting there (quadratic nearly there)
7 Matching Treated Units: T = {t 1,...,t T } Control Units: C = {c 1,...,c C } Observed Covariates: P = {p 1,...,p P } Covariate Values: x t = x t p p2p, x c = x c p p2p, t 2 T c 2 C
8 Maximum Cardinality Exact Matching max s.t. t2t c2c m t,c m t,c apple 1, 8c 2 C t2t m t,c apple 1, c2c 8t 2 T m t,c =0 8t, c x t 6= x c 0 apple m t,c apple 1 m t,c 2 {0, 1} 8t 2 T, c2 C. Solve time for truncated Lalonde CPS 429 = s Why? Can solve relaxation.
9 Maximum Cardinality Marginal Means t2t max c2c m t,c s.t. m t,c apple 1, t2t m t,c apple 1, c2c t2t 8c 2 C 8t 2 T x c pm t,c 8p 2 P x t pm t,c = c2c c2c t2t m t,c 2 {0, 1} 8t 2 T, c2 C. Solve time for truncated Lalonde CPS 429 = 444 s Why? One reason = relaxation has fractions.
10 Maximum Cardinality Fine Balance, Take 1 max s.t. t2t c2c m t,c t2t mp t,c apple 1, 8c 2 C, p 2 P c2c mp t,c apple 1, 8t 2 T, p 2 P m p t,c =0 8t, c x t p 6= x c p, p 2 P t2t mp t,c = t2t mq t,c c2c mp t,c = c2c mq t,c 8c 2 C, p, q 2 P 8t 2 T, p, q 2 P m p t,c 2 {0, 1} 8t 2 T,c2 C, p 2 P. Solve time for truncated Lalonde CPS 429 = 0.81 s Why? One reason = relaxation has fewer fractions.
11 Basic Branch-and-Bound x (1) FRAC z = 1.71 x 1 apple 0 x (2) INT z = 0 x 2 2 (3) FRAC z = 1.5 x 1 apple x 1 3x 1 +2x 2 apple 6 2x 1 + x 2 apple 0 x 1,x 2 0 (4) INT z = 1 Linear Programming (LP) Relaxation (5) INF x 1,x 2 2 Z
12 Stronger Formulations = Faster Solves? x x (1) FRAC z = 1.71 x 1 apple 0 x (2) INT z = 0 (3) INT z = x x 1 3x 1 +2x 2 apple 6 2x 1 + x 2 apple 0 x 1,x 2 0 x 1,x 2 2 Z x 1 +x 2 apple 2 2x 1 +x 2 apple 0 x 1,x 2 0 x 1,x 2 2 Z
13 Cutting Plane Example: Chátal-Gomory Cuts ( ) P := x 2 R 2 x 1 + x 2 apple 3, : 5x 1 3x 2 apple 3 3 H := x 2 R 2 :4x 1 +3x 2 apple 10.5 {z } 2 Z if x 2 Z Valid for H \ Z
14 Branch-and-Bound and Cuts (Branch-and-Cut) x 2 3 (1) FRAC z = (2) FRAC z = 1 x 1 apple 0 x x 1 3x 1 +2x 2 apple 6 (3) INT z = 0 (4) INT z = 1 2x 1 + x 2 apple 0 x 1,x 2 0 x 1,x 2 2 Z
15 Partial Solves and GAP x 2 3 (1) FRAC z = 1.71 x 1 apple 0 x (2) INT z = 0 (3) FRAC z = x 1 GAP =100 bestnode bestinteger bestinteger BB GAP = = BB + GAP = = 50%
16 Maximum Cardinality Marginal Means t2t max c2c m t,c s.t. m t,c apple 1, t2t m t,c apple 1, c2c t2t 8c 2 C 8t 2 T x c pm t,c 8p 2 P x t pm t,c = c2c c2c t2t m t,c 2 {0, 1} 8t 2 T, c2 C. Solve time for truncated Lalonde CPS 429 = 444 s Optimal = 17, LP Relaxation = 19.35
17 Maximum Cardinality Marginal Means Root relaxation: objective e+01, 1947 iterations, 0.16 seconds Nodes Current Node Objective Bounds Work Expl Unexpl Obj Depth IntInf Incumbent BestBd Gap It/Node Time Black Magic ( s s s s s H % s... * % s Cutting planes: Gomory: 1 Flow cover: 1 Explored nodes ( simplex iterations) in seconds Optimal solution found (tolerance 1.00e-04) Best objective e+01, best bound e+01, gap 0.0% Better Feasible Solution Optimal Solution Bound
18 Maximum Cardinality Fine Balance, Take 1 max s.t. t2t c2c m t,c t2t mp t,c apple 1, 8c 2 C, p 2 P c2c mp t,c apple 1, 8t 2 T, p 2 P m p t,c =0 8t, c x t p 6= x c p, p 2 P t2t mp t,c = t2t mq t,c c2c mp t,c = c2c mq t,c 8c 2 C, p, q 2 P 8t 2 T, p, q 2 P m p t,c 2 {0, 1} 8t 2 T,c2 C, p 2 P. Solve time for truncated Lalonde CPS 429 = 0.81 s Optimal = 10, LP Relaxation = 11
19 Maximum Cardinality Fine Balance, Take 1 Root relaxation: objective e+01, 1490 iterations, 0.03 seconds Nodes Current Node Objective Bounds Work Expl Unexpl Obj Depth IntInf Incumbent BestBd Gap It/Node Time s s s H % - 0s % - 0s Cutting planes: Zero half: 2 Explored 8 nodes (9933 simplex iterations) in 0.78 seconds Optimal solution found (tolerance 1.00e-04) Best objective e+01, best bound e+01, gap 0.0% Better Feasible Solution Optimal Solution Bound
20 Maximum Cardinality Fine Balance, Take 1 max s.t. t2t c2c m t,c t2t mp t,c apple 1, 8c 2 C, p 2 P c2c mp t,c apple 1, 8t 2 T, p 2 P m p t,c =0 8t, c x t p 6= x c p, p 2 P t2t mp t,c = t2t mq t,c c2c mp t,c = c2c mq t,c 8c 2 C, p, q 2 P 8t 2 T, p, q 2 P m p t,c 2 {0, 1} 8t 2 T,c2 C, p 2 P. Solve time for truncated Lalonde CPS 429 = 0.81 s Scalability? Size = T C P
21 Maximum Cardinality Fine Balance, Take 2 max s.t. t2t p,k t2t c2c m t,c m t,c apple 1, t2t m t,c apple 1, c2c c/2c p,k m t,c = t/2t p,k Solve time for truncated Lalonde CPS 429 = 0.61 s Scalability? Size = c2c p,k m t,c 8c 2 C 8t 2 T 8p 2 P,k 2 K(p) m t,c 2 {0, 1} 8t 2 T, c 2 C. T C + P K(p) ={x c p} c2p [ {x t p} t2t C p,k = {c 2 C : x c p = k} T p,k = {t 2 T : x t p = k}
22 Take 1 Revisited max s.t. t2t c2c m t,c t2t mp t,c apple 1, 8c 2 C, p 2 P c2c mp t,c apple 1, 8t 2 T, p 2 P m p t,c =0 8t, c x t p 6= x c p, p 2 P t2t mp t,c = t2t mq t,c c2c mp t,c = c2c mq t,c 8c 2 C, p, q 2 P 8t 2 T, p, q 2 P m p t,c 2 {0, 1} 8t 2 T,c2 C, p 2 P. Solve time for truncated Lalonde CPS 429 = 0.81 s Scalability? Size = T C P
23 Take 1 Revisited max s.t. t2t x t m p t,c = y c, 8p 2 P, k 2 K(p), c 2 C p,k t2t p,k m p t,c = x t, 8S 2 F, k 2 K(p), t 2 T p,k c2c p,k m p t,c 2 {0, 1} 8p 2 P, k 2 K(p), t 2 T p,k, c 2 C p,k x t 2 {0, 1} 8t 2 T y c 2 {0, 1} 8c 2 C.
24 Maximum Cardinality Fine Balance, Take 3 max s.t. t2t t2t x t t2t p,k x t = x t = c2c y c, Solve time for truncated Lalonde CPS 429 = s Scalability? Size = c2c p,k y c, 8p 2 P, k 2 K(p) x t 2 {0, 1} 8t 2 T y c 2 {0, 1} 8c 2 C. P ( T + C ) K(p) ={x c p} c2p [ {x t p} t2t C p,k = {c 2 C : x c p = k} T p,k = {t 2 T : x t p = k}
25 Simple Formulation for Univariate Functions f(d 3 ) f(d 2 ) f(d 5 ) f(d 1 ) f(d 4 ) z = f(x) d 1 d 2 d 3 d 4 d 5 Size = O (# of segments) Non-Ideal: Fractional Extreme Points x z = 5 1= 5 j=1 j=1 dj f(d j ) j j, j 0
26 Advanced Formulation for Univariate Functions f(d 3 ) f(d 2 ) f(d 5 ) f(d 1 ) f(d 4 ) z = f(x) x z d 1 d 2 d 3 d 4 d 5 Size = O (log 2 # of segments) Ideal: Integral Extreme Points = 5 1= 5 j=1 j=1 y 2 {0, 1} 2 dj f(d j ) j j, j 0
27 Extended Formulation for PWL Functions S = gr (f) = f(d 3 ) f(d 2 ) f(d 1 ) f(d 4 ) k[ i=1 ( (x, z) 2 R 2 : 0 d 1 d 2 d 3 d 4 d i apple x apple d i+1 m i x + c i = z ) d i y i apple x i apple d i+1 y i m i x i + c i y i = z i k i=1 xi = x k i=1 zi = z k i=1 y i =1 MC Formulation: y 2 {0, 1} k 8i 2 [k] 8i 2 [k]
28 Abstracting Univariate Functions P z = f(x) = i := 2 5 : j =0 8j 5/2 T i T j=1 i := {i, i +1} i 2 {1,...,4} f(d 3 ) 5 = f(d 2 ) f(d 5 ) f(d 1 ) f(d 4 ) x z 1= 5 j=1 dj f(d j ) j j, j 0 d 1 d 2 d 3 d 4 d 5 T [ 4 i=1 P i 5
29 Abstraction Works for Multivariate Functions P i := { 2 m : j =0 8v j /2 T i } f(x,y) v m y x 2 [ n i=1 P i m
30 Complete Abstraction V := n 2 R V + : v2v o v =1, P i = 2 V : v =0 8v /2 T i T i = cliques of a graph
31 From Cliques to (Complement) Conflict Graph SOS
32 From Conflict Graph to Bi-clique Cover SOS
33 From Bi-clique Cover to Formulation SOS2 0 apple apple 1 y 1 0 apple 3 apple y 1 0 apple apple 1 y 2 0 apple apple y
34 Ideal Formulation from Bi-clique Cover Conflict Graph G =(V,E) E = {(u, v) :u, v 2 V, u 6= s.t. u, v 2 T i } Bi-clique cover A j, B j t j=1, Aj, B j V 8{u, v} 2 E 9j s.t. u 2 A j ^ v 2 B j Formulation v2a j v apple 1 y j 8j 2 [t] v2b j v apple y j 8j 2 [t] y 2 {0, 1} t
35 Recursive Construction of Cover for SOS2, Step 1 Base case n=2 1 : Step 1 recursion : Reflect Graph / Cover Stick Graph / Cover Repeat for all bi-cliques from 2 k-1 to cover all edges within first and last half of conflict graph
36 Recursive Construction of Cover for SOS2, Step 2 Only edges missing are those between first and last half of conflict graph Step 2 : Add one more bi-clique 1 2 n/2 n/2 n/2 + n n Cover has log 2 n bi-cliques. For non-power of two just delete extra nodes.
37 Grid Triangulations: Step 1 = SOS2 for Inter-Box Covers all arcs between boxes
38 Grid Triangulations: Step 2 = Ad-hoc Intra-Box Covers all arcs within boxes Sometimes 1 additional cover
39 Grid Triangulations: Step 2 = Ad-hoc Intra-Box Sometimes 2 additional covers Sometimes more, but always less than 9 Simple rules to get (near) optimal in Fall 16
Advanced Mixed Integer Programming (MIP) Formulation Techniques
Advanced Mixed Integer Programming (MIP) Formulation Techniques Juan Pablo Vielma Massachusetts Institute of Technology Center for Nonlinear Studies, Los Alamos National Laboratory. Los Alamos, New Mexico,
More informationAdvanced Mixed Integer Programming Formulations for Non-Convex Optimization Problems in Statistical Learning
Advanced Mixed Integer Programming Formulations for Non-Convex Optimization Problems in Statistical Learning Juan Pablo Vielma Massachusetts Institute of Technology 2016 IISA International Conference on
More informationMixed Integer Programming (MIP) for Daily Fantasy Sports, Statistics and Marketing
Mixed Integer Programming (MIP) for Daily Fantasy Sports, Statistics and Marketing Juan Pablo Vielma Massachusetts Institute of Technology AM/ES 121, SEAS, Harvard. Boston, MA, November, 2016. MIP & Daily
More informationMixed Integer Programming Approaches for Experimental Design
Mixed Integer Programming Approaches for Experimental Design Juan Pablo Vielma Massachusetts Institute of Technology Joint work with Denis Saure DRO brown bag lunch seminars, Columbia Business School New
More information18 hours nodes, first feasible 3.7% gap Time: 92 days!! LP relaxation at root node: Branch and bound
The MIP Landscape 1 Example 1: LP still can be HARD SGM: Schedule Generation Model Example 157323 1: LP rows, still can 182812 be HARD columns, 6348437 nzs LP relaxation at root node: 18 hours Branch and
More informationEncodings in Mixed Integer Linear Programming
Encodings in Mixed Integer Linear Programming Juan Pablo Vielma Sloan School of usiness, Massachusetts Institute of Technology Universidad de hile, December, 23 Santiago, hile. Mixed Integer er inary Formulations
More informationIntroduction to Integer Programming
Lecture 3/3/2006 p. /27 Introduction to Integer Programming Leo Liberti LIX, École Polytechnique liberti@lix.polytechnique.fr Lecture 3/3/2006 p. 2/27 Contents IP formulations and examples Total unimodularity
More informationThe Traveling Salesman Problem: An Overview. David P. Williamson, Cornell University Ebay Research January 21, 2014
The Traveling Salesman Problem: An Overview David P. Williamson, Cornell University Ebay Research January 21, 2014 (Cook 2012) A highly readable introduction Some terminology (imprecise) Problem Traditional
More informationMixed Integer Programming Models for Non-Separable Piecewise Linear Cost Functions
Mixed Integer Programming Models for Non-Separable Piecewise Linear Cost Functions Juan Pablo Vielma H. Milton Stewart School of Industrial and Systems Engineering Georgia Institute of Technology Joint
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 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 informationInteger Linear Programming
Integer Linear Programming Solution : cutting planes and Branch and Bound Hugues Talbot Laboratoire CVN April 13, 2018 IP Resolution Gomory s cutting planes Solution branch-and-bound General method Resolution
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 informationInteger programming (part 2) Lecturer: Javier Peña Convex Optimization /36-725
Integer programming (part 2) Lecturer: Javier Peña Convex Optimization 10-725/36-725 Last time: integer programming Consider the problem min x subject to f(x) x C x j Z, j J where f : R n R, C R n are
More informationto work with) can be solved by solving their LP relaxations with the Simplex method I Cutting plane algorithms, e.g., Gomory s fractional cutting
Summary so far z =max{c T x : Ax apple b, x 2 Z n +} I Modeling with IP (and MIP, and BIP) problems I Formulation for a discrete set that is a feasible region of an IP I Alternative formulations for the
More informationTHE ART OF INTEGER PROGRAMMING - RELAXATION. D. E. Bell
THE ART OF INTEGER PROGRAMMING - RELAXATION D. E. Bell November 1973 WP-73-8 Working Papers are not intended for distribution outside of IIASA, and are solely for discussion and information purposes. The
More information21. Set cover and TSP
CS/ECE/ISyE 524 Introduction to Optimization Spring 2017 18 21. Set cover and TSP ˆ Set covering ˆ Cutting problems and column generation ˆ Traveling salesman problem Laurent Lessard (www.laurentlessard.com)
More informationIP Cut Homework from J and B Chapter 9: 14, 15, 16, 23, 24, You wish to solve the IP below with a cutting plane technique.
IP Cut Homework from J and B Chapter 9: 14, 15, 16, 23, 24, 31 14. You wish to solve the IP below with a cutting plane technique. Maximize 4x 1 + 2x 2 + x 3 subject to 14x 1 + 10x 2 + 11x 3 32 10x 1 +
More informationSolving LP and MIP Models with Piecewise Linear Objective Functions
Solving LP and MIP Models with Piecewise Linear Obective Functions Zonghao Gu Gurobi Optimization Inc. Columbus, July 23, 2014 Overview } Introduction } Piecewise linear (PWL) function Convex and convex
More informationInteger Linear Programming (ILP)
Integer Linear Programming (ILP) Zdeněk Hanzálek, Přemysl Šůcha hanzalek@fel.cvut.cz CTU in Prague March 8, 2017 Z. Hanzálek (CTU) Integer Linear Programming (ILP) March 8, 2017 1 / 43 Table of contents
More informationMVE165/MMG631 Linear and integer optimization with applications Lecture 8 Discrete optimization: theory and algorithms
MVE165/MMG631 Linear and integer optimization with applications Lecture 8 Discrete optimization: theory and algorithms Ann-Brith Strömberg 2017 04 07 Lecture 8 Linear and integer optimization with applications
More informationMVE165/MMG630, Applied Optimization Lecture 6 Integer linear programming: models and applications; complexity. Ann-Brith Strömberg
MVE165/MMG630, Integer linear programming: models and applications; complexity Ann-Brith Strömberg 2011 04 01 Modelling with integer variables (Ch. 13.1) Variables Linear programming (LP) uses continuous
More informationAnalyzing the computational impact of individual MINLP solver components
Analyzing the computational impact of individual MINLP solver components Ambros M. Gleixner joint work with Stefan Vigerske Zuse Institute Berlin MATHEON Berlin Mathematical School MINLP 2014, June 4,
More informationRecent Advances in Mixed Integer Programming Modeling and Computation
Recent Advances in Mixed Integer Programming Modeling and Computation Juan Pablo Vielma Massachusetts Institute of Technology MIT Center for Transportation & Logistics. Cambridge, Massachusetts, June,
More informationBranch-and-Bound. Leo Liberti. LIX, École Polytechnique, France. INF , Lecture p. 1
Branch-and-Bound Leo Liberti LIX, École Polytechnique, France INF431 2011, Lecture p. 1 Reminders INF431 2011, Lecture p. 2 Problems Decision problem: a question admitting a YES/NO answer Example HAMILTONIAN
More informationThe quest for finding Hamiltonian cycles
The quest for finding Hamiltonian cycles Giang Nguyen School of Mathematical Sciences University of Adelaide Travelling Salesman Problem Given a list of cities and distances between cities, what is the
More informationIntroduction to Integer Linear Programming
Lecture 7/12/2006 p. 1/30 Introduction to Integer Linear Programming Leo Liberti, Ruslan Sadykov LIX, École Polytechnique liberti@lix.polytechnique.fr sadykov@lix.polytechnique.fr Lecture 7/12/2006 p.
More informationOutline. Outline. Outline DMP204 SCHEDULING, TIMETABLING AND ROUTING. 1. Scheduling CPM/PERT Resource Constrained Project Scheduling Model
Outline DMP204 SCHEDULING, TIMETABLING AND ROUTING Lecture 3 and Mixed Integer Programg Marco Chiarandini 1. Resource Constrained Project Model 2. Mathematical Programg 2 Outline Outline 1. Resource Constrained
More informationResource Constrained Project Scheduling Linear and Integer Programming (1)
DM204, 2010 SCHEDULING, TIMETABLING AND ROUTING Lecture 3 Resource Constrained Project Linear and Integer Programming (1) Marco Chiarandini Department of Mathematics & Computer Science University of Southern
More informationDiscrete (and Continuous) Optimization WI4 131
Discrete (and Continuous) Optimization WI4 131 Kees Roos Technische Universiteit Delft Faculteit Electrotechniek, Wiskunde en Informatica Afdeling Informatie, Systemen en Algoritmiek e-mail: C.Roos@ewi.tudelft.nl
More informationIntroduction to Mathematical Programming IE406. Lecture 21. Dr. Ted Ralphs
Introduction to Mathematical Programming IE406 Lecture 21 Dr. Ted Ralphs IE406 Lecture 21 1 Reading for This Lecture Bertsimas Sections 10.2, 10.3, 11.1, 11.2 IE406 Lecture 21 2 Branch and Bound Branch
More informationScheduling and Optimization Course (MPRI)
MPRI Scheduling and optimization: lecture p. /6 Scheduling and Optimization Course (MPRI) Leo Liberti LIX, École Polytechnique, France MPRI Scheduling and optimization: lecture p. /6 Teachers Christoph
More informationModeling with Integer Programming
Modeling with Integer Programg Laura Galli December 18, 2014 We can use 0-1 (binary) variables for a variety of purposes, such as: Modeling yes/no decisions Enforcing disjunctions Enforcing logical conditions
More informationThe Surprisingly Complicated Case of [Convex] Quadratic Optimization
The Surprisingly Complicated Case of [Convex] Quadratic Optimization Robert Fourer 4er@ampl.com AMPL Optimization Inc. www.ampl.com +1 773-336-AMPL U.S.-Mexico Workshop on Optimization and Its Applications
More informationLarge-scale optimization and decomposition methods: outline. Column Generation and Cutting Plane methods: a unified view
Large-scale optimization and decomposition methods: outline I Solution approaches for large-scaled problems: I Delayed column generation I Cutting plane methods (delayed constraint generation) 7 I Problems
More informationInteger Programming ISE 418. Lecture 8. Dr. Ted Ralphs
Integer Programming ISE 418 Lecture 8 Dr. Ted Ralphs ISE 418 Lecture 8 1 Reading for This Lecture Wolsey Chapter 2 Nemhauser and Wolsey Sections II.3.1, II.3.6, II.4.1, II.4.2, II.5.4 Duality for Mixed-Integer
More informationMathematics for Decision Making: An Introduction. Lecture 8
Mathematics for Decision Making: An Introduction Lecture 8 Matthias Köppe UC Davis, Mathematics January 29, 2009 8 1 Shortest Paths and Feasible Potentials Feasible Potentials Suppose for all v V, there
More informationIndicator Constraints in Mixed-Integer Programming
Indicator Constraints in Mixed-Integer Programming Andrea Lodi University of Bologna, Italy - andrea.lodi@unibo.it Amaya Nogales-Gómez, Universidad de Sevilla, Spain Pietro Belotti, FICO, UK Matteo Fischetti,
More informationLogic, Optimization and Data Analytics
Logic, Optimization and Data Analytics John Hooker Carnegie Mellon University United Technologies Research Center, Cork, Ireland August 2015 Thesis Logic and optimization have an underlying unity. Ideas
More informationMIXED INTEGER PROGRAMMING APPROACHES FOR NONLINEAR AND STOCHASTIC PROGRAMMING
MIXED INTEGER PROGRAMMING APPROACHES FOR NONLINEAR AND STOCHASTIC PROGRAMMING A Thesis Presented to The Academic Faculty by Juan Pablo Vielma Centeno In Partial Fulfillment of the Requirements for the
More informationSplit Cuts for Convex Nonlinear Mixed Integer Programming
Split Cuts for Convex Nonlinear Mixed Integer Programming Juan Pablo Vielma Massachusetts Institute of Technology joint work with D. Dadush and S. S. Dey S. Modaresi and M. Kılınç Georgia Institute of
More informationA note on : A Superior Representation Method for Piecewise Linear Functions
A note on : A Superior Representation Method for Piecewise Linear Functions Juan Pablo Vielma Business Analytics and Mathematical Sciences Department, IBM T. J. Watson Research Center, Yorktown Heights,
More informationTopics in Complexity Theory
Topics in Complexity Theory Announcements Final exam this Friday from 12:15PM-3:15PM Please let us know immediately after lecture if you want to take the final at an alternate time and haven't yet told
More informationMINLP: Theory, Algorithms, Applications: Lecture 3, Basics of Algorothms
MINLP: Theory, Algorithms, Applications: Lecture 3, Basics of Algorothms Jeff Linderoth Industrial and Systems Engineering University of Wisconsin-Madison Jonas Schweiger Friedrich-Alexander-Universität
More informationHot-Starting NLP Solvers
Hot-Starting NLP Solvers Andreas Wächter Department of Industrial Engineering and Management Sciences Northwestern University waechter@iems.northwestern.edu 204 Mixed Integer Programming Workshop Ohio
More informationOutline. Relaxation. Outline DMP204 SCHEDULING, TIMETABLING AND ROUTING. 1. Lagrangian Relaxation. Lecture 12 Single Machine Models, Column Generation
Outline DMP204 SCHEDULING, TIMETABLING AND ROUTING 1. Lagrangian Relaxation Lecture 12 Single Machine Models, Column Generation 2. Dantzig-Wolfe Decomposition Dantzig-Wolfe Decomposition Delayed Column
More informationThe traveling salesman problem
Chapter 58 The traveling salesman problem The traveling salesman problem (TSP) asks for a shortest Hamiltonian circuit in a graph. It belongs to the most seductive problems in combinatorial optimization,
More informationis called an integer programming (IP) problem. model is called a mixed integer programming (MIP)
INTEGER PROGRAMMING Integer Programming g In many problems the decision variables must have integer values. Example: assign people, machines, and vehicles to activities in integer quantities. If this is
More informationBCOL RESEARCH REPORT 07.04
BCOL RESEARCH REPORT 07.04 Industrial Engineering & Operations Research University of California, Berkeley, CA 94720-1777 LIFTING FOR CONIC MIXED-INTEGER PROGRAMMING ALPER ATAMTÜRK AND VISHNU NARAYANAN
More informationDisconnecting Networks via Node Deletions
1 / 27 Disconnecting Networks via Node Deletions Exact Interdiction Models and Algorithms Siqian Shen 1 J. Cole Smith 2 R. Goli 2 1 IOE, University of Michigan 2 ISE, University of Florida 2012 INFORMS
More informationLecture 8: Column Generation
Lecture 8: Column Generation (3 units) Outline Cutting stock problem Classical IP formulation Set covering formulation Column generation A dual perspective Vehicle routing problem 1 / 33 Cutting stock
More informationNew Integer Programming Formulations of the Generalized Travelling Salesman Problem
American Journal of Applied Sciences 4 (11): 932-937, 2007 ISSN 1546-9239 2007 Science Publications New Integer Programming Formulations of the Generalized Travelling Salesman Problem Petrica C. Pop Department
More informationApproximation Algorithms
Approximation Algorithms Announcements Problem Set 9 due right now. Final exam this Monday, Hewlett 200 from 12:15PM- 3:15PM Please let us know immediately after lecture if you want to take the final at
More informationExtended Formulations, Lagrangian Relaxation, & Column Generation: tackling large scale applications
Extended Formulations, Lagrangian Relaxation, & Column Generation: tackling large scale applications François Vanderbeck University of Bordeaux INRIA Bordeaux-Sud-Ouest part : Defining Extended Formulations
More informationPartial Path Column Generation for the Vehicle Routing Problem with Time Windows
Partial Path Column Generation for the Vehicle Routing Problem with Time Windows Bjørn Petersen & Mads Kehlet Jepsen } DIKU Department of Computer Science, University of Copenhagen Universitetsparken 1,
More informationFrom structures to heuristics to global solvers
From structures to heuristics to global solvers Timo Berthold Zuse Institute Berlin DFG Research Center MATHEON Mathematics for key technologies OR2013, 04/Sep/13, Rotterdam Outline From structures to
More informationSection #2: Linear and Integer Programming
Section #2: Linear and Integer Programming Prof. Dr. Sven Seuken 8.3.2012 (with most slides borrowed from David Parkes) Housekeeping Game Theory homework submitted? HW-00 and HW-01 returned Feedback on
More information16.410/413 Principles of Autonomy and Decision Making
6.4/43 Principles of Autonomy and Decision Making Lecture 8: (Mixed-Integer) Linear Programming for Vehicle Routing and Motion Planning Emilio Frazzoli Aeronautics and Astronautics Massachusetts Institute
More informationA Fast Heuristic for GO and MINLP
A Fast Heuristic for GO and MINLP John W. Chinneck, M. Shafique, Systems and Computer Engineering Carleton University, Ottawa, Canada Introduction Goal: Find a good quality GO/MINLP solution quickly. Trade
More informationDiscrete Optimization 2010 Lecture 7 Introduction to Integer Programming
Discrete Optimization 2010 Lecture 7 Introduction to Integer Programming Marc Uetz University of Twente m.uetz@utwente.nl Lecture 8: sheet 1 / 32 Marc Uetz Discrete Optimization Outline 1 Intro: The Matching
More informationWhat is an integer program? Modelling with Integer Variables. Mixed Integer Program. Let us start with a linear program: max cx s.t.
Modelling with Integer Variables jesla@mandtudk Department of Management Engineering Technical University of Denmark What is an integer program? Let us start with a linear program: st Ax b x 0 where A
More informationCutting Plane Separators in SCIP
Cutting Plane Separators in SCIP Kati Wolter Zuse Institute Berlin DFG Research Center MATHEON Mathematics for key technologies 1 / 36 General Cutting Plane Method MIP min{c T x : x X MIP }, X MIP := {x
More informationMDD-based Postoptimality Analysis for Mixed-integer Programs
MDD-based Postoptimality Analysis for Mixed-integer Programs John Hooker, Ryo Kimura Carnegie Mellon University Thiago Serra Mitsubishi Electric Research Laboratories Symposium on Decision Diagrams for
More informationTeaching Integer Programming Formulations Using the Traveling Salesman Problem
SIAM REVIEW Vol. 45, No. 1, pp. 116 123 c 2003 Society for Industrial and Applied Mathematics Teaching Integer Programming Formulations Using the Traveling Salesman Problem Gábor Pataki Abstract. We designed
More informationNP-Complete Problems. More reductions
NP-Complete Problems More reductions Definitions P: problems that can be solved in polynomial time (typically in n, size of input) on a deterministic Turing machine Any normal computer simulates a DTM
More informationGestion de la production. Book: Linear Programming, Vasek Chvatal, McGill University, W.H. Freeman and Company, New York, USA
Gestion de la production Book: Linear Programming, Vasek Chvatal, McGill University, W.H. Freeman and Company, New York, USA 1 Contents 1 Integer Linear Programming 3 1.1 Definitions and notations......................................
More informationOptimization Exercise Set n.5 :
Optimization Exercise Set n.5 : Prepared by S. Coniglio translated by O. Jabali 2016/2017 1 5.1 Airport location In air transportation, usually there is not a direct connection between every pair of airports.
More informationUnit 1A: Computational Complexity
Unit 1A: Computational Complexity Course contents: Computational complexity NP-completeness Algorithmic Paradigms Readings Chapters 3, 4, and 5 Unit 1A 1 O: Upper Bounding Function Def: f(n)= O(g(n)) if
More informationIncremental Formulations for SOS1 Variables
Incremental Formulations for SOS Variables Juan Pablo Vielma Massachusetts Institute of Technology joint work with Sercan Yıldız arnegie Mellon University INFORMS nnual Meeting, October 22 Phoenix, rizona
More informationStructured Problems and Algorithms
Integer and quadratic optimization problems Dept. of Engg. and Comp. Sci., Univ. of Cal., Davis Aug. 13, 2010 Table of contents Outline 1 2 3 Benefits of Structured Problems Optimization problems may become
More informationMath Models of OR: Traveling Salesman Problem
Math Models of OR: Traveling Salesman Problem John E. Mitchell Department of Mathematical Sciences RPI, Troy, NY 12180 USA November 2018 Mitchell Traveling Salesman Problem 1 / 19 Outline 1 Examples 2
More informationColumn Generation. ORLAB - Operations Research Laboratory. Stefano Gualandi. June 14, Politecnico di Milano, Italy
ORLAB - Operations Research Laboratory Politecnico di Milano, Italy June 14, 2011 Cutting Stock Problem (from wikipedia) Imagine that you work in a paper mill and you have a number of rolls of paper of
More informationStochastic Decision Diagrams
Stochastic Decision Diagrams John Hooker CORS/INFORMS Montréal June 2015 Objective Relaxed decision diagrams provide an generalpurpose method for discrete optimization. When the problem has a dynamic programming
More informationNP-Completeness. Andreas Klappenecker. [based on slides by Prof. Welch]
NP-Completeness Andreas Klappenecker [based on slides by Prof. Welch] 1 Prelude: Informal Discussion (Incidentally, we will never get very formal in this course) 2 Polynomial Time Algorithms Most of the
More informationSEQUENTIAL AND SIMULTANEOUS LIFTING IN THE NODE PACKING POLYHEDRON JEFFREY WILLIAM PAVELKA. B.S., Kansas State University, 2011
SEQUENTIAL AND SIMULTANEOUS LIFTING IN THE NODE PACKING POLYHEDRON by JEFFREY WILLIAM PAVELKA B.S., Kansas State University, 2011 A THESIS Submitted in partial fulfillment of the requirements for the degree
More informationLow-Complexity Algorithms for Sequencing Jobs with a Fixed Number of Job-Classes
Low-Complexity Algorithms for Sequencing Jobs with a Fixed Number of Job-Classes Jack A.A. van der Veen Nijenrode University - The Netherlands School of Business Straatweg 25, 3621 BG Breukelen The Netherlands
More informationCHAPTER 3: INTEGER PROGRAMMING
CHAPTER 3: INTEGER PROGRAMMING Overview To this point, we have considered optimization problems with continuous design variables. That is, the design variables can take any value within a continuous feasible
More informationParallel PIPS-SBB Multi-level parallelism for 2-stage SMIPS. Lluís-Miquel Munguia, Geoffrey M. Oxberry, Deepak Rajan, Yuji Shinano
Parallel PIPS-SBB Multi-level parallelism for 2-stage SMIPS Lluís-Miquel Munguia, Geoffrey M. Oxberry, Deepak Rajan, Yuji Shinano ... Our contribution PIPS-PSBB*: Multi-level parallelism for Stochastic
More informationMixed Integer Programming Solvers: from Where to Where. Andrea Lodi University of Bologna, Italy
Mixed Integer Programming Solvers: from Where to Where Andrea Lodi University of Bologna, Italy andrea.lodi@unibo.it November 30, 2011 @ Explanatory Workshop on Locational Analysis, Sevilla A. Lodi, MIP
More informationInteger Programming for Bayesian Network Structure Learning
Integer Programming for Bayesian Network Structure Learning James Cussens Helsinki, 2013-04-09 James Cussens IP for BNs Helsinki, 2013-04-09 1 / 20 Linear programming The Belgian diet problem Fat Sugar
More informationSOLVING INTEGER LINEAR PROGRAMS. 1. Solving the LP relaxation. 2. How to deal with fractional solutions?
SOLVING INTEGER LINEAR PROGRAMS 1. Solving the LP relaxation. 2. How to deal with fractional solutions? Integer Linear Program: Example max x 1 2x 2 0.5x 3 0.2x 4 x 5 +0.6x 6 s.t. x 1 +2x 2 1 x 1 + x 2
More information15.081J/6.251J Introduction to Mathematical Programming. Lecture 24: Discrete Optimization
15.081J/6.251J Introduction to Mathematical Programming Lecture 24: Discrete Optimization 1 Outline Modeling with integer variables Slide 1 What is a good formulation? Theme: The Power of Formulations
More informationSensitivity Analysis for Intuitionistic Fuzzy Travelling Salesman Problem
Volume 119 No. 9 2018, 143-154 ISSN: 1311-8080 (printed version); ISSN: 1314-3395 (on-line version) url: http://www.ijpam.eu ijpam.eu Sensitivity Analysis for Intuitionistic Fuzzy Travelling Salesman Problem
More informationA note on : A Superior Representation Method for Piecewise Linear Functions by Li, Lu, Huang and Hu
A note on : A Superior Representation Method for Piecewise Linear Functions by Li, Lu, Huang and Hu Juan Pablo Vielma, Shabbir Ahmed and George Nemhauser H. Milton Stewart School of Industrial and Systems
More informationOperations Research Methods in Constraint Programming
Operations Research Methods in Constraint Programming John Hooker Carnegie Mellon University Prepared for Lloret de Mar, Spain June 2007 2007 Slide 1 CP and OR Have Complementary Strengths CP: Inference
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 informationSection Notes 9. Midterm 2 Review. Applied Math / Engineering Sciences 121. Week of December 3, 2018
Section Notes 9 Midterm 2 Review Applied Math / Engineering Sciences 121 Week of December 3, 2018 The following list of topics is an overview of the material that was covered in the lectures and sections
More information19. Fixed costs and variable bounds
CS/ECE/ISyE 524 Introduction to Optimization Spring 2017 18 19. Fixed costs and variable bounds ˆ Fixed cost example ˆ Logic and the Big M Method ˆ Example: facility location ˆ Variable lower bounds Laurent
More informationInteger Programming Chapter 15
Integer Programming Chapter 15 University of Chicago Booth School of Business Kipp Martin November 9, 2016 1 / 101 Outline Key Concepts Problem Formulation Quality Solver Options Epsilon Optimality Preprocessing
More informationAsymmetric Traveling Salesman Problem (ATSP): Models
Asymmetric Traveling Salesman Problem (ATSP): Models Given a DIRECTED GRAPH G = (V,A) with V = {,, n} verte set A = {(i, j) : i V, j V} arc set (complete digraph) c ij = cost associated with arc (i, j)
More informationCombinatorial optimization problems
Combinatorial optimization problems Heuristic Algorithms Giovanni Righini University of Milan Department of Computer Science (Crema) Optimization In general an optimization problem can be formulated as:
More informationApplied Integer Programming: Modeling and Solution
Applied Integer Programming: Modeling and Solution Chen, Batson, Dang Section 6. - 6.3 Blekinge Institute of Technology April 5, 05 Modeling Combinatorical Optimization Problems II Traveling Salesman Problem
More informationDecision Diagrams: Tutorial
Decision Diagrams: Tutorial John Hooker Carnegie Mellon University CP Summer School Cork, Ireland, June 2016 Decision Diagrams Used in computer science and AI for decades Logic circuit design Product configuration
More informationScheduling Home Hospice Care with Logic-Based Benders Decomposition
Scheduling Home Hospice Care with Logic-Based Benders Decomposition John Hooker Carnegie Mellon University Joint work with Aliza Heching Ryo Kimura Compassionate Care Hospice CMU Lehigh University October
More informationDecision Diagrams for Discrete Optimization
Decision Diagrams for Discrete Optimization Willem Jan van Hoeve Tepper School of Business Carnegie Mellon University www.andrew.cmu.edu/user/vanhoeve/mdd/ Acknowledgments: David Bergman, Andre Cire, Samid
More informationExercises NP-completeness
Exercises NP-completeness Exercise 1 Knapsack problem Consider the Knapsack problem. We have n items, each with weight a j (j = 1,..., n) and value c j (j = 1,..., n) and an integer B. All a j and c j
More informationInteger Programming Duality
Integer Programming Duality M. Guzelsoy T. K. Ralphs July, 2010 1 Introduction This article describes what is known about duality for integer programs. It is perhaps surprising that many of the results
More informationPart 4. Decomposition Algorithms
In the name of God Part 4. 4.4. Column Generation for the Constrained Shortest Path Problem Spring 2010 Instructor: Dr. Masoud Yaghini Constrained Shortest Path Problem Constrained Shortest Path Problem
More informationTractable & Intractable Problems
Tractable & Intractable Problems We will be looking at : What is a P and NP problem NP-Completeness The question of whether P=NP The Traveling Salesman problem again Programming and Data Structures 1 Polynomial
More informationOptimization in Process Systems Engineering
Optimization in Process Systems Engineering M.Sc. Jan Kronqvist Process Design & Systems Engineering Laboratory Faculty of Science and Engineering Åbo Akademi University Most optimization problems in production
More information