A bounded feasible standard-form LP attains its minimum at a basic feasible solution
ProvedPolyhedral.lp_min_attained_basicExistence of an optimal basic feasible solution. Consider the linear program in standard form
If it is feasible and its objective is bounded below on the feasible set, then the minimum is attained at a point whose support columns are linearly independent: the columns with form a linearly independent family.
This is the fundamental structural theorem of linear programming — the statement that optimisation may be restricted to basic feasible solutions, of which there are only finitely many. It underlies the simplex method, the finiteness of the set of candidate optima, and uniform bounds on optimal solutions as the right-hand side varies.
Formalization note. A basic feasible solution is described here directly by the linear independence of its support columns, LinearIndepOn \u211d (fun j => fun i => W i j) {j | y0 j \u2260 0}, rather than through a choice of basis matrix; this avoids assuming that has full row rank. Optimality is stated pointwise against all feasible points, and boundedness below as the existence of a single real lower bound for the objective on the feasible set.
import Mathlib open Matrix
theorem Polyhedral.lp_min_attained_basic {m n : ℕ} (W : Matrix (Fin m) (Fin n) ℝ)
(q : Fin n → ℝ) (d : Fin m → ℝ)
(hfeas : ∃ y : Fin n → ℝ, (∀ j, 0 ≤ y j) ∧ W.mulVec y = d)
(hbdd : ∃ beta : ℝ, ∀ y : Fin n → ℝ, (∀ j, 0 ≤ y j) → W.mulVec y = d →
beta ≤ q ⬝ᵥ y) :
∃ y0 : Fin n → ℝ, (∀ j, 0 ≤ y0 j) ∧ W.mulVec y0 = d ∧
(∀ y : Fin n → ℝ, (∀ j, 0 ≤ y j) → W.mulVec y = d → q ⬝ᵥ y0 ≤ q ⬝ᵥ y) ∧
LinearIndepOn ℝ (fun j => (fun i => W i j)) {j | y0 j ≠ 0} := by sorry