Solutions

MATP6640/ISYE6770
Linear Programming
Spring 2016

Midterm Exam, Friday, April 8, 2016.

Please do all three problems. Show all work. No books or calculators allowed. You may use any result from class, the homeworks, or the texts, except where stated. You may use one sheet of handwritten notes. The exam lasts 110 minutes.

Q1       /30



Q2       /35



Q3       /35






Total      /100
  1. (30 points; each part is worth 10 points)
    Let A be an n × n matrix satisfying:

    Show the following:

    1. Let y n. Assume the smallest component of y is y k, so yk yi for i = 1,,n. Show i=1na ikyi 0.
    2. Show the linear program
      maxy           0
subject to  AT y  ≤  - e
               y  ∈  ℝn

      is infeasible, where e denotes the vector of ones.

    3. What can you say about the set {x n : Ax = 0,x 0}?

    Solutions:

    1. We have
      ∑n            ∑n
    aikŻyi ≥       aikŻyk   since aik ≥ 0 and Żyi ≥ Żyk for i ⁄= k
 i=1           i=1
                ∑ n
          =   Żyk    aik
                 i=1
          =   0.
    2. For any y n, let y k be the smallest component of y. We have
      (    )      ∑ n
 AT y    =      aikŻyi
      k      i=1
         ≥  0   from  part (a)
      so (AT y) k ⁄≤-1, so AT y ⁄≤-e.
    3. The dual LP to the infeasible problem in part (b) is
                    T
minx ∈ℝn    - e x
subject to    Ax   =   0
                x  ≥   0.

      This problem is feasible (take x = 0), so it is unbounded. So there exists a nonzero x in the given set.

  2. (35 points)
    In the single commodity network flow problem in the figure, nodes s and t are supply nodes and nodes u and v are demand nodes. The supplies available are bs = 3 and bt = 2. The demands are du = 1 and dv = 4. The unit shipping costs for each edge are indicated in the figure.

    An initial basic feasible solution can be obtained by taking xsu, xut, and xtv to be basic, and xst, xuv to be nonbasic.

    1. (5 points) What is the initial basic feasible solution?
    2. (15 points) What are the primal and dual linear programs? Find a dual solution that satisfies complementary slackness. Show your dual solution is not dual feasible.
    3. (15 points) Make a simplex pivot to improve your primal solution. What is the updated dual solution? Is the new solution optimal?

    Solutions:

    1. xsu = 3, xut = 2, xtv = 4, xst = xuv = 0.
    2. Total supply and total demand are both equal to 5, so we can use equality constraints in the primal LP. Primal LP:
      minx        5xsu  +   9xst  +   5xut +   7xuv  +   3xtv
subject to    xsu +    xst                               =     3

            - xsu           +    xut +    xuv            =   - 1
                  -    xst  -    xut           +    xtv  =     2
                                     -    xuv  -    xtv  =   - 4
                                 xij ≥      0      ∀ edges (i,j).

      Dual LP:

      maxy        3ys  -  yu  +   2yt  -  4yv
subject to   ys  -  yu                    ≤  5
             ys         -    yt           ≤  9
                    y   -    y            ≤  5
                     u        t
                    yu           -    yv  ≤  7
                             yt  -    yv  ≤  3

      From complementary slackness, need y to satisfy

      ys  -   yu                 =   5
        y   -   y          =   5
         u       t
                yt -   yv  =   3

      One solution is yv = 0, yt = 3, yu = 8, ys = 13. This violates the second and fourth dual constraints, each by one unit.

    3. We use network simplex to make the pivot. We choose edge (s,t) to enter the basis (using the edge (u,v) instead is similar). The fundamental cycle also includes edges (s,u) and (u,t). The values of the variables with different choices Δ for the incoming basic variable is indicated in the figure:

      Largest possible value of Δ = 2, and edge (u,t) leaves the basis. Get new bfs as indicated in the figure:

      From complementary slackness, need y to satisfy:

      ys  -   yu                 ≤   5
ys          -   yt         ≤   9
                yt -   yv  ≤   3

      One solution is yv = 0, yt = 3, ys = 12, yu = 7, which is dual feasible. Therefore, the new bfs is primal optimal, and the given dual solution is dual optimal.

  3. (35 points)
    The numbers on the edges in the following graph represent the probability of failure of an edge. Edges fail independently of one another.

    For a price, we can reinforce an edge so that its probability of failure drops to zero. We want to ensure that a path exists between s and t with high probability. More precisely, we want to ensure that, for each cutset,

    Probability {every edge in cutset fails} ≤ 0.01.
    (1)

    A cutset consists of all edges leading from P to Q, where P and Q form a partition of the vertices with s P and t Q, so

    P ∪ Q  = {s,t,u,v,w },  P  ∩ Q =  ∅,  s ∈ P,   t ∈ Q.
    (2)

    Two possible partitions are:

    We introduce variables

         {
x =     1  if edge e is reinforced
 e      0  otherwise.
    (3)

    In order to ensure the solution x satisfies the constraint (1) for the cutset given by P1 and Q1, we can impose the constraint

    xsu + xsv ≥ 1.

    Similarly, we can impose the constraint

    xvt + xwt ≥ 1

    corresponding to the other cutset. For a general graph, there may be a large number of possible cutsets, so we add the corresponding constraints as needed. The initial LP relaxation is as follows, where the objective function coefficients correspond to the costs of reinforcing a particular edge:

    min    7    3x    +   2x    +  7x     +  6x     +  7x     +   5x   +   3x
    x∈ℝ       su        sv        uv        uw        vw        vt       wt
subject to   xsu  +    xsv                                                   ≥   1
                                                               xvt +    xwt  ≥   1
              xe  ≥      0  for all edges e
    (4)

    1. (5 points) Show that the solution to the initial relaxation is xsv = xwt = 1, with all other xe = 0.
    2. (10 points) Find a cutset that leads to a valid inequality that is violated by the solution in part 3a.
    3. (5 points) Show that xsu = xwt = 1 with all other xe = 0 is feasible after the new valid inequality is added to (4).
    4. (15 points) Construct the dual LP to the new problem obtained by adding the constraint to (4). Find a dual solution that satisfies complementary slackness. Is your solution optimal?

    Solution:

    1. Trivial.
    2. Choose P = {s,v}, Q = {t,u,w}. The edges going from P to Q are (s,u), (v,t) and (v,w).

      The product of the probabilities of failure of the three edges in the cutset is

      P (all edges from P to Q fail) = 0.2 × 0.3 × 0.2 = 0.012 > 0.01,

      so we must reinforce at least one edge in this cutset in order to satisfy constraint (1). This gives the valid constraint

      xsu + xvt + xvw ≥ 1,
      (5)

      which is violated by the solution given in part (a).

    3. The updated primal LP is
      minx ∈ℝ7    3xsu  +   2xsv  +  7xuv   +  6xuw   +  7xvw   +  5xvt  +   3xwt
subject to   xsu  +    xsv                                                   ≥   1
                                                               x   +    x    ≥   1
                                                                vt       wt
             xsu                      +   xvw             +    xvt           ≥   1
              xe  ≥      0  for all edges e
      (6)

      and the given solution is feasible in this LP, with value 6.

    4. The dual LP is
      max         y1  +  y2  +   y3
subject to  y          +   y   ≤  3
             1              3
            y1                 ≤  2
                            0  ≤  7
                           y3  ≤  6
                            0  ≤  7

                   y2  +   y3  ≤  5
                   y2          ≤  3
                     y1,y2,y3  ≥  0.

      From complementary slackness, we need

      y1        +  y3  =   3
      y          =   3
       2

      This system has multiple solutions, including y = (2, 3, 1) and (1, 3, 2), both of which are feasible in the dual problem. Since we have primal and dual feasibility as well as complementary slackness, we are optimal.

      The general set of dual optimal solutions has the form y = (1 + z, 3, 2 - z) for 0 z 1.