Abstract

When consumers faced with the choice of competitive chain facilities that offer exclusive services, current rules cannot describe these customers’ behaviors very well. So we propose a partially proportional rule to represent this kind of customer behavior. In addition, the exact demands of customers in many real-world environments are often difficult to determine. This is contradicting to the assumption in most studies of the competitive facility location problem. For the competitive facility location problem with the partially proportional rule, we establish a robust optimization model to handle the uncertainty of customers’ demands. We propose two methods to solve the robust model by studying the properties of the counterpart problem. The first method is presented by solving a mixed-integer optimization model of the counterpart problem directly. The second method is given by embedding a sorting subalgorithm into the simulated annealing framework, in which the sorting subalgorithm can easily solve the subproblem. The effects of the budget and the robust control parameter to the location scheme are analyzed in a quasi-real example. The result shows that changes in the robust control parameter can affect the customer demands that were captured by the new entrants, thereby changing the optimal solution for facility location. In addition, there is a threshold of the robust control parameter for any given budget. Only when the robust control parameter is larger than this threshold, the market share captured by the new entering firm increases with the increases of this parameter. Finally, numerical experiments show the superiority of the algorithm in large-scare competitive facility location problems.

1. Introduction

Competitive facility location (CFL) is the problem of locating new facilities in competitive markets with the aim of maximizing market share [1, 2]. A lot of competitive facility location models have been proposed in recent years because these models are very useful in many competitive situations [3]. There are various competitive facility location models depending on the components to be considered. These components usually include competitive type, location space, and customer behavior (see survey paper [4]).

Suppose that there are different facilities offering similar goods or services, customer behavior refers to the way that a customer how to spend his buying power on these facilities. It is obvious that customer behavior is a key ingredient in the competitive facility location problem [5]. The two most common customer behavior rules employed in literature are the binary rule and the proportional rule [6]. The binary rule dates back to the duopoly model proposed by Hotelling [7]. It is traditionally assumed that customers patronize the nearest facility to be served, but other characteristics of the facilities can also be taken into account [810]. In fact, if we use the parameter (attraction) to uniformly represent the different characteristics considered by competitive facilities, then the essence of the binary rule is that the customers always patronize the most attractive facility. The proportional rule is first proposed by Huff [11], which assumes that the customers patronize all facilities in proportion to their attractions [1216]. Many researchers have studied different customer behavior rules, but most of them can be regarded as variants of the above two rules. Such as in [17], the authors studied the customer behavior with minimal attraction requirement. Drezner et al. [18] investigated the competitive facility location problem based on the concept of “radius of influence.” In their model, although the calculation of the attraction is based on the facility’s cover distance, the division of the market share still relies on the proportional rule. Qi et al. [19] observed that customers usually patronize facilities within a range that they feel is convenient. They proposed a proportional rule with a limited range to study the competitive facility location problem. Fernández et al. [20] studied the customer behavior that customer only patronizes those facilities that the attraction greater than or equal to a threshold value.

According to the classification of Suáez-Vega et al. [21], there is another basic rule of customer behavior called the partially binary rule. Following this rule, the customer first selects the most attractive facility from each firm and then splits his demand among those facilities proportionally to their attractions [22]. The partially binary rule is usually studied together with the binary rule and/or the proportional rule. For example, Suárez-Vega et al. [21] and Biesinger et al. [23] studied the competitive facility location problems in the network space and the discrete space, respectively. Six scenarios have been considered in both papers, they are combinations of two service types (essential and unessential) and three customer behavior rules (binary, proportional, and partially binary). Fernández et al. [24] proposed two new heuristic algorithms to solve competitive facility location problems with the binary rule and the partially binary rule, respectively. Fernández et al. [25] studied a continuous competitive single-facility location and design problem with the partially binary rule. The comparison with the proportional rule reveals their location results are quite different.

Most of the customer behaviors in competitive facility location problems can be described by these three rules. However, when consumers faced with the choice of competitive chain facilities that offer exclusive services, we find that these three rules cannot describe customers’ behaviors well. For example, when somebody wants to apply for a bank account for deposit/withdrew or other services, because each bank has several business service points (including manual service points and ATMs), a customer usually uses the total attraction of all business service points as the evaluation value of the bank. The typical customer behavior is as follows: the customer chooses the most attractive bank to apply a bank account, and then he will patronize the business service points of this bank in proportion to their attractions. From the typical behavior of how customers choose competitive chain facilities that provide exclusive services, we can find that neither of the binary rule, the proportional rule, nor the partially binary rule adequately describes this kind of customer behavior. In this paper, we propose a new kind of customer behavior rule to describe this kind of customer behavior, and we call this rule as the partially proportional rule.

In most studies of the competitive facility location problem, it is usually assumed that the demands of customers are known. However, due to different reasons such as lack of historical data of customers or the high-cost to conduct a comprehensive market survey, the exact demands of customers in many real-world environments are often difficult to determine. To the best of our knowledge, only two literatures involved the uncertainty of the demands of customers. Shiode and Drezner [26] assumed that the weights of demand points are stochastic, and then they studied the competitive facility location problem in the Stackelberg game framework. Beresnev and Melnikov [27] studied the competitive facility location problem in the situation of there are several alternative demand scenarios. The uncertainty of demands in the above two papers is actually expressed by random variables, but it is quite difficult to determine the corresponding random distribution in many cases. Therefore, we use interval data to express the uncertainty of demands and then adopt a robust optimization method to deal with this uncertainty.

The remainder of this paper is organized as follows: Section 2 is devoted to present the definition of the partially proportional rule and the robust optimization method. Section 3 consists the robust model of competitive facility location problem with the partially proportional rule. Section 4 includes two solution methods and to solve the presented model. The effects of the budget and the robust control parameter in the location scheme are presented in Section 5, and the superiority of is also shown in large-scale numerical tests. Finally, some conclusions are presented in Section 6.

2. Problem Definition

2.1. Partially Proportional Rule

In this section, we first introduce a suitable rule to describe the behavior of how customers choose competitive chain facilities that provide exclusive services. Suppose there are several firms in the competitive market, each firm has several chain facilities that offer exclusive services. A customer only chooses one of these services due to the exclusive of these services. If the customer chooses one of these firms, any of the firm’s chain facilities can provide the corresponding exclusive service to the customer. In order to describe well this kind of customer behavior, the partially proportional rule is defined as follows: a customer chooses the most attractive firm from all firms at first and then splits his buy powers on the selected firm’s facilities in proportion to their attractions.

In order to illustrate the difference between the binary rule, the proportional rule, the partially binary rule, and the partially proportional rule, we give an example to show how customer will patronize facilities under different behavior rules. Suppose that there are three firms A, B, and C in the competing market. Each firm has several chain facilities, firm A has facilities 1, 4, 5, and 9, firm B has facilities 2, 6, and 7, and firm C has facilities 3, 8, and 10. For simplicity, suppose that for a customer, the attraction of the facility k is equal to its serial number, i.e., . The patronizing patterns under different customer behavior rules are shown in Figure 1.

If the customer follows the binary rule, he only patronizes the facility 10 because this facility has the largest attraction among all facilities. If the customer follows the proportional rule, he patronizes all the facilities and splits his demands to facility k with the proportion , where 55 is the total attraction of all facilities. If the customer follows the partially binary rule, he selects the most attractive facility from each firm at first, that is, facility 9 of firm A, facility 7 of firm B, and facility 10 of firm C. Then he patronizes facilities 9, 7, and 10 and splits his demands in proportion to 9/26, 7/26, and 10/26, respectively. If the customer follows the partially proportional rule, he compares the total attraction of each firm at first. In this example, the total attraction of firms A, B, and C is 19, 15, and 21, respectively. The most attractive firm is C; hence, he patronizes facilities 3, 8, and 10 of firm C and splits his demands in proportion to 3/21, 8/21, and 10/21, respectively.

Note that for a customer, if there are multiple firms that are most attractive, then the customer patronizes all the facilities of these firms and splits his demands in proportion to their attractions.

2.2. Robust Optimization

A robust model with interval data uncertainty was first proposed by Soyster [28], but his method usually leads to a solution that is too conservative because of always considering the worst case. Ben-Tal and Nemirovsky [29] addressed the over-conservatism by applying an ellipsoid uncertainty set, but the robust counterpart model cannot conserve the complexity of the original problem. Bertsims and Sim [30] presented a method with the flexibility to adjust the degree of conservativeness while maintaining the linear robust counterparts. Before establishing a robust model, we first briefly describe the robust model introduced in [30]. The authors considered the following problem:where x is the n-dimensional decision variable, and is the j-th component of x. Some parameters of the coefficient matrix A are uncertain. For a particular row i of A, let denote the set of coefficients that are uncertain. Suppose each entry takes values in . For every i, a parameter is introduced as the robust control parameter. The goal is to be protected against all cases that up to of these coefficients are allowed to change, and one coefficient changes by . The role of the parameter is to adjust the robustness against the level of conservatism of the solution. Then the robust model of the original problem iswhere

3. Proposed Model

Suppose there are several firms in a region, each providing exclusive services through their chain facilities. A new entering firm wants to set up several chain facilities to offer its exclusive service. The problem of the new entering firm is how to select the locations of its chain facilities to gain the largest market share under a given budget. There is a set of discrete customers that can aggregate to demand points. The demand of each customer is uncertain, but the range of its value is known. The following notations are used:K: set of existing firms, indexed by kI: set of demand points (customers), indexed by i: set of firm k’s chain facilities, k is an existing firmJ: set of existing facilities, i.e., : set of potential facility locations, indexed by j: demand of customer , it is uncertain within the interval : opening cost at potential facility location : distance between customer i and facility j: attraction that customer i feels from facility jG: budget of the new entering firm: binary variable, if facility j is select to open and 0 otherwise

The goal of the new entrant is to maximize the market share. Hence, the competitive facility location problem for the new entering firm is as follows:where indicates the proportion of demand that the new entering firm captured from customer i. , let , and let . If customers follow the partially proportional rule, then can be expressed as

By adding binary variables , and let ϵ be a positive real number that is small enough, condition (5) can be expressed with some linear constraints. Hence, the above competitive facility location model can be transformed into an integer optimization model as follows:

Note that for any , if , then the optimal solution is . In this case, the new entering firm gets the total demand of customer i, i.e., . If , then the optimal solution is , and the proportion of demand that the new entering firm get from the customer i is . If , then the optimal solution is , that is, . The role of ϵ is to compare the size of and ; hence, any positive number that is less than can be used as ϵ.

Let denote the solution set of (7), and , where is a parameter which is used to control the total uncertainty of . Naturally, the uncertainty of can be expressed as with . Similar to the robust model that is proposed in [30], the robust model of the above optimization problem can be expressed as follows:

4. Solution Method

In this section, we propose two methods to solve the robust optimization problem . The first method is based on directly solving an integer model that is equivalent to the origin robust model .

4.1. Method

Given , let us consider the suboptimization problem :

The dual problem of is as follows:

So the robust optimization model can be converted to as follows:

Since the equivalence model is a mixed integer optimization model, the first method is to solve it directly using some commerce software such as CPLEX and GUROBI. We denote this method as Method .

4.2. Method

The second method to solve is based on the simulated annealing framework. The key point is that we found the structure of the optimal solution to the subproblem. Now we first present a Lemma 1.

Lemma 1. For let . Consider the optimization problem : . Let be the largest integer less than or equal to . If is defined as follows, it is the optimal solution of :

Proof. Firstly, we prove that is a necessary condition for the optimal solution of problem . Assuming this conclusion is wrong, i.e., there is . Let . We can construct a new solution as follows: set , and . It is easy to know that , but . This is a contradiction.
Secondly, we prove that there exists an optimal solution, where at most one component takes the value , and the other components are 0 or 1. Assuming this conclusion is wrong, i.e., for every optimal solution, there are at least two components and such that . Without loss of generality, we can assume that . Now we can construct a new solution. Set , and as follows:(1)If , then (2)If , then Obviously, , for case (1) we know that . For case (2), we know that . So is also an optimal solution. This is a contradiction.
Lastly, we prove that the solution defined in (11) is an optimal solution of . Therefore, we know that the optimal solution of (OP) must satisfy , and at most one component takes the value . So we can conclude that there are components equal to 1, and one component equals and all other components are 0. Therefore, according to Chebyshev’s sum inequality [31], we know that defined by (11) is the optimal solution of the problem .
By applying the above Lemma 1 to the problem , we can get the following Theorem 1.

Theorem 1. For given , let . Let . Then is an optimal solution of the problem if it is defined as follows:(1)If , , and (2)If , record the sequence of in the ascending order as , where , then , and for , defined as follows:where .

Proof. The problem can be rewritten as follows:(1)If , from , and , we can get that . Obviously, the objective function corresponding to is 0. Hence, we can conclude that is an optimal solution of in this case.(2)If , firstly, we prove that the optimal solution of satisfies . Assuming this conclusion is wrong, i.e., there is a such that . Let . It is easy to construct a new solution as follows:Set . Find any one such that , let and for other . We can get that . This is a contradiction. Hence, there must be , . We can conclude that is an optimal solution of according to Lemma 1.
Given a facility location solution y, we propose a sorting-based algorithm to solve as follows, this algorithm can be used to calculate the market share that was captured by the new entering firm:SortA:(1)Define , as follows:And define , in a similar way (just substitute by ).(2)Record the sequence of in the ascending order as . Take into account the uncertainty of , the market share of the new entering firm isThe simulated annealing algorithm was introduced by Kirkpatrick et al. [32]. By embedding the sorting-based algorithm (SortA) into the simulated annealing framework, we present the solution method in Algorithm 1.
The process of CreateInitialSolution (ProblemSize) is as follows:
Calculate , by sorting in the descend order. Record the rearranged index sequence as . Find the maximal index such that , then let for and for to generate an initial feasible solution , where .
The process of CreateNeighborSolution () is as follows:

Let , calculate the remaining budget after opening new facilities . Let . Randomly choose two subsets and such that , and then let , , and for other j.

Input: ProblemSize, , , r, L
Output:
 CreateInitialSolution (ProblemSize)
 algorithm SortA
While ()
 For ( To L)
   CreateNeighborSolution ()
   algorithm SortA
  If ()
   
   If ()
    
   End
  Else If (Exp() Rand())
   
  End
 End
End
Return ()

5. Numerical Example

5.1. A Quasi-Real Example

Consider the 49-node data set described in Daskin [33] consisting of the capitals of the continental United States plus Washington, DC. The lower and upper bounds of the customer demands () are proportional to the population of the state (i) in 1890 and 1990, respectively. Hence . Now, there are two firms (A and B) that offer similar exclusive services through their chain facilities. Firm A has 5 facilities located in the states of Alabama, Arizona, Iowa, Maine, and Nebraska, respectively. Firm B has 5 facilities located in the states of New Hampshire, Ohio, South Dakota, Utah, and Wyoming, respectively. The attraction that the customer i feels from facility j is calculated as . The opening cost at the potential facility locations is also taken from the data set of Daskin. Before the new firm C entered the competitive market, the market was completely divided into two parts by firm A and firm B. The division pattern of the market share is shown in Figure 2.

Now firm C wants to enter this market by opening several chain facilities. Any customer point without an existing facility can be regarded as a potential location for new facilities. The goal of firm C is to maximize the market share captured by itself under the constraint of a limited budget G. The robust control parameter can take any value in the interval . For simplicity, takes integer values in our experiments. If the value of is 0 or 49, then is equal to the lower limit or the upper limit , respectively. This means that the demand of each customer has been determined. Hence, we limit the parameter to the integers in the interval .

For , when the parameter , the optimal decision of firm C is to open two facilities in Pennsylvania and Mississippi, respectively. The opening cost of these two facilities is 93000. The location of facilities and the customers that captured by each firm are shown in Figure 3.

However, when , the optimal locations of firm C’s facilities are Pennsylvania and Michigan, respectively. The opening cost of these two facilities is 86,800. Similarly, all information is shown in Figure 4.

For simplicity, we refer to location plans corresponding to and as Plan 1 and Plan 2, respectively. Why does firm C’s location plan change when changes from 42 to 43? In order to find the cause of this change, we list the customer demand ranges captured by firm C in Table 1. The overall contents of Table 1 are divided into two parts: Plan 1 and Plan 2. Take the part of Plan 1 as an example, the first column is the customers captured by firm C, the second column is the demand intervals of customers. The column denotes the optimal solution of μ in when , and the definition of is similar.

Given a facility location plan y, denote the set of customers captured by firm C as . Then the minimal demands captured by firm C can be expressed as . According to the Lemma 1 proposed in Section 4.2, the values of and are easy to obtain. When , the demands captured by firm C with Plan 1 and Plan 2 are 246.0530 and 240.3947, respectively. So in this case, Plan 1 is better than Plan 2. However, when , the demands captured by firm C with Plan 1 and Plan 2 are 283.4438 and 285.7085, respectively. Obviously, Plan 2 is better than Plan 1 in this case. Therefore, we can get a conclusion that changes in the robust control parameter can affect the customer demands that were captured by the new entrants, thereby changing the optimal solution for facility location.

For , the relationships between the demands captured by firm C and the values of are shown in Figure 5. From Figure 5, it is obvious that a larger budget makes the entering firm capture more demands. This is in line with our intuition, as more budgets expand the feasible range and the value of the objective function increases accordingly. For a given budget, there is a threshold for the parameter . Only when the parameter is larger than this threshold, the market share that was captured by the new entering firm increases with the increase of . In addition, the threshold of the robust control parameter decreases as the budget increases. Therefore, uncertainty will have a greater impact on the location of competitive facilities when the budget is larger.

Let . For each , generate scenarios as follows: randomly select number of customers, let their demands to be the corresponding upper bounds, and the demands of the remaining customers to be the corresponding lower bounds. For each , we can solve the deterministic model of the competitive facility location problem to get the optimal location scheme for firm C. Denote the demands captured by firm C in the worst case as . Let be the average of and record the customer demands obtained by firm C based on the robust model as . For each value of , the comparison of and is listed in Table 2.

From Table 2, for all values of . In fact, for the uncertain demands expressed by interval data, the solution obtained by the robust model is better than the solution obtained by other models in the worst case.

The optimal solutions of firm C’s competitive facility location problem are presented in Table 3. The first column denotes the budget G of the firm C, the second column denotes the range of the parameter , the third column denotes the facilities’ opening cost of the firm C, the fourth column denotes the optimal facility locations, and the last column denotes the customers captured by the firm C.

5.2. Performances of Two Methods

To evaluate the performances of the two methods proposed in this paper, we generated the following 12 experimental examples. There are already 3 firms in the competitive market, each with 5 chain facilities. M is the number of customers, N is the number of potential facility locations, and G is the budget of the new entering firm. The attraction are generated randomly within . The lower bound of demand and the ranges of uncertain demand are generated randomly within . The opening costs are generated randomly within . The robust control parameter takes value as . The method is solved by CPLEX 12.5 to get the exact solutions, and the method is coded on the platform of Matlab R2015b. Parameters in the method are . The settings of these parameters are determined based on the results of multiple calculations, and they ensure that the algorithm can find an approximate optimal solution in a short time. Due to the randomness of the algorithm , each example is calculated 20 times. The results include the objective values, CPU times for both methods, and the gaps between objective values and CPU times are presented in Table 4. Here the gaps of two indicators between two methods are defined as follows: Gap (Obj, Time) = [(Obj, Time)(MIP) − (Ave, Time)(SAS)]/(Ave, Time)(SAS).

From Table 4, it is obvious that only for small-scale problems, method is better than method both in objective value and CPU time. When the scale of problems increases, the CPU time required for the method increases dramatically, while the CPU time required for the method remains stable. At the same time, the objective values obtained by the method are only slightly worse than that obtained by the method . Therefore, the method is a good choice to solve large-scale competitive facility location problems.

6. Conclusion

When consumers faced with the choice of competitive chain facilities that offer exclusive services, all current rules cannot describe these customers’ behaviors very well. In this paper, the partially proportional rule is proposed to describe this kind of customer behavior. After expressing the uncertainty of customer demand with interval data, we use robust optimization method to study the competitive facility location problem with uncertain demand. Two methods are proposed to solve the presented robust optimization model. One is an exact method , and it is based on directly solving an integer model equivalent to the original robust model. The other is a heuristic method , which is a combination of the simulated annealing framework and a sorting-based algorithm. A quasi-real example is analyzed in detail. One result that is found through this example is the budget of the new entering firm is an important factor for capturing more market share. Another result is for any given budget, there is a threshold of the robust control parameter . If the parameter is larger than this threshold, the market share captured by the new entering firm will increase as this parameter increases. Random example tests show that for small scale problems, the method is a good choice. But for large-scale problems, the method is better than because the former only requires very little computation time, while maintaining a fairly good calculation result.

For customer behavior of choosing competitive facilities that provide exclusive services, the partially proportional rule is more suitable to describe it than other rules. The proposed model can be used to solve the competitive facility location problem with different customer behavior rules, as long as there is uncertainty in customer demand. We are looking forward to studying the competitive facility problem with more uncertain parameters in the future.

Data Availability

The data supporting this study have been provided.

Conflicts of Interest

The author declares that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

The National Social Science Fund of China (17BGL182) supported this work.

Supplementary Materials

The data set of the 49-node example in Section 5.1 is provided as supplementary materials. The first two columns of the data are latitudes and longitudes of the capitals of the continental United States plus Washington, DC. The third column denotes the existing facilities. The fourth and fifth columns are the lower demand limits and upper demand limits , respectively. The sixth column is the opening cost at the potential facility locations . The last column is the name of the states. (Supplementary Materials)