Heuristics in MILP. Group 1 D. Assouline, N. Molyneaux, B. Morén. Supervisors: Michel Bierlaire, Andrea Lodi. Zinal 2017 Winter School

Save this PDF as:
 WORD  PNG  TXT  JPG

Size: px
Start display at page:

Download "Heuristics in MILP. Group 1 D. Assouline, N. Molyneaux, B. Morén. Supervisors: Michel Bierlaire, Andrea Lodi. Zinal 2017 Winter School"

Transcription

1 Heuristics in MILP Group 1 D. Assouline, N. Molyneaux, B. Morén Supervisors: Michel Bierlaire, Andrea Lodi Zinal 2017 Winter School 0 / 23

2 Primal heuristics Original paper: Fischetti, M. and Lodi, A. (2011). Heuristics in mixed integer programming. Wiley Encyclopedia of Operations Research and Management Science LP-based heuristics Feasibility pump Pivoting methods Rounding and diving OCTANE MILP-based heuristics Local branching RINS Evolutionary algorithms 1 / 23

3 Pivoting methods 1 Pivoting methods 2 Relaxation Induced Neighborhood Search (RINS) 3 An evolutionary algorithm within MILP 4 Conclusion 2 / 23

4 Pivoting methods Pivoting Methods Balas, E. and Martin, C. H. (1980). Pivot and complement a heuristic for 0-1 programming. Management Science MILP with binary and continuous variables, no integer variables MILP formulation min x c T x s. t. Ax b x j binary, j B x j continuous, j C 3 / 23

5 Pivoting methods Pivoting Methods Given a solution to the LP-relaxation Non-basic variables: at lower or upper bound, x i = l i or x i = u i Basic variables: between lower and upper bound, l i < x i < u i (nondegeneracy) Binary variables Either on lower or on upper bound If all binary variables are non-basic, solution is feasible Pivoting methods tries to pivot binary variables out of basis 4 / 23

6 Pivoting methods Pivoting Methods - Pivot steps Pivot types 1 Pivot binary basic variable out of basis, maintaining primal feasibility. 2 Reduce the binary infeasibility i B min{x i x i, x i x i } 3 Pivot binary basic variable out of basis, allowing primal infeasibility. 5 / 23

7 Pivoting methods Pivoting Methods - Extension Pivot and shift - integer variables and improving solution Combined with local branching Combined with tabu search Pivot, cut, and dive Pivot and Gomory cut 6 / 23

8 Relaxation Induced Neighborhood Search (RINS) 1 Pivoting methods 2 Relaxation Induced Neighborhood Search (RINS) 3 An evolutionary algorithm within MILP 4 Conclusion 7 / 23

9 Relaxation Induced Neighborhood Search (RINS) Relaxation Induced Neighborhood Search (RINS) Part of the MILP-based heuristics. Improvement heuristic: start with reference feasible solution x of an MILP, and try to find a better one. Main idea: Uses the local neighborhood search principle, using information from the LP relaxation solution to build the neighborhood. Observation: often the incumbent solution x of an MILP and the optimal LP-relaxation x have a few of variables set to the same values. These values form a partial solution of good objective value. Neighborhood exploration is then formulated as a sub-milp. 8 / 23

10 Relaxation Induced Neighborhood Search (RINS) RINS algorithm At each node of the branch-and-cut tree, we have the incumbent x, and x from the LP-relaxation at the node, and do the following: 1 Fix the integer variables that have the same value in x and x ; 2 Set an objective cutoff based on the objective value of the current incumbent; 3 Solve the associated sub-milp on the remaining variables (a simplified and smaller in size MILP!); 4 If incumbent solution is improved, x is updated in the rest of the tree. Sub-MILP in RINS min x c T x s. t. Ax b c T x (1 ɛ) c T x for some ɛ > 0 x j = x j j I, with x j = x j x j integer, j I 9 / 23

11 Relaxation Induced Neighborhood Search (RINS) RINS in practice Each sub-milp might be large and difficult to solve: the neighborhood exploration if often truncated, by setting a limited number of nodes nl. RINS can be invoked at every node of the branch-and-cut tree: very suitable for integration within solvers. (The LP relaxation is potentially different at each node, so neighborhood diversification is well performed) Yet, neighborhoods from consecutive nodes might be similar: we don t invoke RINS at each node, but every f node (f 1). Another way to control the frequence of RINS use is to impose a minimum ratio of variables mfr to fix allowing to invoke RINS. It is embedded in CPLEX 12.6! 10 / 23

12 Relaxation Induced Neighborhood Search (RINS) RINS results and discussion (a) Small spread Figure: Taken from [Danna et al., 2005] 1 (b) Large spread 1 Danna, E., Rothberg, E., and Le Pape, C. (2005). Exploring relaxation induced neighborhoods to improve mip solutions. Mathematical Programming, 102(1): / 23

13 Relaxation Induced Neighborhood Search (RINS) RINS discussion with LB The main difference between Local Branching and RINS is that if the former uses soft-fixing (with addition of a linear inequality), the latter uses hard-fixing (by considering the LP relaxation into account). That constitutes the biggest advantage of RINS: the sub-milp are smaller. Also, RINS can be invoked more often since x changes between the nodes. RINS theoretically leads to faster improvements, but lacks the flexibility of LB / 23

14 Relaxation Induced Neighborhood Search (RINS) RINS variants DINS: Distance Induced Neighborhood Search Idea: the most promising solutions are those close to x. The neighborhood is defined by the distance inequality: j I xj xj j I xj xj If xj xj <0.5, xj = x j, otherwise we set the bounds of the variables so that the distance with x cannot increase. RENS: Relaxation Enforced Neighborhood Search Idea: define a reduced MILP as the set of all integer solutions which can be obtained by rounding a relaxed solution x. It is basically a sophisticated rounding heuristic, and does not require an incumbent solution to work with. 13 / 23

15 An evolutionary algorithm within MILP 1 Pivoting methods 2 Relaxation Induced Neighborhood Search (RINS) 3 An evolutionary algorithm within MILP 4 Conclusion 14 / 23

16 An evolutionary algorithm within MILP Evolutionary algorithm within MILP Evolutionary algorithms are the parents of genetic algorithms (GA) The four important steps of GA are found in this heuristic: Population: a pool of feasible solutions (hopefully good ones) Combination: take two (or more) solutions and merge into a child Mutation: randomly change part of the solution Selection: procedure for selecting solutions to perform previous actions Unlike DNA, simply changing "bits" of an MILP solution will not generate a new feasible solution need to define a method for adressing this issue 15 / 23

17 An evolutionary algorithm within MILP Evolutionary algorithm: ensuring feasibility Combination and mutation steps follow a similar procedure: select some variables to change (permute if binary) solve the sub-milp using RINS (or other heuristics) For mutation, an adaptive approach for choosing the fraction of variables to change is given. (a) Combination process (b) Mutation process Figure: Taken from [Rothberg, 2007] 2 2 Rothberg, E. (2007). An evolutionary algorithm for polishing mixed integer programming solutions. INFORMS Journal on Computing, 19(4): / 23

18 An evolutionary algorithm within MILP Evolutionary algorithm: selection process When selecting solutions for combination & mutation a preference should be given to "good" solutions, but selecting systematically the best solutions will drastically reduce the diversity Suggested selection process: randomly pick first solution from pool, then second solution is selected from sub-pool of better solutions 17 / 23

19 An evolutionary algorithm within MILP Evolutionary algorithm: results Computationally expensive use for improving already acceptable solution (found with RINS, local branching and guided diving) Compare results from (i) default CPLEX, (2) RINS, local branching and guided diving and (3) polishing (evolutionary algorithm) Figure: Taken from [Rothberg, 2007] 18 / 23

20 An evolutionary algorithm within MILP Evolutionary algorithm: discussion Improves solutions for some problems, but not all lack of diversity or neighborhood size General framework applicable to many different problems Computationally expensive: solves many (60 for example) sub-milps cannot rely sollely on this Implemented in CPLEX 12.4 (and following) 19 / 23

21 Conclusion 1 Pivoting methods 2 Relaxation Induced Neighborhood Search (RINS) 3 An evolutionary algorithm within MILP 4 Conclusion 20 / 23

22 Conclusion Conclusion Various heuristics used for different stages of a problem solving need to choose one/some methods All heuristics have pros and cons, and should be used for an adapted problem. We are just scratching the surface" in terms of all existing heuristics. 21 / 23

23 Conclusion Thank you for you attention Questions? 22 / 23

24 Conclusion References Balas, E. and Martin, C. H. (1980). Pivot and complement a heuristic for 0-1 programming. Management Science. Danna, E., Rothberg, E., and Le Pape, C. (2005). Exploring relaxation induced neighborhoods to improve mip solutions. Mathematical Programming, 102(1): Fischetti, M. and Lodi, A. (2011). Heuristics in mixed integer programming. Wiley Encyclopedia of Operations Research and Management Science. Rothberg, E. (2007). An evolutionary algorithm for polishing mixed integer programming solutions. INFORMS Journal on Computing, 19(4): / 23

Primal Heuristics in SCIP

Primal 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 information

Column Generation Based Primal Heuristics

Column 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 information

Primal Heuristics for Branch-and-Price Algorithms

Primal Heuristics for Branch-and-Price Algorithms 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

More information

Restrict-and-relax search for 0-1 mixed-integer programs

Restrict-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 information

15.083J Integer Programming and Combinatorial Optimization Fall Enumerative Methods

15.083J Integer Programming and Combinatorial Optimization Fall Enumerative Methods 5.8J Integer Programming and Combinatorial Optimization Fall 9 A knapsack problem Enumerative Methods Let s focus on maximization integer linear programs with only binary variables For example: a knapsack

More information

3 INTEGER LINEAR PROGRAMMING

3 INTEGER LINEAR PROGRAMMING 3 INTEGER LINEAR PROGRAMMING PROBLEM DEFINITION Integer linear programming problem (ILP) of the decision variables x 1,..,x n : (ILP) subject to minimize c x j j n j= 1 a ij x j x j 0 x j integer n j=

More information

Branch-and-cut implementation of Benders decomposition Matteo Fischetti, University of Padova

Branch-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 information

Randomized Rounding: A Primal Heuristic for General Mixed Integer Programming Problems

Randomized Rounding: A Primal Heuristic for General Mixed Integer Programming Problems Randomized Rounding: A Primal Heuristic for General Mixed Integer Programming Problems Mahdi Namazifar, Robin Lougee-Heimer, Andrew Miller, and John Forrest August 2009 Abstract. We propose an algorithm

More information

A hard integer program made easy by lexicography

A hard integer program made easy by lexicography Noname manuscript No. (will be inserted by the editor) A hard integer program made easy by lexicography Egon Balas Matteo Fischetti Arrigo Zanette February 16, 2011 Abstract A small but notoriously hard

More information

lpsymphony - Integer Linear Programming in R

lpsymphony - Integer Linear Programming in R lpsymphony - Integer Linear Programming in R Vladislav Kim October 30, 2017 Contents 1 Introduction 2 2 lpsymphony: Quick Start 2 3 Integer Linear Programming 5 31 Equivalent and Dual Formulations 5 32

More information

SUBSTITUTING GOMORY CUTTING PLANE METHOD TOWARDS BALAS ALGORITHM FOR SOLVING BINARY LINEAR PROGRAMMING

SUBSTITUTING GOMORY CUTTING PLANE METHOD TOWARDS BALAS ALGORITHM FOR SOLVING BINARY LINEAR PROGRAMMING ASIAN JOURNAL OF MATHEMATICS AND APPLICATIONS Volume 2014, Article ID ama0156, 11 pages ISSN 2307-7743 http://scienceasia.asia SUBSTITUTING GOMORY CUTTING PLANE METHOD TOWARDS BALAS ALGORITHM FOR SOLVING

More information

Improving CP-based Local Branching via Sliced Neighborhood Search

Improving CP-based Local Branching via Sliced Neighborhood Search Improving CP-based Local Branching via Sliced Neighborhood Search Fabio Parisini D.E.I.S., University of Bologna, Italy fabio.parisini@unibo.it Michela Milano D.E.I.S., University of Bologna, Italy michela.milano@unibo.it

More information

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

Introduction 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 information

Pure Cutting Plane Methods for ILP: a computational perspective

Pure Cutting Plane Methods for ILP: a computational perspective Pure Cutting Plane Methods for ILP: a computational perspective Matteo Fischetti, DEI, University of Padova Rorschach test for OR disorders: can you see the tree? 1 Outline 1. Pure cutting plane methods

More information

Stochastic Separable Mixed-Integer Nonlinear Programming via Nonconvex Generalized Benders Decomposition

Stochastic Separable Mixed-Integer Nonlinear Programming via Nonconvex Generalized Benders Decomposition Stochastic Separable Mixed-Integer Nonlinear Programming via Nonconvex Generalized Benders Decomposition Xiang Li Process Systems Engineering Laboratory Department of Chemical Engineering Massachusetts

More information

SUBSTITUTING GOMORY CUTTING PLANE METHOD TOWARDS BALAS ALGORITHM FOR SOLVING BINARY LINEAR PROGRAMMING

SUBSTITUTING GOMORY CUTTING PLANE METHOD TOWARDS BALAS ALGORITHM FOR SOLVING BINARY LINEAR PROGRAMMING Bulletin of Mathematics Vol. 06, No. 0 (20), pp.. SUBSTITUTING GOMORY CUTTING PLANE METHOD TOWARDS BALAS ALGORITHM FOR SOLVING BINARY LINEAR PROGRAMMING Eddy Roflin, Sisca Octarina, Putra B. J Bangun,

More information

M2 ORO: Advanced Integer Programming. Part IV. Solving MILP (1) easy IP. Outline. Sophie Demassey. October 10, 2011

M2 ORO: Advanced Integer Programming. Part IV. Solving MILP (1) easy IP. Outline. Sophie Demassey. October 10, 2011 M2 ORO: Advanced Integer Programming Sophie Demassey Part IV Solving MILP (1) Mines Nantes - TASC - INRIA/LINA CNRS UMR 6241 sophie.demassey@mines-nantes.fr October 10, 2011 Université de Nantes / M2 ORO

More information

The goal of this paper is to develop models and methods that use complementary

The goal of this paper is to develop models and methods that use complementary for a Class of Optimization Problems Vipul Jain Ignacio E. Grossmann Department of Chemical Engineering, Carnegie Mellon University, Pittsburgh, Pennsylvania, 15213, USA Vipul_Jain@i2.com grossmann@cmu.edu

More information

Cutting Planes for Some Nonconvex Combinatorial Optimization Problems

Cutting Planes for Some Nonconvex Combinatorial Optimization Problems Cutting Planes for Some Nonconvex Combinatorial Optimization Problems Ismael Regis de Farias Jr. Department of Industrial Engineering Texas Tech Summary Problem definition Solution strategy Multiple-choice

More information

Discrete Optimization with Decision Diagrams

Discrete Optimization with Decision Diagrams Discrete Optimization with Decision Diagrams J. N. Hooker Joint work with David Bergman, André Ciré, Willem van Hoeve Carnegie Mellon University Australian OR Society, May 2014 Goal Find an alternative

More information

A Local Dominance Procedure for Mixed-Integer Linear Programming

A Local Dominance Procedure for Mixed-Integer Linear Programming A Local Dominance Procedure for Mixed-Integer Linear Programming Matteo Fischetti ( ) and Domenico Salvagnin ( ) ( ) DEI, University of Padova, Italy ( ) DMPA, University of Padova, Italy e-mail: matteo.fischetti@unipd.it,

More information

Cloud Branching MIP workshop, Ohio State University, 23/Jul/2014

Cloud 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 information

A Steady-State Genetic Algorithm for Traveling Salesman Problem with Pickup and Delivery

A Steady-State Genetic Algorithm for Traveling Salesman Problem with Pickup and Delivery A Steady-State Genetic Algorithm for Traveling Salesman Problem with Pickup and Delivery Monika Sharma 1, Deepak Sharma 2 1 Research Scholar Department of Computer Science and Engineering, NNSS SGI Samalkha,

More information

Column Generation based Primal Heuristics

Column 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 information

Last topic: Summary; Heuristics and Approximation Algorithms Topics we studied so far:

Last topic: Summary; Heuristics and Approximation Algorithms Topics we studied so far: Last topic: Summary; Heuristics and Approximation Algorithms Topics we studied so far: I Strength of formulations; improving formulations by adding valid inequalities I Relaxations and dual problems; obtaining

More information

Evolutionary Computation Algorithms for Cryptanalysis: A Study

Evolutionary Computation Algorithms for Cryptanalysis: A Study Evolutionary Computation Algorithms for Cryptanalysis: A Study Poonam Garg Information Technology and Management Dept. Institute of Management Technology Ghaziabad, India pgarg@imt.edu Abstract The cryptanalysis

More information

Basic Concepts of Constraint Integer Programming

Basic Concepts of Constraint Integer Programming Basic Concepts of Constraint Integer Programming Ambros Gleixner Zuse Institute Berlin September 30, 2015 Outline SCIP Solving Constraint Integer Programs 4 methodologies in optimization An integrated

More information

A NEW MILP APPROACH FOR THE FACILITY LAYOUT DESIGN PROBLEM WITH RECTANGULAR AND L/T SHAPED DEPARTMENTS

A NEW MILP APPROACH FOR THE FACILITY LAYOUT DESIGN PROBLEM WITH RECTANGULAR AND L/T SHAPED DEPARTMENTS A NEW MILP APPROACH FOR THE FACILITY LAYOUT DESIGN PROBLEM WITH RECTANGULAR AND L/T SHAPED DEPARTMENTS Yossi Bukchin Michal Tzur Dept. of Industrial Engineering, Tel Aviv University, ISRAEL Abstract In

More information

An Introduction to Evolutionary Algorithms

An Introduction to Evolutionary Algorithms An Introduction to Evolutionary Algorithms Karthik Sindhya, PhD Postdoctoral Researcher Industrial Optimization Group Department of Mathematical Information Technology Karthik.sindhya@jyu.fi http://users.jyu.fi/~kasindhy/

More information

LaGO. Ivo Nowak and Stefan Vigerske. Humboldt-University Berlin, Department of Mathematics

LaGO. Ivo Nowak and Stefan Vigerske. Humboldt-University Berlin, Department of Mathematics LaGO a Branch and Cut framework for nonconvex MINLPs Ivo Nowak and Humboldt-University Berlin, Department of Mathematics EURO XXI, July 5, 2006 21st European Conference on Operational Research, Reykjavik

More information

A NEW SEQUENTIAL CUTTING PLANE ALGORITHM FOR SOLVING MIXED INTEGER NONLINEAR PROGRAMMING PROBLEMS

A NEW SEQUENTIAL CUTTING PLANE ALGORITHM FOR SOLVING MIXED INTEGER NONLINEAR PROGRAMMING PROBLEMS EVOLUTIONARY METHODS FOR DESIGN, OPTIMIZATION AND CONTROL P. Neittaanmäki, J. Périaux and T. Tuovinen (Eds.) c CIMNE, Barcelona, Spain 2007 A NEW SEQUENTIAL CUTTING PLANE ALGORITHM FOR SOLVING MIXED INTEGER

More information

A Genetic Algorithm Framework

A Genetic Algorithm Framework Fast, good, cheap. Pick any two. The Project Triangle 3 A Genetic Algorithm Framework In this chapter, we develop a genetic algorithm based framework to address the problem of designing optimal networks

More information

Welcome to the Webinar. What s New in Gurobi 7.5

Welcome to the Webinar. What s New in Gurobi 7.5 Welcome to the Webinar What s New in Gurobi 7.5 Speaker Introduction Dr. Tobias Achterberg Director of R&D at Gurobi Optimization Formerly a developer at ILOG, where he worked on CPLEX 11.0 to 12.6 Obtained

More information

On Truck dock assignment problem with operational time constraint within cross docks

On Truck dock assignment problem with operational time constraint within cross docks On Truck dock assignment problem with operational time constraint within cross docks Shahin Gelareh a,b,, Gilles Goncalves a,b, Rahimeh N. Monemi c,1, a UArtois, LGI2A, F-62400, Béthune, France b Univ

More information

/ Approximation Algorithms Lecturer: Michael Dinitz Topic: Linear Programming Date: 2/24/15 Scribe: Runze Tang

/ Approximation Algorithms Lecturer: Michael Dinitz Topic: Linear Programming Date: 2/24/15 Scribe: Runze Tang 600.469 / 600.669 Approximation Algorithms Lecturer: Michael Dinitz Topic: Linear Programming Date: 2/24/15 Scribe: Runze Tang 9.1 Linear Programming Suppose we are trying to approximate a minimization

More information

On the Optimization of CPLEX Models

On the Optimization of CPLEX Models International Research Journal of Applied and Basic Sciences 3 Available online at www.irjabs.com ISSN 5-838X / Vol, 4 (9): 8-86 Science Explorer Publications On the Optimization of CPLEX Models Mohamad

More information

Minimum Weight Constrained Forest Problems. Problem Definition

Minimum Weight Constrained Forest Problems. Problem Definition Slide 1 s Xiaoyun Ji, John E. Mitchell Department of Mathematical Sciences Rensselaer Polytechnic Institute Troy, NY, USA jix@rpi.edu, mitchj@rpi.edu 2005 Optimization Days Montreal, Canada May 09, 2005

More information

A simulated annealing algorithm for the vehicle routing problem with time windows and synchronization constraints

A simulated annealing algorithm for the vehicle routing problem with time windows and synchronization constraints A simulated annealing algorithm for the vehicle routing problem with time windows and synchronization constraints Sohaib Afifi 1, Duc-Cuong Dang 1,2, and Aziz Moukrim 1 1 Université de Technologie de Compiègne

More information

MINLP applications, part II: Water Network Design and some applications of black-box optimization

MINLP applications, part II: Water Network Design and some applications of black-box optimization MINLP applications, part II: Water Network Design and some applications of black-box optimization Claudia D Ambrosio CNRS & LIX, École Polytechnique dambrosio@lix.polytechnique.fr 5th Porto Meeting on

More information

The Mixed Integer Linear Programming Solver

The Mixed Integer Linear Programming Solver SAS/OR 14.3 User s Guide Mathematical Programming The Mixed Integer Linear Programming Solver This document is an individual chapter from SAS/OR 14.3 User s Guide: Mathematical Programming. The correct

More information

Binary decision diagrams for computing the non-dominated set

Binary decision diagrams for computing the non-dominated set Binary decision diagrams for computing the non-dominated set July 13, 2015 Antti Toppila and Ahti Salo 27th European Conference on Operational Research, 12-15 July 2015, University of Strathclyde, Glasgow,

More information

A Branch-and-Cut Algorithm for the Partition Coloring Problem

A Branch-and-Cut Algorithm for the Partition Coloring Problem A Branch-and-Cut Algorithm for the Partition Coloring Problem Yuri Frota COPPE/UFRJ, Programa de Engenharia de Sistemas e Otimização Rio de Janeiro, RJ 21945-970, Brazil abitbol@cos.ufrj.br Nelson Maculan

More information

Optimization of Complex Systems with OptQuest

Optimization of Complex Systems with OptQuest Optimization of Complex Systems with OptQuest MANUEL LAGUNA Graduate School of Business Administration University of Colorado, Boulder, CO 80309-0419 Manuel.Laguna@Colorado.EDU Latest revision: April 8,

More information

L2: Algorithms: Knapsack Problem & BnB

L2: Algorithms: Knapsack Problem & BnB L2: Algorithms: Knapsack Problem & BnB This tutorial covers the basic topics on creating a forms application, common form controls and the user interface for the optimization models, algorithms and heuristics,

More information

Critical Node Detection Problem. Panos Pardalos Distinguished Professor CAO, Dept. of Industrial and Systems Engineering, University of Florida

Critical Node Detection Problem. Panos Pardalos Distinguished Professor CAO, Dept. of Industrial and Systems Engineering, University of Florida Critical Node Detection Problem ITALY May, 2008 Panos Pardalos Distinguished Professor CAO, Dept. of Industrial and Systems Engineering, University of Florida Outline of Talk Introduction Problem Definition

More information

Manpower Planning: Task Scheduling. Anders Høeg Dohn

Manpower Planning: Task Scheduling. Anders Høeg Dohn : Task Scheduling Anders Høeg Dohn Scope During these lectures I will: Go over some of the practical problems encountered in manpower planning. Rostering Task Scheduling Propose models that can be used

More information

A Parallel Architecture for the Generalized Travelling Salesman Problem: Project Proposal

A Parallel Architecture for the Generalized Travelling Salesman Problem: Project Proposal A Parallel Architecture for the Generalized Travelling Salesman Problem: Project Proposal Max Scharrenbroich, maxfs at umd.edu Dr. Bruce Golden, R. H. Smith School of Business, bgolden at rhsmith.umd.edu

More information

The SYMPHONY Callable Library for Mixed-Integer Linear Programming

The SYMPHONY Callable Library for Mixed-Integer Linear Programming The SYMPHONY Callable Library for Mixed-Integer Linear Programming Ted Ralphs and Menal Guzelsoy Industrial and Systems Engineering Lehigh University INFORMS Computing Society Conference, Annapolis, MD,

More information

Integer Programming Chapter 9

Integer Programming Chapter 9 Integer Programming Chapter 9 University of Chicago Booth School of Business Kipp Martin October 25, 2017 1 / 40 Outline Key Concepts MILP Set Monoids LP set Relaxation of MILP Set Formulation Quality

More information

Algorithms for the Precedence Constrained Generalized Travelling Salesperson Problem

Algorithms for the Precedence Constrained Generalized Travelling Salesperson Problem MASTER S THESIS Algorithms for the Precedence Constrained Generalized Travelling Salesperson Problem RAAD SALMAN Department of Mathematical Sciences Division of Mathematics CHALMERS UNIVERSITY OF TECHNOLOGY

More information

GAUSSIAN VARIABLE NEIGHBORHOOD SEARCH FOR THE FILE TRANSFER SCHEDULING PROBLEM

GAUSSIAN VARIABLE NEIGHBORHOOD SEARCH FOR THE FILE TRANSFER SCHEDULING PROBLEM Yugoslav Journal of Operations Research 26 (2016), Number 2, 173 188 DOI: 10.2298/YJOR150124006D GAUSSIAN VARIABLE NEIGHBORHOOD SEARCH FOR THE FILE TRANSFER SCHEDULING PROBLEM Zorica DRAŽIĆ Faculty of

More information

Cutting Stock with Binary Patterns: Arc-flow Formulation with Graph Compression

Cutting Stock with Binary Patterns: Arc-flow Formulation with Graph Compression Cutting Stock with Binary Patterns: Arc-flow Formulation with Graph Compression Filipe Brandão INESC TEC and Faculdade de Ciências, Universidade do Porto, Portugal fdabrandao@dcc.fc.up.pt arxiv:1502.02899v1

More information

Hill Climbing. Assume a heuristic value for each assignment of values to all variables. Maintain an assignment of a value to each variable.

Hill Climbing. Assume a heuristic value for each assignment of values to all variables. Maintain an assignment of a value to each variable. Hill Climbing Many search spaces are too big for systematic search. A useful method in practice for some consistency and optimization problems is hill climbing: Assume a heuristic value for each assignment

More information

Minimal Infeasible Subsystems and Benders cuts

Minimal Infeasible Subsystems and Benders cuts Mathematical Programming manuscript No. (will be inserted by the editor) Matteo Fischetti Domenico Salvagnin Arrigo Zanette Minimal Infeasible Subsystems and Benders cuts Received: date / Revised Sept.

More information

Implementing Constraint Handlers in SCIP

Implementing Constraint Handlers in SCIP Implementing Constraint Handlers in SCIP Gregor Hendel Ambros Gleixner, Felipe Serrano September 30 Outline Why use constraints? Motivation Callback: Default Constraint Handlers of SCIP Implementation

More information

Computing Crossing-Free Configurations with Minimum Bottleneck

Computing Crossing-Free Configurations with Minimum Bottleneck Computing Crossing-Free Configurations with Minimum Bottleneck Sándor P. Fekete 1 and Phillip Keldenich 1 1 Department of Computer Science, TU Braunschweig, Germany {s.fekete,p.keldenich}@tu-bs.de Abstract

More information

Rounding and Propagation Heuristics for Mixed Integer Programming

Rounding 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 information

Operations Research and Optimization: A Primer

Operations Research and Optimization: A Primer Operations Research and Optimization: A Primer Ron Rardin, PhD NSF Program Director, Operations Research and Service Enterprise Engineering also Professor of Industrial Engineering, Purdue University Introduction

More information

Hybrid Enumeration Strategies for Mixed Integer Programming

Hybrid Enumeration Strategies for Mixed Integer Programming Hybrid Enumeration Strategies for Mixed Integer Programming João Pedro Pedroso Technical Report Series: DCC-2004-8 Departamento de Ciência de Computadores Faculdade de Ciências & Laboratório de Inteligência

More information

A Row-and-Column Generation Method to a Batch Machine Scheduling Problem

A Row-and-Column Generation Method to a Batch Machine Scheduling Problem The Ninth International Symposium on Operations Research and Its Applications (ISORA 10) Chengdu-Jiuzhaigou, China, August 19 23, 2010 Copyright 2010 ORSC & APORC, pp. 301 308 A Row-and-Column Generation

More information

IEEE TRANSACTIONS ON AUTOMATION SCIENCE AND ENGINEERING, VOL. 5, NO. 4, OCTOBER

IEEE TRANSACTIONS ON AUTOMATION SCIENCE AND ENGINEERING, VOL. 5, NO. 4, OCTOBER IEEE TRANSACTIONS ON AUTOMATION SCIENCE AND ENGINEERING, VOL. 5, NO. 4, OCTOBER 2008 573 Hybrid Nested Partitions and Mathematical Programming Approach and Its Applications Liang Pi, Student Member, IEEE,

More information

Pre-requisite Material for Course Heuristics and Approximation Algorithms

Pre-requisite Material for Course Heuristics and Approximation Algorithms Pre-requisite Material for Course Heuristics and Approximation Algorithms This document contains an overview of the basic concepts that are needed in preparation to participate in the course. In addition,

More information

THE objective of the redundancy allocation problem is to

THE objective of the redundancy allocation problem is to IEEE TRANSACTIONS ON RELIABILITY, VOL. 55, NO. 3, SEPTEMBER 2006 551 Multiple Weighted Objectives Heuristic for the Redundancy Allocation Problem David W. Coit and Abdullah Konak Abstract A new heuristic

More information

Finding an optimal seating arrangement for employees traveling to an event

Finding an optimal seating arrangement for employees traveling to an event Croatian Operational Research Review 419 CRORR 6(2015), 419 427 Finding an optimal seating arrangement for employees traveling to an event Ninoslav Čerkez 1, Rebeka Čorić 2, Mateja Dumić 2, and Domagoj

More information

Solving the Euclidean Steiner Tree Problem in n-space

Solving the Euclidean Steiner Tree Problem in n-space Solving the Euclidean Steiner Tree Problem in n-space Marcia Fampa (UFRJ), Jon Lee (U. Michigan), and Wendel Melo (UFRJ) January 2015 Marcia Fampa, Jon Lee, Wendel Melo Solving the Euclidean Steiner Tree

More information

(Duality), Warm Starting, and Sensitivity Analysis for MILP

(Duality), Warm Starting, and Sensitivity Analysis for MILP (Duality), Warm Starting, and Sensitivity Analysis for MILP Ted Ralphs and Menal Guzelsoy Industrial and Systems Engineering Lehigh University INFORMS Annual Conference, Denver, CO, Tuesday, October 26,

More information

COVERING POINTS WITH AXIS PARALLEL LINES. KAWSAR JAHAN Bachelor of Science, Bangladesh University of Professionals, 2009

COVERING POINTS WITH AXIS PARALLEL LINES. KAWSAR JAHAN Bachelor of Science, Bangladesh University of Professionals, 2009 COVERING POINTS WITH AXIS PARALLEL LINES KAWSAR JAHAN Bachelor of Science, Bangladesh University of Professionals, 2009 A Thesis Submitted to the School of Graduate Studies of the University of Lethbridge

More information

Introduction. Chapter 15. Optimization Modeling: Applications. Integer Programming. Manufacturing Example. Three Types of ILP Models

Introduction. Chapter 15. Optimization Modeling: Applications. Integer Programming. Manufacturing Example. Three Types of ILP Models Chapter 5 Optimization Modeling: Applications Integer Programming Introduction When one or more variables in an LP problem must assume an integer value we have an Integer Linear Programming (ILP) problem.

More information

An Evolutionary Algorithm for the Multi-objective Shortest Path Problem

An Evolutionary Algorithm for the Multi-objective Shortest Path Problem An Evolutionary Algorithm for the Multi-objective Shortest Path Problem Fangguo He Huan Qi Qiong Fan Institute of Systems Engineering, Huazhong University of Science & Technology, Wuhan 430074, P. R. China

More information

Chapter 14: Optimal Winner Determination Algorithms

Chapter 14: Optimal Winner Determination Algorithms Chapter 14: Optimal Winner Determination Algorithms Tuomas Sandholm 1 Introduction This chapter discusses optimal winner determination algorithms for combinatorial auctions (CAs). We say the auctioneer

More information

Review of Mixed-Integer Nonlinear and Generalized Disjunctive Programming Methods

Review of Mixed-Integer Nonlinear and Generalized Disjunctive Programming Methods Carnegie Mellon University Research Showcase @ CMU Department of Chemical Engineering Carnegie Institute of Technology 2-2014 Review of Mixed-Integer Nonlinear and Generalized Disjunctive Programming Methods

More information

Package Rcplex. February 15, 2013

Package Rcplex. February 15, 2013 Version 0.3-1 Date 2013-02-13 Title R interface to CPLEX Package Rcplex February 15, 2013 Author Hector Corrada Bravo, with contributions from Stefan Theussl and Kurt Hornik Maintainer Hector Corrada Bravo

More information

AN EXACT ALGORITHM FOR A LOCATION-ROUTING PROBLEM. Martine Labbé 1, Inmaculada Rodríguez Martín 2, Juan José Salazar González 2

AN EXACT ALGORITHM FOR A LOCATION-ROUTING PROBLEM. Martine Labbé 1, Inmaculada Rodríguez Martín 2, Juan José Salazar González 2 27 Congreso Nacional de Estadística e Investigación Operativa Lleida, 8 11 de abril de 2003 AN EXACT ALGORITHM FOR A LOCATION-ROUTING PROBLEM Martine Labbé 1, Inmaculada Rodríguez Martín 2, Juan José Salazar

More information

IMPROVING A GREEDY DNA MOTIF SEARCH USING A MULTIPLE GENOMIC SELF-ADAPTATING GENETIC ALGORITHM

IMPROVING A GREEDY DNA MOTIF SEARCH USING A MULTIPLE GENOMIC SELF-ADAPTATING GENETIC ALGORITHM Proceedings of Student/Faculty Research Day, CSIS, Pace University, May 4th, 2007 IMPROVING A GREEDY DNA MOTIF SEARCH USING A MULTIPLE GENOMIC SELF-ADAPTATING GENETIC ALGORITHM Michael L. Gargano, mgargano@pace.edu

More information

Introduction to Linear Programming. Algorithmic and Geometric Foundations of Optimization

Introduction to Linear Programming. Algorithmic and Geometric Foundations of Optimization Introduction to Linear Programming Algorithmic and Geometric Foundations of Optimization Optimization and Linear Programming Mathematical programming is a class of methods for solving problems which ask

More information

Assignment 3b: The traveling salesman problem

Assignment 3b: The traveling salesman problem Chalmers University of Technology MVE165 University of Gothenburg MMG631 Mathematical Sciences Linear and integer optimization Optimization with applications Emil Gustavsson Assignment information Ann-Brith

More information

DIT411/TIN175, Artificial Intelligence. Peter Ljunglöf. 23 January, 2018

DIT411/TIN175, Artificial Intelligence. Peter Ljunglöf. 23 January, 2018 DIT411/TIN175, Artificial Intelligence Chapters 3 4: More search algorithms CHAPTERS 3 4: MORE SEARCH ALGORITHMS DIT411/TIN175, Artificial Intelligence Peter Ljunglöf 23 January, 2018 1 TABLE OF CONTENTS

More information

Artificial Intelligence

Artificial Intelligence Artificial Intelligence Local Search Vibhav Gogate The University of Texas at Dallas Some material courtesy of Luke Zettlemoyer, Dan Klein, Dan Weld, Alex Ihler, Stuart Russell, Mausam Systematic Search:

More information

What is GOSET? GOSET stands for Genetic Optimization System Engineering Tool

What is GOSET? GOSET stands for Genetic Optimization System Engineering Tool Lecture 5: GOSET 1 What is GOSET? GOSET stands for Genetic Optimization System Engineering Tool GOSET is a MATLAB based genetic algorithm toolbox for solving optimization problems 2 GOSET Features Wide

More information

Solving Traveling Salesman Problem for Large Spaces using Modified Meta- Optimization Genetic Algorithm

Solving Traveling Salesman Problem for Large Spaces using Modified Meta- Optimization Genetic Algorithm Solving Traveling Salesman Problem for Large Spaces using Modified Meta- Optimization Genetic Algorithm Maad M. Mijwel Computer science, college of science, Baghdad University Baghdad, Iraq maadalnaimiy@yahoo.com

More information

Informed Search and Exploration

Informed Search and Exploration Informed Search and Exploration Chapter 4 (4.1-4.3) CS 2710 1 Introduction Ch.3 searches good building blocks for learning about search But vastly inefficient eg: Can we do better? Breadth Depth Uniform

More information

CSc 545 Lecture topic: The Criss-Cross method of Linear Programming

CSc 545 Lecture topic: The Criss-Cross method of Linear Programming CSc 545 Lecture topic: The Criss-Cross method of Linear Programming Wanda B. Boyer University of Victoria November 21, 2012 Presentation Outline 1 Outline 2 3 4 Please note: I would be extremely grateful

More information

CS5401 FS2015 Exam 1 Key

CS5401 FS2015 Exam 1 Key CS5401 FS2015 Exam 1 Key This is a closed-book, closed-notes exam. The only items you are allowed to use are writing implements. Mark each sheet of paper you use with your name and the string cs5401fs2015

More information

SCIP. 1 Introduction. 2 Model requirements. Contents. Stefan Vigerske, Humboldt University Berlin, Germany

SCIP. 1 Introduction. 2 Model requirements. Contents. Stefan Vigerske, Humboldt University Berlin, Germany SCIP Stefan Vigerske, Humboldt University Berlin, Germany Contents 1 Introduction.................................................. 673 2 Model requirements..............................................

More information

MINIZSAT. A semi SAT-based pseudo-boolean solver. Master thesis by Rogier Poldner

MINIZSAT. A semi SAT-based pseudo-boolean solver. Master thesis by Rogier Poldner MINIZSAT A semi SAT-based pseudo-boolean solver Master thesis by Rogier Poldner MINIZSAT A semi SAT-based pseudo-boolean solver Master thesis by Rogier Poldner committee: Dr. H. van Maaren Dr. M.J.H.

More information

Local search heuristic for multiple knapsack problem

Local search heuristic for multiple knapsack problem International Journal of Intelligent Information Systems 2015; 4(2): 35-39 Published online February 14, 2015 (http://www.sciencepublishinggroup.com/j/ijiis) doi: 10.11648/j.ijiis.20150402.11 ISSN: 2328-7675

More information

ALGORITHMS AND SOFTWARE FOR CONVEX MIXED INTEGER NONLINEAR PROGRAMS

ALGORITHMS AND SOFTWARE FOR CONVEX MIXED INTEGER NONLINEAR PROGRAMS ALGORITHMS AND SOFTWARE FOR CONVEX MIXED INTEGER NONLINEAR PROGRAMS PIERRE BONAMI, MUSTAFA KILINÇ, AND JEFF LINDEROTH Abstract. This paper provides a survey of recent progress and software for solving

More information

Automatic Generation of Prime Factorization Algorithms Using Genetic Programming

Automatic Generation of Prime Factorization Algorithms Using Genetic Programming 1 Automatic Generation of Prime Factorization Algorithms Using Genetic Programming David Michael Chan Department of Computer Science Post Office Box 12587 Stanford, California 94309 dmchan@stanford.edu

More information

AMPL Integer Linear Programming And Sensitivity Analysis for LP. ORLAB Operations Research Laboratory. Borzou Rostami. Politecnico di Milano, Italy

AMPL Integer Linear Programming And Sensitivity Analysis for LP. ORLAB Operations Research Laboratory. Borzou Rostami. Politecnico di Milano, Italy AMPL Integer Linear Programming And Sensitivity Analysis for LP ORLAB Operations Research Laboratory Borzou Rostami Politecnico di Milano, Italy December 6, 2012 Integer Programming: Many linear programming

More information

Genetic programming. Lecture Genetic Programming. LISP as a GP language. LISP structure. S-expressions

Genetic programming. Lecture Genetic Programming. LISP as a GP language. LISP structure. S-expressions Genetic programming Lecture Genetic Programming CIS 412 Artificial Intelligence Umass, Dartmouth One of the central problems in computer science is how to make computers solve problems without being explicitly

More information

Mechanical Component Design for Multiple Objectives Using Elitist Non-Dominated Sorting GA

Mechanical Component Design for Multiple Objectives Using Elitist Non-Dominated Sorting GA Mechanical Component Design for Multiple Objectives Using Elitist Non-Dominated Sorting GA Kalyanmoy Deb, Amrit Pratap, and Subrajyoti Moitra Kanpur Genetic Algorithms Laboratory (KanGAL) Indian Institute

More information

SBB: A New Solver for Mixed Integer Nonlinear Programming

SBB: A New Solver for Mixed Integer Nonlinear Programming SBB: A New Solver for Mixed Integer Nonlinear Programming Michael R. Bussieck GAMS Development Corp. Arne Drud ARKI Consulting & Development A/S Overview Introduction: The MINLP Model The B&B Algorithm

More information

The Branch & Move algorithm: Improving Global Constraints Support by Local Search

The Branch & Move algorithm: Improving Global Constraints Support by Local Search Branch and Move 1 The Branch & Move algorithm: Improving Global Constraints Support by Local Search Thierry Benoist Bouygues e-lab, 1 av. Eugène Freyssinet, 78061 St Quentin en Yvelines Cedex, France tbenoist@bouygues.com

More information

A List Heuristic for Vertex Cover

A List Heuristic for Vertex Cover A List Heuristic for Vertex Cover Happy Birthday Vasek! David Avis McGill University Tomokazu Imamura Kyoto University Operations Research Letters (to appear) Online: http://cgm.cs.mcgill.ca/ avis revised:

More information

Using Database Storage to Improve Explorative Optimization of Form Critical Structures

Using Database Storage to Improve Explorative Optimization of Form Critical Structures 15 to 19 September 2014, Brasilia, Brazil Reyolando M.L.R.F. BRASIL and Ruy M.O. PAULETTI (eds.) Using Database Storage to Improve Explorative Optimization of Form Critical Structures VON BUELOW, Peter

More information

Construction of Minimum-Weight Spanners Mikkel Sigurd Martin Zachariasen

Construction 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 information

arxiv: v2 [cs.dm] 3 Dec 2014

arxiv: v2 [cs.dm] 3 Dec 2014 The Student/Project Allocation problem with group projects Aswhin Arulselvan, Ágnes Cseh, and Jannik Matuschke arxiv:4.035v [cs.dm] 3 Dec 04 Department of Management Science, University of Strathclyde,

More information

Mixed Integer Programming Class Library (MIPCL)

Mixed 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 information

Outline. Optimales Recycling - Tourenplanung in der Altglasentsorgung

Outline. Optimales Recycling - Tourenplanung in der Altglasentsorgung 1 Optimales Recycling - Ruhr-Universität Bochum, 15.02.2013 2 1. Introduction and Motivation 2. Problem Definition 3. Literature Review 4. Mathematical Model 5. Variable Neighborhood Search 6. Numerical

More information

Using the OPTMODEL Procedure in SAS/OR to Solve Complex Problems

Using the OPTMODEL Procedure in SAS/OR to Solve Complex Problems Paper SAS1502-2015 Using the OPTMODEL Procedure in SAS/OR to Solve Complex Problems Rob Pratt, SAS Institute Inc. ABSTRACT Mathematical optimization is a powerful paradigm for modeling and solving business

More information