Discrete Optimization
A k-product uncapacitated facility location problem

https://doi.org/10.1016/j.ejor.2007.01.010Get rights and content

Abstract

A k-product uncapacitated facility location problem can be described as follows. There is a set of demand points where clients are located and a set of potential sites where facilities of unlimited capacities can be set up. There are k different kinds of products. Each client needs to be supplied with k kinds of products by a set of k different facilities and each facility can be set up to supply only a distinct product with a non-negative fixed cost determined by the product it intends to supply. There is a non-negative cost of shipping goods between each pair of locations. These costs are assumed to be symmetric and satisfy the triangle inequality. The problem is to select a set of facilities to be set up and their designated products and to find an assignment for each client to a set of k facilities so that the sum of the setup costs and the shipping costs is minimized. In this paper, an approximation algorithm within a factor of 2k+1 of the optimum cost is presented. Assuming that fixed setup costs are zero, we give a 2k-1 approximation algorithm for the problem. In addition we show that for the case k=2, the problem is NP-complete when the cost structure is general and there is a 2-approximation algorithm when the costs are symmetric and satisfy the triangle inequality. The algorithm is shown to produce an optimal solution if the 2-product uncapacitated facility location problem with no fixed costs happens to fall on a tree graph.

Introduction

In the classical simple plant location or 1-product uncapacitated facility location problem, we have to select a set of facilities to be set up and a set of clients for each facility to service so as to minimize the total cost of setting up the facilities and servicing the clients. In the last few years, a number of constant factor approximation algorithms have been proposed for this problem when the service cost is assumed to be in the metric space by Aardal et al., 1999, Charikar et al., 1999, Charikar and Guha, 1999, Guha and Khuller, 1999, Jain and Vazirani, 1999, Shmoys et al., 1997. That is, the service cost is assumed to be symmetric and satisfy the triangle inequality.

The first heuristic algorithm with a performance guarantee of 31-e-3 was given by Shmoys et al. (1997), which is based on a linear program (LP) rounding algorithm extended from the filter technique of Lin and Vitter, 1992a, Lin and Vitter, 1992b. Coupling with a local search phase with the LP rounding, Guha and Khuller (1999) improved the factor to be 2.408. Later Chudak and Shmoys (2003) further strengthened the LP rounding approach to obtain a 1+2e-approximation algorithm. This idea has since been extended by Aardal et al. (1999) for the k-level facility location problem. Another interesting and elegant approach to obtain a constant factor approximation is via the primal–dual algorithm proposed by Jain and Vazirani (1999). This idea is extended to the two-level facility location problem where a performance guarantee of 1.77 was obtained by Zhang (2004). Jain et al. (2002) gave a simple greedy algorithm for the 1-product facility location problem with a performance guarantee of 1.61. This ratio was improved to 1.52 by Mahdian et al. (2002), which is close to the lower bound 1.463 proved by Guha and Khuller (1999).

In practice, a client usually requires more than one product to be supplied by several sources. This problem can be referred as multiproduct uncapacitated facility location problem or k-product uncapacitated facility location problem (k-PUFLP). It was mentioned as a special case of the submodular set function by Fisher et al. (1978). Various optimal and heuristic algorithms for the k-PUFLP have previously been considered in the general space by Klincewicz et al., 1986, Klincewicz and Luss, 1987.

For the multiproduct capacitated facility location problem, it has been studied by Lee, 1991, Lee, 1993, Mazzola and Neebe, 1999.

In this paper we consider the k-PUFLP and present heuristic algorithms for a variety of k-PUFLPs in the metric space. In Section 2, we formulate the k-PUFLP. In Section 3, we consider the k-product uncapacitated facility location problem with no fixed costs (k-PUFLPN). We show that the 2-PUFLPN is NP-complete when the cost structure is assumed to be general without the metric space property and a heuristic algorithm with a performance guarantee of 2 is presented when the costs are assumed to satisfy the metric space property. It is shown that the algorithm produces an optimal solution if the 2-PUFLPN happens to fall on a tree graph. For the k-PUFLPN (k3), we present a heuristic algorithm with a performance guarantee of 2k-1. In Section 4, we give a heuristic algorithm with a performance guarantee of 2k+1 for the k-PUFLP when k is fixed. In Section 5 we conclude the paper by suggesting some questions for future research.

Section snippets

The formulation of the k-PUFLP

Let D be the set of clients and F be the set of potential facilities. There are k kinds of products, pl,l=1,2,,k. Each facility iF may be set up to provide at most one of the products. The cost of setting up a facility i to supply product pl is fil, iF,1lk. The cost of shipping between any two points i,jFD is equal to cij. Each client jD must be supplied with k products by a set of k facilities. In other words, split sourcing is not allowed for a given product.

Throughout this paper, we

Algorithms for the 2-PUFLPN and the k-PUFLPN (k3)

Suppose that all of the fixed costs are zero, i.e., fil=0,iF,l{1,2,,k}. We refer to it as a k-PUFLPN.

It is clear that all facilities can be assumed to be set up as no cost will be incurred. Thus constraints (4) become l=1kyil=1. The following is the underlying LP of this problem and its optimal value provides a lower bound for the problem.(Pk)minl=1kiFjDcijxijl,subject toiFxijl=1,jD,l=1,2,,k;xijlyil,iF,jD,l=1,2,,k;l=1kyil=1,iF,xijl0,iF,jD,l{1,2,,k},yil0,iF,l{1,2,,

Algorithm for the k-PUFLP

In the following, we consider the k-PUFLP. In this section, we suggest another formulation for k-PUFLP. For convenience of presentation we use s to represent a sequence of k different facilities ilF,l=1,2,,k, and refer to s as a feasible sequence of facilities. The set of all possible feasible sequences is denoted by S. Hence each client jD must be assigned to exactly one feasible sequence s=(i1,i2,,ik)S. For simplicity, a client jD is said to be assigned to a feasible sequence s=(i1,i2,,

Discussion

In this paper we consider the k-PUFLP and suggest a few approximation algorithms. One interesting question that remains open is whether our Algorithm 1 for the 2-PUFLPN on tree graphs can be extended to the case when k3. Another interesting question for future research is whether there exist approximation algorithms for the k-PUFLP and the k-PUFLPN with constant performance guarantee independent of k.

References (21)

There are more references available in the full text version of this article.

Cited by (15)

  • Location selection of city logistics centers under sustainability

    2015, Transportation Research Part D: Transport and Environment
    Citation Excerpt :

    Since the 1970s, the literature is replete with work on location theory. We name a few recent key works for completeness, for example, the hybrid multi-start heuristic (Resende and Werneck, 2006), second-order cone programming (Wagner et al., 2009), approximation algorithms (Huang and Li, 2008; Li, 2013), greedy heuristic and fix-and-optimize heuristic (Ghaderi and Jabalameli, 2013), Lagrangian relaxation heuristic (Nezhad et al., 2013), mixed integer linear programming model (Kratica et al., 2014), discrete variant of unconscious search (Ardjmand et al., 2014), multi-objective optimization model (Tang et al., 2013), and the weighted Dantzig–Wolfe decomposition and path-relinking combined method (Li et al., 2014), which have been presented for solving an uncapacitated facility location problem. Also, some algorithms and methods have been proposed for solving the capacitated facility location problem to optimality such as the mixed integer programming formulation (Melkote and Daskin, 2001; Aros-Vera et al., 2013; Rosa et al., 2014), branch-and-price algorithm (Klose and Görtz, 2007), Lagrangian heuristic algorithm (Wu et al., 2006; Elhedhli and Merrick, 2012), kernel search heuristic (Guastaroba and Speranza, 2014), Lagrangian Heuristic and Ant Colony System (Chen and Ting, 2008), Lagrangian relaxation algorithm (Yun et al., 2014), a fix-and-optimize heuristic based on the evolutionary fire-fly algorithm (Rahmaniani and Ghaderi, 2013), hybrid Firefly-Genetic Algorithm (Rahmani and MirHassani, 2014), swarm intelligence based on sample average approximation (Aydin and Murat, 2013), modified Clarke and Wright savings heuristic algorithm (Li et al., 2015), iterated tabu search heuristic (Ho, 2015), improved approximation algorithm (Aardal et al., 2015), two-stage robust models and algorithms (An et al., 2014), and the evolutionary multi-objective optimization approach (Rakas et al., 2004; Harris et al., 2014).

  • A Cost-Sharing Scheme for the k-Product Facility Location Game with Penalties

    2024, Lecture Notes in Computer Science (including subseries Lecture Notes in Artificial Intelligence and Lecture Notes in Bioinformatics)
View all citing articles on Scopus
View full text