<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Archiving and Interchange DTD v1.0 20120330//EN" "JATS-archivearticle1.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink">
  <front>
    <journal-meta />
    <article-meta>
      <title-group>
        <article-title>Anna Zykina, Olga Kaneva</article-title>
      </title-group>
      <contrib-group>
        <aff id="aff0">
          <label>0</label>
          <institution>Omsk State Technical University</institution>
          ,
          <addr-line>Omsk</addr-line>
          ,
          <country country="RU">Russia</country>
        </aff>
      </contrib-group>
      <fpage>251</fpage>
      <lpage>255</lpage>
      <abstract>
        <p>The paper deals with the transportation logistics problems in the conditions of incomplete information. The research includes several formulations of the stochastic optimization problem for different variants of the relationship "Resources reserves - Resources consumption". For the solvability of these problems we propose two-stage scheme for solving stochastic transportation problem. The novelty of the two-stage problem stochastic programming formulation has been contained in the statement of the second stage problem. Choosing of compensation plan is determined from the solution of the linear complementarity problem. Transportation problem; two-stage stochastic problems; linear complementarity problem.</p>
      </abstract>
    </article-meta>
  </front>
  <body>
    <sec id="sec-1">
      <title>Introduction</title>
      <p>another word stochastic problem with compensation of residuals [7]. The process of solving the problem
can be divided into two stages: on the first stage we select the preliminary plan from deterministic
conditions, on the second stage we implement compensation discrepancies, that have been identified after
the implementation of random events. This approach can be used for stochastic problems where a
preliminary decision should be taken and put into the implementation before we have know the value of
random parameters. For the transportation problem with random demand the preliminary decision can be
determined by the distribution of materials supplies taking into account the determined reserves of the
materials.</p>
      <p>In general the difficulties with the analysis of the two-stage stochastic transport problems are
determined by the need to choose the best preliminary plan of the original problem, which would guarantee
the existence of residual compensation for all implementations of parameters of uncertainty.</p>
      <p>In this paper a two-stage stochastic transport problem is considered, where the choice of
compensation plan subjects to the terms which are determined by a linear complementarity problem.
Statement of the Linear Complementarity Problem</p>
      <p>It is important to consider the linear complementarity problem in the form [8]:</p>
      <p>v  Bz  q, v j  0, z j  0, v j z j  0, j  1,..., p.</p>
      <p>There B – given a square matrix with size р, (v j , z j ) – is a couple of additional variables.
Condition vj z j  0, j  1,..., p, analogous to the condition of complementarity in the duality theory for
inequalities Bz  q  0 and z  0.This means that in a pair of conjugate inequalities at least one should
be implanted as equality.</p>
      <p>
        Non-negative definiteness of the matrix B ensures the solvability of problem (
        <xref ref-type="bibr" rid="ref1">1</xref>
        ). When the
positive definiteness of the matrix B there is a unique solution z of problem (
        <xref ref-type="bibr" rid="ref1">1</xref>
        ).
      </p>
      <p>
        For the solution building of problem (
        <xref ref-type="bibr" rid="ref1">1</xref>
        ) algorithm can be used as an additionalconversion Lemke
[8], it can be demonstrated as an analogue of the simplex algorithms for linear programming problem.
Mathematical Model of the Transport Problem
      </p>
      <p>It is necessary to consider the classical transport problem: minimize the total cost</p>
      <p>m n
for specified volumes of transportation of a homogeneous product X  {xij }, i  1,...,m, j  1,...,n, with
restrictions:
cij xij
i1 j1
n
 xij  ai , i  1,...,m,
j1
m
 xij  bj , j  1,...,n,
i1
xij  0, i  1,...,m, j  1,...,n.</p>
      <p>
        Conditions (
        <xref ref-type="bibr" rid="ref3">3</xref>
        ) determine the distribution of transportation X
ai , i 1,...,m.
      </p>
      <p>
        Conditions (
        <xref ref-type="bibr" rid="ref4">4</xref>
        ) ensure the fulfillment of the demand bj , j  1,...,n .
      </p>
      <p>It is necessary to introduce the vector representation of the transportation plan x and vector c
with the matrix representation of the transportation plan X and matrices C , by sticking rows of matrix to
vector
in accordance with reserves
xij  x((i1)n j) , cij  c((i1)n j).</p>
      <p>
        Then a group of equations (
        <xref ref-type="bibr" rid="ref3">3</xref>
        ) is denoted as Ax  a , where vector a  (a1 ,...,am ) , and matrix A
corresponds to constraints (
        <xref ref-type="bibr" rid="ref3">3</xref>
        ).
      </p>
      <p>
        Two-stage StochasticTransport Problem
(
        <xref ref-type="bibr" rid="ref1">1</xref>
        )
(
        <xref ref-type="bibr" rid="ref2">2</xref>
        )
(
        <xref ref-type="bibr" rid="ref3">3</xref>
        )
(
        <xref ref-type="bibr" rid="ref4">4</xref>
        )
(
        <xref ref-type="bibr" rid="ref5">5</xref>
        )
It is important to consider the stochastic formulation of the transportation problem(
        <xref ref-type="bibr" rid="ref2">2</xref>
        )–(
        <xref ref-type="bibr" rid="ref5">5</xref>
        ).
      </p>
      <p>Let the demand b  b( ) is a random variable. Also C  C( ) – the cost of transportation of
product is also random.</p>
      <p>
        ~
Define by bj , j  1,...,n
– some implication of a random variable
bj , using
~x , i  1,..., m, j  1,..., n, – a transportation plan that satisfies to the deterministic conditions (
        <xref ref-type="bibr" rid="ref3">3</xref>
        )
ij
and (
        <xref ref-type="bibr" rid="ref5">5</xref>
        ).
      </p>
      <p>Case 1. For some j  J  ,J   1,...,n following inequality holds:</p>
      <p>
        Then there are two cases of the choice of compensation plan from the terms and conditions
determined by the linear problem of complementarity (
        <xref ref-type="bibr" rid="ref1">1</xref>
        ).
      </p>
      <p>m ~
 ~xij  b j .</p>
      <p>
        i1
Condition (
        <xref ref-type="bibr" rid="ref6">6</xref>
        ) means that at the point j  J

      </p>
      <p>the demand is not satisfied.</p>
      <p>For such constraintsit is necessary to introduce a penalty s j , j  J


for one unit of
compensation plan of the product deficit and carry out compensation of obtained discrepancy as follows:</p>
      <p>E  ~x  b~   M  y .</p>
      <p>
Where E – is a matrix and b</p>
      <p>
        
correspond to set J , M  – is positive definite matrix and y
– a vector which were composed by the restrictions (
        <xref ref-type="bibr" rid="ref4">4</xref>
        ) which

– is a nonnegative vector of the
(
        <xref ref-type="bibr" rid="ref6">6</xref>
        )
(
        <xref ref-type="bibr" rid="ref7">7</xref>
        )
(
        <xref ref-type="bibr" rid="ref8">8</xref>
        )
(
        <xref ref-type="bibr" rid="ref9">9</xref>
        )
(11)
(12)
compensation plan of corresponding dimensions.
      </p>
      <p>
        Labeling
we will receive linear complementarity problem (
        <xref ref-type="bibr" rid="ref1">1</xref>
        ) for conditions (
        <xref ref-type="bibr" rid="ref7">7</xref>
        ).
      </p>
      <p>Case 2. For some j  J  , J   1,...,n the following inequality is satisfied
z  y , B  M  , q  E ~x  b~  ,
i1
m ~ ~
 x  b .</p>
      <p>
        ij j
The condition (
        <xref ref-type="bibr" rid="ref8">8</xref>
        ) means that at the point j  J
      </p>
      <p>
For such constraints we introduce a penalty s j , j  J

excess product and hold compensation of obtained discrepancy by following:</p>
      <p>E  ~x  b~   M  y .</p>
      <p>there is a need to store excess product.</p>
      <p></p>
      <p>for one unit of compensation plan

Where E – is a matrix and b</p>
      <p>
        
correspond to set J , M  – is positive definite matrix and y
– a vector which were composed by the restrictions (
        <xref ref-type="bibr" rid="ref4">4</xref>
        ) which

– is a nonnegative vector of the
compensation plan of corresponding dimensions.
      </p>
      <p>By analogy in case 1, when</p>
      <p>
        z  y  , B  M  , q  E  ~x  b~ ,
we will receive linear complementarity problem (
        <xref ref-type="bibr" rid="ref1">1</xref>
        ) for conditions (
        <xref ref-type="bibr" rid="ref9">9</xref>
        ).
      </p>
      <p>Transportation problem with random demand can be presented as the following two-stage
stochastic programming problem:</p>
      <p>
         nm   
M  k1 ck xk  ym ,iny  (s  y   s y )  mxin (
        <xref ref-type="bibr" rid="ref10">10</xref>
        )
~ 
~ 
v  M  y  E x  b~ , v  0, y   0, v y   0, j  J  .
      </p>
      <p>j j</p>
      <p>For the solvability of the second stage with all implementations of a random variable b and with
any preliminary plan x it is necessary and sufficient that the matrix M  and M  was positive definite.
Under such conditions, there is a unique solution y (x,b) and y (x,b) for each complementarity problem
where y (x,b),</p>
      <p>y (x,b) – solutions of the corresponding complementarity problems, the objective
space, and the feasible set D  x | Ax  a, x  0 is a bounded and convex linear set.
function is an expectation of a random function depending on a random vector from a certain probability</p>
      <p>For solving the problem (14) – (15) can use the methods of stochastic approximation [9, 10].</p>
    </sec>
    <sec id="sec-2">
      <title>Algorithm</title>
      <p>Monte-Carlo estimators of the objective function and d
projection of the gradient estimate to the  -feasible set).</p>
      <p>Step 2. Go to the next point</p>
      <p>Let the initial approximation of the solution x0  D and some initial Monte-Carlo sample size
N 0 be given. Put k  0 and move on to the main stage.</p>
      <p>Step 1. Let the vector x is known. We generate N k values of a random variable b andcalculated
k
k k
is an  -feasible direction at the point x
(i.e.,
xk 1   x ( xk  k d k ).</p>
      <p>N k 1 
k  d
Where C is a certain constant k   x k ( d k ) , d k
projection of the gradient estimate to the  -feasible set).</p>
    </sec>
    <sec id="sec-3">
      <title>Conclusion</title>
      <p>is an  -feasible direction at the point x
k
(i.e.,
of the problem for the positive semi definite matrix M
preliminary plan x .</p>
      <p>The novelty of the two-stage schemes stochastic transportation problem formulation contains in
the statement of the second stage problem. Two-stage scheme can be applied for solving the stochastic
transportation problem in following case. Conditions (15) give the distribution of the known deterministic
reserves. This is a preliminary plan. After the values of the random demand becomes known, we have a
possibility of the occurring residuals.</p>
      <p>On the first stage we look for a preliminary plan without taking into account the random
parameters. On the second stage the residual vector is searched for correction. The residual is
parameterized with using of a compensation matrix in new formulation which is proposed in the article,.
Basically it allows us to interpret the correcting process of the residuals in the following way – with the
emerging shortage of resources. An additional purchase is carried out without excess.</p>
      <p>The construction of the linear complementarity problem (12) and (13) for the second stage in a
new production of non-linear two-stage schemes stochastic transportation problem ensures the solvability
 and M  under all implementations b and any</p>
      <p>If matrix M  and M  is positive defined, then there is a unique solution of the linear
complementarity problem (12) and (13) in all implementations of the uncertainty parameters and
preliminary plan.
This research was supported by the Russian Foundation for Basic Research (project no.15-41-04436).</p>
      <p>Литература</p>
      <p>Владимировна, заведующий кафедрой прикладной математики и фундаментальной Омского
государственного технического университета, доктор физико-математических наук, профессор,
avzykina@mail.ru;</p>
    </sec>
  </body>
  <back>
    <ref-list>
      <ref id="ref1">
        <mixed-citation>
          <string-name>
            <surname>1 Hitchcock F. L.</surname>
          </string-name>
          <article-title>The distribution of product from several sources to numerous localities //</article-title>
          <source>Journal of Mathematical Physics</source>
          .
          <year>1941</year>
          . p.
          <fpage>224</fpage>
          -
          <lpage>230</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref2">
        <mixed-citation>
          2 Koopmans T.C.
          <article-title>Uptimum utilization of the transportation system</article-title>
          // Econometrica.
          <year>1949</year>
          . p.
          <fpage>3</fpage>
          -
          <lpage>4</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref3">
        <mixed-citation>
          <string-name>
            <surname>3 Williams A.C.</surname>
          </string-name>
          <article-title>A stochastic transportation</article-title>
          problem // Operations Research.
          <year>1963</year>
          . p.
          <fpage>759</fpage>
          -
          <lpage>770</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref4">
        <mixed-citation>
          4 Anholcer M.
          <article-title>Algorithm for the stochastic generalized transportation problem // Operations research</article-title>
          and decisions.
          <source>2012</source>
          . p.
          <fpage>9</fpage>
          -
          <lpage>20</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref5">
        <mixed-citation>
          <article-title>5 Anholcer М</article-title>
          .
          <article-title>Stochastic generalized transportation problem with discrete distribution of demand //Operations research</article-title>
          and decisions.
          <source>2013</source>
          . p.
          <fpage>9</fpage>
          -
          <lpage>19</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref6">
        <mixed-citation>
          <string-name>
            <surname>6 Biswal</surname>
            <given-names>M. P.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Samal . H. K. Stochastic Transportation</surname>
          </string-name>
          <article-title>Problem with Cauchy Random Variables</article-title>
          and Multi Choice Parameters //Physical Sciences.
          <year>2013</year>
          . p.
          <fpage>117</fpage>
          -
          <lpage>130</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref7">
        <mixed-citation>
          <string-name>
            <surname>7 Judin D.B</surname>
          </string-name>
          .
          <article-title>Matematicheskie metody upravlenija v uslovijah nepolnoj informacii. Zadachi i metody stohasticheskogo programmirovanija</article-title>
          . - M.:
          <string-name>
            <surname>Krasand</surname>
          </string-name>
          ,
          <year>2010</year>
          . - 400 p.
          <article-title>(In Russ)</article-title>
          .
        </mixed-citation>
      </ref>
      <ref id="ref8">
        <mixed-citation>
          <string-name>
            <surname>8 Bazara</surname>
            <given-names>M.</given-names>
          </string-name>
          <article-title>Nelinejnoe programmirovanie</article-title>
          .
          <source>Teorija i algoritmy</source>
          . - M. :
          <string-name>
            <surname>Mir</surname>
          </string-name>
          ,
          <year>1982</year>
          . - 583 p.
          <article-title>(In Russ)</article-title>
          .
        </mixed-citation>
      </ref>
      <ref id="ref9">
        <mixed-citation>
          9 Kaneva
          <string-name>
            <surname>O.N.</surname>
          </string-name>
          <article-title>Dvuhjetapnaja zadacha nelinejnogo stohasticheskogo programmirovanija s determinirovannoj matricej korrekcii //Matematicheskie struktury i modelirovanie</article-title>
          .
          <year>2005</year>
          . p.
          <fpage>25</fpage>
          -
          <lpage>33</lpage>
          . (In Russ).
        </mixed-citation>
      </ref>
      <ref id="ref10">
        <mixed-citation>
          10 Sakalauskas L.
          <article-title>Application of the Monte-Carlo Method to Nonlinear Stochastic Optimization with Linear Constraints /</article-title>
          /INFORMATICA.
          <year>2004</year>
          . №. 2. p.
          <fpage>271</fpage>
          -
          <lpage>282</lpage>
          .
        </mixed-citation>
      </ref>
    </ref-list>
  </back>
</article>