key: cord-0446378-tayd034o authors: Chen, Zhen; Zhang, Ren-qian title: A multi-period multi-product stochastic inventory problem with order-based loan date: 2020-12-09 journal: nan DOI: nan sha: 697a984412c1af6a48f0459cd0954884261023e7 doc_id: 446378 cord_uid: tayd034o This paper investigates a multi-product stochastic inventory problem in which a cash-constrained online retailer can adopt order-based loan provided by some Chinese e-commerce platforms to speed up its cash recovery for deferred revenue. We first build deterministic models for the problem and then develop the corresponding stochastic programming models to maximize the retailers' expected profit over the planning horizon. The uncertainty of customer demand is represented by scenario trees, and a scenario reduction technique is used to solve the problem when the scenario trees are too large. We conduct numerical tests based on real data crawling from an online store. The results show that the stochastic model outperforms the deterministic model, especially when the retailer is less cash-constrained. Moreover, the retailer tends to choose using order-based loan when its initial available cash is small or facing long receipt delay length. Cash availability is very important to small-to-medium-sized enterprises (SMEs), including start-up companies and small retailers. They usually lack sufficient capital to absorb large losses, and it is difficult for them to obtain external financing such as loans or venture capital compared with large companies due to production capacity and sales scale. Research from the Association of Chartered Certified Accountants (ACCA) and the Institute of Management Accountants (IMA) (Doove et al., 2014) found that access to growth capital was one of the two main factors affecting the businesses of SMEs in the UK over the past two years. A report by research firm CB Insights found that 29% of start-ups failed because of cash crises (CBInsights, 2018) . Moreover, SMEs tend to be the first to feel the effects of financial crises such as the COVID-19 pandemic (Caniato et al., 2020) . On Chinese e-commerce platforms such as Taobao.com of Alibaba Group and Jd.com of JD Group, there are many small online retailers operated by only a few people or even a single individual. Generally, the maximum number of products they can purchase depends on their available cash in the current period, which are similar to nanostores in emerging markets (Boulaksil and van Wijk, 2018) . Many forms of supply chain financing services appear in the market to help SMEs alleviate the cash shortage problem. Thomas Olsen of Bain, a consulting firm, reckons that the supply chain finance market is expanding by 15-25% per year in the Americas and by 30-50% in Asia (Economist, October 2017) . Some Chinese e-commerce platforms have also started supply chain financing services for their online retailers in recent years. Alibaba Group opened its own online bank, mybank.cn in 2010, which mainly provides loans for online stores on its three ecommerce websites: alibaba.com, tmall.com and taobao.com. Its main competitor in China-jd.com-also set up a supply chain finance department in 2012 and provides similar financing services. Some Chinese e-commerce platforms provide a financing service called order-based loans that has seldom been addressed in the literature. The aim of this service is to speed up the retailer's deferred revenue, which is common in Chinese online shopping. In the real transactions of Chinese e-commerce platforms, payments from customers are first obtained by the platform. The platform then transfers the payment to the retailer's account after customers confirm the receipt of the order with additional processing time. Therefore, the retailer actually receives the customer's payment after a time lag (usually at least 2 weeks). The transaction process on Taobao.com or Tmall.com is shown in Figure 1 , where Alipay is the online payment system of Alibaba Group. The circled numbers in the figure represent the sequence of the transaction process. Order-based loan is a kind of financing services intended to help speed up retailers' cash recovery. Once it is adopted for an order, the platform immediately makes payments to the retailer on behalf of customers after the order is marked as delivered. The retailer repays the loan with some interest to the platform after an agreed time limit. The retailer does not need to wait for its payment, which is a good method to alleviate the retailer's cash shortage. The payment delay length is zero when applying order-based loan to an order. Given the above background, the aim of this paper is to investigate the impacts of order-based loan on an online retailer's operational decisions. We make the following contribution to the literature. • We consider a supply chain financing service: order-based loan in a stochastic inventory problem, which has been used in business practice but seldom addressed in the inventory literature. • A stochastic multi-product inventory model is built for the problem. Scenario trees are adopted to represent the demand uncertainty, and we use a scenario tree reduction method to solve the problem if the tree size is large. • By crawling real transaction data from an online store and performing numerical analysis, we find that factors such as the initial available cash, receipt delay length, overhead cost and demand fluctuations may influence a retailer's decision to apply for order-based loans. The literature associated with our work can be classified into two main streams: inventory management with financing decisions and scenario programming in operational research. Financial flow management is an important part of supply chain management (Cooper et al., 1997) . In recent years, there have already been many works addressing joint inventory management and financing decisions. Here, we only review some of them for brevity. Buzacott and Zhang (2004) first considered operational and financial decisions by analyzing a cash-constrained retailer with asset-based financing provided by a bank. Dada and Hu (2008) built a Stackelberg gaming model involving a bank, manufacturer and cash-constrained retailer. Raghavan and Mishra (2011) considered a problem with a retailer and manufacturer facing cash constraints. Kouvelis and Zhao (2012) investigated the structure of optimal trade credit contracts between a supplier and bank. Moussawi-Haidar and Jaber (2013) developed a model considering delayed payments for a cash-constrained retailer. Tunca and Zhu (2018) discussed the role and efficiency of buyer intermediation in supplier financing through a game-theoretic approach. Wu et al. (2019) proposed a trade credit model with asymmetric competing retailers in which one weak retailer was capital-constrained. Jin et al. (2019) discussed noncollaborative and collaborative financing. Xu and Fang (2020) investigated a supply chain financing system with one supplier and one emission-dependent and capital-constrained manufacturer. Yuan et al. (2020) integrated supply risk into a cash constraint problem. Zhou et al. (2020) investigated manufacturer guarantor financing and third-party logistics guarantor financing in a four-party supply chain game. The abovementioned works mainly built single-period inventory models to analyze problems. In terms of the multi-period setting, Chao et al. (2008) investigated a multi-period self-financing newsvendor problem and proved that the optimal ordering pattern is an approximate base stock policy. Gong et al. (2014) further considered shortterm financing. Katehakis et al. (2016) showed that the optimal ordering policy is characterized by a pair of threshold parameters for this problem. Boulaksil and van Wijk (2018) formulated a cash-constrained stochastic inventory model with consumer loans and supplier credits and obtained managerial insights by simulating numerical cases. Chen and Zhang (2018) incorporated cash flow constraints into a multi-period lot-sizing problem with trade credit. Chen and Zhang (2019) extended this problem with credit-based loans. Bi et al. (2020) addressed short-term financing in an inventory problem with multi-item products. Li et al. (2020) considered a two-stage inventory model with financing and demand updating using a Bayesian approach. Chen and Rossi (2020) developed an (s, C(x), S ) policy for a stochastic inventory problem with fixed costs and cash constraints. In decision making under uncertainty, one method is scenario programming in which the uncertainty is represented by scenario trees. Høyland and Wallace (2001) presented a method to generate a limited number of scenarios that satisfy specified statistical properties. Haugen et al. (2001) adopted a progressive hedging algorithm for each scenario to solve a stochastic lot-sizing problem. Heitsch and Römisch (2003) developed two scenario reduction algorithms: new versions of forward and backward type algorithms. Brandimarte (2006) proposed a fix-and-relax strategy to solve a scenario programming problem with multi-item capacitated lot sizing. Beraldi et al. (2006) also made use of this method to solve a lot-sizing problem with uncertain processing times. Helber et al. (2013) took a scenario approach to solve a capacitated lot-sizing problem with random demand and dynamic safety stocks. Feng and Ryan (2013) used the scenario generation and reduction method for a stochastic power generation expansion planning problem. Hu and Hu (2016) and Hu and Hu (2018) adopted similar methods to solve two-stage and multi-stage stochastic problems, respectively. Hnaien and Afsar (2017) considered a minmax robust lot-sizing problem with uncertain lead time. Fattahi et al. (2018) used the forward scenario construction technique and Benders decomposition algorithm to solve a multi-period supply chain network redesign problem. To the best of our knowledge, there is no work in the literature that considers order-based loans in a multi-period stochastic inventory problem, which inspired us to investigate this topic. For convenience, some notations adopted in this paper are listed in Table 1 . Other relevant notations will be introduced as needed. Table 1 : Some notations adopted in the paper. t Index of a period, t = 1, . . . , T . n Index of an item, n = 1, 2, . . . , N. Operational parameters: I n,0 Initial inventory for item n. d n,t Demand for item n in period t, which is uncertain v n Unit variable ordering cost of the retailer for item n. Financial parameters: C 0 Initial cash at the beginning of planning horizon. Overhead costs in each period (e.g., wages or rents) L Receipt delay length. p n Selling price of item n. v n Unit variable ordering cost of item n. α Discount rate for unreceived revenues at the end of planning horizon. The amount of order-based loan. r o Interest rate of order-based loan. Decision variables: Q n,t Ordering quantity for item n in period t. The quantity of order-based loan used for product n in period t. R n,t Retailer's revenue from item n in period t. End-of-period inventory level for item n in period t. C t End-of-period cash balance in period t. FC Final cash at period T after discounting unreceived revenues. Our study is about a cash-constrained online retailer selling multiple products on Tmall.com, which is a Chinese e-commerce platform owned by the Alibaba Group. In the problem, the retailer has initial cash C 0 before making operational decisions. There are N products for sale, and its planning horizon length is T . The initial inventory for product n at the beginning of the planning horizon is I n,0 . Overhead costs H t , such as wages or rents, may have to be paid in each period, which are incurred irrespective of the retailer's operation activity. The retailer purchases products from its suppliers, and the unit ordering cost for product n is v n . The retailer faces uncertain demand from customers. The demand for product n in period t is represented by d n,t . In each period, the retailer makes decisions about Q n,t , which is the purchasing quantity for product n in period t. Excess stock is transferred to the next period as inventory, and unmet demand is lost without paying penalty costs. The end-ofperiod inventory for product n in period t is I n,t . The delivery lead time of products is assumed to be zero. The inventory flow equation for product n for two consecutive periods is given below, where (x) + denotes max{x, 0}. Assume the time lag length of receipt for product n is L. Define R n,t as the realized revenue from product n in period t. R n,t can be expressed by Eq. (2), where L is the receipt delay length and min{I n,t−L−1 + Q n,t−L , d n,t−L } is the realized demand for item n in period t. Based on Eq. (1), R n,t can also be written as p n (I n, (2) means that the sales income from period t will be received in period t + L. Assume that there are no unreceived payments at the beginning of the planning horizon. The cash flow equation is: In Eq. (3), when t > 0, one period's end-of-period cash balance equals its initial cash in this period (C t−1 ) plus this period's total revenues ( N n=1 R n,t ), minus this period's total ordering costs ( N n=1 v n Q n,t ) and overhead costs (H t ). At the end of planning horizon T , there are still some unreceived payments. For ease of analysis, we assume that the retailer discounts those payments to the final period T with discounting factor α. According to the above description, the final cash FC at the end of period T after discounting unreceived payments is shown by Eq. (4), where C T is the cash balance at period T before discounting the unreceived payments. The total amount of unreceived payments is N n=1 L k=1 R n,T +k . The retailer faces cash constraints in each period, which is common in real transactions for such small retailers. The cash constraint means that the retailer's maximum purchasing quantity is bounded by its available cash, which can be represented by the following formula: 4 Mathematical models In this section, several mathematical models are formulated to address different financing solutions provided by ecommerce platforms. We first provide the deterministic models and then build the corresponding stochastic models. In the deterministic situation, the uncertain demand for product n in period t is predicted to be the mean value of its demand distribution, which is represented as d n,t . We build two models below that are distinguished by whether order-based loan is used. For many retailers, if operating without supply chain financing, the cash inflow is mainly initially owned cash and revenues from sales. The whole linear deterministic model for the situation of operating with self-owned cash (2), (3), (5) The objective function (4) is to maximize the expected final cash after discounting the unreceived payments. Constraint (2) is the revenue equation, while Constraint (3) is the cash flow equation. Constraint (5) is the cash constraint. Constraints (6)-(10) are the linear expressions for the inventory flow equation (1). Constraint (11) shows the nonnegativity/binarity of end-of-period inventory, ordering quantity and δ t . The nonnegativity of I t is the result of the lost sales assumption. In the presence of order-based loan, when the retailer's application for order-based loan is approved, the platform sets a credit limit for the retailer. Assume that the limit of order-based loans is B U . The retailer can use this kind of loan repeatedly for many orders over a time limit as long as the total payment subject to order-based loans is less The interest rate of the order-based loan is R o . A new decision variable g n,t represents the quantity of order-based loans employed for product n in period t. For ease of analysis, assume that the retailer repays an order-based loan after a receipt delay of length L, and the time limit for the order-based loan is the planning horizon length T . The deterministic model for order-based loans (OL-D) is given below: (3), (5) − (11) g n,t ≤ I n,t−1 + Q n,t − I n,t The differences between Model OL-D and Model SO-D are the constraints related to g n,t . Constraint (12) means that the order-based loan employed for a product should be less than its realized demand (I n,t−1 + Q n,t − I n,t ) in this period. Constraint (13) gives the expression for revenue in the presence of order-based loans: when t ≤ L, the retailer's revenue from product n is p n g n,t since it can use order-based loan g n,t for this product; when t > L, I n,t−L−1 + Q n,t−L − I n,t−L − g n,t−L is the quantity of demand that does not use order-based loan in period t − L, and this portion of revenue comes at period t; and p n g n,t−L (1 + R o ) L ) is the total principle and interest that the retailer needs to repay in period t for period t − L's order-based loans; when t > T , the retailer cannot use order-based loan. Note that because of receipt delaying, the retailer receives some revenue when t > T . Constraint (14) suggests that the total employed order-based loan over the planning horizon should be no more than the order-based loan limit. Constraint (15) ensures the nonnegativity of g n,t . In real transactions, customer demand is uncertain, which is represented by demand scenarios in this paper. The demand for product n in period t under scenario s is represented by d s n,t , and the probability of scenario s occurring is Pr s . The set of all scenarios is S . Each scenario consists of the discrete values of all demands over the whole planning horizon. In the scenario programming model, the decision variables are Q s n,t , R s n,t , I s n,t , g s n,t , C s n,t , FC s , and δ s n,t . The stochastic model when the retailer is operating with self-owned cash (SO-S) is formulated below: I s n,t ≥ 0, Q s n,t ≥ 0, δ s n,t ∈ {0, 1} The objective function (16) is to maximize the expected final cash for all scenarios after discounting the unreceived payments. Constraints (17)- (24) and (28) represents the sets of scenarios that share the same history as scenario s before period t. Non-anticipativity means that all the scenarios with the same history before a given period should result in the same values of decision variables until that period. There are no non-anticipativity constraints for R s n,t and C s n,t because they are basically the expressions of I s n,t , Q s n,t and δ s n,t . In a similar way, the stochastic model when the retailer is operating with order-based loan (OL-S) is formulated below: Pr s ′ g s n,t = g s n,t s ′ ∈J(s,t) Pr s ′ (32) g s n,t ≥ 0 (33) In this section, we describe the data used in the numerical cases and the techniques employed for scenario generation and reduction in the study. A case study is derived from an online store that sells computer peripherals such as keyboards, mice, headsets, etc. on the Chinese e-commerce platforms Tmall.com and Jd.com. Since it is difficult to obtain internal sales data from retailers, we crawl customer comments from the retailer' online stores. Assume that 10% of customers leave comments after buying the products and that there is an approximately one week delay between the date they order and the date they comment, so we can estimate the associated demand. Note that although this seems to be a rought estimate for customers' demand, the scenario generation and reduction methods adopted in this paper can be applied to more exact demand values (e.g. Feng and Ryan, 2013; Hu and Hu, 2018) . For ease of analysis, two weeks is deemed a period in the planning horizon. Because demand typically booms in some months, such as the Single's Day shopping event in November and the semester-opening promotion events in March and September, we classify the planning horizon into two situations: booming-demand periods and normaldemand periods. We fit distribution parameters using a maximum likelihood approach. The empirical study revealed that a log-normal distribution generally provides a good fit for our demand samples. The resulting fitted values of the three main products (keyboard, mouse, headset) are shown in Table 2 ; the respective p-values from Kolmogorov-Smirnov tests are also reported in the table. Since all the p-values are above 0.05, we hold that the demand fits a log-normal distribution at the 95% confidence level. Because the length of order-based loans provided by the e-commerce platform usually is at most 3 months (6 two-week periods), the length of the planning horizon in our problem is 6 periods. A limited number of scenarios can be generated to satisfy the specified statistical properties of stochastic distributions. This paper adopts the moment matching method to generate scenarios. The main idea behind this method is to minimize the distances between the specified statistical properties of the generated scenarios (Høyland and Wallace, 2001) . We briefly summarize the moment matching technique proposed by Høyland and Wallace (2001) . Recall that the probability of scenario s is Pr s and that the set of all scenarios is S . Let ψ represent the specified stochastic property of the demand distribution and Ψ be the set of all specified properties. Let V ψ be the value for the stochastic property ψ and f (s, Pr s ) be the mathematical expression for V ψ . In the moment matching method, the decision variables are scenario realizations s and the corresponding possibility Pr s . For example, to generate scenarios for stochastic demand with a normal distribution, the stochastic property ψ can be the mean or variance. V ψ is the given The objective function (34) The symbol "∼" denotes "close to"; D indicates the dimension of each scenario node vector, and y represents the number of branches from each node. For example, with regard to our three-product, six-period problem, D is 18 for a problem with 3 products and 6 periods. By specifying three moments, the mean, variance and skewness, the number of specifications is 18 × 3 = 54. Therefore, by the above formula, we choose to create 3 realizations in each period because (18 + 1) × 3 − 1 ≥ 54. The fmincon function in MATLAB is adopted to solve this nonlinear optimization model (34)-(36). Multiple starting points are tested, and we select the solution with the smallest objective value sufficiently close to zero. To examine how uncertainty can affect decision making, 3 scenario trees are constructed. The size of each scenario tree is 3 6 . A summary of the demand scenarios for the first period in each scenario tree is shown in Table 3 . Recall that since demand is independent in each period, demand realizations and corresponding possibilities are the same for each period. For example, in the first period of scenario tree 1, the probability of demand realizations of 133, 246, and 87 for the three products is 0.1. In the latter periods, the probability of the same realizations is still 0.1. In each scenario tree, there are 3 6 scenarios. With a larger planning horizon or more products, the problem may become computationally intractable. For example, if using weekly demand, the length of the planning horizon will be 12, and there will be 3 12 scenarios. Therefore, we also make use of the scenario reduction method pro- posed by Heitsch and Römisch (2003) to reduce the computational complexity. With regard to the two reduction methods: fast forward selection method (FFS) and backward selection method (BS), we make use of FFS because it outperforms BS when the number of selected scenarios is no more than 25% of the number of original scenarios (Heitsch and Römisch, 2003; Hu and Hu, 2018) . The main idea behind FFS is described as follows: in the first step, select the scenario with the smallest total weighted distance; then, update the distances and select the scenario among the remaining unselected scenarios; and after selecting enough scenarios as required, add the probabilities of unselected scenarios to the probability of their closest selected scenarios. The pseudocode of FFS is given by Algorithm 1 in the appendix. In this section, we first discuss the stability test results of the scenario generation and reduction method, then make comparisons between deterministic models and stochastic models, and finally investigate the impacts of different initial cash, overhead cost, receipt delay length, and demand fluctuations on the retailer's choice of using order-based loans. Selling prices and unit purchasing costs for the three products are listed in Table 4 . Assume that the initial cash balance of the retailer is 20,000. The discount rate α is 1%, and the interest rate r 0 for order-based loan is 1.5%. The selling prices and unit purchasing costs of the three products are listed in Table 4 . Other parameter values are as follows: receipt delay length L is 2; order-loan limit B u is 10,000; overhead cost H t is 2000 in each period; and the demand patterns are [0, 0, 0, 0, 1, 1], where 0 means normal demand in this period and 1 means booming demand. The scenario sample sizes are generated from 20 to 180 in steps of 20. As Figure 2 shows, there is a slightly decreasing trend in the objective value as the scenario sample size increases. When the scenario size is 120 or 140, the objective value tends to be stable for all the scenario trees in the two situations. Therefore, we decide to select 140 as the scenario sample size. In addition, the objective values have very small fluctuations across different sample sizes. Next, we analyze the in-sample stability and out-of-sample stability of the scenario tree. In-sample stability means that the objective values should not vary substantially across different scenario trees of the same size, while out-of-sample stability indicates that the solutions from one scenario tree should exhibit close objective values in another scenario tree. Table 5 lists the details of the stability test for the two scenario models. For instance, 46092 in row 1 and column 2 means that substituing the solution from Tree 1 into Tree 2 gets the value 46092. Tree 3 46117 46126 46131 46247 46258 46263 As Table 5 shows, since the gaps between different solutions in different scenario trees are all less than 5%, we claim that the stability test for these scenario generation results is valid. By fixing other parameter values and allowing one parameter to vary, we attempt to determine whether it is worth considering demand uncertainty in different parameter settings. The parameter variations are given by Table 6 . Receipt delay length L 0 1 2 3 4 Overhead Cost H t 0 1,000 2,000 3,000 4,000 Let In Figure 3 , VSS grows with increasing initial cash but decreases with increasing overhead cost and receipt delay length. Since a larger overhead cost or receipt delay length would make the retailer more cash-constrained during operation, this means that the stochastic model has a greater advantage when the retailer is less cash-constrained. With regard to EVPI, it grows as the initial cash increases and decreases as the overhead cost increases. However, as the receipt delay length increases, the EVPI does not show a clear monotonous pattern. Moreover, both the VSS and EVPI values when using order-based loans are larger than those when not using loans. Therefore, it is much better to use the stochastic model when using order-based loans. We also want to know which factors affect the retailer's choice of using order-based loans, namely, in which parameter settings the order-based loan model (Model OL-S) achieves higher expected profit than the model of not using loan (SO-S). This is shown by Figure 4 , where the y axis is the expected profit gap between using order-based loan and not using loan. We also compare several demand patterns in Figure 4( The above figures show that the retailer would prefer using order-based loans when the receipt delay length is long or its initial cash is small. When the overhead cost is large, using order-based loans tends to yield more profits than not using them. However, when the overhead cost is too large, there is no advantage in using order-based loans because the retailer cannot make profits. With regard to demand patterns, it seems that it is better for the retailer to use order-based loans when there expected demand is high. Many small online retailers encounter cash shortages during operation, and some Chinese e-commerce platforms have provided financing services to help retailers alleviate this problem. In this paper, we discuss order-based loans by building a multi-period multi-product model with demand uncertainty. Scenario programming and reduction techniques are used to solve this problem. Numerical examples show that the retailer tends to use order-based loans when the receipt delay length is long or initial cash is small. When the overhead cost in each period is large but not too large to make the retailer unprofitable, the retailer would also prefer using order-based loans. Future research may consider several directions: one possible direction is to consider demand substitution in which one product can be substituted by others if it is out of stock; a second possible direction is to investigate other supply chain financing behaviors in a multi-period stochastic problem. Scenario-based planning for lotsizing and scheduling with uncertain processing times Dynamic mixed-item inventory control with limited capital and short-term financing A cash-constrained stochastic inventory model with consumer loans and supplier credits: the case of nanostores in emerging markets Multi-item capacitated lot-sizing with demand uncertainty Inventory management with asset-based financing A financial crisis is looming for smaller suppliers Dynamic inventory management with cash flow constraints A capital flow-constrained lot-sizing problem with trade credit A dynamic ordering policy for a stochastic inventory problem with cash constraints A cash-constrained dynamic lot-sizing problem with loss of goodwill and creditbased loan. International Transactions in Operational Research Supply chain management: more than a new name for logistics. The international journal of logistics management Financing newsvendor inventory Technology is revolutionising supply-chain finance. The Economist A multi-stage stochastic program for supply chain network redesign problem with price-dependent uncertain demands Scenario construction and reduction applied to stochastic power generation expansion planning Dynamic inventory control with limited capital and short-term financing Progressive hedging as a meta-heuristic applied to stochastic lot-sizing Scenario reduction algorithms in stochastic programming Dynamic capacitated lot sizing with random demand and dynamic safety stocks Robust single-item lot-sizing problems with discrete-scenario lead time Generating scenario trees for multistage decision problems A two-stage stochastic programming model for lot-sizing and scheduling under uncertainty A multi-stage stochastic programming for lot-sizing and scheduling under demand uncertainty Non-collaborative and collaborative financing in a bilateral supply chain with capital constraints Cash-flow based dynamic inventory management. Production & Operations Management Financing the newsvendor: supplier vs. bank, and the structure of optimal trade credit contracts Two-stage inventory management with financing under demand updates A joint model for cash and inventory management for a retailer under delay in payments Short-term financing in a cash-constrained supply chain Buyer intermediation in supplier finance A trade credit model with asymmetric competing retailers. Production and Operations Management Partial credit guarantee and trade credit in an emission-dependent supply chain with capital constraint Supply chain with random yield and financing Guarantor financing in a four-party supply chain game with leadership influence The research presented in this paper is supported by Chongqing Social Science Planning Fund under Project Appendix A FFS algorithm for scenario reduction s i and s j represent scenario i and scenario j (i, j ∈ S ). d [k] i, j is the distance of s i and s j in step k and WD [k] i is the overall weighted distance of scenario i in step k. The target number of selected scenarios is K. In step k, the set of selected scenarios is represented by Ω [k] while the set of unselected scenarios is represented by J [k] . L(i) represents the set of scenarios whose closest scenario is scenario i (i ∈ Ω [K] ) in the final step. Data: the set of all scenarios S , target number of selected scenarios K.Result: the selected K scenarios. , ∀ j ∈ J [K] , where L(i) = {s j | arg min i∈Ω [K] , j∈J [K] {d [K] i, j } = i}.