Optimal Replenishment for Perishable Products with Inventory-Dependent Demand and Backlogging under Continuous and Discrete Progressive Payments

Managing material and cash flows attracts concerns of physical and financial departments in most companies. We, therefore, focus on optimizing the replenishment policy for a channel with stock-dependent demand considering item deterioration and order backlogging under two financial schemes of progressive trade credit periods. We take both the continuous payment regime (CPR) and discrete payment regime (DPR) into account in the progressive trading process, which generates ten distinct scenarios. We show that the profit functions may not necessarily be concave and accordingly give a corresponding computing algorithm, which relaxes the convexity assumptions of objective functions in the existing literature and consequently enrich the research. We address the formulation characterization and the logic of pursuing global optimization from models arising in all settings. Computational studies and simulations are conducted to illustrate the effect of various parameters on the optimal replenishing policy and profit. Numerical experiments show that the CPR scheme is dominantly prior to DPR for long replenishment time intervals, whereas it is exactly the opposite for short time intervals. We also examine the impact of the shortage cost and deteriorating rate on optimum ordering policy and channel performance. Finally, future research directions are addressed in the end.


Introduction
Managing material and cash flows is the common concern for both physical and financial departments both in most of business organizations [1].In this sense, it is usually viewed as the interesting interface of financing and operations to seek trade-off between inventory control and accounting (including purchasing payment).For the transaction in supply chains, it is common in practice for the supplier to offer the downstream buyer different paying schemes after the transfer of the goods occurs, namely the Progressive Payment Scheme that has been extensively studied.This kind of strategy of payment delay is offered as an incentives of inspiring buyers to enlarge their purchasing amount, which is expected to incur similar effects resulting from a variety of incentive contracts [2] for enhancing supply chain performance.Aggarwal and Jaggi [3], Sarker et al. [4] and Soni and Shah [5] discuss the optimal ordering policy under the progressive payment scheme.In recent years, considering the progressive payment scheme, Pal [6] and Sen and Saha [7] also investigate the optimal decisions of supply chain.From here we see that the progressive payment scheme is one of the most significant factors in ordering decisions, which is also the key factor for the research in this paper.
For the characteristics of market demand, practitioners in practice and academia have both observed the phenomenon that the demand rate is usually affected by the inventory on hand, which has been addressed by dozens of studies since Levin et al. [8].Shukla et al. [9] developed two different mathematical models in terms of the different types of holding cost function with stock-level dependent demand for obtaining the optimal order quantity.Stock-dependent products is described as a key factor influencing the optimization problem in Fu and Chen [10], in which two functions of minimizing total cost with or without shortage are developed.Thus, it is important for a firm confronting the stock-dependent demand to find the optimum ordering policy associating its stock-on-hand evolving path in order to maximize (or minimize) the total profit (cost) in the long run.
On the other hand, some products may deteriorate depending on the storage conditions.The deteriorating items are classified into two categories in Li et al. [11], the first category are the decayed items, such as meat, fruits, vegetables, the other category are items that lose value, such as computer chips and mobile phone.The deteriorating products in this paper focus on the first category.As we all know, the deterioration of items shows a large impact on inventory decisions of firms.Therefore, how to control the items deterioration to mitigate its influence on inventory is a considerable question, which has gained much attention in academia such as Wang [12], Dye [13], Majumder et al. [14] and Ardak et al. [15].
We investigate the optimal reordering problem of deteriorating items with inventory-dependent demand in the scenario with a plausible delay payment scheme newly proposed.For this purpose, we clearly define two different delay payment regimes, Continuous Payment Regime (hereinafter referred as CPR) and Discrete Payment Regime (hereinafter referred as DPR) as follows: (1) Continuous Payment Regime: over some specific period, the buyer pays the supplier the money owed using the instantaneous revenue continuously until all purchasing cost has been paid back completely; (2) Discrete Payment Regime: the buyer accumulates revenue by depositing it in bank and pay money at some discrete specific time instants until all purchasing cost are cleared off.
The reminder of this paper is organized as follows.In Section 2, the relevant literatures are reviewed.In Section 3, we formulate the problem as well as define notations.In Section 4, we develop models in a variety of scenarios considering CPR and DPR in terms of the variation of permissible delay time, respectively.In Section 5, we firstly give a unified computational algorithm and then conduct computational experiments and sensitivity analyses to obtain some managerial insights.Future research and conclusions are offered in Section 6.

Literature Review
This study is typically featured in three aspects, namely, deterioration, stock-dependent demand and financially progressive payment credits.We therefore give a brief review over the literature from these three tracks.An increasing body of literature has focused on inventory control problems with deterioration in different settings.Wang [12], Dye [13] and Chang and [16] discussed the replenishment policy for deteriorating item considering backlogging.Yang et al. [17] formulate the inventory system for non-instantaneous deteriorating items with the price-dependent demand, aiming to find the optimal selling price and replenishment cycle time to maximize the total profit.Khanra et al. [18] assume the deterioration rate is constant, and developed the economic ordering quantity (EOQ) model for deteriorating item.Further, Singh and Pattnayak [19] explore the EOQ model by extending the aforesaid paper, suppose the deterioration rate is variable rather than constant.More studies related to the effect of deteriorating rate on inventory, and one can refer to Majumder et al. [14], Ardak et al. [15] and Pal and Samanta [20].
There is also a number of studies on modeling progressive payment schemes.Aggarwal and Jaggi [3], Sarker et al. [4] and Soni and Shah [5] are among early researches to develop some optimal ordering/quantity policies for system under progressive payment scheme.Chung [21] developed a new solution procedure to overcome the shortcoming of that in Goyal et al. [22] in terms of the classical EOQ model considering payment delay.Mohan et al. [23] explored the optimal replenishment policy for multi-item in three different situations: the permissible delay is negative, is less than the replenishment interval and is greater than the replenishment interval.Recently, Chen and Teng [24] proposed a simple arithmetic-geometric method to solve the extended EOQ model under a conditionally permissible delay in payments.The same year, Li et al. [25] describe the inventory games where the retailer's delay in payments is permitted by the supplier, and proved that the core of inventory games is nonempty under a permissible delay.For more research on the replenishment considering the permissible delay in payment also see Li et al. [1], Vandana and Sharma [26], Sen and Saha [7] and Pal [6].
For the track of stock-dependent demand, Min and Zhou [27] addressed a deterministic inventory model with demand backlogging permitted and exploited a solution procedure for optimum replenishment policy in a similar way as Lee and Dye [28].In a more complicated situation, Min et al. [29] considered sequential delayed payments between the retailer with her supplier and customers, respectively and focused on optimizing retailer's average profit.Although Sana [30] and Sana [31] also both considered stock dependent demand, they moved to study multiple products considering sales efforts' effect on demand as well.However, these two studies do not take delay payment into consideration while only showing the optimal replenishment policies.Extending the single-firm problem to the supply chain situation, Zhou et al. [32] analyzed the interactions of ordering and paying issues in a non-cooperative game theory framework.Their work opens the door to study stock-dependent demand featured ordering and trading policies with intertwined decision-makings in supply chain scenario.There are also a series of EOQ-based models with stock-dependent demand, such as Ray and Chaudhuri [33], Teng and Chang [34], Teng and Ouyang [35], Tyagi et al. [36], Shukla et al. [9], Ghosh et al. [37] and Fu and Chen [10].
Among the most relevant research to and our study, Teng et al. [38] provide a comprehensive extensions by allowing a non-zero ending inventory, bounded inventory capacity, and deterioration based on Soni and Shah [5].Glock et al. [39] correct some formulation errors in Soni and Shah [5] and research a new scenario in which the retailer's credit interest rate is greater than the supplier's charging interest rate.However, they focus on variant problems using a similar framework of modeling and a solving method with that in Soni and Shah [5].Different from the above studies, we not only simultaneously consider three factors of deterioration, stock-dependent demand and financially progressive payment credits, but also incorporate two delayed payment schemes and permit order backlogging, which is closer to reality yet seldom examined together.
Compared with the aforementioned existing literatures only adopting a DPR-like policy, our study considers CPR and DPR simultaneously in two main aspects: (1) under CPR, we initially formulate the reordering policy optimization problem, develop models and exploit associated solving processes; (2) under DPR, we reconstruct the formulation and modeling and moreover we permit the time point for clearing off purchasing cost (including charging interest) either in CPR or in DPR, which does not necessarily just hit the time instant T exactly as in Soni and Shah [5] and Teng et al. [38].Namely, for Z(T) the inventory level at time T equals to zero does not mean that the retailer needs to settle all unpaid money just at T. Relaxating the constraint of the clearing off time point at T is significant and will substantially affect the problem formulation, modeling and solving process.Second, when modeling and solving our problem, we remove the pre-conditions and assumptions that all objective functions are convex ensured by putting constraints upon the decision variable T which has been shown frequently in Soni and Shah [5], Teng et al. [38], Glock et al. [39] and so on.This relaxation is managerially significant again due to its impact on modeling and solving process and of course the optimal solution since the objective function should not be confined or necessarily convex (or concave) from the point of practice.

Preliminaries and Notations
Consider a retailer purchases finished goods from the upstream supplier in repeated cycles.The market demand facing the retailer is dependent on the timely stock level on hand, which is commonly assumed extensively in the literature and found in reality.Suppose the items are perishable goods that can deteriorate over time.Therefore, the inventory is depleted at three factors of stock level, deteriorating and constant consumption rates.Moreover, assume the retailer permits the order backlogging in each cycle and the unfilled orders will be served in the succeeding period.Hence, the schematic curve of inventory level over time can be plotted as shown in Figure 1, where the retailer receives Q + V units goods from the supplier at the beginning of each replenishment period (t = 0) and immediately put V units of products to meet the backlogging orders inherited from the preceding period, and accordingly Q units are held as the initial inventory of the present period.In each replenishment cycle, the retailer can earn interest income from depositing sales revenue in the bank until the deadline of the permitted credit periods.The time instant T c starts the shortage interval implying that the retailer needs to decide suitable V for backlogging all orders generated in this meantime.Consider a retailer purchases finished goods from the upstream supplier in repeated cycles.The market demand facing the retailer is dependent on the timely stock level on hand, which is commonly assumed extensively in the literature and found in reality.Suppose the items are perishable goods that can deteriorate over time.Therefore, the inventory is depleted at three factors of stock level, deteriorating and constant consumption rates.Moreover, assume the retailer permits the order backlogging in each cycle and the unfilled orders will be served in the succeeding period.Hence, the schematic curve of inventory level over time can be plotted as shown in Figure 1, where the retailer receives Q V  units goods from the supplier at the beginning of each replenishment period ( 0 t  ) and immediately put V units of products to meet the backlogging orders inherited from the preceding period, and accordingly Q units are held as the initial inventory of the present period.In each replenishment cycle, the retailer can earn interest income from depositing sales revenue in the bank until the deadline of the permitted credit periods.The time instant c T starts the shortage interval implying that the retailer needs to decide suitable V for backlogging all orders generated in this meantime.In addition to the ordering and replenishment problem mentioned above, we also incorporate the issue of payment in instalment for goods into the study, which interwines inventory ordering policy to substantially impact the system profit since the payment by instalment can also bring penalty and/or income.Therefore, considering the perishable goods with deterioration, we intend to investigate the system optimal ordering policy with shortage and stock-dependent demand under two trading regimes and credits.Specifically, several cases are discussed according to the variation of credit periods.
For the necessity of problem formulating, the following assumptions are adopted throughout the whole paper: (i) The demand rate, , where a, b and  are intrinsic demand, sensitivity parameter and deteriorating rate to inventory level, respectively, with and replenishment usually occurs at the beginning of each period.Furthermore, assume the supplier's capacity is infinite without defective item.(iii) as for the payment schedule, the retailer will charge additional fees for the unpaid balances as follows: nothing in the case of payment clearance before time M, at a rate of

Mathematical Models
According to the illustrations and assumptions in Section 2, the depletion of inventory is due to three aspects: the constant demand a , the stock-depended demand and the deterioration in the products.Thereby, we use the differential equations below to govern the instantaneous state of inventory at time t In addition to the ordering and replenishment problem mentioned above, we also incorporate the issue of payment in instalment for goods into the study, which interwines inventory ordering policy to substantially impact the system profit since the payment by instalment can also bring penalty and/or income.Therefore, considering the perishable goods with deterioration, we intend to investigate the system optimal ordering policy with shortage and stock-dependent demand under two trading regimes and credits.Specifically, several cases are discussed according to the variation of credit periods.
For the necessity of problem formulating, the following assumptions are adopted throughout the whole paper: (i) The demand rate, R(t) = a + bZ(t) + θZ(t), is a function of the inventory level Z(t), where a, b and θ are intrinsic demand, sensitivity parameter and deteriorating rate to inventory level, respectively, with a > b > 0, 0 ≤ θ < 1 and R(t) = a when T c ≤ t ≤ T. (ii) Lead-time is neglected and replenishment usually occurs at the beginning of each period.Furthermore, assume the supplier's capacity is infinite without defective item.(iii) as for the payment schedule, the retailer will charge additional fees for the unpaid balances as follows: nothing in the case of payment clearance before time M, at a rate of I c 1 supposing payment clearance during time interval [M, N] and a rate of I c 2 given payment clearance after N with I c 2 > I c 1 > I e and I e denoting interest the rate of banks.

Mathematical Models
According to the illustrations and assumptions in Section 2, the depletion of inventory is due to three aspects: the constant demand a, the stock-depended demand and the deterioration in the products.Thereby, we use the differential equations below to govern the instantaneous state of inventory at time t Then the stock level can be obtained below in view of the boundary conditions Z(T c ) = 0 and Combining the condition Z(T) = −V and Function (4) yields Thus, the inventory level in period [0, T] can be summarized as below Additionally, the retailer's other costs per unit time after value assignment include ordering cost A/T, shortage cost c 0 V/T, and holding cost In view of diverse credit periods with different penalty interests charged by supplier as mentioned previously, four major cases with subcases will be discussed as follows, after which we will compare all cases under different progressive payment schemes to recognize the optimum replenishing policy.

Case 1 T c < T ≤ M < N
We first consider an extreme but simplest situation where the replenishment cycle is less than the first credit time M (shown in Figure 2).In this case, assume the retailer intends to exploit the advantage of no penalty incurred in the case of payment clearance before the first credit time M. Thus, the rational retailer will definitely clear all outstanding payment at M since she can deposit the sales revenue in the bank to earn interest at the rate of I e .Considering the costs incurred by ordering, wholesaling, shortage and inventory holding, the retailer's problem is to determine suitable replenishing cycle time and shortage to maximize her annual profit, which can be expressed as: The first term, the sum of second and third terms and the rest in the curly braces of above Equation ( 7) represent sales revenue, interest earned and manifold aforesaid costs, respectively.
We now consider the situation in which the replenishment time interval can be divided into two parts: [0, M] and [M, T] (shown in Figure 3).The profit of retailer includes both sales revenue and interest earned during the credit period [0, M].While penalty interest is charged at rate of I c 1 during in the period [M, T], the wise retailer will repay the wholesale cost at the time M to avoid the penalty interest charged by the supplier.The interest earned during the credit period [0, M] takes the form:   [0, ] 0 ( )( ) The remaining deposit of the retailer after repaying the wholesale cost to the supplier at time M can be expressed as follows: The first two terms of the Equation ( 9) represents the sales revenue, the sum of the third and fourth terms denotes the earned interest of the retailer by depositing the sales revenue in the bank, and the last term of the function represents the wholesale cost.
It should be noted that whether the revenue including sales revenue and earned interest is sufficient to repay all the wholesale cost.Considering the problem, two subcases are discussed below.

Subcase 2-1 The Revenue Being Sufficient to Repay the Wholesale Cost at Time M
The retailer repays all the wholesale cost at time M , according to the Equation ( 9), immediately, the annual profit of retailer in the whole replenishment time interval takes the form:   [0, ] 0 ( )( ) The remaining deposit of the retailer after repaying the wholesale cost to the supplier at time M can be expressed as follows: The first two terms of the Equation ( 9) represents the sales revenue, the sum of the third and fourth terms denotes the earned interest of the retailer by depositing the sales revenue in the bank, and the last term of the function represents the wholesale cost.
It should be noted that whether the revenue including sales revenue and earned interest is sufficient to repay all the wholesale cost.Considering the problem, two subcases are discussed below.

Subcase 2-1 The Revenue Being Sufficient to Repay the Wholesale Cost at Time M
The retailer repays all the wholesale cost at time M , according to the Equation ( 9), immediately, the annual profit of retailer in the whole replenishment time interval takes the form: The remaining deposit of the retailer after repaying the wholesale cost to the supplier at time M can be expressed as follows: The first two terms of the Equation ( 9) represents the sales revenue, the sum of the third and fourth terms denotes the earned interest of the retailer by depositing the sales revenue in the bank, and the last term of the function represents the wholesale cost.
It should be noted that whether the revenue including sales revenue and earned interest is sufficient to repay all the wholesale cost.Considering the problem, two subcases are discussed below.

Subcase 2-1 The Revenue Being Sufficient to Repay the Wholesale Cost at Time M
The retailer repays all the wholesale cost at time M, according to the Equation ( 9), immediately, the annual profit of retailer in the whole replenishment time interval takes the form:

Subcase 2-2 The Revenue Is Insufficient to Repay the Wholesale Cost at Time M
In this subcase, we consider that the revenue is insufficient to repay the wholesale cost in the permissible credit period.For the wise retailer, all revenue including sales revenue and earned interest will be spent on repaying the wholesale cost at time M to avoid paying higher penalty interest.The unpaid money accordingly at time M is shown as: Obviously, the unpaid money must be paid by retailer and the penalty interest will be incurred beyond the credit period M.Meanwhile, the supplier provides two ways for the retailer regarding how to pay the remaining balance: continuous and discrete payment regimes.

Subcase 2-2a Continuous Payment Regime (CPR)
With the continuous payment regime, the retailer continuously uses the instantaneous revenue to repay the debt until all wholesale costs are completely paid off.Denote t a the time at which the unpaid balance is cleared off under the continuous payment regime through the paper, M < t a ≤ T c .Especially, the profit of the retailer in the period (T c < t a ) is negative, to facilitate the research, the situation is excluded in this model.
The unpaid money during the [M, T c ] can be calculated from Combining the Equations ( 12) and ( 13), the unpaid money takes the form: The clearance point t a can be obtained by letting U(t) = 0. Furthermore, the earned interest in the interval [t a , T c ] is Immediately, the annual profit of the retailer during the whole cycle period can be expressed as: Subcase 2-2b Discrete Payment Regime (DPR) Different from the continuous payment regime, the retailer will pay all the remaining debt at the specific time t b with the discrete payment regime.Similarly, the case with T < t b is not discussed in the paper.
The charged interest during the period [M, t b ] can be written as: The earned interest during the period [M, t b ] can be written as: The revenue during the period [t b , T c ] can be written as: Considering the retailer settles all unpaid balance at time t b , the following equation holds The left-hand side (LHS) of above equation represents the sales revenue and interest earned in the interval [M, t b ], and the right-hand side (RHS) denotes the unpaid balance and interest charged.Solving the equation below can obtain t b where φ(M) < 0, φ(t b ) = 0, φ(T c ) > 0. Subsequently, the annual profit of the retailer during the whole cycle period is shown as: Consider the different penalty interest rates over progressive delays in payments with the relation M < N < T c < T (shown in Figure 4).In this case, for necessity of analysis the whole cycle period is divided to three intervals [0, M], [M, N] and [N, T], each of which corresponds to a subcase as discussed below.Similarly, whether the retailer can settle all unpaid balance at time M generates the following two subcases.

Subcase 3-1 The Revenue Is Sufficient to Repay the Wholesale Cost at Time M
This case is no different from the Case 2-1 except for the credit period, which will not be reiterated here.

Subcase 3-2 The Revenue Is Insufficient to Repay the Wholesale Cost at Time M
Consider the subcase in which total purchasing cost cannot be settled at time M.As mentioned previously, assume the supplier offers the retailer two progressive payment regimes, i.e. continuous and discrete payment regimes, whose essential difference reflects the way the retailer behaves continuously or discretely for doing the payment.Then we explore the subcases below.With the continuous payment regime, the retailer continuously transfers the instantaneous income to defray the debt until all balance is cleared off at a certain time t a .However, there are different penalty interests incurred when the clearance point t a belongs to the intervals [M, N] or [N, T c ] due to the principles described in supplier's installment payment schemes.Hence, the Subcase 3-2a can be further decomposed into two separate parts below.Subcase 3-2a' M < t a ≤ N This subcase is similar to the Subcase 2-2a as described previously, so it will not be addressed here to avoid the repetition.
In this subcase, the retailer is not able to settle all unpaid balance in the two progressive periods.Under continuous payment regime, the retailer will pay the penalty interest at the rate I c 1 during the time interval [M, N] and I c 2 during the period [N, T c ]. Then the differential equations below hold According to the Equation ( 12) and conditions U(N) > 0 and U(T c ) ≤ 0, the unpaid balance is governed by the equations below Hence, retailer's annual profit in whole period in Subcase 3-2a" can be expressed as: Subcase 3-2b Discrete Payment Regime (DPR) Now we pause here to discuss the subcase under the other progressive payment regime, i.e. discrete payment regime.As the name suggests, DPR is just on the opposite of CPR, namely, the retailer deposits money in some discrete times until all unpaid balance is cleared off.at time t b .Similarly, depending on intervals [M, N] or [N, T c ], to which the point t b belongs, we further discuss the problem in two separate situations, respectively.Subcase 3-2b' M < t b ≤ N Similarly, this subcase is the same as the Subcase 2-2b, so we again omit the deriving process.

Subcase 3-2b" N < t b ≤ T c
Considering the penalty interest charged by the supplier, the retailer tries to settle the unpaid balance by time N and then he pays off all the debt at the time t b under the discrete payment regime.Hence, the unpaid balance at time N can be written as: Accordingly, the interest charged and the interest earned during the credit period [N, t b ] are shown respectively as: Note that the retailer can square up all unpaid balance at the time t b only when her accumulated income equals the sum of unpaid money and charged interest, that is the equation below holds where t b is determined by the following function as well as its boundary conditions with Given the above, the retailer's average annual profit can be expressed as follows

Max
(T,V) Sustainability 2018, 10, x FOR PEER REVIEW 10 of 24 Considering the penalty interest charged by the supplier, the retailer tries to settle the unpaid balance by time N and then he pays off all the debt at the time b t under the discrete payment regime.Hence, the unpaid balance at time N can be written as: Accordingly, the interest charged and the interest earned during the credit period [ , ] b N t are shown respectively as: Note that the retailer can square up all unpaid balance at the time b t only when her accumulated income equals the sum of unpaid money and charged interest, that is the equation below holds where b t is determined by the following function as well as its boundary conditions x x with. ( ) 0, ( ) 0, ( ) Given the above, the retailer's average annual profit can be expressed as follows

Case 4 ( ) c c c T M T T M N T or T M T N        
Now we shed the light on the case in which the relationship c T M T   holds, which can be 5).In view of the retailer's rationality, she will return a large amount money at time M while she earns income from sales and deposit interest.Then retailer's remaining bank savings after payment at time M can be expressed as follows:

Case 4 T c <
Now we shed the light on the case in which the relationship T c < M < T holds, which can be decomposed into sub-cases T c < M < N < T or T c < M < T < N (shown in Figure 5).In view of the retailer's rationality, she will return a large amount money at time M while she earns income from sales and deposit interest.Then retailer's remaining bank savings after payment at time M can be expressed as follows: where the first term in RHS of Equation ( 32) represents the sales revenue, the last term the total purchasing cost and the second term the earned interest of the retailer by depositing the sales revenue in the bank.Furthermore, the remaining savings of retailer is continuously deposited in bank to earn the interest during the credit period [M, T], and it takes the form Hence, the average annual profit of retailer takes can be formulated as follows Sustainability 2018, 10, x FOR PEER REVIEW 11 of 24 where the first term in RHS of Equation ( 32) represents the sales revenue, the last term the total purchasing cost and the second term the earned interest of the retailer by depositing the sales revenue in the bank.Furthermore, the remaining savings of retailer is continuously deposited in bank to earn the interest during the credit period [ , ] M T , and it takes the form 4 ( )( ) Hence, the average annual profit of retailer takes can be formulated as follows Figure 5. Inventory level over time in case 4.

The Characterization and Logic of Global Optimization
Given all the formulation above, there exist, in total, ten subcases explored in terms of two progressive payment approaches and a comparison of trade credits and the replenishing cycle, which is shown below:

The Characterization and Logic of Global Optimization
Given all the formulation above, there exist, in total, ten subcases explored in terms of two progressive payment approaches and a comparison of trade credits and the replenishing cycle, which is shown below: Essentially, it is a process to pursue global optimization based on multiple piecewise nonlinear optimizations.Its logic is to first solve each subcase by getting the local optimum solution, and then compare all of these ten local optimum values to determine the global optimum solution.In view of the complication and nonlinearity of the whole problem, we are not able to give the elegantly analytical solution, so an effective way is to conduct numerical experiments to investigate the problem's characterization as well as managerial insights.Note that the detailed exposition of formulation and representations for all ten cases are addressed in the Appendix.

Numerical Examples
To illustrate the previously developed models, the paper utilizes Particle Swarm Optimization (PSO), an extensively used intelligent optimization algorithm, to solve the problem.For more detailed information about the PSO, one can refer to Shi and Eberhart [40] and Eberhart and Shi [41].Here some computational experiments are proposed under given parameters below, A = 200, c = 20, c 0 = 3, p = 30, h = 6, a = 1000, b = 0.25, M = 17/365, N = 35/365, θ = 0.1, I e = 12%, I c 1 = 13% and I c 2 = 20%.In addition, assume the upper bounds of shortage quantity and replenishment time interval are V 0 = 500 and T 0 = 2, respectively.The scale of particle swarm is 100, the iteration times 300 and the calculation time is 10.
Example 1.The annual profit (maximum between CPR and DPR) of the retailer is increasing in the replenishment time interval originally, however, with the growth of the replenishment time interval, the annual profit presents a downtrend.This is the same as our intuition that the annual profit is decreasing in shortage quantity.In addition, the annual profit is zero in some situations as shown in Figure 6.The reason for rendering the phenomenon is that the constraint condition T c = T − V/a ≥ 0 is not satisfied.In view of the situation, we define the annual profit equals to zero when T c = T − V/a < 0. Observing Figure 7, when there is a short replenishment time interval, the profit of the retailer in the discrete payment regime (DPR) is greater than that in the continuous payment regime (CPR), while CPR is superior to DPR with the long replenishment time interval.The result reveals that it is beneficial for the retailer to choose DPR with a short replenishment time interval, otherwise, the CPR will be the best choice.
Essentially, it is a process to pursue global optimization based on multiple piecewise nonlinear optimizations.Its logic is to first solve each subcase by getting the local optimum solution, and then compare all of these ten local optimum values to determine the global optimum solution.In view of the complication and nonlinearity of the whole problem, we are not able to give the elegantly analytical solution, so an effective way is to conduct numerical experiments to investigate the problem's characterization as well as managerial insights.Note that the detailed exposition of formulation and representations for all ten cases are addressed in the Appendix.

Numerical Examples
To illustrate the previously developed models, the paper utilizes Particle Swarm Optimization (PSO), an extensively used intelligent optimization algorithm, to solve the problem.For more detailed information about the PSO, one can refer to Shi and Eberhart [40] and Eberhart and Shi [41].Here some computational experiments are proposed under given parameters below, =200 A , 20 c  , . In addition, assume the upper bounds of shortage quantity and replenishment time interval are 0 500 V  and 0 2 T  , respectively.The scale of particle swarm is 100, the iteration times 300 and the calculation time is 10.Example 1.The annual profit (maximum between CPR and DPR) of the retailer is increasing in the replenishment time interval originally, however, with the growth of the replenishment time interval, the annual profit presents a downtrend.This is the same as our intuition that the annual profit is decreasing in shortage quantity.In addition, the annual profit is zero in some situations as shown in Figure 6.The reason for rendering the phenomenon is that the constraint condition 0 c T T V a    is not satisfied.In view of the situation, we define the annual profit equals to zero when 0 c T T V a    .Observing Figure 7, when there is a short replenishment time interval, the profit of the retailer in the discrete payment regime (DPR) is greater than that in the continuous payment regime (CPR), while CPR is superior to DPR with the long replenishment time interval.The result reveals that it is beneficial for the retailer to choose DPR with a short replenishment time interval, otherwise, the CPR will be the best choice.Example 2. Figures 8 and 9 depict how the unit shortage cost and deteriorating rate influence on the retailer's profit in the two different payment regimes, and some results can be obtained from the figures.First, the profit is decreasing in unit shortage cost (the same result described aforesaid) and deteriorating rate.Second, the time hitting the profit peak shifts right with the increase of shortage quantity, that is, the retailer obtains the maximum profit in a short replenishment time interval with a small shortage quantity.Third, originally, the profit of retailer in DPR is higher than that in CPR, while, as the growth of replenishment time interval, the CPR is superior to the DPR.The result illustrates when there is a short replenishment time interval, and it is beneficial to the retailer with DPR, otherwise, the CPR will be chosen by the retailer (the same result as obtained in Example 1).Example 2. Figures 8 and 9 depict how the unit shortage cost and deteriorating rate influence on the retailer's profit in the two different payment regimes, and some results can be obtained from the figures.First, the profit is decreasing in unit shortage cost (the same result described aforesaid) and deteriorating rate.Second, the time hitting the profit peak shifts right with the increase of shortage quantity, that is, the retailer obtains the maximum profit in a short replenishment time interval with a small shortage quantity.Third, originally, the profit of retailer in DPR is higher than that in CPR, while, as the growth of replenishment time interval, the CPR is superior to the DPR.The result illustrates when there is a short replenishment time interval, and it is beneficial to the retailer with DPR, otherwise, the CPR will be chosen by the retailer (the same result as obtained in Example 1).Example 2. Figures 8 and 9 depict how the unit shortage cost and deteriorating rate influence on the retailer's profit in the two different payment regimes, and some results can be obtained from the figures.First, the profit is decreasing in unit shortage cost (the same result described aforesaid) and deteriorating rate.Second, the time hitting the profit peak shifts right with the increase of shortage quantity, that is, the retailer obtains the maximum profit in a short replenishment time interval with a small shortage quantity.Third, originally, the profit of retailer in DPR is higher than that in CPR, while, as the growth of replenishment time interval, the CPR is superior to the DPR.The result illustrates when there is a short replenishment time interval, and it is beneficial to the retailer with DPR, otherwise, the CPR will be chosen by the retailer (the same result as obtained in Example 1).   1, and some information is reflected in this phenomenon.With the increases of the unit shortage cost, the retailer decreases the replenishment time interval and order quantity by increasing the replenishment times for reducing the shortage cost.Nevertheless, the profit of the retailer is still decreasing in the shortage cost.In addition, the shortage quantity closes to the maximum shortage quantity with a low shortage cost.However, there is no shortage with the high shortage cost.  2 and 3, some conclusions can be found.As the increase of stockdependent demand parameter b and the credit periods M and N , the annual profit of retailer is increasing, while it is decreasing in penalty interest rates that when the demand is more dependent on the stock-level or there is a long credit period, the retailer should increase the ordering quantity and extend the replenishment time interval for getting high profits.However, considering the penalty interest, the retailer must shorten the replenishment time interval and reduce the ordering quantity to mitigate the negative effect of penalty interest.In reality, firms should focus on two issues for getting high profit.First, investigate the sensitivity of the consumer on inventory to evaluate the value of b , then decide to increase or reduce the order quantity and replenishment cycle interval for profiting more in term of the value of b .Second, retailers should strengthen the bargaining power with upstream companies to get the long credit periods and low penalty interest rate.Example 3. Replenishment time interval, order quantity and profit of the retailer are decreasing in unit shortage cost, as shown in Table 1, and some information is reflected in this phenomenon.With the increases of the unit shortage cost, the retailer decreases the replenishment time interval and order quantity by increasing the replenishment times for reducing the shortage cost.Nevertheless, the profit of the retailer is still decreasing in the shortage cost.In addition, the shortage quantity closes to the maximum shortage quantity with a low shortage cost.However, there is no shortage with the high shortage cost.  2 and 3, some conclusions can be found.As the increase of stock-dependent demand parameter b and the credit periods M and N, the annual profit of retailer is increasing, while it is decreasing in penalty interest rates I c 1 and I c 2 .The phenomenon indicates that when the demand is more dependent on the stock-level or there is a long credit period, the retailer should increase the ordering quantity and extend the replenishment time interval for getting high profits.However, considering the penalty interest, the retailer must shorten the replenishment time interval and reduce the ordering quantity to mitigate the negative effect of penalty interest.In reality, firms should focus on two issues for getting high profit.First, investigate the sensitivity of the consumer on inventory to evaluate the value of b, then decide to increase or reduce the order quantity and replenishment cycle interval for profiting more in term of the value of b.Second, retailers should strengthen the bargaining power with upstream companies to get the long credit periods and low penalty interest rate.

Concluding Remarks
In this research, we investigated the optimizing ordering policy for a perishable product supply chain with stock-dependent demand facing items deterioration and order backlogging under progressive payment schemes.We considered the continuous payment regime (CPR) and discrete payment regime (DPR) both in the progressive trading process to generate piecewise nonlinear optimization problem of ten distinct scenarios.We demonstrated the characterization of problem formulation as well as the logic of exploring global optimization from all piecewise local optimizations.We further conducted computational studies and numerical simulation to illustrate the impact of parameters on the optimal ordering policy and channel performance.
Some key findings of this study are summarized as follows: (i) we decompose the problem into an equivalent set of multiple piecewise nonlinear optimizations and elaborate its formulation characterization and logic, which provides an integral approach for exploring the optimal replenishing policy in this specific problem.(ii) comparing the profits incurred by two progressive payment schemes implies that the relatively long replenishment cycle will let the retailer better off under CPR, whereas the relatively short cycle works better under DPR.Therefore, for firms, it is best choice to repay the money by CPR when there is long replenishment cycle, otherwise, the DPR will be the best way for getting more profits.(iii) reducing order backlogging quantity may benefit retailer's pursuing greater profit in short time range.Thus, the one who wants to maximize the profit in a short time would better reduce order backlogging quantity.(iv) the increasing of stock-dependent demand parameter and the length of the progressive credit period will push up the ordering quantity and replenishment cycle for getting more profit.The result implies two directions that the firms should pay much attention: first, increasing the sensitivity of the consumer on inventory, and second, enhancing bargaining power to obtain a long progressive credit period.
Our research contributes to the track of managing material and cash flows.The proposed model can be extended in several ways.For example, the model may extend the constant deterioration rate to time-varying deterioration.Furthermore, retailers are also likely to provide the permissible delay in payment to consumers in day-to-day deal; then the model will become more complex, and more significant results may be explored.

Figure 1 .
Figure 1.The change of the inventory level over time.

1 cI 2 cI
supposing payment clearance during time interval [M, N] and a rate of given payment clearance after N with the rate of banks.

Figure 1 .
Figure 1.The change of the inventory level over time.

)I
Sustainability 2018, 10, x FOR PEER REVIEW 6 of 24 interest earned during the credit period [0, ] M .While penalty interest is charged at rate of 1 c during in the period [ , ] M T , the wise retailer will repay the wholesale cost at the time M to avoid the penalty interest charged by the supplier.The interest earned during the credit period [0, ] M takes the form:

Figure 2 .
Figure 2. Inventory level over time in case 1.

Figure 3 .
Figure 3. Inventory level over time in case 2.

Figure 2 .
Figure 2. Inventory level over time in case 1.

Figure 2 .
Figure 2. Inventory level over time in case 1.

Figure 3 .
Figure 3. Inventory level over time in case 2.

Figure 3 .
Figure 3. Inventory level over time in case 2.

Figure 4 .
Figure 4. Inventory level over time in case 3.

Figure 4 .
Figure 4. Inventory level over time in case 3.

Figure 5 .
Figure 5. Inventory level over time in case 4.

Figure 6 .
Figure 6.The annual profit of retailer.Figure 6.The annual profit of retailer.

Figure 6 .
Figure 6.The annual profit of retailer.Figure 6.The annual profit of retailer.

Figure 7 .
Figure 7.The difference of the annual profit between CPR and DPR.

Figure 7 .
Figure 7.The difference of the annual profit between CPR and DPR.

Figure 7 .
Figure 7.The difference of the annual profit between CPR and DPR.

Table 1 .
Effect of 0 c on decision parameters.

Table 2 .
Effect of b , M and

Table 1 .
Effect of c 0 on decision parameters.

Table 2 .
Effect of b, M and I c 1 on decision parameters.

Table 3 .
Effect of b, N and I c 1 on decision parameters.