A Particle Swarm Optimization Algorithm for Solving Pricing and Lead Time Quotation in a Dual-Channel Supply Chain with Multiple Customer Classes

The combination of traditional retail channel with direct channel adds a new dimension of competition to manufacturers’ distribution system. In this paper, we consider a make-to-order manufacturer with two channels of sale, sale through retailers and online direct sale. The customers are classiﬁed into diﬀerent classes, based on their sensitivity to price and due date. The orders of traditional retail channel customers are fulﬁlled in the same period of ordering. However, price and due date are quoted to the online customers based on the available capacity as well as the other orders in the pipeline. We develop two diﬀerent structures of the supply chain: centralized and decentralized dual-channel supply chain which are formulated as bilevel binary nonlinear models. The Particle Swarm Optimization algorithm is also developed to obtain a satisfactory near-optimal solution and compared to a genetic algorithm. Through various numerical analyses, we investigate the eﬀects of the customers’ preference of a direct channel on the model’s variables.


Introduction
e rapidly expanding Internet provides an opportunity for organizations to distribute their products via both direct channel and traditional retail channel. Some personal computer manufacturers (like Dell company), apparel retailers, and automotive industries (like General Motors) are examples of companies that use hybrid of both direct and retailer channels.
In the direct channel, the firm interacts with consumers directly through Internet. ere are a number of benefits from direct channel distribution such as controlling the distribution and pricing directly, providing a broader product selection, and improving firms' visibility [1].
Despite hybrid channel's benefits such as capturing a larger share of the market, combining the retail distribution channel with direct channel may pose some challenges, including pricing policies, distribution strategies, and conflicting demands placed on internal company resources such as capital, personal, products, and technology by multiple channels [2].
To overcome these challenges, supply chain's members negotiate the retail as well as wholesale price to cooperate with each other. In addition to the product price, there are other factors such as product availability and service quality that contribute to consumer preference of the direct channel [3].
Delivery lead time is one of the important factors that can affect customers' demand [3]. is is the reason many e-retailers, such as Amazon.com, BestBuy.com, Walmart. com, and FYE.com, try to offer competitively quoted lead times [4].
is paper is focused on competitive pricing strategy as well lead time decisions in a supply chain for a manufacturer that sells the products through two channels. One channel is the traditional one in which the firm uses an intermediary to reach final consumers, while the other is a direct channel in which the customer places direct orders through the Internet. It is assumed that there are multiple classes of customers in the market. e demand of each class depends on the price and lead time. A finite planning horizon is considered, where the production capacity in each period is finite but varies. e manufacturer has to respond to the customers quickly, based on the available capacity and the orders in pipeline. Due to the finite production capacity and the demand sensitivity to price and lead time changes, the manufacturer needs to decide on the selling price and the lead time of direct customers, the contract terms with the retailer, and the wholesale price as well as the production schedule to maximize his own profit.
We consider two different cases, centralized and decentralized dual-channel supply chain. We propose a model to determine suitable lead time and price, simultaneously for each case. In a decentralized dual-channel supply chain, the two members interact within the framework of Stackelberg game. Under this framework, the manufacturer, as the leader, determines the wholesale price for retailer and also a price and lead time for direct sale in each period. e intermediary then reacts by choosing a retail price to maximize its own profit. On the other hand, for a centralized dual-channel supply chain, the direct sale price, the traditional retail price, and the quoted lead time in the direct channel are determined by a vertically integrated manufacturer.
e resulting models happen to be binary nonlinear programs. ese kinds of problems with large dimensions are not usually solved by exact algorithms. erefore, an alternative method for solving the models is to use metaheuristic algorithms to obtain a near-optimal solution with reasonable computational time.
us, an algorithm of Particle Swarm Optimization (PSO) is developed to solve the models. e results will be compared, in terms of solution quality and computational time, with genetic algorithm (GA). e remainder of this paper is structured as follows. Section 2 reviews the more relevant literature. We formulate the models in Section 3. e principles of the PSO metaheuristic and our algorithm are explained in Section 4. Section 5 presents some numerical examples of our model to examine the effects of customer preference of a direct channel on the pricing strategies and lead time decisions. Our conclusions are summarized in Section 6.

Literature Review
is paper focuses on competition in a dual-channel supply chain. A comprehensive review of multichannel models can be found in [1,5].
Several researchers and practitioners have focused on dual-channel supply chains during the last decade. Different aspects of this chain are investigated. Most of the papers that formulated the dual-channel supply chain focused on the competition context in the supply chain and pricing optimization issues in each channel [6][7][8][9]. e researchers combined the pricing issues with other aspects of supply chain. In this regard, we can refer to sales effort determining and service management in each channel [10][11][12][13], contract optimization [6,14], disruption management [15][16][17], product variety in supply chain [18,19], and multiperiod model [20].
Xu et al. [32] considered the unit delivery cost (m/t) with m > 0 if the product is to be delivered with lead time t. Batarfi et al. [33] considered the combination of two production approaches: make-to-order and make-to-stock in direct online channel and indirect offline channel, respectively. e demand depends on the prices, the quoted delivery lead time, and the product differentiation.
Yang et al. [4] modeled delivery lead time optimization for perishable products in a dual-channel supply chain. Modak and Kelle [34] considered that the demand not only is dependent on price and delivery time but also is stochastic. e lead times in the above papers are determined exogenously (determined by the sales department, without knowing the actual production schedule). In addition, the production capacity is unlimited. Table 1 illustrates the major literature review with our paper included.
Based on the above literature, this paper investigates the joint decision on production, pricing, and lead time in a dual-channel distribution system where the lead times are assigned internally by the scheduling model. Aside from considering the above literature, we address a new model in which the production capacity in each period is limited. We also consider multiple customer classes that differ in their arrival (commitment) times, quantities demanded, and sensitivity to price and lead time.

Problem Description
We consider a dual-channel supply chain in which a manufacturer sells to retailers as well as directly to end customers. e planning horizon is assumed to be finite and divided into periods of equal length. e capacity of the manufacturer is limited but may vary in different periods. Customers can be classified with respect to their sensitivity to price and lead time. Furthermore, the attributes of customers such as arrival (commitment) times, quantities demanded, unit production, and holding costs are different.
In each period, the manufacturer must set a wholesale price for the traditional retail channel as well as setting both 2 Advances in Operations Research price and lead time (or due date) for the customers of direct channel. For the customers of traditional retail channel, the production must be scheduled at the same period, at which the order arrives. However, the manufacturer can quote a lead time (due date) to the customers of direct channel. e production for these customers is scheduled within any period between the arrival time of order and the quoted due date. A holding cost occurs for any order of direct channel which is completed before the quoted due date. As mentioned before, in order to evaluate the effects of the delivery lead time and customer's preference of the direct channel on the pricing decisions of the manufacturer and retailer, we consider two different dual-channel supply chains, centralized and decentralized systems.
3.1. Notation. We use the following notation. i,t : production cost of customer of class i ∈ Ψ in the direct channel in period t ∈ T Cp 2 i,t : production cost of customer of class i ∈ Ψ in the retail channel in period t ∈ T Cr i : operational cost in the retail channel for customer of class i ∈ Ψ Ce i : operational cost in the direct channel for customer of class i ∈ Ψ K t : production capacity available in period t ∈ T D r i,j : demand (in terms of production capacity units) for customer of class i ∈ Ψ in the retail channel, quoted due date j ∈ [e(i), . . . , T] D s i,j : demand (in terms of production capacity units) for customer of class i ∈ Ψ in the direct channel, quoted due date j ∈ [e(i), . . . , T]

Decision Variables
Z i,j : 1 if the due date j, j ∈ [e(i), . . . , T], is selected (quoted) for customer of class i ∈ Ψ in direct channel; 0 otherwise P r i : price in the retail channel for customer of class i ∈ Ψ P s i : price in the direct channel for customer of class i ∈ Ψ W i : wholesale price charged for customer of class i ∈ Ψ in the retail channel Y i,t : total production (in units of capacity) for customer of class i ∈ Ψ in the direct channel in period t ∈ T H i : total inventory of customer of class i ∈ Ψ

Demand Function.
In line with Kurata et al. [35], Cai et al. [6], and Hua et al. [3], we assume that the demand is a linear function of the self-and cross-price and lead time, as follows: where D s i and D r i denote the consumer demand to the manufacturer and the consumer demand to the retailer, respectively, a i denotes the base level of industry demand or the demand rate for customer of class i, and θ i (0 < θ i < 1) represents the initial portion of customers of class i that prefer the direct channel if the price and lead time are zero. us, 1 − θ i represents the portion of customers of class i that prefer purchasing from the retailer. e coefficients b s i and b r i are the coefficients of the price elasticity in the direct channel and retail channel demand functions, respectively. e cross-price sensitivities α r i and α s i reflect the substitution's degree of the two channels, and c s i is the lead time sensitivity of the demand in the direct channel. If the lead time L s i increases by one unit, c r i units of demand will transfer to the retail channel and c s i units of direct channel's demand will be decreased. e total demand of the two channels should have a negative slope with respect to the retailer's price, the direct sale price, and the quoted lead time. us, we have α r i < b r i , α s i < b s i , and c s i > c r i . Following Hua et al. [3], it is assumed that the manufacturer uses dual channels to sell their goods, and the base level of industry demand or demand rate in both channels is very large; i.e., θ i should not be unreasonably small or large.

Model of the Decentralized Dual-Channel Supply Chain.
In this section, we study a decentralized dual-channel supply chain. In this model both the manufacturer and the retailer make their own decisions separately to maximize their profits. e manufacturer, as the Stackelberg leader, first determines the wholesale price W i , the direct sale price P s i , and the direct sale quoted lead time L s i . en, the retailer chooses his own optimal retail price P r i based on the manufacturer's decisions.

Manufacturer's Best Response.
e goal of the manufacturer as a leader is to maximize his profit, considering the capacity constraints and demand constraints in the systems and constraints determined by the retailer optimization problem. e problem is formulated within the framework of bilevel programming (BLP), first level (the manufacturer model called the leader) and second level (the retailer model called a follower). In BLP model, each decision maker tries to optimize its own objective function without considering the objective of the other party. However, the decision of each party affects the objective value of the other one as well as the decision space.
T j�e(i)
e first and second terms in the first-level objective represent the total revenue and production cost of direct sales. e third term represents the total profit obtained by sales to the retailer and the fourth one is the carrying cost.
Constraint (4) indicates that the wholesale price cannot be higher than the direct channel price; otherwise, the retailer may purchase from the direct channel at a lower price. Constraints (5) ensure that only one due date is chosen for each direct channel customer order. Constraint (6) ensures that if due date j is selected for customer of class i in the direct channel, then exact D s i,j units must be produced and delivered, where D s i,j depends on the selected price and due date. Constraint (7) is a capacity constraint that ensures that the production capacity in each period is not exceeded. An order of a retail channel customer is produced in the customer's arrival time period and delivered instantaneously in the same period, whereas an order of a direct channel customer can be produced in any period between its arrival period and the quoted due date. e first term in the left hand side of constraint (7) indicates the total production for direct channel customers and the second term is the total production for retail channel customers. Constraint (8) indicates that an order of a direct channel customer can be produced in any period between its arrival period and the quoted due date. e required inventory for orders is scheduled in any period prior to its commitment, and the negotiated due date is controlled by constraint (9), where M is a sufficiently large number. Constraint (10) is demand functions for customer orders in the direct channel. e term j − e(i) + 1 is a time interval between the arrival time of an order and the quoted due date, which is called the lead time.
e objective function of second-level optimization problem for the retailer channel is represented by (11). Constraints (12) and (13) control feasibility as well as demand restrictions.
Without considering constraints (12)- (14) for retailer, its best response to the wholesale price can be defined in Proposition 1.

Proposition 1.
e retailer's best response to the wholesale price W i , the direct sale price P s i , and the direct quoted lead time L s i set by the manufacturer is as follows: e Proof of Proposition 1 as well as the other propositions is given in the appendix.
e BLP model (4)- (14) can be formulated as a single level mixed binary problem. is is achieved by replacing the lower level problem (12)- (14) with its Kuhn-Tucker conditions, which we name Model I, as follows.

Advances in Operations Research
Constraints (24)- (27) indicate the optimality condition for the retailer to convert the bilevel problem into a single level one. Proposition 2. For any given due date j, the wholesale price W i , and the direct sale price P s i for customer class i, the optimal retail price P r i is always positive. As a result of Proposition 2, it can be concluded that ] i is always equal to zero, because ] i · P r i � 0 in Model I.

Proposition 3.
Consider an unlimited production capacity in each period for a decentralized dual-channel supply chain.
i is always greater than zero. us, in this case, the manufacturer profit Π s is always strictly jointly concave in W i and P s i . us, for the sake of analytical tractability, we will assume that the cross-price effects are symmetric and α s i � α r i .

Proposition 4.
Let the production capacity be unlimited in each period for a decentralized dual-channel supply chain and suppose there is Advances in Operations Research en, the optimal lead time L s i ∈ [L l i , L u i ] can be found from L l i , L u i , and L 0 i by comparing their Π s values; the one with the largest Π s value is the optimal lead time.

Model of the Centralized Dual-Channel Supply Chain.
In a centralized dual-channel supply chain, a vertically integrated manufacturer controls the wholesale price W i , the direct sale price P s i , the direct sale quoted lead time L s i , and the retailer sale price P r i . e model, called Model II, is as follows.

Model II.
Max e objective function (30) maximizes the total profit of the centralized supply chain, which is the sum of manufacturer's profit and retailer's profit. e remaining constraints are the same as those in Model I.

Proposition 5.
Consider an unlimited production capacity in each period for a centralized dual-channel supply chain. If can be found from L l i and L u i by comparing their Π c values; the largest Π c value is the optimal lead time.

Solution Approach
Models I and II are 0-1 nonlinear programs with continuous price and production decision variables and binary variables for due date selection. It is well known that the problem of nonlinear 0-1 programming is NP-hard [36]. us, these problems cannot be solved easily, and exact algorithms can be time consuming when the number of 0-1 variables is large. An alternative for solving these models is to use approximate algorithms or heuristic methods with reasonable computational times.
In this paper, PSO metaheuristic is used and compared to the GA for modeling and optimization of the pricing and due date setting problems in a dual-channel supply chain.
PSO and GA can handle the whole MINLP easily and naturally, and it is easy to apply it to various problems for comparison with conventional methods.

Advances in Operations Research
If we fix the set of binary variables Z in the models, then both models are converted to quadratic programs that can be solved by the nonlinear programming solver LINGO.

Particle Swarm Optimization (PSO).
e PSO algorithm was first proposed by Kennedy and Eberhart [37]. is method is based on the information shared among members of a species and then used for evolution. Update the velocity and position for every particle.

Advances in Operations Research
Simple structure, ease of implementation, speed of acquiring solutions, and robustness are the advantages of the PSO that persuade us to use it in solving the models.
Recently, PSO algorithms were successfully applied to a wide range of applications. A comprehensive survey of PSO algorithms and applications can be found in the paper by Kennedy et al. [38].
In the proposed PSO algorithm, particle k is represented as L k � L k,1 , L k,2 , . . . , L k,N , which denotes the due dates quoted to N customers. For each L k,i in N-dimensional space, the Z iL k,i is set to one in Models I and II, and the other variables Z i,j , j � e(i), . . . T, and j ≠ L k,i , are fixed to zero for customer order i ∈ Ψ.
is section will present the application of the PSO algorithm for Models I and II. e PSO algorithm for solving these models is illustrated in Figure 1.
e key point in a constrained optimization process is to deal with the constraints. We must modify the original PSO method for constrained optimization. Many methods were proposed for handling constraints, such as methods based on preserving the feasibility of solutions, methods based on penalty functions, methods that make a clear distinction between feasible and infeasible solutions, and other hybrid methods [39]. e most straightforward method is the one based on preserving the feasibility of solutions. In this method, each particle searches the whole space but only keeps track of the feasible solutions. To accelerate the process, we derive an upper bound on the search range by finding the largest lead times such that both demand in retail and direct channel will be nonnegative; i.e., D s i,j ≥ 0 D r i,j ≥ 0. e upper bounds on lead times are given in the following proposition.

Proposition 6.
e best due date for each customer i ∈ Ψ is restricted to a range L l i ≤ L i ≤ L u i . e lower bound of the range is the arrival time of the customer, i.e., e(i), and the upper bound of the range is

Genetic Algorithm (GA)
. GA is a well-known metaheuristic optimization technique based on Darwin's theory of the "survival of the fittest," proposed by Holland [40]. e GA starts with a group of individuals created randomly. e individuals in the population are then evaluated using a fitness value. Two individuals are then selected based on their fitness. ese individuals "reproduce" to create one or more offspring by crossover and mutation operators. e one-point crossover and swap mutation are used in this article. According to the description provided in PSO algorithm, each solution can be represented by an integer string with length N (number of customers). Each gene represents the due date quoted to each customer. e GA algorithm for solving these models is illustrated in Figure 2:

Numerical Studies
In this section, we investigate the relation between the wholesale price, the direct sale price, and the retail price in the decentralized and centralized supply chain. e pricing strategies and quoted lead time decisions are compared for the two models using numerical experiments. e parameters of the proposed PSO algorithm are selected based on parameter tuning. With testing different values for PSO' parameters, the ones selected are those that gave the best results. e best values of the computational experiments are as follows: (1) the values of a population of 30 individuals are used for both GA and PSO, (2) the initial inertia weight is set to 0.9, and (3) the values of the acceleration constants c1 and c2 are fixed to 0.9. e maximum velocity is set as the difference between the upper and the lower bounds, which ensures that the particles are able to fly across the problem-specific constraints region. e other parameters used in GA are crossover rate of 0.80 and mutation rate of 0.3.
Each instance was run for 50 iterations, and 30 replications were conducted for each instance. e performance of PSO algorithm was compared with MINLP solver GAMS and GA algorithm by testing on a set of 8 small size instances in the decentralized (Model I) and centralized (Model II) supply chain.

Advances in Operations Research
GAMS solutions and the presented PSO algorithm. e sum of manufacturer's and retailer's profit is shown in Table 2. All the results show that PSO algorithm can achieve better results than those obtained by GAMS solver.
Comparing the results of the GA and PSO methods, we can see that there is no conclusive winner. In some cases, the GA method results in a better solution and in some cases the PSO method does. Due to the simplicity of the PSO structure, we use it for our sensitivity analysis.
To analyze the effects of models' parameters on prices and profits, we considered 18 different problem groups to carry out this numerical study. For each group of problems, the production capacity of each period (K) is categorized as high and low capacity. e ratio of operational costs in direct channel and retailer (Ce/Cr) is categorized as high, medium, and low. We can have single class customer problems where the price and lead time sensitivities are the same across the customers and multiclass customer problems based on the variability of these sensitivities. With this assumption, the additional categories are based on price sensitivity variability (PV) and lead time sensitivity variability (LV). Table 3 presents the characteristics of different groups of problems. e other parameters are considered as follows:       Figure 5: Comparison of the retail price, direct sale price, and wholesale price in decentralized dual-channel supply chain for different instances.
e holding costs per unit of time for each customer of class i are Ch i � 5, and the production cost in each period t for each customer of class i is Cp i t � 10. e demand rate for customer of class i, a i , is generated randomly from a i ∈ [500, 3000]. We considered a planning horizon with 12 time periods and 6 customer classes arriving during the planning horizon for all problems. Two types of arrival time distributions are considered: the arrival times near the beginning of planning horizon (B) (e(i) ∼ uniform [1,4]) and the uniform distribution along planning horizon (U) (e(i) ∼ uniform [1,12]) [29,33].
In all instances, we investigated the effect of customer preference of the direct channel, θ, on pricing and lead time decisions. e results of some carried instances are summarized in Figures 3-9 and Tables 4 and 5. Tables 4 and 5 show the differences between average values of retailer prices (P As we can see from Figures 5 and 8, in both the centralized and decentralized dual-channel supply chain, when customer preference of the direct channel θ is below a certain level, the retail price is higher than the direct sale price; conversely, when θ exceeds that level, the direct price becomes higher than the retail sale price. is result shows that if the base level of demand or demand rate in one channel is relatively higher than a certain threshold, the sale price in that channel should be set higher than the one in the other channel.
From Figure 5, in some examples, when θ is relatively low, the direct sale price should be set equal to the wholesale price, and when θ is relatively high, the direct sale price should be set higher than the wholesale price. In other words, when θ is lower than a certain threshold, the base level of demand in the retail channel is high. us, the retail price can also be set to be high, but according to the low base demand in the direct channel, the direct sale price is set to be low. e wholesale price must be less than or equal to the direct sale price. erefore, the direct sale price is equal to the wholesale price, and the retail price is higher than the wholesale price; conversely, when θ exceeds the threshold, the base level of demand in the retail channel is    low, and the retail price is generally also set to be low. We know that the wholesale price must be lower than the retail price. erefore, with increasing θ, the retail price and consequently the wholesale price decrease. As we can see from Figures 5 and 8, for some values of θ, both the retail price and the direct price in the decentralized dual-channel supply chain should be set higher than those in the centralized one, while the total profit in the decentralized dualchannel supply chain is lower than that in the centralized one. e negative value of gap in price differences in Table 4 and the positive value of gap for profits in Table 5 show this fact. is shows that when each partner in the decentralized dual-channel supply chain maximizes his own profit, this leads to double marginalization. Double marginalization means that the retailer and the manufacturer independently set their price to maximize their profit margins; as a result, the price is higher and the sales volume and profits are lower than those of a vertically integrated channel [3].
In instances with low (Ce/Cr) or with high Cr, the direct sale price in the centralized setting tends to equal the direct sale price in decentralized settings. is trend is also true for retail prices in the two settings. Moreover, when Cr is very large, the retail channel has a small impact on the manufacturer's profits, and in both settings, the manufacturer has a main role in determining prices. erefore, the prices defined in the two settings become close.
As illustrated in Figure 7, in the centralized dual-channel supply chain, when customer preference of the direct channel θ is lower than a certain threshold, the total profit decreases with increasing θ; conversely, when θ exceeds the threshold, the total profit increases with increasing θ. is threshold differs from one instance to another. In examples with high direct channel operational cost, the threshold is near 1, meaning the function seems as decreasing function. To compare the profitability changes of the firm according to different direct and retail channel operational costs, we carried some more examples by varying these two parameters. Figure 10 shows the results. As we can see from Figure 10, under high value of direct channel operational cost, the profit decreases with increasing θ. erefore, it is more profitable to encourage customers buying from retail channel.
As we can see from Figures 6 and 9, the high lead time changes depending on θ in two centralized and decentralized setting show that the delivery lead time has a strong effect on the manufacturer's and the retailer's pricing strategies and profits. We can see the results of Propositions 4 and 5 in instances 2 and 16, where the unlimited capacity is considered. e quoted lead times for the customers are equal to the lower bound or the upper bound of the lead time.
We perform a numerical study to compare how the lead time flexibility affects the profitability of the firm.
We consider two policies of lead time flexibility: P1 (lead time flexibility), P2 (no lead time flexibility). With lead time flexibility, we quote different lead times to different customers. When there is no lead time flexibility, a single (fixed) lead time is quoted to all customers.
We have generated the problems in three situations: the customers' sensitivity to price in direct channel is 1-greater,   2-equal, and 3-lower than customers' sensitivity to price in retail channel. We have considered the same price sensitivities for all customers. e categorization for production capacity and variability of customers' lead time sensitivity is shown in Table 3. e percentage of profit increases over P2 policy with no lead time flexibility under each policy and category is represented in Table 6. e main conclusions one can draw from Table 6 are as follows: (1) e effect of production capacity: comparing instances 2 and 4 (instance 4 with high production capacity), we can see that more capacity makes the manufacturer able to charge equal lead times to all customers, and manufacturers usually determine a lead time equal to that of the first period for all customers.
(2) e effect of variability of customers' sensitivity to lead time: comparing instances 1 and 2 (instance 1 with high variability in customers' sensitivity), we can see that if there is high variability in customers' lead time sensitivities, the manufacturer can obtain more profit from lead time flexibility. e manufacturer can charge the high lead time for customer with lower sensitivity to lead time and reserve the capacity of the first period for customers with higher sensitivity to lead time.

Conclusion
In this paper, we presented a pricing and due date setting model for a manufacturer with a dual sale channel, an online direct channel, and a traditional retail channel. We see that a large number of e-retailers, such as Amazon.com, BestBuy. com, and Walmart.com, attempt to offer a proper delivery lead time (Yang et al. [4]). Also the model can be developed for 's auto industry which is the second largest industry in Iran after oil and gas. Currently, Iran's auto industry has an online sale and the proposed list price is sometimes based on delivery time.
We assumed a finite planning horizon divided into periods of equal length. We considered the limit production capacity in each period and multiple classes of customers arriving during these periods. We developed a Stackelberg game to model the decentralized dual-channel supply chain. Under this game framework, the supplier, as the leader, announces the wholesale price to an intermediary in addition to the direct channel sales price and lead time. e intermediary then reacts by choosing a retail price to maximize its own profit. Bilevel programming was developed to model this situation. For a centralized dual-channel supply chain, the traditional retail price, the direct sale price, and the quoted lead time in the direct channel are determined in an integrated manner. Considering the lead time selection for each customer order, mixed binary integer nonlinear programming models for the design of the centralized and decentralized dual-channel supply chains were used.
Because exact algorithms for solving the proposed models can be expensive and time consuming for instances with large numbers of 0-1 variables, the PSO algorithm solving model is adopted to get a satisfactory near-optimal solution efficiently.
rough numerical analyses, we have examined the effects of the customer preference of a direct channel on the manufacturers' and retailers' pricing behaviors. We also compared the optimal lead times, prices, and profits in the two settings of a centralized and decentralized dual-channel supply chain.
Numerical examples have shown that the retail price is decreasing with increasing customer preference of the direct channel and that the direct sale price is increasing with increasing customer preference of the direct channel. erefore, when the customer preference of the direct channel is low, the retail price is higher than the direct sale price; conversely, when the customer preference of the direct channel is high, the retail price is lower than the direct sale price. Comparing the two settings of dual-channel supply chains, we found that when θ is relatively low and, thus, the base level demand in the retail channel is high, both the retail price and the direct price in the decentralized dual-channel supply chain should be set to be higher than those in the centralized supply chain, while the total profit in the decentralized dual-channel supply chain is lower than that in the centralized supply chain. is shows double marginalization in the decentralized dual-channel supply chain which each member maximizes his own profit. e high changes in the lead time as a function of θ in the centralized and decentralized settings show that the delivery lead time strongly influences the manufacturers' and the retailers' pricing strategies and profits. erefore, considering a constant value for the lead time in all scenarios leads to a decrease in profits.
is research can be extended in several directions in future work. First, our model is deterministic, and we can study a model with stochastic demand. Second, we can consider competition among several manufacturers and retailers. ird, it is worth investigating the coordination of the dual-channel supply chain by contracts when the price and lead time are both considered in the models.
Obviously, Π r is a concave quadratic function of P r i . Using the partial first-order condition, we get erefore, the ] i is always negative which is a contradiction.

Proof of Proposition 3.
e retailer's best response to the wholesale price W i , the direct sale price P s i , and the direct sale quoted lead time L s i set by the manufacturer is given by (3).
According to unlimited capacity in each period, the manufacturer's profit is determined by Substituting (15) into (A.7) and simplifying it, we get Taking the second-order partial derivatives of Π s with respect to P s i and W i , we have the Hessian matrix as follows: From assumption 8b r i · b s i − (α s i ) 2 − (α r i ) 2 − 6α s i .α r i > 0, Π s is strictly jointly concave in P s i and W i , since (z 2 Π s /z(p s i ) 2 ) � − 2b s i + (α r i α s i /b r i ) < 0 and Proof of Proposition 4. As expressed in Proposition 3, under unlimited production capacity, the manufacturer's best price is obtained from e second-order partial derivative of Π s with respect to L s i is (z 2 Π s /z(L s i ) 2 ) � 0 which indicates that the Π s is a linear function of L s i . erefore, in this case, the optimal value of L s i will be equal to its upper bound (L u i ) or lower bound (L 0 i ).
Substituting P s i � 0 into the above inequality, we can derive (39) which is upper bound for L s i .

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that they have no conflicts of interest.