Name:

MATP6640/ISYE6770
Linear Programming
Spring 2018

Midterm Exam, Tuesday, March 27, 2018.
Solutions

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       /45



Q2       /20



Q3       /35






Total      /100
  1. (45 points) We have a two stage stochastic program where the nonnegative first stage variables x 2 must satisfy x 1 + 3x2 = 6. The first stage costs are cT x = 2x 1 + 7x2. There are two scenarios.
    1. (10 points) Assume x satisfies the first stage constraints. Show that in each scenario, there exists a recourse decision y that is feasible in the second stage problem.
    2. (10 points) Assume the cost of the second stage variables is 3y. Formulate the stochastic program in extended form as an LP. (Hint: your problem should have three equality constraints and four variables.)
    3. (10 points) Find all the basic feasible solutions to the LP formulation. Find the optimal solution.
    4. (5 points) What is the dual to the LP formulation?
    5. (10 points) Use complementary slackness to find the set of optimal solutions to the dual problem.

    Solution:

      • In Scenario 1: y = 13 - 2x1 - 5x2 13 - 2x1 - 6x2 = 1, so y 0.
      • In Scenario 2: y = 22 - 3x1 - 11x2 22 -11
 3x1 - 11x2 = 0, so y 0.
    1. Stochastic program is:
      minx,y     2x1   +   7x2  +   y1  +   2y2

subject to  x1   +   3x2                   =  6
           2x1   +   5x2  +   y1           =  13
           3x1   +  11x2          +    y2  =  22
                             x1, x2,y1,y2  ≥  0

    2. We know from a homework question that this problem has at most two extreme points. The extreme points are found by setting x1 = 0 or x2 = 0. They are:
      • x = (0, 2), y = (3, 0), with value 17.
      • x = (6, 0), y = (1, 4), with value 21.

      Hence, the optimal solution is x = (0, 2), with recourse variables y = (3, 0).

    3. Dual problem is
      max π       6π1  +   13π2  +   22π3
subject to   π1  +    2π2  +    3π3  ≤   2
            3π   +    5π   +   11π   ≤   7
              1         2         3
                       π2            ≤   1
                                 π3  ≤   2
                          π1, π2,π3      free

    4. The optimal primal solution is degenerate, so there may be multiple optimal dual solutions. From complementary slackness and dual feasibility, we need:
       π   +   2π   +   3 π   ≤  2
  1        2         3
3π1  +   5π2  +  11 π3  =  7
          π2            =  1
                    π3  ≤  2

      We can use the two equalities to set π2 = 1 and then express π1 = 1
3(2 - 11π3). The two inequalities then become:

      1-(2 - 11 π ) + 3π   ≤  0
3         3      3
               π3   ≤  2
      or
      1 ≤ π3 ≤ 2

      Thus, the set of dual optimal solutions is the line segment:

                               (  20     )
π = λ (- 3,1, 1) + (1 - λ) ----,1,2  ,     for 0 ≤ λ ≤  1.
                             3

  2. (20 points) Let A be a real p × n matrix and B be a real q × n matrix, with p, q, and n all positive integers. Show that exactly one of the following systems has a solution:

    System 1: Ax < 0, Bx = 0 for some x n.
    System 2: AT u + BT v = 0 for some u p, v q, with u 0 and u0.

    Solution:

    Let e denote the vector of ones. Consider the primal-dual pair of LPs:

    maxx          0                       minu,v      - eTu
subject to  Ax   ≤   - e    (P )      subject to   AT u  +   BT v  =   0      (D )

            Bx   =     0                              u            ≥   0
              x      free                                        v      free

    1. Assume System 1 is consistent: Any feasible solution to System 1 can be scaled to give a feasible solution to (P). Thus, the optimal value of (P) is 0, so (D) has optimal value 0, so System 2 is inconsistent.
    2. Assume System 1 is inconsistent: Then (P) is infeasible. Problem (D) is always feasible: take u = 0 and v = 0. Thus, (D) must have unbounded optimal value. So there exists u, v with AT u + BT v = 0, u 0, eT u > 0, so System 2 is consistent.
  3. (35 points) In a vehicle routing problem, a set of cities {1,,n} must be visited, starting from a base city 0. More than one vehicle can be used, with each vehicle traversing a subtour that starts and ends at the base but only visits a subset of the cities. This problem can be formulated as an integer program with each variable corresponding to an allowable subtour. Let T be the set of allowable subtours. Let ct be the cost of subtour t T. Let at ∈{0, 1}n be the incidence vector of subtour t T, so ait = 1 if and only if tour t visits city i for i ∈{1,,n}. The LP relaxation of the vehicle routing problem can then be expressed as
                ∑
min         ∑ t∈T ctxt
subject to    t∈T atixt  ≥   1  for i = 1,...,n      (CP )
                   xt  ≥   0  ∀t ∈ T

    Let n = 4. Assume the cost ct of a tour is given by the sum of the edge lengths in the tour, and the edge lengths are as follows:

    city01234






    04794
    14468
    27475
    39677
    44857

    No tour is allowed to visit more than three cities (plus the base city). A tour can visit just one city and the base city. So, for example, valid tours include 0 - 4 - 0 with cost 8 and 0 - 1 - 2 - 3 - 0 with cost 24.

    1. (5 points) What is the dual LP to (CP)?
    2. (10 points) Let the current dual solution be y = (6, 8, 10, 7)T . Find a violated dual constraint.
    3. (10 points) There are 14 feasible tours for this problem. Enumerate them all and hence show that ŷ = (4, 8, 12, 5)T is dual feasible.
    4. (10 points) What does complementary slackness imply about the primal solution if the dual solution is ŷ? Show there is a primal feasible solution satisfying complementary slackness with each positive xt taking the same value. What do you conclude?

    Solution:

    1. Dual problem is
                  ∑
maxy          ni=1yi
subject to   (at)Ty  ≤   ct ∀t ∈ T     (CD )

                  y  ≥   0

    2. Need a subset of the cities with the length ct smaller than the sum of the y values for the cities in the subset. Take cities {1, 2, 4}. Then ct = 17, but (at)T y = 21. The violated constraint is
      y  + y  + y  ≤ 17.
  1   2    4

    3. Enumerate the various subsets:

      subset t ityict(at)T y(     )
{  <  }
   =
(  >  )ct




      1 4 8 <
      2 8 14 <
      3 12 18 <
      4 5 8 <
      1,2 12 15 <
      1,3 16 19 <
      1,4 9 16 <
      2,3 20 23 <
      2,4 13 16 <
      3,4 17 20 <
      1,2,3 24 24 =
      1,2,4 17 17 =
      1,3,4 21 21 =
      2,3,4 25 25 =
    4. From complementary slackness, the only primal xt variables that can be positive are those corresponding to the subsets of size 3. Also, since each ŷi > 0, all the primal constraints must hold at equality. So the primal variables must satisfy:
      x123  +  x124  +   x134            =  1
x123  +  x124           +   x234  =  1
x123           +   x134  +   x234  =  1

         x124  +   x134  +   x234  =  1

      The solution to this system is x123 = x124 = x134 = x234 = 13, with all other xt equal to zero. We have primal feasibility, dual feasibility and complementary slackness, so we have optimality.

      Check objective function values:

      • Dual value: 4 + 8 + 12 + 5 = 29.
      • Primal value: 1
3(24 + 17 + 21 + 25) = 1
3(87) = 29.

      So primal and dual objective function values agree.