A column-and-row generation approach for the flying sidekick travelling salesman problem

https://doi.org/10.1016/j.trc.2020.102913Get rights and content

Highlights

Abstract

The flying sidekick traveling salesman problem (FS-TSP) is a parcel delivery problem arising in the last-mile logistics, where the distribution plan of a driver-operated truck assisted by a drone (unmanned aerial vehicle, UAV) has to be defined. The FS-TSP is a variant of the TSP where routing decisions are integrated with customer-to-drone and customer-to-truck assignment decisions and truck-and-drone synchronization constraints. The objective is the minimization of the time required to serve all the customers, taking into account drone payload capacity and battery power constraints. In this work we provide a new representation of the FS-TSP based on the definition of an extended graph. This representation allows to model the problem by a new and compact integer linear programming formulation, where the synchronization issue is tackled in a column generation fashion, thus avoiding the usage of big-M constraints, representing one of the main drawbacks of the models present in literature. The proposed formulation has been solved by an exact approach which combines a Branch-and-Cut algorithm and a column generation procedure, strengthened by variable fixing strategies and new valid inequalities specifically defined for the problem. The proposed method has been experienced on a large set of benchmark instances. Computational results show that the proposed approach either is competitive or outperforms the best exact approach present in literature for the FS-TSP. Indeed, it is able to provide the optimal solution for all small size instances with 10 customers and for several medium size instances with 20 customers, some of them never solved before.

Introduction

In the last years, the attention of the scientific research community towards the usage of drones has continuously increased (Otto et al., 2018). The motivations behind this growing research interest can be mainly found in the continuous developments of drone technologies and in the subsequent exploration of new drone applications, driven by large corporations such as Amazon, Google, UPS, and DHL (Bezos, 2013, Cary and Bose, 2016, Reuters staff, 2018, Reuters staff, 2015, UPS staff, 2017). Among the others, great attention has been paid to the usage of drones in surveillance and security, disaster management and transportation fields (Shakhatreh et al., 2019). This has led to the definition of new and complex decision and optimization problems for which operations research methodologies represent a valuable support tool.

In this context, several studies have been made on the usage of drones in the last-mile logistic, showing the benefits in terms of emissions and completion time reduction that can be achieved in parcel delivery operations (Kille et al., 2019, Rojas Viloria et al., 2020).

The most studied drone assisted logistic system consists of a single truck and a single drone (or unmanned aerial vehicle, UAV) operating in tandem for the parcel delivery to the customers. In a nutshell, the system performs as follows: the truck and the drone must depart and return to a single depot, either in tandem or independently; each customer must be served once either by the truck or by the drone; the truck has an infinite capacity and acts as a mobile depot for the drone, replenishing its batteries and providing the parcels to be delivered. A more formal and detailed description of the system features, assumptions and operations will be provided in Section 3.

To the best of authors’ knowledge, the first work dealing with this hybrid working unit (truck-and-drone) is the one presented in Murray and Chu (2015), where it is referred to as Flying Sidekick Traveling Salesman Problem (FS-TSP). The FS-TSP is a variant of the TSP where routing decisions are integrated with customer-to-drone and customer-to-truck assignment decisions and truck-and-drone synchronization constraints. The objective is the minimization of the time required to serve all the customers (completion time), taking into account drone payload capacity and battery power constraints.

Several contributions have appeared in the last years providing integer and mixed-integer linear programming (ILP and MILP) formulations for the FS-TSP and modifications of it, where different objective functions and/or operating conditions are considered. However, most of these formulations suffer from dimensional drawbacks which make their solution impracticable even on small instances. This is mainly due to the synchronization issue which is generally modeled by the usage of big-M constraints. Indeed, these constraints make the usage of off-the-shelf optimization software ineffective and they would require the development of ad-hoc exact row and/or column generation solution approaches.

In this context, as will be discussed in the review provided in Section 2, even if the majority of the literature contributions present formulations for the specific truck-and-drone problem under investigation, then they mainly focus on heuristic solving approaches able to find good sub-optimal solutions with an acceptable computational burden.

The aims of this work is filling this gap by proposing an exact Branch-and-Cut algorithm coupled with a column generation procedure for the FS-TSP presented in Murray and Chu (2015).

Our method exploits a new representation of the FS-TSP based on the definition of an extended graph. This representation allows to model the problem by a new and compact ILP formulation, where the synchronization issue is tackled in a column generation fashion, thus avoiding the usage of big-M constraints. Moreover, the proposed method is strengthened by the usage of variable fixing strategies and new valid inequalities specifically defined for the problem.

To the best of our knowledge, the state of the art of exact approaches for the FS-TSP is represented by the work of Dell’Amico et al. (2019b). Our Branch-and-Cut algorithm has been experienced on the same set of benchmark instances used in Dell’Amico et al. (2019b) and the obtained computational results show that our approach either is competitive or outperforms it. Indeed, it is able to provide the optimal solution for all small size instances with 10 customers and several medium size instances with 20 customers, some of them never solved before. We underline that the maximum size of the tackled instances is coherent with the potential usage of the described truck-and-drone system. Indeed, such system can easily find real applications in low density inhabited (rural) areas due to the present flight regulation. In these areas the time required to travel distances represents the main part of the courier shift, so reducing the number of customers that can be served per day. In Karcz and Slusarczyk (2016) a real example of parcel delivery in rural areas is described and the same order of magnitude of customers is considered.

The paper is organized as follows: in Section 2, we discuss upon the most relevant studies focusing on exact solution methods for delivery problems with drones; in Section 3, we provide a detailed description of the FS-TSP and we present its extended graph representation; Section 4 is devoted to the presentation of the proposed ILP formulation; Section 5 introduces several variable fixing strategies and families of valid inequalities to strengthen the proposed model; Section 6 describes our exact approach to solve the problem; Section 7 proves the effectiveness of the proposed method showing the obtained computational results; finally, conclusions are given and future perspectives are discussed in Section 8.

Section snippets

Literature review

A large body of recent transportation science and operations research literature has focused on the civil applications of drones, with particular interest in logistics and parcel delivery fields. For a comprehensive survey of the main research issues and deriving optimization problems, we address the reader to the survey works by Chung et al., 2020, Macrina et al., 2020, Otto et al., 2018 and Rojas Viloria et al. (2020).

In this section, consistently with the scope of this paper, we discuss only

Problem description

The FS-TSP can be seen as a particular variant of the TSP. Solving a TSP instance requires to decide the order according to which the customers will be served by a truck with the aim of minimizing the total route length. On the other hand, solving an FS-TSP instance requires to decide:

  • the route to be performed by the truck;

  • the set of clients served by the drone taking into account its payload capacity and battery power;

  • the nodes to be used as launch and pick-up (rendezvous) nodes when the drone

Problem formulation

On the basis of the notation previously introduced, we can model the FS-TSP by an ILP formulation where the following three sets of binary variables are used:

  • yij,(i,j)AT, is equal to 1 if the arc (i,j) belongs to the truck path, 0 otherwise.

  • xij,(i,j)AD, is equal to 1 if the arc (i,j) belongs to the drone path, 0 otherwise.

  • zPijk,kD,i,jT{s,t}andPijΠij such that WPijk>0, is equal to 1 if the drone travels the path i-k-j and the truck travels the path Pij, 0 otherwise. In other words, we

Variable fixing and valid inequalities for the FS-TSP

In this section we present several variable fixings and valid inequalities which will be integrated in our solution approach in order to improve its effectiveness. The variable fixings are implemented to reduce the solution space. The valid inequalities are introduced to strengthen the proposed formulation.

A branch-and-cut approach

We developed a row-and-column generation procedure and we embedded it into a Branch-and-Cut framework to optimally solve the FS-TSP. In particular, the column generation procedure allows us to avoid a complete enumeration of the z-variables, which would be impracticable. Likewise, as widely known, the row generation procedure is required for all the constraints which are exponential in the number and for which a separation procedure is required.

The Branch-and-Cut algorithm considers the

Computational results

In this section we present and discuss the computational results of the experimentation performed to evaluate and validate the proposed Branch-and-Cut algorithm (B&C). The experiments have been performed on an Intel(R) Core(TM) i7-6500U, 2.50 GHz, 8.00 GB of RAM. The B&C algorithm has been coded in C language using Cplex 12.7 with default setting as MILP solver, imposing a computation time limit of 1 h. Two test beds have been considered:

  • The set of 72 instances proposed in Murray and Chu (2015)

Conclusions

In this study, we proposed a new and compact ILP formulation for the FS-TSP, which unlike the ones present in the literature, overcome the drawbacks of the big-M constraints required to model the synchronization issue. We developed a row-and-column generation procedure and we embedded it into a Branch-and-Cut algorithm. Moreover, we presented variable fixing strategies and original set of valid inequalities which further strengthen the algorithm. Computational results confirmed the

CRediT authorship contribution statement

Maurizio Boccia: Conceptualization, Methodology, Software, Validation, Writing - original draft, Writing - review & editing. Adriano Masone: Conceptualization, Methodology, Software, Validation, Writing - original draft, Writing - review & editing. Antonio Sforza: Conceptualization, Methodology, Software, Validation, Writing - original draft, Writing - review & editing. Claudio Sterle: Conceptualization, Methodology, Software, Validation, Writing - original draft, Writing - review & editing.

References (42)

  • C.C. Murray et al.

    The flying sidekick traveling salesman problem: Optimization of drone-assisted parcel delivery

    Transp. Res. Part C: Emerg. Technol.

    (2015)
  • C.C. Murray et al.

    The multiple flying sidekicks traveling salesman problem: Parcel delivery with multiple drones

    Transp. Res. Part C: Emerg. Technol.

    (2020)
  • S. Poikonen et al.

    Multi-visit drone routing problem

    Comput. Oper. Res.

    (2020)
  • D. Sacramento et al.

    An adaptive large neighborhood search metaheuristic for the vehicle routing problem with drones

    Transp. Res. Part C: Emerg. Technol.

    (2019)
  • M. Salama et al.

    Joint optimization of customer location clustering and drone-based routing for last-mile deliveries

    Transp. Res. Part C: Emerg. Technol.

    (2020)
  • D. Schermer et al.

    A hybrid vns/tabu search algorithm for solving the vehicle routing problem with drones and en route operations

    Comput. Oper. Res.

    (2019)
  • Z. Wang et al.

    Vehicle routing problem with drones

    Transp. Res. Part B: Methodol.

    (2019)
  • N. Agatz et al.

    Optimization approaches for the traveling salesman problem with drone

    Transp. Sci.

    (2018)
  • Bezos, J., 2013. Interview with Charlie Rose. broadcasted on: 60 minutes, cbs. (transcript)....
  • P. Bouman et al.

    Dynamic programming approaches for the traveling salesman problem with drone

    Networks

    (2018)
  • Cary, N., Bose, N., 2016. UPS, FedEx and Amazon gather flight data to prove drone safety....
  • Cited by (0)

    View full text