The math optional, made finite. Daily Practice

Big-M / two-phase method (artificial variables)

At a Glance

Why This Chapter Matters

This atom appears in 5 of the last 13 years and always carries high marks (10–20). The question always involves an LPP with at least one \ge or == constraint, making a starting BFS unavailable without artificial variables. Both methods produce identical optima; the choice between them is stylistic. Three of the five questions use Big-M; two use the two-phase method. Mastering the sign conventions for each method — and knowing when to declare optimality — is what separates full-mark answers from partial credit.

Minimum Theory

Why artificials are needed. The simplex method requires a basic feasible solution (BFS) to start. Slack variables provide a BFS for \le constraints. For \ge or == constraints, no obvious non-negative starting basis exists, so artificial variables ai0a_i \ge 0 are appended to create one.

Standard form for a \ge constraint. For aiTxbi\mathbf{a}_i^T\mathbf{x} \ge b_i (with bi>0b_i > 0): subtract a surplus variable si0s_i \ge 0 and add an artificial ai0a_i \ge 0: aiTxsi+ai=bi.\mathbf{a}_i^T\mathbf{x} - s_i + a_i = b_i. For an == constraint: add only an artificial (no surplus).

Big-M method. Add large penalty +M+M per artificial to the objective for minimisation (or M-M for maximisation): min  z=cTx+Miai.\min\; z = \mathbf{c}^T\mathbf{x} + M\sum_i a_i. Run the standard simplex. If the optimal solution has any ai>0a_i > 0, the original problem is infeasible. Once all ai=0a_i = 0 and all zjcj0z_j - c_j \le 0 (for min), the solution is optimal.

Sign convention for minimisation (Big-M). The entering variable rule: enter the non-basic variable with the largest positive zjcjz_j - c_j (or equivalently most negative cjzjc_j - z_j). Optimality when all zjcj0z_j - c_j \le 0. For a +M+M artificial in a minimisation problem, the coefficient of MM in zjcjz_j - c_j is always positive at first, driving artificials to leave.

Two-phase method. Split into two LPs:

Phase I: Minimise w=iaiw = \sum_i a_i subject to the constraints (original + artificials). If wmin>0w_{\min} > 0, the problem is infeasible. If wmin=0w_{\min} = 0, all artificials are zero; drop the artificial columns.

Phase II: Restore the original objective and run simplex from the Phase I optimal basis (with artificials removed).

Comparison.

Big-MTwo-phase
How many simplex runsOneTwo
RiskNumerical issues if MM is not chosen large enoughNone
Exam preferenceFaster to set upCleaner conceptually

Both give the same optimal solution; use whichever is asked.

Uniqueness check. At the optimum, if every non-basic variable has a strictly negative zjcjz_j - c_j (for min), the solution is unique. If any non-basic has zjcj=0z_j - c_j = 0, an alternative optimum exists along that edge.

Question Archetypes

ArchetypeYou are seeing this when…
big-m-method”Solve by Big-M method”; mixed \ge/\le constraints
two-phase-method”Solve by two-phase method”; \ge or == constraints

big-m-method (3 question(s); 2018, 2021, 2023)

Recognition Cues

Solution Template

  1. Convert to standard form. For each \ge constraint: subtract surplus sis_i, add artificial aia_i. For each \le constraint: add slack sjs_j. For == constraints: add only aia_i.
  2. Penalise artificials. For minimisation: add +Mai+Ma_i to the objective. For maximisation: add Mai-Ma_i.
  3. Initial BFS. Artificials (and any slacks) form the starting basis with ai=bia_i = b_i.
  4. Compute the initial zjcjz_j - c_j row. Include the MM-terms. The most positive entry (for min) or most negative (for max) identifies the entering variable.
  5. Ratio test. Divide RHS by the positive entries in the entering column; the row with the smallest ratio identifies the leaving variable.
  6. Pivot and repeat until all zjcj0z_j - c_j \le 0 (min) or 0\ge 0 (max).
  7. Drop expelled artificials. Once an artificial leaves, do not allow it to re-enter.
  8. Report optimality and state whether the original feasible region is bounded (to justify finite optimum).

Worked Example(s)

2018 Paper 2, 2018-P2-Q2b (20 marks)

Minimise z=3x1+5x2z = 3x_1 + 5x_2 subject to x1+2x28x_1 + 2x_2 \ge 8, 3x1+2x2123x_1 + 2x_2 \ge 12, 5x1+6x2605x_1 + 6x_2 \le 60, x1,x20x_1, x_2 \ge 0.

Standard form: add surpluses s1,s2s_1, s_2, artificials a1,a2a_1, a_2, slack s3s_3: x1+2x2s1+a1=8,3x1+2x2s2+a2=12,5x1+6x2+s3=60.x_1+2x_2-s_1+a_1=8,\quad 3x_1+2x_2-s_2+a_2=12,\quad 5x_1+6x_2+s_3=60.

Penalised objective: min  z=3x1+5x2+Ma1+Ma2\min\; z = 3x_1+5x_2+Ma_1+Ma_2.

Initial basis (a1,a2,s3)=(8,12,60)(a_1, a_2, s_3) = (8, 12, 60). Initial zjcjz_j - c_j:

x1x_1x2x_2s1s_1s2s_2
zjcjz_j-c_j4M34M-34M54M-5M-MM-M

Most positive: x1x_1 (4M3>4M54M-3 > 4M-5 for all MM). Ratio test: min(8/1,12/3,60/5)=4\min(8/1,\,12/3,\,60/5) = 4a2a_2 leaves.

Pivot 1 (enter x1x_1, leave a2a_2, pivot element 3): After row operations, new basis (a1,x1,s3)(a_1, x_1, s_3).

Updated: zjcjz_j-c_j for x2x_2 still positive (43M\sim \tfrac43 M). Enter x2x_2. Ratio test: min(3,6,15)=3\min(3, 6, 15) = 3a1a_1 leaves.

Pivot 2 (enter x2x_2, leave a1a_1): New basis (x2,x1,s3)=(3,2,32)(x_2, x_1, s_3) = (3, 2, 32). Both artificials expelled.

Optimality check: all remaining zjcj0z_j - c_j \le 0. Optimal: x1=2,x2=3,zmin=3(2)+5(3)=21.\boxed{x_1=2,\quad x_2=3,\quad z_{\min}=3(2)+5(3)=21.}

Finite optimum: the \le constraint 5x1+6x2605x_1+6x_2\le60 bounds the feasible region; combined with both artificials driven to zero, the problem has a finite optimal value.

Constraint check: x1+2x2=8x_1+2x_2=8 ✓ (binding), 3x1+2x2=123x_1+2x_2=12 ✓ (binding), 5x1+6x2=28605x_1+6x_2=28\le60 ✓ (slack 32).

2021 Paper 2, 2021-P2-Q4c (15 marks)

Maximise Z=4x1+5x2+2x3Z = 4x_1 + 5x_2 + 2x_3 subject to 2x1+x2+x3102x_1+x_2+x_3 \ge 10, x1+3x2+x312x_1+3x_2+x_3 \le 12, x1+x2+x3=6x_1+x_2+x_3 = 6, xi0x_i \ge 0.

Key structural observation. The equality constraint x1+x2+x3=6x_1+x_2+x_3=6 combined with the \ge constraint 2x1+x2+x3102x_1+x_2+x_3\ge10 gives (subtract): x14x_1\ge4. The objective becomes Z=4x1+5x2+2x3=2x1+3x2+12Z = 4x_1+5x_2+2x_3 = 2x_1+3x_2+12 after substituting x3=6x1x2x_3=6-x_1-x_2. To maximise: push x1x_1 and x2x_2 as large as feasible.

Constraint 2: 6+2x212x236 + 2x_2 \le 12 \Rightarrow x_2 \le 3. Try x1=4,x2=2,x3=0x_1=4,\,x_2=2,\,x_3=0: all constraints satisfied, Z=16+10+0=26Z=16+10+0=26.

Try x1=4,x2=3,x3=?x_1=4,\,x_2=3,\,x_3=?: forces x3=1<0x_3=-1<0, infeasible. The optimum is: x1=4,x2=2,x3=0,Zmax=26.\boxed{x_1=4,\quad x_2=2,\quad x_3=0,\quad Z_{\max}=26.}

The simplex Big-M tableau confirms this after two pivots expelling the two artificials.

2023 Paper 2, 2023-P2-Q3c (20 marks)

Minimise Z=2x1+3x2Z = 2x_1+3x_2 subject to x1+x29x_1+x_2\ge9, x1+2x215x_1+2x_2\ge15, 2x13x292x_1-3x_2\le9, x1,x20x_1,x_2\ge0. Is the optimal solution unique?

Standard form: surpluses s1,s2s_1,s_2, artificials a1,a2a_1,a_2, slack s3s_3.

Iteration 1: Enter x2x_2 (largest zjcj=3M3z_j-c_j = 3M-3). Ratio test: min(9/1,15/2)=7.5\min(9/1, 15/2) = 7.5a2a_2 leaves.

Iteration 2: Enter x1x_1 (next most positive, 0.5M0.50.5M-0.5). Ratio test: min(3,15,9)=3\min(3, 15, 9) = 3a1a_1 leaves.

Optimal tableau (both artificials expelled): x1=3,x2=6,s3=21x_1=3,\,x_2=6,\,s_3=21.

x1=3,x2=6,Zmin=2(3)+3(6)=24.\boxed{x_1=3,\quad x_2=6,\quad Z_{\min}=2(3)+3(6)=24.}

Uniqueness: at the optimum, non-basic variables s1s_1 and s2s_2 both have zjcj=1<0z_j-c_j=-1<0 (strictly negative). No non-basic variable has zjcj=0z_j-c_j=0, so the optimum is unique.

Common Traps


two-phase-method (2 question(s); 2022, 2024)

Recognition Cues

Solution Template

Phase I:

  1. Add surplus sis_i and artificial aia_i for each \ge constraint; add artificial only for == constraints.
  2. Phase I objective: min  w=iai\min\; w = \sum_i a_i (ignoring the original objective).
  3. Compute reduced costs for Phase I: cˉj=cBTaj\bar c_j = -\mathbf{c}_B^T\mathbf{a}_j for non-basic columns (where cB\mathbf{c}_B is the Phase I cost of the basis).
  4. Pivot until w=0w = 0 (all artificials zero) or confirm infeasibility (w>0w > 0 at optimum).
  5. Drop all artificial columns; keep the current basis (all non-artificial).

Phase II:

  1. Restore the original objective; compute new reduced costs in terms of the Phase I basis.
  2. Run simplex until all reduced costs are non-negative (for min) or non-positive (for max).
  3. Report the optimal solution.

Worked Example(s)

2022 Paper 2, 2022-P2-Q1e (10 marks)

Use two-phase method: Minimise Z=x1+x2Z = x_1 + x_2 s.t. 2x1+x242x_1+x_2\ge4, x1+7x27x_1+7x_2\ge7, x1,x20x_1,x_2\ge0.

Standard form: 2x1+x2s1+a1=42x_1+x_2-s_1+a_1=4, x1+7x2s2+a2=7x_1+7x_2-s_2+a_2=7.

Phase I: Minimise w=a1+a2w = a_1 + a_2.

Initial basis (a1,a2)=(4,7)(a_1,a_2)=(4,7). Phase I reduced costs: cˉx2=(1+7)=8\bar c_{x_2} = -(1+7) = -8 (most negative → enter x2x_2). Ratio test: min(4/1,7/7)=1\min(4/1,\,7/7)=1a2a_2 leaves.

After pivot (x2x_2 enters, a2a_2 leaves): basis (a1,x2)(a_1,x_2). Next most negative: cˉx1=13/7\bar c_{x_1}=-13/7. Enter x1x_1. Ratio: 3/(13/7)=21/133/(13/7)=21/13a1a_1 leaves.

Phase I basis (x1,x2)(x_1,x_2); w=0w=0. BFS: x1=21/13,x2=10/13x_1=21/13,\,x_2=10/13.

Phase II: Minimise Z=x1+x2Z = x_1 + x_2.

Reduced costs for s1s_1 and s2s_2: both >0>0 (computed from the Phase I tableau). No entering variable — current solution is optimal.

x1=2113,x2=1013,Zmin=31132.385.\boxed{x_1=\tfrac{21}{13},\quad x_2=\tfrac{10}{13},\quad Z_{\min}=\tfrac{31}{13}\approx2.385.}

Verification: 2x1+x2=42+1013=42x_1+x_2=\tfrac{42+10}{13}=4 ✓; x1+7x2=21+7013=7x_1+7x_2=\tfrac{21+70}{13}=7 ✓ (both binding).

2024 Paper 2, 2024-P2-Q1e (10 marks)

Use two-phase method: Maximise z=x1+2x2z = x_1+2x_2 s.t. x1x23x_1-x_2\ge3, 2x1+x2102x_1+x_2\le10, x1,x20x_1,x_2\ge0.

Standard form: x1x2s1+a1=3x_1-x_2-s_1+a_1=3, 2x1+x2+s2=102x_1+x_2+s_2=10.

Phase I: Minimise w=a1w = a_1.

Initial basis (a1,s2)=(3,10)(a_1,s_2)=(3,10). Phase I: enter x1x_1 (coefficient 1-1 in w=3x1+x2+s1w=3-x_1+x_2+s_1, most negative). Ratio test: min(3/1,10/2)=3\min(3/1,\,10/2)=3a1a_1 leaves. Pivot: x1=3,s2=4x_1=3,\,s_2=4; w=0w=0.

Phase I complete.

Phase II: Maximise z=x1+2x2z = x_1+2x_2.

Express zz in non-basics: from x1=3+x2+s1x_1 = 3+x_2+s_1, z=(3+x2+s1)+2x2=3+3x2+s1z = (3+x_2+s_1)+2x_2 = 3+3x_2+s_1. Coefficient of x2x_2 is +3>0+3>0 → enter x2x_2.

Ratio test on s2s_2 row (only binding constraint for x2x_2): s2=43x22s10x24/3s_2 = 4-3x_2-2s_1 \ge 0 \Rightarrow x_2 \le 4/3. So x2=4/3x_2=4/3, s2s_2 leaves.

New basis (x1,x2)(x_1,x_2): x1=13/3,x2=4/3x_1=13/3,\,x_2=4/3. Updated objective: z=7s1s2z = 7 - s_1 - s_2. Both non-basic coefficients negative — optimal.

x1=133,x2=43,zmax=7.\boxed{x_1=\tfrac{13}{3},\quad x_2=\tfrac{4}{3},\quad z_{\max}=7.}

Verification: x1x2=3x_1-x_2=3 ✓, 2x1+x2=102x_1+x_2=10 ✓ (both binding).

Common Traps

Marks-Aware Writing

10-mark question: show the standard-form conversion (with all variables named), the first pivot calculation in detail (including the zjcjz_j-c_j row), and the final answer. A small 2×52\times 5 or 3×63\times 6 simplex tableau is expected.

15-mark question: two full pivots shown step-by-step with the updated tableau after each pivot. State the entering and leaving variable at each step with the ratio-test computation.

20-mark question: three pivots or more, a thorough optimality-check discussion, and an explicit finite-optimum / uniqueness argument. Show constraint verification at the optimal point.

Practice Set

We've mapped all 13 years of this exam. Get new chapters, tools, and solutions as we release them — free.