A Capacitated Location-Inventory Model with Demand Selection

In this paper, we consider an integrated supply chain network design problem, which incorporates inventory and pricing decisions into the capacitated facility locationmodel.We assume that eachwarehouse has a capacity limitation that limits the average demand flowing through the warehouse and that the supplier can choose whether to satisfy each potential retailer’s demand. We formulate the problem as a nonlinear integer programming model and solve the model via a Lagrangian relaxation based approach. We develop an efficient algorithm to solve the subproblem that arises from the Lagrangian relaxation procedure. Finally, we conduct extensive computational experiments to test the performance of the algorithms proposed in this paper and provide the managerial insights based on the computational results.


Introduction
Network design is among the most important supply chain decisions, as their implications are significant and long lasting.Thus, in designing a supply chain distribution network, the supply chain drivers, such as facilities, transportation, inventory, and pricing, should be considered together to support the competitive strategy of a company and maximize the supply chain profits [1].Traditionally, decisions associated with these drivers are made separately.In this paper, we study an integrated supply chain distribution network design problem, which incorporates inventory and pricing decisions into the capacitated facility location model.We consider that a supplier (company) owns the production facility and produces one type of product to serve a set of potential retailers, each facing an uncertain demand.The supplier replenishes the retailers via the warehouses that serve as intermediate facilities for the shipment of the product.The objective of the supplier is to determine how many warehouses to locate and where to locate them, which retailers to assign to each warehouse, when and how many to reorder and what level of the safety stock to maintain at warehouses, and how much the sale price of the product at each region associated with a particular warehouse to maximize the total profits, while ensuring that the capacity constraints for warehouses are not exceeded and maintaining a specific service level.The total profits for the supplier are equal to the total revenues subtracting the transportation cost, the working inventory and safety stock holding cost at the warehouses, and the locating and operating cost of the warehouses.
In the supply chain distribution network design problems, most papers assume that all the demands should be served by the supplier, and then the objective is to minimize the system-wide costs.In practice, it is advisable that the supplier may lose some potential retailers to competitors in order to increase the profits because either the costs of serving those retailers can be very costly [2,3] or the warehouses serving those retailers have capacity limitations, possibly both.Therefore, in this paper, we consider that the supplier can choose whether to satisfy each potential retailer's demand and then explore the benefits that the supplier will achieve by conducting this market strategy.The supplier determines which retailers to serve via the pricing decisions.In the pricing decisions, each retailer holds a reserve price for the product, and if the per unit cost charged for a retailer by the supplier is higher than the retailer's reserve price, the supplier will lose the retailer.The per unit cost charged for a retailer contains the wholesale price of the product and the per unit transportation cost from a specific warehouse to the retailer.Thus, it is important for the supplier to determine the wholesale price of the product at each region associated with a particular warehouse, which directly affects the supplier in terms of the level of profit achieved as well as the demand profile that the supplier attempts to serve.
Furthermore, in this paper, we assume that each warehouse is subject to a capacity limitation, such as in the capacitated fixed charged facility location problem, which limits the average demand that flows through the warehouse.This situation might occur, for example, in contexts where the plant has a finite annual production capacity, or the warehouse has only a fixed number of square feet of storage space.The warehouse capacity constraints considered in this paper are different with those considered in Ozsen et al. [4], which limit the maximum possible inventory accumulation at the warehouses.
In summary, the distribution network design problem studied in this paper has the following features.(i) Each warehouse has a capacity limitation that limits the average demand flowing through the warehouse; (ii) each retailer faces an uncertain demand; (iii) the supplier can choose which potential retailers to serve; (iv) the objective of the supplier is to maximize the total profits.The main contribution of this paper is that, in designing the supply chain distribution network, we consider these features simultaneously and introduce an algorithm to solve the problem effectively.
The rest of this paper is organized as follows.In Section 2, we provide a review of the literature.In Section 3, we present a nonlinear integer programming model for the capacitated distribution network design problem.In Section 4, we introduce a Lagrangian relaxation based approach to solve the capacitated distribution network design problem and an efficient algorithm to solve the subproblem that arises from the Lagrangian relaxation procedure.In Section 5, we present the computational results.Finally, we conclude our discussions in Section 6.
In this paper, we only review literature on the locationinventory models that aim to minimize the sum of facility location costs, transportation costs, and inventory-related costs by choosing the optimal warehouse locations among given candidate sites to serve a set of retailers and at the same time determining the optimal inventory replenishment policies for the warehouses.Nozick and Turnquist [27] study a distribution network design problem, in which the inventory policy adopted at the distribution centers is an (s-1, s) continuous review policy, and the safety stock holding costs are approximated by a linear function of the number of distribution centers.They formulate the problem as a fixed-charge facility location model with inventory costs included in the fixed location cost as a constant.Shen et al. [28] study an uncapacitated location-inventory problem, in which the mean-to-variance ratio of the demand at each retailer is identical for all retailers, and the risk-pooling benefits are considered.They formulate the problem as a setcovering model and then make use of a column generation based approach to solve the set-covering model.Daskin et al. [29] formulate the same problem as a nonlinear integer programming model, and they develop a Lagrangian relaxation procedure to solve their model.Shu et al. [30] extend the works of Shen et al. [28] by proposing an effective algorithm for the pricing problem that arises from the column generation for all demand patterns.Sourirajan et al. [31] study a distribution network design problem that focuses on the trade-off between the safety stocks and the leading times at the DCs, in which they assume the leading time at a DC is directly related to the total demands assigned to the DC.They consider the facility location costs, the pipeline inventory costs between the production facility and the DCs, and the safety stock holding costs at the DCs, but they ignore the transportation costs between the DCs and the retailers in their model.They solve their model by a Lagrangian relaxation based algorithm.Berman et al. [32] study a coordinated location-inventory problem in which they consider a periodic review (, ) inventory policy at distribution centers and two types of coordination for the distribution network system, i.e., partial coordination and full coordination.They show that the full coordination increases the location and inventory costs, but it is likely to reduce the overall costs of running the distribution system.
The models we discuss above deal with the problems in which the facilities are uncapacitated.Miranda and Garrido [33] propose a distribution network design model by incorporating inventory decisions into the capacitated facility location model, which also captures the risk-pooling effects.They develop a heuristic algorithm based on Lagrangian relaxation for solving the model by decomposing the corresponding Lagrangian dual problem into three subproblems.Ozsen et al. [4] study a capacitated logistic network design problem with risk-pooling, in which they consider the capacity constraint on the maximum possible inventory accumulation at a given warehouse instead of the capacity constraint that limits the average quantity of products flowing through the warehouse.They formulate the capacitated logistic network design problem as a nonlinear integer programming model and solve it by a Lagrangian relaxation based algorithm.Ozsen et al. [9] study a similar problem such that they consider multisourcing retailers; i.e., each retailer can be supplied by two or more warehouses.Vidyarthi et al. [34] and Park et al. [35] study the multistage distribution network design problems that need to determine the number and locations of suppliers or plants.Vidyarthi et al. [34] assume that both DCs and suppliers have finite capacity, and Park et al. [35] assume that only DCs have finite capacity.
However, the literature on supply chain distribution network design problems with integrating the pricing decisions is rather limited.Zhang [36] proposes a profit-maximizing location model, where a firm, facing a set of potential retailers, needs to determine where to locate its center warehouse and the price for its product.He assumes that each retailer holds a reserve price for the product and if the actual price charged by the firm is higher than the reserve price, then the retailer will be lost from the firm.But he ignores the inventory cost in his model.Shen [3] studies a location-inventory problem with integrating the pricing decisions.He considers that each warehouse has no capacity limitation and each retailer faces a deterministic demand.He formulates the problem as a set-packing model and then solves it via column generation.Shu et al. [37] propose a twoechelon location-inventory model with a profit-maximizing objective.They assume that both warehouses and retailers carry inventory, and the supplier is in charge of the warehouse-retailer inventory replenishment.They formulate the problem as a set-packing model and then introduce a column generation based approach for solving the setpacking model.Ahmadi-Javid and Hoseinpour [38] study a profit-maximization multicommodity supply chain network design problem.They assume that the customer demands are price-sensitive and formulate the problem as a mixedinteger nonlinear programming model for the uncapacitated and capacitated distribution centers, respectively.They use Lagrangian relaxation based algorithm to solve their model.
In this paper, we present a capacitated location-inventory model by incorporating inventory and pricing decisions into the capacitated facility location model.The problem we studied in this paper is closely related to Miranda and Garrido [33], Shen [3], Ozsen et al. [4], Shu et al. [37], and Ahmadi-Javid and Hoseinpour [38].Comparing with Shen [3] and Shu et al. [37], the feature of our paper is that we consider each warehouse has a capacity limitation, and each retailer faces an uncertain demand; comparing with Miranda and Garrido [33], the feature of our paper is that we integrate the pricing decisions into the supply chain network design; comparing with Ozsen et al. [4], the feature of our paper is that we consider different warehouse capacity constraints; comparing with Ahmadi-Javid and Hoseinpour [38], the feature of our paper is that we consider that the supplier can choose which retailers to serve, and each retailer faces an uncertain demand.We formulate this capacitated distribution network design problem as a nonlinear integer programming model and then introduce a Lagrangian relaxation based approach to solve this model.The subproblem that arises from the Lagrangian relaxation procedure is a nonlinear zero-one integer programming model, which maximizes a convex objective function subject to a single knapsack constraint.We also use the Lagrangian relaxation to solve the subproblem.Finally, we conduct the computational experiments to test the efficiency and effectiveness of the algorithms proposed in this paper.

Model Formulation
We consider a distribution network consisting of a production facility, a set of candidate warehouses, and a set of potential retailers, where the location of the production facility and the retailers is known, and the production rate of production facility is infinite.We assume that a supplier owns the production facility and produces one type of product, where we ignore the production cost.After production, the supplier transports the product to several warehouses for distribution.We also assume that the leading time from the production facility to each candidate warehouse is deterministic, and the demand at each retailer is uncorrelated across retailers, and that if the supplier chooses to serve a retailer, the supplier should satisfy all the demand of the retailer.In this paper, we consider single source strategy, which means a retailer is only served by one warehouse.The notations we used throughout the paper are defined as follows.
Sets : set of potential retailers; : set of candidate warehouse locations;

Inputs and Parameters
: mean demand per year at retailer , for each  ∈ .We calculate   using the formula   =    , where    is the mean demand per day at retailer  and  is a constant used to convert daily demand into annual demand, e.g.,  = 300;   : variance of daily demand at retailer , for each  ∈ ;  0 : per unit transportation cost of delivering product from the supplier to warehouse , for each  ∈ ;   : per unit transportation cost of delivering product from warehouse  to retailer , for each  ∈ ;   : fixed cost of placing an order from warehouse  to the supplier, for each  ∈ ; ℎ  : per unit inventory holding cost rate at warehouse  per year, for each  ∈ ;   : leading time in days for deliveries from the supplier to warehouse , for each  ∈ ; : the desired percentage of retailers orders satisfied;   : standard normal deviate such that ( ≤   ) = ; V  : the reserve price of retailer , for each  ∈ ;   (⋅): the (annual) cost of locating and operating warehouse , which is a concave and nondecreasing function of the annual demand flows through warehouse , for each  ∈ ;   : the capacity of warehouse , for each  ∈ ;

Decision Variables
: per unit wholesale price of the product at the region served by warehouse , for each  ∈ ;   = 1, if we locate a warehouse at candidate location , and 0 otherwise, for each  ∈ ;   = 1, if the demand at retailer  is served by warehouse , and 0 otherwise, for each  ∈ ,  ∈ .
Let   (  ,   ) denote the per unit gross profits of serving retailer  using warehouse  for the supplier, and we define   (  ,   ) in the following.We know that when the wholesale price plus the per unit transportation cost charged to a retailer is higher than the retailer's reserve price, this retailer will turn to some other suppliers; i.e., the supplier will lose this retailer.Then we have [3,36] We know that when   +   ≤ V  , retailer  can be served by the supplier for the per unit cost retailer  has to pay lower than his reserve price, and thus the supplier can get revenues from serving retailer ; when   +   > V  , retailer  cannot be served by the supplier for the per unit cost, retailer  has to pay higher than his reserve price, and thus the supplier cannot get any revenue from retailer .Therefore, the pricing decisions determine which retailers can be served by the supplier, which is important for maximizing the supplier's profits.
In this paper, we focus on Type I service level [39], i.e., probability of no stock out per order cycle, which is one of the most commonly used service levels in inventory management [40,41], and for calculating the safety stocks at warehouses, we assume that the normal approximation for demands is appropriate [28].Since retailer demands are independent and the normal approximation for demands is appropriate, the warehouse demands can be represented by the Normal distribution.Thus, for each warehouse , the leading time demand at warehouse  is Normal distribution with mean   √∑ ∈      and variance   ∑ ∈  2    , and then the safety stock required at warehouse  with the probability  of no stock out per order cycle is   √  ∑ ∈  2    .Using the above notations, we formulate the capacitated distribution network design problem as follows: ∈ {0, 1} , ∀ ∈ .
The first term of the objective function ( 2) is the gross profits of serving retailers using warehouses; the second term describes the working inventory holding cost determined by the EOQ model with ordering cost   , per unit holding cost ℎ  , and demand ∑ ∈     ; the third term represents the safety stocks holding cost; the fourth term is the cost of locating and operating warehouses, which is associated with how many demands flow through the warehouses.Constraint (3) states that a retailer may not be assigned to any warehouse.Constraint (4) stipulates that retailers can be only assigned to opened warehouses.Constraint (5) stipulates that the total demands served by a warehouse cannot exceed the capacity limit of the warehouse.Lastly, constraints (6) and ( 7) are standard integrality constraints, and constraint (6) also requires that each retailer is assigned to exactly one warehouse.
Note that the above model formulated is an extension of the classic capacitated facility location problem, which is already NP-hard, and furthermore the objective function (2) of the above model contains nonlinear cost terms as well.Thus any commercial solver cannot solve the above model effectively.In what follows, we develop a Lagrangian relaxation based approach to solve the above nonlinear integer programming model.

Solution Approach
Lagrangian relaxation method has been used successfully to solve many integer programming problems, such as facility location problems and supply chain network design problems [4,9,29,34,42,43].Fisher [44] gives an excellent discussion of the Lagrangian relaxation method for solving the integer programming problems.In this section, we develop a Lagrangian relaxation procedure to solve P. We first obtain the Lagrangian dual problem of P by relaxing the assignment constraint (3), and then, for fixed Lagrangian multipliers, the Lagrangian dual problem is a nonlinear zero-one integer programming model, which maximizes a convex objective function subject to a single knapsack constraint.We also use the Lagrangian relaxation method to solve the Lagrangian dual problem.We describe the details of the Lagrangian relaxation procedure as follows.
Note that we solve the Lagrangian dual problem for fixed values of Lagrange multipliers,   , to obtain an upper bound on optimal objective function value of the capacitated distribution network design problem.Clearly, we want to maximize the Lagrangian function ( 8) over decision variables,   ,   , and   , and to minimize the function over the Lagrange multipliers,   .Minimizing the Lagrangian function (8) over the Lagrange multipliers,   , can be done using subgradient optimization [44].Therefore, we focus on solving the problem for fixed values of the Lagrange multipliers,   .
For fixed values of Lagrange multipliers,   , we solve the following subproblem for each candidate warehouse .
We note that P  is a 0-1 integer programming problem with a nonlinear convex objective function and a single knapsack constraint.Let   be the optimal objective function value of P  , and we easily know that ∑ :  >0,∈   +∑ ∈   is an upper bound for the optimal objective function value of P. Since   is a decision variable in the objective function of P  , we first show how to solve P  for a fixed   , and then we show how to solve P  when   is a decision variable as well.
. . .Solution Approach for P  with a Fixed   .For a fixed   , we solve P  by the Lagrangian relaxation method as well.By relaxing the capacity constraint (5), the Lagrangian dual problem of P  is then formulated as follows: min where   is the Lagrange multiplier associated with the capacity constraint (5), and   ≥ 0.
For a fixed value of the Lagrange multiplier   , we want to maximize (11) over the assignment variables   .Thus, for a fixed   , the Lagrangian relaxation dual problem of P  reduces to the following problem: Let   be the optimal objective function value of L   , and then   +     is an upper bound for the optimal objective function value of P  .Furthermore, we can rewrite L   as a different but equivalent formulation as follows: where  is a subset of .Let   denote the optimal solution of ( 15), and we have the following observation.
Proof.We prove it by contradiction.Suppose that there exists  ∈   such that   (  ,   ) ≥ 0. Then we let Ŝ =   \ {}, and note that the objective function value of (15) for Ŝ is less than that for   , which contradicts the optimality of   .We therefore establish that if  ∈   , then   (  ,   ) < 0.
For a fixed   , ( 15) can be solved in ( 2 log ) time [23,30], where  is the number of retailers.In this paper, we solve (15) by a dual algorithm [23], and we depict the steps of the dual algorithm as follows.
Step .For each pair of (,   ),  ∈  −  ,   ∈  −  , we solve the following linear equations: Let (   ,    ) denote the solution of ( 16), and we choose the solution if both    and    are positive.
Step .For each  = 0, 1, . . ., , Step .From the feasible solutions of (15), we select the one with minimum objective function value of ( 15), and then this solution is the optimal solution of (15).
Having solved the Lagrangian relaxation problem L   , we need to minimize the Lagrangian function (11) over the Lagrange multiplier   .We do so using a subgradient optimization procedure.
We then show how to obtain a lower bound for the optimal objective function value of P  , and it is sufficient to find a feasible solution for P  .At each iteration of the Lagrangian procedure, we depict the details of the step for obtaining a feasible solution for P  as follows: Let Step .If ∑ ∈ 1      ≤   , then the current solution for the Lagrangian relaxation problem is feasible for P  , and for a retailer  ∈  0  , we also let   = 1 for which the capacity constraint is satisfied and that increases the objective function value to the maximum based on the assignments made so far.
Step .If ∑ ∈ 1      ≥   , for each  ∈  1  , we first let   = 0 for which   = 1 and that decreases the objective function value to the minimum based on the assignments made so far until the capacity constraint is satisfied.Secondly, we update the set  0  , and then for a retailer  ∈  0  , we also let   = 1 for which the capacity constraint is satisfied and that increases the objective function value to the maximum based on the assignments made so far.
. . .Solution Approach for P  When   Is a Decision Variable.In this subsection, we discuss how to solve P  as well as decide the optimal sale price   when the price   charged to retailers is a decision variable.We let  *  denote the optimal solution of P  .
Let   = { |   < V  ,  ∈ }.Note that   is a set of potential retailers that can be served by warehouse  when the per unit sale price of the product is zero.To solve P  and decide the optimal sale price   , we first give the following theorem.

Theorem 2. Let 𝑝 *
be the optimal price of P  and    = V  −   ,  ∈   .When the optimal solution  *  of P  is a nonempty set, the optimal price  *  is in the set {   ,  ∈   }.
Proof.For the notational convenience, we relabel the retailers in set   as {1, 2, . . ., |  |} so that Then we have   ( *  ,   )  −   ≤ 0, ∀ ∈ , and it is obvious that the optimal solution  *  of P  is an empty set.This is to say, no retailers can be served by warehouse  because of the high price, and at this time the optimal price has no practical significance.Thus we have 0 <  *  ≤ This would imply that  *  is not the optimal price of P  and, therefore, a contradiction exists.We establish that  *  ∈ {   ,  ∈   }.
Theorem 2 shows that we can attain the optimal price  *  for P  within the set {   ,  ∈   }.Based on Theorem 2, we describe how to solve P  when   is a decision variable as follows: Step .For each  ∈   , we first solve P  at   =    by the algorithm proposed in Section 4.1.1.
Step .Then, we select the one with maximum objective value, and the corresponding solution and price are the optimal solution and price of P  .
. .Obtaining a Lower Bound.In this subsection, we show how to obtain a feasible solution of P, and, therefore, we obtain a lower bound on optimal objective function value for P. We depict the details of obtaining a feasible solution of P as follows.At each iteration of the Lagrangian procedure, we first fix the warehouse locations at those locations for which   = 1 in the current Lagrangian solution.Then we assign the retailers to the fixed warehouses in a two-step heuristic process.
Step .For retailer , if ∑ ∈   ≥ 1, then we assign retailer  to warehouse  for which   = 1 and that increases the profit to the maximum based on the assignments made so far.
Step .For retailer , if ∑ ∈   = 0, then we assign retailer  to warehouse  for which the capacity constraint is satisfied and that increases the total profit to the maximum based on the assignments made so far.

Computational Results
In this section, we summarize our computational experience with the algorithms outlined in the above section.We conduct two sets of experiments as follows.The first set of experiments is designed to solve P  that arises from the Lagrangian relaxation procedure, and the second set of experiments is designed to solve P. We solve P  via the algorithm proposed in Section 4.1 and solve P by the Lagrangian relaxation procedure.All the algorithms are coded in C++.We conduct all the experiments on a PC with a dual-core CPU of 2.27 GHz and 4G RAM running the Windows 7 operating system.We randomly generate the potential retailer locations and the candidate warehouse locations on the plane [0, 1000] ×[0, 1000], and we denote the point (0, 0) as the supplier or plant.Without loss of generality, we let  = 300.
All the parameters for the problem instances are generated as follows.For each retailer , we randomly generate the daily mean demand    in [50,100], the variance of daily demand   in [10,50], and the reserve price V  in [10,50].The fixed ordering cost   is set to 1000, the per unit holding cost rate ℎ  is set to 5, the lead time   is set to 10, and the standard normal deviate   is set to 1.96 corresponding to a 97.5% service level.Table 1 summarizes the values for the parameters.
Let  be the total potential retailer demands, i.e.,  = ∑ ∈   .In the first experiment, we fix a potential warehouse, say , and solve the corresponding subproblem P  .We denote the capacity of the fixed warehouse  as 0.5, 0.7, and 0.9, respectively, which allows for the test problem instances with small capacity, medium capacity, and large capacity relative to the total potential retailer demands.In the second experiment, for each candidate warehouse , we randomly generate the capacity limit   in [0.3, 0.5], [0.5, 0.7], and [0.7, ], respectively, which allows for the test problem instances with small capacity, medium capacity, and large capacity relative to the total potential retailer demands.We let the locating and operating cost for warehouse  be 1000 × √  , where   is the total demands that flow through warehouse ,  ∈ .
Let  0 be the Euclidean distance between the supplier and warehouse  and   be the Euclidean distance between warehouse  and retailer .Then, we define the per unit transportation cost from the supplier to warehouse ,  0 , and the per unit transportation cost from warehouse  to retailer ,   , in the following:  number of warehouses and the number of potential retailers in the test problem instances, respectively; the columns titled "#WHs Opened", "#Ret Served", "CPU(s)", and "Gap(%)" represent the number of warehouses opened, the number of the retailers served by the opened warehouses, the number of CPU seconds taken by the Lagrangian procedure, and the percentage gap between upper bound and lower bound solutions, where Gap (%) = 100 × (Upper Bound -Lower Bound)/Lower Bound; the column titled "Δ(%)" in Tables 5 and 6 represents the percentage gap between the flexible model (with the constraints ∑ ∈   ≤ 1, ∀ ∈ ) and the corresponding serve-all model (with the constraints ∑ ∈   = 1, ∀ ∈ ), where Δ(%) = 100 ×( 1 −  2 )/ 2 ,  1 is the best objective value of the flexible model we obtained, and  2 is the best objective value of the corresponding serveall model we obtained.We solve the serve-all model via the Lagrangian relaxation algorithm as well.Note that the serveall model means that all the potential retailers must be served by the opened warehouses, which requires that ∑ ∈   ≥ , i.e., the problem instances with medium and large capacity constraints.
In Tables 4-6, we observe that the average gap is 3.64% and the largest gap is no more than 8%, which shows that the algorithm proposed in this paper for solving P is practical and efficient.As shown in Tables 5 and 6, the percentage of extra profits gained from the flexible model is 10% on average, and the largest percentage of extra profits gained from the flexible model is 14.7%, which shows that the flexible model is more profitable than the serve-all model.From Tables 5 and  6, we also obtain that the percentage of extra profits gained from the flexible model decreases with increasing of the capacities of the warehouses and increases with increasing of the number of potential retailers.These results show that the flexible model proposed in this paper can gain more profits when the problem instances are large-scale and when the capacities of the warehouses are small relative to the retailer's demands.Based on the above analyses, we conclude that the demand selection strategy is more profitable for the supplier in practice.

Conclusions
In this paper, we study a supply chain distribution network design problem in which each warehouse has a capacity limitation that limits the average demand flowing through the warehouse and each retailer faces an uncertain demand.

Journal of Advanced Transportation
We consider that the supplier can choose whether to satisfy each potential retailer's demand for the warehouse capacity limitation, or the high cost of serving the retailer, and the supplier determines which retailers to serve via the pricing decisions.The one contribution of this paper to the literature is that in designing the supply chain distribution network, we consider warehouses have capacity limits, retailers face uncertain demands, and the supplier has flexibility in determining which retailers to serve simultaneously.We formulate this problem as a nonlinear integer programming model and solve the nonlinear integer programming model via a Lagrangian relaxation based approach.The subproblem that arises from the Lagrangian relaxation procedure is a 0-1 integer programming problem, which maximizes a nonlinear convex objective function subject to a single knapsack constraint.The other contribution of this paper to the literature is that we introduce an efficient algorithm to solve the subproblem.Finally, we conduct two sets of computational experiments to test the efficiency of the algorithms proposed in this paper.From the computational experiments, we obtain the following results: (1) the algorithm proposed in this paper for solving the subproblem P  is efficient, and, in most cases, the algorithm can solve the P  to optimality; (2) the Lagrangian relaxation based algorithm proposed in this paper for solving the master problem P is efficient, where the average gap is 3.64%, and the largest gap is no more than 8%; (3) the flexible model is more profitable than the serve-all model, which means the demand selection strategy is more profitable for the supplier; (4) when the problem instances are largescale, and when the capacities of the warehouses are tight, the demand selection strategy can make more profits for the supplier.
The problem studied in this paper can be extended in many ways.First, in this paper, we consider the capacity of each warehouse is exogenous; i.e., the capacity of each warehouse is regarded as a fixed input, and, in practice, we want to determine the best capacity level to install at each warehouse.Thus, an extension of the problem studied in this paper is to allow the capacity of each warehouse to be a decision variable [45,46].Second, we only deal with medium size problem instances, and, in the future, we want to explore how to design the efficient approximation algorithms for the problem [6,[47][48][49][50].

Table 4 :
The computational results for P with small capacity.

Table 5 :
The computational results for P with medium capacity.

Table 6 :
The computational results for P with large capacity.