Theorem 2, (20)–(22) — weak duality: and bound the reduced relaxation (5) from below
ProvedL0BnB.Duality.weak_dualitydual-boundslagrangian-dualityp2o-batch-pfp2bp2o-gran-per-chapterp2o-plan-paperp2o-v1weak-duality
Let , , , and let be the objective of the reduced relaxation (5), minimized over .
- If , then for every , every and every with ,
- If , then for every with for all , and every with ,
Here is the objective of (20)–(21) and that of (22). This is the "a dual of Problem (5) is given by" half of Theorem 2: every dual feasible point gives a valid lower bound on the relaxation, which is how the branch-and-bound algorithm prunes nodes from inexact primal solutions.
Formalization Note The dual (20) is unconstrained, so part 1 quantifies over all ; in (22) is not sign-constrained.
Preamble
import Mathlib import Definitions.Def_L0BnB_Duality_Dual
Formal statement
namespace L0BnB.Duality
/-- Theorem 2, (20)–(22), p. 13 — weak duality. If `√(λ₀/λ₂) ≤ M`, then `h₁(α, γ) ≤ L0BnB.Reduced.F(β)` for
every `α ∈ ℝⁿ`, every `γ ∈ ℝᵖ` (the dual (20) is unconstrained) and every `β` feasible for (5).
If `√(λ₀/λ₂) > M`, then `h₂(ρ, μ) ≤ L0BnB.Reduced.F(β)` for every `(ρ, μ)` satisfying the constraint of (22)
and every `β` feasible for (5). -/
theorem weak_duality {n p : ℕ} (X : Matrix (Fin n) (Fin p) ℝ) (y : Fin n → ℝ)
(lam0 lam2 M : ℝ) (hlam0 : 0 < lam0) (hlam2 : 0 < lam2) (hM : 0 < M) :
(Real.sqrt (lam0 / lam2) ≤ M →
∀ (α : Fin n → ℝ) (γ : Fin p → ℝ), ∀ β ∈ L0BnB.Reduced.box p M,
h1 X y lam0 lam2 M α γ ≤ L0BnB.Reduced.F X y lam0 lam2 M β) ∧
(M < Real.sqrt (lam0 / lam2) →
∀ (ρ : Fin n → ℝ) (μ : Fin p → ℝ), Feas22 X lam0 lam2 M ρ μ → ∀ β ∈ L0BnB.Reduced.box p M,
h2 y M ρ μ ≤ L0BnB.Reduced.F X y lam0 lam2 M β) := by sorry
end L0BnB.Duality
Source
Hazimeh, Mazumder, Saab, Sparse Regression at Scale: Branch-and-Bound rooted in First-Order Optimization, arXiv:2004.06152v2, p. 13, Theorem 2, (20)–(22); p. 31, Proof of Theorem 2, (48)–(49)
Human review
Confirmed by the mission captain (proposal self-audit).
Confirmed by the moderator at approval.