1.INTRODUCTION
Supply chain is defined as an integrated process where different business entities such as suppliers, manufacturers, distributors, and retailers work together to deliver products to customers. This chain is concerned with two distinct flows: a forward flow of materials and a backward flow of information (Karakitsiou and Migdalas, 2008). The competition forces companies to look at their operations from a supplychain perspective and to seek improvements from better coordination and communication across the supply chain.
In this paper, we consider the problem of selecting the appropriate distribution strategy for delivering a family of products from a set of suppliers to a set of retailers so that the total transportation and inventory costs are minimized. There were some research works about this issue. Harris (1990) and Burns (1985) have presented a model to deal with these problems in an integrated manner. Viswanathan (1997) have considered a twostage model for multiproduct capacity planning of supply chain networks. Ye and Du (2003) formulated the minimization of total cost as the objective function, and established the optimal model of the supplier’s maximum inventory and distribution periods. However, these models do not take any uncertain factors into account.
Bertazzi and Cherubini (2013) researched an inventory transportation system with stochastic demand. They studied a logistic system in which a supplier has to deliver a set of products to a set of retailers to face a random demand over a given time horizon. Agrawal and Smith (2013) optimized inventory management for a retail chain with diverse store demands. Accounting for these variations appropriately in inventory management can significantly improve retailers’ profits. Pyke (1993) have developed a stochastic programming model that can be used to determine the optimal replenishment policies as well as the optimal dispatch strategy simultaneously. Shu (2005) dealt with the problem of minimizing the sum of transportation and inventory costs for shipping several products on a common link where the demand is a random variable. Zhang (2005) presented a model to design simple inventory policies and transportation strategies to satisfy uncertainty demands over a finite horizon, while minimizing system wide cost by taking advantage of quantity discount in the transportation cost structures. Hu et al. (2014) studied on the influence of supply chain joint decision from the perspective of supply chain finance. Based on the research for Xi’an equipment manufacturing and biomedicine industry enterprises, structural equation model is established. We know that, by improving the use frequency of supply chain finance and enhancing trust, communication, cooperation, small and mediumsized enterprises have a significant impact. Fu (2007) discussed an analytical model for coordinating inventory and transportation decisions in a vendormanaged inventory system and constructed the adaptive genetic algorithm to solve the model. Liu (2009) presented a model to consolidate inbound freight at transshipment point from multiple sources to a single destination to minimize the overall transportation and inventory costs, which is subject to demand uncertainty, at the same time, heuristic algorithm also was used to solve the two subproblems of inventory and transportation. As the new research area in supply chain management, financial supply chain focuses that cash flow should correspond with information and logistic. Such as Liu (2016) and Lambert (2016), through communication, information resource sharing, efficient and low cost supply chain formation, enterprises can better cope with the fierce market competition.
From the latest literatures, we can know that the research of supply chain still has certain significance and value. However, these references only consider single period constraint model, didn’t put the periodic into account. Similarly, we also focus on minimizing the transportation and inventory costs in a supply chain where the demand is a random variable. But, further, we extend the model to the multiperiod with random variables, and then, translate the stochastic programming (Panos et al., 2015) into a deterministic model. Based on the genetic algorithm and stochastic simulation, a hybrid intelligent algorithm (HIA) having been explored can help us to find a good system solution.
The rest of this paper is organized as follows. Section 2 presents a singleperiod and a multiperiod stochastic programming model where the demands are random variables. Section 3 presents HIA to handle the multi period stochastic programming model. Numerical examples and computational experiments have been reported in Section 4. Section 5 concludes the paper with a summary and direction for future work.
2.MODEL ESTABLISH
In this section, we will give the description of transportation inventory strategy and how to construct our model.
2.1.SinglePeriod Model
In this section we present the singleperiod model that serves as the building block for modeling a multi period stochastic programming. The model will do the following assumptions:

1. Only consider the same kind of products;

2. The random demand for each product at each retailer form each supplier is any positive number;

3. The transportation and inventory costs are only determined by the quantity, regardless of the weight;

4. The customer demand is a random variable, and the demands of different retailers are independent of each other;

5. The retailers’ inventory supplement take single period (T,S) replenishment policy (According to the predetermined interval T , each retailer need to check the stock, put forward the order, and then, replenish the target inventory stock S );

6. The delivery from each supplier to each retailer can be completed immediately;

7. The transportation costs includes order costs.
Secondly, we introduce the following notation to present the mathematical formulation throughout the remainder of this paper.
Indices, Sets, and Parameters.

I : set of suppliers, I = {1, 2, ..., m} .

J : set of retailers, J = {1, 2, ..., n} .

i : index for set of suppliers I.

j : index for set of retailers J.

a_{i} : output of supplier i.

c_{ij} : transportation cost of shipping one unit of product from supplier i to retailer j.

c_{j} : inventory cost of one unit of product for retailer j.

x_{ij} : decision variable, transport volume of products from supplier i to retailer j.

y_{j} : random variable, demand of products of retailer j.

μ_{j} : the expected demand volume of retailer j.

σ_{j} : the demand standard deviation of retailer j.

s_{j} : retailer j inventory volume of products at the end of this period.

s_{j0} : retailer j inventory volume of products at the beginning of this period.

α_{j} : confidence level of retailer j.

Notice: Amount of product units: ten thousand.

The cost unit: $
we consider the following objective function as a chance constraint of our model. In addition, the objective function is to minimize the inventory and transportation costs. Then, we take available resources into account as other constraints to construct the stochastic programming for transportation inventory model as follows:
subject to
In the model, (1) is the objective function of minimizing the inventorytransportation costs. Constraints (2) is the significance that the transport volume of the supplier i to all the number of retailers can not exceed the total output of this supplier. Constraints (3) is a chance constraint, which indicates that the total storage capacity of retailer j at the periodic initial state is greater than or equal to the probability of the retailer’s current demand that is not less than the confidence level of retailer j. Constraints (4) meets the condition of decision variable is non negative.
Obviously, our model above is a stochastic chanceconstrained programming model (SCCPM). We assume that customer demands is subject to normal distribution, the expected demand volume of retailer j is μ_{j} , and the demand standard deviation of retailer j is σ_{j} . Density function and distribution function are, respectively, f(y_{j}) and F(y_{j}). After the equivalent transformation, the above model will be equivalent to the following deterministic model:(5)
This is a deterministic nonlinear programming model, and we can call directly the MATLAB optimizer to solve.
2.2.MultiPeriod Model
In real strategy decisionmaking process, the multi period situation is more common and in this situation the decisionmaker can revise his production plan when more information about the customer demand becomes available. Many latest papers have addressed such situation, for example, Masoumeh et al. (2009) proposed a multiproduct, multiperiod planning problem under quality of raw materials. And then, we will construct the multi period model which is an extension of the singeperiod stochastic programming in Section 2.1. At first, we add the following assumptions:

1. The customer demand is a random variable, and the demands of different retailers and the different periods of the same retailers are independent of each other.

2. The retailers’ inventory supplement take multiple period replenishment policy and each replenishment period must satisfy the constraint conditions.

3. It is assumed that the moment of replenishment in each period is the periodic initial state.
Secondly, we introduce the following notation to present the mathematical formulation throughout the remainder of this paper.
Indices, Sets, and Parameters.

I : set of suppliers, I = {1, 2, ..., m} .

J : set of retailers, J = {1, 2, ..., n} .

T : set of periods, T = {1, 2, ..., p} .

i : index for set of suppliers I .

j : index for set of retailers J .

t : index for set of periods T .

a_{it} : output of supplier i in period t.

c_{ijt} : transportation cost of shipping one unit of product from supplier i to retailer j in period t.

c_{jt} : inventory cost of one unit of product for retailer j in period t.

x_{ijt} : decision variable, transport volume of products from supplier i to retailer j in period t.

y_{jt} : random variable, demand of products of retailer j in period t.

μ_{jt} : the expected demand volume of retailer j in period t.

σ_{jt} : the demand standard deviation of retailer j in period t.

s_{jt} : retailer j inventory volume of products at the end of the period t.

s_{j0} : retailer j inventory volume of products at the beginning of the period.

α_{jt} : confidence level of retailer j in period t.

Notice: Amount of product units: ten thousand.

The cost unit: $
Suppose that it is a multiperiod delivering system from m suppliers to n retailers. Output of supplier i in period t is a_{it}. Suppliers delivery products at the periodic initial state for retailers, and each period’s total transport volume is $\sum _{i\in I}{x}_{ijt}$ from all suppliers to the retailer j . Set up the first period of the retailer j ’s inventory is s_{j0}. Obviously, the inventory at the end of the last period is the initial inventory of the next period. We consider the following objective function as a chance constraint of our model. In addition, the objective function is to minimize the cost of transportationinventory. Then, we take available resources into account as other constraints to construct the model as follows:
subject to
The objective function (9) represents the sum of the total cost of transportation and inventory in all periods, and it needs to be minimized. Constraints (10) are expressed the number goods from supplier i to all retailers in period t can not exceed the total output of the supplier i. Constraints (11) is chance constrained conditions, which indicates that the total storage capacity ($T{S}_{jt}={s}_{j,\text{\hspace{0.17em}}t1}+{\displaystyle \sum _{i\in I}{x}_{ijt}}$) of retailers j , at the periodic initial state of period t , is greater than or equal to the probability of the demand of the retailer that’s not less than the confidence level α_{jt} . Constraint condition (12) is a non negative condition of decision variable. (13) indicates that the stocks at the beginning and end of the whole period is 0.
It is a multiperiod dynamic programming model. Hossein (2015) and Nan (2014)have studied the similar problems. In the case of customer demand is a random variable, the number of replenishment for any retailer is not only affected by the amount of preinventory, but also depends on the random demand in this period. In view of the characteristics of multiperiod model, we can not translate the stochastic programming problem into a deterministic model. We may not solve it by classical optimization methods or intelligent algorithms. The next work is to explore the HIA  based on genetic algorithm and stochastic simulation, so as to handle the model.
3.SOLUTION METHOD
The traditional method of solving the chance constrained programming is to convert the chance constrains into their determinate equivalent forms, and then to solve the determinate model (such as singleperiod model). However, this traditional method isn’t applied to multi period model above. Because the multiperiod model is essentially a dynamic programming model and the variables of different period influence each other. So, in this paper, we explore HIA (Jose et al., 2013) which was mentioned by Liu (2003). This algorithm can solve the general chanceconstrained programming. The main thought of this algorithm is to convert the stochastic model to the certain one. First, we adopt stochastic simulation technologies to produce the demands as the input data, and then calculate the output data according to a formulation we have designed. These inputs and outputs are used as the training samples for the back propagation neural network (BP neural network). Second, we train a BP neural network as the fitting function in the next stage. Finally, we use the improvedgenetic algorithm to find the best solution. The flowchart for HIA is showed in the following Figure 1:
3.1.Stochastic Simulation Technologies
The SCCPM is different from the traditional model which contains some random variables. The probabilities of random events are taken as constrains. Thus, these probabilities should be obtained first. We use stochastic simulation method to compute them.
3.1.1.The Principle of Monte Carlo Technology Approaching Probability
We consider $P\left\{f\left(x,\text{\hspace{0.17em}}\xi \right)\le b\right\}$, where f (x, ξ) is a function of stochastic vector. First of all, N independent random vectors ${\xi}_{n},\text{\hspace{0.17em}}\left(m=1\text{\hspace{0.17em}},\text{\hspace{0.17em}}2,\text{\hspace{0.17em}}\cdots ,\text{\hspace{0.17em}}N\right)$ are generated from the probability distribution φ (ξ) of random vector ξ.
Denote the number of f_{n} which satisfy the condition f (x, ξ) ≤ b , as M. According to the Law of Large Numbers, when N is sufficiently large, $P=\frac{M}{N}$ is used as an approximation of $P\left\{f\left(x,\text{\hspace{0.17em}}\xi \right)\le b\right\}$.
We know that the maximum y_{jt} which satisfies the inequality:
$P\left\{{s}_{j,\text{\hspace{0.17em}}t1}+{\displaystyle \sum _{i\in I}{x}_{ijt}\ge {y}_{jt}}\right\}\ge {\alpha}_{jt}$
should be achieved at the equal sign, that is
$P\left\{{s}_{j,\text{\hspace{0.17em}}t1}+{\displaystyle \sum _{i\in I}{x}_{ijt}\ge {y}_{jt}}\right\}={\alpha}_{jt}$
Here, we denote
$f\left(\xi \right)={s}_{j,\text{\hspace{0.17em}}t1}+{\displaystyle \sum _{i\in I}{x}_{ijt}}$, the sequence
$f\left({\xi}^{1}\right),\text{\hspace{0.17em}}f\left({\xi}^{2}\right),\text{\hspace{0.17em}}\cdots ,\text{\hspace{0.17em}}f\left({\xi}^{N}\right)$
is produced after N simulations.
Define $h\left({\xi}^{k}\right),$
When $\begin{array}{l}f\left({\xi}^{k}\right)\ge {y}_{jt},\text{\hspace{0.17em}}\left(k=1,\text{\hspace{0.17em}}2,\text{\hspace{0.17em}}\cdots ,\text{\hspace{0.17em}}N\right)\\ h\left({\xi}^{k}\right)=1;\end{array}$
otherwise, $h\left({\xi}^{k}\right)=0$
According to the law of large number, when N tends to infinity the value of
$\frac{{\displaystyle \sum _{k\in N}h\left({\xi}^{k}\right)}}{N}$
tend to α_{jt} which presents the probability of f (ξ) ≥ y_{jt}. On this basis, when the sequence
$f\left({\xi}^{1}\right),\text{\hspace{0.17em}}f\left({\xi}^{2}\right),\text{\hspace{0.17em}}\cdots ,\text{\hspace{0.17em}}f\left({\xi}^{N}\right)$
is sorted in ascending order, the N^{’th} element in the arrangement gives the required value, where N ’is the integer part of ${\alpha}_{jt}\times N$.
3.1.2.Monte Carlo Technology to Produce the Simulation Data
Now, note α_{jt} = 0.95 and define the following function U(x) .
Where
$U\left(x\right)=\text{min{}{y}_{jt}P({s}_{j,\text{\hspace{0.17em}}t1}+{\displaystyle \sum _{i\in I}{x}_{ijt}\ge {y}_{jt}})\ge 0.95\text{}}$
The steps for simulating U(x) are as follows:

Step 1. Initialize the parameters: number of simulation N , confidential level α_{jt} = 0.95.

Step 2. Generate (i, j, t) randomly, and then, determine x_{ijt}.

Step 3. Generate ξ, randomly. Note that the value generated by the random function should be treated, so as to ensure nonnegative.

Step 4. Calculate the value of ${s}_{j,\text{\hspace{0.17em}}t1}+{\displaystyle \sum _{i\in I}{x}_{ijt}}$.

Step 5. Repeat steps 24 N times and then get N values of ${s}_{j,\text{\hspace{0.17em}}t1}+{\displaystyle \sum _{i\in I}{x}_{ijt}}$. Sort these values in ascending order, the N ’th is the maximum y_{jt} needed.
Thus, give an x , there is a U(x) (the maximum y_{jt} ) through N times’ simulation. Repeat M times like this, we obtain M training samples. By these samples, train a neural network which is used to approximate the uncertain function U.
3.2.Training for Neural Network
Train a neural network (the number of input neuron is m, that of hidden layer neurons is y, and the number of output neurons is 2) by the input and output data to approximate the uncertain function. The uncertain function is used as the evaluation function in the next layer algorithms. Where y is determined by the formula,
$y<\sqrt{p+q}+a$
Where p , q represent the number of input neurons and output neurons respectively, a is a constant, 1 < a < 10 (Li Xiaoyan, 2009).
3.3.Genetic Algorithm to Solve the Model
The general steps of genetic algorithm:

Step 1. Initialize the chromosomes, compute the fitness of chromosomes by the trained neural network.

Step 2. Cross, mutate, select, and compute the fitness of the offspring chromosomes by the trained neural network.

Step 3. Pick up chromosomes by roulette wheel selection.

Step 4. Repeat steps until the completion of the given number of cycles.

Step 5. Obtain the best chromosome, adopt this chromosome as the optimal solution.
In the following the given example, assuming that algorithm parameters as follows:
Population size is 30;
Using twopoint crossover and crossover probability P_{c} = 0.3;
Mutation probability P_{m} = 0.2;
The parameter in the evaluation function based on the sequence is α = 0.05 .
4.ILLUSTRATIVE EXAMPLE
In order to test the model of transportation and inventory joint optimization and the solution algorithms applied in the actual situation, we perform numerical tests based on the domestic operation department with reasonable assumptions.
4.1.Illustrative Example of SinglePeriod Model
We consider two suppliers, I = {1, 2} , noting A_{1}, A_{2}, And three retailers, J = {1, 2, 3} , notes B_{1}, B_{2}, B_{3}. This is a two level supply chain system that needs deal with the singleperiod transportation and inventory optimization problem under the situation of demand is a random variable. Confidence level is α_{j} = 0.95, so, Φ^{1} (α_{j}) = 1.65. c_{ij} is the transportation cost of shipping one unit of product from supplier i to retailer j; c_{j} is the inventory cost of one unit of product for retailer j; a_{i} is the total output of supplier i. The random demand of each retailer is subject to the normal distribution $N\left({\mu}_{j},\text{\hspace{0.17em}}{\sigma}_{jt}{}^{2}\right)$,
and the density function is $f\left({y}_{j}\right)=\frac{\text{exp}\frac{{\left({y}_{j}{\mu}_{j}\right)}^{2}}{2{\sigma}_{j}^{2}}}{\sqrt{2\pi}{\sigma}_{j}}$. The original schedule is Table 1 and Table 2 as follows: (14)
As shown above, the singleperiod model is obtained according to the known data. Calling directly MATLAB optimizer to solve, then, we can get the optimal solution making the total cost of transportation inventory joint optimization model least (Minimum value is $17679) that has been showed in the following Table 3:
4.2.Numerical Experiment of MultiPeriod Model
In order to verify the correctness of the model and demonstrate how to use the proposed algorithm to solve the multiperiod model, we assume that the product P will be delivered from two suppliers A_{1}, A_{2} to three retailers B_{1}, B_{2}, B_{3}, and the period T which is made up of three periods t_{1}, t_{2}, t_{3}. Table 4 and Table 5 give the value of the parameters and others.
The confidence levels are α_{j} = 0.95, ( j =1, 2, 3, t = 1, 2, 3 ), and ${\Phi}_{jt}^{1}$(0.95) = 1.65. So, we have the following model about this problem:(18)
According to the specific algorithm parameters and process in section 3.3, we can make use of MATLAB to program and calculate the procedure. The results can be seen in Table 6. The minimum cost is $48483.
We consider a small problem consisting of two suppliers and three retailers whose computational time is only a few seconds. Even if we consider a largescale case, computation time will not be changed. At the same time, due to space limitations, we don’t need to study largescale problems.
5.CONCLUSIONS
This paper considers the singleperiod and multi period twolevel supply network and establishes the stochastic chanceconstrained programming models. In these models, we consider the whole costs of transportation and inventory as the objective function, and the probabilities of random events are taken as constrains. In order to solve models above, we propose HIA based on genetic algorithm and stochastic simulation. In the end of this paper, numerical examples was given and showed the effectiveness of the method.
A major contribution of this paper is that we provide a comprehensive framework for supply chain under uncertain conditions. Much work still needs to be done to improve on the current framework. One obvious direction is that reducing assumptions and restriction conditions, and studying the transportation of multiproducts. Thus, it will be more comprehensive to deal with the issue of supply chain.