A Compact Linearisation of Euclidean Single Allocation Hub Location Problems
|
|
- Sara Rogers
- 6 years ago
- Views:
Transcription
1 A Compact Linearisation of Euclidean Single Allocation Hub Location Problems J. Fabian Meier 1,2, Uwe Clausen 1 Institute of Transport Logistics, TU Dortmund, Germany Borzou Rostami 1, Christoph Buchheim 1 Fakultät für Mathematik, TU Dortmund, Germany Abstract Hub location problems are strategic network planning problems. They formalise the challenge of mutually exchanging shipments between a large set of depots. The aim is to choose a set of hubs (out of a given set of possible hubs) and connect every depot to a hub so that the total transport costs for exchanging shipments between the depots are minimised. In classical hub location problems, the unit cost for transport between hubs is proportional to the distance between the hubs. Often these distances are Euclidean distances: Then is possible to replace the quadratic cost term for hub-hub-transport in the objective function by a linear term and a set of linear inequalities. The resulting model can be solved by a row generation scheme. The strength of the method is shown by solving all AP instances to optimality. Keywords: Hub Location, Euclidean Distance, Row Generation 1 This research has been funded by the German Research Foundation (DFG) within the project Lenkung des Güterflusses in durch Gateways gekoppelten Logistik-Service- Netzwerken mittels quadratischer Optimierung (CL 318/14 and BU 2313/2) 2 meier@itl.tu-dortmund.de
2 1 Introduction O Kelly [5] introduced the Uncapacitated Single Allocation p-hub Median Problem (USApHMP) in 1987: A set of p hubs is chosen from n possible hub locations and every depot is connected to exactly one hub. The depots mutually exchange shipments: A shipment is sent from the source depot to the assigned hub, then to the assigned hub of the sink depot and then finally to the sink depot. The aim is to choose hubs and assignments with minimal overall transport costs. We state the problem in the original quadratic form: Consider a complete graph G = (V, E) with V = n. Each node i V of the graph corresponds to origins, destinations and possible hub locations. Let C ij be the transport cost per unit of flow from node i to node j, and W ij be the amount of flow from node i to node j (the shipment from i to j). The cost per unit of flow for each path i k l j from an origin node i to a destination node j which passes hubs k and l respectively, is β 1 C ik + αc kl + β 2 C lj, where β 1, α, and β 2 are the collection, transfer and distribution costs respectively. We define the binary variable x ik to indicate the allocation of node i to the hub located at node k. If node i is assigned to itself, then node i is a hub. To ease the argumentation in the following sections, we define d kl = αc kl K ik = β 1 C ik W ij + β 2 C ki W ji. The quadratic 0 1 formulation of the USApHMP is stated as follows: min W ij d kl x ik x jl i V k V K ik x ik + i V j V j V k V l V j V s.t. x ik = 1 i V (1) k V x ik x kk i, k V (2) x kk = p (3) k V x ik {0, 1} i, k V. (4) Constraints (1) indicate that node i is allocated to precisely one hub node. The inequalities (2) make sure that a node i can only be allocated to a hub node. Constraints (3) force the number of selected hubs to be p.
3 The main difficulty of the problem lies in the quadratic structure of the objective function. Different attempts have been made to linearise the objective function. The MILP formulations of Skorin-Kapov et al. [6] and Ernst and Krishnamoorthy [2] result in O ( n 3) or O ( n 4) variables. Other authors successfully applied metaheuristic techniques to construct good solutions for large instances, recently even for more than 100 nodes [3,4]. In many problem instances (also in those from the usual benchmark instances CAB and AP [1]) d kl is given by the Euclidean distance of the hubs (or a constant multiple of it). Assuming this kind of structure, we can construct a linearisation with only O ( n 2) variables, but O ( n 4) additional constraints. The resulting linearized problem can be solved very efficiently in practice by a row generation procedure. The next section will explain our new approach to the linearisation of the quadratic term in the objective function, while Sect. 3 gives computational evidence. Section 4 gives an outlook for further applications. 2 The Euclidean Hub Location Problem To reduce the number of linearization variables to O ( n 2), we will make the assumption that all distances d kl are calculated by a vector space norm derived from an inner product,, i.e. every hub k is represented by a vector w k in a Hilbert space W and we have d kl = w k w l. Although it is nearly always possible to realise this in an n 1 dimensional vector space, we are more interested in the situation in which costs are proportional to the Euclidean distance in a 2-dimensional landscape. We call this problem Euclidean Hub Location Problem. The Euclidean distance assumption is only relevant for the distance between possible hubs, i.e. the coefficient of x ik x jl in the objective function of (USApHMP). Let us introduce the variables y ij k V W ij d kl x ik x jl (5) l V which represent the total cost of any hub-hub flow induced by the shipment (i, j). The quadratic term of the objective function of (USApHMP) is then equal to i V j V y ij. As the equation (5) is quadratic itself, it is not sensible to add it directly to the problem. Instead, we want to derive a set of linear inequalities which forces y ij to fulfil the constraint (5). In W, we can define the orthogonal projection P w (u) of a vector u W to another vector w W. We know that P w (u 1 ) P w (u 2 ) u 1 u 2 for all
4 u 1, u 2 W. Furthermore, P w (u) = λ u w w w with λ u w = u, w w R (6) and P w (u 1 ) P w (u 2 ) = λ u 1 w λ u 2 w. (7) The value λ u w can be understood as co-ordinate in direction w. We see that it can be directly computed from the knowledge of the involved vectors. This implies, that for any k, l V and w W, w 0, we have: d kl = w k w l P w (w k ) P w (w l ) λ w k w If w = w k w l we have equality everywhere. This implies that y ij = W ij d kl x ik x jl k V l V ( ) W ij λ w k w λ w l w xik x jl k V l V = W ij λ w k w x ik x jl W ij λ w l w x ik x jl k V l V k V l V where we can reduce the sums due to (1): λ w l w (8) = k V W ij λ w k w x ik l V W ij λ w l w x jl (9) Let λ k mh be the co-ordinate associated with the projection of w k to w m w h (see also Fig. 1). Then we know that y ij k V W ij λ k mh x ik l V W ij λ l mh x jl i, j, m, h V (I mh ij ) If x im = x jh = 1, then (Iij mh ) is an equality. This means that the constraints (Iij mh ) completely determine y ij if the values of x ik are binary and chosen according to the constraint (1), so that we can insert y ij into the objective function and get an equivalent MIP to USApHMP which we call EUSApHMP. The LP formulation of EUSApHMP may be weaker than in the original formulation. The number of variables of
5 EUSApHMP is O ( n 2). Unfortunately, the number of inequalities resulting from our construction is O ( n 4). Although we can halve the number of variables and inequalities by exploiting the symmetry in the indices of y ij, it is still not sensible to add all these inequalities at once. Hence we use a row generation scheme and add the inequalities step by step. From a theoretical point of view, an efficient separation algorithm is given trivially, as the number of constraints to consider is polynomial. However, from a practical point of view, checking O ( n 4) inequalities in every call of a separation algorithm is still too expensive. We thus decided to adopt the following heuristic separation method in order to limit both iterations and added inequalities, splitting the solution process into two rounds: In the first round, we iteratively solve the LP relaxation and determine a solution for USApHMP by rounding. We start by considering only the inequalities ( ) I ij ij for all i j V. After every iteration, we add those inequalities that are violated by the rounded solution, i.e. those ( ) Iij kl for ˆx ik = ˆx jl = 1, where {ˆx ik } ik is the solution of the USApHMP. In this way, we restrict the number of new inequalities to at most n 2. We stop the first round when there are no more inequalities to add. Note that this method does not guarantee that the final LP solution solves the LP relaxation of EUSApHMP; but the rounded solutions can be used as initial values for the second round. In the second round, we use a MIP solver. We start by adding all inequalities we used so far in the LP relaxation. Then we proceed as before: In every step we determine a MIP solution {ˆx ik } ik for USApHMP from the given solution and add all inequalities for ˆx ik = ˆx jl = 1. If we have no more inequalities to add, we stop, because this proves the solution is optimal. As shown by the computational results presented in the following section, this approach is able to effectively limit the number of added inequalities while at the same time leading to a very small number of iterations Fig. 1. Illustration of the orthogonal projection of all hub positions onto the vector w 5 w 4. The distance of any two vectors exceeds the difference of the λ values.
6 3 Computational Results We code the procedures in C# using the Gurobi 5.6 solver for the LPs and MIPs. We used a 3.4 GHz computer with 16GB RAM and four threads. The famous AP data served as our test data set [1]; we only give results for instances starting at 50 nodes. We see in Table 1 and 2 that all instances are solved to optimality. After stating the size of V and the number of hubs, we give two specific iterations of the calculation: (i) The first iteration which reaches less than 1% gap between the best found MIP solution and the last lower bound. (ii) The iteration after which optimality is proven. Iterations are written as L/M where L and M are the iterations of the first and second round. For each iteration we also give the total time in seconds and the total number of added inequalities. Table 1 AP results up to 90 nodes. n p Gap 1% Optimality Optimal Value It s (Iij mh ) It s (Iij mh ) / / / / / / / / / / / / / / / / / / / /
7 Table 2 AP results for more than 90 nodes. No published optimal solutions for these instances could be found. n p Gap 1% Optimality Optimal Value It s (Iij mh ) It s (Iij mh ) / / / / / / / / / / / / / / / / / / / / / / / / / / / / / / / / / / / / /
8 We see that at most 4n 2 of the possible n 4 /2 inequalities of (Iij mh )-type are used. Often, rounding one of the first LP iterations gives the optimal value, so that the optimal value is reached long before optimality is proven. In [3] we find the best known values for 100 and 200 nodes, found by a metaheuristic method: Our optimal results are identical to these best known results, proving them to be optimal. 4 Conclusion Constructing linearisations from the Euclidean structure that is present in many instances of hub location problems leads to compact formulations with few variables. These can be solved even for a large number of nodes. As future work, we plan to apply this method to even larger instances and other variants of single allocation problems, like capacitated or stochastic problems. Furthermore, the method can be extended to the case where the distances are not Euclidean but just fulfil the triangle inequality, i.e., the metric case. Setting λ k mh = d kh leads to the same kind of model, but turned out to be weaker in preliminary numerical experiments. Nevertheless it is worth exploring these techniques for more general cases. References [1] Beasley, J. E., OR library (2012). URL [2] Ernst, A. T. and M. Krishnamoorthy, Efficient algorithms for the uncapacitated single allocation p-hub median problem, Location science 4 (1996), pp [3] Ilić, A., D. Urošević, J. Brimberg and N. Mladenović, A general variable neighborhood search for solving the uncapacitated single allocation p-hub median problem, European Journal of Operational Research 206 (2010), pp [4] Kratica, J., Z. Stanimirović, D. Tošić and V. Filipović, Two genetic algorithms for solving the uncapacitated single allocation p-hub median problem, European Journal of Operational Research 182 (2007), pp [5] O Kelly, M. E., A quadratic integer program for the location of interacting hub facilities, European Journal of Operational Research 32 (1987), pp [6] Skorin-Kapov, D., J. Skorin-Kapov and M. O Kelly, Tight linear programming relaxations of uncapacitated p-hub median problems, European Journal of Operational Research 94 (1996), pp
The Uncapacitated Single Allocation p-hub Median Problem with Stepwise Cost Function
The Uncapacitated Single Allocation p-hub Median Problem with Stepwise Cost Function Borzou Rostami a, Christoph Buchheim b a Department of Mathematics and Industrial Engineering, Polytechnique Montreal
More informationSolving Classical and New Single Allocation Hub Location Problems on Euclidean Data
Solving Classical and New Single Allocation Hub Location Problems on Euclidean Data J. Fabian Meier a, Uwe Clausen a a Institute of Transport Logistics, TU Dortmund, Leonhard-Euler-Str. 2, 44227 Dortmund,
More informationReliable single allocation hub location problem under hub breakdowns
Reliable single allocation hub location problem under hub breakdowns Borzou Rostami 1, Nicolas Kämmerling 2, Christoph Buchheim 3, and Uwe Clausen 2 1 École de Technologie Supérieure de Montréal and Interuniversity
More informationA Hub Location Problem with Fully Interconnected Backbone and Access Networks
A Hub Location Problem with Fully Interconnected Backbone and Access Networks Tommy Thomadsen Informatics and Mathematical Modelling Technical University of Denmark 2800 Kgs. Lyngby Denmark tt@imm.dtu.dk
More informationSome Numerical Studies for a Complicated Hub Location Problem
Some Numerical Studies for a Complicated Hub Location Problem J. Fabian Meier and Uwe Clausen Abstract We consider a complicated hub location problem which includes multiallocation, different hub sizes
More informationA LAGRANGIAN RELAXATION FOR CAPACITATED SINGLE ALLOCATION P-HUB MEDIAN PROBLEM WITH MULTIPLE CAPACITY LEVELS
A LAGRANGIAN RELAXATION FOR CAPACITATED SINGLE ALLOCATION P-HUB MEDIAN PROBLEM WITH MULTIPLE CAPACITY LEVELS Ching-Jung Ting Department of Industrial Engineering and Management, Yuan Ze University Kuo-Rui
More informationOptimization Exercise Set n.5 :
Optimization Exercise Set n.5 : Prepared by S. Coniglio translated by O. Jabali 2016/2017 1 5.1 Airport location In air transportation, usually there is not a direct connection between every pair of airports.
More informationLinear Programming-Based Algorithms for the Fixed-Hub Single Allocation Problem
Linear Programming-Based Algorithms for the Fixed-Hub Single Allocation Problem Dongdong Ge, Yinyu Ye, Jiawei Zhang March 20, 2007 Abstract This paper discusses the fixed-hub single allocation problem.
More informationOptimization Exercise Set n. 4 :
Optimization Exercise Set n. 4 : Prepared by S. Coniglio and E. Amaldi translated by O. Jabali 2018/2019 1 4.1 Airport location In air transportation, usually there is not a direct connection between every
More informationConic optimization under combinatorial sparsity constraints
Conic optimization under combinatorial sparsity constraints Christoph Buchheim and Emiliano Traversi Abstract We present a heuristic approach for conic optimization problems containing sparsity constraints.
More informationA heuristic algorithm for the Aircraft Landing Problem
22nd International Congress on Modelling and Simulation, Hobart, Tasmania, Australia, 3 to 8 December 2017 mssanz.org.au/modsim2017 A heuristic algorithm for the Aircraft Landing Problem Amir Salehipour
More informationA Mixed-Integer Linear Program for the Traveling Salesman Problem with Structured Time Windows
A Mixed-Integer Linear Program for the Traveling Salesman Problem with Structured Time Windows Philipp Hungerländer Christian Truden 5th January 2017 Abstract In this extended abstract we introduce the
More informationMETAHEURISTICS FOR HUB LOCATION MODELS
Clemson University TigerPrints All Dissertations Dissertations 8-2011 METAHEURISTICS FOR HUB LOCATION MODELS Ornurai Sangsawang Clemson University, osangsawang@yahoo.com Follow this and additional works
More informationDiscrete PSO for the Uncapacitated Single Allocation Hub Location Problem
Brock University Department of Computer Science Discrete PSO for the Uncapacitated Single Allocation Hub Location Problem Alexander Bailey, Beatrice Ombuki-Berman, and Stephen Asobiela Technical Report
More informationThe tree-of-hubs location problem: A comparison of formulations
1 The tree-of-hubs location problem: A comparison of formulations Iván Contreras, Elena Fernández Departament d Estadística i Investigació Operativa. Universitat Politècnica de Catalunya, Barcelona, Spain.
More informationStabilized Branch-and-cut-and-price for the Generalized Assignment Problem
Stabilized Branch-and-cut-and-price for the Generalized Assignment Problem Alexandre Pigatti, Marcus Poggi de Aragão Departamento de Informática, PUC do Rio de Janeiro {apigatti, poggi}@inf.puc-rio.br
More informationA Capacity Scaling Procedure for the Multi-Commodity Capacitated Network Design Problem. Ryutsu Keizai University Naoto KATAYAMA
A Capacity Scaling Procedure for the Multi-Commodity Capacitated Network Design Problem Ryutsu Keizai University Naoto KATAYAMA Problems 2006 1 Multi-Commodity Network Design Problem The basic model for
More informationBelo Horizonte, Minas Gerais, Brasil. Universidade Federal de Minas Gerais Belo Horizonte, Minas Gerais, Brasil
Detailed computational results for the paper entitled The k-cardinality Tree Problem: Reformulations and Lagrangian Relaxation, under revision for Discrete Applied Mathematics Frederico P. Quintão a,b
More informationA New Compact Formulation for Discrete p-dispersion
Gutenberg School of Management and Economics & Research Unit Interdisciplinary Public Policy Discussion Paper Series A New Compact Formulation for Discrete p-dispersion David Sayah and Stefan Irnich November
More informationInderjit Dhillon The University of Texas at Austin
Inderjit Dhillon The University of Texas at Austin ( Universidad Carlos III de Madrid; 15 th June, 2012) (Based on joint work with J. Brickell, S. Sra, J. Tropp) Introduction 2 / 29 Notion of distance
More informationEasy and not so easy multifacility location problems... (In 20 minutes.)
Easy and not so easy multifacility location problems... (In 20 minutes.) MINLP 2014 Pittsburgh, June 2014 Justo Puerto Institute of Mathematics (IMUS) Universidad de Sevilla Outline 1 Introduction (In
More informationWeek 4. (1) 0 f ij u ij.
Week 4 1 Network Flow Chapter 7 of the book is about optimisation problems on networks. Section 7.1 gives a quick introduction to the definitions of graph theory. In fact I hope these are already known
More informationIntroduction to integer programming III:
Introduction to integer programming III: Network Flow, Interval Scheduling, and Vehicle Routing Problems Martin Branda Charles University in Prague Faculty of Mathematics and Physics Department of Probability
More informationAn Adaptive Partition-based Approach for Solving Two-stage Stochastic Programs with Fixed Recourse
An Adaptive Partition-based Approach for Solving Two-stage Stochastic Programs with Fixed Recourse Yongjia Song, James Luedtke Virginia Commonwealth University, Richmond, VA, ysong3@vcu.edu University
More informationON THE HUB-AND-SPOKE MODEL WITH ARC CAPACITY CONATRAINTS
Journal of the Operations Research Society of Japan 2003, Vol. 46, No. 4, 409-428 2003 The Operations Research Society of Japan ON THE HUB-AND-SPOKE MODEL WITH ARC CAPACITY CONATRAINTS Mihiro Sasaki Nanzan
More informationIntroduction into Vehicle Routing Problems and other basic mixed-integer problems
Introduction into Vehicle Routing Problems and other basic mixed-integer problems Martin Branda Charles University in Prague Faculty of Mathematics and Physics Department of Probability and Mathematical
More informationTotally unimodular matrices. Introduction to integer programming III: Network Flow, Interval Scheduling, and Vehicle Routing Problems
Totally unimodular matrices Introduction to integer programming III: Network Flow, Interval Scheduling, and Vehicle Routing Problems Martin Branda Charles University in Prague Faculty of Mathematics and
More informationThe Multiple Checkpoint Ordering Problem
The Multiple Checkpoint Ordering Problem Philipp Hungerländer Kerstin Maier November 19, 2017 Abstract The multiple Checkpoint Ordering Problem (mcop) aims to find an optimal arrangement of n one-dimensional
More informationParallel PIPS-SBB Multi-level parallelism for 2-stage SMIPS. Lluís-Miquel Munguia, Geoffrey M. Oxberry, Deepak Rajan, Yuji Shinano
Parallel PIPS-SBB Multi-level parallelism for 2-stage SMIPS Lluís-Miquel Munguia, Geoffrey M. Oxberry, Deepak Rajan, Yuji Shinano ... Our contribution PIPS-PSBB*: Multi-level parallelism for Stochastic
More informationAn Optimization-Based Heuristic for the Split Delivery Vehicle Routing Problem
An Optimization-Based Heuristic for the Split Delivery Vehicle Routing Problem Claudia Archetti (1) Martin W.P. Savelsbergh (2) M. Grazia Speranza (1) (1) University of Brescia, Department of Quantitative
More informationThe Multiple Traveling Salesperson Problem on Regular Grids
Philipp Hungerländer Anna Jellen Stefan Jessenitschnig Lisa Knoblinger Manuel Lackenbucher Kerstin Maier September 10, 2018 Abstract In this work we analyze the multiple Traveling Salesperson Problem (mtsp)
More informationA maritime version of the Travelling Salesman Problem
A maritime version of the Travelling Salesman Problem Enrico Malaguti, Silvano Martello, Alberto Santini May 31, 2015 Plan 1 The Capacitated TSP with Pickup and Delivery 2 The TSPPD with Draught Limits
More information3.7 Cutting plane methods
3.7 Cutting plane methods Generic ILP problem min{ c t x : x X = {x Z n + : Ax b} } with m n matrix A and n 1 vector b of rationals. According to Meyer s theorem: There exists an ideal formulation: conv(x
More informationDesigning Survivable Networks: A Flow Based Approach
Designing Survivable Networks: A Flow Based Approach Prakash Mirchandani 1 University of Pittsburgh This is joint work with Anant Balakrishnan 2 of the University of Texas at Austin and Hari Natarajan
More informationModels and Cuts for the Two-Echelon Vehicle Routing Problem
Models and Cuts for the Two-Echelon Vehicle Routing Problem Guido Perboli Roberto Tadei Francesco Masoero Department of Control and Computer Engineering, Politecnico di Torino Corso Duca degli Abruzzi,
More informationNotes on Dantzig-Wolfe decomposition and column generation
Notes on Dantzig-Wolfe decomposition and column generation Mette Gamst November 11, 2010 1 Introduction This note introduces an exact solution method for mathematical programming problems. The method is
More information- Well-characterized problems, min-max relations, approximate certificates. - LP problems in the standard form, primal and dual linear programs
LP-Duality ( Approximation Algorithms by V. Vazirani, Chapter 12) - Well-characterized problems, min-max relations, approximate certificates - LP problems in the standard form, primal and dual linear programs
More informationComputational testing of exact separation for mixed-integer knapsack problems
Computational testing of exact separation for mixed-integer knapsack problems Pasquale Avella (joint work with Maurizio Boccia and Igor Vasiliev ) DING - Università del Sannio Russian Academy of Sciences
More informationAn Exact Algorithm for the Steiner Tree Problem with Delays
Electronic Notes in Discrete Mathematics 36 (2010) 223 230 www.elsevier.com/locate/endm An Exact Algorithm for the Steiner Tree Problem with Delays Valeria Leggieri 1 Dipartimento di Matematica, Università
More informationReconnect 04 Introduction to Integer Programming
Sandia is a multiprogram laboratory operated by Sandia Corporation, a Lockheed Martin Company, Reconnect 04 Introduction to Integer Programming Cynthia Phillips, Sandia National Laboratories Integer programming
More informationXLVI Pesquisa Operacional na Gestão da Segurança Pública
A linear formulation with O(n 2 ) variables for the quadratic assignment problem Serigne Gueye and Philippe Michelon Université d Avignon et des Pays du Vaucluse, Laboratoire d Informatique d Avignon (LIA),
More informationBenders Decomposition for the Uncapacitated Multicommodity Network Design Problem
Benders Decomposition for the Uncapacitated Multicommodity Network Design Problem 1 Carlos Armando Zetina, 1 Ivan Contreras, 2 Jean-François Cordeau 1 Concordia University and CIRRELT, Montréal, Canada
More informationin specific rooms. An objective is to minimise the geographical distance between the allocated rooms and the requested rooms of the events. A highly p
Room Allocation Optimisation at the Technical University of Denmark Niels-Christian Fink Bagger Jesper Larsen Thomas Stidsen Keywords University Course Timetabling Room Allocation Mathematical Programming
More informationInstitut für Numerische und Angewandte Mathematik
Institut für Numerische und Angewandte Mathematik When closest is not always the best: The distributed p-median problem Brimberg, J.; Schöbel, A. Nr. 1 Preprint-Serie des Instituts für Numerische und Angewandte
More informationStochastic Hub and Spoke Networks
Stochastic Hub and Spoke Networks Edward Eric Hult Homerton College Judge Business School University of Cambridge A thesis submitted for the degree of Doctor of Philosophy 2011 I would like to dedicate
More informationTHE LATEST ARRIVAL HUB LOCATION PROBLEM FOR CARGO DELIVERY SYSTEMS WITH STOPOVERS. Hande YAMAN Bahar Y. KARA Barbaros Ç. TANSEL.
THE LATEST ARRIVAL HUB LOCATION PROBLEM FOR CARGO DELIVERY SYSTEMS WITH STOPOVERS Hande YAMAN Bahar Y. KARA Barbaros Ç. TANSEL July, 2006 Department of Industrial Engineering, Bilkent University, Bilkent
More informationLagrangean Decomposition for Mean-Variance Combinatorial Optimization
Lagrangean Decomposition for Mean-Variance Combinatorial Optimization Frank Baumann, Christoph Buchheim, and Anna Ilyina Fakultät für Mathematik, Technische Universität Dortmund, Germany {frank.baumann,christoph.buchheim,anna.ilyina}@tu-dortmund.de
More informationDS-GA 1002 Lecture notes 0 Fall Linear Algebra. These notes provide a review of basic concepts in linear algebra.
DS-GA 1002 Lecture notes 0 Fall 2016 Linear Algebra These notes provide a review of basic concepts in linear algebra. 1 Vector spaces You are no doubt familiar with vectors in R 2 or R 3, i.e. [ ] 1.1
More informationLOWER BOUNDS FOR THE UNCAPACITATED FACILITY LOCATION PROBLEM WITH USER PREFERENCES. 1 Introduction
LOWER BOUNDS FOR THE UNCAPACITATED FACILITY LOCATION PROBLEM WITH USER PREFERENCES PIERRE HANSEN, YURI KOCHETOV 2, NENAD MLADENOVIĆ,3 GERAD and Department of Quantitative Methods in Management, HEC Montréal,
More informationDiscrete (and Continuous) Optimization WI4 131
Discrete (and Continuous) Optimization WI4 131 Kees Roos Technische Universiteit Delft Faculteit Electrotechniek, Wiskunde en Informatica Afdeling Informatie, Systemen en Algoritmiek e-mail: C.Roos@ewi.tudelft.nl
More informationInteger Programming Formulations for the Minimum Weighted Maximal Matching Problem
Optimization Letters manuscript No. (will be inserted by the editor) Integer Programming Formulations for the Minimum Weighted Maximal Matching Problem Z. Caner Taşkın Tınaz Ekim Received: date / Accepted:
More informationMulti-objective Quadratic Assignment Problem instances generator with a known optimum solution
Multi-objective Quadratic Assignment Problem instances generator with a known optimum solution Mădălina M. Drugan Artificial Intelligence lab, Vrije Universiteit Brussel, Pleinlaan 2, B-1050 Brussels,
More informationLecture 15 (Oct 6): LP Duality
CMPUT 675: Approximation Algorithms Fall 2014 Lecturer: Zachary Friggstad Lecture 15 (Oct 6): LP Duality Scribe: Zachary Friggstad 15.1 Introduction by Example Given a linear program and a feasible solution
More information3.4 Relaxations and bounds
3.4 Relaxations and bounds Consider a generic Discrete Optimization problem z = min{c(x) : x X} with an optimal solution x X. In general, the algorithms generate not only a decreasing sequence of upper
More information16.1 Min-Cut as an LP
600.469 / 600.669 Approximation Algorithms Lecturer: Michael Dinitz Topic: LPs as Metrics: Min Cut and Multiway Cut Date: 4//5 Scribe: Gabriel Kaptchuk 6. Min-Cut as an LP We recall the basic definition
More informationA strongly polynomial algorithm for linear systems having a binary solution
A strongly polynomial algorithm for linear systems having a binary solution Sergei Chubanov Institute of Information Systems at the University of Siegen, Germany e-mail: sergei.chubanov@uni-siegen.de 7th
More informationILP-Based Reduced Variable Neighborhood Search for Large-Scale Minimum Common String Partition
Available online at www.sciencedirect.com Electronic Notes in Discrete Mathematics 66 (2018) 15 22 www.elsevier.com/locate/endm ILP-Based Reduced Variable Neighborhood Search for Large-Scale Minimum Common
More informationRevisiting the Hamiltonian p-median problem: a new formulation on directed graphs and a branch-and-cut algorithm
Revisiting the Hamiltonian p-median problem: a new formulation on directed graphs and a branch-and-cut algorithm Tolga Bektaş 1, Luís Gouveia 2, Daniel Santos 2 1 Centre for Operational Research, Management
More informationA Column Generation Based Heuristic for the Dial-A-Ride Problem
A Column Generation Based Heuristic for the Dial-A-Ride Problem Nastaran Rahmani 1, Boris Detienne 2,3, Ruslan Sadykov 3,2, François Vanderbeck 2,3 1 Kedge Business School, 680 Cours de la Libération,
More informationTime Aggregation for Network Design to Meet Time-Constrained Demand
20th International Congress on Modelling and Simulation, Adelaide, Australia, 1 6 December 2013 www.mssanz.org.au/modsim2013 Time Aggregation for Network Design to Meet Time-Constrained Demand N. Boland
More informationConvex and Semidefinite Programming for Approximation
Convex and Semidefinite Programming for Approximation We have seen linear programming based methods to solve NP-hard problems. One perspective on this is that linear programming is a meta-method since
More informationACO Comprehensive Exam October 14 and 15, 2013
1. Computability, Complexity and Algorithms (a) Let G be the complete graph on n vertices, and let c : V (G) V (G) [0, ) be a symmetric cost function. Consider the following closest point heuristic for
More informationEarliest Arrival Flows in Networks with Multiple Sinks
Electronic Notes in Discrete Mathematics 36 (2010) 607 614 www.elsevier.com/locate/endm Earliest Arrival Flows in Networks with Multiple Sinks Melanie Schmidt 2 Lehrstuhl II, Fakultät für Informatik Technische
More informationInteger Equal Flows. 1 Introduction. Carol A. Meyers Andreas S. Schulz
Integer Equal Flows Carol A Meyers Andreas S Schulz Abstract We examine an NP-hard generalization of the network flow problem known as the integer equal flow problem The setup is the same as a standard
More informationValid Inequalities and Restrictions for Stochastic Programming Problems with First Order Stochastic Dominance Constraints
Valid Inequalities and Restrictions for Stochastic Programming Problems with First Order Stochastic Dominance Constraints Nilay Noyan Andrzej Ruszczyński March 21, 2006 Abstract Stochastic dominance relations
More informationA cutting plane algorithm for the capacitated connected facility location problem
Noname manuscript No. (will be inserted by the editor) A cutting plane algorithm for the capacitated connected facility location problem Stefan Gollowitzer Bernard Gendron Ivana Ljubić Received: date /
More informationAggregation and Mixed Integer Rounding to Solve MILPs (Marchand and Wolsey)
Aggregation and Mixed Integer Rounding to Solve MILPs (Marchand and Wolsey) SAS Institute - Analytical Solutions Lehigh University - Department of Industrial and Systems Engineering July 7, 2005 Classical
More informationHow does data quality in a network affect heuristic
Working papers in transport, tourism, information technology and microdata analysis How does data quality in a network affect heuristic solutions? Ev. underrubrik ska stå i Arial, 14pt, Fet stil Mengjie
More informationPresolve Reductions in Mixed Integer Programming
Zuse Institute Berlin Takustr. 7 14195 Berlin Germany TOBIAS ACHTERBERG, ROBERT E. BIXBY, ZONGHAO GU, EDWARD ROTHBERG, AND DIETER WENINGER Presolve Reductions in Mixed Integer Programming This work has
More informationSequential pairing of mixed integer inequalities
Sequential pairing of mixed integer inequalities Yongpei Guan, Shabbir Ahmed, George L. Nemhauser School of Industrial & Systems Engineering, Georgia Institute of Technology, 765 Ferst Drive, Atlanta,
More informationChapter 7 Network Flow Problems, I
Chapter 7 Network Flow Problems, I Network flow problems are the most frequently solved linear programming problems. They include as special cases, the assignment, transportation, maximum flow, and shortest
More informationExploring the Numerics of Branch-and-Cut for Mixed Integer Linear Optimization
Zuse Institute Berlin Takustr. 7 14195 Berlin Germany MATTHIAS MILTENBERGER 1, TED RALPHS 2, DANIEL E. STEFFY 3 Exploring the Numerics of Branch-and-Cut for Mixed Integer Linear Optimization 1 Zuse Institute
More informationOn the Approximate Linear Programming Approach for Network Revenue Management Problems
On the Approximate Linear Programming Approach for Network Revenue Management Problems Chaoxu Tong School of Operations Research and Information Engineering, Cornell University, Ithaca, New York 14853,
More informationMVE165/MMG631 Linear and integer optimization with applications Lecture 8 Discrete optimization: theory and algorithms
MVE165/MMG631 Linear and integer optimization with applications Lecture 8 Discrete optimization: theory and algorithms Ann-Brith Strömberg 2017 04 07 Lecture 8 Linear and integer optimization with applications
More informationNetwork Flows. 6. Lagrangian Relaxation. Programming. Fall 2010 Instructor: Dr. Masoud Yaghini
In the name of God Network Flows 6. Lagrangian Relaxation 6.3 Lagrangian Relaxation and Integer Programming Fall 2010 Instructor: Dr. Masoud Yaghini Integer Programming Outline Branch-and-Bound Technique
More informationMATHEMATICAL ENGINEERING TECHNICAL REPORTS. The Symmetric Quadratic Semi-Assignment Polytope
MATHEMATICAL ENGINEERING TECHNICAL REPORTS The Symmetric Quadratic Semi-Assignment Polytope Hiroo SAITO (Communicated by Kazuo MUROTA) METR 2005 21 August 2005 DEPARTMENT OF MATHEMATICAL INFORMATICS GRADUATE
More informationMulti-Skill Resource-Constrained Project Scheduling: Formulation and Inequalities
Multi-Skill Resource-Constrained Project Scheduling: Formulation and Inequalities Isabel Correia, Lídia Lampreia Lourenço and Francisco Saldanha-da-Gama CIO Working Paper 17/2008 Multi-Skill Resource-Constrained
More informationInteger Linear Programs
Lecture 2: Review, Linear Programming Relaxations Today we will talk about expressing combinatorial problems as mathematical programs, specifically Integer Linear Programs (ILPs). We then see what happens
More informationMILP reformulation of the multi-echelon stochastic inventory system with uncertain demands
MILP reformulation of the multi-echelon stochastic inventory system with uncertain demands Axel Nyberg Åbo Aademi University Ignacio E. Grossmann Dept. of Chemical Engineering, Carnegie Mellon University,
More informationExtended Formulations, Lagrangian Relaxation, & Column Generation: tackling large scale applications
Extended Formulations, Lagrangian Relaxation, & Column Generation: tackling large scale applications François Vanderbeck University of Bordeaux INRIA Bordeaux-Sud-Ouest part : Defining Extended Formulations
More informationComputational Linear Algebra
Computational Linear Algebra PD Dr. rer. nat. habil. Ralf Peter Mundani Computation in Engineering / BGU Scientific Computing in Computer Science / INF Winter Term 2017/18 Part 2: Direct Methods PD Dr.
More informationMinimum Linear Arrangements
Minimum Linear Arrangements Rafael Andrade, Tibérius Bonates, Manoel Câmpelo, Mardson Ferreira ParGO - Research team in Parallel computing, Graph theory and Optimization Department of Statistics and Applied
More informationON SUM OF SQUARES DECOMPOSITION FOR A BIQUADRATIC MATRIX FUNCTION
Annales Univ. Sci. Budapest., Sect. Comp. 33 (2010) 273-284 ON SUM OF SQUARES DECOMPOSITION FOR A BIQUADRATIC MATRIX FUNCTION L. László (Budapest, Hungary) Dedicated to Professor Ferenc Schipp on his 70th
More informationThe Strength of Multi-Row Relaxations
The Strength of Multi-Row Relaxations Quentin Louveaux 1 Laurent Poirrier 1 Domenico Salvagnin 2 1 Université de Liège 2 Università degli studi di Padova August 2012 Motivations Cuts viewed as facets of
More information3.10 Lagrangian relaxation
3.10 Lagrangian relaxation Consider a generic ILP problem min {c t x : Ax b, Dx d, x Z n } with integer coefficients. Suppose Dx d are the complicating constraints. Often the linear relaxation and the
More informationTravelling Salesman Problem
Travelling Salesman Problem Fabio Furini November 10th, 2014 Travelling Salesman Problem 1 Outline 1 Traveling Salesman Problem Separation Travelling Salesman Problem 2 (Asymmetric) Traveling Salesman
More informationCAPACITATED LOT-SIZING PROBLEM WITH SETUP TIMES, STOCK AND DEMAND SHORTAGES
CAPACITATED LOT-SIZING PROBLEM WITH SETUP TIMES, STOCK AND DEMAND SHORTAGES Nabil Absi,1 Safia Kedad-Sidhoum Laboratoire d Informatique d Avignon, 339 chemin des Meinajariès, 84911 Avignon Cedex 09, France
More informationA comparison of sequencing formulations in a constraint generation procedure for avionics scheduling
A comparison of sequencing formulations in a constraint generation procedure for avionics scheduling Department of Mathematics, Linköping University Jessika Boberg LiTH-MAT-EX 2017/18 SE Credits: Level:
More informationSolving Railway Track Allocation Problems
Konrad-Zuse-Zentrum für Informationstechnik Berlin Takustraße 7 D-14195 Berlin-Dahlem Germany RALF BORNDÖRFER THOMAS SCHLECHTE Solving Railway Track Allocation Problems Supported by the Federal Ministry
More informationStrengthened Benders Cuts for Stochastic Integer Programs with Continuous Recourse
Strengthened Benders Cuts for Stochastic Integer Programs with Continuous Recourse Merve Bodur 1, Sanjeeb Dash 2, Otay Günlü 2, and James Luedte 3 1 Department of Mechanical and Industrial Engineering,
More informationModels and valid inequalities to asymmetric skill-based routing problems
EURO J Transp Logist (2013) 2:29 55 DOI 10.1007/s13676-012-0012-y RESEARCH PAPER Models and valid inequalities to asymmetric skill-based routing problems Paola Cappanera Luis Gouveia Maria Grazia Scutellà
More informationThinning out facilities: a Benders decomposition approach for the uncapacitated facility location problem with separable convex costs
Thinning out facilities: a Benders decomposition approach for the uncapacitated facility location problem with separable convex costs Matteo Fischetti 1, Ivana Ljubić 2, Markus Sinnl 2 1 Department of
More informationA Branch-and-Cut Algorithm for the Stochastic Uncapacitated Lot-Sizing Problem
Yongpei Guan 1 Shabbir Ahmed 1 George L. Nemhauser 1 Andrew J. Miller 2 A Branch-and-Cut Algorithm for the Stochastic Uncapacitated Lot-Sizing Problem December 12, 2004 Abstract. This paper addresses a
More informationFundamental Domains for Integer Programs with Symmetries
Fundamental Domains for Integer Programs with Symmetries Eric J. Friedman Cornell University, Ithaca, NY 14850, ejf27@cornell.edu, WWW home page: http://www.people.cornell.edu/pages/ejf27/ Abstract. We
More information1 Column Generation and the Cutting Stock Problem
1 Column Generation and the Cutting Stock Problem In the linear programming approach to the traveling salesman problem we used the cutting plane approach. The cutting plane approach is appropriate when
More informationRecognizing single-peaked preferences on aggregated choice data
Recognizing single-peaked preferences on aggregated choice data Smeulders B. KBI_1427 Recognizing Single-Peaked Preferences on Aggregated Choice Data Smeulders, B. Abstract Single-Peaked preferences play
More information(K + L)(c x) = K(c x) + L(c x) (def of K + L) = K( x) + K( y) + L( x) + L( y) (K, L are linear) = (K L)( x) + (K L)( y).
Exercise 71 We have L( x) = x 1 L( v 1 ) + x 2 L( v 2 ) + + x n L( v n ) n = x i (a 1i w 1 + a 2i w 2 + + a mi w m ) i=1 ( n ) ( n ) ( n ) = x i a 1i w 1 + x i a 2i w 2 + + x i a mi w m i=1 Therefore y
More informationChapter 3: Discrete Optimization Integer Programming
Chapter 3: Discrete Optimization Integer Programming Edoardo Amaldi DEIB Politecnico di Milano edoardo.amaldi@polimi.it Sito web: http://home.deib.polimi.it/amaldi/ott-13-14.shtml A.A. 2013-14 Edoardo
More informationVariable Objective Search
Variable Objective Search Sergiy Butenko, Oleksandra Yezerska, and Balabhaskar Balasundaram Abstract This paper introduces the variable objective search framework for combinatorial optimization. The method
More informationCSC Linear Programming and Combinatorial Optimization Lecture 12: The Lift and Project Method
CSC2411 - Linear Programming and Combinatorial Optimization Lecture 12: The Lift and Project Method Notes taken by Stefan Mathe April 28, 2007 Summary: Throughout the course, we have seen the importance
More information