On Aug 4, 10:16 am, Henrik Ziegeldorf
<henrik.ziegeld
...@googlemail.com> wrote:
> > The MODI method is really just the pricing portion of the simplex
> > algorithm, tailored to the specific form of the Hitchcock problem.
> > /Paul
> Then, hopefully, my final question is:
> Why does the reduced costs compute to c_ij - u_i - v_j
> I guess the answer is somewhere close to the following considerations:
> The vector of reduced costs is
> c^T - c_B^T B^-1 A
> for every solution x to the primal problem (optimal or not)
> and we somehow can interprete c_B^T B^-1 as the values for the dual
> variables y
> so we can compute y , e.g. the u_i, v_j values, by solving
> c_ij - u_i - v_j = 0 for basisvariables x_ij (as they have
> reduced costs = 0)
> I know that for an optimal Basis B, the dual solution is just c_B^T
> B^-1
> but why does the interpretation of y = c_B^T B^-1 make sense for a non
> optimal basis?
> Or in other words, which is the meaning of y = c_B B^-1 for a non
> optimal basis?
> Thanks for helping
I am not sure I quite understand your question, or where you are
having trouble, so let me just summarize what happens. For the
transportation problem in equality-constraint form, we can get the
dual solution corresponding to any primal basis just by requiring that
each basic variable have a reduced cost equal to zero. That is, the
dual solution corresponding to a primal basis is determined by
complementary slackness: x_ij basic <---> dual constraint ij is an
equality. Why is this the case? The simplex method implies that at any
stage, the new objective equation = original objective equation + sum
w_i*(constraint equation i) for some constants w_i; in the TP case
this implies that new c_ij = c_ij - u_i - v_j for some constants u_i
(supply) and v_j (demand). We need c_ij = u_i + v_j for all basic ij,
so u_i and v_j are determined (up to an overall constant k, because
u_i + k and v_j - k gives the same result). Then, for the nonbasic ij
the reduced cost is c_ij - u_i - v_j. If these are all >= 0 we are
done: the simplex method stops at an optimal solution. However, if
some of these are < 0 we do not yet have an optimal basis although, in
the degenerate case, we might have an optimal solution. For
simplicity, assume no degenerate bases are encountered (since
degeneracy just complicates the fundamental issue at this point). I f
nonbasic variable x)ij ahas negative reduced cost, letting x_ij enter
the basis and letting some other currently-basic variable x_pq leave
the basis will give a strict total cost improvement. Prior to
optimality, you are working with primal feasible solutions
corresponding dual solutions are infeasible in the sense that some of
the dual constraints are violated. Once you arrive at primal
optimality, the dual solution becomes feasible as well.
R.G. Vickson