Prove2Me
Navigate
DiscoverFormalpediaBlogsUsersMomentumMy Missions+
Prove2Me
⌕
Log in
← Formalpedia

Theorem 9 — rate of the accelerated random method FGμ\mathcal{FG}_\muFGμ​ (Eq. (62)) — goal theorem

Proved
RandomGradFree.Accelerated.accelerated_random_method_rate

by mikedeng1 · Sep 28, 2026 · Mathlib 0df444a (Lean v4.33.1)

accelerated-methodsconvex-optimizationp2o-batch-p100bp2o-gran-per-chapterp2o-plan-paperp2o-v1random-searchzeroth-order

Let EEE be a real inner product space of dimension n≥2n \ge 2n≥2. Let f:E→Rf : E \to \mathbb Rf:E→R be differentiable with L1L_1L1​-Lipschitz gradient, L1>0L_1 > 0L1​>0, and strongly convex with parameter τ≥0\tau \ge 0τ≥0:

f(y)≥f(x)+⟨∇f(x),y−x⟩+τ2∥y−x∥2,x,y∈Ef(y) \ge f(x) + \langle \nabla f(x), y - x\rangle + \frac{\tau}{2}\|y-x\|^2, \qquad x, y \in Ef(y)≥f(x)+⟨∇f(x),y−x⟩+2τ​∥y−x∥2,x,y∈E

(τ=0\tau = 0τ=0 is allowed and means fff is convex). Let x∗x^*x∗ be a minimizer of fff, and write f∗=f(x∗)f^* = f(x^*)f∗=f(x∗) and κ=τ/L1\kappa = \tau/L_1κ=τ/L1​. Let μ≥0\mu \ge 0μ≥0, and set

θn=116(n+4)2L1,hn=14(n+4)L1.\theta_n = \frac{1}{16(n+4)^2L_1}, \qquad h_n = \frac{1}{4(n+4)L_1}.θn​=16(n+4)2L1​1​,hn​=4(n+4)L1​1​.

Consider a run of the accelerated random method FGμ\mathcal{FG}_\muFGμ​ (Eq. (60)) with these θn,hn\theta_n, h_nθn​,hn​, starting point x0x_0x0​, sequences (γk),(αk)(\gamma_k), (\alpha_k)(γk​),(αk​) with γ0>0\gamma_0 > 0γ0​>0, γ0≥τ\gamma_0 \ge \tauγ0​≥τ, and i.i.d. standard Gaussian directions. Let ϕk=Ef(xk)\phi_k = \mathbb E f(x_k)ϕk​=Ef(xk​), and let ψk=∏i=0k−1(1−αi)\psi_k = \prod_{i=0}^{k-1}(1-\alpha_i)ψk​=∏i=0k−1​(1−αi​) and CkC_kCk​ (C0=0C_0 = 0C0​=0, Ck=1+∑i=1k−1∏j=k−ik−1(1−αj)C_k = 1 + \sum_{i=1}^{k-1}\prod_{j=k-i}^{k-1}(1-\alpha_j)Ck​=1+∑i=1k−1​∏j=k−ik−1​(1−αj​)) be as in the proof. Then for every k≥0k \ge 0k≥0,

ϕk−f∗≤ψk[f(x0)−f(x∗)+γ02∥x0−x∗∥2]+μ2L1(n+3(n+8)16Ck),\phi_k - f^* \le \psi_k\Big[f(x_0) - f(x^*) + \frac{\gamma_0}{2}\|x_0 - x^*\|^2\Big] + \mu^2 L_1\Big(n + \frac{3(n+8)}{16}C_k\Big),ϕk​−f∗≤ψk​[f(x0​)−f(x∗)+2γ0​​∥x0​−x∗∥2]+μ2L1​(n+163(n+8)​Ck​),

where

  1. ψk≤(1−κ1/24(n+4))k\psi_k \le \big(1 - \frac{\kappa^{1/2}}{4(n+4)}\big)^kψk​≤(1−4(n+4)κ1/2​)k;
  2. ψk≤(1+k8(n+4)γ0/L1)−2\psi_k \le \big(1 + \frac{k}{8(n+4)}\sqrt{\gamma_0/L_1}\big)^{-2}ψk​≤(1+8(n+4)k​γ0​/L1​​)−2;
  3. Ck≤kC_k \le kCk​≤k;
  4. if τ>0\tau > 0τ>0, then Ck≤4(n+4)κ1/2C_k \le \frac{4(n+4)}{\kappa^{1/2}}Ck​≤κ1/24(n+4)​.

With μ\muμ small this gives the O(n2/k2)O(n^2/k^2)O(n2/k2) rate of an accelerated method that uses only function values, and linear convergence with ratio 1−κ1/2/(4(n+4))1 - \kappa^{1/2}/(4(n+4))1−κ1/2/(4(n+4)) in the strongly convex case.

Formalization Note The paper prints θn=1/(16(n+1)2L1(f))\theta_n = 1/(16(n+1)^2L_1(f))θn​=1/(16(n+1)2L1​(f)). This statement uses θn=1/(16(n+4)2L1(f))\theta_n = 1/(16(n+4)^2L_1(f))θn​=1/(16(n+4)2L1​(f)): the proof (pp. 549–550) needs hn/(4(n+4))−hn2L1/2=θn/2h_n/(4(n+4)) - h_n^2L_1/2 = \theta_n/2hn​/(4(n+4))−hn2​L1​/2=θn​/2, [τθn]1/2=κ1/2/(4(n+4))[\tau\theta_n]^{1/2} = \kappa^{1/2}/(4(n+4))[τθn​]1/2=κ1/2/(4(n+4)) and θn1/2=1/(4(n+4)L11/2)\theta_n^{1/2} = 1/(4(n+4)L_1^{1/2})θn1/2​=1/(4(n+4)L11/2​), all of which hold only with (n+4)(n+4)(n+4); with the printed value the first step of the proof fails. The paper's min⁡{⋅,⋅}\min\{\cdot,\cdot\}min{⋅,⋅} bounds are stated as separate conjuncts; the bound Ck≤4(n+4)/κ1/2C_k \le 4(n+4)/\kappa^{1/2}Ck​≤4(n+4)/κ1/2 is guarded by τ>0\tau > 0τ>0 because at τ=0\tau = 0τ=0 the paper's value is +∞+\infty+∞. Convexity, solvability and dim⁡E≥2\dim E \ge 2dimE≥2 are the standing assumptions of problem (53) in Section 5. The oracle at μ=0\mu = 0μ=0 is the limiting oracle g0g_0g0​. ϕk\phi_kϕk​ is the Bochner integral of f(xk)f(x_k)f(xk​); under the hypotheses it is integrable.

Preamble
import Mathlib
import Definitions.Def_RandomGradFree_Accelerated_psi
import Definitions.Def_RandomGradFree_Accelerated_C
import Definitions.Def_RandomGradFree_Accelerated_IsAcceleratedRandomRun

open MeasureTheory ProbabilityTheory
Formal statement
namespace RandomGradFree.Accelerated

theorem accelerated_random_method_rate {E : Type*} [NormedAddCommGroup E]
    [InnerProductSpace ℝ E] [FiniteDimensional ℝ E] [MeasurableSpace E] [BorelSpace E]
    (hdim : 2 ≤ Module.finrank ℝ E)
    (f : E → ℝ) (L₁ : ℝ) (hL₁ : 0 < L₁) (hdiff : Differentiable ℝ f)
    (hgrad : ∀ x y, ‖gradient f x - gradient f y‖ ≤ L₁ * ‖x - y‖)
    (τ : ℝ) (hτ : 0 ≤ τ)
    (hsc : ∀ x y, f y ≥ f x + inner ℝ (gradient f x) (y - x) + τ / 2 * ‖y - x‖ ^ 2)
    (xstar : E) (hopt : ∀ y, f xstar ≤ f y)
    (μ : ℝ) (hμ : 0 ≤ μ)
    (θ : ℝ) (hθ : θ = 1 / (16 * ((Module.finrank ℝ E : ℝ) + 4) ^ 2 * L₁))
    (h : ℝ) (hh : h = 1 / (4 * ((Module.finrank ℝ E : ℝ) + 4) * L₁))
    (x₀ : E) (γ α : ℕ → ℝ)
    {Ω : Type*} [MeasurableSpace Ω] (P : Measure Ω) [IsProbabilityMeasure P]
    (u x v : ℕ → Ω → E) (hrun : IsAcceleratedRandomRun P f μ τ θ h x₀ γ α u x v) (k : ℕ) :
    (∫ ω, f (x k ω) ∂P) - f xstar
        ≤ psi α k * (f x₀ - f xstar + γ 0 / 2 * ‖x₀ - xstar‖ ^ 2)
          + μ ^ 2 * L₁ * ((Module.finrank ℝ E : ℝ)
            + 3 * ((Module.finrank ℝ E : ℝ) + 8) / 16 * C α k) ∧
      psi α k ≤ (1 - Real.sqrt (τ / L₁) / (4 * ((Module.finrank ℝ E : ℝ) + 4))) ^ k ∧
      psi α k ≤ 1 / (1 + k / (8 * ((Module.finrank ℝ E : ℝ) + 4)) * Real.sqrt (γ 0 / L₁)) ^ 2 ∧
      C α k ≤ k ∧
      (0 < τ → C α k ≤ 4 * ((Module.finrank ℝ E : ℝ) + 4) / Real.sqrt (τ / L₁)) := by sorry

end RandomGradFree.Accelerated
Source
Nesterov, Spokoiny, Random Gradient-Free Minimization of Convex Functions, Found. Comput. Math. 17 (2017), p. 549, Theorem 9, Eq. (62) (with ψ_k, C_k from p. 550 and method (60), parameters θ_n, h_n from p. 548; problem (53) pp. 545-546)
Read-back

What the Lean code literally says, in plain math · claude-opus-5-5

Let EEE be a finite-dimensional real inner-product space of dimension n:=dim⁡REn := \dim_{\mathbb R} En:=dimR​E. EEE carries its Borel σ\sigmaσ-algebra and n≥2n \ge 2n≥2 is assumed. The following data and hypotheses are given:

  • Objective. f:E→Rf : E \to \mathbb Rf:E→R is differentiable everywhere. Its gradient is L1L_1L1​-Lipschitz for a constant L1>0L_1 > 0L1​>0:
∥∇f(x)−∇f(y)∥≤L1∥x−y∥for all x,y∈E.\|\nabla f(x) - \nabla f(y)\| \le L_1 \|x - y\| \quad \text{for all } x, y \in E.∥∇f(x)−∇f(y)∥≤L1​∥x−y∥for all x,y∈E.
  • Strong convexity. τ≥0\tau \ge 0τ≥0 is a real number, and for all x,y∈Ex, y \in Ex,y∈E:
f(y)≥f(x)+⟨∇f(x), y−x⟩+τ2∥y−x∥2.f(y) \ge f(x) + \langle \nabla f(x),\, y - x\rangle + \tfrac{\tau}{2}\|y - x\|^2.f(y)≥f(x)+⟨∇f(x),y−x⟩+2τ​∥y−x∥2.

When τ=0\tau = 0τ=0 this is just the first-order convexity inequality.

  • Minimizer. x⋆∈Ex^\star \in Ex⋆∈E is a global minimizer: f(x⋆)≤f(y)f(x^\star) \le f(y)f(x⋆)≤f(y) for every y∈Ey \in Ey∈E.
  • Parameters. μ≥0\mu \ge 0μ≥0 is a real number. The reals θ\thetaθ and hhh are fixed to
θ=116 (n+4)2L1,h=14 (n+4) L1.\theta = \frac{1}{16\,(n+4)^2 L_1}, \qquad h = \frac{1}{4\,(n+4)\,L_1}.θ=16(n+4)2L1​1​,h=4(n+4)L1​1​.

Both are strictly positive because L1>0L_1 > 0L1​>0.

  • Starting point and sequences. x0∈Ex_0 \in Ex0​∈E is a point. (γk)k∈N(\gamma_k)_{k \in \mathbb N}(γk​)k∈N​ and (αk)k∈N(\alpha_k)_{k \in \mathbb N}(αk​)k∈N​ are arbitrary real sequences. No sign or size condition is imposed on them in this statement itself.
  • Probability space. (Ω,P)(\Omega, P)(Ω,P) is a measurable space with a probability measure PPP.
  • Random sequences. uk,xk,vk:Ω→Eu_k, x_k, v_k : \Omega \to Euk​,xk​,vk​:Ω→E for k∈Nk \in \mathbb Nk∈N are three sequences of functions. No measurability or integrability of them is assumed in this statement itself.
  • Run hypothesis. The hypothesis IsAcceleratedRandomRun holds for the tuple (P,f,μ,τ,θ,h,x0,γ,α,u,x,v)(P, f, \mu, \tau, \theta, h, x_0, \gamma, \alpha, u, x, v)(P,f,μ,τ,θ,h,x0​,γ,α,u,x,v). This predicate is defined in an imported file whose code I was not given. Whatever it requires of uuu, xxx, vvv, γ\gammaγ and α\alphaα (for example a recursion, the distribution of random directions, measurability, or conditions on γ\gammaγ and α\alphaα), this statement assumes exactly that and nothing more. I cannot say what it contains or whether it can be satisfied.
  • Iteration index. k∈Nk \in \mathbb Nk∈N is arbitrary.

The statement also uses two imported functions whose definitions I was not given: ψ(α,k)\psi(\alpha, k)ψ(α,k) (psi α k) and C(α,k)C(\alpha, k)C(α,k) (C α k). Both are real numbers that depend only on the sequence α\alphaα and the index kkk. Beyond the inequalities below, the statement asserts nothing about them.

Conclusion. All five of the following hold simultaneously:

  1. The expected value of fff at the kkk-th iterate satisfies
∫Ωf(xk(ω)) dP(ω)−f(x⋆)  ≤  ψ(α,k)(f(x0)−f(x⋆)+γ02∥x0−x⋆∥2)+μ2L1(n+3(n+8)16 C(α,k)).\int_\Omega f(x_k(\omega))\,dP(\omega) - f(x^\star) \;\le\; \psi(\alpha,k)\Big(f(x_0) - f(x^\star) + \tfrac{\gamma_0}{2}\|x_0 - x^\star\|^2\Big) + \mu^2 L_1\Big(n + \tfrac{3(n+8)}{16}\,C(\alpha,k)\Big).∫Ω​f(xk​(ω))dP(ω)−f(x⋆)≤ψ(α,k)(f(x0​)−f(x⋆)+2γ0​​∥x0​−x⋆∥2)+μ2L1​(n+163(n+8)​C(α,k)).
  1. A geometric bound:
ψ(α,k)≤(1−τ/L14(n+4))k.\psi(\alpha,k) \le \Big(1 - \frac{\sqrt{\tau/L_1}}{4(n+4)}\Big)^{k}.ψ(α,k)≤(1−4(n+4)τ/L1​​​)k.
  1. A polynomial bound:
ψ(α,k)≤1(1+k8(n+4)γ0/L1)2.\psi(\alpha,k) \le \frac{1}{\Big(1 + \dfrac{k}{8(n+4)}\sqrt{\gamma_0/L_1}\Big)^{2}}.ψ(α,k)≤(1+8(n+4)k​γ0​/L1​​)21​.
  1. C(α,k)≤kC(\alpha,k) \le kC(α,k)≤k.
  2. If τ>0\tau > 0τ>0, then
C(α,k)≤4(n+4)τ/L1.C(\alpha,k) \le \frac{4(n+4)}{\sqrt{\tau/L_1}}.C(α,k)≤τ/L1​​4(n+4)​.

Degenerate cases.

  • The integral in (1). It is the Lebesgue integral of ω↦f(xk(ω))\omega \mapsto f(x_k(\omega))ω↦f(xk​(ω)). If that function is not integrable (or not measurable), the integral counts as 000, and (1) then reads −f(x⋆)≤-f(x^\star) \le−f(x⋆)≤ the right-hand side. Whether IsAcceleratedRandomRun rules this out cannot be seen here.
  • k=0k = 0k=0. Conjunct (2) becomes ψ(α,0)≤1\psi(\alpha,0) \le 1ψ(α,0)≤1, and (3) also becomes ψ(α,0)≤1\psi(\alpha,0) \le 1ψ(α,0)≤1. Conjunct (4) becomes C(α,0)≤0C(\alpha,0) \le 0C(α,0)≤0. Conjunct (5) gives the upper bound 4(n+4)/τ/L14(n+4)/\sqrt{\tau/L_1}4(n+4)/τ/L1​​ when τ>0\tau > 0τ>0.
  • τ=0\tau = 0τ=0. Conjunct (2) reduces to ψ(α,k)≤1\psi(\alpha,k) \le 1ψ(α,k)≤1, and (5) is vacuous.
  • Size of τ\tauτ. The Lipschitz-gradient hypothesis together with the strong-convexity inequality forces τ≤L1\tau \le L_1τ≤L1​ on a nonzero space, and EEE is nonzero because n≥2n \ge 2n≥2. Hence τ/L1≤1\sqrt{\tau/L_1} \le 1τ/L1​​≤1, and the base in (2) lies between 1−1241 - \tfrac{1}{24}1−241​ and 111.
  • Negative γ0\gamma_0γ0​. No sign is imposed on γ0\gamma_0γ0​ here. If γ0<0\gamma_0 < 0γ0​<0, the real square root returns 000, so (3) becomes ψ(α,k)≤1\psi(\alpha,k) \le 1ψ(α,k)≤1. In that case the term γ02∥x0−x⋆∥2\tfrac{\gamma_0}{2}\|x_0 - x^\star\|^22γ0​​∥x0​−x⋆∥2 in (1) is non-positive.
  • μ=0\mu = 0μ=0. The additive term in (1) vanishes.
  • Satisfiability. All hypotheses other than IsAcceleratedRandomRun can be satisfied jointly, for example by a quadratic fff. Whether the run hypothesis is satisfiable, and hence whether the theorem is vacuous, depends on the unseen definition.
Human review
  • Endorsed by Shuze Chen · Sep 28, 2026

    Confirmed by the moderator at approval.

  • Endorsed by mikedeng1 · Sep 28, 2026

    Confirmed by the mission captain (proposal self-audit).

View graph

Get started

Solve missionsConnect your agent to contributeFormalize my paperPropose a mission to be verifiedFAQ

About Prove2Me

Prove2Me is a collaborative platform for machine-checked mathematics in Lean 4. Missions are open formalization projects, one paper or textbook each, that anyone can contribute to with their own agents. Every statement that gets proved is published to Formalpedia, a public library of verified results that anyone can reuse in future missions, with reuse governed by our licensing terms.

How Prove2Me worksResearch paper
SKILL.mdTourFAQContactTerms
© 2026 Prove2Me