Time Dependent Traveling Salesman Problem with Time Windows: Properties and an Exact Algorithm
|
|
- Miranda Stevens
- 5 years ago
- Views:
Transcription
1 Time Dependent Traveling Salesman Problem with Time Windows: Properties and an Exact Algorithm Anna Arigliano, Gianpaolo Ghiani, Antonio Grieco, Emanuela Guerriero Dipartimento di Ingegneria dell Innovazione, Università del Salento Via per Monteroni, Lecce, Italy Abstract In this paper, we deal with the Time-Dependent Traveling Salesman Problem with Time Windows (T DT SP T W ). Firstly, we prove that under special conditions the T DT SP T W can be solved as an Asymmetric Traveling Salesman Problem with Time Windows (AT SP T W ), with suitable-defined time windows and (constant) travel times. Secondly, we show that, if the special conditions do not hold, the AT SP T W optimal solution provides both a lower bound and (eventually) an upper bound with a worst-case guarantee for the original T DT SP T W. Finally an integer linear programming model is presented and valid inequalities are embedded into a branch-and-cut algorithm. Computational results show that the proposed algorithm is able to solve instances with up to 40 vertices. Keywords: traveling salesman problem, time dependence, time windows, lower and upper bounds, branch-and-cut 1. Introduction This paper presents an exact algorithm for solving the time-windowed version of the Time-Dependent Traveling Salesman Problem defined as follows. Let G = (V {0}, A) be a complete directed graph, where V = {1,..., n} is the set of n customers to be visited, vertex 0 is the depot, and A is the set of arcs. To each arc (i, j) A is associated a non-negative travel distance L ij. The depot is the starting point of the vehicle. Let τ ij (t) be the travel duration of arc (i, j) A when the vehicle leaves vertex i at time t. Moreover, emanuela.guerriero@unisalento.it Preprint submitted to European Journal of Operational Research February 10, 2014
2 at each i V {0}, are associated an earliest and a latest service time r i and d i, respectively. The interval [r i, d i ] is called the time window of node i. The TDTSPTW consists in finding a least-duration Hamiltonian circuit on G, starting from the depot at time t = 0, such that for every node i V the arrival time at node i V is lower than or equal to d i. One may arrive at node i V earlier than r i and it is allowed to wait until the node becomes available for service at time r i. Finally, it is supposed that the time window of node depot is relaxed. The earliest papers on the T SP T W focused on exact approaches are [1] and [2]. Later, in [3] Dumas et al. extended earlier dynamic programming approaches by using state space reduction techniques that allowed them to solve larger problem instances. Ascheuer et al. [4] developed a branch-andcut algorithm in which they applied techniques tailored for the asymmetric T SP T W. Balas and Simonetti [5] presented a dynamic programming algorithm for several TSP variants with precedence constraints, including the TSPTW. Constraint programming is another exact approach that has been recently applied to the TSPTW in [6] and [7]. The main limit of these contributions is due to hypothesis of constant travel times. Indeed, due to the road congestion, which varies over the day, solutions derived from time-independent models to the T SP T W could be infeasible when applied in real life. Despite numerous publications deal with exact algorithms for the time independent vehicle routing problems, few address the inherent time-dependent nature of this problems. Recently, Ropke et al. in [8] have proposed a Branch-and-Price algorithm for the Time-Dependent vehicle routing problem with time windows. With respect to vehicle routing problems with a single vehicle including time-dependent travel times, recently Cordeau et. al. [9] exploited a number of properties and devised an exact algorithm. As far as we know, the only exact algorithm for the TDTSPTW is given in [10]. Here the problem is firstly transformed into an Asymmetric Generalized TSP and then into an Asymmetric Graphical TSP. In this way, the authors can apply a known exact algorithm for the Mixed General Routing Problem. This approach can solve instances with up to 60 vertices, with very tight time windows, resulting in relatively small Asymmetric Graphical TSP instances. The contribution of this paper is threefold. Firstly, we prove that under special conditions the T DT SP T W can be solved as an Asymmetric Traveling Salesman Problem with Time Windows (AT SP T W ), with suitable-defined time windows and (constant) travel times. Secondly, we show that, if the 2
3 special conditions do not hold, the AT SP T W optimal solution provides both a lower bound and (eventually) an upper bound with a worst-case guarantee for the original T DT SP T W. Finally an integer linear programming model is formulated for TDTSPTW. Valid inequalities are embedded into a branchand-cut algorithm, that was able to solve instances with up to 40 vertices. The remainder of the paper is organized as follows. The next section introduces the notation used throughout the paper. Section 3 illustrates some properties of the T DT SP T W. In Section 4 a procedure to compute lower and upper bounds on the optimal solution value is described. It is shown that the lower bound corresponds to the optimal solution value when all arcs share a common congestion pattern. Section 5 introduces an integer linear programming formulation of the problem. Section 6 is devoted to valid inequalities that strengthen the linear programming relaxation of the formulation given in Section 5. These inequalities are then used within a branch-and-cut algorithm which is described in Section 7. This is followed by computational experiments in Section 8, and by the conclusion. 2. Notation Let [0, T ] be the time horizon (typically a single day) within which a tour may be completed. This horizon is partitioned into H subintervals [T h, T h+1 ] (h = 0,..., H 1). We assume that, for each arc (i, j) A, traffic information systems provide the average value of travel speed v ijh during time interval [T h, T h+1 ] with h = 0,..., H 1. The time-dependent travel value τ ij (t) on arc (i, j) A at time t can be computed as proposed by Ichoua et al. in [11]. As stated in [12], following a generic piecewise linear travel time model allows to express the travel speeds as where: v ijh = δ ijh b h u ij (1) u ij is the maximum travel speed across arc (i, j) A during [0, T ],i.e. u ij = max v ijh; h=0,...,h 1 b h belongs to [0, 1] and is the best (i.e. lightest) congestion factor during interval [T h, T h+1 ] on the complete graph, i.e. b h = max v ijh /u ij ; (i,j) A 3
4 δ ij belongs to [0, 1] and represents the degradation of the congestion factor of arc (i, j) in interval [T h, T h+1 ] with respect to the less congested arc in [T h, T h+1 ]. Let = min ijh (δ ijh) be the heaviest degradation of the traffic jam factor of any arc (i, j) A over [0, T ]. If = 1, then all arcs (i, j) A share a common traffic jam factor b h during interval [T h, T h+1 ] (h = 1,..., H 1). Hence the relationship (1) can be expressed as: v ijh = b h u ij. (2) Given an arc (i, j) A, under speed variation law (2), it results that: L ij u ij = t+τij (t) t b(µ)dµ, (3) where b(t) = b h t [T h, T h + 1[. (4) 3. Properties Let p i be a simple path {i 0, i 1,..., i m } on G with i 0 = 0 and i m = i. Given a circuit c, let σ(0), σ(1),..., σ(n) denote the sequence of customers visited in positions 0, 1,..., n respectively. When the vehicle traverses the path p σ(i) under the speed variation law (1) starting from the depot at time t = 0, y(p σ(i) ) indicates the arrival time in σ(i), with i = 1,..., n. Since the vehicle might arrive before r σ(i), we also define the visiting time as z(p σ(i) ) = max{r σ(i), y(p σ(i) )}. A Hamiltonian circuit c is said to be feasible under speed variation law (1) if each node i V is visited within the time window [r i, d i ], with i = 1,..., n. We denote with z(c) the time needed to traverse a Hamiltonian circuit c on G under speed variation law (1). Moreover, let C denote the set of all Hamiltonian circuits on G feasible under speed variation law (1). Similarly, z(c) and C will indicate the duration of c and the set of feasible Hamiltonian circuits when the speed variation law (2) holds. Of course it results that:. C C 4
5 If all v ijh are further increased to u ij, travel speeds become constant and the problem reduces to the ATSP with Time Windows. In particular we denote with Ω the set of Hamiltonian circuits on G such that if c is traversed with constant speeds u ij, then each node i V is visited within the time window [R i, D i ], where R i = ri 0 b(µ)dµ, D i = di 0 b(µ)dµ. (5) Under the hypothesis of constant speed values, we denote with z(c) the duration of c Ω. Finally, as we did for the duration of a circuit, we also define z(p σ(i) ), y(p σ(i) ), z(p σ(i) ) and y(p σ(i) ). If all time windows [r i, d i ] are relaxed then any Hamiltonian circuit c becomes feasible under all three considered speed variation laws. This is the hypothesis underlying the TDTSP, for which Cordeau et al. in [9] demonstrated that, given a circuit c, the relationship (6) holds. z(p σ(i) ) = z(pσ(i) ) 0 b(µ)dµ. (6) The following theorem extends this result to the T DT SP T W. Theorem 1. For each c C, relationship (6) holds. Proof. Let us denote with m the number of customers, that contribute with a not-null waiting time to z(p σ(i) ). We demonstrate the thesis by induction on m. Case (m = 0). In this case, only travel times contribute to the computation of z(p σ(i) ). Since this is the hypothesis underlying the TDTSP, we can assert that Cordeau et al. in [9] have already proved that relationship (6) holds with m = 0. Case (m > 0). Let us suppose that (6) holds for (m 1). Let σ(q) denote the node at which the m th not-null waiting time occurs under speed variation law (2), with q = 1,..., n. Then it results that and y(p σ(q) ) < r σ(q) (7) 5
6 z(p σ(q) ) = r σ(q). We observe that only travel times contribute to the computation of z(p σ(i) ), with i = q + 1,..., n. Then from (5) it follows that the thesis is proved if we demonstrate that y(p σ(q) ) < R σ(q), that is z(p σ(q) ) = R σ(q). We start by observing that and y(p σ(q) ) = z(p σ(q 1) ) + L σ(q 1)σ(q) u σ(q 1)σ(q), y(p σ(q) ) = z(p σ(q 1) ) + τ σ(q 1)σ(q) (z(p σ(q 1) )). Since (6) holds for m 1, from (3) and (7) it follows that : y(p σ(q) ) = which completes the proof. y(pσ(q) ) 0 b(µ)dµ < rσ(q) 0 b(µ)dµ = R σ(q), The following Corollary extends Theorem 1 to durations of circuits. Corollary 1. For each c C, relationship (8) holds. z(c) = z(c) 0 b(µ)dµ. (8) Proof. We start by observing that z(c) and z(c) can be expressed as : and z(c) = z(σ(n)) + τ σ(n),0 (z(σ(n))) (9) z(c) = z(p(σ(n))) + L σ(n) u σ(n),0. (10) The thesis is proved by applying Theorem 1 and (3) to relationships (9) and (10). 6
7 4. Computing Lower and Upper Bounds We start by observing that min(z(c) c C) and min(z(c) c C) are both TDTSPs with time windows, meanwhile min(z(c) c Ω) corresponds to the classical ATSPTW. Let c, c be optimal solutions of the TDTSPTW s under the speed variation law (1) and (2), i.e. c = arg min{z(c) c C} and c = arg min{z(c) c C}. Similarly, let c denote an optimal solution of the ATSPTW defined on Ω, that is c = arg min{z(c) c Ω}. In the following we demonstrate that efficient algorithms developed for the AT SP T W can be used to solve the T DT SP T W when speed variation law (2) holds A Lower Bound In this section we demonstrate both that z(c ) is a lower bound on z(c ) and that if speed variation law (2) holds, then c is optimal for TDTSPTW, i.e. z(c ) = z(c ). Theorem 2. Circuit c is optimal under speed variation law (2). Proof. From Theorem 1 and definition (5), it follows that for each c C: z(p σ(i) ) d pσ(i) z(p σ(i) ) D pσ(i). (11) This implies that circuits belonging to Ω are feasible under speed variation law (2),that is C Ω. Moreover from Corollary 1 we have that, given two generic circuits c and c of C: z(c) z(c ) z(c) z(c ). (12) An implication of (11) and (12) is that z(c ) = min{z(c) c C}. Corollary 2. Given two Hamiltonian circuits c Ω and c C (optimal with respect to z(c)) and z(c), respectively), the following relationship is satisfied: z(c ) z(c ) (13) Proof. Given a generic solution c C, from Theorem 2 and Corollary 1 one has z(c ) z(c) z(c ) z(c). 7
8 Since c remains feasible under speed variation law (2), i.e. C C, it follows that z(c ) z(c ) z(c ) Worst case performance analysis In this subsection we provide a worst case analysis on z(c ). Theorem 3. If c is feasible under speed variation law (1), then z(c ) is an upper bound not greater than 1/ times the optimal value z(c ), that is z(c ) 1 z(c ). Proof. Given the speed variation law (1), during each time interval h the congestion factor may take a value in the interval [ b h, b h ] with h = 0,..., H 1. Therefore, for a generic solution c C, one has Combining (14) and (13) we may assert that z(c) z(c) 1 z(c). (14) z(c ) 1 z(c ) 1 z(c ). 5. Integer Linear Programming Formulation In this section, the T DT SP T W is formulated as an integer linear program for which new valid inequalities are devised. Let x ij be a binary decision variable equal to 1 iff arc (i, j) A is part of the solution. Given a customer i, we denote with P i the set of all simple paths {i 0, i 1,..., i m } on G with i 0 = 0 and i m = i. Let denote by LB p a lower bound on the tour duration when the path p P i, i V, is part of the solution. As in [13], the time window restrictions are modeled by a class of inequalities called unfeasible path elimination constraints (IP ECs). For this purpose, we denote with x(p) the sum of variables corresponding to path p, i.e. x(p) = m 1 l=0 x i l,i l+1. Finally let C p denote the set of circuits, which the path p is part of. The T DT SP T W is formulated as follows. 8
9 M inimize z (15) subject to z LB p 1 + (x ab 1), i V {0}, p P i (16) (a,b) p i V \{j} j V \{i} x ij = 1, j V {0} (17) x ij = 1, i V {0} (18) x ij 1, S V {0}, S 2 (19) i S j / S x(p) m 1 C p C = (20) x ij {0, 1}, (i, j) A. (21) Equation (16) states the relationships between the tour duration variable z and the flow variables x ij : the value of z must be larger than or equal to the longest path belonging to the circuit defined by the x ij variables taking value 1. Equations (17) and (18) are referred to as degree constraints. Constraints (17) impose that a unique arc enters each vertex j V. Similarly, constraints (18) state that a single arc exits each vertex i V. Inequalities (19) are the subtour elimination constraints (SECs). Finally, (20) forbid infeasible paths, i.e. paths violating time windows (IP ECs) Computing LB p Let z(p, t) denote the visiting time of node i, when the vehicle traverses the path p = {i 0 = 0,..., i m = i} under variation law (1) starting from the depot at time t > 0. As it has been made for the starting time t = 0, z(p, t) and z(p, t) are defined. Theorem 4 shows how one can use z(c ) to compute a lower bound on z(c) with c C p C. Given a start travel time t, to ease the explanation, we will denote by τ(l, t) the traversal time of a dummy arc which has a length equal to L and travel speed values equal to b h for h = 0,..., H 1. 9
10 Theorem 4. Given a path p with C p C, if z(c ) z(p) it follows that where c C p. z(c) z(p) + τ(z(c ) z(p), z(p)), (22) Proof. Each Hamiltonian circuit c C p can also be modeled as the concatenation of p and p c, where p c represents the simple path starting from i. Therefore the circuit duration z(c) can be rewritten as: z(c) = z(p) + z(p c, z(p)). Since z(p) is the same for any c C p, a lower bound on z(p c, z(p)) can be determined. From Corollary 1 it can be said that z(p c, z(p c )) = τ(z(p c ), z(p)) z(p c, z(p)) c C p. (23) Since c = arg min{z(c)}, it follows that c Ω z(c ) z(p) z(p c ) c C p. (24) Combining Theorem 1 with (23) and (24) yields 0 z(c ) z(p) z(p c ) τ(z(c ) z(p), z(p)) τ(z(p c ), z(p)) z(p c, z(p)) c Ω p, which completes the proof. We observe that another lower bound on z(p c, z(p)) is r kc + τ kc,0 (r kc ), where k c = arg max(r i i p c ). Therefore we adopt the procedure described in Algorithm 1 to calculate the LB p value. 6. Valid Inequalities The linear programming relaxation of formulation (15) (21) can be strengthened with two families of valid inequalities. 10
11 Algorithm 1 Computing the lower bound LB p if z(c ) < z(p) then LB p = z(p) else LB p = z(p) + max{τ(z(c ) z(p), z(p)), r kc + τ kc,0 (r kc )} end if return LB p 6.1. Lower bounding inequalities We introduce A additional non-negative variables t ij representing the visiting time of node i under constant travel speeds (i.e. v ij (t) = u ij ), whenever x ij = 1. Moreover we introduce a non-negative variable Γ, representing the circuit duration z(c). The following inequalities state that c Ω. t ij + i V i V Γ j V {0} L ij x ij t jk j V (25) u ij k V t j0 + j V {0} L j0 u j0 x j0 (26) R i x ij t ij D i x ij (i, j) A (27) Γ 0 (28) t ij 0 (i, j) A (29) By observing that Theorem 2 also demonstrated that C Ω( C), then we can assert that inequalities (25)-(29) are valid. If 1, it results that z(c) z(c) for any c C. In this case, the following valid inequalities can greatly speed up the solution process: z z(c) c C. (30) We assume that T k z(c) T k+1 and k {0,..., H 1}. Since the function b(µ) is constant stepwise, then the relationship (8) can be rewritten as follows: z(c) = k (T h+1 T h )b h + (z(c) T k )b k. (31) h=0 11
12 The relationship (31) is equivalent to z(c) = H 1 h=0 z h (c)/b h, where, as reported in [9], z 0 (c),..., z H 1 (c) have been recursively defined in the following { } z H 1 (c) = max 0, z(c) H 2 { h 1 z h (c) = max 0, z(c) b l [T l+1 T l ] l=0 z 0 (c) = max { l=0 0, z(c) b l [T l+1 T l ] H 1 l=h+1 H 1 l=1 z l (c) z l (c) } }.. h = 1,..., H 2 Therefore the inequality (30) can be rendered through the following linear inequalities: z H 1 h=0 z 0 Γ h 1 z h Γ b l (T l T l 1 ) l=0 z H 1 Γ H 2 l=0 z h b h (32) H 1 l=1 H 1 l=h+1 z l (33) z l h = 1,.., H 2 (34) b l (T l T l 1 ) (35) H 1 z h b h (T h T h 1 ) y l h = 0,.., H 1 (36) z h b h (T h T h 1 ) l=h H 1 l=h+1 12 y l h = 0,.., H 2 (37)
13 y h = 1 (38) H 1 h=0 z h 0, y h {0, 1} h = 0,..., H 1, (39) Z 0, (40) where the binary variable y h is equal to 1 iff z(c) [T h 1, T h ], with h = 0,..., H 1. Since the number H of periods is typically rather small, no separation procedure has been devised for the (3H+1) valid inequalities (32) (39), which have been directly added as standard constraints to the integer linear programming formulation. Similarly we dealt with the (n(1 + n) + 1) valid inequalities (25) (27) Infeasible Path Elimination Constraints Given a simple path P = (v 1,..., v k ), it results infeasible if time windows are violated. The infeasible paths have to be forbidden. The tournament constraints (41), proposed in [13] for the AT SP T W, are also valid for the T DT SP T W. x([p ]) := k 1 k i=1 j=i+1 7. Branch-and-Cut Algorithm x vi v j k 2 ( P 1) (41) In this section the main elements of a branch-and-cut algorithm based on formulation (15) (21) are described Data preprocessing As for many other combinatorial optimization problems, data preprocessing plays a central role in the efficient implementation of the proposed branch-and-cut algorithm. During preprocessing, the aim is to determine a tighter formulation of the problem without loosing any optimal solution of the original one. Since the set of feasible circuits C is a subset of Ω, we apply some of the preprocessing steps proposed in [4] for the AT SP T W, i.e. min(z(c) c Ω). In particular we considered three main steps: tightening the time windows, constructing precedences among the nodes, and fixing variables permanently. 13
14 It is worth noting that the tightening of the time windows has been carried out into two steps. Firstly, the release date and due-date adjustment procedures proposed in [4] are applied to R i and D i, with i = 1,..., n. By extending the notation introduced in the previous section to the relationships (5), we have that : r(i) = τ(r(i), 0) d(i) = τ(d(i), 0), with i = 1,..., n. Therefore, during the second step, the corresponding adjustments on r i and d i are obtained as output of the travel time computation algorithm proposed in [11] Initialization Because z(c ) represents a lower bound on the optimal value z(c ), the following valid inequality has been added, as initial lower bound, to the model: z z(c ) Furthermore we tried to determine an initial feasible solution as follows. If c C, then a start upper bound is provided by the value z(c ). Otherwise we try to determine an initial upper bound value by running the time-dependent version of the construction and improvement heuristics proposed in [4] for AT SP T W Separation Routines The proposed formulation has an exponential number of constraints (16), (19) and (20). For this reason, procedures to identify violated inequalities, are used. For inequalities (16), we adopted the separation procedure proposed in [9]. For a given solution x, with each arc (i, j) A is associated a cost 1 ˆx ij, where ˆx ij denotes the value of variable x ij in the current fractional solution. We determine the least cost path from vertex 0 to any other vertex with respect to these costs. If these shortest path costs are greater than or equal to 1 then no violated constraint (16) exists. Otherwise, the shortest path corresponds to a violated inequality of the form (16). Conversely, if a violated constraint (16) exists, then there is at least one path with a cost smaller than 1 and it will be detected by our procedure. For the subtour elimination constraints (19), an appropriate separation routine, based on the minimum cut algorithm of [14], is used. Tournments constraints (20) are separated with the simple procedure defined in [4]. Let ˆx be a fractional 14
15 solution. The paths P k for which k 1 k ˆx vi v j k+2 is greater than 0, under i=1 j=i+1 the assumption that ˆx ij = ˆx ij = 1, are detected by enumeration with i V j V backtrack as soon as k 1 k i=1 j=i+1 ˆx vi v j k Computational Results The branch-and-cut algorithm was coded in C++ and run on MacBook Pro with a 2.33-GHz Intel Core 2 Duo processor and 4 GB of memory. The Concert Technology library of Cplex 11 was used as a framework. For our numerical study we use data sets generated in [9] for the T DT SP, where traffic density was characterized with respect to three concentric circular zones C 1, C 2 and C 3 (see Figure 1). For each possible value of n in the set {15, 20, 30, 40} Cordeau et al. [9] generated 30 instances, that is 30 distance matrices. Two distinct traffic density patterns, A and B, have been considered as indicated in Tables 1 and 2. As stated in [9], traffic pattern A is typical of cities where the old town centre C 1 is a limited traffic zone and has a low traffic congestion comparable to the residential zone C 3. In this scenario the value refers to the business centre located in the middle of zone C 2. The traffic pattern B is more representative of cities where the highest traffic congestion, i.e. the value, characterizes the city centre. Under this scenario, the middle zone C 2 experiences medium congestion as the traffic thins out gradually from the city centre C 1 to the residential zone C 3. For each traffic pattern we have considered five different values ranging from 0.70 to For each zone c and each time interval h, the δ ijh values of arcs (i, j) with vertex i in zone c have been set within the intervals reported in Tables 1 and 2. Given a T DT SP instance we generated a feasible T DT SP T W as follows. Firstly n release dates were randomly generated. Then we determined a circuit c, where customers are visited in no-increasing order of release date. The main issue was to assign a time window width to each customer so that c is feasible. Therefore, the duedate value is set to be 30% greater than the visiting time in c, that is d σ(i) = 1.30 z(p σ(i) ) with i = 1,..., n. 15
16 Figure 1: Customer locations Combining the four possible values for n, five values and two traffic patterns yields a total of 1200 instances. Table 1: Intervals containing δ ijh for traffic pattern A h = 1 h = 2 h = 3 C C 2 [0.70; 0.98] 1 [0.70; 0.98] C 3 [0.95; 0.99] 1 [0.95; 0.99] A maximum of 3600 seconds was imposed on the solution of each instance. To determine c we formulated the corresponding AT SP T W instance and we solved it using the separation procedures proposed for (19) and (20) in [14] and [13], respectively. Results are reported in Table 3 and Table 4 under the following headings: OP T : number of instances solved to optimality; T IME: time spent to determine an initial feasible solution; z(c ) z(c ): ratio of the initial upper bound on the best lower bound available at the end of the search; 16
17 Table 2: Intervals containing δ ijh for traffic pattern B h = 1 h = 2 h = 3 C 1 [0.70; 0.98] 1 [0.70; 0.98] C 2 [0.80; 0.99] 1 [0.80; 0.99] C GAP I : the initial optimality gap in percentage, tha is (z(c ) z(c )) z(c ) GAP f : final optimality gap in percentage; NODES: number of nodes in the search tree; (19): total number of subtour elimination constraints (19) added by the algorithm; (20): total number of infeasible path elimination constraints (20) added by the algorithm; (16) : total number of constraints (16) added by the algorithm; T IME: time spent to determine the optimal solution. For column T IME we report averages across all instances. For columns from z(c )/z(c ) to GAP f, we report results on two distinct rows: the first row is the average across instances solved to optimality, while the second row is the average for the remaining instances. For the sake of conciseness, the first or the second row has been omitted whenever none or all instances were solved to optimality. For columns from NODES to T IME we report only averages for instances that were solved to optimality. When the time limit was reached before proving optimality, the ratio z(c )/z(c ) was evaluated with respect to the best lower bound found. Optimality gaps were computed as upper bound lower bound. (42) lower bound The final optimality GAP f is zero if the instance was solved to optimality. Otherwise, GAP f was evaluated with respect to the best lower and upper bounds found within the time limit of 3600 seconds. It is worth noting that, 17
18 at the root node the LP relaxation of (15)-(21), inequalities (25)-(29) and (32)-(39) were never able to improve the lower bound z(c ). Therefore, the value GAP I can also be interpreted as the root node optimality gap. Computational results confirms that the proposed algorithm exhibits a computational behavior similar to that one of the branch-and-cut algorithm proposed in [9] for the T DT SP. The time needed to compute the initial solution is usually below 1 second and the initial gap was below 6%. When traffic pattern A holds, the proposed formulation is tight and the branchand-cut algorithm is capable of solving all instances with 15 customers and most instances containing up to 40 customers. Under this traffic congestion pattern, the percentage of arcs (i, j) A characterized by δ ijh 1 is greater than 60% and z(c ) z(c ).When traffic pattern B holds, the success rate decreases, because a large fraction of the arcs appearing in any solution belong to the highly congested clusters C 1 and C 2. The algorithm was able to solve to optimality only instances with an initial gap value GAP I usually lower than 1% and 2.5% under traffic pattern A and B respectively. Moreover we also observe that for all instances solved to optimality the algorithm closed the optimality gap within the first 60 seconds of computation on average. In all other cases, the deviation between the best lower bound and the heuristic solution is usually below 5%. 9. Conclusions In this research work, an exact algorithm for the T DT SP T W has been proposed. In particular, a new integer programming formulation of the problem has been defined and valid inequalities have been introduced that have been embedded into a branch-and-cut algorithm. Computational experiments show that the algorithm can successfully solve instances derived from the literature on ATSPTW with up to 40 vertices. Future research will be focused on more generic test cases in order to assess the relationships existing between the value and the performance of the proposed algorithm. References [1] N. Christofides, A. Mingozzi, P. Toth, State-space relaxation procedures for the computation of bounds to routing problems, Networks 11 (2) (1981)
19 [2] E. K. Baker, An exact algorithm for the time-constrained traveling salesman problem, Networks 31 (5) (1983) [3] Y. Dumas, J. Desrosiers, E. Gelinas, An optimal algorithm for the traveling salesman problem with time windows, Operations Research 43 (2) (1995) [4] N. Ascheuer, M. Fischetti, G. M, A polyhedral study of the asymmetric traveling salesman problem with time windows, Networks 90 (2000) [5] E. Balas, N. Simonetti, Linear time dynamic-programming algorithms for new classes of restricted tsps: A computational study, INFORMS Journal on Computing 13 (1) (2001) [6] F. Focacci, A. Lodi, M. Milano, A hybrid exact algorithm for the tsptw, INFORMS Journal on Computing (2002) [7] G. Pesant, M. Gendreau, J.-Y. Potvin, J.-M. Rousseau, Logic programming algorithm for the traveling salesman problem with time windows, Transportation Science 32 (1998) [8] S. Dabia, S. Ropke, T. W. andt. G. de Kok, Branch and price for the time-dependent vehicle routing problem with time windows, Transportation Science 47 (3) (2013) [9] J. F. Cordeau, G. Ghiani, E. Guerriero, Analysis and branch-and-cut algorithm for the time-dependent travelling salesman problem, Transportation Science (2011) Published on line before print, December 13, 2012, doi: /trsc [10] J. Albiach, J. Sanchis, D. Soler, An asymmetric TSP with time windows and with time-dependent travel times and costs: An exact solution through a graph transformation, European Journal of Operational Research 189 (2008) [11] I. Soumia, M. Gendreau, J.-Y. Potvin, Vehicle dispatching with timedependent travel times, European Journal of Operational Research 144 (2003)
20 [12] G. Ghiani, E. Guerriero, A note on the ichoua et al. (2003) travel time model, Transportation Science Accepted for publication (2013). [13] N. Ascheuer, M. Fischetti, G. M, Solving the asymmetric travelling salesman problem with time windows by branch-and-cut, Mathematical Programming 36 (2001) [14] H. Nagamochi, T. Ono, T. Ibaraki, Implementing an efficient minimum capacity cut algorithm, Math. Program. 67 (1994)
21 Table 3: Computational results for traffic pattern A V OP T TIME z(c )/z(c ) GAP I GAP f NODES SEC IP EC LB T IME
22 Table 4: Computational results for traffic pattern B V OP T TIME z(c )/z(c ) GAP I GAP f NODES SEC IP EC LB T IME
A Time Bucket Formulation for the TSP with Time Windows
A Time Bucket Formulation for the TSP with Time Windows Sanjeeb Dash, Oktay Günlük IBM Research Andrea Lodi, Andrea Tramontani University of Bologna November 10, 2009 Abstract The Traveling Salesman Problem
More informationSolving Time Dependent Traveling Salesman Problems with Time Windows
Solving Time Dependent Traveling Salesman Problems with Time Windows Duc Minh Vu, Mike Hewitt 1 and Natashia Boland, Martin Savelsbergh 2 1 Department of Information Systems and Supply Chain Management,
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 informationSolving a Production Scheduling Problem as a Time-Dependent Traveling Salesman Problem
Solving a Production Scheduling Problem as a Time-Dependent Traveling Salesman Problem GABRIELLA STECCO Department of Applied Mathematics, University Ca Foscari of Venice, Dorsoduro n. 3825/E, 30123 Venice,
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 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 informationChapter 3: Discrete Optimization Integer Programming
Chapter 3: Discrete Optimization Integer Programming Edoardo Amaldi DEIB Politecnico di Milano edoardo.amaldi@polimi.it Website: http://home.deib.polimi.it/amaldi/opt-16-17.shtml Academic year 2016-17
More informationIntroduction to Mathematical Programming IE406. Lecture 21. Dr. Ted Ralphs
Introduction to Mathematical Programming IE406 Lecture 21 Dr. Ted Ralphs IE406 Lecture 21 1 Reading for This Lecture Bertsimas Sections 10.2, 10.3, 11.1, 11.2 IE406 Lecture 21 2 Branch and Bound Branch
More 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 informationPartial Path Column Generation for the Vehicle Routing Problem with Time Windows
Partial Path Column Generation for the Vehicle Routing Problem with Time Windows Bjørn Petersen & Mads Kehlet Jepsen } DIKU Department of Computer Science, University of Copenhagen Universitetsparken 1,
More 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 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 informationThe Traveling Salesman Problem with Pickup and Delivery. A polyhedral approach. IBM Research - Australia. Irina Dumitrescu
Australia The Traveling Salesman Problem with Pickup and Delivery A polyhedral approach Irina Dumitrescu Jean-Francois Cordeau, Gilbert Laporte, Stefan Ropke The TSP with Pickup and Delivery (TSPPD) Given:
More informationA Branch-and-Cut Algorithm for the Dial-a-Ride Problem
A Branch-and-Cut Algorithm for the Dial-a-Ride Problem JEAN-FRANÇOIS CORDEAU Canada Research Chair in Distribution Management, HEC Montréal 3000, chemin de la Côte-Sainte-Catherine Montréal, Canada H3T
More informationPart III: Traveling salesman problems
Transportation Logistics Part III: Traveling salesman problems c R.F. Hartl, S.N. Parragh 1/74 Motivation Motivation Why do we study the TSP? it easy to formulate it is a difficult problem many significant
More informationBounds on the Traveling Salesman Problem
Bounds on the Traveling Salesman Problem Sean Zachary Roberson Texas A&M University MATH 613, Graph Theory A common routing problem is as follows: given a collection of stops (for example, towns, stations,
More informationNew Integer Programming Formulations of the Generalized Travelling Salesman Problem
American Journal of Applied Sciences 4 (11): 932-937, 2007 ISSN 1546-9239 2007 Science Publications New Integer Programming Formulations of the Generalized Travelling Salesman Problem Petrica C. Pop Department
More informationIntroduction to Integer Programming
Lecture 3/3/2006 p. /27 Introduction to Integer Programming Leo Liberti LIX, École Polytechnique liberti@lix.polytechnique.fr Lecture 3/3/2006 p. 2/27 Contents IP formulations and examples Total unimodularity
More 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 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 informationEXACT ALGORITHMS FOR THE ATSP
EXACT ALGORITHMS FOR THE ATSP Branch-and-Bound Algorithms: Little-Murty-Sweeney-Karel (Operations Research, ); Bellmore-Malone (Operations Research, ); Garfinkel (Operations Research, ); Smith-Srinivasan-Thompson
More informationPart III: Traveling salesman problems
Transportation Logistics Part III: Traveling salesman problems c R.F. Hartl, S.N. Parragh 1/282 Motivation Motivation Why do we study the TSP? c R.F. Hartl, S.N. Parragh 2/282 Motivation Motivation Why
More informationDecision Diagrams for Sequencing and Scheduling
Decision Diagrams for Sequencing and Scheduling Willem-Jan van Hoeve Tepper School of Business Carnegie Mellon University www.andrew.cmu.edu/user/vanhoeve/mdd/ Plan What can MDDs do for Combinatorial Optimization?
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 informationValid Inequalities and Separation for the Symmetric Sequential Ordering Problem
Valid Inequalities and Separation for the Symmetric Sequential Ordering Problem Adam N. Letchford Yanjun Li Draft, April 2014 Abstract The sequential ordering problem (SOP) is the generalisation of the
More informationPricing Routines for Vehicle Routing with Time Windows on Road Networks
Pricing Routines for Vehicle Routing with Time Windows on Road Networks Adam N. Letchford Saeideh D. Nasiri Amar Oukil Published in Computers & Operations Research, July 2014 Abstract Several very effective
More informationAsymmetric Traveling Salesman Problem (ATSP): Models
Asymmetric Traveling Salesman Problem (ATSP): Models Given a DIRECTED GRAPH G = (V,A) with V = {,, n} verte set A = {(i, j) : i V, j V} arc set (complete digraph) c ij = cost associated with arc (i, j)
More informationVehicle Routing with Traffic Congestion and Drivers Driving and Working Rules
Vehicle Routing with Traffic Congestion and Drivers Driving and Working Rules A.L. Kok, E.W. Hans, J.M.J. Schutten, W.H.M. Zijm Operational Methods for Production and Logistics, University of Twente, P.O.
More informationA Branch-and-Cut Algorithm for the Periodic Rural Postman Problem with Irregular Services
A Branch-and-Cut Algorithm for the Periodic Rural Postman Problem with Irregular Services Enrique Benavent Ángel Corberán Demetrio Laganà Francesca Vocaturo DEIO, Universitat de València, Spain DIMEG,
More informationAn Exact Algorithm for the Traveling Salesman Problem with Deliveries and Collections
An Exact Algorithm for the Traveling Salesman Problem with Deliveries and Collections R. Baldacci DISMI, University of Modena and Reggio E., V.le Allegri 15, 42100 Reggio E., Italy E. Hadjiconstantinou
More informationA priori performance measures for arc-based formulations of the Vehicle Routing Problem
A priori performance measures for arc-based formulations of the Vehicle Routing Problem Fernando Ordóñez, Ilgaz Sungur, and Maged Dessouky Industrial and Systems Engineering, University of Southern California
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 informationMemorandum COSOR 95-19, 1995, Eindhoven University of Technology
Memorandum COSOR 95-19, 1995, Eindhoven University of Technology A polyhedral approach to the delivery man problem C.A. van Eijl Abstract We propose a mixed integer programming formulation for the delivery
More informationSolving Elementary Shortest-Path Problems as Mixed-Integer Programs
Gutenberg School of Management and Economics Discussion Paper Series Solving Elementary Shortest-Path Problems as Mixed-Integer Programs Michael Drexl and Stefan Irnich Januar 2012 Discussion paper number
More informationCompact Formulations of the Steiner Traveling Salesman Problem and Related Problems
Compact Formulations of the Steiner Traveling Salesman Problem and Related Problems Adam N. Letchford Saeideh D. Nasiri Dirk Oliver Theis March 2012 Abstract The Steiner Traveling Salesman Problem (STSP)
More informationSolving the Asymmetric Travelling Salesman Problem with time windows by branch-and-cut
Math. Program., Ser. A 90: 475 506 (2001) Digital Object Identifier (DOI) 10.1007/s101070100218 Norbert Ascheuer Matteo Fischetti Martin Grötschel Solving the Asymmetric Travelling Salesman Problem with
More informationBranch and Price for the Vehicle Routing Problem with Discrete Split Deliveries and Time Windows
Branch and Price for the Vehicle Routing Problem with Discrete Split Deliveries and Time Windows Matteo Salani Ilaria Vacca December 24, 2009 Report TRANSP-OR 091224 Transport and Mobility Laboratory Ecole
More informationThe Probabilistic Pickup and Delivery Problem
The Probabilistic Pickup and Delivery Problem Enrique Benavent 1, M. Landete 2, J.J. Salazar 3, G. Tirado 4 1 Departament d'estadística i Investigació Operativa, Universitat de València, Spain 2 Universidad
More informationTRAVELING SALESMAN PROBLEM WITH TIME WINDOWS (TSPTW)
TRAVELING SALESMAN PROBLEM WITH TIME WINDOWS (TSPTW) Aakash Anuj 10CS30043 Surya Prakash Verma 10AE30026 Yetesh Chaudhary 10CS30044 Supervisor: Prof. Jitesh Thakkar TSP Given a list of cities and the distances
More informationStronger Multi-Commodity Flow Formulations of the (Capacitated) Sequential Ordering Problem
Stronger Multi-Commodity Flow Formulations of the (Capacitated) Sequential Ordering Problem Adam N. Letchford Juan-José Salazar-González May 2015 Abstract The sequential ordering problem (SOP) is the generalisation
More informationVehicle Routing and MIP
CORE, Université Catholique de Louvain 5th Porto Meeting on Mathematics for Industry, 11th April 2014 Contents: The Capacitated Vehicle Routing Problem Subproblems: Trees and the TSP CVRP Cutting Planes
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 informationBranching Rules for Minimum Congestion Multi- Commodity Flow Problems
Clemson University TigerPrints All Theses Theses 8-2012 Branching Rules for Minimum Congestion Multi- Commodity Flow Problems Cameron Megaw Clemson University, cmegaw@clemson.edu Follow this and additional
More informationInteger Programming ISE 418. Lecture 8. Dr. Ted Ralphs
Integer Programming ISE 418 Lecture 8 Dr. Ted Ralphs ISE 418 Lecture 8 1 Reading for This Lecture Wolsey Chapter 2 Nemhauser and Wolsey Sections II.3.1, II.3.6, II.4.1, II.4.2, II.5.4 Duality for Mixed-Integer
More informationTechnische Universität München, Zentrum Mathematik Lehrstuhl für Angewandte Geometrie und Diskrete Mathematik. Combinatorial Optimization (MA 4502)
Technische Universität München, Zentrum Mathematik Lehrstuhl für Angewandte Geometrie und Diskrete Mathematik Combinatorial Optimization (MA 4502) Dr. Michael Ritter Problem Sheet 1 Homework Problems Exercise
More informationCombinatorial optimization problems
Combinatorial optimization problems Heuristic Algorithms Giovanni Righini University of Milan Department of Computer Science (Crema) Optimization In general an optimization problem can be formulated as:
More informationPartial Path Column Generation for the Elementary Shortest Path Problem with Resource Constraints
Partial Path Column Generation for the Elementary Shortest Path Problem with Resource Constraints Mads Kehlet Jepsen & Bjørn Petersen Department of Computer Science, University of Copenhagen Universitetsparken
More informationInteger Linear Programming (ILP)
Integer Linear Programming (ILP) Zdeněk Hanzálek, Přemysl Šůcha hanzalek@fel.cvut.cz CTU in Prague March 8, 2017 Z. Hanzálek (CTU) Integer Linear Programming (ILP) March 8, 2017 1 / 43 Table of contents
More informationTransformations of node-balanced routing problems
Transformations of node-balanced routing problems Antonio Martinez-Sykora Tolga Bektaş 1 Southampton Business School Centre for Operational Research, Management Science and Information Systems (CORMSIS)
More informationThe Multiple Traveling Salesman Problem with Time Windows: Bounds for the Minimum Number of Vehicles
The Multiple Traveling Salesman Problem with Time Windows: Bounds for the Minimum Number of Vehicles Snežana Mitrović-Minić Ramesh Krishnamurti School of Computing Science, Simon Fraser University, Burnaby,
More informationEquivalence of an Approximate Linear Programming Bound with the Held-Karp Bound for the Traveling Salesman Problem
Equivalence of an Approximate Linear Programming Bound with the Held-Karp Bound for the Traveling Salesman Problem Alejandro Toriello H. Milton Stewart School of Industrial and Systems Engineering Georgia
More informationLecture 8: Column Generation
Lecture 8: Column Generation (3 units) Outline Cutting stock problem Classical IP formulation Set covering formulation Column generation A dual perspective Vehicle routing problem 1 / 33 Cutting stock
More informationApplied Integer Programming: Modeling and Solution
Applied Integer Programming: Modeling and Solution Chen, Batson, Dang Section 6. - 6.3 Blekinge Institute of Technology April 5, 05 Modeling Combinatorical Optimization Problems II Traveling Salesman Problem
More 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 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 informationExact and Heuristic Algorithms for the Symmetric and Asymmetric Vehicle Routing Problem with Backhauls
Exact and Heuristic Algorithms for the Symmetric and Asymmetric Vehicle Routing Problem with Backhauls Paolo Toth, Daniele Vigo ECCO IX - Dublin 1996 Exact and Heuristic Algorithms for VRPB 1 Vehicle Routing
More informationCapacitated ring arborescence problems with profits
Capacitated ring arborescence problems with profits Alessandro Hill 1, Roberto Baldacci 2, and Edna Ayako Hoshino 3 1 School of Business Universidad Adolfo Ibáñez Diagonal Las Torres 2640, Santiago, Chile
More informationDecomposition-based Methods for Large-scale Discrete Optimization p.1
Decomposition-based Methods for Large-scale Discrete Optimization Matthew V Galati Ted K Ralphs Department of Industrial and Systems Engineering Lehigh University, Bethlehem, PA, USA Départment de Mathématiques
More informationA generalized formulation for vehicle routing problems
Working paper A generalized formulation for vehicle routing problems Pedro Munari a, Twan Dollevoet b, Remy Spliet b a Production Engineering Department, Federal University of São Carlos, Brazil b Econometric
More informationAn Integer Cutting-Plane Procedure for the Dantzig-Wolfe Decomposition: Theory
An Integer Cutting-Plane Procedure for the Dantzig-Wolfe Decomposition: Theory by Troels Martin Range Discussion Papers on Business and Economics No. 10/2006 FURTHER INFORMATION Department of Business
More informationCS/COE
CS/COE 1501 www.cs.pitt.edu/~nlf4/cs1501/ P vs NP But first, something completely different... Some computational problems are unsolvable No algorithm can be written that will always produce the correct
More informationThe traveling salesman problem
Chapter 58 The traveling salesman problem The traveling salesman problem (TSP) asks for a shortest Hamiltonian circuit in a graph. It belongs to the most seductive problems in combinatorial optimization,
More informationVEHICLE ROUTING ON REAL ROAD NETWORKS
VEHICLE ROUTING ON REAL ROAD NETWORKS by Saeideh Dehghan Nasiri SUBMITTED FOR THE DEGREE OF DOCTOR OF PHILOSOPHY STOR-I DOCTORAL TRAINING CENTER LANCASTER UNIVERSITY LANCASTER, UK, SEPTEMBER 2014 c Copyright
More informationto work with) can be solved by solving their LP relaxations with the Simplex method I Cutting plane algorithms, e.g., Gomory s fractional cutting
Summary so far z =max{c T x : Ax apple b, x 2 Z n +} I Modeling with IP (and MIP, and BIP) problems I Formulation for a discrete set that is a feasible region of an IP I Alternative formulations for the
More information16.410/413 Principles of Autonomy and Decision Making
6.4/43 Principles of Autonomy and Decision Making Lecture 8: (Mixed-Integer) Linear Programming for Vehicle Routing and Motion Planning Emilio Frazzoli Aeronautics and Astronautics Massachusetts Institute
More informationWeek Cuts, Branch & Bound, and Lagrangean Relaxation
Week 11 1 Integer Linear Programming This week we will discuss solution methods for solving integer linear programming problems. I will skip the part on complexity theory, Section 11.8, although this is
More informationComplexity of Routing Problems with Release Dates and Deadlines
Complexity of Routing Problems with Release Dates and Deadlines Alan Erera, Damian Reyes, and Martin Savelsbergh H. Milton Stewart School of Industrial and Systems Engineering Georgia Institute of Technology
More informationON COST MATRICES WITH TWO AND THREE DISTINCT VALUES OF HAMILTONIAN PATHS AND CYCLES
ON COST MATRICES WITH TWO AND THREE DISTINCT VALUES OF HAMILTONIAN PATHS AND CYCLES SANTOSH N. KABADI AND ABRAHAM P. PUNNEN Abstract. Polynomially testable characterization of cost matrices associated
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 informationScheduling and Optimization Course (MPRI)
MPRI Scheduling and optimization: lecture p. /6 Scheduling and Optimization Course (MPRI) Leo Liberti LIX, École Polytechnique, France MPRI Scheduling and optimization: lecture p. /6 Teachers Christoph
More informationIntroduction to Integer Linear Programming
Lecture 7/12/2006 p. 1/30 Introduction to Integer Linear Programming Leo Liberti, Ruslan Sadykov LIX, École Polytechnique liberti@lix.polytechnique.fr sadykov@lix.polytechnique.fr Lecture 7/12/2006 p.
More information3.7 Strong valid inequalities for structured ILP problems
3.7 Strong valid inequalities for structured ILP problems By studying the problem structure, we can derive strong valid inequalities yielding better approximations of conv(x ) and hence tighter bounds.
More informationMVE165/MMG630, Applied Optimization Lecture 6 Integer linear programming: models and applications; complexity. Ann-Brith Strömberg
MVE165/MMG630, Integer linear programming: models and applications; complexity Ann-Brith Strömberg 2011 04 01 Modelling with integer variables (Ch. 13.1) Variables Linear programming (LP) uses continuous
More informationModeling with Integer Programming
Modeling with Integer Programg Laura Galli December 18, 2014 We can use 0-1 (binary) variables for a variety of purposes, such as: Modeling yes/no decisions Enforcing disjunctions Enforcing logical conditions
More informationOptimizing departure times in vehicle routes
Optimizing departure times in vehicle routes A.L. Kok, E.W. Hans, J.M.J. Schutten Operational Methods for Production and Logistics, University of Twente, P.O. Box 217, 7500AE, Enschede, Netherlands Abstract
More informationChapter 3 THE TIME-DEPENDENT SHORTEST PAIR OF DISJOINT PATHS PROBLEM: COMPLEXITY, MODELS, AND ALGORITHMS
Chapter 3 THE TIME-DEPENDENT SHORTEST PAIR OF DISJOINT PATHS PROBLEM: COMPLEXITY, MODELS, AND ALGORITHMS 3.1 Introduction The general time-dependent shortest pair of disjoint paths problem (TD-2SP) can
More informationBranch-and-Price for Routing with Probabilistic Customers
Branch-and-Price for Routing with Probabilistic Customers Felipe Lagos 1, Mathias A. Klapp 2, Alejandro Toriello 1 1 H. Milton Stewart School of Industrial and Systems Engineering Georgia Institute of
More informationBranch-and-Bound. Leo Liberti. LIX, École Polytechnique, France. INF , Lecture p. 1
Branch-and-Bound Leo Liberti LIX, École Polytechnique, France INF431 2011, Lecture p. 1 Reminders INF431 2011, Lecture p. 2 Problems Decision problem: a question admitting a YES/NO answer Example HAMILTONIAN
More informationFormulations and Algorithms for Minimum Connected Dominating Set Problems
Formulations and Algorithms for Minimum Connected Dominating Set Problems Abilio Lucena 1 Alexandre Salles da Cunha 2 Luidi G. Simonetti 3 1 Universidade Federal do Rio de Janeiro 2 Universidade Federal
More informationThe Separation Problem for Binary Decision Diagrams
The Separation Problem for Binary Decision Diagrams J. N. Hooker Joint work with André Ciré Carnegie Mellon University ISAIM 2014 Separation Problem in Optimization Given a relaxation of an optimization
More informationPOLYNOMIAL MILP FORMULATIONS
POLYNOMIAL MILP FORMULATIONS Miller-Tucker-Zemlin (J. ACM, 1960); Gavish-Graves (MIT Tech. Report 1978) Fox-Gavish-Graves (Operations Research 1980); Wong (IEEE Conference, 1980); Claus (SIAM J. on Algebraic
More informationCable Trench Problem
Cable Trench Problem Matthew V Galati Ted K Ralphs Joseph C Hartman magh@lehigh.edu Department of Industrial and Systems Engineering Lehigh University, Bethlehem, PA Cable Trench Problem p.1 Cable Trench
More informationDecision Procedures An Algorithmic Point of View
An Algorithmic Point of View ILP References: Integer Programming / Laurence Wolsey Deciding ILPs with Branch & Bound Intro. To mathematical programming / Hillier, Lieberman Daniel Kroening and Ofer Strichman
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 informationTechnische Universität Dresden Institute of Numerical Mathematics
Technische Universität Dresden Institute of Numerical Mathematics An Improved Flow-based Formulation and Reduction Principles for the Minimum Connectivity Inference Problem Muhammad Abid Dar Andreas Fischer
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 informationInternational ejournals
ISSN 0976 1411 Available online at www.internationalejournals.com International ejournals International ejournal of Mathematics and Engineering 2 (2017) Vol. 8, Issue 1, pp 11 21 Optimization of Transportation
More informationWorld Academy of Science, Engineering and Technology International Journal of Computer and Information Engineering Vol:5, No:4, 2011
A New Integer Programming Formulation for the Chinese Postman Problem with Time Dependent Travel Times Jinghao Sun, Guozhen Tan, and Guangjian Hou Abstract The Chinese Postman Problem (CPP) is one of the
More informationInteger Linear Programming
Integer Linear Programming Solution : cutting planes and Branch and Bound Hugues Talbot Laboratoire CVN April 13, 2018 IP Resolution Gomory s cutting planes Solution branch-and-bound General method Resolution
More informationDecision Diagrams for Discrete Optimization
Decision Diagrams for Discrete Optimization Willem Jan van Hoeve Tepper School of Business Carnegie Mellon University www.andrew.cmu.edu/user/vanhoeve/mdd/ Acknowledgments: David Bergman, Andre Cire, Samid
More informationLogic, Optimization and Data Analytics
Logic, Optimization and Data Analytics John Hooker Carnegie Mellon University United Technologies Research Center, Cork, Ireland August 2015 Thesis Logic and optimization have an underlying unity. Ideas
More informationVehicle Routing and Scheduling. Martin Savelsbergh The Logistics Institute Georgia Institute of Technology
Vehicle Routing and Scheduling Martin Savelsbergh The Logistics Institute Georgia Institute of Technology Vehicle Routing and Scheduling Part II: Algorithmic Enhancements Handling Practical Complexities
More informationHybrid Optimization Methods for Time-Dependent Sequencing Problems
Hybrid Optimization Methods for Time-Dependent Sequencing Problems Joris Kinable a,b,, Andre A. Cire c, Willem-Jan van Hoeve a a Tepper School of Business, Carnegie Mellon University 5000 Forbes Ave.,
More informationAn exact algorithm for Aircraft Landing Problem
An exact algorithm for Aircraft Landing Problem Min Wen Jesper Larsen Jens Clausen Informatics and Mathematical Modelling Technical University of Denmark 2800 Kgs. Lyngby Denmark September 4, 2005 Abstract
More informationHow to Relax. CP 2008 Slide 1. John Hooker Carnegie Mellon University September 2008
How to Relax Slide 1 John Hooker Carnegie Mellon University September 2008 Two ways to relax Relax your mind and body. Relax your problem formulations. Slide 2 Relaxing a problem Feasible set of original
More informationModels and Algorithms for Stochastic and Robust Vehicle Routing with Deadlines
Accepted in Transportation Science manuscript (Please, provide the mansucript number!) Authors are encouraged to submit new papers to INFORMS journals by means of a style file template, which includes
More informationMinimizing Total Delay in Fixed-Time Controlled Traffic Networks
Minimizing Total Delay in Fixed-Time Controlled Traffic Networks Ekkehard Köhler, Rolf H. Möhring, and Gregor Wünsch Technische Universität Berlin, Institut für Mathematik, MA 6-1, Straße des 17. Juni
More information2 Notation and Preliminaries
On Asymmetric TSP: Transformation to Symmetric TSP and Performance Bound Ratnesh Kumar Haomin Li epartment of Electrical Engineering University of Kentucky Lexington, KY 40506-0046 Abstract We show that
More informationSeparating Simple Domino Parity Inequalities
Separating Simple Domino Parity Inequalities Lisa Fleischer Adam Letchford Andrea Lodi DRAFT: IPCO submission Abstract In IPCO 2002, Letchford and Lodi describe an algorithm for separating simple comb
More informationA Gossip Algorithm for Heterogeneous Multi-Vehicle Routing Problems
A Gossip Algorithm for Heterogeneous Multi-Vehicle Routing Problems Mauro Franceschelli Daniele Rosa Carla Seatzu Francesco Bullo Dep of Electrical and Electronic Engineering, Univ of Cagliari, Italy (e-mail:
More information