Primal Heuristics for Branch-and-Price Algorithms
|
|
- Georgiana Brooke Freeman
- 6 years ago
- Views:
Transcription
1 Primal Heuristics for Branch-and-Price Algorithms Marco Lübbecke and Christian Puchert Abstract In this paper, we present several primal heuristics which we implemented in the branch-and-price solver GCG based on the SCIP framework. This involves new heuristics as well as heuristics from the literature that make use of the reformulation yielded by the Dantzig-Wolfe decomposition. We give short descriptions of those heuristics and briefly discuss computational results. Furthermore, we give an outlook on current and further development. 1 Introduction A wealth of practical decision and planning problems in operations research can be modeled as mixed integer programs (MIPs) with a special structure; this structure has been recently successfully exploited in branch-and-price algorithms [6] which are based on Dantzig-Wolfe decomposition [5]; we assume the reader to be familiar with both concepts. Primal heuristics are a very important aspect of MIP solving: In practice, the decision maker is often satisfied with a good solution to her problem; besides, such solutions can help accelerate the solution process. We developed and tested heuristics specially tailored for branch-and-price that exploit the fact that two different formulations of the MIP are available via the Dantzig-Wolfe decomposition: the original MIP, and the extended reformulation. We implemented our heuristics in the branch-and-price solver GCG [8] which is based on the well-known SCIP framework [1]. Whereas today s branch-and-price algorithms are almost always tailored to the application and are thus reserved to experts, the goal of the GCG project is to provide a more generally accessible solver Marco Lübbecke RWTH Aachen University, Operations Research, Templergraben 64, Aachen, Germany marco.luebbecke@rwth-aachen.de Christian Puchert RWTH Aachen University, Operations Research, Templergraben 64, Aachen, Germany puchert@or.rwth-aachen.de 1
2 2 Marco Lübbecke and Christian Puchert that is available to a much broader range of users. Our work contributes towards this goal by providing primal heuristics for branch-and-price. In the remainder, we will present three types of such heuristics called extreme point, restricted master, and column selection heuristics and discuss computational results. This work is based on the master s thesis of the second author [11]. 2 Preliminaries We want to solve MIPs of the form min c T x s. t. Ax b x X, (1) with X = {x R n : Dx d,x i Z for all i I {1,...,n}}. Typically, Ax b are some complicated or linking constraints, whereas the set X often has a block structure X = X 1 X K with X k = {x k R n k : D k x k d k,x k i Z for all i I k }, k = 1,...,K, I k {1,...n k }, K k=1 n k = n. Via Dantzig-Wolfe decomposition, (1) is reformulated to a master problem min s. t. K k=1 K k=1 p P k c p λ kp a p λ kp b p P k λ kp = 1 p P k for all k λ k R P k + for all k, (2) where integrality of the x i is enforced by constraints p Pk λ kp x kp = x k for all k and xi k Z for all i I k (convexification) or λ k Z P k + for all k (discretization). The master problem (2) bears too many variables to explicitly work with and column generation is typically applied. The pricing subproblems to generate variables as needed read as min (c k ) T x k (π k ) T A k x k π 0 s. t. x k (3) X k, one for each block k.
3 Primal Heuristics for Branch-and-Price Algorithms 3 3 Extreme Points Crossover Each component x k of the solution x yielded by solving the master problem (2) is a convex combination x k = p Pk λkp x kp, with P k being the (index) set of extreme points obtained by solving the k-th pricing problem. These are integer feasible and satisfy the constraints D k x k d k, while x k satisfies all linear constraints but is usually not integer feasible. This leads to the question: Can we exploit the integer feasibility of the x kp to obtain a solution x satisfying all constraints? This is what extreme point based heuristics try to achieve. The extreme points we are particularly interested in are those with λ kp > 0. It may happen that there is a number of coordinates in which these extreme points are identical. Especially when their values are also shared by x, it seems worth exploring whether there are more integer feasible or feasible points taking the same values. Therefore, Extreme Points Crossover performs a neighborhood search on the extreme points by crossing them: For each block, it takes those x xp with the highest λ kp, fixes the variables in which the selected points are identical and solves a sub-mip. As its name indicates, the heuristic can be compared to the Crossover heuristic discussed in [4,12] and implemented in SCIP, which considers a number of already known feasible solution for crossing. The main advantage of our heuristic against this one is that it does not need any previously found feasible solution. Fig. 1 Extreme Points Crossover: An LP feasible solution x is a convex combination of x p 1,x p 2 conv(x). All the three points have a coordinate in common which is also shared by the encircled feasible point. x p1 x x p2 Subproblem Master problem 4 Restricted Master When we solve the master problem, we dynamically add variables that improve the current solution. As the total number of master variables may be exponentially high compared to the number of original variables, this means that the master problem may soon reach a size where searching for feasible solutions would become very time-consuming. This is what the Restricted Master heuristic (suggested by Joncour et al. [9]) tries to overcome. It searches for an integer feasible master solution by restricting the formulation again to a subset of promising variables, hence regarding a problem which is of considerably smaller size than the current master formulation.
4 4 Marco Lübbecke and Christian Puchert That is, we solve a sub-mip of the master formulation (without adding any new variables to it) where all variables not contained in this subset are fixed to zero. As our subset, we chose those master variables which are part of the master LP solution, i.e., the λ p with λ p > 0. 5 Column Selection Column Selection also suggested by Joncour et al. [9] is a constructive approach to build feasible solutions from scratch, with only the knowledge of extreme points previously generated by the pricing problems. In the discretization approach, integrality is enforced by integrality constraints on the master variables. This in particular means that choosing one x kp k, p k P k, for each block k and building the Cartesian product will lead to a point x = x 1p 1 x K p K that is feasible w.r.t. the integrality and the subproblem constraints; only the master constraints may be violated. The heuristic constructs a product as above by successively choosing x kp k previously generated by the subproblems, while trying to avoid master infeasibility. Technically, this happens by increasing the corresponding master variable λ kpk in the master problem by one. There are different approaches trying to achieve that, we implemented the following two: Greedy Column Selection: Starting with λ = 0, increase a master variable such that the number of currently violated master constraints is decreased as much as possible. As a tie breaker, one can optionally use the objective coefficients of the λ kpk variables; Relaxation Based Column Selection: Starting with a rounded down master LP solution λ, we increase master variables, preferring those which actually have been rounded (i.e., the λ kpk with λ kpk > 0); this method can also be viewed as a rounding heuristic on the master variables. 6 Results We ran GCG with the above mentioned heuristics on four types of combinatorial optimization problems: bin packing, vertex coloring, capacitated p-median and resource allocation. These problems bear a bordered block-diagonal structure or at least can be brought into one which can be easily exploited by GCG. Apart from that, we also tested how several of the default SCIP heuristics [4] such as rounding heuristics, RENS, diving heuristics, and improvement heuristics perform within the branch-and-price scenario (that is, when applied to the original x variables). Previously, no primal heuristics had been implemented in GCG; only the default SCIP heuristics on the master variables were used. Now, the new heuristics as well as the SCIP heuristics on the original variable space are available to GCG.
5 Primal Heuristics for Branch-and-Price Algorithms 5 The heuristics did unfortunately not so much improve the solution time of GCG. This is actually no big surprise since it already performed well on especially the bin packing and the smaller coloring instances before the heuristics had been implemented. However, solutions could be found and the number of branch-and-bound nodes could be reduced on a large number of the problems on which the heuristics were tested; in particular, they were successful on some instances where GCG previously failed. Restricted Master Column Selection Extreme Point BIN Solutions found [% inst.] Best solution [% inst.] Running time [sec, total] Running time [sec, geom. mean] COLEASY Solutions found [% inst.] Best solution [% inst.] Running time [sec, total] Running time [sec, geom. mean] COLHARD Solutions found [% inst.] Best solution [% inst.] Running time [sec, total] Running time [sec, geom. mean] CPMP Solutions found [% inst.] Best solution [% inst.] Running time [sec, total] Running time [sec, geom. mean] RAP Solutions found [% inst.] Best solution [% inst.] Running time [sec, total] Running time [sec, geom. mean] Table 1 Results: For each heuristic and each test set, we indicate in how many percent of the instances the heuristics found a solution, how often it found the best known solution and how much time it spent over all instances in total and on average, respectively. Besides, we also had test runs on some MIPLIB2003 [2] and MIPLIB2010 [10] instances; recently, there has been made progress in detecting structures in such general MIPs even if it is not apparently visible (see [3] for details). Our results there indicate that the heuristics are also of use in that general setting; in particular, we recently succeeded with applying diving heuristics to MIPLIB2003 instances. The GCG versions of Coefficient Diving and Fractional Diving which use the master LP instead of the original LP were able to reduce the computational effort for the instances 10teams, modglob, p2756, pp08acuts, and rout. 7 Conclusion & Outlook The heuristics presented here are able to generate feasible solutions for quite a number of problem instances. Yet, there still lies some potential in them. Currently, we are improving the Extreme Points Crossover heuristic and the GCG diving heuristics. In addition, we are developing a column generation based variant of the Feasibility Pump heuristic by Fischetti et al. [7]. While the original Feasibility Pump decom-
6 6 Marco Lübbecke and Christian Puchert poses the problem into linear and integer constraints, our heuristic distinguishes like the Dantzig-Wolfe decomposition between master and subproblem constraints. That is, by alternately solving modified master problems and subproblems, we try to construct two sequences of points, hopefully converging towards each other and yielding a feasible solution. Often, like in stochastic programs or other multi-stage decision problems, a matrix has a staircase structure, i.e., we have linking variables belonging to multiple blocks. The presence of such variables is undesirable since in order to apply a Dantzig-Wolfe decomposition, one needs to copy these variables and introduce additional master constraints. A way to overcome this issue is to heuristically fix them, which is what fifty-fifty heuristics do. Another possibility is letting diving heuristics branch on those variables, hoping that they will get fixed by doing so. The results show that our heuristics have now become an important part of GCG; a number of them are now included by default. This motivates us to further pursue our research in this important area of computational integer programming. References 1. Achterberg, T.: SCIP: Solving constraint integer programs. Mathematical Programming Computation 1(1), 1 41 (2009) Achterberg, T., Koch, T., Martin, A.: MIPLIB Operations Research Letters 34(4), (2006) Bergner, M., Caprara, A., Furini, F., Lübbecke, M., Malaguti, E., Traversi, E.: Partial convexification of general MIPs by Dantzig-Wolfe reformulation. In: Günlük, O., Woeginger, G.J. (eds.), Integer Programming and Combinatorial Optimization (IPCO 2011), LNCS, 6655, pp Springer, Berlin (2011) 4. Berthold, T.: Primal Heuristics for Mixed Integer Programs. Master s thesis. Technische Universität Berlin (2006) 5. Dantzig, G.B., Wolfe, P.: Decomposition Principle for Linear Programs. Operations Research 8(1), (1960) 6. Desrosiers, J., Lübbecke, M.: A Primer in Column Generation. In: Desaulniers, G., Desrosiers, J., Solomon, M.M. (eds.), Column Generation, pp Springer, US (2005) 7. Fischetti, M., Glover, F., Lodi, A.: The feasibility pump. Mathematical Programming 104, (2005) 8. Gamrath, G., Lübbecke, M.: Experiments with a Generic Dantzig-Wolfe Decomposition for Integer Programs. In: Festa, P. (ed.), Experimental Algorithms, LNCS, 6049, pp Springer, Berlin (2010) 9. Joncour, C., Michel, S., Sadykov, R., Sverdlov, D., and Vanderbeck, F.: Column Generation based Primal Heuristics. Electronic Notes in Discrete Mathematics 3, (2010). ISCO 2010 International Symposium on Combinatorial Optimization. 10. T. Koch, T. Achterberg, E. Andersen, O. Bastert, T. Berthold, R. E. Bixby, E. Danna, G. Gamrath, A. M. Gleixner, S. Heinz, A. Lodi, H. Mittelmann, T. Ralphs, D. Salvagnin, D. E. Steffy, and K. Wolter. MIPLIB Mathematical Programming Computation, 3(2): , Puchert, C.: Primal Heuristics for Branch-and-Price Algorithms. Master s thesis. Technische Universität Darmstadt (2011) 12. Rothberg, E.: An Evolutionary Algorithm for Polishing Mixed Integer Programming Solutions. ILOG Inc. (2005)
Automatic Decomposition and Branch-and-Price A Status Report
Automatic Decomposition and Branch-and-Price A Status Report Marco E. Lübbecke RWTH Aachen University, Operations Research, Kackertstraße 7, D-52072 Aachen, Germany marco.luebbecke@rwth-aachen.de Abstract.
More informationPrimal Heuristics for Mixed Integer Programs with a Staircase Structure
Primal Heuristics for Mixed Integer Programs with a Staircase Structure Marco E. Lübbecke and Christian Puchert Chair of Operations Research, RWTH Aachen University, Kackertstr. 7, 52072 Aachen, Germany
More informationExperiments with a Generic Dantzig-Wolfe Decomposition for Integer Programs
Experiments with a Generic Dantzig-Wolfe Decomposition for Integer Programs Gerald Gamrath 1 and Marco E. Lübbecke 2 1 Zuse Institute Berlin, Takustr. 7, 14195 Berlin, Germany gamrath@zib.de 2 Technische
More informationColumn Generation Based Primal Heuristics
Column Generation Based Primal Heuristics C. Joncour, S. Michel, R. Sadykov, D. Sverdlov, F. Vanderbeck University Bordeaux 1 & INRIA team RealOpt Outline 1 Context Generic Primal Heuristics The Branch-and-Price
More informationThe Heuristic (Dark) Side of MIP Solvers. Asja Derviskadic, EPFL Vit Prochazka, NHH Christoph Schaefer, EPFL
The Heuristic (Dark) Side of MIP Solvers Asja Derviskadic, EPFL Vit Prochazka, NHH Christoph Schaefer, EPFL 1 Table of content [Lodi], The Heuristic (Dark) Side of MIP Solvers, Hybrid Metaheuristics, 273-284,
More informationLearning when to use a decomposition
Learning when to use a decomposition Markus Kruber 1, Marco E. Lübbecke 1, and Axel Parmentier 2 1 Chair of Operations Research, RWTH Aachen University, Kackertstrasse 7, 52072 Aachen, Germany, {kruber,
More informationRestrict-and-relax search for 0-1 mixed-integer programs
EURO J Comput Optim (23) :2 28 DOI.7/s3675-3-7-y ORIGINAL PAPER Restrict-and-relax search for - mixed-integer programs Menal Guzelsoy George Nemhauser Martin Savelsbergh Received: 2 September 22 / Accepted:
More informationThe MIP-Solving-Framework SCIP
The MIP-Solving-Framework SCIP Timo Berthold Zuse Institut Berlin DFG Research Center MATHEON Mathematics for key technologies Berlin, 23.05.2007 What Is A MIP? Definition MIP The optimization problem
More informationHeuristics in MILP. Group 1 D. Assouline, N. Molyneaux, B. Morén. Supervisors: Michel Bierlaire, Andrea Lodi. Zinal 2017 Winter School
Heuristics in MILP Group 1 D. Assouline, N. Molyneaux, B. Morén Supervisors: Michel Bierlaire, Andrea Lodi Zinal 2017 Winter School 0 / 23 Primal heuristics Original paper: Fischetti, M. and Lodi, A. (2011).
More informationHeuristics in Commercial MIP Solvers Part I (Heuristics in IBM CPLEX)
Andrea Tramontani CPLEX Optimization, IBM CWI, Amsterdam, June 12, 2018 Heuristics in Commercial MIP Solvers Part I (Heuristics in IBM CPLEX) Agenda CPLEX Branch-and-Bound (B&B) Primal heuristics in CPLEX
More informationAlternating Criteria Search: A Parallel Large Neighborhood Search Algorithm for Mixed Integer Programs
Alternating Criteria Search: A Parallel Large Neighborhood Search Algorithm for Mixed Integer Programs Lluís-Miquel Munguía 1, Shabbir Ahmed 2, David A. Bader 1, George L. Nemhauser 2, and Yufen Shao 3
More informationColumn Generation based Primal Heuristics
Electronic Notes in Discrete Mathematics 36 (2010) 695 702 www.elsevier.com/locate/endm Column Generation based Primal Heuristics C. Joncour(1,3), S. Michel (2), R. Sadykov (3,1), D. Sverdlov (3), F. Vanderbeck
More informationComparisons of Commercial MIP Solvers and an Adaptive Memory (Tabu Search) Procedure for a Class of 0-1 Integer Programming Problems
Comparisons of Commercial MIP Solvers and an Adaptive Memory (Tabu Search) Procedure for a Class of 0-1 Integer Programming Problems Lars M. Hvattum The Norwegian University of Science and Technology Trondheim,
More informationRounding and Propagation Heuristics for Mixed Integer Programming
Konrad-Zuse-Zentrum für Informationstechnik Berlin Takustraße 7 D-9 Berlin-Dahlem Germany TOBIAS ACHTERBERG TIMO BERTHOLD GREGOR HENDEL Rounding and Propagation Heuristics for Mixed Integer Programming
More informationApplied Mixed Integer Programming: Beyond 'The Optimum'
Applied Mixed Integer Programming: Beyond 'The Optimum' 14 Nov 2016, Simons Institute, Berkeley Pawel Lichocki Operations Research Team, Google https://developers.google.com/optimization/ Applied Mixed
More informationPrimal Heuristics in SCIP
Primal Heuristics in SCIP Timo Berthold Zuse Institute Berlin DFG Research Center MATHEON Mathematics for key technologies Berlin, 10/11/2007 Outline 1 Introduction Basics Integration Into SCIP 2 Available
More informationThe Gurobi Optimizer. Bob Bixby
The Gurobi Optimizer Bob Bixby Outline Gurobi Introduction Company Products Benchmarks Gurobi Technology Rethinking MIP MIP as a bag of tricks 8-Jul-11 2010 Gurobi Optimization 2 Gurobi Optimization Incorporated
More informationPenalty Alternating Direction Methods for Mixed- Integer Optimization: A New View on Feasibility Pumps
Penalty Alternating Direction Methods for Mixed- Integer Optimization: A New View on Feasibility Pumps Björn Geißler, Antonio Morsi, Lars Schewe, Martin Schmidt FAU Erlangen-Nürnberg, Discrete Optimization
More informationComparisons of Commercial MIP Solvers and an Adaptive Memory (Tabu Search) Procedure for a Class of 0 1 Integer Programming Problems
Algorithmic Operations Research Vol.7 (2012) 13 20 Comparisons of Commercial MIP Solvers and an Adaptive Memory (Tabu Search) Procedure for a Class of 0 1 Integer Programming Problems Lars MagnusHvattum
More informationOutline. Column Generation: Cutting Stock A very applied method. Introduction to Column Generation. Given an LP problem
Column Generation: Cutting Stock A very applied method thst@man.dtu.dk Outline History The Simplex algorithm (re-visited) Column Generation as an extension of the Simplex algorithm A simple example! DTU-Management
More informationColumn Generation: Cutting Stock
Column Generation: Cutting Stock A very applied method thst@man.dtu.dk DTU-Management Technical University of Denmark 1 Outline History The Simplex algorithm (re-visited) Column Generation as an extension
More informationTwo-layer Network Design by Branch-and-Cut featuring MIP-based Heuristics
Two-layer Network Design by Branch-and-Cut featuring MIP-based Heuristics Sebastian Orlowski, Zuse Institute Berlin, Takustr. 7, D-14195 Berlin, orlowski@zib.de Arie M.C.A. Koster, Zuse Institute Berlin,
More informationPivot and Gomory Cut. A MIP Feasibility Heuristic NSERC
Pivot and Gomory Cut A MIP Feasibility Heuristic Shubhashis Ghosh Ryan Hayward shubhashis@randomknowledge.net hayward@cs.ualberta.ca NSERC CGGT 2007 Kyoto Jun 11-15 page 1 problem given a MIP, find a feasible
More informationColumn Generation Method for an Agent Scheduling Problem
Column Generation Method for an Agent Scheduling Problem Balázs Dezső Alpár Jüttner Péter Kovács Dept. of Algorithms and Their Applications, and Dept. of Operations Research Eötvös Loránd University, Budapest,
More informationRecursive column generation for the Tactical Berth Allocation Problem
Recursive column generation for the Tactical Berth Allocation Problem Ilaria Vacca 1 Matteo Salani 2 Michel Bierlaire 1 1 Transport and Mobility Laboratory, EPFL, Lausanne, Switzerland 2 IDSIA, Lugano,
More informationGENERAL ASSIGNMENT PROBLEM via Branch and Price JOHN AND LEI
GENERAL ASSIGNMENT PROBLEM via Branch and Price JOHN AND LEI Outline Review the column generation in Generalized Assignment Problem (GAP) GAP Examples in Branch and Price 2 Assignment Problem The assignment
More informationMotivation for Heuristics
MIP Heuristics 1 Motivation for Heuristics Why not wait for branching? Produce feasible solutions as quickly as possible Often satisfies user demands Avoid exploring unproductive sub trees Better reduced
More informationOutline. Modeling. Outline DMP204 SCHEDULING, TIMETABLING AND ROUTING. 1. Models Lecture 5 Mixed Integer Programming Models and Exercises
Outline DMP204 SCHEDULING, TIMETABLING AND ROUTING 1. Lecture 5 Mixed Integer Programming and Exercises Marco Chiarandini 2. 3. 2 Outline Modeling 1. Min cost flow Shortest path 2. Max flow Assignment
More informationCloud Branching MIP workshop, Ohio State University, 23/Jul/2014
Cloud Branching MIP workshop, Ohio State University, 23/Jul/2014 Timo Berthold Xpress Optimization Team Gerald Gamrath Zuse Institute Berlin Domenico Salvagnin Universita degli Studi di Padova This presentation
More informationExperiments On General Disjunctions
Experiments On General Disjunctions Some Dumb Ideas We Tried That Didn t Work* and Others We Haven t Tried Yet *But that may provide some insight Ted Ralphs, Serdar Yildiz COR@L Lab, Department of Industrial
More informationSelected Topics in Column Generation
Selected Topics in Column Generation February 1, 2007 Choosing a solver for the Master Solve in the dual space(kelly s method) by applying a cutting plane algorithm In the bundle method(lemarechal), a
More informationWireless frequency auctions: Mixed Integer Programs and Dantzig-Wolfe decomposition
Wireless frequency auctions: Mixed Integer Programs and Dantzig-Wolfe decomposition Laszlo Ladanyi (IBM T.J. Watson Research Center) joint work with Marta Eso (The Hotchkiss School) David Jensen (IBM T.J.
More informationExact solutions to mixed-integer linear programming problems
Exact solutions to mixed-integer linear programming problems Dan Steffy Zuse Institute Berlin and Oakland University Joint work with Bill Cook, Thorsten Koch and Kati Wolter November 18, 2011 Mixed-Integer
More informationColumn Generation based Primal Heuristics
Column Generation based Primal Heuristics Cédric Joncour(1), Sophie Michel (2), Ruslan Sadykov (3,1), Dimitri Sverdlov (3,1), François Vanderbeck (1,3) (1) Université Bordeaux 1, Institut de Mathématiques
More informationModern Benders (in a nutshell)
Modern Benders (in a nutshell) Matteo Fischetti, University of Padova (based on joint work with Ivana Ljubic and Markus Sinnl) Lunteren Conference on the Mathematics of Operations Research, January 17,
More informationParallel Branch & Bound
Parallel Branch & Bound Bernard Gendron Université de Montréal gendron@iro.umontreal.ca Outline Mixed integer programming (MIP) and branch & bound (B&B) Linear programming (LP) based B&B Relaxation and
More informationAccelerating Column Generation for a Survivable Network Design Problem
Accelerating Column Generation for a Survivable Network Design Problem Markus Leitner, Günther R. Raidl Ulrich Pferschy Department of Telematics and Network Engineering, Carinthia University of Applied
More information5.3 Cutting plane methods and Gomory fractional cuts
5.3 Cutting plane methods and Gomory fractional cuts (ILP) min c T x s.t. Ax b x 0integer feasible region X Assumption: a ij, c j and b i integer. Observation: The feasible region of an ILP can be described
More informationExploiting Degeneracy in MIP
Exploiting Degeneracy in MIP Tobias Achterberg 9 January 2018 Aussois Performance Impact in Gurobi 7.5+ 35% 32.0% 30% 25% 20% 15% 14.6% 10% 5.7% 7.9% 6.6% 5% 0% 2.9% 1.2% 0.1% 2.6% 2.6% Time limit: 10000
More informationSolving lexicographic multiobjective MIPs with Branch-Cut-Price
Solving lexicographic multiobjective MIPs with Branch-Cut-Price Marta Eso (The Hotchkiss School) Laszlo Ladanyi (IBM T.J. Watson Research Center) David Jensen (IBM T.J. Watson Research Center) McMaster
More information15.082J and 6.855J. Lagrangian Relaxation 2 Algorithms Application to LPs
15.082J and 6.855J Lagrangian Relaxation 2 Algorithms Application to LPs 1 The Constrained Shortest Path Problem (1,10) 2 (1,1) 4 (2,3) (1,7) 1 (10,3) (1,2) (10,1) (5,7) 3 (12,3) 5 (2,2) 6 Find the shortest
More informationRENS. The optimal rounding. Timo Berthold
Math. Prog. Comp. (2014) 6:33 54 DOI 10.1007/s12532-013-0060-9 FULL LENGTH PAPER RENS The optimal rounding Timo Berthold Received: 25 April 2012 / Accepted: 2 October 2013 / Published online: 1 November
More informationAdvanced Use of GAMS Solver Links
Advanced Use of GAMS Solver Links Michael Bussieck, Steven Dirkse, Stefan Vigerske GAMS Development 8th January 2013, ICS Conference, Santa Fe Standard GAMS solve Solve william minimizing cost using mip;
More informationGenerating Hard Instances for Robust Combinatorial Optimization
Generating Hard Instances for Robust Combinatorial Optimization Marc Goerigk 1 and Stephen J. Maher 2 1 Network and Data Science Management, University of Siegen, Germany 2 Department of Management Science,
More informationReceived: 27 October 2008 / Accepted: 1 September 2009 / Published online: 17 September 2009 Springer and Mathematical Programming Society 2009
Math. Prog. Comp. (2009) 1:201 222 DOI 10.1007/s12532-009-0007-3 FULL LENGTH PAPER Feasibility pump 2.0 Matteo Fischetti Domenico Salvagnin Received: 27 October 2008 / Accepted: 1 September 2009 / Published
More informationPart 4. Decomposition Algorithms Dantzig-Wolf Decomposition Algorithm
In the name of God Part 4. 4.1. Dantzig-Wolf Decomposition Algorithm Spring 2010 Instructor: Dr. Masoud Yaghini Introduction Introduction Real world linear programs having thousands of rows and columns.
More informationBranching rules revisited
Operations Research Letters 33 (2005) 42 54 Operations Research Letters www.elsevier.com/locate/dsw Branching rules revisited Tobias Achterberg a;, Thorsten Koch a, Alexander Martin b a Konrad-Zuse-Zentrum
More informationA Feasibility Pump heuristic for general Mixed-Integer Problems
A Feasibility Pump heuristic for general Mixed-Integer Problems Livio Bertacco, Matteo Fischetti, Andrea Lodi Department of Pure & Applied Mathematics, University of Padova, via Belzoni 7-35131 Padova
More informationCOLUMN GENERATION IN LINEAR PROGRAMMING
COLUMN GENERATION IN LINEAR PROGRAMMING EXAMPLE: THE CUTTING STOCK PROBLEM A certain material (e.g. lumber) is stocked in lengths of 9, 4, and 6 feet, with respective costs of $5, $9, and $. An order for
More informationThe Gurobi Solver V1.0
The Gurobi Solver V1.0 Robert E. Bixby Gurobi Optimization & Rice University Ed Rothberg, Zonghao Gu Gurobi Optimization 1 1 Oct 09 Overview Background Rethinking the MIP solver Introduction Tree of Trees
More informationCrew Scheduling Problem: A Column Generation Approach Improved by a Genetic Algorithm. Santos and Mateus (2007)
In the name of God Crew Scheduling Problem: A Column Generation Approach Improved by a Genetic Algorithm Spring 2009 Instructor: Dr. Masoud Yaghini Outlines Problem Definition Modeling As A Set Partitioning
More informationInteger Programming ISE 418. Lecture 7. Dr. Ted Ralphs
Integer Programming ISE 418 Lecture 7 Dr. Ted Ralphs ISE 418 Lecture 7 1 Reading for This Lecture Nemhauser and Wolsey Sections II.3.1, II.3.6, II.4.1, II.4.2, II.5.4 Wolsey Chapter 7 CCZ Chapter 1 Constraint
More informationIntroduction to Mathematical Programming IE406. Lecture 20. Dr. Ted Ralphs
Introduction to Mathematical Programming IE406 Lecture 20 Dr. Ted Ralphs IE406 Lecture 20 1 Reading for This Lecture Bertsimas Sections 10.1, 11.4 IE406 Lecture 20 2 Integer Linear Programming An integer
More informationA Generalized Wedelin Heuristic for Integer Programming
Online Supplement to A Generalized Wedelin Heuristic for Integer Programming Oliver Bastert Fair Isaac, Leam House, 64 Trinity Street, Leamington Spa, Warwicks CV32 5YN, UK, OliverBastert@fairisaac.com
More informationHandling first-order linear constraints with SCIP
Handling first-order linear constraints with SCIP James Cussens, University of York KU Leuven, 2015-02-16 James Cussens, University of York FO with SCIP KU Leuven, 2015-02-16 1 / 18 MIP Mixed integer programs
More informationAn Extension of the Multicut L-Shaped Method. INEN Large-Scale Stochastic Optimization Semester project. Svyatoslav Trukhanov
An Extension of the Multicut L-Shaped Method INEN 698 - Large-Scale Stochastic Optimization Semester project Svyatoslav Trukhanov December 13, 2005 1 Contents 1 Introduction and Literature Review 3 2 Formal
More informationComputational Integer Programming. Lecture 12: Branch and Cut. Dr. Ted Ralphs
Computational Integer Programming Lecture 12: Branch and Cut Dr. Ted Ralphs Computational MILP Lecture 12 1 Reading for This Lecture Wolsey Section 9.6 Nemhauser and Wolsey Section II.6 Martin Computational
More informationMVE165/MMG630, Applied Optimization Lecture 8 Integer linear programming algorithms. Ann-Brith Strömberg
MVE165/MMG630, Integer linear programming algorithms Ann-Brith Strömberg 2009 04 15 Methods for ILP: Overview (Ch. 14.1) Enumeration Implicit enumeration: Branch and bound Relaxations Decomposition methods:
More informationAlgorithms for Integer Programming
Algorithms for Integer Programming Laura Galli November 9, 2016 Unlike linear programming problems, integer programming problems are very difficult to solve. In fact, no efficient general algorithm is
More informationBranch-and-cut implementation of Benders decomposition Matteo Fischetti, University of Padova
Branch-and-cut implementation of Benders decomposition Matteo Fischetti, University of Padova 8th Cargese-Porquerolles Workshop on Combinatorial Optimization, August 2017 1 Mixed-Integer Programming We
More informationThe SCIP Optimization Suite 4.0
Zuse Institute Berlin Takustrasse 7 D-14195 Berlin-Dahlem Germany STEPHEN J. MAHER, TOBIAS FISCHER, TRISTAN GALLY, GERALD GAMRATH, AMBROS GLEIXNER, ROBERT LION GOTTWALD, GREGOR HENDEL, THORSTEN KOCH, MARCO
More informationAlgorithms for Decision Support. Integer linear programming models
Algorithms for Decision Support Integer linear programming models 1 People with reduced mobility (PRM) require assistance when travelling through the airport http://www.schiphol.nl/travellers/atschiphol/informationforpassengerswithreducedmobility.htm
More informationExact Algorithms for Mixed-Integer Bilevel Linear Programming
Exact Algorithms for Mixed-Integer Bilevel Linear Programming Matteo Fischetti, University of Padova (based on joint work with I. Ljubic, M. Monaci, and M. Sinnl) Lunteren Conference on the Mathematics
More informationBenders Decomposition
Benders Decomposition Using projections to solve problems thst@man.dtu.dk DTU-Management Technical University of Denmark 1 Outline Introduction Using projections Benders decomposition Simple plant location
More informationImplementing a B&C algorithm for Mixed-Integer Bilevel Linear Programming
Implementing a B&C algorithm for Mixed-Integer Bilevel Linear Programming Matteo Fischetti, University of Padova 8th Cargese-Porquerolles Workshop on Combinatorial Optimization, August 2017 1 Bilevel Optimization
More informationThe Three Phases of MIP Solving
Zuse Institute Berlin Takustrasse 7 D-14195 Berlin-Dahlem Germany TIMO BERTHOLD, GREGOR HENDEL, AND THORSTEN KOCH The Three Phases of MIP Solving The work for this article has been conducted within the
More informationMVE165/MMG631 Linear and integer optimization with applications Lecture 9 Discrete optimization: theory and algorithms
MVE165/MMG631 Linear and integer optimization with applications Lecture 9 Discrete optimization: theory and algorithms Ann-Brith Strömberg 2018 04 24 Lecture 9 Linear and integer optimization with applications
More informationBranch-Cut-and-Price solver for Vehicle Routing Problems
1 / 28 Branch-Cut-and-Price solver for Vehicle Routing Problems Ruslan Sadykov 1,2 Issam Tahiri 1,2 François Vanderbeck 2,1 Remi Duclos 1 Artur Pessoa 3 Eduardo Uchoa 3 1 Inria Bordeaux, France 2 Université
More informationA COMPUTATIONAL STUDY OF THE CUTTING PLANE TREE ALGORITHM FOR GENERAL MIXED-INTEGER LINEAR PROGRAMS
A COMPUTATIONAL STUDY OF THE CUTTING PLANE TREE ALGORITHM FOR GENERAL MIXED-INTEGER LINEAR PROGRAMS BINYUAN CHEN, SİMGE KÜÇÜKYAVUZ, AND SUVRAJEET SEN Abstract. The cutting plane tree (CPT) algorithm provides
More informationDecomposition in Integer Linear Programming
Decomposition in Integer Linear Programming T.K. Ralphs M.V. Galati December 3, 00 Abstract Both cutting plane methods and traditional decomposition methods are procedures that compute a bound on the optimal
More informationSolving a Challenging Quadratic 3D Assignment Problem
Solving a Challenging Quadratic 3D Assignment Problem Hans Mittelmann Arizona State University Domenico Salvagnin DEI - University of Padova Quadratic 3D Assignment Problem Quadratic 3D Assignment Problem
More informationBenders in a nutshell Matteo Fischetti, University of Padova
Benders in a nutshell Matteo Fischetti, University of Padova ODS 2017, Sorrento, September 2017 1 Benders decomposition The original Benders decomposition from the 1960s uses two distinct ingredients for
More informationSurrogate Gradient Algorithm for Lagrangian Relaxation 1,2
Surrogate Gradient Algorithm for Lagrangian Relaxation 1,2 X. Zhao 3, P. B. Luh 4, and J. Wang 5 Communicated by W.B. Gong and D. D. Yao 1 This paper is dedicated to Professor Yu-Chi Ho for his 65th birthday.
More informationIntroduction to Mathematical Programming IE496. Final Review. Dr. Ted Ralphs
Introduction to Mathematical Programming IE496 Final Review Dr. Ted Ralphs IE496 Final Review 1 Course Wrap-up: Chapter 2 In the introduction, we discussed the general framework of mathematical modeling
More information56:272 Integer Programming & Network Flows Final Examination -- December 14, 1998
56:272 Integer Programming & Network Flows Final Examination -- December 14, 1998 Part A: Answer any four of the five problems. (15 points each) 1. Transportation problem 2. Integer LP Model Formulation
More informationImproving the heuristic performance of Benders decomposition
Improving the heuristic performance of Benders decomposition Stephen J. Maher Department of Management Science, Lancaster University, Bailrigg, Lancaster LA1 4YX, UK Abstract A general enhancement of the
More informationNOTATION AND TERMINOLOGY
15.053x, Optimization Methods in Business Analytics Fall, 2016 October 4, 2016 A glossary of notation and terms used in 15.053x Weeks 1, 2, 3, 4 and 5. (The most recent week's terms are in blue). NOTATION
More informationTABU SEARCH FOR MIXED INTEGER PROGRAMMING
Chapter 11 TABU SEARCH FOR MIXED INTEGER PROGRAMMING João Pedro Pedroso Departamento de Ciência de Computadores, Faculdade de Ciências da Universidade do Porto, Rua do Campo Alegre, 823 4150-180 Porto,
More informationUsing Multiple Machines to Solve Models Faster with Gurobi 6.0
Using Multiple Machines to Solve Models Faster with Gurobi 6.0 Distributed Algorithms in Gurobi 6.0 Gurobi 6.0 includes 3 distributed algorithms Distributed concurrent LP (new in 6.0) MIP Distributed MIP
More informationConstruction of Minimum-Weight Spanners Mikkel Sigurd Martin Zachariasen
Construction of Minimum-Weight Spanners Mikkel Sigurd Martin Zachariasen University of Copenhagen Outline Motivation and Background Minimum-Weight Spanner Problem Greedy Spanner Algorithm Exact Algorithm:
More informationA Generic Separation Algorithm and Its Application to the Vehicle Routing Problem
A Generic Separation Algorithm and Its Application to the Vehicle Routing Problem Presented by: Ted Ralphs Joint work with: Leo Kopman Les Trotter Bill Pulleyblank 1 Outline of Talk Introduction Description
More informationMixed Integer Programming Class Library (MIPCL)
Mixed Integer Programming Class Library (MIPCL) Nicolai N. Pisaruk Belarus State University, Faculty of Economy, Nezavisimosty Av., 4, 220088 Minsk, Belarus April 20, 2016 Abstract The Mixed Integer Programming
More informationAn Integer Programming Approach to Packing Lightpaths on WDM Networks 파장분할다중화망의광경로패킹에대한정수계획해법. 1. Introduction
Journal of the Korean Institute of Industrial Engineers Vol. 32, No. 3, pp. 219-225, September 2006. An Integer Programming Approach to Packing Lightpaths on WDM Networks Kyungsik Lee 1 Taehan Lee 2 Sungsoo
More informationUniversità degli Studi di Milano. Branch-and-price algorithms for partitioning problems
Università degli Studi di Milano Facoltà di Scienze Matematiche, Fisiche e Naturali Dipartimento di Tecnologie dell Informazione Dottorato di Ricerca in Informatica, XVIII Ciclo Settore scientifico/disciplinare
More informationDecomposition in Integer Linear Programming
Decomposition in Integer Linear Programming M.V. Galati Advanced Analytics - Operations R & D, SAS Institute, Chesterbrook, PA 9 Ted K. Ralphs Department of Industrial and System Engineering, Lehigh University,
More informationAlgorithms II MIP Details
Algorithms II MIP Details What s Inside Gurobi Optimizer Algorithms for continuous optimization Algorithms for discrete optimization Automatic presolve for both LP and MIP Algorithms to analyze infeasible
More informationHow to use your favorite MIP Solver: modeling, solving, cannibalizing. Andrea Lodi University of Bologna, Italy
How to use your favorite MIP Solver: modeling, solving, cannibalizing Andrea Lodi University of Bologna, Italy andrea.lodi@unibo.it January-February, 2012 @ Universität Wien A. Lodi, How to use your favorite
More informationUsing the primal-dual interior point algorithm within the branch-price-and-cut method
Using the primal-dual interior point algorithm within the branch-price-and-cut method Pedro Munari a,1, Jacek Gondzio b a Instituto de Ciências Matemáticas e de Computação, University of São Paulo, Av.
More informationNew developments in the primal-dual column generation technique
New developments in the primal-dual column generation technique Jacek Gondzio a, Pablo González-Brevis a,1,, Pedro Munari b,2 a School of Mathematics, The University of Edinburgh, James Clerk Maxwell Building,
More informationPrimal Heuristics for Branch-and-Price: the assets of diving methods
Primal Heuristics for Branch-and-Price: the assets of diving methods Ruslan Sadykov, François Vanderbeck, Artur Pessoa, Issam Tahiri, Eduardo Uchoa To cite this version: Ruslan Sadykov, François Vanderbeck,
More information3 No-Wait Job Shops with Variable Processing Times
3 No-Wait Job Shops with Variable Processing Times In this chapter we assume that, on top of the classical no-wait job shop setting, we are given a set of processing times for each operation. We may select
More informationIn this paper we address the problem of finding a feasible solution of a generic MIP problem of the form
Mathematical Programming manuscript No. (will be inserted by the editor) Matteo Fischetti Fred Glover Andrea Lodi The feasibility pump Revised version Abstract. In this paper we consider the NP-hard problem
More informationThe feasibility pump. Matteo Fischetti Fred Glover Andrea Lodi. 1. Introduction
Math. Program., Ser. A 104, 91 104 (2005) Digital Object Identifier (DOI) 10.1007/s10107-004-0570-3 Matteo Fischetti Fred Glover Andrea Lodi The feasibility pump Received: May 17, 2004 / Accepted: November
More informationImproved Gomory Cuts for Primal Cutting Plane Algorithms
Improved Gomory Cuts for Primal Cutting Plane Algorithms S. Dey J-P. Richard Industrial Engineering Purdue University INFORMS, 2005 Outline 1 Motivation The Basic Idea Set up the Lifting Problem How to
More informationInteger Programming Theory
Integer Programming Theory Laura Galli October 24, 2016 In the following we assume all functions are linear, hence we often drop the term linear. In discrete optimization, we seek to find a solution x
More informationParallel and Distributed Optimization with Gurobi Optimizer
Parallel and Distributed Optimization with Gurobi Optimizer Our Presenter Dr. Tobias Achterberg Developer, Gurobi Optimization 2 Parallel & Distributed Optimization 3 Terminology for this presentation
More informationConflict Analysis in Mixed Integer Programming
Konrad-Zuse-Zentrum für Informationstechnik Berlin Takustraße 7 D-14195 Berlin-Dahlem Germany TOBIAS ACHTERBERG Conflict Analysis in Mixed Integer Programming URL: http://www.zib.de/projects/integer-optimization/mip
More informationIntroduction to Stochastic Combinatorial Optimization
Introduction to Stochastic Combinatorial Optimization Stefanie Kosuch PostDok at TCSLab www.kosuch.eu/stefanie/ Guest Lecture at the CUGS PhD course Heuristic Algorithms for Combinatorial Optimization
More informationHybrid Constraint Programming and Metaheuristic methods for Large Scale Optimization Problems
Hybrid Constraint Programming and Metaheuristic methods for Large Scale Optimization Problems Fabio Parisini Tutor: Paola Mello Co-tutor: Michela Milano Final seminars of the XXIII cycle of the doctorate
More informationGraph Coloring via Constraint Programming-based Column Generation
Graph Coloring via Constraint Programming-based Column Generation Stefano Gualandi Federico Malucelli Dipartimento di Elettronica e Informatica, Politecnico di Milano Viale Ponzio 24/A, 20133, Milan, Italy
More information