Theorem 10.6 — strictly complementary feasible solutions
OpenVanderbeiLP.StrictComp.strictly_complementary_feasibledualitylinear-programmingp2o-batch-b23ap2o-gran-per-chapterp2o-plan-bookp2o-v1strict-complementarity
Consider the linear program with explicit slacks
and its dual
with , , . If both the primal and the dual have feasible solutions, then there exist a primal feasible solution and a dual feasible solution such that
where means that every component of the vector is strictly positive. Since all four vectors are nonnegative, this says that for each at least one of , is positive, and for each at least one of , is positive.
This is the feasible-solution version of strict complementarity; the Strict Complementary Slackness Theorem strengthens it to optimal solutions.
Formalization Note The slacks are determined by the solutions: and .
Preamble
import Mathlib import Definitions.Def_VanderbeiLP_StrictComp_PrimalDualPair open Matrix
Formal statement
namespace VanderbeiLP.StrictComp
/-- **Vanderbei, Theorem 10.6 (p. 148).** If both the primal (10.9) and the dual (10.10) have
feasible solutions, then there are a primal feasible `x̄` (slack `w̄ = b - Ax̄`) and a dual
feasible `ȳ` (slack `z̄ = Aᵀȳ - c`) with `x̄ + z̄ > 0` and `ȳ + w̄ > 0` componentwise. -/
theorem strictly_complementary_feasible {m n : ℕ} (A : Matrix (Fin m) (Fin n) ℝ)
(b : Fin m → ℝ) (c : Fin n → ℝ)
(hP : ∃ x : Fin n → ℝ, PrimalFeasible A b x) (hD : ∃ y : Fin m → ℝ, DualFeasible A c y) :
∃ (xbar : Fin n → ℝ) (ybar : Fin m → ℝ),
PrimalFeasible A b xbar ∧ DualFeasible A c ybar ∧
(∀ j, 0 < xbar j + dualSlack A c ybar j) ∧
(∀ i, 0 < ybar i + primalSlack A b xbar i) := by sorry
end VanderbeiLP.StrictComp
Source
Vanderbei, Linear Programming: Foundations and Extensions, 4th ed., Springer 2014, p. 148, Theorem 10.6 and footnote 2 (PDF p. 161); Eqs. (10.9)–(10.10), p. 147
Human review
Confirmed by the mission captain (proposal self-audit).
Confirmed by the moderator at approval.