Facility Location-Routing-Scheduling Problem: LRS

Size: px
Start display at page:

Download "Facility Location-Routing-Scheduling Problem: LRS"

Transcription

1 Facility Location-Routing-Scheduling Problem: LRS Zeliha Akca Lehigh University Rosemary Berger Lehigh University ➀ Definition of LRS ➁ Branch and Price Algorithm ➂ Mater and Pricing Problem ➃ Solution of Pricing Problem ➄ What We Have Done ➅ What We Will Do June 15,

2 LRS Location Routing and Scheduling Problem: 3 dependent problems: ➀ locate facilities ➁ construct routes for vehicles ➂ assign routes to vehicles capacitated facilities capacitated vehicles time restriction for the vehicles RS 2

3 Location-Routing and Scheduling Problem: In literature, heuristic solution for LRS problem (no IP formulation) Exact solutions for RS and LR We choose, Branch and Price Algorithm: IP formulation includes many constraints (s.t. sub tour elimination constraints) Can be written in set partitioning problem easily Easy to think routes in terms of columns With set partitioning formulation, many possible columns Other methods to solve: Lagrangian Relaxation Branch and bound and cut Heuristic design? RS 3

4 Problems in Literature Facility Location -too many Vehicle Routing -too many Routing Scheduling Location Routing Location Routing and Scheduling routing: Location one-to-one relation btw routes and vehicles. Location scheduling: not necessarily one-to-one relation assignment of one vehicle to many paths. RS 4

5 Objective: Minimize total cost. IP FORMULATION Constraints: TotalCost Fixed cost of Facility and Vehicle Operating cost of Vehicles ➀ Each demand node should be served once ➁ # of a vehicle entering a node must be equal to # of the vehicle leaves this node ➂ Capacity restriction for facility ➃ Capacity restriction for vehicles ➄ Flow balance equations (to satisfy demand and eliminate the subtours) ➅ Time restriction to the routes FORMULATION 5

6 Set Partitioning Model: ALTERNATE FORMULATION Pairing Concept: Set of routes assigned to a vehicle and can be served within the given time limit. PAIR Customer facility LTERNATE FORMULATION

7 #! $%=!$:set of customers; &!! Set Partitioning Model: ALTERNATE FORMULATION Variables for set partitioning based on pairing concept: 1 if pairing "otherwise is chosen for facility,, :set of feasible pairs of facility j 1 if facility j is open, "otherwise :set of facilities; LTERNATE FORMULATION 7

8 5 of " A 7 # /5 # 7 4 / =1 if node? # = )8 ) ;: 5 <( + >= + )* ) ' )* )* )3! (' #,+-.' / Set Partitioning Model: (1) s.t. $(2),!,!,!(3) (4) (5) 'is in pairing. LTERNATE FORMULATION 8

9 BRANCH AND PRICE ALGORITHM Restricted Master Prblem (RMP) Initial Columns Solve LP LP LOWER BD IP UPPER BD. Dual Variables info Until no new column new columns Pricing Problem Find columns with ( ) reduced cost ROOT NODE If LP is not integral Branching Modified RPM Cut is added Solve LP modified Pricing Problem RANCH AND PRICE ALGORITHM 9

10 Restricted Master Problem: Initial pairs are formed pair pair 3 Y13 Y12 2 pair Y14 4 pair Y11 fac pair Y21 fac.2 pair Y22 Each pair represent a column in set partitioning formulation Restricted-since includes set of columns, not all columns RANCH AND PRICE ALGORITHM 10

11 B Pricing problem: Create pair : a column for If 3rd const changed to: )3 4 5 = / / F B D E C ' #!and $() We have:,, dual variables Reduced Cost for HG (7) / )8 )8 )8 D 5 C 5 5 I ;: 5 <( Operating cost of the vehicle ( Jtravel time) + Fixed Cost of a vehicle Independent pricing problem for each facility RANCH AND PRICE ALGORITHM 11

12 K ELEMENTARY SHOR THEST PATH WITH RESOURCE CONSTRAINT Pricing Problem = ESPRC If: Set up a network, including all customers and a source and sink nodes Arc costs: LMON QSRT LM\[ e ^ f em (8) d P UVW ] M^ ba Z cy _ `R XY Z find minimum cost path to the sink in our case allow visits more than once to sink If Total cost of path + Vehicle fixed Cost gih, add the column to restricted master problem stop when the shorthest path does not give negative cost column LEMENTARY SHORTHEST PATH WITH RESOURCE CONSTRAINT 12

13 j 7 ELEMENTARY SHOR THEST PATH WITH RESOURCE CONSTRAINT What is an elementary path? Each node can be visited at most once. Why elementary instead of walks? Trade of between more difficult pricing problem and more depth in branch and bound tree In our case: # of visits to sink In each visit to sink, current truck load is set to zero Adapt the Labelling Algorithm for ESPRC by Feillet, Dejax, Gendreau, Geuguen. LEMENTARY SHORTHEST PATH WITH RESOURCE CONSTRAINT 13

14 ESPRC Problem: too many feasible paths Keep resource consumptions, visited nodes, and cost Keep unreachable nodes for each label A node may be unreachable from other if not enough resource or is already visited. Eliminate dominated labels with respect to resource consumption and unreachable nodes. SPRC 14

15 CURRENTLY What we have done Design Master Problem and Pricing Problem Adapted ESPRC algorithm to solve Pricing Problem Do the column generation Solve the root node URRENTLY 15

16 MINTO MINTO:Mixed INTeger Optimizer MINTO uses LP solver and do branch and bound algorithm MINTO can do many applications such as preprocessing, constraint generation, primal heurisitcs MINTO allows user to write own algorithm (for column generation, constraint generation, heuristics,..) specific to the problem Prof. Linderoth supports MINTO in our University WHAT ELSE? COIN-BCP:(Common Optimization INterface) and SYMPHONY Open source allows parallellization in branch and bound tree supported by Prof. Ralphs HAT ELSE? 1

17 NEXT What we will do More implementation Test problems Determine the right number of columns to be generated in each time Different LP algorithms, to find better reduced costs See how well the root node solution Create column pool Branching strategies parallelization Alternate solution: 2-sub problems approach IP formulation Focus on the pricing problem EXT 17

18 o n v p q rw u o n k Q p q t e q x e q m u m 2-SUB PROBLEM APPROACH Master problem includes 3 set of variables: Location variables 1 if pairing is chosen for facility, Qsr, p r el N hotherwise 1 if path wis chosen for facility,, p r el N hotherwise 2 nested sub problems: generating paths,and combining these paths as pairs. SP1: Generating paths: vehicle routing problem or elementary shortest path with 2 resources SP2: Combining paths: knapsack problem -SUB PROBLEM APPROACH 18

19 2-SUB PROBLEM ALGORITHM Solve LP RMP Get dual variables for each customer and facility Solve SP1: New Paths with ( ) cost Form pairs heuristically add to MP With new paths and old paths Solve Master Problem Solve SP2: generate pairs with ( ) cost Get dual variables -SUB PROBLEM ALGORITHM 19

20 THANKS... Any Questions? HANKS... 20

Solution Methods for the Multi-trip Elementary Shortest Path Problem with Resource Constraints

Solution Methods for the Multi-trip Elementary Shortest Path Problem with Resource Constraints Solution Methods for the Multi-trip Elementary Shortest Path Problem with Resource Constraints Zeliha Akca Ted K. Ralphs Rosemary T. Berger December 31, 2010 Abstract We investigate the multi-trip elementary

More information

Modeling and Solving Location Routing and Scheduling Problems

Modeling and Solving Location Routing and Scheduling Problems Modeling and Solving Location Routing and Scheduling Problems Z. Akca R.T. Berger T.K Ralphs October 13, 2008 Abstract This paper studies location routing and scheduling problems, a class of problems in

More information

A Branch-and-Price Algorithm for Combined Location and Routing Problems Under Capacity Restrictions

A Branch-and-Price Algorithm for Combined Location and Routing Problems Under Capacity Restrictions A Branch-and-Price Algorithm for Combined Location and Routing Problems Under Capacity Restrictions Z. Akca R.T. Berger T.K Ralphs September 17, 2008 Abstract We investigate the problem of simultaneously

More information

56:272 Integer Programming & Network Flows Final Examination -- December 14, 1998

56: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 information

A Bucket Graph Based Labelling Algorithm for the Resource Constrained Shortest Path Problem with Applications to Vehicle Routing

A Bucket Graph Based Labelling Algorithm for the Resource Constrained Shortest Path Problem with Applications to Vehicle Routing 1 / 23 A Bucket Graph Based Labelling Algorithm for the Resource Constrained Shortest Path Problem with Applications to Vehicle Routing Ruslan Sadykov 1,2 Artur Pessoa 3 Eduardo Uchoa 3 1 Inria Bordeaux,

More information

Introduction to Mathematical Programming IE496. Final Review. Dr. Ted Ralphs

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

Algorithms for Integer Programming

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

DETERMINISTIC OPERATIONS RESEARCH

DETERMINISTIC OPERATIONS RESEARCH DETERMINISTIC OPERATIONS RESEARCH Models and Methods in Optimization Linear DAVID J. RADER, JR. Rose-Hulman Institute of Technology Department of Mathematics Terre Haute, IN WILEY A JOHN WILEY & SONS,

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

Solutions for Operations Research Final Exam

Solutions for Operations Research Final Exam Solutions for Operations Research Final Exam. (a) The buffer stock is B = i a i = a + a + a + a + a + a 6 + a 7 = + + + + + + =. And the transportation tableau corresponding to the transshipment problem

More information

GENERAL ASSIGNMENT PROBLEM via Branch and Price JOHN AND LEI

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

Branch-Cut-and-Price solver for Vehicle Routing Problems

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

Selected Topics in Column Generation

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

Implementing Scalable Parallel Search Algorithms for Data-Intensive Applications

Implementing Scalable Parallel Search Algorithms for Data-Intensive Applications Implementing Scalable Parallel Search Algorithms for Data-Intensive Applications Ted Ralphs Industrial and Systems Engineering Lehigh University http://www.lehigh.edu/~tkr2 Laszlo Ladanyi IBM T.J. Watson

More information

Using Decomposition Techniques for Solving Large-Scale Capacitated Hub Location Problems with Single Assignment

Using Decomposition Techniques for Solving Large-Scale Capacitated Hub Location Problems with Single Assignment Using Decomposition Techniques for Solving Large-Scale Capacitated Hub Location Problems with Single Assignment Ivan Contreras*, Elena Fernández Department of Statistics and Operations Research Technical

More information

Branch-price-and-cut for vehicle routing. Guy Desaulniers

Branch-price-and-cut for vehicle routing. Guy Desaulniers Guy Desaulniers Professor, Polytechnique Montréal, Canada Director, GERAD, Canada VeRoLog PhD School 2018 Cagliari, Italy, June 2, 2018 Outline 1 VRPTW definition 2 Mathematical formulations Arc-flow formulation

More information

Column Generation and its applications

Column Generation and its applications Column Generation and its applications Murat Firat, dept. IE&IS, TU/e BPI Cluster meeting Outline Some real-life decision problems Standard formulations Basics of Column Generation Master formulations

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 Generic Separation Algorithm and Its Application to the Vehicle Routing Problem

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

The Size Robust Multiple Knapsack Problem

The Size Robust Multiple Knapsack Problem MASTER THESIS ICA-3251535 The Size Robust Multiple Knapsack Problem Branch and Price for the Separate and Combined Recovery Decomposition Model Author: D.D. Tönissen, Supervisors: dr. ir. J.M. van den

More information

A Recovery Algorithm for a Disrupted Airline Schedule

A Recovery Algorithm for a Disrupted Airline Schedule 1 A Recovery Algorithm for a Disrupted Airline Schedule Niklaus Eggenberg Matteo Salani and Prof. M. Bierlaire In collaboration with APM Technologies Index 2 Index Airline Scheduling in general The Disrupted

More information

Algorithms for Decision Support. Integer linear programming models

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

MVE165/MMG630, Applied Optimization Lecture 8 Integer linear programming algorithms. Ann-Brith Strömberg

MVE165/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 information

The MIP-Solving-Framework SCIP

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

Parallel Branch & Bound

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

56:272 Integer Programming & Network Flows Final Exam -- December 16, 1997

56:272 Integer Programming & Network Flows Final Exam -- December 16, 1997 56:272 Integer Programming & Network Flows Final Exam -- December 16, 1997 Answer #1 and any five of the remaining six problems! possible score 1. Multiple Choice 25 2. Traveling Salesman Problem 15 3.

More information

Integer Programming and Network Modeis

Integer Programming and Network Modeis H.A. Eiselt C.-L. Sandblom Integer Programming and Network Modeis With Contributions by K. Spielberg, E. Richards, B.T. Smith, G. Laporte, B.T. Boffey With 165 Figures and 43 Tables &m Springer CONTENTS

More information

Integer and Combinatorial Optimization

Integer and Combinatorial Optimization Integer and Combinatorial Optimization GEORGE NEMHAUSER School of Industrial and Systems Engineering Georgia Institute of Technology Atlanta, Georgia LAURENCE WOLSEY Center for Operations Research and

More information

Disjunctive cuts in branch-and-but-and-price algorithms Application to the capacitated vehicle routing problem

Disjunctive cuts in branch-and-but-and-price algorithms Application to the capacitated vehicle routing problem Disjunctive cuts in branch-and-but-and-price algorithms Application to the capacitated vehicle routing problem Stefan Ropke Technical University of Denmark, Department of Transport (DTU Transport) Column

More information

Column Generation Method for an Agent Scheduling Problem

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

Computational Integer Programming. Lecture 12: Branch and Cut. Dr. Ted Ralphs

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

7KH9HKLFOH5RXWLQJSUREOHP

7KH9HKLFOH5RXWLQJSUREOHP 7K9KO5RXWJSUREOP Given a set of vehicles with a certain capacity located at a depot and a set of customers with different demands at various locations, the vehicle routing problem (VRP) is how to satisfy

More information

Outline. Column Generation: Cutting Stock A very applied method. Introduction to Column Generation. Given an LP problem

Outline. 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 information

Column Generation: Cutting Stock

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

Fundamentals of Integer Programming

Fundamentals of Integer Programming Fundamentals of Integer Programming Di Yuan Department of Information Technology, Uppsala University January 2018 Outline Definition of integer programming Formulating some classical problems with integer

More information

A Parallel Macro Partitioning Framework for Solving Mixed Integer Programs

A Parallel Macro Partitioning Framework for Solving Mixed Integer Programs This research is funded by NSF, CMMI and CIEG 0521953: Exploiting Cyberinfrastructure to Solve Real-time Integer Programs A Parallel Macro Partitioning Framework for Solving Mixed Integer Programs Mahdi

More information

A Note on the Separation of Subtour Elimination Constraints in Asymmetric Routing Problems

A Note on the Separation of Subtour Elimination Constraints in Asymmetric Routing Problems Gutenberg School of Management and Economics Discussion Paper Series A Note on the Separation of Subtour Elimination Constraints in Asymmetric Routing Problems Michael Drexl March 202 Discussion paper

More information

Methods and Models for Combinatorial Optimization Exact methods for the Traveling Salesman Problem

Methods and Models for Combinatorial Optimization Exact methods for the Traveling Salesman Problem Methods and Models for Combinatorial Optimization Exact methods for the Traveling Salesman Problem L. De Giovanni M. Di Summa The Traveling Salesman Problem (TSP) is an optimization problem on a directed

More information

Implementing Custom Applications with CHiPPS. INFORMS Annual Conference, San Diego, CA October 12, 2009

Implementing Custom Applications with CHiPPS. INFORMS Annual Conference, San Diego, CA October 12, 2009 Implementing Custom TED RALPHS LEHIGH UNIVERSITY YAN XU SAS INSTITUTE LASZLO LADÁNYI IBM ILOG MATTHEW SALTZMAN CLEMSON UNIVERSITY INFORMS Annual Conference, San Diego, CA October 12, 2009 Thanks: Work

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

An Algorithm for a 3-Level Location-Routing Problem

An Algorithm for a 3-Level Location-Routing Problem An Algorithm for a 3-Level Location-Routing Problem Leonardo Souza Ribeiro Production Engineering Dept., COPPE/UFRJ Cidade Universitária, CT, Bloco F, Sala 03, RJ, Brasil Petrobras Avenida Almirante Barroso

More information

An Introduction to Dual Ascent Heuristics

An Introduction to Dual Ascent Heuristics An Introduction to Dual Ascent Heuristics Introduction A substantial proportion of Combinatorial Optimisation Problems (COPs) are essentially pure or mixed integer linear programming. COPs are in general

More information

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

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

Linear Programming. Course review MS-E2140. v. 1.1

Linear Programming. Course review MS-E2140. v. 1.1 Linear Programming MS-E2140 Course review v. 1.1 Course structure Modeling techniques Linear programming theory and the Simplex method Duality theory Dual Simplex algorithm and sensitivity analysis Integer

More information

MVE165/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 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 information

DM515 Spring 2011 Weekly Note 7

DM515 Spring 2011 Weekly Note 7 Institut for Matematik og Datalogi Syddansk Universitet May 18, 2011 JBJ DM515 Spring 2011 Weekly Note 7 Stuff covered in Week 20: MG sections 8.2-8,3 Overview of the course Hints for the exam Note that

More information

2 is not feasible if rounded. x =0,x 2

2 is not feasible if rounded. x =0,x 2 Integer Programming Definitions Pure Integer Programming all variables should be integers Mied integer Programming Some variables should be integers Binary integer programming The integer variables are

More information

sbb COIN-OR Simple Branch-and-Cut

sbb COIN-OR Simple Branch-and-Cut sbb COIN-OR Simple Branch-and-Cut John Forrest Lou Hafer lou@cs.sfu.ca CORS/INFORMS Joint Meeting Banff, May 2004 p.1/12 Outline Basic Branch-and-Bound Branching Node Selection Cuts and Heuristics Using

More information

Benders in a nutshell Matteo Fischetti, University of Padova

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

A Recovery Algorithm for a Disrupted Airline Schedule

A Recovery Algorithm for a Disrupted Airline Schedule 1 A Recovery Algorithm for a Disrupted Airline Schedule Niklaus Eggenberg Matteo Salani and Prof. M. Bierlaire In collaboration with APM Technologies Index 2 Index Airline Scheduling in general The Disrupted

More information

A Computational Study of Bi-directional Dynamic Programming for the Traveling Salesman Problem with Time Windows

A Computational Study of Bi-directional Dynamic Programming for the Traveling Salesman Problem with Time Windows A Computational Study of Bi-directional Dynamic Programming for the Traveling Salesman Problem with Time Windows Jing-Quan Li California PATH, University of California, Berkeley, Richmond, CA 94804, jingquan@path.berkeley.edu

More information

Mathematical Tools for Engineering and Management

Mathematical Tools for Engineering and Management Mathematical Tools for Engineering and Management Lecture 8 8 Dec 0 Overview Models, Data and Algorithms Linear Optimization Mathematical Background: Polyhedra, Simplex-Algorithm Sensitivity Analysis;

More information

Graphs and Network Flows IE411. Lecture 13. Dr. Ted Ralphs

Graphs and Network Flows IE411. Lecture 13. Dr. Ted Ralphs Graphs and Network Flows IE411 Lecture 13 Dr. Ted Ralphs IE411 Lecture 13 1 References for Today s Lecture IE411 Lecture 13 2 References for Today s Lecture Required reading Sections 21.1 21.2 References

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

Recursive column generation for the Tactical Berth Allocation Problem

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

Outline. 1 The matching problem. 2 The Chinese Postman Problem

Outline. 1 The matching problem. 2 The Chinese Postman Problem Outline The matching problem Maximum-cardinality matchings in bipartite graphs Maximum-cardinality matchings in bipartite graphs: The augmenting path algorithm 2 Let G = (V, E) be an undirected graph.

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

Final Exam Spring 2003

Final Exam Spring 2003 .8 Final Exam Spring Name Instructions.. Please answer all questions in the exam books that are provided.. Please budget your time carefully. It is often a good idea to read the entire exam first, so that

More information

February 19, Integer programming. Outline. Problem formulation. Branch-andbound

February 19, Integer programming. Outline. Problem formulation. Branch-andbound Olga Galinina olga.galinina@tut.fi ELT-53656 Network Analysis and Dimensioning II Department of Electronics and Communications Engineering Tampere University of Technology, Tampere, Finland February 19,

More information

Interior Point Stabilization for Column Generation

Interior Point Stabilization for Column Generation Interior Point Stabilization for Column Generation LOUIS-MARTIN ROUSSEAU, MICHEL GENDREAU CRT, Université de Montréal, C.P. 6128, succursale centre-ville, Montréal, Canada H3C 3J7 {louism,michelg}@crt.umontreal.ca

More information

Outline: Finish uncapacitated simplex method Negative cost cycle algorithm The max-flow problem Max-flow min-cut theorem

Outline: Finish uncapacitated simplex method Negative cost cycle algorithm The max-flow problem Max-flow min-cut theorem Outline: Finish uncapacitated simplex method Negative cost cycle algorithm The max-flow problem Max-flow min-cut theorem Uncapacitated Networks: Basic primal and dual solutions Flow conservation constraints

More information

Parallel Branch, Cut, and Price for Large-Scale Discrete Optimization

Parallel Branch, Cut, and Price for Large-Scale Discrete Optimization Parallel Branch, Cut, and Price for Large-Scale Discrete Optimization T. K. Ralphs L. Ladányi M. J. Saltzman April 4, 2003 Abstract In discrete optimization, most exact solution approaches are based on

More information

Implementing Scalable Parallel Search Algorithms for Data-intensive Applications

Implementing Scalable Parallel Search Algorithms for Data-intensive Applications Implementing Scalable Parallel Search Algorithms for Data-intensive Applications L. Ladányi 1, T. K. Ralphs 2, and M. J. Saltzman 3 1 Department of Mathematical Sciences, IBM T. J. Watson Research Center,

More information

Noncommercial Software for Mixed-Integer Linear Programming

Noncommercial Software for Mixed-Integer Linear Programming Noncommercial Software for Mixed-Integer Linear Programming J. T. Linderoth T. K. Ralphs December 23, 2004 Abstract We present an overview of noncommercial software tools for the solution of mixed-integer

More information

Martin W.P. Savelsbergh. Georgia Institute of Technology. School of Industrial and Systems Engineering

Martin W.P. Savelsbergh. Georgia Institute of Technology. School of Industrial and Systems Engineering A MINTO short course 1 Martin W.P. Savelsbergh George L. Nemhauser Introduction Georgia Institute of Technology School of Industrial and Systems Engineering Atlanta, GA 30332-0205 (February 1, 1995) MINTO

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

Integer Programming. Xi Chen. Department of Management Science and Engineering International Business School Beijing Foreign Studies University

Integer Programming. Xi Chen. Department of Management Science and Engineering International Business School Beijing Foreign Studies University Integer Programming Xi Chen Department of Management Science and Engineering International Business School Beijing Foreign Studies University Xi Chen (chenxi0109@bfsu.edu.cn) Integer Programming 1 / 42

More information

Solving small VRPTWs with Constraint Programming Based Column Generation

Solving small VRPTWs with Constraint Programming Based Column Generation Proceedings CPAIOR 02 Solving small VRPTWs with Constraint Programming Based Column Generation Louis-Martin Rousseau, Michel Gendreau, Gilles Pesant Center for Research on Transportation Université de

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

Computational Experience with Parallel Integer Programming using the CHiPPS Framework. INFORMS Annual Conference, San Diego, CA October 12, 2009

Computational Experience with Parallel Integer Programming using the CHiPPS Framework. INFORMS Annual Conference, San Diego, CA October 12, 2009 Computational Experience with Parallel Integer Programming using the CHiPPS Framework TED RALPHS LEHIGH UNIVERSITY YAN XU SAS INSTITUTE LASZLO LADÁNYI IBM ILOG MATTHEW SALTZMAN CLEMSON UNIVERSITY INFORMS

More information

Outline. Modeling. Outline DMP204 SCHEDULING, TIMETABLING AND ROUTING. 1. Models Lecture 5 Mixed Integer Programming Models and Exercises

Outline. 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 information

to the Traveling Salesman Problem 1 Susanne Timsj Applied Optimization and Modeling Group (TOM) Department of Mathematics and Physics

to the Traveling Salesman Problem 1 Susanne Timsj Applied Optimization and Modeling Group (TOM) Department of Mathematics and Physics An Application of Lagrangian Relaxation to the Traveling Salesman Problem 1 Susanne Timsj Applied Optimization and Modeling Group (TOM) Department of Mathematics and Physics M lardalen University SE-721

More information

Department of Mathematics Oleg Burdakov of 30 October Consider the following linear programming problem (LP):

Department of Mathematics Oleg Burdakov of 30 October Consider the following linear programming problem (LP): Linköping University Optimization TAOP3(0) Department of Mathematics Examination Oleg Burdakov of 30 October 03 Assignment Consider the following linear programming problem (LP): max z = x + x s.t. x x

More information

COLUMN GENERATION. Diego Klabjan, Northwestern University

COLUMN GENERATION. Diego Klabjan, Northwestern University PARALLEL PRIMAL-DUAL WITH COLUMN GENERATION Diego Klabjan, Northwestern University Primal-Dual Algorithm Started with Dantzig, Ford, and Fulkerson in 1956. Primal-dual algorithms Primal step: Solve a primal

More information

LaGO - A solver for mixed integer nonlinear programming

LaGO - A solver for mixed integer nonlinear programming LaGO - A solver for mixed integer nonlinear programming Ivo Nowak June 1 2005 Problem formulation MINLP: min f(x, y) s.t. g(x, y) 0 h(x, y) = 0 x [x, x] y [y, y] integer MINLP: - n

More information

A Bi-directional Resource-bounded Dynamic Programming Approach for the Traveling Salesman Problem with Time Windows

A Bi-directional Resource-bounded Dynamic Programming Approach for the Traveling Salesman Problem with Time Windows Submitted manuscript A Bi-directional Resource-bounded Dynamic Programming Approach for the Traveling Salesman Problem with Time Windows Jing-Quan Li California PATH, University of California, Berkeley,

More information

Two models of the capacitated vehicle routing problem

Two models of the capacitated vehicle routing problem Croatian Operational Research Review 463 CRORR 8(2017), 463 469 Two models of the capacitated vehicle routing problem Zuzana Borčinová 1, 1 Faculty of Management Science and Informatics, University of

More information

Branch-and-Price-and-Cut for the Vehicle Routing and Truck Driver Scheduling Problem

Branch-and-Price-and-Cut for the Vehicle Routing and Truck Driver Scheduling Problem Branch-and-Price-and-Cut for the Vehicle Routing and Truck Driver Scheduling Problem Christian Tilk,a a Chair of Logistics Management, Johannes Gutenberg University Mainz, Jakob-Welder-Weg 9, D-5512 Mainz,

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

Building a Custom Solver with the COIN-OR Branch, Cut, and Price Frameworks

Building a Custom Solver with the COIN-OR Branch, Cut, and Price Frameworks Building a Custom Solver with the COIN-OR Branch, Cut, and Price Frameworks Ted Ralphs and Menal Guzelsoy Lehigh University László Ladányi IBM T. J. Watson Research Center Matthew Saltzman Clemson University

More information

LEAST COST ROUTING ALGORITHM WITH THE STATE SPACE RELAXATION IN A CENTRALIZED NETWORK

LEAST COST ROUTING ALGORITHM WITH THE STATE SPACE RELAXATION IN A CENTRALIZED NETWORK VOL., NO., JUNE 08 ISSN 896608 00608 Asian Research Publishing Network (ARPN). All rights reserved. LEAST COST ROUTING ALGORITHM WITH THE STATE SPACE RELAXATION IN A CENTRALIZED NETWORK Y. J. Lee Department

More information

Doing It in Parallel (DIP) with The COIN-OR High-Performance Parallel Search Framework (CHiPPS)

Doing It in Parallel (DIP) with The COIN-OR High-Performance Parallel Search Framework (CHiPPS) Doing It in Parallel (DIP) with The COIN-OR High-Performance Parallel Search Framework (CHiPPS) TED RALPHS LEHIGH UNIVERSITY YAN XU SAS INSTITUTE LASZLO LADÁNYI IBM T.J. WATSON RESEARCH CENTER MATTHEW

More information

Advanced Operations Research Prof. G. Srinivasan Department of Management Studies Indian Institute of Technology, Madras

Advanced Operations Research Prof. G. Srinivasan Department of Management Studies Indian Institute of Technology, Madras Advanced Operations Research Prof. G. Srinivasan Department of Management Studies Indian Institute of Technology, Madras Lecture - 35 Quadratic Programming In this lecture, we continue our discussion on

More information

6 ROUTING PROBLEMS VEHICLE ROUTING PROBLEMS. Vehicle Routing Problem, VRP:

6 ROUTING PROBLEMS VEHICLE ROUTING PROBLEMS. Vehicle Routing Problem, VRP: 6 ROUTING PROBLEMS VEHICLE ROUTING PROBLEMS Vehicle Routing Problem, VRP: Customers i=1,...,n with demands of a product must be served using a fleet of vehicles for the deliveries. The vehicles, with given

More information

A. Atamturk. G.L. Nemhauser. M.W.P. Savelsbergh. Georgia Institute of Technology. School of Industrial and Systems Engineering.

A. Atamturk. G.L. Nemhauser. M.W.P. Savelsbergh. Georgia Institute of Technology. School of Industrial and Systems Engineering. A Combined Lagrangian, Linear Programming and Implication Heuristic for Large-Scale Set Partitioning Problems 1 A. Atamturk G.L. Nemhauser M.W.P. Savelsbergh Georgia Institute of Technology School of Industrial

More information

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

EXERCISES SHORTEST PATHS: APPLICATIONS, OPTIMIZATION, VARIATIONS, AND SOLVING THE CONSTRAINED SHORTEST PATH PROBLEM. 1 Applications and Modelling SHORTEST PATHS: APPLICATIONS, OPTIMIZATION, VARIATIONS, AND SOLVING THE CONSTRAINED SHORTEST PATH PROBLEM EXERCISES Prepared by Natashia Boland 1 and Irina Dumitrescu 2 1 Applications and Modelling 1.1

More information

Search Algorithms. IE 496 Lecture 17

Search Algorithms. IE 496 Lecture 17 Search Algorithms IE 496 Lecture 17 Reading for This Lecture Primary Horowitz and Sahni, Chapter 8 Basic Search Algorithms Search Algorithms Search algorithms are fundamental techniques applied to solve

More information

Recoverable Robust Optimization for (some) combinatorial problems

Recoverable Robust Optimization for (some) combinatorial problems Recoverable Robust Optimization for (some) combinatorial problems Marjan van den Akker joint work with 1 Han Hoogeveen, Paul Bouman, Denise Tönissen, Judith Stoef Recoverable robustness Method to cope

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

Graph Coloring via Constraint Programming-based Column Generation

Graph 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

Accelerating Column Generation for a Survivable Network Design Problem

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

Algorithms for the Bin Packing Problem with Conflicts

Algorithms for the Bin Packing Problem with Conflicts Algorithms for the Bin Packing Problem with Conflicts Albert E. Fernandes Muritiba *, Manuel Iori, Enrico Malaguti*, Paolo Toth* *Dipartimento di Elettronica, Informatica e Sistemistica, Università degli

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 Menal Guzelsoy and Ted Ralphs Industrial and Systems Engineering Lehigh University INFORMS Annual Meeting, San Francisco, CA, November

More information

The Gurobi Optimizer. Bob Bixby

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

Decomposition in Integer Linear Programming

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

Decomposition in Integer Linear Programming

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

DECOMPOSITION METHODS FOR INTEGER LINEAR PROGRAMMING

DECOMPOSITION METHODS FOR INTEGER LINEAR PROGRAMMING DECOMPOSITION METHODS FOR INTEGER LINEAR PROGRAMMING by Matthew Galati Presented to the Graduate and Research Committee of Lehigh University in Candidacy for the Degree of Doctor of Philosophy in Industrial

More information

Algorithms for Graph Visualization Layered Layout

Algorithms for Graph Visualization Layered Layout Algorithms for Graph Visualization Layered Layout INSTITUT FÜR THEORETISCHE INFORMATIK FAKULTÄT FÜR INFORMATIK Tamara Mchedlidze 13.12.2017 1 Dr. Tamara Mchedlidze Algorithmen zur Visualisierung von Graphen

More information