Prove2Me
Navigate
DiscoverCollectionsFormalpediaBlogsUsersMomentumMy Missions+
Prove2Me
⌕
Log in

Get started

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

Machine Learning

653 missions · 436 completed

The science of systems that learn from data and experience. Its scope runs from the statistical and mathematical foundations of learning, including generalization, expressivity, and computational limits, through the design of learning algorithms, deep learning, reinforcement learning, and probabilistic methods, to the empirical study of large models and the trustworthiness, interpretability, and societal impact of learned systems.

Missions

Open217Completed436All653
🏆Completed
Convex OptimizationOptimization·Captain: mikedeng1

A Modern Introduction to Online Learning XIII: Online Mirror Descent Has Dynamic Regret (B_ψ(u_T; x_1) + Q·P(u_1, …, u_T))/η + (η/(2λ)) Σ_t ‖g_t‖⋆²Textbook

Why moving comparators matter

In online learning, a learner chooses a point before seeing each round's loss. Ordinary regret compares the learner with one fixed point chosen after all rounds. That comparison can be too coarse when a useful action changes over time. A portfolio allocation, a prediction rule, or an operating point may need to move after the environment shifts. Dynamic regret compares the learner with a specified sequence of feasible points, one for each round. The sequence is allowed to move, but its movement appears explicitly in the bound. Chapter 14 of Orabona's Online Learning: A Modern Introduction Using Convex Optimization develops this comparison, then studies strongly adaptive and tracking regret. This mission isolates its first capstone, Theorem 14.2.

The result extends the book's earlier analysis of online mirror descent to changing comparators. As the chapter explains, a moving comparator can achieve a smaller cumulative loss than any single fixed point, so a dynamic bound answers a stricter performance question. The bound also makes its own limit clear: frequent or large comparator changes cost more. The chapter relates this framework to earlier work on drift in online learning; the formal target here is the precise statement in the pinned arXiv version, with its explicit constants and final Bregman term. Source, Chapter 14, pp. 227–229.

Setting and notation

Let EEE be a finite-dimensional real normed space, V⊆X⊆EV\subseteq X\subseteq EV⊆X⊆E, and VVV a nonempty closed convex feasible set. On round t=1,…,Tt=1,\ldots,Tt=1,…,T, the learner plays xt∈Vx_t\in Vxt​∈V, receives a real-valued loss ℓt\ell_tℓt​, and selects a subgradient gtg_tgt​ at xtx_txt​. A comparison sequence u1,…,uTu_1,\ldots,u_Tu1​,…,uT​ consists of arbitrary points in VVV. The dynamic regret is

DRegret⁡T(u1,…,uT)=∑t=1T(ℓt(xt)−ℓt(ut)).\operatorname{DRegret}_T(u_1,\ldots,u_T) =\sum_{t=1}^{T}\bigl(\ell_t(x_t)-\ell_t(u_t)\bigr).DRegretT​(u1​,…,uT​)=t=1∑T​(ℓt​(xt​)−ℓt​(ut​)).

The path length of the comparison sequence is the total distance between successive points:

P(u1,…,uT)=∑t=2T∥ut−ut−1∥.P(u_1,\ldots,u_T)=\sum_{t=2}^{T}\lVert u_t-u_{t-1}\rVert.P(u1​,…,uT​)=t=2∑T​∥ut​−ut−1​∥.

It is zero for a constant sequence. The path starts at round two, so neither u0u_0u0​ nor an artificial first movement enters the claim. The norm on EEE is arbitrary. Its dual norm on continuous linear functionals is written ∥⋅∥∗\lVert\cdot\rVert_*∥⋅∥∗​.

Online mirror descent uses a differentiable regularizer ψ\psiψ and its Bregman divergence

Bψ(a;b)=ψ(a)−ψ(b)−⟨∇ψ(b),a−b⟩.B_\psi(a;b)=\psi(a)-\psi(b)-\langle\nabla\psi(b),a-b\rangle.Bψ​(a;b)=ψ(a)−ψ(b)−⟨∇ψ(b),a−b⟩.

At each round, the next point minimizes ⟨gt,z⟩+Bψ(z;xt)/η\langle g_t,z\rangle+B_\psi(z;x_t)/\eta⟨gt​,z⟩+Bψ​(z;xt​)/η over z∈Vz\in Vz∈V, with one constant step size η>0\eta>0η>0. The regularizer is closed and λ\lambdaλ-strongly convex, with λ>0\lambda>0λ>0, on the stated domain. Its derivative must be evaluated at interior points of XXX. For a horizon T≥1T\ge1T≥1, the theorem measures the range of the gradients seen along the run by the finite maximum

Q=max⁡1≤t≤T∥∇ψ(xt)−∇ψ(x1)∥∗.Q=\max_{1\le t\le T}\lVert\nabla\psi(x_t)-\nabla\psi(x_1)\rVert_*.Q=1≤t≤Tmax​∥∇ψ(xt​)−∇ψ(x1​)∥∗​.

The maximum includes x1x_1x1​, at which the difference is zero, and stops at xTx_TxT​. These conventions are the ones in Definition 14.1 and Theorem 14.2.

Formalization targets

One-step mirror descent estimate

The first milestone is the inequality part of Lemma 6.10. For each feasible comparator, one OMD update controls the loss difference by a decrease in Bregman divergence plus a quadratic dual-norm term. The milestone records the book's whole inequality chain, including the intermediate term that exposes the change in the iterate. It is the local estimate used at every round of Theorem 14.2. Lemma 6.10, p. 67.

Summed comparison with an extra endpoint

The second milestone is the display in the proof of Theorem 14.2 immediately before the final choice uT+1=uTu_{T+1}=u_TuT+1​=uT​. It keeps that additional comparator free and shows the gradient-difference sum explicitly. This intermediate statement is useful for inspecting the sign and time index of each comparator shift. Proof of Theorem 14.2, pp. 228–229.

Goal: dynamic regret of online mirror descent

The goal is the full bound of Theorem 14.2:

DRegret⁡T(u1,…,uT)≤Bψ(uT;x1)+QP(u1,…,uT)η+η2λ∑t=1T∥gt∥∗2−Bψ(uT;xT+1)η.\operatorname{DRegret}_T(u_1,\ldots,u_T) \le \frac{B_\psi(u_T;x_1)+Q P(u_1,\ldots,u_T)}{\eta} +\frac{\eta}{2\lambda}\sum_{t=1}^{T}\lVert g_t\rVert_*^2 -\frac{B_\psi(u_T;x_{T+1})}{\eta}.DRegretT​(u1​,…,uT​)≤ηBψ​(uT​;x1​)+QP(u1​,…,uT​)​+2λη​t=1∑T​∥gt​∥∗2​−ηBψ​(uT​;xT+1​)​.

It applies to every feasible comparator sequence, including a constant one. The negative terminal divergence is part of the target. Retaining it makes the statement match the source and preserves information for later use. Theorem 14.2, p. 228.

What the result gives

The theorem quantifies how much a comparator may move before its performance becomes hard to match. The cost of movement is linear in the path length and scaled by the regularizer-gradient range QQQ. With a bounded dual norm for the observed subgradients and a suitable constant step size, the chapter derives a square-root dependence on the horizon for slowly moving sequences. A static comparator has path length zero and recovers the familiar mirror-descent setting. The bound also applies to arbitrary norms and regularizers, so it is not limited to Euclidean projection. Discussion after Theorem 14.2, p. 229.

The mathematical result is proved in the book. The remaining task is to give its exact assumptions, update rule, intermediate inequality, and final inequality machine-checked proofs. The already published Bregman definition supplies a reusable representation of BψB_\psiBψ​; the new run, path-length, and gradient-maximum definitions supply the chapter's missing interface. Those definitions can support later formalizations of tracking regret and other moving-comparator results. No proof of the present chapter result is claimed by these draft statements.

Difficulty

For a fixed comparator, the Bregman terms from successive rounds line up directly. With utu_tut​ changing, they do not: each transition changes the first argument of BψB_\psiBψ​. The resulting correction depends on differences of regularizer gradients, and its time index must agree with the path-length sum. A simple replacement of the static comparator by utu_tut​ in a static regret theorem therefore leaves an unaccounted term. The terminal point xT+1x_{T+1}xT+1​ is also essential to the full telescoping statement. The mission asks for the book's exact constant and terminal term rather than only an asymptotic rate.

Formalization scope

The Lean model uses a finite-dimensional real normed space. Losses and the regularizer are real-valued on their evaluated domains; the feasible-set and domain predicates represent the book's extended-real convention outside those domains. Closedness of the regularizer is represented by a closed epigraph over XXX. Subgradients are continuous linear functionals, so their operator norm is the dual norm. The run predicate represents every admissible subgradient and argmin choice for Algorithm 6.1; it also carries the required membership of x1,…,xT+1x_1,\ldots,x_{T+1}x1​,…,xT+1​ in V∩int⁡XV\cap\operatorname{int}XV∩intX. The run is a predicate on sequences, not a selected canonical sequence.

The book's round labels start at one. Index zero is unused. The horizon is at least one, and η,λ\eta,\lambdaη,λ are positive. These conditions keep the finite maximum nonempty and every denominator meaningful. The subgradient condition is relative to VVV, which is sufficient for the theorem's comparisons and slightly weaker than the book's full-space subgradient convention. This choice is disclosed in the item statements. The comparator sequence is universally quantified through its feasibility condition; it is never restricted to a constant sequence or assumed to satisfy the conclusion. The definitions and theorem statements invite proofs of the one-step inequality, the summed estimate, and the final bound, as well as reusable results about Bregman divergence on arbitrary normed spaces.

Selected references

  • Francesco Orabona, Online Learning: A Modern Introduction Using Convex Optimization, arXiv:1912.13213v10, 2026. arXiv
7 thms3 active usersReviewed
🏆Completed
Linear OptimizationStatistics·Captain: mikedeng1

Fast Rates for Contextual Linear Optimization 3: Under the Margin Condition, a Plug-In Policy Whose Estimate Has Pointwise Sub-Gaussian Tails at Rate aₙ Has Regret O(B·aₙ^(−(1+α)/2))Research Paper

Why estimate-then-optimize

Many operational decisions are linear programs whose cost vector is uncertain but predictable from side information: routing with travel times that depend on weather and time of day, procurement with prices that depend on market signals, assignment with costs that depend on observed features. Contextual linear optimization (CLO) formalizes this: observe features XXX, choose a decision zzz in a polytope Z\mathcal ZZ, pay Y⊤zY^\top zY⊤z for an unobserved random cost vector YYY.

Two families of data-driven methods compete. Estimate-then-optimize (ETO) fits a regression f^(x)≈E[Y∣X=x]\hat f(x) \approx \mathbb E[Y \mid X = x]f^​(x)≈E[Y∣X=x] and then solves the linear program with the estimated costs. Integrated methods (for example the SPO+ approach of Elmachtoub and Grigas, 2022) choose a decision rule by minimizing the empirical decision cost directly over a policy class, on the argument that a regression fit ignores how its errors affect decisions. Hu, Kallus and Mao (arXiv:2011.03030v3, Management Science 2022) show that, under a margin condition on the problem, the comparison goes the other way: ETO attains faster regret rates than any method that only uses a policy class.

This mission formalizes their ETO upper bound, Theorem 8. It is the third of a three-mission series on the paper; the first formalizes the upper bound for empirical risk minimization over a policy class (Theorem 6) and the second the matching lower bound for policy-class methods (Theorem 7).

Setting

Write Rk\mathbb R^kRk for Euclidean space with inner product u⊤vu^\top vu⊤v and norm ∥⋅∥\|\cdot\|∥⋅∥.

  • Decisions. The polytope Z={z∈Rd:Az≤b}\mathcal Z = \{z \in \mathbb R^d : Az \le b\}Z={z∈Rd:Az≤b} is nonempty and satisfies ∥z∥≤B\|z\| \le B∥z∥≤B for all z∈Zz \in \mathcal Zz∈Z, with B>0B > 0B>0. Its set of extreme points is Z∠\mathcal Z^\angleZ∠, a finite set.
  • Data. Features X∈RpX \in \mathbb R^pX∈Rp have law PX\mathbb P_XPX​; costs Y∈RdY \in \mathbb R^dY∈Rd satisfy ∥Y∥≤1\|Y\| \le 1∥Y∥≤1. The regression function is f∗(x)=E[Y∣X=x]f^*(x) = \mathbb E[Y \mid X = x]f∗(x)=E[Y∣X=x]. The data D\mathcal DD are nnn i.i.d. copies of (X,Y)(X, Y)(X,Y).
  • Optimal decisions. Z∗(x)=arg⁡min⁡z∈Zf∗(x)⊤z\mathcal Z^*(x) = \arg\min_{z \in \mathcal Z} f^*(x)^\top zZ∗(x)=argminz∈Z​f∗(x)⊤z.
  • Gap. If Z∗(x)≠Z\mathcal Z^*(x) \ne \mathcal ZZ∗(x)=Z, the gap is
Δ(x)=min⁡z∈Z∠∖Z∗(x)f∗(x)⊤z−min⁡z∈Zf∗(x)⊤z,\Delta(x) = \min_{z \in \mathcal Z^\angle \setminus \mathcal Z^*(x)} f^*(x)^\top z - \min_{z \in \mathcal Z} f^*(x)^\top z,Δ(x)=z∈Z∠∖Z∗(x)min​f∗(x)⊤z−z∈Zmin​f∗(x)⊤z,

the extra cost of the best suboptimal vertex; otherwise Δ(x)=0\Delta(x) = 0Δ(x)=0.

  • Noise condition (Assumption 2). For some α,γ≥0\alpha, \gamma \ge 0α,γ≥0,
PX(0<Δ(X)≤δ)≤(γδ/B)αfor all δ>0.\mathbb P_X\bigl(0 < \Delta(X) \le \delta\bigr) \le (\gamma\delta/B)^\alpha \quad \text{for all } \delta > 0.PX​(0<Δ(X)≤δ)≤(γδ/B)αfor all δ>0.

It always holds with α=0\alpha = 0α=0; larger α\alphaα means less mass near dual degeneracy.

  • Plug-in policy. For an estimate fff, a plug-fff-in policy πf\pi_fπf​ chooses at every xxx an extreme point in arg⁡min⁡z∈Zf(x)⊤z\arg\min_{z \in \mathcal Z} f(x)^\top zargminz∈Z​f(x)⊤z. The ETO policy is πf^\pi_{\hat f}πf^​​, where f^=f^D\hat f = \hat f_{\mathcal D}f^​=f^​D​ is computed from the data.
  • Regret. Regret(π^)=EDEX[f∗(X)⊤π^D(X)−min⁡z∈Zf∗(X)⊤z]\mathrm{Regret}(\hat\pi) = \mathbb E_{\mathcal D}\mathbb E_X\bigl[f^*(X)^\top \hat\pi_{\mathcal D}(X) - \min_{z \in \mathcal Z} f^*(X)^\top z\bigr]Regret(π^)=ED​EX​[f∗(X)⊤π^D​(X)−minz∈Z​f∗(X)⊤z].

In Lean the polytope is FastCLO.ETO.Polytope, the distribution is FastCLO.ETO.Instance (the law μ\muμ of XXX and the conditional law κ\kappaκ of YYY given XXX), and optVal, Zstar, gap, NoiseCond, IsPlugIn and regret are the objects above.

Formalization targets

Goal: Theorem 8

Fix α,γ≥0\alpha, \gamma \ge 0α,γ≥0 and C1,C2>0C_1, C_2 > 0C1​,C2​>0. There is a constant CCC depending only on these such that, for every instance satisfying the noise condition, every nnn, every an>0a_n > 0an​>0 and every estimator f^\hat ff^​ with

PD(∥f^(x)−f∗(x)∥≥δ)≤C1exp⁡(−C2anδ2)for every δ>0 and almost every x,\mathbb P_{\mathcal D}\bigl(\|\hat f(x) - f^*(x)\| \ge \delta\bigr) \le C_1 \exp(-C_2 a_n \delta^2) \quad \text{for every } \delta > 0 \text{ and almost every } x,PD​(∥f^​(x)−f∗(x)∥≥δ)≤C1​exp(−C2​an​δ2)for every δ>0 and almost every x,

the plug-in policy satisfies

Regret(πf^)≤C B an−(1+α)/2.\mathrm{Regret}(\pi_{\hat f}) \le C\,B\,a_n^{-(1+\alpha)/2}.Regret(πf^​​)≤CBan−(1+α)/2​.

The constant is existential, as on the page; the proof gives an explicit value, but the goal does not commit to it.

Milestones

The milestones are the displayed steps of the paper's proof (Appendix A.5, p. 30), in order:

  1. Pointwise bound. The plug-in decision's excess cost at a point is at most 2B∥f∗(x)−f^(x)∥2B\|f^*(x) - \hat f(x)\|2B∥f∗(x)−f^​(x)∥.
  2. Implication display. If the plug-in decision is strictly suboptimal at xxx and ∥f∗(x)−f^(x)∥≤ε\|f^*(x) - \hat f(x)\| \le \varepsilon∥f∗(x)−f^​(x)∥≤ε, then 0<Δ(x)≤2Bε0 < \Delta(x) \le 2B\varepsilon0<Δ(x)≤2Bε.
  3. Peeled bound. For every δ>0\delta > 0δ>0,
Regret(πf^)≤2Bδ P(0<Δ(X)≤2Bδ)+Bδ∑r≥12r+1 P(∥f∗(X)−f^(X)∥>2r−1δ, 0<Δ(X)≤2r+1Bδ).\mathrm{Regret}(\pi_{\hat f}) \le 2B\delta\,\mathbb P(0 < \Delta(X) \le 2B\delta) + B\delta\sum_{r \ge 1} 2^{r+1}\,\mathbb P\bigl(\|f^*(X) - \hat f(X)\| > 2^{r-1}\delta,\ 0 < \Delta(X) \le 2^{r+1}B\delta\bigr).Regret(πf^​​)≤2BδP(0<Δ(X)≤2Bδ)+Bδr≥1∑​2r+1P(∥f∗(X)−f^​(X)∥>2r−1δ, 0<Δ(X)≤2r+1Bδ).
  1. Bound after iterating expectations. Under the hypotheses of Theorem 8, for every δ>0\delta > 0δ>0,
Regret(πf^)≤Bγα(2δ)α+1+BγαC1δα+1∑r≥12(r+1)(α+1)e−C2an(2r−1δ)2.\mathrm{Regret}(\pi_{\hat f}) \le B\gamma^\alpha(2\delta)^{\alpha+1} + B\gamma^\alpha C_1\delta^{\alpha+1}\sum_{r \ge 1} 2^{(r+1)(\alpha+1)} e^{-C_2 a_n (2^{r-1}\delta)^2}.Regret(πf^​​)≤Bγα(2δ)α+1+BγαC1​δα+1r≥1∑​2(r+1)(α+1)e−C2​an​(2r−1δ)2.

Significance

For well-specified least squares regression, ana_nan​ is of order nnn up to a logarithmic factor (Corollary 4 of the paper), and Theorem 8 gives regret of order n−(1+α)/2n^{-(1+\alpha)/2}n−(1+α)/2 up to logarithms. Theorem 7 of the same paper shows that no algorithm that only uses a policy class of Natarajan dimension η\etaη can beat n−(1+α)/(2+α)n^{-(1+\alpha)/(2+\alpha)}n−(1+α)/(2+α) in the worst case over such instances. Since (1+α)/2>(1+α)/(2+α)(1+\alpha)/2 > (1+\alpha)/(2+\alpha)(1+α)/2>(1+α)/(2+α) for every α≥0\alpha \ge 0α≥0, the two results separate ETO from integrated methods when the regression model is well specified. This runs against the common intuition behind end-to-end learning. The noise condition is the CLO analogue of the Tsybakov margin condition in classification (Tsybakov, 2004): the gap Δ\DeltaΔ plays the role of the distance of the regression function from the decision boundary.

The result is proved in the paper; as far as is known it has not been machine-checked. Formalizing it fixes the exact hypotheses (which quantifiers the tail condition carries, what "universal" means for C1,C2C_1, C_2C1​,C2​, and the measurability the regret integral needs) and yields reusable pieces: the extreme-point argument behind the implication display and the peeling bound in [0,∞][0, \infty][0,∞].

Difficulty

The obvious argument bounds regret by 2B E∥f∗(X)−f^(X)∥2B\,\mathbb E\|f^*(X) - \hat f(X)\|2BE∥f∗(X)−f^​(X)∥, which gives only the rate an−1/2a_n^{-1/2}an−1/2​ and ignores the noise condition (this is Theorem 4 of the paper). The faster rate needs two things the naive bound does not use. First, a wrong decision at xxx is possible only where the gap Δ(x)\Delta(x)Δ(x) is positive and at most a multiple of the estimation error. This uses that plug-in decisions are extreme points, so a strictly suboptimal decision is at least Δ(x)\Delta(x)Δ(x) worse than optimal; it requires finiteness of Z∠\mathcal Z^\angleZ∠, which must be proved for an H-polytope. Second, the error must be controlled pointwise in xxx, with sub-Gaussian tails over the data, and combined with the noise condition shell by shell, through a Fubini exchange between the data and the feature.

Formalization scope

  • Representation. Vectors live in EuclideanSpace ℝ (Fin k). The polytope bundles AAA, bbb, B>0B > 0B>0, nonemptiness and the bound ∥z∥≤B\|z\| \le B∥z∥≤B. The distribution is a probability measure μ\muμ on features and a Markov kernel κ\kappaκ with κ(x){∥y∥>1}=0\kappa(x)\{\|y\| > 1\} = 0κ(x){∥y∥>1}=0; f∗(x)=∫y κ(x)(dy)f^*(x) = \int y\,\kappa(x)(dy)f∗(x)=∫yκ(x)(dy), and the data law is the nnn-fold product of μ⊗κ\mu \otimes \kappaμ⊗κ.
  • Conventions. Optimal values and gaps are real infima, taken over nonempty bounded (indeed finite or compact) sets. The regret is measured against the optimal value, which equals the paper's f∗(X)⊤π∗(X)f^*(X)^\top\pi^*(X)f∗(X)⊤π∗(X). Probabilities are [0,∞][0, \infty][0,∞]-valued measures; (γδ/B)α(\gamma\delta/B)^\alpha(γδ/B)α is a real power with 00=10^0 = 100=1.
  • Disclosed additions. Plug-in policies are any selection from the argmin with values in Z∠\mathcal Z^\angleZ∠ (the paper's consistent tie-breaking is dropped, which strengthens the upper bound). The estimator and the policy are jointly measurable in the data and the feature. C1,C2>0C_1, C_2 > 0C1​,C2​>0 and an>0a_n > 0an​>0. The constant CCC is chosen after α,γ,C1,C2\alpha, \gamma, C_1, C_2α,γ,C1​,C2​ and before the instance, the sample size, the estimator and ana_nan​.
  • Ruled out. The tail hypothesis bounds the probability over the data for each fixed xxx, for every δ>0\delta > 0δ>0; it is not a probability over XXX and not a single δ\deltaδ. The conclusion is uniform: one constant for all instances, estimators and ana_nan​. The peeled bound is stated in [0,∞][0, \infty][0,∞], because a divergent real series would sum to 000.
  • Infrastructure needed. Finiteness of the extreme points of a bounded H-polytope and attainment of linear minima at extreme points; measurability of x↦Δ(x)x \mapsto \Delta(x)x↦Δ(x); Fubini–Tonelli on the product of the data law and μ\muμ; summability of ∑r2(r+1)(α+1)e−c 4r\sum_r 2^{(r+1)(\alpha+1)} e^{-c\,4^r}∑r​2(r+1)(α+1)e−c4r. The polytope facts are reusable well beyond this mission.

Selected references

  • Y. Hu, N. Kallus, X. Mao, Fast Rates for Contextual Linear Optimization, arXiv:2011.03030v3, 2021; Management Science, 2022. https://arxiv.org/abs/2011.03030v3
  • A. N. Elmachtoub, P. Grigas, Smart "Predict, then Optimize", Management Science 68(1), 2022. https://doi.org/10.1287/mnsc.2020.3922
  • A. B. Tsybakov, Optimal aggregation of classifiers in statistical learning, Annals of Statistics 32(1), 2004. https://doi.org/10.1214/aos/1079120131
  • P. Massart, É. Nédélec, Risk bounds for statistical learning, Annals of Statistics 34(5), 2006. https://doi.org/10.1214/009053606000000786
9 thms3 active usersReviewed
🏆Completed
Convex OptimizationOptimization·Captain: mikedeng1

A Modern Introduction to Online Learning V: Online Mirror Descent with a λ-Strongly Convex Regularizer Has Regret max_t B_ψ(u; x_t)/η_T + Σ_t η_t‖g_t‖⋆²/(2λ)Textbook

Motivation

Online learning studies repeated decisions made before their losses are known. A learner chooses a point, observes a loss function, and then uses that observation to choose the next point. The natural performance measure is regret: the learner's cumulative loss minus the cumulative loss of one fixed competitor chosen after the sequence is known. A regret bound makes a claim for every loss sequence, including a sequence adapted to the learner's earlier decisions. Orabona's Chapter 6 develops online mirror descent as a way to obtain such bounds in geometries other than Euclidean distance Orabona, Chapter 6.

The geometry matters because the decision set and its loss gradients can have very different scales under different norms. On a simplex, for example, a Euclidean description treats coordinate changes differently from an entropy-based description. Theorem 6.11 expresses the guarantee in a primal norm chosen for the regularizer and the corresponding dual norm for the observed subgradients. It is the chapter's general regret statement for the basic mirror-descent algorithm Orabona, Theorem 6.11, p. 67.

Setting

Let EEE be a finite-dimensional real normed space, representing the book's Rd\mathbb R^dRd with a chosen norm ∥⋅∥\|\cdot\|∥⋅∥. The learner chooses from a nonempty closed convex set V⊆EV\subseteq EV⊆E. A regularizer ψ\psiψ has domain X⊇VX\supseteq VX⊇V; it is strictly convex on XXX and differentiable on the interior of XXX. Its Bregman divergence is

Bψ(y;x)=ψ(y)−ψ(x)−⟨∇ψ(x),y−x⟩,B_\psi(y;x)=\psi(y)-\psi(x)-\langle\nabla\psi(x),y-x\rangle,Bψ​(y;x)=ψ(y)−ψ(x)−⟨∇ψ(x),y−x⟩,

defined when y∈Xy\in Xy∈X and x∈int⁡Xx\in\operatorname{int}Xx∈intX. It need not be symmetric or satisfy the triangle inequality. The regularizer is λ\lambdaλ-strongly convex on VVV with respect to the chosen norm, where λ>0\lambda>0λ>0. The norm of a gradient ggg, written ∥g∥∗\|g\|_*∥g∥∗​, is the norm dual to ∥⋅∥\|\cdot\|∥⋅∥. These are the conventions of Definition 6.4 and Lemma 6.10 Orabona, pp. 63–67.

At round ttt, beginning with t=1t=1t=1, the learner has xt∈Vx_t\in Vxt​∈V, then observes a loss ℓt\ell_tℓt​ and a selected subgradient gtg_tgt​ at xtx_txt​. With step size ηt>0\eta_t>0ηt​>0, Algorithm 6.1 chooses any minimizer

xt+1∈arg⁡min⁡y∈V{⟨gt,y⟩+1ηtBψ(y;xt)}.x_{t+1}\in\arg\min_{y\in V}\left\{\langle g_t,y\rangle+\frac{1}{\eta_t}B_\psi(y;x_t)\right\}.xt+1​∈argy∈Vmin​{⟨gt​,y⟩+ηt​1​Bψ​(y;xt​)}.

The book gives two ways to keep each update in the interior of XXX. The statements assume either of them, (6.5) or (6.6), V⊆int⁡XV\subseteq\operatorname{int}XV⊆intX, so every iterate is a valid reference point for the next Bregman divergence Orabona, Algorithm 6.1 and (6.6), p. 65. Regret is measured against one fixed u∈Vu\in Vu∈V:

RT(u)=∑t=1T(ℓt(xt)−ℓt(u)).R_T(u)=\sum_{t=1}^{T}\bigl(\ell_t(x_t)-\ell_t(u)\bigr).RT​(u)=t=1∑T​(ℓt​(xt​)−ℓt​(u)).

Formalization targets

For T≥1T\ge1T≥1 and nonincreasing positive step sizes, Theorem 6.11 gives the chapter's general bound

RT(u)≤max⁡1≤t≤TBψ(u;xt)ηT+12λ∑t=1Tηt∥gt∥∗2(u∈V).R_T(u)\le\frac{\max_{1\le t\le T}B_\psi(u;x_t)}{\eta_T} +\frac{1}{2\lambda}\sum_{t=1}^{T}\eta_t\|g_t\|_*^2\qquad(u\in V).RT​(u)≤ηT​max1≤t≤T​Bψ​(u;xt​)​+2λ1​t=1∑T​ηt​∥gt​∥∗2​(u∈V).

When every step equals the same η>0\eta>0η>0, it also gives the more specific bound

RT(u)≤Bψ(u;x1)η+η2λ∑t=1T∥gt∥∗2.R_T(u)\le\frac{B_\psi(u;x_1)}{\eta} +\frac{\eta}{2\lambda}\sum_{t=1}^{T}\|g_t\|_*^2.RT​(u)≤ηBψ​(u;x1​)​+2λη​t=1∑T​∥gt​∥∗2​.

The decreasing-step bound is the goal; the constant-step conclusion is a companion theorem. The milestone list records the book's three-point identity (Lemma 6.7), the unique-minimizer theorem for a closed strongly convex function (Theorem 6.9), and the existence and full three-part inequality of the one-step mirror-descent lemma (Lemma 6.10). Theorem 6.11 follows these results in the source Orabona, pp. 64–67.

Significance

The first bound makes the cost of the update explicit: one term measures the divergence from the competitor to the iterates, and the other accumulates squared dual norms of the selected subgradients. The constant-step form depends only on the initial divergence, which can be substantially easier to estimate. These forms let a later analysis choose a regularizer and norm pair suited to its feasible set and gradients, as the chapter does in its subsequent examples Orabona, §§6.6–6.7.

The theorem and its proof are known mathematically from the textbook. The work here is to express the algorithm, its admissible oracle choices, and the exact bounds as Lean statements and then prove them. A formal proof of the one-step inequality would also be reusable for later dynamic-regret and optimistic mirror-descent results. The published Bregman-divergence definition is reused, but the statements in this mission remain proof targets; a related published three-point identity assumes continuous differentiability, stronger than the book's hypothesis.

Difficulty

The usual squared Euclidean distance has a symmetric algebraic expansion. Bregman divergence is asymmetric and its second argument must lie where the regularizer is differentiable. Thus a direct replacement of squared distance by BψB_\psiBψ​ does not by itself establish that every update is defined or that the terms in a regret sum have the right orientation. Theorem 6.11 also allows changing step sizes, so the divergence terms carry different denominators from one round to the next. The one-step lemma packages the exact inequality needed to make those issues precise Orabona, Lemma 6.10 and Theorem 6.11, p. 67.

Formalization scope

The Lean development uses an arbitrary finite-dimensional real normed space; continuous linear functionals represent subgradients, and their operator norm is ∥⋅∥∗\|\cdot\|_*∥⋅∥∗​. The regularizer is a total real-valued Lean function, but only its values on XXX represent the book's function. Closedness is expressed by the closed epigraph of the extension equal to +∞+\infty+∞ outside XXX. Bregman divergence uses the published definition with the derivative at the second argument. The run predicate records the initial point, a selected subgradient, and a minimizer over VVV on each round; it imposes no tie-breaking rule. The loss functions are real valued on VVV, and subgradient inequalities are stated relative to VVV, the only points compared in the regret result.

All sums use rounds 1,…,T1,\ldots,T1,…,T; index zero carries no data. The maximum is a finite maximum with the explicit condition T≥1T\ge1T≥1. Positive ηt\eta_tηt​ and λ\lambdaλ protect every division in the displayed bounds. Every drafted statement assumes the book's "(6.5) or (6.6)", with (6.5) read as a one-sided limit r→0+r\to0^+r→0+. The update is required to minimize over VVV itself, and uuu ranges over VVV; allowing points outside that set would change the theorem. The minimizer result, the full one-step chain, and both regret bounds are welcome proof targets. The domain-restricted closed-function representation and the finite-round run predicate can be reused in related online convex optimization developments.

Selected references

  • Francesco Orabona, Online Learning: A Modern Introduction Using Convex Optimization, arXiv preprint arXiv:1912.13213v10, 2026. arXiv
7 thms3 active usersReviewed
🏆Completed
ProbabilityStatistics·Captain: mikedeng1

Robustness May Be at Odds with Accuracy 1: Every Classifier with Standard Accuracy ≥ 1 − δ Has ℓ∞ Robust Accuracy ≤ pδ/(1 − p) Against ε ≥ 2ηResearch Paper

Motivation

Classifiers trained to be robust against small, adversarially chosen input perturbations routinely lose standard accuracy. Adversarial training (Madry et al., 2018) improves worst-case accuracy under an ℓp\ell_pℓp​-bounded attack, but on MNIST, CIFAR-10 and ImageNet the resulting models are measurably less accurate on unperturbed inputs, and the loss grows with the attack radius. A natural reading is that this is an artefact of training procedures or of finite data.

Tsipras, Santurkar, Engstrom, Turner and Mądry (ICLR 2019) argue otherwise: they exhibit a simple distribution on which every classifier faces the trade-off, so that the tension is a property of the data distribution itself and survives infinite data, including for the Bayes-optimal classifier. The argument isolates a distinction that has since been influential: robust features, moderately correlated with the label and stable under small perturbations, versus non-robust features, individually weak but jointly highly predictive and easily flipped. This mission formalizes the paper's first theorem, Theorem 2.1, and its proof in Appendix C.

Setting

Fix d∈Nd \in \mathbb Nd∈N, a probability ppp with 12≤p<1\tfrac12 \le p < 121​≤p<1, and η≥0\eta \ge 0η≥0. The distribution D\mathcal DD of an input-label pair (x,y)(x, y)(x,y), x=(x1,x2,…,xd+1)∈Rd+1x = (x_1, x_2, \dots, x_{d+1}) \in \mathbb R^{d+1}x=(x1​,x2​,…,xd+1​)∈Rd+1, is the data model (3):

y∼Unif{−1,+1},x1={+yw.p. p−yw.p. 1−p,x2,…,xd+1∼i.i.d.N(ηy,1),y \sim \mathrm{Unif}\{-1,+1\}, \qquad x_1 = \begin{cases} +y & \text{w.p. } p\\ -y & \text{w.p. } 1-p\end{cases}, \qquad x_2, \dots, x_{d+1} \overset{\text{i.i.d.}}{\sim} \mathcal N(\eta y, 1),y∼Unif{−1,+1},x1​={+y−y​w.p. pw.p. 1−p​,x2​,…,xd+1​∼i.i.d.N(ηy,1),

with x1x_1x1​ and (x2,…,xd+1)(x_2, \dots, x_{d+1})(x2​,…,xd+1​) independent given yyy. Write G+G_+G+​ and G−G_-G−​ for the law N(±η,1)⊗d\mathcal N(\pm\eta, 1)^{\otimes d}N(±η,1)⊗d of the weak features given y=+1y = +1y=+1 and y=−1y = -1y=−1.

A classifier is a measurable map f:Rd+1→{−1,+1}f : \mathbb R^{d+1} \to \{-1, +1\}f:Rd+1→{−1,+1}. Its standard accuracy is Pr⁡(x,y)∼D[f(x)=y]\Pr_{(x,y)\sim\mathcal D}[f(x) = y]Pr(x,y)∼D​[f(x)=y], the objective (1) of the paper with the 0–1 loss. Its ℓ∞\ell_\inftyℓ∞​ robust accuracy at radius ε\varepsilonε is

accε(f)=Pr⁡(x,y)∼D[f(x+δ)=y for every δ∈Rd+1 with ∥δ∥∞≤ε],\mathrm{acc}_\varepsilon(f) = \Pr_{(x,y)\sim\mathcal D}\big[f(x + \delta) = y \text{ for every } \delta \in \mathbb R^{d+1} \text{ with } \|\delta\|_\infty \le \varepsilon\big],accε​(f)=(x,y)∼DPr​[f(x+δ)=y for every δ∈Rd+1 with ∥δ∥∞​≤ε],

the objective (2) with the 0–1 loss and perturbation set Δ={∥δ∥∞≤ε}\Delta = \{\|\delta\|_\infty \le \varepsilon\}Δ={∥δ∥∞​≤ε}. In Lean these are stdAcc d p η f and robustAcc d p η f ε, each an average with weight 12\tfrac1221​ of the two class-conditional probabilities. The proof also uses pijp_{ij}pij​, the probability that fff predicts +1+1+1 when x1=ix_1 = ix1​=i and the weak features follow GjG_jGj​ (pij), and the shift adversary xi↦xi−2ηyx_i \mapsto x_i - 2\eta yxi​↦xi​−2ηy for i≥2i \ge 2i≥2 (shiftAdv, accuracy shiftAcc).

Formalization targets

Goal: Theorem 2.1 (p. 5)

For every measurable classifier fff and every δ∈R\delta \in \mathbb Rδ∈R,

Pr⁡[f(x)=y]≥1−δ ⟹ accε(f)≤p1−p δfor every ε≥2η.\Pr[f(x) = y] \ge 1 - \delta \ \Longrightarrow\ \mathrm{acc}_\varepsilon(f) \le \frac{p}{1-p}\,\delta \qquad \text{for every } \varepsilon \ge 2\eta.Pr[f(x)=y]≥1−δ ⟹ accε​(f)≤1−pp​δfor every ε≥2η.

The statement quantifies over all classifiers, not a parametric family, and over every radius ε≥2η\varepsilon \ge 2\etaε≥2η.

Milestones (App. C, pp. 14–16)

  1. The shift zi↦zi−2ηyz_i \mapsto z_i - 2\eta yzi​↦zi​−2ηy pushes N(ηy,1)⊗d\mathcal N(\eta y, 1)^{\otimes d}N(ηy,1)⊗d forward to N(−ηy,1)⊗d\mathcal N(-\eta y, 1)^{\otimes d}N(−ηy,1)⊗d: it changes G+G_+G+​ to G−G_-G−​ and vice versa.
  2. For ε≥2η≥0\varepsilon \ge 2\eta \ge 0ε≥2η≥0, robust accuracy is at most the accuracy against the shift adversary.
  3. Pr⁡[f(x)=y]=12(p(1+p++−p−−)+(1−p)(1+p−+−p+−))\Pr[f(x) = y] = \tfrac12\big(p(1 + p_{++} - p_{--}) + (1-p)(1 + p_{-+} - p_{+-})\big)Pr[f(x)=y]=21​(p(1+p++​−p−−​)+(1−p)(1+p−+​−p+−​)).
  4. Pr⁡[f(xadv)=y]=12(p(1+p+−−p−+)+(1−p)(1+p−−−p++))\Pr[f(x_{\mathrm{adv}}) = y] = \tfrac12\big(p(1 + p_{+-} - p_{-+}) + (1-p)(1 + p_{--} - p_{++})\big)Pr[f(xadv​)=y]=21​(p(1+p+−​−p−+​)+(1−p)(1+p−−​−p++​)).
  5. With a=1−p+++p−−≥0a = 1 - p_{++} + p_{--} \ge 0a=1−p++​+p−−​≥0 and b=1−p−++p+−b = 1 - p_{-+} + p_{+-}b=1−p−+​+p+−​: 1−12(pa+(1−p)b)≥1−δ1 - \tfrac12(pa + (1-p)b) \ge 1 - \delta1−21​(pa+(1−p)b)≥1−δ implies 12((1−p)a+pb)≤p1−pδ\tfrac12((1-p)a + pb) \le \tfrac{p}{1-p}\delta21​((1−p)a+pb)≤1−pp​δ.

Companion statements (§2.1, pp. 4–5)

  • Tightness. sign(x1)\mathrm{sign}(x_1)sign(x1​) has standard accuracy ppp and robust accuracy ppp for 0≤ε<10 \le \varepsilon < 10≤ε<1; at δ=1−p\delta = 1-pδ=1−p the bound of Theorem 2.1 is attained.
  • Standard classification is easy. For favg(x)=sign(1d∑i≥2xi)f_{\mathrm{avg}}(x) = \mathrm{sign}(\tfrac1d\sum_{i\ge2} x_i)favg​(x)=sign(d1​∑i≥2​xi​), Pr⁡[favg(x)=y]=Pr⁡[N(η,1/d)>0]\Pr[f_{\mathrm{avg}}(x) = y] = \Pr[\mathcal N(\eta, 1/d) > 0]Pr[favg​(x)=y]=Pr[N(η,1/d)>0], which exceeds 99%99\%99% when η≥3/d\eta \ge 3/\sqrt dη≥3/d​.

Significance

Theorem 2.1 converts the empirical observation into a statement about a distribution: on D\mathcal DD, standard accuracy 1−δ1 - \delta1−δ caps robust accuracy at p1−pδ\frac{p}{1-p}\delta1−pp​δ, so as δ→0\delta \to 0δ→0 robust accuracy falls to 000, while the adversary's budget ε=2η=Θ(1/d)\varepsilon = 2\eta = \Theta(1/\sqrt d)ε=2η=Θ(1/d​) is small compared with the feature scale ±1\pm 1±1. Together with the companion statements it shows that the model admits a classifier with accuracy above 99%99\%99% and a classifier with robust accuracy ppp, but no classifier with both. It underlies the paper's conclusion that robust training must give up non-robust features, and the later line of work on robust versus non-robust features.

The result is proved in the paper, with an elementary argument. As far as the platform record and a search of the literature show, it has not been machine-checked. The formalization adds an exact, convention-explicit version of the model and of robust accuracy (worst case over a closed ℓ∞\ell_\inftyℓ∞​ ball, outer measure of a possibly non-measurable event), a Lean statement of the shift-invariance property of Gaussian products that the argument rests on, and a reusable template for 0–1-loss adversarial objectives that later missions on robustness in Gaussian models can follow.

Difficulty

The obvious attempt is to bound robust accuracy for a specific classifier, or for linear classifiers, where the adversary's best response is explicit. That does not reach the theorem: the claim is about every measurable fff, for which no best response is computable and the robust event need not even be measurable. On the formal side, the statements involve the pushforward of a product of Gaussians under a coordinatewise shift, the decomposition of the class-conditional law into the two-point law of x1x_1x1​ and the product G±G_\pmG±​, and a measure applied to a possibly non-measurable robust event; none of these is available as a ready-made lemma for this model.

Formalization scope

Inputs live in Fin (d + 1) → ℝ; coordinate 0 is x1x_1x1​ and coordinate i.succ is xi+2x_{i+2}xi+2​. The Mathlib norm on this type is the sup norm, so ‖δ‖ ≤ ε is exactly the ℓ∞\ell_\inftyℓ∞​ ball. Labels and predictions are the reals ±1\pm1±1. The class-conditional law is ((x1Law p y).prod (weakLaw d η y)).map (Matrix.vecCons ·.1 ·.2) with weakLaw d η y = Measure.pi (fun _ => gaussianReal (y * η) 1). Accuracies are real numbers in [0,1][0,1][0,1] computed with toReal. The robust event is measured by its outer measure, i.e. under the completion of D\mathcal DD; no measurability of that event is assumed.

The standing assumption p≥12p \ge \tfrac12p≥21​ of model (3) is a hypothesis of every statement. Four hypotheses implicit on the page are explicit: p<1p < 1p<1 in the goal (the bound divides by 1−p1-p1−p, and Lean's p/0=0p/0 = 0p/0=0 would make the statement false at p=1p = 1p=1); η≥0\eta \ge 0η≥0 (the page takes η\etaη large, Θ(1/d)\Theta(1/\sqrt d)Θ(1/d​); with η<0\eta < 0η<0 the case ε=0\varepsilon = 0ε=0 has no adversary and the statement fails); measurability of fff; and f(x)∈{±1}f(x) \in \{\pm1\}f(x)∈{±1} for every xxx.

A trivializing formalization is ruled out: the goal is stated for every measurable {±1}\{\pm1\}{±1}-valued classifier and for the worst case over the full ℓ∞\ell_\inftyℓ∞​ ball of every radius ε≥2η\varepsilon \ge 2\etaε≥2η, not for linear classifiers, not for the fixed shift adversary, and not with the pijp_{ij}pij​ among its hypotheses; those appear only in milestones.

The development needs Mathlib's gaussianReal, Measure.pi, product measures and pushforwards. Welcome contributions: the Gaussian shift identity (reusable for any Gaussian location model), the disintegration of condLaw along x1x_1x1​, and the two accuracy formulas.

Selected references

  • D. Tsipras, S. Santurkar, L. Engstrom, A. Turner, A. Mądry, Robustness May Be at Odds with Accuracy, ICLR 2019; arXiv:1805.12152v5 (2019). https://arxiv.org/abs/1805.12152
  • A. Madry, A. Makelov, L. Schmidt, D. Tsipras, A. Vladu, Towards Deep Learning Models Resistant to Adversarial Attacks, ICLR 2018. https://arxiv.org/abs/1706.06083
  • L. Schmidt, S. Santurkar, D. Tsipras, K. Talwar, A. Mądry, Adversarially Robust Generalization Requires More Data, NeurIPS 2018. https://arxiv.org/abs/1804.11285
7 thms3 active usersReviewed
Linear algebraNumerical AnalysisProbability+2·Captain: mikedeng1

An Introduction to Matrix Concentration Inequalities III: A Matrix Sampling Estimator Satisfies E‖R̄ₙ − B‖ ≤ √(2m₂(R)log(d₁ + d₂)/n) + 2L log(d₁ + d₂)/(3n)Textbook

Motivation

Matrix approximation replaces a large or complicated matrix with a random matrix that is easier to store or compute. Averaging independent realizations improves accuracy while retaining the structure of each realization. Examples in Tropp's monograph include sparse approximations and randomized matrix products. The practical question is how the number of samples, the size of one sample, and its second moment control the spectral norm error. Tropp, §6.2

The book proves a matrix Bernstein inequality in Theorem 6.1.1, then applies it to an empirical matrix average in Corollary 6.2.1. Theorem 6.1.1 is already posed as TroppMatrixConcentration.matrix_bernstein on Prove2Me and is currently Open. This mission records the application with its own hypotheses and constants. The capstone concerns the estimator that a reader would construct, rather than a sum of abstract centered summands. Tropp, §§6.1–6.2

Setting

Let BBB be a fixed complex d1×d2d_1\times d_2d1​×d2​ matrix. A sample matrix RRR has the same size, has expectation BBB, and has spectral norm at most LLL almost surely. The equality ER=B\mathbb ER=BER=B makes one sample an unbiased estimate of BBB. The almost-sure bound limits the size of every sample except on a probability-zero set. The dimensions d1,d2d_1,d_2d1​,d2​ are positive integers, and nnn is a positive sample count.

Take independent copies R1,…,RnR_1,\ldots,R_nR1​,…,Rn​ with the distribution of RRR. Their matrix sampling estimator is

R‾n=1n∑k=1nRk.\overline R_n=\frac1n\sum_{k=1}^nR_k.Rn​=n1​k=1∑n​Rk​.

For any rectangular matrix AAA, A∗A^*A∗ denotes its conjugate transpose and ∥A∥\|A\|∥A∥ its spectral norm. The per-sample second moment is

m2(R)=max⁡{∥E(RR∗)∥,∥E(R∗R)∥}.m_2(R)=\max\{\|\mathbb E(RR^*)\|,\|\mathbb E(R^*R)\|\}.m2​(R)=max{∥E(RR∗)∥,∥E(R∗R)∥}.

Both terms matter for a rectangular matrix: RR∗RR^*RR∗ and R∗RR^*RR∗R act on different coordinate spaces. All expectations of matrices are entrywise expectations, represented by Bochner integrals in Lean. The almost-sure norm bound guarantees the moments used here are finite. Tropp, equations (6.2.4)–(6.2.6)

Formalization targets

Empirical matrix approximation

Corollary 6.2.1 gives an expectation estimate and a tail estimate for the same estimator. For every real t≥0t\ge0t≥0,

E∥R‾n−B∥≤2m2(R)log⁡(d1+d2)n+2Llog⁡(d1+d2)3n,\mathbb E\|\overline R_n-B\|\le \sqrt{\frac{2m_2(R)\log(d_1+d_2)}{n}}+ \frac{2L\log(d_1+d_2)}{3n},E∥Rn​−B∥≤n2m2​(R)log(d1​+d2​)​​+3n2Llog(d1​+d2​)​, P{∥R‾n−B∥≥t}≤(d1+d2)exp⁡(−nt2/2m2(R)+2Lt/3).\mathbb P\{\|\overline R_n-B\|\ge t\}\le (d_1+d_2)\exp\left(\frac{-nt^2/2}{m_2(R)+2Lt/3}\right).P{∥Rn​−B∥≥t}≤(d1​+d2​)exp(m2​(R)+2Lt/3−nt2/2​).

The goal includes both bounds. The expectation estimate controls average error; the tail estimate controls the chance of an error above a chosen threshold. The dimension enters only through log⁡(d1+d2)\log(d_1+d_2)log(d1​+d2​) in the expectation bound and the factor d1+d2d_1+d_2d1​+d2​ in the tail bound. These are the book's constants, including the factors 222 and 333 in the linear terms. Tropp, Corollary 6.2.1, p. 81

Supporting statements

The milestone list follows Corollary 6.1.2 for uncentered summands, the bound ∥n−1(Rk−B)∥≤2L/n\|n^{-1}(R_k-B)\|\le2L/n∥n−1(Rk​−B)∥≤2L/n, and the variance estimate v(R‾n−B)≤m2(R)/nv(\overline R_n-B)\le m_2(R)/nv(Rn​−B)≤m2​(R)/n. Equation (6.2.7) is a companion: a stated lower bound on nnn implies E∥R‾n−B∥≤2ε\mathbb E\|\overline R_n-B\|\le2\varepsilonE∥Rn​−B∥≤2ε for ε>0\varepsilon>0ε>0. Tropp, pp. 77, 82

Significance

The result turns a bound on one random matrix into a quantitative guarantee for an empirical approximation. It lets a user choose a sample count for a target error and evaluate both average error and failure probability. The same framework is instantiated later in the monograph for matrix sparsification and multiplication; those applications depend on controlling LLL and m2(R)m_2(R)m2​(R) for their particular sampling laws. Tropp, §§6.2–6.4

This is a known mathematical result proved in the book. Its Lean statement is an open formalization target. The mission adds a reusable empirical-mean definition and application statements built on the already published matrix spectral norm and rectangular second-moment definitions. A proof of the general matrix Bernstein theorem remains separate; the mission can draw on it when that platform theorem is proved.

Difficulty

The estimator error is a centered sum, but the sampling assumptions are given for an uncentered random matrix RRR. Bounding the centered summands changes the uniform bound from LLL to 2L/n2L/n2L/n. The variance parameter must be related to the second moment of one uncentered sample, with the correct 1/n1/n1/n scaling. Dropping the common-distribution condition would leave no relation between the samples' second moments and m2(R)m_2(R)m2​(R). Applying scalar concentration to entries would not directly give the stated spectral norm bound or its dimension dependence. Tropp, proof of Corollary 6.2.1, p. 82

Formalization scope

Matrices are Matrix (Fin d₁) (Fin d₂) ℂ. The spectral norm and rectSecondMoment are imported from tc's published TroppMatrixConcentration definitions. An arbitrary probability space carries a measurable reference matrix R0R_0R0​ and a finite family RkR_kRk​; iIndepFun expresses independence and IdentDistrib expresses equality in law with R0R_0R0​. The estimator is the finite sum scaled by n−1n^{-1}n−1. No particular sampling scheme or deterministic choice of RRR is fixed: the formal goal requires only the source's unbiasedness and boundedness conditions.

The Lean statements explicitly require d1,d2≥1d_1,d_2\ge1d1​,d2​≥1, n>0n>0n>0, measurability, and the almost-sure form of ∥R∥≤L\|R\|\le L∥R∥≤L. These make the logarithm and inverse sample count meaningful and make the random matrix laws well behaved. For the uncentered Bernstein corollary, integrability of each SkS_kSk​ is explicit because a bounded deviation from an unspecified mean alone does not establish that the mean exists. The goal also asserts integrability of its error norm, preventing a Bochner integral from silently taking Lean's default value for a nonintegrable function. The bound on RRR supplies that integrability. When the denominator in a tail formula is zero, the displayed total real division gives a loose bound at the degenerate case; the source's probabilistic claim remains valid.

Contributions that formalize the centered-summand calculation, the variance estimate, or the transfer from the general matrix Bernstein theorem are within scope. The imported matrix and probability definitions are reusable in other matrix concentration missions. The local empirical-mean definition is reusable whenever a matrix estimator averages finitely many random samples.

Selected references

  • Joel A. Tropp, An Introduction to Matrix Concentration Inequalities, Foundations and Trends in Machine Learning 8(1–2), 2015, arXiv:1501.01571v1, especially §§6.1–6.4.
7 thms3 active usersReviewed
🏆Completed
Algorithmic Game Theory·Captain: mikedeng1

Intrinsic Robustness of the Price of Anarchy II: Outcome Sequences with Vanishing Average Regret Have Average Cost at Most the Robust Price of Anarchy Times OptimalResearch Paper

Motivation

In a game, each participant chooses a strategy while the cost of that choice depends on what everyone else does. A network user may choose a route, for example, but its travel time depends on the routes chosen by other users. The price of anarchy compares the cost of self-interested outcomes with the best joint outcome. A guarantee restricted to a Nash equilibrium leaves open what happens when players repeatedly adapt their choices without settling on any equilibrium. Roughgarden's extension theorem addresses sequences produced by players whose accumulated costs are asymptotically competitive with fixed alternatives. Its bound depends on a property of the game, rather than on a particular learning algorithm or on convergence of play. Roughgarden, §§1 and 3.2.

Setting

There is a finite set of players iii, and player iii has a strategy set SiS_iSi​. These sets may be infinite. An outcome sss specifies one strategy si∈Sis_i\in S_isi​∈Si​ for every player. The real number Ci(s)C_i(s)Ci​(s) is player iii's cost, and the joint cost is C(s)=∑iCi(s)C(s)=\sum_i C_i(s)C(s)=∑i​Ci​(s). The paper treats this objective as nonnegative. If player iii alone changes to the strategy si′s'_isi′​ of another outcome s′s's′, its cost is written Ci(si′,s−i)C_i(s'_i,s_{-i})Ci​(si′​,s−i​); every other player's strategy stays as in sss. Roughgarden, §2.1.

A game is (λ,μ)(\lambda,\mu)(λ,μ)-smooth if every pair s,s′s,s's,s′ satisfies ∑iCi(si′,s−i)≤λC(s′)+μC(s)\sum_i C_i(s'_i,s_{-i})\leq\lambda C(s')+\mu C(s)∑i​Ci​(si′​,s−i​)≤λC(s′)+μC(s). The robust price of anarchy ρ(G)\rho(G)ρ(G) is the infimum of λ/(1−μ)\lambda/(1-\mu)λ/(1−μ) over all such smoothness pairs with μ<1\mu<1μ<1. This number packages every bound obtainable from the same smoothness condition. It is taken in the extended reals, with value +∞+\infty+∞ when no pair exists; the mission assumes an admissible pair exists, so it is below +∞+\infty+∞. Roughgarden, Definitions 2.1–2.2.

For a played outcome sts^tst at time ttt and a comparison outcome s′s's′, define the hypothetical improvement δi(st;s′)=Ci(st)−Ci(si′,s−it)\delta_i(s^t;s')=C_i(s^t)-C_i(s'_i,s^t_{-i})δi​(st;s′)=Ci​(st)−Ci​(si′​,s−it​). A player has vanishing average external regret when its average cost, up to an error tending to zero, is no greater than the average cost it would have paid by using any one fixed strategy throughout the first TTT rounds. The comparison strategy is held fixed across time, while the other players' strategies may change on every round. Each player may have its own error function. Roughgarden, (22), (24).

Formalization targets

Additive bound

The paragraph following equation (25) asserts a bound with a vanishing additive error. For any comparison outcome s′s's′ and every horizon T≥1T\geq1T≥1, there is a function η(T)→0\eta(T)\to0η(T)→0 with

1T∑t=1TC(st)≤ρ(G)C(s′)+η(T).\frac1T\sum_{t=1}^{T}C(s^t)\leq\rho(G)C(s')+\eta(T).T1​t=1∑T​C(st)≤ρ(G)C(s′)+η(T).

This target includes the case C(s′)=0C(s')=0C(s′)=0. The paper states the conclusion for an optimal outcome; the same displayed relations also apply to any comparison outcome. Relations (23) and (25) are separate milestones: the former bounds one outcome using the players' hypothetical improvements, and the latter averages that bound over time. Roughgarden, §3.2, pp. 13–14.

Multiplicative repeated-play bound

The mission goal is Theorem 3.3's multiplicative presentation. For a comparison outcome with C(s′)>0C(s')>0C(s′)>0, there is a function η(T)→0\eta(T)\to0η(T)→0 such that

1T∑t=1TC(st)≤(ρ(G)+η(T))C(s′)(T≥1).\frac1T\sum_{t=1}^{T}C(s^t)\leq\bigl(\rho(G)+\eta(T)\bigr)C(s')\qquad(T\geq1).T1​t=1∑T​C(st)≤(ρ(G)+η(T))C(s′)(T≥1).

The positive-cost condition is needed to make the printed statement true. Consider one player with strategies of costs 111 and 000, playing the costly strategy once and the zero-cost strategy thereafter. Its average regret tends to zero, but its average cost is 1/T>01/T>01/T>0 at every finite horizon. Taking the zero-cost strategy as s′s's′ makes the right-hand side of the printed multiplicative bound zero for any η\etaη. The additive milestone preserves the valid zero-cost case. Roughgarden, Theorem 3.3 and the paragraph preceding it.

Significance

The result transfers a bound defined by the game's smoothness condition to repeated play even when the sequence does not converge. It therefore gives a welfare guarantee for the behavior captured by vanishing external regret without identifying which learning rule generated the sequence. The additive statement also records exactly what remains true when an optimal outcome has zero cost. Roughgarden, §3.2.

Formalizing the result makes the quantifiers behind o(1)o(1)o(1) explicit: for each fixed game, sequence, and comparison outcome, one error sequence works across all positive horizons and tends to zero. The development also distinguishes the paper's valid additive consequence from its multiplicative display at zero cost. The paper proves the result mathematically; the Lean files in this mission are theorem statements for solvers and do not yet contain machine-checked proofs.

Difficulty

For a fixed smoothness pair, a finite-time inequality controls average joint cost by the comparison cost and players' average hypothetical improvements. Regret controls those improvements asymptotically, but the robust price of anarchy takes an infimum over all admissible pairs. Passing from bounds for individual pairs to that infimum requires the sign and boundedness conditions on the joint cost to be tracked correctly. A direct replacement of the fixed pair's ratio by ρ(G)\rho(G)ρ(G) does not itself justify the limit. The zero-cost example also shows why an additive error cannot simply be displayed as a multiplicative o(1)o(1)o(1) term in every case. Roughgarden, (23)–(25) and Theorem 3.3.

Formalization scope

Lean uses a finite type of players, arbitrary strategy types, and real-valued player costs. It defines the joint objective by the finite sum ∑iCi(s)\sum_i C_i(s)∑i​Ci​(s); it does not assume that each individual cost is nonnegative. The goal and additive milestone assume that the joint objective is nonnegative at every outcome, matching the paper's standing convention. They also assume at least one smoothness pair with μ<1\mu<1μ<1, the case in which ρ(G)=+∞\rho(G)=+\inftyρ(G)=+∞ and the page's bound is vacuous. ρ(G)\rho(G)ρ(G) is an extended-real infimum and the inequalities are compared in the extended reals. For a nonnegative game with an admissible pair, ρ(G)\rho(G)ρ(G) is a real number at least 1 when some outcome has positive cost; if all outcomes cost zero, the target cost bounds are zero. Roughgarden, §§1, 2.1, 3.2.

Time indices are 1,…,T1,\ldots,T1,…,T, with T≥1T\geq1T≥1. The Lean sequence is infinite so that the asymptotic assertion concerns prefixes of one sequence; its value at index zero is unused. The regret definition quantifies over every fixed alternative strategy rather than assuming a minimum is attained when a strategy set is infinite. A separate vanishing error sequence may be chosen for each player. The conclusion chooses one vanishing error sequence after the game, play sequence, and comparison outcome, and uses it at every horizon. These choices make the o(1)o(1)o(1) terms genuine limits, rather than unconstrained values chosen separately for each TTT. Roughgarden, (24), Theorem 3.3.

The game objective, smoothness predicate, robust price of anarchy, deviation improvement, and regret predicate are the reusable definitions. Solvers can contribute proofs of the one-step bound, its averaged form, the additive conclusion, and the main theorem. Mixed-outcome sequences from Remark 3.4 are outside this mission's core; their probability and integrability conventions would require another definition layer.

Selected references

  • Tim Roughgarden, Intrinsic Robustness of the Price of Anarchy, Journal of the ACM 62(5), Article 32, 2015. DOI: 10.1145/2806883. Statements and page numbers follow the author's final version dated July 14, 2015.
7 thms3 active usersReviewed
🏆Completed
ProbabilityStatistics·Captain: mikedeng1

A Modern Introduction to Online Learning XV: Online Learning Bounds Rademacher Complexity, Giving Massart's Lemma G∞√(2 ln d)/√TTextbook

Why random signs matter in learning

Learning theory asks whether a rule that fits an observed sample will also perform well on future observations. A class with many functions can fit accidental patterns in the sample, so its size and geometry affect the reliability of empirical risk estimates. Empirical Rademacher complexity measures this ability to fit noise: attach an independent random sign to each observed point, let the class choose its best function after seeing all signs, and average the resulting signed score. A smaller value supports tighter uniform generalization bounds. Orabona, §16.1, pp. 270–271

Online learning studies a different protocol. A learner chooses a prediction before each round's loss is revealed, and its regret compares its cumulative loss with that of a fixed competitor chosen in hindsight. Chapter 16 of Orabona's book connects these two quantities. The connection matters because an online algorithm's regret estimate can be used to control the statistical complexity of a class of predictors, even though the online game and the random-sign experiment are defined differently. Orabona, Theorem 16.4, p. 272

Setting and notation

Fix a horizon T≥1T\ge1T≥1 and vectors z1,…,zT∈Rdz_1,\ldots,z_T\in\mathbb R^dz1​,…,zT​∈Rd. For each ttt, let εt\varepsilon_tεt​ be an independent Rademacher sign, equal to 111 or −1-1−1 with probability one half each. Given a nonempty bounded set V⊆RdV\subseteq\mathbb R^dV⊆Rd, define its linear class by

Hlin={z↦⟨x,z⟩:x∈V}.H_{\mathrm{lin}}=\{z\mapsto\langle x,z\rangle:x\in V\}.Hlin​={z↦⟨x,z⟩:x∈V}.

The empirical complexity on the fixed sample S=(z1,…,zT)S=(z_1,\ldots,z_T)S=(z1​,…,zT​) is

R^S(Hlin)=1TEε[sup⁡x∈V∑t=1Tεt⟨x,zt⟩].\widehat{\mathfrak R}_S(H_{\mathrm{lin}}) =\frac1T\mathbb E_\varepsilon\left[ \sup_{x\in V}\sum_{t=1}^T\varepsilon_t\langle x,z_t\rangle \right].RS​(Hlin​)=T1​Eε​[x∈Vsup​t=1∑T​εt​⟨x,zt​⟩].

This is a finite expectation over the 2T2^T2T possible sign assignments. The supremum is over the same comparison set VVV for every assignment; the chosen xxx may depend on the signs. The book's general Definition 16.1 allows an arbitrary class of real-valued functions, while this mission uses its linear specialization from the discussion immediately preceding Theorem 16.4. Orabona, Definition 16.1 and §16.1, pp. 270–272

For the final bound, VVV is the probability simplex Δd−1={x∈Rd:xi≥0, ∑ixi=1}\Delta^{d-1}=\{x\in\mathbb R^d:x_i\ge0,\ \sum_i x_i=1\}Δd−1={x∈Rd:xi​≥0, ∑i​xi​=1}. The supremum of a linear score over this simplex is the maximum of its ddd coordinate scores. Write G∞G_\inftyG∞​ for a common bound on the infinity norm of the sample vectors, so ∣zt,i∣≤G∞|z_{t,i}|\le G_\infty∣zt,i​∣≤G∞​ for every round and coordinate. Orabona, Corollary 16.7, p. 273

Formalization targets

The central target is Corollary 16.7, called Massart's Lemma in the book:

1TEε[max⁡1≤i≤d∑t=1Tεtzt,i]≤G∞2ln⁡dT.\frac1T\mathbb E_\varepsilon\left[ \max_{1\le i\le d}\sum_{t=1}^T\varepsilon_tz_{t,i} \right] \le \frac{G_\infty\sqrt{2\ln d}}{\sqrt T}.T1​Eε​[1≤i≤dmax​t=1∑T​εt​zt,i​]≤T​G∞​2lnd​​.

Its constant and dependence on ddd and TTT are part of the target. The result includes d=1d=1d=1: in that case the expectation and the bound are both zero. Orabona, Corollary 16.7, p. 273

Two statements supply the chapter's route to this target. Theorem 16.4 says that if a nonanticipating online linear algorithm has worst-case regret at most R(T)R(T)R(T) for all loss vectors in a symmetric set ZZZ, then R^S(Hlin)≤R(T)/T\widehat{\mathfrak R}_S(H_{\mathrm{lin}})\le R(T)/TRS​(Hlin​)≤R(T)/T for any sample in ZZZ. The Exponentiated Gradient regret bound of §6.6 gives, on the simplex, ∑t(ℓt(xt)−ℓt(u))≤(ln⁡d)/η+(η/2)∑t∥gt∥∞2\sum_t(\ell_t(x_t)-\ell_t(u))\le(\ln d)/\eta+(\eta/2)\sum_t\|g_t\|_\infty^2∑t​(ℓt​(xt​)−ℓt​(u))≤(lnd)/η+(η/2)∑t​∥gt​∥∞2​ for every simplex competitor uuu. The milestone list quotes both source statements. Orabona, Theorem 16.4, p. 272; §6.6, p. 81

What the result provides

The goal supplies a sample-dependent complexity bound using only a coordinate bound, the number of coordinates, and the sample size. Its square-root dependence on ln⁡d\ln dlnd explains why a simplex-indexed class can have much smaller complexity than a crude bound that scales directly with ddd. In applications where a generalization estimate contains empirical Rademacher complexity, this result gives a concrete quantity to insert for a finite coordinate class. Orabona, §16.1, pp. 271–273

The mathematical statements are proved in the book. This mission asks for machine-checked proofs of the exact finite-sign formulation of the goal, the regret-to-complexity reduction, and the source's Exponentiated Gradient regret estimate. The local proposal currently contains statements with proof placeholders; it does not claim that these three items have already been proved in Lean. The resulting finite-sign and bounded-supremum definitions can also support other empirical-process arguments.

Where the work is difficult

The reduction has to respect the order of the online game. A bound for a learner allowed to inspect the current random sign would not imply the same complexity estimate: it could choose its round-ttt point using information that a genuine online learner lacks. The algorithm therefore needs a precise causality condition, as well as a regret guarantee that holds for every admissible gradient sequence. Orabona, proof of Theorem 16.4, p. 272

The simplex result also depends on the exact regret coefficient. A generic O(T)O(\sqrt T)O(T​) estimate does not recover the displayed 2ln⁡d\sqrt{2\ln d}2lnd​ constant. At the boundaries, real division and logarithms in Lean are total, so the statements must say when TTT and ddd are positive. The one-dimensional case is meaningful and must stay in the theorem.

Formalization scope

Vectors are functions Fin d → ℝ, representing Rd\mathbb R^dRd with its coordinate pairing. Samples, predictions and gradients use rounds 1,…,T1,\ldots,T1,…,T; index zero is unused. The random-sign expectation is an exact average over Boolean assignments, with no measure-theoretic integrability convention. Both Rademacher averages require T≥1T\ge1T≥1. The signed linear supremum is real-valued and receives evidence that VVV is nonempty and bounded. Theorem 16.4 also records the book's closedness hypothesis.

Theorem 16.4 formalizes deterministic online algorithms. The source permits internal randomization as well; that extension requires an independent probability model for the algorithm's random choices. The Exponentiated Gradient item represents losses by real-valued functions and requires the subgradient inequality on the simplex, where the book evaluates them. Its run predicate requires η>0\eta>0η>0 and includes the uniform initial point and coordinate update of Algorithm 6.2. The goal itself has no algorithm hypothesis. Its infinity-norm assumption is expressed as coordinate absolute-value bounds.

The definitions exclude an empty comparison class or an unbounded real supremum being silently assigned a default Lean value. A complete development needs finite-sign averaging identities, properties of real suprema over bounded sets, simplex geometry, and the regret analysis of Exponentiated Gradient. Contributions to these reusable components, to the two milestones, and to the goal are in scope. The PAC-Bayes and time-uniform concentration parts of Chapter 16 are reserved for separate work because they need different probability and measure-theory objects.

Selected references

  • Francesco Orabona, Online Learning: A Modern Introduction Using Convex Optimization, arXiv:1912.13213v10, 2026, Chapter 16 and §6.6. arXiv
5 thms3 active usersReviewed
🏆Completed
Convex OptimizationProbability·Captain: mikedeng1

A Modern Introduction to Online Learning X: The Weighted Average Algorithm with a Uniform Prior Has Regret (d/α)(1 + ln(T/d + 1)) on α-Exp-Concave LossesTextbook

Motivation

Online prediction compares the losses of a sequence of decisions with the losses of one fixed decision chosen after the sequence is known. A learner would like a regret guarantee that remains valid for every such competitor, even when the losses arrive in an arbitrary order. Exponentially weighted averaging is a way to make each decision from all feasible points rather than committing to one candidate early. The method gives the Weighted Average Algorithm (WAA) in Chapter 11 of Orabona's text. It is particularly effective when the exponential of each negative loss is concave, a condition called exp-concavity.

The chapter first compares WAA with a competing probability distribution. Theorem 11.5 turns that guarantee into a bound against an individual feasible point. That distinction matters when the benchmark is a concrete portfolio, parameter vector, or other fixed decision, rather than a distribution over decisions. The same chapter later extends the weighted-average viewpoint to mixable losses through the Aggregating Algorithm; this mission focuses on the earlier deterministic-comparator result.

Setting

The feasible set VVV is a closed, bounded, convex, full-dimensional subset of Rd\mathbb R^dRd, with d≥1d\geq1d≥1. Its uniform distribution assigns a measurable subset probability equal to its Euclidean volume divided by Vol⁡(V)\operatorname{Vol}(V)Vol(V). Because VVV is bounded and full-dimensional, its volume is finite and positive. Rounds are t=1,…,Tt=1,\ldots,Tt=1,…,T, with T≥1T\geq1T≥1. In round ttt, the learner predicts xtx_txt​, then observes a loss ℓt\ell_tℓt​ and pays ℓt(xt)\ell_t(x_t)ℓt​(xt​). A competitor u∈Vu\in Vu∈V is fixed across all rounds. Its regret is

RT(u)=∑t=1T(ℓt(xt)−ℓt(u)).R_T(u)=\sum_{t=1}^{T}\bigl(\ell_t(x_t)-\ell_t(u)\bigr).RT​(u)=t=1∑T​(ℓt​(xt​)−ℓt​(u)).

For a parameter α>0\alpha>0α>0, a real-valued loss ℓt\ell_tℓt​ is α\alphaα-exp-concave on VVV if x↦e−αℓt(x)x\mapsto e^{-\alpha\ell_t(x)}x↦e−αℓt​(x) is concave on VVV. Definition 7.31 in the source also permits a loss value of +∞+\infty+∞, with e−α(+∞)=0e^{-\alpha(+\infty)}=0e−α(+∞)=0. The mission defines that extended-valued concept for Proposition 7.34 and uses real-valued losses in its WAA target. The weighted average in round ttt is

xt=∫Vx exp⁡ ⁣(−α∑s=1t−1ℓs(x)) dx∫Vexp⁡ ⁣(−α∑s=1t−1ℓs(x)) dx.x_t=\frac{\int_V x\,\exp\!\left(-\alpha\sum_{s=1}^{t-1}\ell_s(x)\right)\,dx}{\int_V\exp\!\left(-\alpha\sum_{s=1}^{t-1}\ell_s(x)\right)\,dx}.xt​=∫V​exp(−α∑s=1t−1​ℓs​(x))dx∫V​xexp(−α∑s=1t−1​ℓs​(x))dx​.

Thus x1x_1x1​ is the centroid of VVV. This formula is Algorithm 11.1 with the uniform initial distribution and 1/λt=α1/\lambda_t=\alpha1/λt​=α. The existing published EWOO point definition uses precisely this one-based formula, so the mission imports that object. The integrals are with respect to Euclidean Lebesgue measure.

Formalization targets

The main target is Theorem 11.5, equation (11.6): for every u∈Vu\in Vu∈V,

RT(u)≤dα(1+ln⁡(Td+1)).R_T(u)\leq\frac d\alpha\left(1+\ln\left(\frac Td+1\right)\right).RT​(u)≤αd​(1+ln(dT​+1)).

The target fixes the exact dimension, horizon, and exp-concavity dependence from the source. It applies to each fixed competitor, with no distribution chosen in the statement. The milestone list records three source claims in reading order. Proposition 7.34 says that exp-concavity implies convexity, including the book's extended-real convention. Two displays in the proof of Theorem 11.5 describe the translated, contracted set V′={T/dT/d+1u+1T/d+1w:w∈V}V'=\{\frac{T/d}{T/d+1}u+\frac{1}{T/d+1}w:w\in V\}V′={T/d+1T/d​u+T/d+11​w:w∈V}: the exponential losses there are bounded against those at uuu, and its uniform-prior mass is 1/(T/d+1)d1/(T/d+1)^d1/(T/d+1)d. These are the chapter's explicit geometric and loss-comparison targets.

Significance

The theorem converts a distributional WAA guarantee into a bound for each deterministic competitor. Its dimension factor is ddd, and its logarithm contains T/d+1T/d+1T/d+1. This exact expression is sharper than replacing that logarithm by ln⁡(T+1)\ln(T+1)ln(T+1) for d>1d>1d>1; the distinction is visible in the target statement. A related EWOO result with the weaker (d/α)(1+ln⁡(T+1))(d/\alpha)(1+\ln(T+1))(d/α)(1+ln(T+1)) form is already represented in the platform's library, but it does not establish the bound selected here. The source theorem is proved in the book; the remaining work is a machine-checked proof of this statement and the geometric milestones, with the source's constants intact.

The formal development can also support later results on exponentially weighted prediction. Its reusable pieces are the definition of exp-concavity, the extended-real epigraph convention, and the volume identity for a translated scaling of a convex body. The WAA prediction itself is represented by a previously published definition, giving future results a common algorithmic object rather than another copy of the same integral formula.

Difficulty

The general WAA statement compares the algorithm with a probability distribution, whereas Theorem 11.5 compares it with a point uuu. A point mass at uuu is unsuitable as a comparator to a uniform continuous prior: it has no density with respect to that prior, so its relative entropy is infinite. The finite bound depends on identifying a region around uuu that has positive uniform mass while retaining control of every exp-concave loss on that region. A formal treatment must also ensure that the weighted vector integral denotes an actual mean and that the normalization denominator is positive, rather than inheriting a default value from a total integral or division operator.

Formalization scope

The state space is EuclideanSpace ℝ (Fin d), and a loss is indexed from round one; index zero is unused. The prior is uniform on VVV, and the point played is the published one-based EWOO weighted mean. The hypotheses d,T≥1d,T\geq1d,T≥1 make the printed ratios meaningful. Closedness, boundedness, convexity, and positive volume express the source's full-dimensional convex-body setting. The theorem additionally records almost-everywhere measurability of each loss on VVV, and integrability of each exponential weight and its vector moment; these are explicit analytic conditions for the real and Bochner integrals in Lean. They are not algorithmic restrictions or assumptions of the regret inequality's conclusion. Positive volume and integrability rule out a vacuous zero denominator or Lean's zero value for a non-integrable integral.

The source allows extended-real losses, but the goal uses real-valued losses on VVV. Proposition 7.34 separately retains the full extended-real statement through WithTop ℝ: +∞+\infty+∞ maps to zero under the exponential, and convexity is expressed by a convex epigraph. A complete proof needs the exp-concavity and Jensen facts, the weighted-mean properties on convex VVV, the measure and volume properties of V′V'V′, and the distributional WAA inequality from the chapter. Contributions that supply these lemmas under the stated conventions are within scope. The prediction cannot be replaced by an arbitrary point of VVV: the exact uniform-prior weighted mean is part of the goal.

Selected references

  • Francesco Orabona, Online Learning: A Modern Introduction Using Convex Optimization, arXiv:1912.13213v10, 2026. Chapter 11 and Proposition 7.34.
7 thms3 active usersReviewed
🏆Completed
Algorithmic Game TheoryConvex OptimizationOptimization·Captain: mikedeng1

A Modern Introduction to Online Learning XIV: Optimistic FTRL Self-Play on a Smooth Convex-Concave Function Has Duality Gap at Most a Constant over TTextbook

Motivation

A saddle-point problem asks one player to choose a decision that keeps a payoff small while another chooses a decision that makes it large. Such problems include two-player zero-sum games and constrained optimization written through a Lagrangian. Online convex optimization supplies a way to construct approximate saddle points from repeated play: each player responds to losses generated by the other's current decision. Chapter 15 of Orabona's textbook develops this connection and studies the resulting averages. The chapter's earlier online-to-saddle-point result gives a duality-gap certificate from the two players' regrets; §15.5 asks how the certificate improves when the payoff is smooth and both players use predictions of their next gradients.

The minimax equality for continuous convex-concave payoffs on compact convex sets is already a known mathematical result. Theorem 15.11 in the book is a case of Sion's minimax theorem, and a proved Sion theorem already exists in the platform library. This mission centers on the algorithmic rate in Theorem 15.26, which has a separate statement and a different role.

Setting

Let XXX and YYY be nonempty closed convex sets in finite-dimensional real normed spaces, with norms ∥⋅∥X\|\cdot\|_X∥⋅∥X​ and ∥⋅∥Y\|\cdot\|_Y∥⋅∥Y​. The payoff f(x,y)f(x,y)f(x,y) is real-valued, convex as a function of xxx for fixed yyy, and concave as a function of yyy for fixed xxx. The XXX-player tries to minimize it; the YYY-player tries to maximize it. A saddle point would be a pair whose XXX coordinate minimizes the payoff against its YYY coordinate and whose YYY coordinate maximizes it against its XXX coordinate.

The duality gap of a pair (x,y)(x,y)(x,y), when the extrema are attained, is max⁡v∈Yf(x,v)−min⁡u∈Xf(u,y)\max_{v\in Y}f(x,v)-\min_{u\in X}f(u,y)maxv∈Y​f(x,v)−minu∈X​f(u,y). It is nonnegative and measures how far the pair is from the simultaneous best-response condition. For iterates indexed from 111 through TTT, the book evaluates this gap at the averages xˉT=T−1∑t=1Txt\bar x_T=T^{-1}\sum_{t=1}^{T}x_txˉT​=T−1∑t=1T​xt​ and yˉT=T−1∑t=1Tyt\bar y_T=T^{-1}\sum_{t=1}^{T}y_tyˉ​T​=T−1∑t=1T​yt​.

Optimistic follow-the-regularized-leader, or optimistic FTRL, chooses each point by minimizing a regularizer plus all earlier linearized losses and a hint for the next loss. Algorithm 15.6 sets the hint equal to the previous observed partial gradient, starting with zero at round 000. The XXX-player's observed linear functional is gX,t=∇xf(xt,yt)g_{X,t}=\nabla_x f(x_t,y_t)gX,t​=∇x​f(xt​,yt​); the YYY-player's is gY,t=−∇yf(xt,yt)g_{Y,t}=-\nabla_y f(x_t,y_t)gY,t​=−∇y​f(xt​,yt​). Their regularizers ψX\psi_XψX​ and ψY\psi_YψY​ have strong-convexity parameters λX\lambda_XλX​ and λY\lambda_YλY​ with respect to their respective norms. The norm of each partial gradient is the corresponding dual norm.

Smoothness here means that the partial gradients exist on an open neighborhood of X×YX\times YX×Y and satisfy four Lipschitz inequalities. Changes in the XXX coordinate affect the XXX gradient by at most LXXL_{XX}LXX​ times the XXX distance; changes in the YYY coordinate affect it by at most LXYL_{XY}LXY​ times the YYY distance. For the YYY gradient the cross constant is again LXYL_{XY}LXY​, and the own-coordinate constant is LYYL_{YY}LYY​, exactly as in equations (15.7)–(15.10) of the textbook.

Formalization targets

The goal is the exact finite-horizon inequality of Theorem 15.26. Fix α>0\alpha>0α>0 and parameters satisfying λX≥22(LXX+LXYα)\lambda_X\ge2\sqrt2(L_{XX}+L_{XY}\alpha)λX​≥22​(LXX​+LXY​α) and λY≥22(LYY+LXY/α)\lambda_Y\ge2\sqrt2(L_{YY}+L_{XY}/\alpha)λY​≥22​(LYY​+LXY​/α). For any minimizer xT′x'_TxT′​ of f(⋅,yˉT)f(\cdot,\bar y_T)f(⋅,yˉ​T​) and maximizer yT′y'_TyT′​ of f(xˉT,⋅)f(\bar x_T,\cdot)f(xˉT​,⋅), the target is

max⁡v∈Yf(xˉT,v)−min⁡u∈Xf(u,yˉT)≤ψX(xT′)−ψX(x1)+ψY(yT′)−ψY(y1)+∥gX,1∥X,∗2/λX+∥gY,1∥Y,∗2/λYT.\max_{v\in Y}f(\bar x_T,v)-\min_{u\in X}f(u,\bar y_T) \le \frac{\psi_X(x'_T)-\psi_X(x_1)+\psi_Y(y'_T)-\psi_Y(y_1)+\|g_{X,1}\|_{X,*}^{2}/\lambda_X+\|g_{Y,1}\|_{Y,*}^{2}/\lambda_Y}{T}.v∈Ymax​f(xˉT​,v)−u∈Xmin​f(u,yˉ​T​)≤TψX​(xT′​)−ψX​(x1​)+ψY​(yT′​)−ψY​(y1​)+∥gX,1​∥X,∗2​/λX​+∥gY,1​∥Y,∗2​/λY​​.

The milestones are the linear-loss, fixed-regularizer instance of Theorem 7.47 used in the chapter, Theorem 15.12 relating player regrets to the averaged duality gap, and the explicit summed-regret inequality on page 262. Together they specify the algorithmic input, the conversion to a saddle-point certificate, and the exact stability coefficients in the goal.

Significance

The bound certifies the quality of the averaged pair after a specified number of rounds. When the regularizers are bounded on the feasible sets, its numerator is uniformly bounded and the gap is O(1/T)O(1/T)O(1/T), as stated in the paragraph before Theorem 15.26. Without that further boundedness, the displayed inequality remains the theorem's finite-horizon statement, but its numerator can depend on TTT through xT′x'_TxT′​ and yT′y'_TyT′​. The mission name's “constant over TTT” describes the bounded-regularizer setting; the formal goal preserves the more general printed bound.

The book proves the theorem mathematically. The work here is a machine-checked representation of Algorithm 15.6, its smoothness assumptions, the exact parameter thresholds, and the finite-horizon inequality. The run and smoothness definitions can be reused for other optimistic methods in the chapter, including Theorem 15.28. Existing library results about minimax equality address existence and equality of optimal values; they do not supply the quantitative guarantee for this run.

Difficulty

A standard regret bound for each player leaves variation terms involving successive iterates. With two players, each player's gradient variation depends on both players' movement and on different norms. A bound that discards this coupling loses the displayed coefficients and the rate in Theorem 15.26. The formal statement therefore keeps both sequences, both dual norms, and all four smoothness inequalities visible. It also has to treat the first round separately: the previous-gradient hint is zero there, while later rounds compare two observed gradients.

Formalization scope

The Lean development uses finite-dimensional real normed spaces; a partial gradient is a continuous linear functional and its operator norm is the book's dual norm. The payoff and regularizers are represented by real-valued functions on the ambient spaces, restricted to XXX and YYY when optimized. Smoothness includes an open set containing the feasible product. Strong convexity uses Mathlib's set-relative StrongConvexOn. The run predicate enforces both argmin updates, and index 000 is used only for the zero hint. It constrains exactly the rounds 1,…,T1,\ldots,T1,…,T. In the optimistic FTRL milestone, the look-ahead point xT+1x_{T+1}xT+1​ is the minimizer of the regularized loss sum with a zero hint. The book's proof uses the hint of round TTT there instead, which makes that milestone's first inequality false.

The goal assumes T≥1T\ge1T≥1, positive α,λX,λY\alpha,\lambda_X,\lambda_Yα,λX​,λY​, and nonnegative smoothness constants. These record the book's parameter conventions and keep every denominator meaningful. The extrema in the conclusion are represented by arbitrary attained minimizer and maximizer witnesses, matching the theorem's nonempty-argmin and nonempty-argmax hypotheses. No boundedness or compactness of XXX or YYY is imposed. An arbitrary sequence of feasible points cannot satisfy the goal: the optimistic FTRL argmin updates are part of the run hypothesis. Contributions that establish the two-player regret estimate, the regret-to-gap step, or the final inequality against this run model directly support the mission.

Selected references

  • Francesco Orabona, Online Learning: A Modern Introduction Using Convex Optimization, arXiv:1912.13213v10, 2026, Chapter 15, pp. 244–269; Theorem 7.47, pp. 129–130. arXiv
  • Maurice Sion, On general minimax theorems, Pacific Journal of Mathematics 8 (1958), pp. 171–176. DOI
5 thms3 active usersReviewed
🏆Completed
Convex OptimizationOptimizationProbability·Captain: mikedeng1

Linear Convergence of Gradient and Proximal-Gradient Methods Under the Polyak-Łojasiewicz Condition VI: Randomized Proximal Coordinate Descent Has Expected Rate (1 − μ/(dL))^kResearch Paper

Motivation

Many problems in machine learning and statistics have the composite form min⁡xF(x)=f(x)+g(x)\min_x F(x)=f(x)+g(x)minx​F(x)=f(x)+g(x), where fff is a smooth loss and ggg is a non-smooth regularizer such as λ∥x∥1\lambda\|x\|_1λ∥x∥1​. When ggg is separable, g(x)=∑igi(xi)g(x)=\sum_i g_i(x_i)g(x)=∑i​gi​(xi​), the method of choice at scale is often randomized proximal coordinate descent: at each step one coordinate is drawn at random and a one-dimensional proximal problem is solved in that coordinate only. Each step is cheap, and the method underlies standard solvers for ℓ1\ell_1ℓ1​-regularized regression and for the dual of support vector machines.

Classical linear-convergence guarantees for this method assume strong convexity of fff or of FFF (Richtárik and Takáč, 2014), or rely on the Kurdyka–Łojasiewicz property with constants that are hard to make explicit. Karimi, Nutini and Schmidt (arXiv:1608.04636, ECML PKDD 2016) introduced the proximal-PL inequality, a generalization of the Polyak–Łojasiewicz condition to composite problems, and showed that it alone yields linear rates for proximal-gradient methods. Theorem 6 of that paper, the subject of this mission, extends the analysis to randomized proximal coordinate descent.

Timeline. Polyak (1963) and Łojasiewicz (1963) introduced the gradient-domination inequality for smooth problems. Nesterov (2012) analysed randomized coordinate descent for smooth convex problems; Richtárik and Takáč (2014) extended it to composite separable problems under strong convexity. Karimi, Nutini and Schmidt (2016) obtained the rate (1−μ/(dL))k(1-\mu/(dL))^k(1−μ/(dL))k for the proximal coordinate method under proximal-PL alone (Theorem 6, p. 10, proof in Appendix H, p. 23).

Setting

Write Rd\mathbb R^dRd for the Euclidean space with inner product ⟨⋅,⋅⟩\langle\cdot,\cdot\rangle⟨⋅,⋅⟩ and norm ∥⋅∥\|\cdot\|∥⋅∥; xix_ixi​ is the iii-th coordinate of xxx, eie_iei​ the iii-th basis vector, and ∇if(x)\nabla_i f(x)∇i​f(x) the iii-th coordinate of the gradient. The data are:

  • a differentiable f:Rd→Rf:\mathbb R^d\to\mathbb Rf:Rd→R whose gradient is LLL-Lipschitz, ∥∇f(x)−∇f(y)∥≤L∥x−y∥\|\nabla f(x)-\nabla f(y)\|\le L\|x-y\|∥∇f(x)−∇f(y)∥≤L∥x−y∥, with L>0L>0L>0;
  • convex univariate functions g1,…,gd:R→Rg_1,\dots,g_d:\mathbb R\to\mathbb Rg1​,…,gd​:R→R and the separable g(x)=∑igi(xi)g(x)=\sum_i g_i(x_i)g(x)=∑i​gi​(xi​);
  • F=f+gF=f+gF=f+g, which has a minimizer x∗x^*x∗; write F∗=F(x∗)F^*=F(x^*)F∗=F(x∗).

For α>0\alpha>0α>0 the quantity

Dg(x,α)=−2αmin⁡y[⟨∇f(x),y−x⟩+α2∥y−x∥2+g(y)−g(x)]\mathcal D_g(x,\alpha)=-2\alpha\min_y\Big[\langle\nabla f(x),y-x\rangle+\frac\alpha2\|y-x\|^2+g(y)-g(x)\Big]Dg​(x,α)=−2αymin​[⟨∇f(x),y−x⟩+2α​∥y−x∥2+g(y)−g(x)]

measures the best decrease of a quadratic upper model of FFF at xxx. FFF satisfies the proximal-PL inequality with constant μ>0\mu>0μ>0 if

12 Dg(x,L) ≥ μ (F(x)−F∗)for all x.\tfrac12\,\mathcal D_g(x,L)\ \ge\ \mu\,\big(F(x)-F^*\big)\qquad\text{for all }x.21​Dg​(x,L) ≥ μ(F(x)−F∗)for all x.

The coordinate-wise proximal-gradient method starts at x0x_0x0​. At step kkk it draws iki_kik​ uniformly from {1,…,d}\{1,\dots,d\}{1,…,d}, independently of the past, and sets xk+1=xk+α⋆eikx_{k+1}=x_k+\alpha^\star e_{i_k}xk+1​=xk​+α⋆eik​​, where α⋆\alpha^\starα⋆ minimizes

α ∇ikf(xk)+L2α2+gik(xk,ik+α)−gik(xk,ik).(16)\alpha\,\nabla_{i_k}f(x_k)+\frac L2\alpha^2+g_{i_k}(x_{k,i_k}+\alpha)-g_{i_k}(x_{k,i_k}).\tag{16}α∇ik​​f(xk​)+2L​α2+gik​​(xk,ik​​+α)−gik​​(xk,ik​​).(16)

Formalization targets

Goal: Theorem 6 (p. 10)

For every k≥0k\ge0k≥0,

E[F(xk)−F∗] ≤ (1−μdL)k[F(x0)−F∗].(17)\mathbb E\big[F(x_k)-F^*\big]\ \le\ \Big(1-\frac{\mu}{dL}\Big)^k\big[F(x_0)-F^*\big].\tag{17}E[F(xk​)−F∗] ≤ (1−dLμ​)k[F(x0​)−F∗].(17)

The constant μ/(dL)\mu/(dL)μ/(dL) is the paper's; no further constant is fixed.

Milestones (Appendix H, p. 23)

  1. (62), the coordinate upper model: F(x+tei)−F(x)≤t∇if(x)+L2t2+gi(xi+t)−gi(xi)F(x+te_i)-F(x)\le t\nabla_i f(x)+\frac L2t^2+g_i(x_i+t)-g_i(x_i)F(x+tei​)−F(x)≤t∇i​f(x)+2L​t2+gi​(xi​+t)−gi​(xi​).
  2. The expected-improvement identity: 1d∑imin⁡t[t∇if(x)+L2t2+gi(xi+t)−gi(xi)]=−12LdDg(x,L)\frac1d\sum_i\min_{t}\big[t\nabla_i f(x)+\frac L2t^2+g_i(x_i+t)-g_i(x_i)\big]=-\frac1{2Ld}\mathcal D_g(x,L)d1​∑i​mint​[t∇i​f(x)+2L​t2+gi​(xi​+t)−gi​(xi​)]=−2Ld1​Dg​(x,L).
  3. (63), the expected descent: E[F(xk+1)∣xk=x]≤F(x)−12LdDg(x,L)\mathbb E[F(x_{k+1})\mid x_k=x]\le F(x)-\frac1{2Ld}\mathcal D_g(x,L)E[F(xk+1​)∣xk​=x]≤F(x)−2Ld1​Dg​(x,L).

A companion statement (not a milestone, not in the paper) asserts that the one-dimensional step (16) exists for convex gig_igi​ and L>0L>0L>0.

Significance

The result. Theorem 6 gives a linear rate for proximal coordinate descent on every composite problem satisfying proximal-PL, a class that contains ℓ1\ell_1ℓ1​-regularized least squares and other problems that are not strongly convex. The rate matches that of plain randomized coordinate descent under PL (Theorem 3 of the paper) with the composite term added, and the factor ddd is the usual price of updating one coordinate per step.

Formalizing it. The result is proved on paper; no machine-checked proof is known. A formal proof requires the coordinate descent lemma, the exchange of a separable minimum with a sum, and an exact accounting of the expectation over random index sequences. The definitions of Dg\mathcal D_gDg​, proximal-PL and the random-coordinate run built here are reusable for other proximal and coordinate methods.

Difficulty

The one-step argument is short once the expected-improvement identity is available; the difficulty is in that identity. Its left side averages over coordinates the minima of one-dimensional models, while Dg\mathcal D_gDg​ is a minimum over the whole space. Equating them requires that the full-space minimum split across coordinates, which is where separability of ggg is used, and that each minimum be attained, which requires convexity of the gig_igi​ and L>0L>0L>0. A second point is the expectation: the iterate after kkk steps is a function of the whole index word (i0,…,ik−1)(i_0,\dots,i_{k-1})(i0​,…,ik−1​), and the one-step expected descent must be propagated through the tower of conditional expectations to the kkk-step bound.

Formalization scope

  • The space is EuclideanSpace ℝ (Fin d) with d≥1d\ge1d≥1 (0 < d); the gradient is Mathlib's gradient f, and ∇if(x)\nabla_i f(x)∇i​f(x) is gradient f x i.
  • The Lipschitz gradient is the literal hypothesis ∥∇f(x)−∇f(y)∥≤L∥x−y∥\|\nabla f(x)-\nabla f(y)\|\le L\|x-y\|∥∇f(x)−∇f(y)∥≤L∥x−y∥ with L>0L>0L>0 and Differentiable ℝ f, as in the setup of Theorem 5. Appendix H uses only coordinate-wise Lipschitz continuity, which follows.
  • X∗≠∅X^*\neq\emptysetX∗=∅ is a point x∗x^*x∗ minimizing FFF; F∗F^*F∗ is F(x∗)F(x^*)F(x∗), never a free real.
  • The gig_igi​ are real-valued convex functions. Extended-valued gig_igi​ (box indicators, used for the SVM dual in §4.4) are outside this formalization.
  • The minimum in Dg\mathcal D_gDg​ is a real infimum; it is attained, hence equal to the paper's minimum, for convex gig_igi​ and L>0L>0L>0, and every statement using it carries these hypotheses.
  • The update is any step selection sel(x,i)=x+t⋆ei\mathrm{sel}(x,i)=x+t^\star e_isel(x,i)=x+t⋆ei​ with t⋆t^\start⋆ a minimizer of (16); the selection is a hypothesis-constrained parameter, never an arbitrary map. Equation (16) as printed assigns the minimizing step to xk+1x_{k+1}xk+1​; the iterate xk+α⋆eikx_k+\alpha^\star e_{i_k}xk​+α⋆eik​​ of Appendix H is what is encoded.
  • E[F(xk)−F∗]\mathbb E[F(x_k)-F^*]E[F(xk​)−F∗] is the exact average over all dkd^kdk index words, each of probability d−kd^{-k}d−k, which is the law of i.i.d. uniform coordinates.
  • A formalization that leaves the update unconstrained, takes F∗F^*F∗ as a free constant, or computes Dg\mathcal D_gDg​ as a junk infimum would make the goal false or vacuous; these are ruled out by the conventions above.

Contributions welcome: proofs of the three milestones and the goal, the existence of the coordinate step, and a version with extended-valued gig_igi​.

Selected references

  • H. Karimi, J. Nutini, M. Schmidt, Linear Convergence of Gradient and Proximal-Gradient Methods Under the Polyak-Łojasiewicz Condition, ECML PKDD 2016; arXiv:1608.04636v4, 2020. https://arxiv.org/abs/1608.04636
  • P. Richtárik, M. Takáč, Iteration complexity of randomized block-coordinate descent methods for minimizing a composite function, Mathematical Programming 144, 2014. https://doi.org/10.1007/s10107-012-0614-z
  • Yu. Nesterov, Efficiency of coordinate descent methods on huge-scale optimization problems, SIAM J. Optim. 22(2), 2012. https://doi.org/10.1137/100802001
  • B. T. Polyak, Gradient methods for minimizing functionals, USSR Comput. Math. Math. Phys. 3(4), 1963. https://doi.org/10.1016/0041-5553(63)90382-3
7 thms3 active usersReviewed
Bandit AlgorithmsProbability·Captain: mikedeng1

A Modern Introduction to Online Learning VIII: Tsallis-INF Has Adversarial Regret 32L∞√((d − 1)T) and Stochastic Pseudo-Regret 256L∞² Σ_i (1 + ln T)/Δ_iTextbook

Why bandit feedback matters

In online learning with bandit feedback, the learner chooses one action and sees only that action's loss. This is the information pattern of sequential trials in which testing one option does not reveal the outcomes of the others. The learner must keep its cumulative loss close to that of a good fixed action while learning which actions are good. A full-information method can update from every arm's loss each round; here the unseen losses have to be handled through the observations the learner actually receives. Chapter 9 of Orabona's text develops this distinction for multi-armed bandits.

Two environments are treated in the chapter. With an oblivious adversary, all losses are fixed before play, although the learner does not know them. With a stochastic environment, each arm has a fixed loss distribution and new losses are drawn across rounds. An algorithm with a guarantee in both environments does not need to be told which generated its data. Theorem 9.17 gives such a pair of guarantees for the FTRL version of Tsallis-INF, using the precise constants stated in the pinned version of the book.

Arms, losses, and the algorithm

There are ddd arms and TTT rounds, numbered 1,…,T1,\ldots,T1,…,T. In round ttt, the learner forms a probability vector xt=(xt,1,…,xt,d)x_t=(x_{t,1},\ldots,x_{t,d})xt​=(xt,1​,…,xt,d​) on the simplex Δd−1\Delta^{d-1}Δd−1, draws an arm AtA_tAt​ with probability xt,ix_{t,i}xt,i​ of choosing iii, and observes only gt,Atg_{t,A_t}gt,At​​. The loss gt,ig_{t,i}gt,i​ of an unchosen arm is not revealed. For an oblivious adversary, the entire table (gt,i)(g_{t,i})(gt,i​) is deterministic and the expectation is over the learner's draws.

Algorithm 9.6 uses the importance-weighted estimate g~t,i=gt,i/xt,i\widetilde g_{t,i}=g_{t,i}/x_{t,i}g​t,i​=gt,i​/xt,i​ if At=iA_t=iAt​=i, and 000 otherwise. Its Tsallis-INF FTRL prediction minimizes, over x∈Δd−1x\in\Delta^{d-1}x∈Δd−1,

⟨∑s=1t−1g~s,x⟩−4L∞t∑i=1dxi.\left\langle\sum_{s=1}^{t-1}\widetilde g_s,x\right\rangle -4L_\infty\sqrt t\sum_{i=1}^{d}\sqrt{x_i}.⟨s=1∑t−1​g​s​,x⟩−4L∞​t​i=1∑d​xi​​.

The minimizer has positive coordinates. Thus each denominator xt,ix_{t,i}xt,i​ used by the estimator is positive. In the stochastic environment, arm iii has mean μi\mu_iμi​; if jjj is a best arm, its gap is Δi=μi−μj\Delta_i=\mu_i-\mu_jΔi​=μi​−μj​. The pull count ST,iS_{T,i}ST,i​ is the number of rounds through TTT in which arm iii is selected. Stochastic pseudo-regret compares the learner's expected loss with TμjT\mu_jTμj​, the expected loss of always playing the best arm.

Formalization targets

The mission's main theorem is the adversarial claim of Theorem 9.17, p. 165. For every fixed comparator arm kkk and losses in [0,L∞][0,L_\infty][0,L∞​],

E ⁣[∑t=1T(gt,At−gt,k)]≤32L∞(d−1)T.\mathbb E\!\left[\sum_{t=1}^{T}(g_{t,A_t}-g_{t,k})\right] \le 32L_\infty\sqrt{(d-1)T}.E[t=1∑T​(gt,At​​−gt,k​)]≤32L∞​(d−1)T​.

Taking kkk to be the best arm in the fixed loss table yields the stated adversarial regret bound. The companion item formalizes the second claim of the same theorem: if the losses are stochastic with a unique best arm and T≥1T\ge1T≥1, then

P ⁣-⁡RegretT≤256L∞2∑i:Δi≠01+ln⁡TΔi.\mathrm{P\!\operatorname{-}Regret}_T \le256L_\infty^2\sum_{i:\Delta_i\ne0}\frac{1+\ln T}{\Delta_i}.P-RegretT​≤256L∞2​i:Δi​=0∑​Δi​1+lnT​.

The milestone list contains the chapter's pseudo-regret decomposition (Lemma 9.10), its sub-Gaussian sample-mean bound (Corollary 9.9, reused from an existing formal statement), and the probability-sensitive Tsallis-INF inequality (9.8). These are distinct mathematical claims rather than alternate names for the final bound.

What the result provides

The adversarial bound grows like T\sqrt{T}T​ and remains valid for any fixed bounded loss table. The stochastic bound instead measures the difficulty of distinguishing arms through the positive gaps Δi\Delta_iΔi​ and grows logarithmically with the horizon for fixed gaps. Both refer to the same FTRL update and the same importance-weighted estimator. The finite-horizon constants 323232 and 256256256 are part of the target, not unspecified factors hidden in asymptotic notation.

Both claims are proved in the book. The work here is to give their hypotheses, random experiment, and algorithm a machine-checkable Lean interface, then formalize the known arguments. At this drafting stage the theorem bodies are open Lean problems; a successful build checks that their statements are well formed, not that the regret bounds have been proved. The resulting finite-history bandit protocol and stochastic pseudo-regret definitions can also support other loss-based bandit algorithms. The existing formalized sub-Gaussian tail bound can be reused for concentration results elsewhere in the chapter.

Where the proof is difficult

The learner observes only one coordinate, while regret is measured against losses across all arms. The importance-weighted estimate corrects the missing observations in expectation, but its size can grow when a selected arm has small probability. Consequently, a bound based only on a uniform estimate of its magnitude loses the probability dependence needed for the theorem's constants. The stochastic claim adds another difficulty: the learner's choices depend on past observed losses, so a current loss must be related to a random pull count without letting the choice depend on that current loss.

There are two tempting but invalid shortcuts. Permitting xtx_txt​ to inspect gtg_tgt​ would change the bandit protocol, and taking an expectation under an arbitrary path measure would no longer describe the arm draws of Algorithm 9.6. The formal statements tie each prediction to its FTRL objective and each path probability to the predictions made along that history.

Formalization scope

Arms are Fin d; the main theorem permits d=1d=1d=1, where both its regret and bound are zero. Rounds retain the book's 1,…,T1,\ldots,T1,…,T indexing. A finite path has a dummy coordinate zero, fixed to an arbitrary anchor arm; this coordinate is never played or scored. The adversarial losses are a fixed table, selected before any draw. A path has probability equal to the product of the algorithm's conditional arm probabilities, and regret against a fixed arm is a finite weighted sum over paths. The FTRL run predicate requires simplex membership, strict positivity, minimization against every simplex point, and dependence only on past arm choices. The last two properties specify the intended update on every possible history.

The stochastic model first samples a loss table from a probability measure and then draws an arm path conditionally using the same Algorithm 9.6. Current-round loss vectors are independent across rounds, while coordinates within one round need not be independent. Each arm has the same mean in every round. The policy is a measurable function of the observed history (past arms and their losses); its expected loss is explicitly integrable. Loss bounds are imposed on every outcome in the chosen sample space, which can be taken as the support of bounded distributions. The stochastic statement adds T≥1T\ge1T≥1 to make the logarithmic formula's domain explicit, and the unique-best-arm assumption ensures every denominator Δi\Delta_iΔi​ in the displayed sum is positive. Lemma 9.10 also spells out the integrability of its means and pull counts.

The development needs finite probability laws, real integrals, simplex geometry, and properties of the square-root regularizer. The independent tail bound for sample means is already available as a referenced theorem. Contributions that establish existence and interiority of the FTRL minimizer, the normalization of the path law, or the stated regret inequalities fit this scope. A choice rule that sees the current loss, a path law unrelated to xtx_txt​, or a zero integrability default does not satisfy the intended bandit model.

Selected references

  • Francesco Orabona, Online Learning: A Modern Introduction Using Convex Optimization, arXiv:1912.13213v10, 2026, Chapter 9, pp. 151–169. arXiv
6 thms3 active usersReviewed
🏆Completed
Convex OptimizationProbability·Captain: mikedeng1

A Modern Introduction to Online Learning IV: Every Deterministic Online Linear Optimization Algorithm Has Regret at Least √2·LD√T/4Textbook

Why lower bounds matter

Online learning studies decisions made before their costs are known. An algorithm selects a point, observes a loss, and repeats. Regret compares its accumulated loss with the loss of one fixed point chosen after the sequence is known. Upper bounds show what a particular algorithm can guarantee; a lower bound shows what every algorithm must sometimes suffer. The distinction matters when deciding whether a rate such as T\sqrt TT​ can be improved by changing the update rule. Chapter 5 of Orabona's Online Learning: A Modern Introduction Using Convex Optimization addresses this question for online linear optimization.

The chapter begins with bounded decision sets and then examines unbounded decision spaces. On bounded sets, the main result gives a quantitative obstruction with an explicit constant for every deterministic algorithm. The book compares that order with the regret upper bound for online subgradient descent: the square-root dependence on the number of rounds is optimal up to a constant factor under the same loss-scale and diameter parameters (Orabona, Theorem 5.1 and discussion, pp. 50–51). The same chapter shows a different phenomenon for unprojected subgradient descent with steps that decay too slowly: its regret can have order T2−αT^{2-\alpha}T2−α for 0<α<10<\alpha<10<α<1 (Orabona, Theorem 5.4, p. 52).

The online linear optimization setting

Let VVV be a nonempty subset of Euclidean space Rd\mathbb R^dRd. On round ttt, a deterministic online algorithm selects xt∈Vx_t\in Vxt​∈V using only the loss vectors from rounds 1,…,t−11,\ldots,t-11,…,t−1. It then observes gt∈Rdg_t\in\mathbb R^dgt​∈Rd and pays the linear loss ⟨gt,xt⟩\langle g_t,x_t\rangle⟨gt​,xt​⟩. The adversary may select any loss vectors satisfying ∥gt∥2≤L\|g_t\|_2\le L∥gt​∥2​≤L. This is the online linear optimization protocol used in the chapter; it also supplies lower bounds for online convex optimization because linear functions are convex (Orabona, §5.1, p. 50).

For a fixed competitor u∈Vu\in Vu∈V, define

Regret⁡T(u)=∑t=1T⟨gt,xt⟩−∑t=1T⟨gt,u⟩.\operatorname{Regret}_T(u) =\sum_{t=1}^{T}\langle g_t,x_t\rangle -\sum_{t=1}^{T}\langle g_t,u\rangle.RegretT​(u)=t=1∑T​⟨gt​,xt​⟩−t=1∑T​⟨gt​,u⟩.

The competitor is fixed across all TTT rounds, although it may be selected after the loss sequence. The algorithm must make its own choices causally. These two quantifier positions are part of the mathematical claim. In particular, the hard sequence and competitor can depend on the algorithm, while the algorithm cannot use a future loss vector to choose its current prediction.

The bounded setting also uses the diameter D=max⁡v,w∈V∥v−w∥2D=\max_{v,w\in V}\|v-w\|_2D=maxv,w∈V​∥v−w∥2​. The theorem assumes VVV is closed, bounded, and convex, and that D>0D>0D>0. In finite-dimensional Euclidean space, closedness and boundedness make VVV compact, so this maximum exists. The positive diameter excludes a singleton feasible set, where the algorithm and the only comparator coincide. These hypotheses and the diameter convention are stated in Theorem 5.1, pp. 50–51.

Formalization targets

Worst-case regret on every bounded convex set

The goal is the chapter's Theorem 5.1. For every deterministic algorithm on such a set VVV, every positive integer TTT, and every nonnegative loss scale LLL, there are bounded loss vectors and a competitor u∈Vu\in Vu∈V for which

Regret⁡T(u)≥24LDT.\operatorname{Regret}_T(u)\ge \frac{\sqrt2}{4}LD\sqrt T.RegretT​(u)≥42​​LDT​.

The book prints the exact factor 2/4\sqrt2/42​/4. The scale LLL appears in its conclusion and in the restriction ∥gt∥2≤L\|g_t\|_2\le L∥gt​∥2​≤L but is not introduced before that sentence, so the formal statement quantifies it explicitly. At L=0L=0L=0, the bound reads zero; for L>0L>0L>0 it gives a positive lower bound when D>0D>0D>0 and T≥1T\ge1T≥1 (Orabona, Theorem 5.1, pp. 50–51).

Supporting targets

Three milestones support and contextualize the goal. The first is the real-coefficient, p=1p=1p=1 lower-bound case of Khintchine's inequality, which controls the mean absolute value of a weighted sum of independent random signs (Orabona, Theorem A.14, p. 296). The second is the finite first-moment principle: a random sequence whose expected loss is at least KKK has an outcome with loss at least KKK (Orabona, §5.1, p. 50). The third records the lower-bound example for unprojected online subgradient descent with step size t−αt^{-\alpha}t−α:

Regret⁡T(0)≥12ϕ(α)T2−α,ϕ(α)=12−α+(1/2)1−α−11−α.\operatorname{Regret}_T(0)\ge\frac12\phi(\alpha)T^{2-\alpha}, \qquad \phi(\alpha)=\frac1{2-\alpha}+\frac{(1/2)^{1-\alpha}-1}{1-\alpha}.RegretT​(0)≥21​ϕ(α)T2−α,ϕ(α)=2−α1​+1−α(1/2)1−α−1​.

This last target uses the book's threshold T≥2/((1−α)ϕ(α))T\ge2/((1-\alpha)\phi(\alpha))T≥2/((1−α)ϕ(α)) and 0<α<10<\alpha<10<α<1 (Orabona, Theorem 5.4, p. 52).

What the results establish

The goal rules out an o(T)o(\sqrt T)o(T​) worst-case regret guarantee for deterministic online linear optimization under a fixed positive loss scale and positive diameter. Because the result applies to every nonempty bounded closed convex VVV with D>0D>0D>0, it does not rely on a symmetric feasible set or a chosen cube. The explicit constant permits direct comparison with concrete algorithmic guarantees. Theorem 5.4 separately identifies a failure mode for unprojected online subgradient descent on an unbounded domain. Both conclusions are proved in the book; the remaining task here is to give their statements and supporting objects machine-checkable Lean meanings, then prove them in Lean. Hazan's related lower bound uses a hypercube instance and an unspecified absolute constant, so it is related context rather than the same formal target.

The mathematical difficulty

A hard sequence must be chosen after the algorithm is fixed and must respect the algorithm's information at every round. A sequence selected independently of that quantifier order need not be hard: an algorithm tailored to a known sequence can choose its best fixed point on every round. The bound must also hold on any convex body with positive diameter, including sets whose geometry is asymmetric. The proof therefore needs a way to extract a uniform numerical lower bound from only the diameter and the loss norm bound, while retaining the algorithm's causal dependence on earlier losses (Orabona, §5.1, pp. 50–51).

Formalization scope

Lean represents Rd\mathbb R^dRd by EuclideanSpace ℝ (Fin d), with vectors gtg_tgt​ and decisions xtx_txt​ indexed from round one; index zero is unused. A deterministic algorithm is a function of a full loss sequence together with an explicit non-anticipation condition. Its outputs are required to lie in VVV. Regret is a finite sum against one fixed competitor. The positive horizon, positive diameter, closedness, boundedness, convexity, and nonemptiness are explicit. This rules out trivializing the goal by using a singleton domain, an anticipatory algorithm, or a fixed hard sequence that precedes the algorithm's quantifier.

The Khintchine milestone uses the uniform distribution on a finite Boolean cube; each Boolean coordinate represents an independent Rademacher sign. Its expectation is an ordinary finite average and is always defined. The finite first-moment milestone uses arbitrary nonnegative weights summing to one. Theorem 5.4 uses linear losses of norm at most one in arbitrary positive Euclidean dimension, which are convex and 111-Lipschitz. Its auxiliary function ϕ\phiϕ is used only for 0<α<10<\alpha<10<α<1, where the printed formula's denominators are nonzero. The separate limiting claim about ϕ\phiϕ and the later coin-betting lower bounds are outside this mission's current milestones. Contributions establishing the Rademacher inequality, the geometric lower bound, the causal regret argument, and the unprojected descent example all fit this scope.

Selected references

  • Francesco Orabona, Online Learning: A Modern Introduction Using Convex Optimization, arXiv:1912.13213v10, 2026. arXiv.
  • Elad Hazan, Introduction to Online Convex Optimization, Foundations and Trends in Optimization, 2016; arXiv:1909.05207. arXiv.
6 thms3 active usersReviewed
🏆Completed
Convex OptimizationOptimization·Captain: mikedeng1

A Modern Introduction to Online Learning VI: Follow-the-Regularized-Leader with a µ-Strongly Convex Regularizer Has Regret ((ψ(u) − min ψ)/α + α)·L√T/√µTextbook

Motivation

In online convex optimization, a learner chooses a point before seeing the loss assigned to that round. Its performance is measured against a fixed point that could have been chosen after the full sequence became known. A learner that merely minimizes accumulated past losses can react sharply when successive losses favor different points. Follow-the-Regularized-Leader (FTRL) changes that decision by adding a regularizer to the past losses. Chapter 7 of Orabona's Online Learning: A Modern Introduction Using Convex Optimization develops the resulting regret guarantees for a general norm, time-varying regularizers, and convex losses. This mission targets its first explicit square-root regret bound, Corollary 7.6, together with the equality and stability statements used to analyze the algorithm.

The general-norm formulation matters because different feasible sets make different norms natural. A Euclidean constraint set and a simplex, for example, need not give useful bounds under the same geometry. The chapter states the central inequality using the norm of the learner's decisions and the corresponding dual norm of loss subgradients. Its regularizer may vary with time, which lets the same result describe a learning-rate schedule chosen before each loss is observed.

Setting

Let VVV be a nonempty, closed, convex subset of a finite-dimensional real normed space EEE. In round t=1,…,Tt=1,\ldots,Tt=1,…,T, the learner outputs xt∈Vx_t\in Vxt​∈V, then receives a loss function ℓt\ell_tℓt​ and incurs ℓt(xt)\ell_t(x_t)ℓt​(xt​). For a fixed comparator u∈Vu\in Vu∈V, regret is ∑t=1T(ℓt(xt)−ℓt(u))\sum_{t=1}^T(\ell_t(x_t)-\ell_t(u))∑t=1T​(ℓt​(xt​)−ℓt​(u)). The comparison point stays the same across all rounds; it is not a round-by-round optimum.

FTRL supplies a sequence of regularizers ψt\psi_tψt​. At the start of round ttt, it selects any minimizer of the regularized cumulative loss

Ft(z)=ψt(z)+∑i=1t−1ℓi(z),xt∈argmin⁡z∈VFt(z).F_t(z)=\psi_t(z)+\sum_{i=1}^{t-1}\ell_i(z),\qquad x_t\in\operatorname*{argmin}_{z\in V}F_t(z).Ft​(z)=ψt​(z)+i=1∑t−1​ℓi​(z),xt​∈z∈Vargmin​Ft​(z).

The sum contains only past losses. Its empty value at t=1t=1t=1 means x1x_1x1​ minimizes ψ1\psi_1ψ1​. The analysis also names xT+1x_{T+1}xT+1​ and ψT+1\psi_{T+1}ψT+1​; that extra iterate occurs after the last charged round. A loss subgradient gtg_tgt​ at xtx_txt​ is a continuous linear functional with ℓt(y)≥ℓt(xt)+gt(y−xt)\ell_t(y)\ge\ell_t(x_t)+g_t(y-x_t)ℓt​(y)≥ℓt​(xt​)+gt​(y−xt​) for every y∈Ey\in Ey∈E. Its operator norm is the dual norm ∥gt∥∗\|g_t\|_*∥gt​∥∗​. The norm in the strong-convexity and stability statements is the norm on EEE.

Formalization targets

The goal is Corollary 7.6 on printed page 102. Let ψ\psiψ be closed and μ\muμ-strongly convex on VVV, with μ>0\mu>0μ>0, and let m=min⁡z∈Vψ(z)m=\min_{z\in V}\psi(z)m=minz∈V​ψ(z). The round-ttt regularizer is ψt(z)=(ψ(z)−m)/ηt−1\psi_t(z)=(\psi(z)-m)/\eta_{t-1}ψt​(z)=(ψ(z)−m)/ηt−1​, with ψT+1=ψT\psi_{T+1}=\psi_TψT+1​=ψT​. If the losses are subdifferentiable on VVV, the first target is the bound, for every u∈Vu\in Vu∈V and every choice of subgradients gt∈∂ℓt(xt)g_t\in\partial\ell_t(x_t)gt​∈∂ℓt​(xt​),

∑t=1T(ℓt(xt)−ℓt(u))≤ψ(u)−mηT−1+∑t=1Tηt−12μ∥gt∥∗2+∑t=1T−1(1ηt−1−1ηt)(ψ(xt+1)−m).\sum_{t=1}^T\bigl(\ell_t(x_t)-\ell_t(u)\bigr) \le \frac{\psi(u)-m}{\eta_{T-1}} +\sum_{t=1}^T\frac{\eta_{t-1}}{2\mu}\|g_t\|_*^2 +\sum_{t=1}^{T-1}\left(\frac1{\eta_{t-1}}-\frac1{\eta_t}\right) \bigl(\psi(x_{t+1})-m\bigr).t=1∑T​(ℓt​(xt​)−ℓt​(u))≤ηT−1​ψ(u)−m​+t=1∑T​2μηt−1​​∥gt​∥∗2​+t=1∑T−1​(ηt−1​1​−ηt​1​)(ψ(xt+1​)−m).

The last sum records the change in the regularizer; it is part of the stated result. The second target uses losses that are LLL-Lipschitz on an open neighborhood of VVV and the specified schedule ηt−1=αμ/(Lt)\eta_{t-1}=\alpha\sqrt\mu/(L\sqrt t)ηt−1​=αμ​/(Lt​), for α,L>0\alpha,L>0α,L>0:

∑t=1Tℓt(xt)−∑t=1Tℓt(u)≤(ψ(u)−mα+α)LTμ.\sum_{t=1}^T\ell_t(x_t)-\sum_{t=1}^T\ell_t(u) \le\left(\frac{\psi(u)-m}{\alpha}+\alpha\right) \frac{L\sqrt T}{\sqrt\mu}.t=1∑T​ℓt​(xt​)−t=1∑T​ℓt​(u)≤(αψ(u)−m​+α)μ​LT​​.

The goal records both conclusions. Its milestones are the exact regret identity of Lemma 7.1, the stability inequality and minimizer assertion of Lemma 7.5, the proximal-regularizer variant of Lemma 7.7, and the linear-perturbation stability result of Lemma 7.49. The last result extends the chapter's stability theme beyond the immediate proof of Corollary 7.6.

Significance

The second inequality gives an explicit T\sqrt TT​ bound with a comparator-dependent regularization term, without imposing a finite diameter on VVV. The first inequality retains the complete dependence on the chosen step sizes and on changes in the regularizer. A reader can therefore see what a schedule must control before specializing to the square-root choice. These are results already proved in the cited textbook; the open work in this mission is their Lean formalization, along with proofs of the stated milestones.

Formalizing the chapter would provide reusable statements about minimizers of strongly convex functions, subgradients in a dual normed space, and the stability of regularized minimization. The identities and stability bounds also connect FTRL's update rule to regret estimates without tying the interface to a single Euclidean norm or a single learning-rate choice. The definitions are intended to be usable by later chapters that employ linearized or optimistic FTRL, although those later results are outside this mission's goal.

Difficulty

The algebraic regret identity is exact, but it still contains the losses on both sides of the comparison. It does not by itself provide an explicit upper bound. The challenge is to control the change between successive minimizers while the objective and regularizer both change, then retain the correct index shift in the regularizer terms. In the general-norm setting, a subgradient is a dual vector; treating it as a primal vector would change the inequality's meaning. The corollary also has two quantifier patterns: its first bound covers every admissible subgradient sequence, while its second uses a particular learning-rate schedule and a Lipschitz condition.

Formalization scope

The Lean development uses a finite-dimensional real normed space. Losses and regularizers are real-valued functions on that space, and minimization is restricted to VVV. This is the chapter's finite-valued regime of its broader extended-real notation. A closed function on VVV is represented by the closed epigraph of its restriction. Strong convexity is Mathlib's StrongConvexOn; a subgradient is a continuous linear functional and its operator norm is the book's dual norm. Subdifferentiability is required at every feasible point. The Lipschitz condition is stated on an open set containing VVV, as on the page.

Indices start at 111 for rounds and at 000 for step sizes. The theorem requires T≥1T\ge1T≥1, positive step sizes, and positive values for every quantity used as a denominator. The minimum mmm is supplied by a point attaining it, and an FTRL run is a sequence of actual minimizers through round T+1T+1T+1. This rules out a sequence of arbitrary feasible decisions or a regularizer schedule unrelated to η\etaη. The first inequality retains the regularizer-change sum even when step sizes are not monotone. The book's one-line proof remarks that this change is nonpositive under an unstated monotonicity condition; the printed statement itself does not impose that condition.

The definition layer and all four milestone statements are in scope. Contributions proving any of them, or developing general lemmas about closed strongly convex objectives and dual-space subgradients that they need, fit this mission. The only simplification is replacing extended-real losses by real-valued losses on the feasible region, so no infinite value or undefined extended-real subtraction enters a displayed bound.

Selected references

  • Francesco Orabona, Online Learning: A Modern Introduction Using Convex Optimization, arXiv:1912.13213v10, 2026, Chapter 7, pp. 99–140. arXiv
7 thms3 active usersReviewed
🏆Completed
Convex OptimizationProbability·Captain: mikedeng1

A Modern Introduction to Online Learning II: Online-to-Batch Conversion Turns a Regret Bound into a Risk Bound with Deviation 2√(2 ln(2/δ)/T)Textbook

Motivation

An online learning algorithm chooses a decision before seeing each round's data and is judged by regret: how much more loss it incurs than a fixed decision chosen in hindsight. Statistical learning asks a different question. After seeing data drawn from a distribution, how small is the decision's expected loss on another draw? Chapter 3 of Orabona's textbook connects these settings. It shows that a deterministic online rule with a pathwise regret guarantee also controls expected loss when its input data are independent samples. The chapter's high-probability result is the target here.

The distinction matters when an online guarantee is already available. A regret theorem holds for every loss sequence, including adversarial sequences, yet it does not by itself state what a learner achieves on future data from a distribution. The conversion supplies that statistical reading. It applies to the sequence of predictions produced along one sample path, with no need to select an empirical-risk minimizer from the full decision set. Chapter 3 also studies finite-set validation, in which a separate sample selects one candidate from a fixed collection Orabona, Chapter 3.

Setting

The decision set VVV is a subset of Euclidean space Rd\mathbb R^dRd. A data point zzz belongs to a measurable space DDD and is drawn according to a probability distribution ρ\rhoρ. The loss f(x,z)f(x,z)f(x,z) lies in [0,1][0,1][0,1] for every x∈Vx\in Vx∈V and z∈Dz\in Dz∈D, and is measurable as a function of zzz when xxx is fixed. Its true risk is F(x)=Ez∼ρ[f(x,z)]F(x)=\mathbb E_{z\sim\rho}[f(x,z)]F(x)=Ez∼ρ​[f(x,z)]. For a sample ξ1,…,ξT\xi_1,\ldots,\xi_Tξ1​,…,ξT​ drawn independently from ρ\rhoρ, round ttt presents the loss function ℓt(x)=f(x,ξt)\ell_t(x)=f(x,\xi_t)ℓt​(x)=f(x,ξt​).

A deterministic online algorithm produces predictions x1,…,xT+1x_1,\ldots,x_{T+1}x1​,…,xT+1​ in VVV. Its prediction xtx_txt​ uses only ξ1,…,ξt−1\xi_1,\ldots,\xi_{t-1}ξ1​,…,ξt−1​. For a fixed competitor u∈Vu\in Vu∈V, its regret is Regret⁡T(u)=∑t=1Tℓt(xt)−∑t=1Tℓt(u)\operatorname{Regret}_T(u)=\sum_{t=1}^T\ell_t(x_t)-\sum_{t=1}^T\ell_t(u)RegretT​(u)=∑t=1T​ℓt​(xt​)−∑t=1T​ℓt​(u). A proposed regret bound R(u,T)R(u,T)R(u,T) is required to hold on every sample path and for every competitor uuu. The risk comparison uses a minimizer of the full expression F(u)+R(u,T)/TF(u)+R(u,T)/TF(u)+R(u,T)/T, whose existence the theorem assumes.

Formalization targets

The goal is Theorem 3.10, with T≥1T\ge1T≥1 and 0<δ<10<\delta<10<δ<1:

P ⁣{1T∑t=1TF(xt)≤min⁡u∈V(F(u)+R(u,T)T)+22ln⁡(2/δ)T}≥1−δ.\mathbb P\!\left\{\frac1T\sum_{t=1}^T F(x_t) \le \min_{u\in V}\left(F(u)+\frac{R(u,T)}{T}\right) +2\sqrt{\frac{2\ln(2/\delta)}{T}}\right\}\ge1-\delta.P{T1​t=1∑T​F(xt​)≤u∈Vmin​(F(u)+TR(u,T)​)+2T2ln(2/δ)​​}≥1−δ.

The minimum includes both risk and normalized regret. The theorem controls the average risk of the produced predictions; it does not require FFF to be convex or claim the risk of their average decision. This makes the result applicable to decision sets and losses beyond the convex setting of the chapter's expectation result.

The first milestone is the two-sided bounded-martingale inequality of Theorem 3.9. A second milestone is the displayed concentration estimate in the proof of Theorem 3.10, which compares average true risk with average observed loss for the adaptive predictions. The final milestone is Theorem 3.17: if a fixed nonempty finite set S⊆VS\subseteq VS⊆V is evaluated on fresh data and x^\widehat xx minimizes its empirical risk, then

F(x^)≤min⁡x∈SF(x)+22ln⁡(2∣S∣/δ)TF(\widehat x)\le\min_{x\in S}F(x) +2\sqrt{\frac{2\ln(2|S|/\delta)}{T}}F(x)≤x∈Smin​F(x)+2T2ln(2∣S∣/δ)​​

with probability at least 1−δ1-\delta1−δ Orabona, pp. 27–32.

Significance

The goal gives a direct way to read an online regret certificate as a distributional risk certificate. The bound retains the competitor-dependent function R(u,T)R(u,T)R(u,T) inside the minimization: a decision with slightly larger risk but substantially smaller regret bound may give the better guarantee. Its deviation term has the exact constant and confidence dependence stated in the textbook. The finite-set theorem supplies a distinct validation guarantee for choosing one candidate from a fixed collection; it applies even when the risk is not convex.

Formalizing these results adds the probability layer absent from a deterministic regret statement: independent samples, a prediction rule constrained by past observations, true risk as an integral, and probability events for the adaptive sequence. Related online-to-batch statements already exist in the platform's developments of Hazan's online convex optimization text and Mohri, Rostamizadeh and Talwalkar's Foundations of Machine Learning. They use different algorithm and loss models and do not state this competitor-dependent bound. The present mission records the precise Chapter 3 formulation as an open formal proof target; the local Lean statements compile, but their proofs are not supplied.

Difficulty

The observed loss f(xt,ξt)f(x_t,\xi_t)f(xt​,ξt​) uses a prediction chosen from earlier observations, while the event in the goal concerns its expected loss F(xt)F(x_t)F(xt​). Treating the predictions as fixed would ignore their dependence on the sample history. Allowing a prediction to inspect the current sample would break the comparison. The required concentration statement therefore has to retain the exact information available at each round. The comparator in the final bound is fixed before the samples are observed, while the online predictions remain sample dependent Orabona, proof of Theorem 3.10.

For the finite-set result, the selected predictor depends on the validation sample. A concentration estimate for a single fixed predictor does not immediately give the asserted guarantee for that selected predictor. The finite set must remain fixed before sampling, and the statement must account for every possible selected member.

Formalization scope

The Lean decision space is EuclideanSpace ℝ (Fin d). Samples are functions N→D\mathbb N\to DN→D under the infinite product measure of ρ\rhoρ; their values at index zero are unused. The algorithm is a function of the full sample path constrained by a non-anticipation predicate: changing samples from round ttt onward cannot change xtx_txt​. Regret is against a fixed uuu, with sums over rounds 1,…,T1,\ldots,T1,…,T. The pathwise regret guarantee quantifies over every sample path, so it cannot be replaced by an expected-regret assumption.

The integral defining FFF is used with the book's measurable [0,1][0,1][0,1] losses on VVV, avoiding an undefined-integral convention. Measurability of each prediction, realized loss, and realized risk is stated explicitly where a probability event uses them. These regularity conditions make the probability statements meaningful and are the only extra conditions beyond the textbook's mathematical setup. A minimizer u⋆∈Vu_\star\in Vu⋆​∈V represents the assumed minimum in Theorem 3.10. The finite-set minimum is a real infimum whose input is nonempty, finite, and bounded. The positive horizon guards division by TTT; 0<δ<10<\delta<10<δ<1 guards the logarithms. In Theorem 3.9, B>0B>0B>0 and ε≥0\varepsilon\ge0ε≥0 make its printed tail formula valid, including at the threshold boundary.

A complete development needs Mathlib's infinite product measure, filtration and martingale interfaces, conditional expectation or conditional moment bounds, and the sub-Gaussian concentration library. The risk, regret, non-anticipation, and empirical-risk definitions can be reused for later stochastic learning results. Contributions that establish the concentration milestone, connect it to the regret identity, or prove the finite-set validation statement all fit this scope.

Selected references

  • Francesco Orabona, Online Learning: A Modern Introduction Using Convex Optimization, arXiv:1912.13213v10, 2026. arXiv
  • Elad Hazan, Introduction to Online Convex Optimization, Foundations and Trends in Optimization, 2016; revised arXiv version, 2019. arXiv
5 thms3 active usersReviewed
🏆Completed
OptimizationProbability·Captain: mikedeng1

Linear Convergence of Gradient and Proximal-Gradient Methods Under the Polyak-Łojasiewicz Condition III: Randomized Coordinate Descent Has Expected Rate (1 − μ/(dL))^kResearch Paper

Motivation

Coordinate descent changes one component of a parameter vector at a time. This matters when evaluating or updating a single component is cheaper than handling the full vector, for example in problems with sparse features. The coordinate chosen at each step can be random or selected from the current gradient. Karimi, Nutini, and Schmidt study these choices for objectives that satisfy a Polyak–Łojasiewicz (PL) inequality, a condition that controls objective error by the size of the gradient without assuming convexity. Their Theorem 3 gives an explicit expected rate for uniform random coordinates, while §3.2 gives a rate for the greedy Gauss–Southwell choice. Karimi, Nutini, and Schmidt (2016), §3

The question is how much of the familiar full-gradient linear rate survives when an update sees only one coordinate. The theorem answers this with a dimension factor: each uniform draw accesses one of ddd coordinates, and the rate depends on μ/(dL)\mu/(dL)μ/(dL). The same paper relates the greedy rule to earlier work on coordinate descent and boosting; this mission treats its stated greedy rate as a companion result, using the same coordinate update. Karimi, Nutini, and Schmidt (2016), §§3.1–3.2

Setting

Let d≥1d\ge1d≥1, and let f:Rd→Rf:\mathbb R^d\to\mathbb Rf:Rd→R be differentiable. Write ∇if(x)\nabla_i f(x)∇i​f(x) for coordinate iii of its gradient and eie_iei​ for the corresponding standard basis vector. The paper's coordinate-wise upper bound (5), with L>0L>0L>0, is

f(x+αei)≤f(x)+α∇if(x)+L2α2for every x∈Rd, i∈{1,…,d}, α∈R.f(x+\alpha e_i)\le f(x)+\alpha\nabla_i f(x)+\frac L2\alpha^2 \quad\text{for every }x\in\mathbb R^d,\ i\in\{1,\ldots,d\},\ \alpha\in\mathbb R.f(x+αei​)≤f(x)+α∇i​f(x)+2L​α2for every x∈Rd, i∈{1,…,d}, α∈R.

Here iii ranges over the ddd coordinates. The bound concerns motion along one coordinate. It is the exact smoothness condition used in Theorem 3. Let x∗x^*x∗ be a global minimizer and put f∗=f(x∗)f^*=f(x^*)f∗=f(x∗); thus the solution set is nonempty and f∗f^*f∗ is an attained optimal value. The PL inequality (3) with μ>0\mu>0μ>0 says

12∥∇f(x)∥2≥μ(f(x)−f∗)for every x∈Rd.\frac12\|\nabla f(x)\|^2\ge\mu\bigl(f(x)-f^*\bigr) \quad\text{for every }x\in\mathbb R^d.21​∥∇f(x)∥2≥μ(f(x)−f∗)for every x∈Rd.

The randomized coordinate-descent step (6) draws iki_kik​ uniformly from the ddd coordinates, independently at each time, and sets xk+1=xk−L−1∇ikf(xk)eikx_{k+1}=x_k-L^{-1}\nabla_{i_k}f(x_k)e_{i_k}xk+1​=xk​−L−1∇ik​​f(xk​)eik​​. The initial point x0x_0x0​ is arbitrary. For kkk draws, expectation means the average over all dkd^kdk words (i0,…,ik−1)(i_0,\ldots,i_{k-1})(i0​,…,ik−1​), each assigned weight d−kd^{-k}d−k. In particular, at k=0k=0k=0 the empty word has weight one and the iterate is x0x_0x0​.

The Gauss–Southwell rule instead chooses a coordinate with the largest ∣∇if(xk)∣|\nabla_i f(x_k)|∣∇i​f(xk​)∣; ties may be resolved in any way. Its rate uses a PL inequality in the maximum-coordinate norm, ∥∇f(x)∥∞=max⁡i∣∇if(x)∣\|\nabla f(x)\|_\infty=\max_i|\nabla_i f(x)|∥∇f(x)∥∞​=maxi​∣∇i​f(x)∣, with constant μ1>0\mu_1>0μ1​>0. Karimi, Nutini, and Schmidt (2016), pp. 5–6

Formalization targets

Theorem 3: uniform random coordinates

For every k≥0k\ge0k≥0, the goal is the paper's expected suboptimality bound

E[f(xk)−f∗]≤(1−μdL)k(f(x0)−f∗).\mathbb E[f(x_k)-f^*] \le\left(1-\frac{\mu}{dL}\right)^k\bigl(f(x_0)-f^*\bigr).E[f(xk​)−f∗]≤(1−dLμ​)k(f(x0​)−f∗).

The target includes the specific update (6), the coordinate condition (5), and the PL inequality (3). The three milestones are the paper's displayed one-coordinate decrease, its uniform expected decrease involving ∥∇f(x)∥2\|\nabla f(x)\|^2∥∇f(x)∥2, and the resulting conditional expected gap contraction. They are indexed to the proof of Theorem 3 on page 5. Karimi, Nutini, and Schmidt (2016), Theorem 3 and proof

Greedy coordinate choice

The companion target is the rate stated in §3.2 when the maximum-norm PL inequality holds:

f(xk)−f∗≤(1−μ1L)k(f(x0)−f∗).f(x_k)-f^* \le\left(1-\frac{\mu_1}{L}\right)^k\bigl(f(x_0)-f^*\bigr).f(xk​)−f∗≤(1−Lμ1​​)k(f(x0​)−f∗).

This bound applies to any run that chooses a maximizing coordinate at every step. The different PL modulus is part of the statement; it need not equal the modulus μ\muμ in Theorem 3. Karimi, Nutini, and Schmidt (2016), §3.2

Significance

Theorem 3 gives a quantitative guarantee for a method that reads one gradient component per update. The guarantee depends on the attained optimum, the coordinate smoothness constant, the PL modulus, and the number of coordinates. It applies to objectives satisfying PL even when convexity is unavailable. The companion result records how the stated rate changes when the coordinate is selected by the largest gradient component and the PL condition is measured in the corresponding dual norm. Karimi, Nutini, and Schmidt (2016), §3

The paper proves these claims, but this mission asks for machine-checked proofs of their exact hypotheses and constants. A complete development will make the finite random experiment explicit and give reusable definitions for coordinate updates and expectations of finite coordinate-choice histories. The PL predicate is already published in the theorem library; the coordinate model and the results here remain to be formalized.

Difficulty

A bound for a single coordinate update depends on one component of the gradient. An arbitrary coordinate can have a zero component even when the full gradient is large, so a pointwise decrease with the full gradient norm cannot be asserted for every choice. The theorem instead concerns the expectation over a uniform choice and then over the complete history of choices. The main formal difficulty is keeping those two levels of expectation aligned with a state xkx_kxk​ that itself depends on earlier draws. For the greedy companion, the selected coordinate may not be unique; the statement must cover every maximizing choice without assuming a fixed tie-breaking rule.

Formalization scope

The vector space is EuclideanSpace ℝ (Fin d). The shared definitions give (5), the step (6), the iterate produced by a word Fin k → Fin d, and the expected gap as a real finite sum. The published ShiOptRates.PolyakLojasiewicz predicate represents (3) after multiplying both sides by two. Every theorem using uniform draws assumes d>0d>0d>0; every step with 1/L1/L1/L assumes L>0L>0L>0. The goal also states differentiability, an explicit global minimizer, and μ>0\mu>0μ>0. No convexity is assumed. A free real number called f∗f^*f∗ would change the PL condition and could invalidate the rate, so the formal statement uses f(x∗)f(x^*)f(x∗) throughout.

The index-word average is an exact finite representation of independent uniform sampling, including k=0k=0k=0. It does not sample coordinates without replacement. The companion theorem represents the maximum-norm PL condition by requiring its inequality at every maximizing coordinate, which is equivalent for d>0d>0d>0, and quantifies over every run satisfying the Gauss–Southwell update relation. Useful contributions include proofs of the three displayed milestones, the passage from one draw to the word average, and the greedy companion. The paper's optional remark on sampling proportional to coordinate-specific constants LiL_iLi​ and its exact-coordinate-minimization remark are outside this proposal.

Selected references

  • H. Karimi, J. Nutini, and M. Schmidt, Linear Convergence of Gradient and Proximal-Gradient Methods Under the Polyak–Łojasiewicz Condition, ECML PKDD, 2016; arXiv:1608.04636v4.
7 thms3 active usersReviewed
🏆Completed
Convex OptimizationOptimization·Captain: mikedeng1

Linear Convergence of Gradient and Proximal-Gradient Methods Under the Polyak-Łojasiewicz Condition II: SC → ESC → WSC → RSI → EB ≡ PL → QG, and RSI ≡ EB ≡ PL ≡ QG for Convex fResearch Paper

Motivation

Gradient descent on a smooth function f:Rd→Rf:\mathbb R^d\to\mathbb Rf:Rd→R converges at a linear rate (the optimality gap shrinks by a constant factor each step) under strong convexity. Many problems in machine learning and statistics are not strongly convex: least squares with a rank-deficient design, logistic regression on separable or degenerate data, and some non-convex objectives. Over roughly 25 years several weaker conditions were proposed that still give linear convergence. Examples are the error bound of Luo and Tseng, restricted strong convexity, essential strong convexity, weak strong convexity, and quadratic growth (also called optimal or semi-strong convexity). Each came with its own proof, usually assuming convexity. How these conditions relate to one another was rarely made explicit.

Karimi, Nutini and Schmidt (arXiv:1608.04636, ECML PKDD 2016) put the Polyak–Łojasiewicz (PL) inequality (Polyak, 1963) at the centre. Their Theorem 2 orders seven conditions by strength for functions with a Lipschitz gradient, without convexity. It shows that PL is equivalent to the error bound and implies quadratic growth. For convex functions it shows that the four weakest conditions coincide. The equivalence of EB and PL is a special case of a result of Bolte, Nguyen, Peypouquet and Suter (arXiv:1510.08234, Theorem 5). Zhang (arXiv:1606.00269) independently related RSI, EB, PL and QG.

Setting

Let f:Rd→Rf:\mathbb R^d\to\mathbb Rf:Rd→R be differentiable with gradient ∇f\nabla f∇f. The paper's standing assumptions (§2, p. 2) are:

  1. ∇f\nabla f∇f is LLL-Lipschitz: ∥∇f(x)−∇f(y)∥≤L∥x−y∥\|\nabla f(x)-\nabla f(y)\|\le L\|x-y\|∥∇f(x)−∇f(y)∥≤L∥x−y∥ for all x,yx,yx,y, with L>0L>0L>0;
  2. the solution set X∗={x: f(x)≤f(y) ∀y}\mathcal X^*=\{x:\ f(x)\le f(y)\ \forall y\}X∗={x: f(x)≤f(y) ∀y} is non-empty. Write f∗f^*f∗ for the optimal value.

For every xxx, xpx_pxp​ denotes a projection of xxx onto X∗\mathcal X^*X∗: a point of X∗\mathcal X^*X∗ at distance dist⁡(x,X∗)\operatorname{dist}(x,\mathcal X^*)dist(x,X∗) from xxx. Each of the following conditions involves a constant μ>0\mu>0μ>0, which may differ between conditions (Appendix A, pp. 11–12):

  • SC (strong convexity): f(y)≥f(x)+⟨∇f(x),y−x⟩+μ2∥y−x∥2f(y)\ge f(x)+\langle\nabla f(x),y-x\rangle+\frac\mu2\|y-x\|^2f(y)≥f(x)+⟨∇f(x),y−x⟩+2μ​∥y−x∥2 for all x,yx,yx,y.
  • ESC (essential strong convexity): the same inequality for all x,yx,yx,y with xp=ypx_p=y_pxp​=yp​.
  • WSC (weak strong convexity): f∗≥f(x)+⟨∇f(x),xp−x⟩+μ2∥xp−x∥2f^*\ge f(x)+\langle\nabla f(x),x_p-x\rangle+\frac\mu2\|x_p-x\|^2f∗≥f(x)+⟨∇f(x),xp​−x⟩+2μ​∥xp​−x∥2 for all xxx.
  • RSI (restricted secant inequality): ⟨∇f(x),x−xp⟩≥μ∥xp−x∥2\langle\nabla f(x),x-x_p\rangle\ge\mu\|x_p-x\|^2⟨∇f(x),x−xp​⟩≥μ∥xp​−x∥2 for all xxx.
  • EB (error bound): ∥∇f(x)∥≥μ∥xp−x∥\|\nabla f(x)\|\ge\mu\|x_p-x\|∥∇f(x)∥≥μ∥xp​−x∥ for all xxx.
  • PL (Polyak–Łojasiewicz): 12∥∇f(x)∥2≥μ(f(x)−f∗)\frac12\|\nabla f(x)\|^2\ge\mu(f(x)-f^*)21​∥∇f(x)∥2≥μ(f(x)−f∗) for all xxx.
  • QG (quadratic growth): f(x)−f∗≥μ2∥xp−x∥2f(x)-f^*\ge\frac\mu2\|x_p-x\|^2f(x)−f∗≥2μ​∥xp​−x∥2 for all xxx.

Formalization targets

Goal: Theorem 2 (p. 3)

Read each condition as "it holds for some μ>0\mu>0μ>0". For fff with a Lipschitz-continuous gradient,

(SC)→(ESC)→(WSC)→(RSI)→(EB)≡(PL)→(QG),(SC)\to(ESC)\to(WSC)\to(RSI)\to(EB)\equiv(PL)\to(QG),(SC)→(ESC)→(WSC)→(RSI)→(EB)≡(PL)→(QG),

and if fff is also convex,

(RSI)≡(EB)≡(PL)≡(QG).(RSI)\equiv(EB)\equiv(PL)\equiv(QG).(RSI)≡(EB)≡(PL)≡(QG).

The goal states no constants, because the constant may change along the chain.

Milestones: the implications of Appendix A, with their constants

The appendix proves a subset of the arrows, each with an explicit constant. The rest follow by transitivity.

SC(μ)⇒ESC(μ),ESC(μ)⇒WSC(μ),WSC(μ)⇒RSI(μ/2),RSI(μ)⇒EB(μ),\text{SC}(\mu)\Rightarrow\text{ESC}(\mu),\quad \text{ESC}(\mu)\Rightarrow\text{WSC}(\mu),\quad \text{WSC}(\mu)\Rightarrow\text{RSI}(\mu/2),\quad \text{RSI}(\mu)\Rightarrow\text{EB}(\mu),SC(μ)⇒ESC(μ),ESC(μ)⇒WSC(μ),WSC(μ)⇒RSI(μ/2),RSI(μ)⇒EB(μ), EB(μ)⇒PL(μ2/L),PL(μ)⇒QG(μ),PL(μ)⇒EB(μ),QG(μ)+convex⇒RSI(μ/2).\text{EB}(\mu)\Rightarrow\text{PL}(\mu^2/L),\quad \text{PL}(\mu)\Rightarrow\text{QG}(\mu),\quad \text{PL}(\mu)\Rightarrow\text{EB}(\mu),\quad \text{QG}(\mu)+\text{convex}\Rightarrow\text{RSI}(\mu/2).EB(μ)⇒PL(μ2/L),PL(μ)⇒QG(μ),PL(μ)⇒EB(μ),QG(μ)+convex⇒RSI(μ/2).

One further milestone is display (20). Under PL, g=f−f∗g=\sqrt{f-f^*}g=f−f∗​ satisfies ∥∇g(x)∥≥μ/2\|\nabla g(x)\|\ge\sqrt{\mu/2}∥∇g(x)∥≥μ/2​ at every x∉X∗x\notin\mathcal X^*x∈/X∗.

Significance

The result. Theorem 2 shows that QG is the weakest of the seven conditions. QG alone does not rule out non-global local minima, so it cannot guarantee that gradient descent reaches a global minimizer. Among the conditions considered, PL and the equivalent EB are therefore the most general that still give linear convergence to a global minimum. The paper's Theorem 1 gives that rate under PL with a one-paragraph proof. For convex fff, RSI (there called restricted strong convexity), EB, PL and QG are one condition. A rate proved under any of them transfers to the others, up to the constants of the milestones.

Formalizing it. The result is proved on paper, and the arrows of Appendix A are short, except PL → QG. A machine-checked version adds three things. It fixes a precise meaning for xpx_pxp​ when X∗\mathcal X^*X∗ is not convex and a nearest point is not unique. It records the correct constant in EB → PL: the printed text says μ/L\mu/Lμ/L, while the displayed argument yields μ2/L\mu^2/Lμ2/L. It puts PL → QG for non-convex fff on a checked footing; the acknowledgments note an earlier error in this proof. No formalization of these implications is known to be on the platform or in Mathlib.

Difficulty

Most arrows are one-line rearrangements or the Cauchy–Schwarz inequality. EB → PL also needs the quadratic upper bound implied by the Lipschitz gradient, and ∇f(xp)=0\nabla f(x_p)=0∇f(xp​)=0 at a global minimizer.

The obstacle is PL → QG without convexity. The natural first idea fails. Integrating along the segment from xxx to xpx_pxp​ uses no information from PL, because the gradient along that segment is unrelated to the gap f−f∗f-f^*f−f∗. The implication is global: it converts a pointwise bound on the gradient into a lower bound on fff in terms of the distance to a possibly non-convex set X∗\mathcal X^*X∗. Any proof has to bound the distance to X∗\mathcal X^*X∗ by the length of some path into it that is controlled by fff. Building such a path in Lean, with its existence, its arrival at X∗\mathcal X^*X∗ in finite time and a bound on its length, needs analysis that is not yet assembled in Mathlib for this purpose.

Formalization scope

  • Space. Rd\mathbb R^dRd is EuclideanSpace ℝ (Fin d), abbreviated E d. Gradients are Mathlib's gradient. The Lipschitz gradient is the literal hypothesis ∀ x y, ‖gradient f x - gradient f y‖ ≤ L * ‖x - y‖ with 0 < L, together with Differentiable ℝ f.
  • Solution set. Non-emptiness is carried by a global minimizer xstar with ∀ y, f xstar ≤ f y, and f∗f^*f∗ is f xstar. In WSC, f∗f^*f∗ is written f(xp)f(x_p)f(xp​).
  • Projection. xpx_pxp​ is an explicit map p with IsProjSel f p: p x lies in X∗\mathcal X^*X∗ and ‖x - p x‖ = Metric.infDist x X*. Every condition, and Theorem 2, refers to an arbitrary such selection. ESC's "xp=ypx_p=y_pxp​=yp​" is p x = p y for that same selection. No Classical.choose is used.
  • PL is the published predicate ShiOptRates.PolyakLojasiewicz f (gradient f) μ fstar (definition ShiOptrates_FunctionClass), i.e. 2μ(f(x)−f∗)≤∥∇f(x)∥22\mu(f(x)-f^*)\le\|\nabla f(x)\|^22μ(f(x)−f∗)≤∥∇f(x)∥2. This is the paper's (3) multiplied by 2.
  • Constants. Each milestone carries the constant of its appendix bullet and the hypothesis μ>0\mu>0μ>0. The goal uses Holds C := ∃ μ, 0 < μ ∧ C μ independently on each side of every arrow. Convexity is ConvexOn ℝ Set.univ f.
  • Dropped standing hypotheses. Milestones whose argument does not use the Lipschitz gradient omit it: SC → ESC, ESC → WSC, WSC → RSI, RSI → EB, QG + convex → RSI, and (20). Each is then a slightly stronger statement.
  • Corrected constant. EB → PL is stated with μ2/L\mu^2/Lμ2/L, which is what the paper's displayed argument proves. The quoted milestone text keeps the printed "μ/L\mu/Lμ/L".
  • No trivialization. Every condition is quantified over all xxx (and yyy). The selection p is constrained by IsProjSel, not free, and the goal's hypotheses are satisfiable: f(x)=12∥x∥2f(x)=\frac12\|x\|^2f(x)=21​∥x∥2 with p≡0p\equiv0p≡0, L=1L=1L=1 satisfies all of them and SC with μ=1\mu=1μ=1. A formalization that fixed one μ\muμ for the whole chain, or replaced xpx_pxp​ by an unconstrained map, would be a different theorem.
  • Reuse. The condition predicates and the projection selection are reusable by any development on linear convergence without strong convexity. A general "PL implies QG" lemma via gradient flow is reusable beyond this mission. Contributions toward the ODE and arc-length infrastructure are welcome.

Selected references

  • H. Karimi, J. Nutini, M. Schmidt, Linear Convergence of Gradient and Proximal-Gradient Methods Under the Polyak-Łojasiewicz Condition, ECML PKDD 2016; arXiv:1608.04636v4 (2020). https://arxiv.org/abs/1608.04636
  • B. T. Polyak, Gradient methods for minimizing functionals, Zh. Vychisl. Mat. Mat. Fiz. 3(4):643–653, 1963 (English translation: USSR Comput. Math. Math. Phys. 3(4):864–878). https://doi.org/10.1016/0041-5553(63)90382-3
  • J. Bolte, T. P. Nguyen, J. Peypouquet, B. W. Suter, From error bounds to the complexity of first-order descent methods for convex functions, Math. Program. 165, 2017. https://arxiv.org/abs/1510.08234
  • H. Zhang, The restricted strong convexity revisited: analysis of equivalence to error bound and quadratic growth, Optim. Lett. 11, 2017. https://arxiv.org/abs/1511.01635
  • H. Zhang, New analysis of linear convergence of gradient-type methods via unifying error bound conditions, Math. Program. 180, 2020. https://arxiv.org/abs/1606.00269
  • Z.-Q. Luo, P. Tseng, Error bounds and convergence analysis of feasible descent methods: a general approach, Ann. Oper. Res. 46, 1993. https://doi.org/10.1007/BF02096261
12 thms3 active usersReviewed
🏆Completed
Optimization·Captain: mikedeng1

Linear Convergence of Gradient and Proximal-Gradient Methods Under the Polyak-Łojasiewicz Condition I: Gradient Descent with Step 1/L Satisfies f(x_k) − f* ≤ (1 − μ/L)^k (f(x_0) − f*)Research Paper

Motivation

Gradient descent is often introduced with strong convexity as the condition that gives a geometric reduction in objective error. That condition excludes useful smooth objectives whose minimizer is not unique and some nonconvex objectives. Karimi, Nutini, and Schmidt place the Polyak–Łojasiewicz inequality at the center of a more general account: it controls the objective gap using the squared gradient norm without requiring convexity. Their paper uses it to analyze gradient, coordinate, stochastic, and proximal gradient methods, and identifies settings including least squares and logistic regression in which the condition is useful. Karimi, Nutini, and Schmidt, §1–2

This mission isolates the paper's first rate theorem and its immediate analytic setting. The question is whether the ordinary gradient step with length 1/L1/L1/L has a global geometric rate when the objective has a Lipschitz gradient, attains its minimum, and satisfies the PL inequality. The claim applies to every initial point, even when the function is nonconvex. It also gives a precise test case for deciding which properties of smoothness and optimality must be represented explicitly in Lean. Karimi, Nutini, and Schmidt, Theorem 1, p. 2

Setting

Let ddd be a natural number and let f:Rd→Rf:\mathbb R^d\to\mathbb Rf:Rd→R be a differentiable objective. Its gradient at xxx is ∇f(x)\nabla f(x)∇f(x), with the Euclidean norm and inner product denoted by ∥⋅∥\|\cdot\|∥⋅∥ and ⟨⋅,⋅⟩\langle\cdot,\cdot\rangle⟨⋅,⋅⟩. An LLL-Lipschitz gradient means that ∥∇f(x)−∇f(y)∥≤L∥x−y∥\|\nabla f(x)-\nabla f(y)\|\le L\|x-y\|∥∇f(x)−∇f(y)∥≤L∥x−y∥ for every x,yx,yx,y, with L>0L>0L>0. The corresponding quadratic upper bound is the paper's inequality (2):

f(y)≤f(x)+⟨∇f(x),y−x⟩+L2∥y−x∥2.f(y)\le f(x)+\langle\nabla f(x),y-x\rangle+\frac L2\|y-x\|^2.f(y)≤f(x)+⟨∇f(x),y−x⟩+2L​∥y−x∥2.

The optimization problem is unconstrained. It has a nonempty set of solutions: choose a global minimizer x∗x^*x∗, so f(x∗)≤f(y)f(x^*)\le f(y)f(x∗)≤f(y) for every yyy, and write f∗=f(x∗)f^*=f(x^*)f∗=f(x∗). This notation ties the reference value to an attained optimum. The PL inequality with constant μ>0\mu>0μ>0 is

12∥∇f(x)∥2≥μ(f(x)−f∗)for every x∈Rd.\frac12\|\nabla f(x)\|^2\ge\mu\bigl(f(x)-f^*\bigr) \qquad\text{for every }x\in\mathbb R^d.21​∥∇f(x)∥2≥μ(f(x)−f∗)for every x∈Rd.

Given any starting point x0x_0x0​, the gradient descent run is the sequence specified by xk+1=xk−L−1∇f(xk)x_{k+1}=x_k-L^{-1}\nabla f(x_k)xk+1​=xk​−L−1∇f(xk​). The term “linear convergence” here refers to a geometric bound on the difference in objective values, not to a claim that the points xkx_kxk​ approach a particular minimizer. Multiple minimizers are allowed. These conventions are those of §2 and Theorem 1 of the source. Karimi, Nutini, and Schmidt, p. 2

Formalization targets

Gradient descent at step size 1/L1/L1/L

The mission goal is the paper's Theorem 1. For every starting point and every k∈Nk\in\mathbb Nk∈N,

f(xk)−f∗≤(1−μL)k(f(x0)−f∗).f(x_k)-f^*\le\left(1-\frac\mu L\right)^k\bigl(f(x_0)-f^*\bigr).f(xk​)−f∗≤(1−Lμ​)k(f(x0​)−f∗).

The target includes the gradient update itself. Its three milestones record the source's analytic path: the quadratic upper bound (2) from the Lipschitz gradient assumption, the one-step decrease displayed in the proof, and the one-step PL contraction displayed immediately afterwards. The latter two are linked to proved platform theorems whose exact-gradient specializations give the source's displays. Karimi, Nutini, and Schmidt, Theorem 1 and proof, p. 2

Nearby consequences

The mission also states the paper's observation that an exact line search over all real step sizes has the same function-value rate. Two companion claims record that a stationary point under PL is a global minimizer, and that μ\muμ-strong convexity implies PL with the same μ\muμ. These claims locate the rate theorem among the conditions discussed in §2.2 and Appendix B; they are separate from the goal and its milestone list. Karimi, Nutini, and Schmidt, pp. 3, 14

Significance

The rate bounds the objective error at each finite iteration, with an explicit factor 1−μ/L1-\mu/L1−μ/L. It explains why uniqueness of the minimizer and convexity are not prerequisites for this particular function-value guarantee. It does not by itself assert a rate for the distance to the solution set, a rate for individual iterates, or convergence to a uniquely determined point. The strong-convexity companion shows that the classical strongly convex case fits inside the same PL framework. Karimi, Nutini, and Schmidt, §§2–2.2 and Appendix B

For formalization, the mission adds the precise bridge from a literal Lipschitz-gradient hypothesis to the quadratic upper-bound form used by the existing step results. It also states the complete source theorem with an attained optimum, a genuine gradient, and the update rule, rather than treating any reference value or gradient-like map as interchangeable with those objects. The two published one-step theorems are already machine checked on Prove2Me; the five source-specific draft statements in this proposal remain targets until their proofs are supplied and verified. This distinction matters because the previously published geometric-rate theorem assumes μ≤L\mu\le Lμ≤L, while Theorem 1 does not state that restriction.

Difficulty

Smoothness controls the change in fff during a gradient step through the size of the gradient. A geometric rate needs the change to be related to the remaining objective gap. That relation does not follow from smoothness alone; the PL inequality supplies the additional control. Another issue is that the factor 1−μ/L1-\mu/L1−μ/L may be negative if μ>L\mu>Lμ>L. The source does not exclude that parameter range, so the formal theorem must handle it without inserting an extra hypothesis. An exact line search raises a distinct issue: its asserted minimum along a line must be attained, and the chosen next iterate must be a point realizing it. Karimi, Nutini, and Schmidt, Theorem 1 and following paragraph, pp. 2–3

Formalization scope

Lean represents Rd\mathbb R^dRd as EuclideanSpace ℝ (Fin d), so d=0d=0d=0 is included and denotes a one-point space. It uses Mathlib's gradient on an explicitly differentiable fff, ordinary real scalars and norms, and a sequence indexed by natural numbers. The fixed step is written with 1/L1/L1/L and requires L>0L>0L>0. The minimizer witness x∗x^*x∗ and its global-minimality property are explicit; the PL predicate uses exactly f(x∗)f(x^*)f(x∗) as its reference value. The PL inequality is represented in the equivalent form 2μ(f(x)−f(x∗))≤∥∇f(x)∥22\mu(f(x)-f(x^*))\le\|\nabla f(x)\|^22μ(f(x)−f(x∗))≤∥∇f(x)∥2. The goal assumes the literal Lipschitz bound on gradients; equation (2) is a milestone conclusion. No convexity and no condition μ≤L\mu\le Lμ≤L are added.

The exact line-search companion quantifies over every real step size and requires an attained minimizing iterate. The stationary-point consequence retains differentiability but does not need the LLL-Lipschitz gradient bound. The strong-convexity companion uses the first inequality in Appendix A and the implication in Appendix B, with the same μ\muμ on both sides. The shared ShiOptRates predicates for the quadratic upper bound and PL inequality are referenced as published definitions. Useful contributions include proofs of the Lipschitz-to-quadratic-bound bridge and the paper-specific results, as well as reusable facts about gradients on Euclidean spaces. A freely chosen f∗f^*f∗ or a default value of gradient at a nondifferentiable point would change the mathematical claim; the goal excludes both.

Selected references

  • H. Karimi, J. Nutini, and M. Schmidt, Linear Convergence of Gradient and Proximal-Gradient Methods Under the Polyak–Łojasiewicz Condition, ECML PKDD, 2016; arXiv:1608.04636v4. This mission cites the pinned preprint's pp. 2–3, 11, and 14.
5 thms3 active usersReviewed
🏆Completed
OptimizationProbability·Captain: mikedeng1

Linear Convergence of Gradient and Proximal-Gradient Methods Under the Polyak-Łojasiewicz Condition IV: Stochastic Gradient Rates LC²/(2kμ²) and (1 − 2μα)^k [f(x_0) − f*] + LC²α/(4μ)Research Paper

Motivation

Stochastic gradient (SG) methods are the default algorithms for training machine-learning models: when an objective is an average over a very large number of training examples, computing its full gradient even a few times is infeasible, while the gradient of a single randomly chosen example is cheap. The classical convergence theory of SG rates its expected suboptimality under strong convexity — the O(1/k)O(1/k)O(1/k) rate of Nemirovski, Juditsky, Lan and Shapiro (2009) — yet common models such as least squares with a rank-deficient design and logistic regression are convex but not strongly convex.

Karimi, Nutini and Schmidt (2016) observed that the much older Polyak–Łojasiewicz (PL) inequality of Polyak (1963) is enough to recover linear rates for gradient descent and its variants without any convexity. Their Theorem 4 transfers this observation to SG: under PL alone, SG with a decreasing step attains the O(1/k)O(1/k)O(1/k) rate known from the strongly convex case, and SG with a constant step converges linearly to a neighbourhood of the optimum whose size is proportional to the step. This mission formalizes that theorem.

Setting

Let d≥0d\ge0d≥0 and let Rd\mathbb R^dRd carry the Euclidean norm ∥⋅∥\|\cdot\|∥⋅∥. An index iii ranges over a measurable space III with a probability measure ν\nuν, and each index has a differentiable component function fi:Rd→Rf_i:\mathbb R^d\to\mathbb Rfi​:Rd→R. The stochastic optimization problem (7) is to minimize

f(x)=Ei∼ν[fi(x)],f(x) = \mathbb E_{i\sim\nu}[f_i(x)],f(x)=Ei∼ν​[fi​(x)],

of which the finite sum f=1n∑i=1nfif = \frac1n\sum_{i=1}^n f_if=n1​∑i=1n​fi​ (8) is the case of the uniform distribution on nnn indices. The component gradients are unbiased: Ei∼ν[∇fi(x)]=∇f(x)\mathbb E_{i\sim\nu}[\nabla f_i(x)] = \nabla f(x)Ei∼ν​[∇fi​(x)]=∇f(x) for every xxx.

The SG method (9) starts from a point x0x_0x0​ and, given step sizes αk>0\alpha_k>0αk​>0 and indices i0,i1,…i_0,i_1,\dotsi0​,i1​,… drawn independently from ν\nuν, iterates

xk+1=xk−αk∇fik(xk).x_{k+1} = x_k - \alpha_k\nabla f_{i_k}(x_k).xk+1​=xk​−αk​∇fik​​(xk​).

The iterate xkx_kxk​ is a random vector, a function of i0,…,ik−1i_0,\dots,i_{k-1}i0​,…,ik−1​.

The standing assumptions on fff (§2 of the paper) are:

  1. ∇f\nabla f∇f is LLL-Lipschitz: ∥∇f(x)−∇f(y)∥≤L∥x−y∥\|\nabla f(x)-\nabla f(y)\|\le L\|x-y\|∥∇f(x)−∇f(y)∥≤L∥x−y∥ for all x,yx,yx,y, with L>0L>0L>0;
  2. the solution set is non-empty: some x∗x^*x∗ minimizes fff, and f∗=f(x∗)f^* = f(x^*)f∗=f(x∗);
  3. the PL inequality (3) with constant μ>0\mu>0μ>0: 12∥∇f(x)∥2≥μ (f(x)−f∗)\tfrac12\|\nabla f(x)\|^2\ge\mu\,(f(x)-f^*)21​∥∇f(x)∥2≥μ(f(x)−f∗) for all xxx.

Theorem 4 adds a second-moment bound along the iterates: Ei∼ν[∥∇fi(xk)∥2]≤C2\mathbb E_{i\sim\nu}\big[\|\nabla f_i(x_k)\|^2\big]\le C^2Ei∼ν​[∥∇fi​(xk​)∥2]≤C2 for every iterate xkx_kxk​.

Formalization targets

Goal: Theorem 4 (p. 6)

Under the assumptions above,

  1. with the decreasing step αk=2k+12μ(k+1)2\alpha_k = \dfrac{2k+1}{2\mu(k+1)^2}αk​=2μ(k+1)22k+1​,
E[f(xk)−f∗]≤LC22kμ2(k≥1);\mathbb E[f(x_k)-f^*]\le\frac{LC^2}{2k\mu^2}\qquad(k\ge1);E[f(xk​)−f∗]≤2kμ2LC2​(k≥1);
  1. with a constant step αk=α\alpha_k=\alphaαk​=α, 0<α<12μ0<\alpha<\dfrac1{2\mu}0<α<2μ1​,
E[f(xk)−f∗]≤(1−2μα)k [f(x0)−f∗]+LC2α4μ(k≥0).\mathbb E[f(x_k)-f^*]\le(1-2\mu\alpha)^k\,[f(x_0)-f^*]+\frac{LC^2\alpha}{4\mu}\qquad(k\ge0).E[f(xk​)−f∗]≤(1−2μα)k[f(x0​)−f∗]+4μLC2α​(k≥0).

The second-moment bound is assumed separately for the run of each scheme, with the same constant CCC.

Milestones (proof of Theorem 4, p. 7)

  1. One-step bound. For any point xxx with Ei∥∇fi(x)∥2≤C2\mathbb E_i\|\nabla f_i(x)\|^2\le C^2Ei​∥∇fi​(x)∥2≤C2 and any step aaa,
Ei[f(x−a∇fi(x))]≤f(x)−a∥∇f(x)∥2+LC2a22.\mathbb E_i[f(x-a\nabla f_i(x))]\le f(x)-a\|\nabla f(x)\|^2+\frac{LC^2a^2}{2}.Ei​[f(x−a∇fi​(x))]≤f(x)−a∥∇f(x)∥2+2LC2a2​.
  1. Recursion (10). Along the run, whenever αk≥0\alpha_k\ge0αk​≥0,
E[f(xk+1)−f∗]≤(1−2αkμ) E[f(xk)−f∗]+LC2αk22.\mathbb E[f(x_{k+1})-f^*]\le(1-2\alpha_k\mu)\,\mathbb E[f(x_k)-f^*]+\frac{LC^2\alpha_k^2}{2}.E[f(xk+1​)−f∗]≤(1−2αk​μ)E[f(xk​)−f∗]+2LC2αk2​​.
  1. Decreasing step rate, the first half of the goal.
  2. Constant step rate, the second half of the goal.

Significance

The result. Theorem 4 places SG under PL on the same footing as SG under strong convexity: the same O(1/k)O(1/k)O(1/k) rate, with constants depending on LLL, μ\muμ and CCC only, and no convexity. Since PL holds for least squares with any design matrix and, over any compact set, for logistic regression (§2.3 of the paper), the theorem gives explicit SG rates for these models. The constant-step bound quantifies the practitioners' rule of halving the step when progress stalls: halving α\alphaα at least halves the noise floor LC2α/(4μ)LC^2\alpha/(4\mu)LC2α/(4μ).

Formalizing it. The result is proved, in a short argument, in the paper. To our knowledge no machine-checked proof exists: the existing Lean formalizations of SG rates on Prove2Me concern stationarity of non-convex objectives under a variance bound, a different noise model with a different conclusion. The remaining work is to formalize the paper's proof, which requires a usable treatment of expectations over an i.i.d. index sequence (conditioning the step from xkx_kxk​ on the past) — infrastructure that every later analysis of a stochastic first-order method can reuse.

Difficulty

The deterministic content is elementary: the descent lemma, unbiasedness, the PL inequality, a telescoping sum and a geometric series. The difficulty is probabilistic bookkeeping. Recursion (10) takes the expectation of a one-step inequality with respect to iki_kik​ alone, holding xkx_kxk​ fixed; this is legitimate because iki_kik​ is independent of i0,…,ik−1i_0,\dots,i_{k-1}i0​,…,ik−1​, of which xkx_kxk​ is a function. Making this precise requires measurability of the iterates in the sample path and a Fubini-type decomposition of the product law into the first kkk coordinates, the kkk-th coordinate and the rest. A second subtlety is that the second-moment bound holds only at the iterates, not at every point, so the one-step bound must be applied pointwise in the path, on a set of full measure, rather than once and for all.

Formalization scope

The space is EuclideanSpace ℝ (Fin d) and ∇\nabla∇ is Mathlib's gradient. The index law is an arbitrary probability measure ν\nuν on a measurable space, and the path law is the infinite product Measure.infinitePi (fun _ => ν); the run is a deterministic function sgRun of the path. The Lipschitz-gradient assumption is stated literally on ∇f\nabla f∇f; the PL inequality is the published predicate ShiOptRates.PolyakLojasiewicz f (gradient f) μ (f xstar) (inequality (3) multiplied by 2); the solution set is a point xstar with ∀ y, f xstar ≤ f y, and f∗f^*f∗ is f xstar, never a free real.

Every expectation of a gap is a lower Lebesgue integral of max⁡(f(xk)−f∗,0)\max(f(x_k)-f^*,0)max(f(xk​)−f∗,0) with values in [0,∞][0,\infty][0,∞], which equals E[f(xk)−f∗]\mathbb E[f(x_k)-f^*]E[f(xk​)−f∗] because f∗f^*f∗ is the minimum; no integrability side condition can make a bound hold by default. The second-moment bound is the lower Lebesgue integral of ∥∇fi(xk)∥2\|\nabla f_i(x_k)\|^2∥∇fi​(xk​)∥2, required for almost every path and every kkk. Unbiasedness and problem (7) carry explicit integrability, so a non-integrable gradient cannot turn them into trivialities.

The second-moment bound is deliberately not imposed at every point xxx. With the bound at every point, ∥∇f∥≤C\|\nabla f\|\le C∥∇f∥≤C everywhere, so f−f∗≤C2/(2μ)f-f^*\le C^2/(2\mu)f−f∗≤C2/(2μ) by PL, while along the gradient-ascent flow f−f∗f-f^*f−f∗ grows exponentially unless it vanishes; hence fff would be constant and the theorem vacuous. Along the iterates the hypothesis is met by non-constant objectives such as f(x)=12∥x∥2f(x)=\tfrac12\|x\|^2f(x)=21​∥x∥2.

Additions relative to the page, all disclosed in the statements: each fif_ifi​ is differentiable; (i,x)↦∇fi(x)(i,x)\mapsto\nabla f_i(x)(i,x)↦∇fi​(x) is jointly measurable, so that the iterates are random variables; the constant step is positive; the decreasing-step rate is stated for k≥1k\ge1k≥1; and "each fff has an LLL-Lipschitz continuous gradient" is read as a condition on fff, the only one the proof uses. The page's proof of the constant-step bound writes the exponent kkk where k+1k+1k+1 is meant; the statements follow the theorem's (correct) form.

Contributions welcome: the measurability of the run, the conditioning step for infinitePi, the deterministic one-step bound, and the two scalar recursions, each of which is reusable for other SG analyses.

Selected references

  • H. Karimi, J. Nutini, M. Schmidt, Linear Convergence of Gradient and Proximal-Gradient Methods Under the Polyak-Łojasiewicz Condition, ECML PKDD 2016; arXiv:1608.04636v4, 2020. https://arxiv.org/abs/1608.04636
  • B. T. Polyak, Gradient methods for minimizing functionals, USSR Computational Mathematics and Mathematical Physics 3(4), 1963. https://doi.org/10.1016/0041-5553(63)90382-3
  • A. Nemirovski, A. Juditsky, G. Lan, A. Shapiro, Robust stochastic approximation approach to stochastic programming, SIAM Journal on Optimization 19(4), 2009. https://doi.org/10.1137/070704277
  • F. Bach, E. Moulines, Non-strongly-convex smooth stochastic approximation with convergence rate O(1/n), NeurIPS 2013. https://arxiv.org/abs/1306.2119
7 thms3 active usersReviewed
Bandit AlgorithmsOperations ResearchStatistics·Captain: mikedeng1

Explore First, Exploit Next: The True Shape of Regret in Bandit Problems V: Symmetric Monotonic Strategies Draw the Suboptimal Arms Collectively a Linear Number of TimesResearch Paper

Motivation

In a stochastic multi-armed bandit, a player repeatedly chooses one of KKK arms and receives a random reward drawn from the chosen arm's unknown distribution. The classical lower bounds of Lai and Robbins (1985) and Burnetas and Katehakis (1996) describe the regret of good strategies when the horizon TTT tends to infinity: it grows like ln⁡T\ln TlnT, with a constant determined by Kullback–Leibler divergences. They say nothing about small horizons.

Garivier, Ménard and Stoltz (Math. Oper. Res., 2019; arXiv:1602.07182v3) show that the regret has a definite shape over the whole range of TTT: it is almost linear at first ("explore first"), and only afterwards logarithmic ("exploit next"). Their §3 gives three non-asymptotic lower bounds for the initial phase, each for a natural class of strategies. This mission is the third of them, Theorem 4 (p. 13): a collective lower bound on the total number of pulls of all suboptimal arms. It is the bound that remains informative when many arms are suboptimal, and it is the one illustrated for large KKK in the paper's §3.4.

Setting

A bandit problem ν‾=(νa)a=1,…,K\underline\nu=(\nu_a)_{a=1,\dots,K}ν​=(νa​)a=1,…,K​ is a family of probability distributions on R\mathbb RR, each with an expectation μa=E(νa)\mu_a=E(\nu_a)μa​=E(νa​). A model D\mathcal DD is a set of such distributions, and ν‾\underline\nuν​ is in D\mathcal DD when every νa∈D\nu_a\in\mathcal Dνa​∈D. Put μ⋆=max⁡aμa\mu^\star=\max_a\mu_aμ⋆=maxa​μa​.

At each round t≥1t\ge1t≥1 a strategy ψ\psiψ chooses an arm AtA_tAt​ from the past rewards (possibly with auxiliary randomization) and observes a reward Yt∼νAtY_t\sim\nu_{A_t}Yt​∼νAt​​. The number of pulls of arm aaa up to round TTT is Nψ,a(T)=∑t=1T1{At=a}N_{\psi,a}(T)=\sum_{t=1}^T\mathbb 1\{A_t=a\}Nψ,a​(T)=∑t=1T​1{At​=a}, and Eν‾\mathbb E_{\underline\nu}Eν​​ is the expectation when ψ\psiψ interacts with ν‾\underline\nuν​.

The set of optimal arms is A⋆(ν‾)={a:μa=μ⋆}\mathcal A^\star(\underline\nu)=\{a:\mu_a=\mu^\star\}A⋆(ν​)={a:μa​=μ⋆}, of cardinality Aν‾⋆A^\star_{\underline\nu}Aν​⋆​, and the set of worst arms is W(ν‾)\mathcal W(\underline\nu)W(ν​), the arms of smallest mean. The paper defines a partial order: ν‾′≼ν‾\underline\nu'\preccurlyeq\underline\nuν​′≼ν​ when νa′=νa\nu'_a=\nu_aνa′​=νa​ for every optimal arm aaa of ν‾\underline\nuν​ and E(νa′)≤E(νa)E(\nu'_a)\le E(\nu_a)E(νa′​)≤E(νa​) for every other arm. Two classes of strategies enter the theorem.

  1. ψ\psiψ is pairwise symmetric for optimal arms on D\mathcal DD (Definition 3, p. 11) if, for every ν‾\underline\nuν​ in D\mathcal DD and every pair of optimal arms a⋆,a⋆a^\star,a_\stara⋆,a⋆​ with νa⋆=νa⋆\nu_{a^\star}=\nu_{a_\star}νa⋆​=νa⋆​​, for all T≥1T\ge1T≥1, the pairs (Nψ,a⋆(T),Nψ,a⋆(T))(N_{\psi,a^\star}(T),N_{\psi,a_\star}(T))(Nψ,a⋆​(T),Nψ,a⋆​​(T)) and (Nψ,a⋆(T),Nψ,a⋆(T))(N_{\psi,a_\star}(T),N_{\psi,a^\star}(T))(Nψ,a⋆​​(T),Nψ,a⋆​(T)) have the same distribution.
  2. ψ\psiψ is monotonic on D\mathcal DD (Definition 4, p. 13) if for all ν‾′≼ν‾\underline\nu'\preccurlyeq\underline\nuν​′≼ν​ in D\mathcal DD, ∑a⋆∈A⋆(ν‾′)Eν‾′[Nψ,a⋆(T)]≥∑a⋆∈A⋆(ν‾)Eν‾[Nψ,a⋆(T)]\sum_{a^\star\in\mathcal A^\star(\underline\nu')}\mathbb E_{\underline\nu'}[N_{\psi,a^\star}(T)]\ge\sum_{a^\star\in\mathcal A^\star(\underline\nu)}\mathbb E_{\underline\nu}[N_{\psi,a^\star}(T)]∑a⋆∈A⋆(ν​′)​Eν​′​[Nψ,a⋆​(T)]≥∑a⋆∈A⋆(ν​)​Eν​​[Nψ,a⋆​(T)]: on an easier problem the optimal arms are pulled at least as often.

With KL\mathrm{KL}KL the Kullback–Leibler divergence, the constant of the theorem is

Kν‾max⁡=min⁡w∈W(ν‾) max⁡a⋆∈A⋆(ν‾)KL(νw,νa⋆).\mathcal K^{\max}_{\underline\nu}=\min_{w\in\mathcal W(\underline\nu)}\ \max_{a^\star\in\mathcal A^\star(\underline\nu)}\mathrm{KL}(\nu_w,\nu_{a^\star}).Kν​max​=w∈W(ν​)min​ a⋆∈A⋆(ν​)max​KL(νw​,νa⋆​).

The Bernoulli divergence is kl(p,q)=pln⁡pq+(1−p)ln⁡1−p1−q\mathrm{kl}(p,q)=p\ln\frac pq+(1-p)\ln\frac{1-p}{1-q}kl(p,q)=plnqp​+(1−p)ln1−q1−p​, (5).

Formalization targets

Goal: Theorem 4 (main display), p. 13

For every model D\mathcal DD, every strategy ψ\psiψ pairwise symmetric for optimal arms and monotonic on D\mathcal DD, every bandit problem ν‾\underline\nuν​ in D\mathcal DD and every T≥1T\ge1T≥1,

∑a∉A⋆(ν‾)Eν‾[Nψ,a(T)] ≥ T(1−Aν‾⋆K−Aν‾⋆2TKν‾max⁡K−2Aν‾⋆TKν‾max⁡K).\sum_{a\notin\mathcal A^\star(\underline\nu)}\mathbb E_{\underline\nu}\big[N_{\psi,a}(T)\big]\ \ge\ T\left(1-\frac{A^\star_{\underline\nu}}K-\frac{A^\star_{\underline\nu}\sqrt{2T\mathcal K^{\max}_{\underline\nu}}}K-\frac{2A^\star_{\underline\nu}T\mathcal K^{\max}_{\underline\nu}}K\right).a∈/A⋆(ν​)∑​Eν​​[Nψ,a​(T)] ≥ T(1−KAν​⋆​​−KAν​⋆​2TKν​max​​​−K2Aν​⋆​TKν​max​​).

Milestones, in the order of the paper's proof

  1. Lemma 2 (p. 14): if (x−α)2≤βx(x-\alpha)^2\le\beta x(x−α)2≤βx with α,β≥0\alpha,\beta\ge0α,β≥0, then x≤α+β+αβx\le\alpha+\beta+\sqrt{\alpha\beta}x≤α+β+αβ​.
  2. Lemma 6 (p. 20): for 0≤p<q≤10\le p<q\le10≤p<q≤1, kl(p,q)≥(p−q)2/(2max⁡x∈[p,q]x(1−x))≥(p−q)2/(2q)\mathrm{kl}(p,q)\ge(p-q)^2/(2\max_{x\in[p,q]}x(1-x))\ge(p-q)^2/(2q)kl(p,q)≥(p−q)2/(2maxx∈[p,q]​x(1−x))≥(p−q)2/(2q).
  3. Monotonicity reduction (p. 13): with w~\tilde ww~ a worst arm and ν‾‾\underline{\underline\nu}ν​​ equal to ν‾\underline\nuν​ on the optimal arms and to νw~\nu_{\tilde w}νw~​ elsewhere, ∑a∉A⋆(ν‾)Eν‾[Nψ,a(T)]≥∑a∉A⋆(ν‾‾)Eν‾‾[Nψ,a(T)]\sum_{a\notin\mathcal A^\star(\underline\nu)}\mathbb E_{\underline\nu}[N_{\psi,a}(T)]\ge\sum_{a\notin\mathcal A^\star(\underline{\underline\nu})}\mathbb E_{\underline{\underline\nu}}[N_{\psi,a}(T)]∑a∈/A⋆(ν​)​Eν​​[Nψ,a​(T)]≥∑a∈/A⋆(ν​​)​Eν​​​[Nψ,a​(T)].
  4. Uniform split (p. 13): if all arms of ν~‾\underline{\tilde\nu}ν~​ carry the same distribution, Eν~‾[Nψ,a(T)]=T/K\mathbb E_{\underline{\tilde\nu}}[N_{\psi,a}(T)]=T/KEν~​​[Nψ,a​(T)]=T/K for every arm.
  5. The kl display (p. 14): TAν‾⋆Kν‾max⁡/K≥kl(Aν‾⋆/K,x)T A^\star_{\underline\nu}\mathcal K^{\max}_{\underline\nu}/K\ge\mathrm{kl}(A^\star_{\underline\nu}/K,x)TAν​⋆​Kν​max​/K≥kl(Aν​⋆​/K,x) with x=1T∑a⋆∈A⋆(ν‾)Eν‾‾[Nψ,a⋆(T)]x=\frac1T\sum_{a^\star\in\mathcal A^\star(\underline\nu)}\mathbb E_{\underline{\underline\nu}}[N_{\psi,a^\star}(T)]x=T1​∑a⋆∈A⋆(ν​)​Eν​​​[Nψ,a⋆​(T)].
  6. Equation (14) (p. 14): x≤Aν‾⋆K(1+2TKν‾max⁡+2TKν‾max⁡)x\le\frac{A^\star_{\underline\nu}}K\big(1+2T\mathcal K^{\max}_{\underline\nu}+\sqrt{2T\mathcal K^{\max}_{\underline\nu}}\big)x≤KAν​⋆​​(1+2TKν​max​+2TKν​max​​).

Significance

Theorem 4 bounds the regret from below through Rψ,ν‾,T≥(min⁡a∉A⋆Δa)∑a∉A⋆Eν‾[Nψ,a(T)]R_{\psi,\underline\nu,T}\ge(\min_{a\notin\mathcal A^\star}\Delta_a)\sum_{a\notin\mathcal A^\star}\mathbb E_{\underline\nu}[N_{\psi,a}(T)]Rψ,ν​,T​≥(mina∈/A⋆​Δa​)∑a∈/A⋆​Eν​​[Nψ,a​(T)], (13). While TTT is small compared with K/(Aν‾⋆Kν‾max⁡)K/(A^\star_{\underline\nu}\mathcal K^{\max}_{\underline\nu})K/(Aν​⋆​Kν​max​), the right-hand side of the theorem is of order T(1−A⋆/K)T(1-A^\star/K)T(1−A⋆/K), so the regret is linear: a reasonable strategy cannot avoid spending a constant fraction of the first rounds on the suboptimal arms. Unlike the per-arm bounds of Theorems 2 and 3, the bound concerns all suboptimal arms at once and does not degrade as KKK grows; §3.4 of the paper shows numerically that it gives the stronger regret bound for many arms.

The theorem is proved in the paper; to our knowledge no part of it has been machine-checked. A formal proof would exercise the published canonical bandit model of the platform (StochasticBandit, BanditPolicy, banditMeasure) on a genuine lower-bound argument with two auxiliary bandit problems, together with a local Pinsker inequality for Bernoulli divergences in [0,+∞][0,+\infty][0,+∞] that is reusable in other bandit lower bounds.

Difficulty

The obvious attempt applies the fundamental inequality (6) to ν‾\underline\nuν​ and an alternative in which one suboptimal arm becomes optimal, as in Theorem 2. That gives a bound arm by arm, which deteriorates with KKK and does not see the number A⋆A^\starA⋆ of optimal arms; it also needs a hypothesis of the type "smarter than uniform", which is not assumed here. Nothing is known about how a general ψ\psiψ behaves on ν‾\underline\nuν​ itself, so the two hypotheses on ψ\psiψ must carry all the information, and both are weak: Definition 3 only constrains pairs of optimal arms with equal laws, and only through the joint distribution of their pull counts under the canonical bandit measure; Definition 4 only compares problems ordered by ≼\preccurlyeq≼. On the analytic side, the classical Pinsker inequality kl(p,q)≥2(p−q)2\mathrm{kl}(p,q)\ge2(p-q)^2kl(p,q)≥2(p−q)2 loses the factor A⋆/KA^\star/KA⋆/K and gives a weaker theorem, so a local refinement is required, and every divergence is valued in [0,+∞][0,+\infty][0,+∞].

Formalization scope

Arms are indexed by Fin K (the paper's arm aaa is index a−1a-1a−1), with K≥1K\ge1K≥1 assumed in the goal. A bandit problem is the published BanditAlgorithm.StochasticBandit K, a strategy is a published BanditPolicy K (a family of Markov kernels from histories to arms; every strategy of the paper with auxiliary uniform randomization induces one with the same law of arms and rewards), and Eν‾[Nψ,a(T)]\mathbb E_{\underline\nu}[N_{\psi,a}(T)]Eν​​[Nψ,a​(T)] is the integral of the pull count against banditMeasure. A model is a Set (Measure ℝ) whose members are probability measures with an integrable identity, the standing assumption of §1.1–§1.2. Definitions 3 and 4 leave TTT implicit and are read "for all T≥1T\ge1T≥1"; in Definition 3 equality of distributions is equality of the image measures of the pair of counts, not of their expectations. KL\mathrm{KL}KL, kl\mathrm{kl}kl and Kmax⁡\mathcal K^{\max}Kmax take values in [0,+∞][0,+\infty][0,+∞]; the goal assumes Kν‾max⁡<+∞\mathcal K^{\max}_{\underline\nu}<+\inftyKν​max​<+∞, since otherwise the printed bound is −∞-\infty−∞ and a real-valued conversion would silently replace +∞+\infty+∞ by 000. The Bernoulli divergence is defined as the divergence between two Bernoulli laws on R\mathbb RR, so it is +∞+\infty+∞ where (5) is, rather than a finite junk value. The auxiliary problems of the proof are not constructed; the milestones quantify over them with their defining properties.

The paper's "In particular" clause, ∑a∉A⋆E[Nψ,a(T)]≥T2(1−A⋆/K)\sum_{a\notin\mathcal A^\star}\mathbb E[N_{\psi,a}(T)]\ge\frac T2(1-A^\star/K)∑a∈/A⋆​E[Nψ,a​(T)]≥2T​(1−A⋆/K) for T≤K/(8A⋆Kmax⁡)T\le K/(8A^\star\mathcal K^{\max})T≤K/(8A⋆Kmax), is not stated: it follows from the main display only when A⋆/K+A⋆/K≤1/2A^\star/K+\sqrt{A^\star/K}\le1/2A⋆/K+A⋆/K​≤1/2 (at A⋆/K=1/2A^\star/K=1/2A⋆/K=1/2 the main bound is negative), and whether it holds in general is not known here.

A formalization that weakens the strategy classes (full permutation symmetry, or equality of expected counts in Definition 3), drops the requirement that the problems lie in the model, or uses a real Bernoulli divergence with finite values at q∈{0,1}q\in\{0,1\}q∈{0,1} does not count.

A complete development needs the fundamental inequality (6) on the canonical bandit model (the subject of mission I of this series; the chain rule for the divergence of bandit measures is already proved on the platform), the symmetry argument on the canonical measure, Lemma 6 and Lemma 2. Contributions of any of these, and of the identity ∑aNψ,a(T)=T\sum_a N_{\psi,a}(T)=T∑a​Nψ,a​(T)=T integrated against banditMeasure, are welcome; Lemma 6 and the symmetry lemma are reusable beyond this mission.

Selected references

  • A. Garivier, P. Ménard, G. Stoltz, Explore First, Exploit Next: The True Shape of Regret in Bandit Problems, Mathematics of Operations Research, 2019; arXiv:1602.07182v3. https://arxiv.org/abs/1602.07182, https://doi.org/10.1287/moor.2017.0928
  • T. L. Lai, H. Robbins, Asymptotically efficient adaptive allocation rules, Advances in Applied Mathematics, 1985. https://doi.org/10.1016/0196-8858(85)90002-8
  • A. N. Burnetas, M. N. Katehakis, Optimal adaptive policies for sequential allocation problems, Advances in Applied Mathematics, 1996. https://doi.org/10.1006/aama.1996.0007
  • T. Lattimore, Cs. Szepesvári, Bandit Algorithms, Cambridge University Press, 2020. https://doi.org/10.1017/9781108571401
13 thms3 active usersReviewed
🏆Completed
Bandit AlgorithmsReinforcement Learning·Captain: mikedeng1

The Statistical Complexity of Interactive Decision Making VII: Policy-Cover Inverse Gap Weighting Certifies dec_γ ≤ 95H³SA/γ for Tabular Reinforcement LearningResearch Paper

Why an explicit exploration policy matters

In episodic reinforcement learning, an agent learns the rewards and dynamics of a system while choosing actions that determine which states it observes. A bound on the number of states and actions alone does not prescribe how to explore: a policy can be valuable because it earns reward, because it visits an uncertain state-action pair, or both. Foster, Kakade, Qian, and Rakhlin introduced the decision-estimation coefficient (DEC) to measure this trade-off in interactive decision problems. Their tabular reinforcement learning example gives an explicit policy distribution with a bound depending on the horizon and the sizes of the state and action sets Foster et al., §5.2.

The target is Proposition 5.6 of that paper. Its policy-cover inverse gap weighting strategy, or PC-IGW, places weight on policies that cover state-action pairs under an estimated model, while accounting for their estimated loss in reward. The paper also gives a posterior-sampling bound for a dual DEC in Proposition 5.4. PC-IGW addresses the frequentist coefficient with an explicit distribution and the constant 959595 Foster et al., Propositions 5.4 and 5.6.

Tabular models and the policy cover

A tabular Markov decision process has a finite state set SSS, a finite action set AAA, and H≥1H\ge1H≥1 decision layers. At layer hhh, state sss, and action aaa, a model MMM supplies a transition distribution PhM(⋅∣s,a)P_h^M(\cdot\mid s,a)PhM​(⋅∣s,a) and a reward distribution RhM(⋅∣s,a)R_h^M(\cdot\mid s,a)RhM​(⋅∣s,a). Every model shares the initial state distribution d1d_1d1​. A randomized nonstationary policy π\piπ supplies an action distribution πh(⋅∣s)\pi_h(\cdot\mid s)πh​(⋅∣s) for each layer and state; all such policies form ΠRNS\Pi_{\mathrm{RNS}}ΠRNS​ Foster et al., Example 1.2.

The trajectory law M(π)M(\pi)M(π) records the sequence of states, actions, and rewards generated by MMM and π\piπ. The final transition reaches a deterministic terminal state and contributes no model-dependent factor to this law. The policy's value is fM(π)=EM,π[∑h=1Hrh]f^M(\pi)=\mathbb E^{M,\pi}[\sum_{h=1}^H r_h]fM(π)=EM,π[∑h=1H​rh​]. Write πM\pi_MπM​ for any value-maximizing policy and dhM,π(s,a)d_h^{M,\pi}(s,a)dhM,π​(s,a) for the probability of visiting (s,a)(s,a)(s,a) at layer hhh. The squared Hellinger distance between two trajectory laws has the normalization DH2(P,Q)=∑x(P(x)−Q(x))2D_H^2(P,Q)=\sum_x(\sqrt{P(x)}-\sqrt{Q(x)})^2DH2​(P,Q)=∑x​(P(x)​−Q(x)​)2, with no factor 1/21/21/2 Foster et al., §2, Eq. (5).

Given an estimated model Mˉ\bar MMˉ and η>0\eta>0η>0, Algorithm 4 chooses, for each (h,s,a)(h,s,a)(h,s,a), a policy πh,s,a\pi_{h,s,a}πh,s,a​ that maximizes

dhMˉ,π(s,a)2HSA+η(fMˉ(πMˉ)−fMˉ(π))over π∈ΠRNS.\frac{d_h^{\bar M,\pi}(s,a)}{2HSA+\eta\bigl(f^{\bar M}(\pi_{\bar M})-f^{\bar M}(\pi)\bigr)} \quad\text{over }\pi\in\Pi_{\mathrm{RNS}}.2HSA+η(fMˉ(πMˉ​)−fMˉ(π))dhMˉ,π​(s,a)​over π∈ΠRNS​.

Let Ψ\PsiΨ be the resulting set of distinct cover policies. The PC-IGW distribution gives each π∈Ψ∪{πMˉ}\pi\in\Psi\cup\{\pi_{\bar M}\}π∈Ψ∪{πMˉ​} weight p(π)=1/[λ+η(fMˉ(πMˉ)−fMˉ(π))]p(\pi)=1/[\lambda+\eta(f^{\bar M}(\pi_{\bar M})-f^{\bar M}(\pi))]p(π)=1/[λ+η(fMˉ(πMˉ​)−fMˉ(π))], where λ\lambdaλ makes the weights sum to one. Proposition 5.7 states that this positive normalizer is unique and lies in [1,2HSA][1,2HSA][1,2HSA] Foster et al., Algorithm 4 and Proposition 5.7.

Formalization targets

The main target is the explicit-policy inequality of Proposition 5.6. For γ>0\gamma>0γ>0 and η=γ/(21H2)\eta=\gamma/(21H^2)η=γ/(21H2), the distribution from Algorithm 4 satisfies, for every normalized tabular model MMM,

Eπ∼p ⁣[fM(πM)−fM(π)−γDH2(M(π),Mˉ(π))]≤95H3SAγ.\mathbb E_{\pi\sim p}\!\left[ f^M(\pi_M)-f^M(\pi)-\gamma D_H^2(M(\pi),\bar M(\pi)) \right]\le\frac{95H^3SA}{\gamma}.Eπ∼p​[fM(πM​)−fM(π)−γDH2​(M(π),Mˉ(π))]≤γ95H3SA​.

The uniform quantifier over MMM is the paper's supremum over its model class. Because this is a bound for a particular distribution ppp, it certifies the stated DEC bound. The Lean goal expresses that explicit distribution and its inequality directly; it does not introduce a separate supremum over all models and infimum over all distributions. The milestone list records the unique normalizer, the exploration-bias display (59), the cover comparison (61), the local and global simulation estimates, the Hellinger information inequalities, and Lemma 5.2's change of measure Foster et al., pp. 36–38 and Appendices A.2, F.1.

What the bound contributes

The result controls the cost of exploration for a finite-horizon model class by H3SA/γH^3SA/\gammaH3SA/γ with a specified constant. Its importance is the constructive policy distribution: the bound applies to the strategy returned by Algorithm 4, rather than asserting only that some distribution exists. In the paper, this bound combines with estimation and the general DEC-to-regret result to obtain a tabular regret guarantee; the authors also use the policy-cover idea as a starting point for their bilinear-class treatment Foster et al., §§5.2.2–5.2.3 and 7.

The paper result is proved. This mission asks for machine-checked proofs of its finite-alphabet statements; the draft theorems currently carry sorry. Its MDP interface includes reward distributions, full trajectory laws, occupancy probabilities, backward values, and local divergences. Those definitions can be reused for other tabular DEC results. Published FoundationsRL definitions already provide the discrete Hellinger and total-variation distances and a proved comparison between them. A separate published MDP interface records reward means, which cannot represent the reward-kernel Hellinger terms in this result Foster and Rakhlin, Foundations of Reinforcement Learning and Interactive Decision Making.

Where the difficulty lies

The cover is chosen using occupancies under Mˉ\bar MMˉ, while the regret is evaluated under any MMM in the class. An ordinary value bound for one fixed policy does not by itself control this mismatch. The page's intermediate expressions also divide by an average occupancy that may be zero. A formal statement must still cover state-action pairs that the policy distribution never visits, without assigning meaning to a zero-denominator ratio. The final inequality uses Hellinger distance between whole trajectories even though local simulation estimates are expressed through transition and reward kernels. Keeping those two levels of comparison aligned is the main mathematical burden Foster et al., pp. 37–38.

There is also a reward-bound issue in the appendix. Lemma F.3 uses continuation values in [0,1][0,1][0,1] at states visited under Mˉ\bar MMˉ, including states that may be unreachable under MMM from the shared initial law. The paper states an almost-sure bound on total episode reward, but that statement alone does not control every such continuation value Foster et al., Lemma F.3, p. 107.

Formalization scope

State, action, and reward alphabets are finite and nonempty. Probability laws and expectations are finite vectors and sums. Policies may randomize; the cover condition compares its choice against every valid randomized nonstationary policy. Layers use zero-based Fin H indices, while the paper uses 1,…,H1,\ldots,H1,…,H. All models share d1d_1d1​. The terminal transition is omitted from the trajectory law and from one-step divergence sums. Duplicate cover policies are merged, as required by the paper's set Ψ∪{πMˉ}\Psi\cup\{\pi_{\bar M}\}Ψ∪{πMˉ​}. The selector πM\pi_MπM​ must maximize value over valid policies, and the algorithm's weights must sum to one; these conditions exclude a free policy choice or an unnormalized pseudo-law.

The paper assumes ∑hrh∈[0,1]\sum_h r_h\in[0,1]∑h​rh​∈[0,1] almost surely. The local simulation lemma's stated bound needs two disclosed additions here: nonnegative stage rewards, and the same total-reward bound on every supported suffix started from any state. RewardsNormalized records them. The change-of-measure and global-simulation milestones use the weaker episode-only condition where sufficient. The cover comparison is cross-multiplied so it remains meaningful when an occupancy is zero. The goal's right side is a real number with γ>0\gamma>0γ>0, H≥1H\ge1H≥1, and nonempty S,AS,AS,A, so its divisions and size factors have the intended meaning.

A complete development needs the finite trajectory law, the policy-cover normalizer, occupancy bounds, the two simulation lemmas, Hellinger comparison inequalities, and the change-of-measure result. Contributions to these foundations, to the exact constant in Proposition 5.6, and to later generalizations of the cover are within scope. The abstract DEC infimum over all randomized policies is a later interface target; this mission's goal is the stronger, explicit-strategy certificate printed first in the proposition.

Selected references

  • Dylan J. Foster, Sham M. Kakade, Jian Qian, and Alexander Rakhlin, The Statistical Complexity of Interactive Decision Making, arXiv preprint, version 3, 2023. arXiv:2112.13487v3.
  • Dylan J. Foster and Alexander Rakhlin, Foundations of Reinforcement Learning and Interactive Decision Making, arXiv preprint, 2023. arXiv:2312.16730.
16 thms3 active usersReviewed
Bandit AlgorithmsOperations ResearchStatistics·Captain: mikedeng1

Explore First, Exploit Next: The True Shape of Regret in Bandit Problems III: A Strategy Smarter Than Uniform Draws Every Arm at Least (T/K)(1 − √(2T K_inf)) TimesResearch Paper

Why early exploration needs a lower bound

A stochastic multi-armed bandit learner repeatedly selects one of finitely many arms and sees only the reward from the selected arm. The learner must gather information about uncertain rewards while trying to select high-mean arms. Expected regret summarizes the lost reward relative to always selecting an optimal arm, but regret alone does not explain how much sampling each particular arm must receive. Garivier, Ménard, and Stoltz study this allocation question directly in Explore First, Exploit Next: The True Shape of Regret in Bandit Problems. Their Section 3 gives lower bounds for small horizons, before the longer-run logarithmic behavior becomes dominant.

The paper's first small-horizon result concerns strategies that allocate at least a uniform share to every optimal arm, on every problem in a specified model. This requirement excludes a strategy that always selects one fixed label and succeeds only when that label happens to be best. Such a requirement still permits many adaptive policies; it constrains their performance on optimal arms rather than prescribing their decisions round by round. The result establishes a compulsory initial allocation to every arm, including one whose mean is below the best mean. It is Theorem 2 on page 10 of the pinned arXiv:1602.07182v3.

Bandit problems, models, and strategies

A bandit problem ν‾=(νa)a=1K\underline\nu=(\nu_a)_{a=1}^Kν​=(νa​)a=1K​ consists of KKK probability laws on R\mathbb RR. At each round, a strategy ψ\psiψ chooses an arm using the observed past and randomization, then receives an independent reward drawn from that arm's law. The mean reward of arm aaa is μa=Eνa[X]\mu_a=\mathbb E_{\nu_a}[X]μa​=Eνa​​[X]. Write μ∗=max⁡aμa\mu^*=\max_a\mu_aμ∗=maxa​μa​, and call an arm optimal when μa=μ∗\mu_a=\mu^*μa​=μ∗. The count Nψ,a(T)N_{\psi,a}(T)Nψ,a​(T) is the number of times arm aaa is chosen during the first TTT rounds. These objects and the information available to a strategy are specified in Section 1.1, pages 2–3.

A model D\mathcal DD is a collection of possible arm laws, each a probability measure with an expectation. The notation ν‾∈D\underline\nu\in\mathcal Dν​∈D means that every arm law νa\nu_aνa​ belongs to D\mathcal DD. A strategy is smarter than uniform on D\mathcal DD when, for every bandit problem in the model, every optimal arm a∗a^*a∗, and every integer T≥1T\ge1T≥1, it satisfies

Eν‾[Nψ,a∗(T)]≥TK.\mathbb E_{\underline\nu}[N_{\psi,a^*}(T)]\ge\frac{T}{K}.Eν​​[Nψ,a∗​(T)]≥KT​.

This is Definition 2 on page 10. It imposes a lower bound on each optimal arm, even when several arms tie. The requirement is evaluated across the whole model, so a strategy cannot be declared smarter based only on its behavior on the one problem used in the theorem.

For a law PPP, a real threshold xxx, and the model D\mathcal DD, the information number Kinf⁡(P,x,D)\mathcal K_{\inf}(P,x,\mathcal D)Kinf​(P,x,D) is the infimum of KL(P,Q)\mathrm{KL}(P,Q)KL(P,Q) over those Q∈DQ\in\mathcal DQ∈D whose mean exceeds xxx. The infimum is +∞+\infty+∞ if there is no such QQQ. This is the definition on page 4. It measures how close an arm law is, in Kullback–Leibler divergence, to another admissible law that would give it a mean above the present optimum.

Formalization targets

The main target is Theorem 2. For a model D\mathcal DD, a strategy ψ\psiψ smarter than uniform on it, a bandit problem ν‾∈D\underline\nu\in\mathcal Dν​∈D, every arm aaa, and every integer T≥1T\ge1T≥1, the finite-Kinf⁡\mathcal K_{\inf}Kinf​ case states

Eν‾[Nψ,a(T)]≥TK(1−2TKinf⁡(νa,μ∗,D)).\mathbb E_{\underline\nu}[N_{\psi,a}(T)]\ge \frac{T}{K}\left(1-\sqrt{2T\mathcal K_{\inf}(\nu_a,\mu^*,\mathcal D)}\right).Eν​​[Nψ,a​(T)]≥KT​(1−2TKinf​(νa​,μ∗,D)​).

The theorem includes optimal and suboptimal arms. Its printed companion says that when T≤1/(8Kinf⁡)T\le 1/(8\mathcal K_{\inf})T≤1/(8Kinf​), the arm receives at least T/(2K)T/(2K)T/(2K) expected selections. In the formal statement this threshold is 8TKinf⁡≤18T\mathcal K_{\inf}\le18TKinf​≤1, which also expresses the intended case Kinf⁡=0\mathcal K_{\inf}=0Kinf​=0. Both claims come from Theorem 2, page 10.

The milestone list follows the source's intermediate claims. Lemma 6, page 20 bounds Bernoulli relative entropy using the maximum of x(1−x)x(1-x)x(1−x) over an interval. The proof of Theorem 2 on page 11 states that Bernoulli relative entropy increases in its second argument on the interval starting at its first argument. The same page applies the paper's fundamental inequality (6) to the fraction of rounds spent on one arm, then states a lower bound against a single alternative arm law. These are the four curated milestones for the main theorem.

What the result supplies

The theorem turns a model-relative performance condition into a numerical constraint on early sampling. If an arm's law is close in divergence to a law that would make it best, then the bound keeps its expected pull count near the uniform allocation until the horizon is large enough for the term under the square root to matter. It identifies the divergence scale 1/Kinf⁡1/\mathcal K_{\inf}1/Kinf​ at which this particular guarantee becomes weak. The paper contrasts this absolute bound with later relative and collective bounds that capture different dependencies on the number of arms in Section 3, pages 10–14.

The mathematical result is proved in the paper. The remaining work in this mission is a machine-checked development of the statement and its source-backed intermediate inequalities in the canonical stochastic-bandit model. The published bandit structure, policy kernels, and information number provide reusable foundations. Shared definitions supply the Bernoulli divergence, model predicate, and expected pull counts; the local strategy class makes the hypothesis of this particular theorem explicit. The Bernoulli inequalities and the pull-count change-of-measure statement can be reused in other lower-bound results.

Where the difficulty lies

A performance condition about optimal arms does not directly constrain an arm that is suboptimal in the problem of interest. The small-horizon bound must remain valid even though the strategy may adapt to observed rewards, and the relevant alternative arm laws can have divergences approaching an infimum without attaining it. The paper's lower bound is therefore sensitive to the exact orientation of the two bandit problems, the expected count multiplying their divergence, and the boundary conventions for Bernoulli relative entropy. Replacing the infimum with a chosen alternative or changing a divergence to an ordinary real with an arbitrary boundary value would alter the statement.

Formalization scope

Lean uses arms Fin K, indexed from zero; the paper indexes them from one. The horizon is a natural number, and counts are integrated against the probability measure on TTT-round histories generated by a stochastic kernel policy. Such kernels represent the arm-and-reward histories induced by the paper's strategies with auxiliary randomization. The model is a set of probability measures on R\mathbb RR whose identity function is integrable. The published StochasticBandit, BanditPolicy, and ConsistentBanditPolicy definitions supply the basic problem, history measure, pull count, means, and Kinf⁡\mathcal K_{\inf}Kinf​.

Bernoulli divergence and Kinf⁡\mathcal K_{\inf}Kinf​ take values in [0,+∞][0,+\infty][0,+∞]. The real square-root formula assumes Kinf⁡<+∞\mathcal K_{\inf}<+\inftyKinf​<+∞: when it is infinite, the paper's right-hand side is −∞-\infty−∞, while converting infinity directly to Lean's real value would produce a false numerical bound. The starting bandit problem is explicitly in D\mathcal DD, a condition needed for the proof's use of Definition 2 on a one-arm modification, although the printed quantifier on page 10 omits it. The arm parameter guarantees K≥1K\ge1K≥1; every target assumes T≥1T\ge1T≥1, so its real divisions by KKK and TTT have the stated meaning.

A faithful development must retain Definition 2's quantification over all model problems and each optimal arm, rather than only over the bandit problem used in the conclusion. It must use the full Bernoulli relative entropy, including infinite endpoint values, and the infimum over precisely those model laws with mean greater than μ∗\mu^*μ∗. The local instance of (6) also assumes integrable arm rewards in both bandit problems, as Section 1.1 requires. Contributions that establish the four milestones, the companion threshold claim, or reusable count-integrability and divergence facts are in scope.

Selected references

  • Aurélien Garivier, Pierre Ménard, and Gilles Stoltz, Explore First, Exploit Next: The True Shape of Regret in Bandit Problems, Mathematics of Operations Research, 2019; arXiv:1602.07182v3.
  • Tor Lattimore and Csaba Szepesvári, Bandit Algorithms, Cambridge University Press, 2020; DOI:10.1017/9781108571401. The published bandit definitions used here follow its canonical framework.
11 thms3 active usersReviewed
Bandit AlgorithmsOperations ResearchStatistics·Captain: mikedeng1

Explore First, Exploit Next: The True Shape of Regret in Bandit Problems VI: Uniformly Super-Fast Convergent Strategies Draw Each Suboptimal Arm ln T / K_inf Times up to Explicit Lower-Order TermsResearch Paper

Why the large-horizon bound matters

In a stochastic multi-armed bandit, a decision maker repeatedly chooses one of several unknown reward sources. Sampling an arm reveals its reward but uses a round that could have been spent on a better arm. A good strategy must therefore reduce sampling of inferior arms, yet it cannot stop sampling them too early: it must still distinguish the current problem from a nearby one in which an apparently inferior arm is best. This tension produces the logarithmic lower bounds that benchmark bandit algorithms.

Garivier, Ménard and Stoltz revisit that benchmark in Explore First, Exploit Next. Their Theorem 5 gives a finite-horizon lower bound with explicit terms beyond the leading logarithm. It concerns every uniformly super-fast convergent strategy on a well-behaved model, not a particular algorithm. The distinction matters when comparing an algorithm's upper bound to what any strategy with that uniform performance guarantee must pay. The theorem is proved in the paper; this mission asks for a machine-checked statement and proof of that established result.

Bandit model and information cost

There are KKK arms. Arm aaa has an unknown probability distribution νa\nu_aνa​ over real rewards, with finite mean μa\mu_aμa​. The vector ν‾=(νa)a=1K\underline\nu=(\nu_a)_{a=1}^Kν​=(νa​)a=1K​ is a bandit problem. Its best mean is μ∗=max⁡aμa\mu^*=\max_a\mu_aμ∗=maxa​μa​, and the gap of arm aaa is Δa=μ∗−μa\Delta_a=\mu^*-\mu_aΔa​=μ∗−μa​. An arm with Δa>0\Delta_a>0Δa​>0 is suboptimal. A strategy ψ\psiψ selects an arm each round using the previous arms and rewards, possibly with randomization. Write Nψ,a(T)N_{\psi,a}(T)Nψ,a​(T) for the number of times it selects arm aaa by round TTT, and Eν‾[Nψ,a(T)]\mathbb E_{\underline\nu}[N_{\psi,a}(T)]Eν​​[Nψ,a​(T)] for the expected count under the problem ν‾\underline\nuν​.

A model D\mathcal DD specifies the probability distributions that an arm may have. All its laws have a finite first moment, and a bandit problem is in D\mathcal DD when every arm law belongs to it. For a law PPP and a target mean xxx, the paper's minimal information cost is

Kinf⁡(P,x,D)=inf⁡{KL(P,Q):Q∈D, E(Q)>x},\mathcal K_{\inf}(P,x,\mathcal D) =\inf\{\mathrm{KL}(P,Q):Q\in\mathcal D,\ \mathrm E(Q)>x\},Kinf​(P,x,D)=inf{KL(P,Q):Q∈D, E(Q)>x},

where an empty infimum is +∞+\infty+∞. This measures how much the law of an arm must change to make its mean exceed xxx. The Bernoulli divergence kl(p,q)\mathrm{kl}(p,q)kl(p,q) used in the intermediate change-of-measure bound keeps its infinite values at the boundary.

Let E(D)E(\mathcal D)E(D) be the interior of the set of expectations of laws in D\mathcal DD. A model is well behaved when positive functions εD\varepsilon_{\mathcal D}εD​ and ωD\omega_{\mathcal D}ωD​ satisfy, whenever P∈DP\in\mathcal DP∈D, x∈E(D)x\in E(\mathcal D)x∈E(D), E(P)<x\mathrm E(P)<xE(P)<x, and 0<ε<εD(x)0<\varepsilon<\varepsilon_{\mathcal D}(x)0<ε<εD​(x),

Kinf⁡(P,x+ε,D)≤Kinf⁡(P,x,D)+εωD(P,x).\mathcal K_{\inf}(P,x+\varepsilon,\mathcal D) \le \mathcal K_{\inf}(P,x,\mathcal D)+\varepsilon\omega_{\mathcal D}(P,x).Kinf​(P,x+ε,D)≤Kinf​(P,x,D)+εωD​(P,x).

A strategy is uniformly super-fast convergent with constant Cψ,DC_{\psi,\mathcal D}Cψ,D​ when Eν‾[Nψ,a(T)]≤Cψ,Dln⁡T/Δa2\mathbb E_{\underline\nu}[N_{\psi,a}(T)]\le C_{\psi,\mathcal D}\ln T/\Delta_a^2Eν​​[Nψ,a​(T)]≤Cψ,D​lnT/Δa2​ for every problem in D\mathcal DD, every suboptimal arm, and every T≥2T\ge2T≥2. The same constant applies throughout. Finally, H(ν‾)=∑a:Δa>0Δa−2H(\underline\nu)=\sum_{a:\Delta_a>0}\Delta_a^{-2}H(ν​)=∑a:Δa​>0​Δa−2​ records the inverse-squared gaps of all suboptimal arms.

Formalization targets

For each suboptimal arm let K=Kinf⁡(νa,μ∗,D)\mathcal K=\mathcal K_{\inf}(\nu_a,\mu^*,\mathcal D)K=Kinf​(νa​,μ∗,D). The goal is Theorem 5, with εT=(ln⁡T)−4\varepsilon_T=(\ln T)^{-4}εT​=(lnT)−4 and

aT=ωD(νa,μ∗)KεT,bT=Cψ,DH(ν‾)ln⁡TT,cT=ln⁡(KCψ,D(ln⁡T)9)ln⁡T.a_T=\frac{\omega_{\mathcal D}(\nu_a,\mu^*)}{\mathcal K}\varepsilon_T, \qquad b_T=C_{\psi,\mathcal D}H(\underline\nu)\frac{\ln T}{T}, \qquad c_T=\frac{\ln(KC_{\psi,\mathcal D}(\ln T)^9)}{\ln T}.aT​=KωD​(νa​,μ∗)​εT​,bT​=Cψ,D​H(ν​)TlnT​,cT​=lnTln(KCψ,D​(lnT)9)​.

For every T≥2T\ge2T≥2 satisfying εT<εD(μ∗)\varepsilon_T<\varepsilon_{\mathcal D}(\mu^*)εT​<εD​(μ∗), aT<1a_T<1aT​<1, bT<1b_T<1bT​<1, and 0≤cT<10\le c_T<10≤cT​<1, the target is

Eν‾[Nψ,a(T)]≥ln⁡TK(1−(aT+bT+cT))−ln⁡2K.\mathbb E_{\underline\nu}[N_{\psi,a}(T)] \ge \frac{\ln T}{\mathcal K}\bigl(1-(a_T+b_T+c_T)\bigr) -\frac{\ln 2}{\mathcal K}.Eν​​[Nψ,a​(T)]≥KlnT​(1−(aT​+bT​+cT​))−Kln2​.

The milestones follow the paper's displayed bounds: the one-alternative inequality (15), the optimized inequality (17), the inverse-information-cost bound after (17), and the product inequality quoted immediately before Theorem 5. They preserve the same constants and the same H(ν‾)H(\underline\nu)H(ν​).

What the result supplies

The leading term ln⁡T/K\ln T/\mathcal KlnT/K is the familiar distribution-dependent sampling barrier. The explicit aTa_TaT​, bTb_TbT​, and cTc_TcT​ quantify the correction at a finite horizon; the paper observes that their combined contribution has order ln⁡ln⁡T\ln\ln TlnlnT under its conditions. This makes the result useful when an algorithm's upper bound includes a second-order term and one wants to know the scale that a lower bound permits. It also identifies precisely which model regularity and uniform strategy guarantee support that comparison.

Formalizing Theorem 5 would add a reusable finite-horizon lower bound to the bandit library. The published definitions of stochastic bandits, policies, pull counts, and Kinf⁡\mathcal K_{\inf}Kinf​ supply the base model; this mission adds the paper's well-behavedness and super-fast convergence predicates and states the intermediate inequalities as separate targets. The paper's proof is established mathematically, while these mission statements remain open Lean theorems until solvers supply machine-checked proofs.

Main difficulty

The information cost is an infimum over alternative reward laws, so the finite-horizon bound must remain meaningful when the feasible set is empty or its divergence is infinite. The displayed formulas also divide by that information cost and use logarithms of quantities involving the strategy's uniform constant. For Theorem 5, the explicit correction factors must all lie in the range where the paper's final product inequality applies. A purely asymptotic statement would discard exactly the constants and horizon conditions that distinguish this theorem.

Formalization scope

Arms are Fin K, so the paper's arm a∈{1,…,K}a\in\{1,\ldots,K\}a∈{1,…,K} has Lean index a−1a-1a−1. A strategy is represented by Markov kernels on observed arm-reward histories; its expected pull count is an integral of a bounded count. The model consists of probability measures on R\mathbb RR with integrable identity. The information cost and Bernoulli divergence live in [0,+∞][0,+\infty][0,+∞]; a finite positive cost is converted to a real number only where a displayed reciprocal requires it. The goal contains a suboptimal arm, which rules out the empty-arm and single-arm cases. The model's witness functions are total in Lean, with requirements only on their paper-defined domains.

The theorem explicitly assumes μ∗∈E(D)\mu^*\in E(\mathcal D)μ∗∈E(D) and 0<K<+∞0<\mathcal K<+\infty0<K<+∞, which are implicit in the paper's use of εD(μ∗)\varepsilon_{\mathcal D}(\mu^*)εD​(μ∗) and 1/K1/\mathcal K1/K. It also assumes cT≥0c_T\ge0cT​≥0: the printed proof invokes a product inequality requiring that sign, although the theorem prints only cT<1c_T<1cT​<1. The version of (17) in the milestone list assumes bT≤1b_T\le1bT​≤1 in the case of a negative logarithm. These conditions and the model's integrability prevent vacuous or totalized-value readings of the lower bound. Contributions toward the four displayed milestones, the required measure-theoretic change of measure, and the finite-horizon arithmetic are in scope.

Selected references

  • Aurélien Garivier, Pierre Ménard, and Gilles Stoltz, Explore First, Exploit Next: The True Shape of Regret in Bandit Problems, Mathematics of Operations Research 44 (2019); arXiv:1602.07182v3, DOI:10.1287/moor.2017.0928.
  • Tor Lattimore and Csaba Szepesvári, Bandit Algorithms, Cambridge University Press (2020), Chapters 4 and 16; DOI:10.1017/9781108571401.
11 thms3 active usersReviewed
🏆Completed
Bandit AlgorithmsReinforcement Learning·Captain: mikedeng1

The Statistical Complexity of Interactive Decision Making VI: Posterior Sampling Bounds the Dual DEC of Tabular Reinforcement Learning by 26H²SA/γResearch Paper

Why tabular reinforcement learning needs a decision-estimation bound

In episodic reinforcement learning, an agent acts repeatedly in an unknown system and must learn which policy earns the largest expected reward. Even when states and actions are finite, the agent's choices change which data it observes. Foster, Kakade, Qian, and Rakhlin use the decision-estimation coefficient (DEC) to connect this information constraint to regret across interactive decision problems. Their tabular example asks for an explicit bound in terms of the horizon, the numbers of states and actions, and an information penalty Foster et al., §5.2.

The result here is the posterior-sampling bound of Proposition 5.4. It gives a uniform constant, rather than a rate stated only with asymptotic notation. Its role is to show that a policy drawn from a distribution of model-optimal policies controls the Bayesian or dual DEC in a tabular class. The same article also studies a frequentist policy-cover strategy with a different coefficient; this mission isolates the posterior-sampling argument and the information inequalities it needs Foster et al., Propositions 5.4 and 5.6.

Models, policies, and trajectory laws

A tabular Markov decision process has a finite state set SSS, a finite action set AAA, and H≥1H\ge1H≥1 decision layers. At layer hhh, state sss, and action aaa, the model MMM supplies a transition distribution PhM(⋅∣s,a)P_h^M(\cdot\mid s,a)PhM​(⋅∣s,a) and a reward distribution RhM(⋅∣s,a)R_h^M(\cdot\mid s,a)RhM​(⋅∣s,a). Every model has the same initial state distribution d1d_1d1​. A randomized nonstationary policy π\piπ supplies an action distribution πh(⋅∣s)\pi_h(\cdot\mid s)πh​(⋅∣s) at each layer and state. This is the policy class ΠRNS\Pi_{\mathrm{RNS}}ΠRNS​ used in Example 1.2, including genuinely randomized policies Foster et al., Example 1.2.

The trajectory law M(π)M(\pi)M(π) records states, actions, and realized rewards over one episode. Its weight on a trajectory is the product of the shared initial probability, the policy probabilities, reward-kernel probabilities, and transitions between consecutive decision layers. The state after layer HHH is a deterministic terminal state. Therefore the last model transition contributes no factor to the recorded law. The reward value of an outcome is given by a fixed map rrr, and the value of a policy is fM(π)=EM,π[∑h=1Hrh]f^M(\pi)=\mathbb E^{M,\pi}[\sum_{h=1}^{H}r_h]fM(π)=EM,π[∑h=1H​rh​]. The occupancy dhM,π(s,a)d_h^{M,\pi}(s,a)dhM,π​(s,a) is the probability of visiting state-action pair (s,a)(s,a)(s,a) at layer hhh Foster et al., §2.

For each model, πM\pi_MπM​ denotes a policy maximizing fMf^MfM over ΠRNS\Pi_{\mathrm{RNS}}ΠRNS​. Let μ\muμ be a prior on models. Posterior sampling here means drawing a model M′∼μM'\sim\muM′∼μ and using πM′\pi_{M'}πM′​. The model whose performance is evaluated is a separate draw M∼μM\sim\muM∼μ; these two draws must remain independent. The reference model is denoted Mˉ\bar MMˉ. Squared Hellinger distance has the paper's normalization DH2(P,Q)=∑x(P(x)−Q(x))2D_H^2(P,Q)=\sum_x(\sqrt{P(x)}-\sqrt{Q(x)})^2DH2​(P,Q)=∑x​(P(x)​−Q(x)​)2, with no factor of one half Foster et al., Eq. (5).

Formalization targets

The main target is the policy-specific inequality of Proposition 5.4. For γ>0\gamma>0γ>0, it asserts

EM∼μ, π∼p[fM(πM)−fM(π)−γDH2(M(π),Mˉ(π))]≤26H2SAγ,p=Law⁡(πM′), M′∼μ.\mathbb E_{M\sim\mu,\,\pi\sim p} \left[f^M(\pi_M)-f^M(\pi)-\gamma D_H^2(M(\pi),\bar M(\pi))\right] \le \frac{26H^2SA}{\gamma}, \qquad p=\operatorname{Law}(\pi_{M'}),\ M'\sim\mu.EM∼μ,π∼p​[fM(πM​)−fM(π)−γDH2​(M(π),Mˉ(π))]≤γ26H2SA​,p=Law(πM′​), M′∼μ.

The claim holds for each finitely supported prior μ\muμ in the finite-alphabet representation. It directly gives the same upper bound when the dual DEC is defined as the supremum over such priors and the infimum over policy distributions. The Lean goal states the inequality for posterior sampling itself; no real supremum is needed to express that stronger, explicit claim. The reference model belongs to the normalized tabular class, which is the case actually proved in §5.2.1. The printed proposition also names reference models outside the class and points to the later general result Foster et al., p. 33.

The milestones are the conditional-to-joint Hellinger inequality of Lemma A.9, both inequalities of Lemma A.11, the global and local simulation lemmas F.2 and F.3, the change-of-measure statement in Lemma 5.2, and the unnumbered posterior-sampling decoupling inequality on pp. 34–35. Together they identify the stated estimates on which the proposition rests Foster et al., Appendices A.2 and F.1.

What the result contributes

The 26H2SA/γ26H^2SA/\gamma26H2SA/γ bound quantifies the price of learning the dynamics and rewards while selecting policies. Its explicit dependence on SSS and AAA is linear, and its dependence on the horizon is quadratic. In the paper, this estimate feeds the DEC-based analysis of tabular reinforcement learning and provides a comparison point for the policy-cover bound and later extensions beyond tabular models Foster et al., §5.2.

This formalization supplies a reusable interface for tabular MDPs with reward distributions: a full trajectory law, policy value, occupancy, backward values, and local divergences. The published FoundationsRL divergence definitions and a proved TV–Hellinger comparison can be reused. None of the seven statements has a machine-checked proof yet. A finite-horizon MDP object already published elsewhere records deterministic reward means and does not express the reward-kernel Hellinger distances required here.

Where the difficulty lies

A value difference is controlled by total variation between complete trajectory laws, but the proposed bound is expressed through local transition and reward discrepancies and ultimately through squared Hellinger distance between whole trajectories. These quantities average under different model-policy laws. Moving between them while retaining the exact constants in Lemma 5.2 is the central obstacle. A direct bound for a policy paired with the same model does not supply the double expectation in the target: the performance model and the policy-generating model are independent under posterior sampling.

The reward convention adds a second obstacle. Section 5.2 sets the reward space to [0,1][0,1][0,1] and also assumes that total episode reward lies in [0,1][0,1][0,1] almost surely. A model's continuation values at states it cannot reach from d1d_1d1​ are not controlled by the episode condition. The local simulation bound is nevertheless asserted under the reference model's occupancy, which may visit such states.

Formalization scope

State, action, and reward alphabets are finite; reward distributions are finite probability vectors. The paper allows general measurable reward spaces, so the finite reward alphabet is a declared scope restriction. Layers are indexed by Fin H, starting at zero, while the source labels them 1,…,H1,\ldots,H1,…,H. The terminal transition is omitted from trajectory products and local discrepancy sums. The initial distribution is a separate parameter shared by all models. Priors and policy laws have finite support, as in the paper's discrete-topology convention for model distributions. Every policy quantified in a theorem is checked to be a probability kernel, and each selected πM\pi_MπM​ is checked to maximize value over valid randomized policies.

RewardsNormalized expresses both standing conditions of §5.2: reward values lie in [0,1][0,1][0,1], and under every randomized policy every trajectory of positive probability has total reward ∑h=1Hrh∈[0,1]\sum_{h=1}^{H}r_h\in[0,1]∑h=1H​rh​∈[0,1]. The proposition and its decoupling milestone use it. The Appendix F simulation lemmas and Lemma 5.2 use EpisodeTotalInUnitInterval, which expresses only the total-reward condition stated in those lemmas.

The target uses the full trajectory law M(π)M(\pi)M(π) in DH2(M(π),Mˉ(π))D_H^2(M(\pi),\bar M(\pi))DH2​(M(π),Mˉ(π)). Replacing it with one-step kernel distances would change Proposition 5.4. A valid model is a family of probability kernels, a valid policy is a family of action distributions, and the model selector must really choose an optimum; these checks prevent the expected inequality from becoming true through zero-weight pseudo-laws or a free, suboptimal policy selector. The definition and Hellinger lemmas can support later tabular DEC missions as well.

Selected references

  • Dylan J. Foster, Sham M. Kakade, Jian Qian, and Alexander Rakhlin, The Statistical Complexity of Interactive Decision Making, arXiv preprint, version 3, 2023. arXiv:2112.13487v3.
  • Dylan J. Foster and Alexander Rakhlin, Foundations of Reinforcement Learning and Interactive Decision Making, arXiv preprint, 2023. arXiv:2312.16730.
11 thms3 active usersReviewed
🏆Completed
Reinforcement LearningStatistics·Captain: mikedeng1

The Statistical Complexity of Interactive Decision Making IV: Log-Loss Regret Controls Hellinger Estimation Error, Est_H ≤ Reg_KL + 2 log(1/δ) with Probability 1−δResearch Paper

Motivation

Interactive decision methods often need an estimate of an unknown outcome law after each round. A learner can be good at decisions only if its model of the observations is accurate where future decisions will sample. Foster, Kakade, Qian and Rakhlin measure this accuracy by a cumulative squared Hellinger distance and use it in their regret bounds for decision making with structured observations. They connect that estimation quantity to online log-loss regret, a standard objective for sequential probability prediction. The connection lets a regret guarantee for a density predictor supply the estimation guarantee needed by a decision method. Foster et al., Appendix A.3 and §4.1.1.

This mission isolates the online estimation result, Lemma A.14, from the rest of the paper. It addresses a learner that predicts a conditional outcome distribution before the next covariate and outcome are revealed. The covariates may respond to the entire preceding history. This matters in interactive decision making, where the learner's earlier choices change what it observes next. Foster et al., pp. 75–76.

Setting

There are TTT rounds. At round ttt, the learner has a history H<t=((xs,ys))s<tH_{<t}=((x_s,y_s))_{s<t}H<t​=((xs​,ys​))s<t​ of past covariates and outcomes. It predicts a probability distribution g^t(⋅∣x;H<t)\widehat g_t(\cdot\mid x;H_{<t})g​t​(⋅∣x;H<t​) for every possible next covariate xxx. Nature draws xtx_txt​ from a history-dependent distribution ct(⋅∣H<t)c_t(\cdot\mid H_{<t})ct​(⋅∣H<t​), then draws yty_tyt​ from the conditional law of a distinguished true expert i⋆i^\stari⋆. An expert class indexed by III supplies laws gi,t(⋅∣x;H<t)g_{i,t}(\cdot\mid x;H_{<t})gi,t​(⋅∣x;H<t​); realizability means gi⋆,tg_{i^\star,t}gi⋆,t​ is the actual outcome law at every round. The prediction g^t\widehat g_tg​t​ is selected from the past history, before xtx_txt​ is observed. Foster et al., Assumption A.1, p. 76.

For a predicted distribution qqq, the log loss on observation yyy is −log⁡q(y)-\log q(y)−logq(y). The learner's log-loss regret RegKL\mathrm{Reg}_{\mathrm{KL}}RegKL​ is its cumulative log loss minus the smallest cumulative loss of an expert in III. This is a realized quantity: both sums use the covariates and outcomes that actually occurred. The squared Hellinger distance is

DH2(P,Q)=∑y(P(y)−Q(y))2.D_{\mathrm H}^2(P,Q)=\sum_y\bigl(\sqrt{P(y)}-\sqrt{Q(y)}\bigr)^2.DH2​(P,Q)=y∑​(P(y)​−Q(y)​)2.

The Hellinger estimation error EstH\mathrm{Est}_{\mathrm H}EstH​ is different from a sum evaluated only at realized covariates. At each round it averages DH2(g^t(⋅∣x),gi⋆,t(⋅∣x))D_{\mathrm H}^2(\widehat g_t(\cdot\mid x),g_{i^\star,t}(\cdot\mid x))DH2​(g​t​(⋅∣x),gi⋆,t​(⋅∣x)) over the conditional law of the next covariate xxx, given H<tH_{<t}H<t​, then sums these conditional errors. Thus the error measures prediction quality in the distribution the learner is about to face. Foster et al., (99)–(100), pp. 75–76.

Formalization targets

Main high-probability bound

The main target is Lemma A.14 (102): for every δ∈(0,1)\delta\in(0,1)δ∈(0,1),

Pr⁡ ⁣{EstH≤RegKL+2log⁡(1/δ)}≥1−δ.\Pr\!\left\{\mathrm{Est}_{\mathrm H} \le \mathrm{Reg}_{\mathrm{KL}}+2\log(1/\delta)\right\}\ge 1-\delta.Pr{EstH​≤RegKL​+2log(1/δ)}≥1−δ.

The coefficient 222 and the absence of a factor 1/21/21/2 in DH2D_{\mathrm H}^2DH2​ are part of the claim. The probability is under the complete history law generated by the adaptive covariate process and the true expert. Foster et al., Lemma A.14, p. 76.

Supporting results

The first milestone is Lemma A.4, an anytime conditional exponential-moment bound. It controls one event on which an adapted inequality fails at any horizon through TTT, without replacing it by a single fixed-horizon claim. The second milestone is the affinity identity used in the proof of Lemma A.14:

∑yP(y)Q(y)=1−12DH2(P,Q).\sum_y\sqrt{P(y)Q(y)}=1-\tfrac12D_{\mathrm H}^2(P,Q).y∑​P(y)Q(y)​=1−21​DH2​(P,Q).

The proposal also states the expectation clause E[EstH]≤E[RegKL]\mathbb E[\mathrm{Est}_{\mathrm H}]\le\mathbb E[\mathrm{Reg}_{\mathrm{KL}}]E[EstH​]≤E[RegKL​] from (101) and the predictable-mask clause (103), which bounds error accumulated only on rounds selected from the past. Foster et al., pp. 70, 76, 79.

Significance

The high-probability bound lets an online density estimator serve as an estimation oracle for the paper's interactive decision bounds. A user can choose a predictor because it has a log-loss regret guarantee, then obtain control of Hellinger estimation error under adaptive observations. The expectation clause supplies the corresponding average guarantee. The masked clause handles applications that evaluate accuracy only on rounds selected before the current observation. Foster et al., Proposition 4.1 and Lemma A.14.

The result is already proved in the source paper. This mission asks for machine-checked proofs of the finite-alphabet statements and a reusable interface for sequential probability kernels, conditional error and log-loss regret. The published FoundationsRL Hellinger definition supplies the same normalization used by the paper; the mission does not propose a new distance. These items have no machine-checked proof in this proposal yet. The source also gives Vovk's aggregation guarantee for finite expert classes, but that algorithmic bound is outside the current item list. Foster et al., Lemma A.15, p. 76.

Difficulty

The quantity to be bounded is conditional on the past. A pathwise distance at the observed xtx_txt​ is a different random quantity and does not give (102). Adaptive covariates add another dependence on history, while the learner must choose its prediction before seeing xtx_txt​. A fixed covariate law would therefore omit a central feature of the theorem. The anytime milestone likewise concerns all prefixes of one process at once; separate claims at each fixed horizon do not preserve its statement. Foster et al., pp. 70, 75–79.

The log-loss formula also has a boundary. A zero predicted probability gives infinite loss in the paper, whereas Lean's real logarithm returns a finite default at zero. The formal statements must keep that default from affecting a positive-probability observation. Squared Hellinger distance has its own normalization convention, so the affinity identity checks the factor needed in the high-probability constant.

Formalization scope

The paper uses measurable covariate and outcome spaces and regular conditional densities. This mission specializes to finite covariates XXX, finite outcomes YYY, and a finite nonempty expert class III. Each kernel is represented by a nonnegative vector summing to one. Complete histories are finite sequences, and the probability of a history is the product of its covariate and true-expert outcome probabilities. Expectations and event probabilities are finite sums. Time is indexed from zero; round ttt in Lean corresponds to round t+1t+1t+1 in the paper.

The learner, experts and covariate kernel may all depend on preceding history. The learner and experts predict a law for every possible covariate while seeing only that history. A Boolean mask for (103) has the same information restriction. The Hellinger error sums conditional averages over the next covariate, and the probability guarantee is expressed as a bound on the total mass of histories where the displayed inequality fails.

The real-valued log loss requires one additional condition: on every complete history of positive probability, the learner and every expert give the observed outcome positive probability at every round. The condition is imposed only where the history law has support. It excludes Lean's default value for log⁡0\log 0log0; it is a scope restriction from the paper's extended-valued formulation. A normalized history law and a distinguished true expert are explicit, so the probability statement cannot be discharged by an empty event space or an unnormalized zero law. The sequential-kernel definitions can be reused in other finite online learning results, and the probability, affinity, expectation and masked bounds are separate proof targets.

Selected references

  • Dylan J. Foster, Sham M. Kakade, Jian Qian, Alexander Rakhlin, The Statistical Complexity of Interactive Decision Making, arXiv preprint arXiv:2112.13487v3, 2023. arXiv.
8 thms3 active usersReviewed
Bandit AlgorithmsReinforcement Learning·Captain: mikedeng1

The Statistical Complexity of Interactive Decision Making V: The Multi-Armed Bandit Has Localized DEC at Least 2⁻⁶·A/γ, via Hard Families of ModelsResearch Paper

Motivation

The Decision-Estimation Coefficient (DEC) of Foster, Kakade, Qian and Rakhlin (arXiv:2112.13487) is a single complexity measure that bounds the minimax regret of interactive decision making (bandits, reinforcement learning, and their structured variants) from above and from below. Its lower bound (Theorem 3.2 of the paper) states that every algorithm suffers expected regret at least a constant times decγ(Mεγ∞(M‾),M‾)⋅T\mathsf{dec}_\gamma(\mathcal{M}^\infty_{\varepsilon_\gamma}(\overline{M}), \overline{M}) \cdot Tdecγ​(Mεγ​∞​(M),M)⋅T, so concrete lower bounds on regret follow from concrete lower bounds on the DEC of localized model classes.

The multi-armed bandit is the first test of this program. Its minimax regret Θ(AT)\Theta(\sqrt{AT})Θ(AT​) is classical (Auer et al., 2002; Audibert and Bubeck, 2010). Section 5.1 of the paper recovers it through the DEC: posterior sampling (Proposition 5.1) and inverse gap weighting (Proposition 5.2) bound the DEC by A/γA/\gammaA/γ from above, and Proposition 5.3 bounds it by 2−6A/γ2^{-6} A/\gamma2−6A/γ from below. The lower bound is proved through a hard family of models (Definition 5.1, Lemma 5.1), the device that the paper reuses for linear, ReLU and gap-dependent lower bounds in §§6–7.

Setting

A decision space Π\PiΠ is a finite set; here Π=[A]={1,…,A}\Pi = [A] = \{1, \dots, A\}Π=[A]={1,…,A}, the arms. A model MMM assigns to every arm π\piπ a distribution M(π)M(\pi)M(π) over rewards in R=[0,1]\mathcal{R} = [0,1]R=[0,1] (observations are trivial). Its mean reward is fM(π)=Er∼M(π)[r]f^M(\pi) = \mathbb{E}_{r \sim M(\pi)}[r]fM(π)=Er∼M(π)​[r], and πM∈arg⁡max⁡πfM(π)\pi_M \in \arg\max_\pi f^M(\pi)πM​∈argmaxπ​fM(π) is an optimal arm. The squared Hellinger distance between distributions P,QP, QP,Q on a finite set is DH2(P,Q)=∑y(P(y)−Q(y))2D^2_{\mathrm{H}}(P, Q) = \sum_y (\sqrt{P(y)} - \sqrt{Q(y)})^2DH2​(P,Q)=∑y​(P(y)​−Q(y)​)2 (no factor 12\tfrac1221​, the paper's (5)).

For a class M\mathcal{M}M of models, a reference model M‾\overline{M}M and γ>0\gamma > 0γ>0, the DEC (2) is

decγ(M,M‾)=inf⁡p∈Δ(Π)sup⁡M∈MEπ∼p[fM(πM)−fM(π)−γ⋅DH2(M(π),M‾(π))].\mathsf{dec}_\gamma(\mathcal{M}, \overline{M}) = \inf_{p \in \Delta(\Pi)} \sup_{M \in \mathcal{M}} \mathbb{E}_{\pi \sim p}\Bigl[ f^M(\pi_M) - f^M(\pi) - \gamma \cdot D^2_{\mathrm{H}}\bigl(M(\pi), \overline{M}(\pi)\bigr) \Bigr].decγ​(M,M)=p∈Δ(Π)inf​M∈Msup​Eπ∼p​[fM(πM​)−fM(π)−γ⋅DH2​(M(π),M(π))].

With gM(π)=fM(πM)−fM(π)g^M(\pi) = f^M(\pi_M) - f^M(\pi)gM(π)=fM(πM​)−fM(π), the L∞L_\inftyL∞​-localized class (12) is Mε∞(M‾)={M∈M:∣gM(π)−gM‾(π)∣≤ε ∀π}\mathcal{M}^\infty_\varepsilon(\overline{M}) = \{ M \in \mathcal{M} : |g^M(\pi) - g^{\overline{M}}(\pi)| \le \varepsilon\ \forall \pi \}Mε∞​(M)={M∈M:∣gM(π)−gM(π)∣≤ε ∀π}.

A reference model M‾∈M\overline{M} \in \mathcal MM∈M and models M1,…,MN∈MM_1, \dots, M_N \in \mathcal MM1​,…,MN​∈M, N≥2N \ge 2N≥2, form an (α,β,δ)(\alpha, \beta, \delta)(α,β,δ)-family (Definition 5.1) if there are ui,vi:Π→[0,1]u_i, v_i : \Pi \to [0,1]ui​,vi​:Π→[0,1] with ∑iui(π)≤N/2\sum_i u_i(\pi) \le N/2∑i​ui​(π)≤N/2, ∑ivi(π)≤1\sum_i v_i(\pi) \le 1∑i​vi​(π)≤1, fMi(πMi)−fMi(π)≥α(1−ui(π))f^{M_i}(\pi_{M_i}) - f^{M_i}(\pi) \ge \alpha(1 - u_i(\pi))fMi​(πMi​​)−fMi​(π)≥α(1−ui​(π)) and DH2(Mi(π),M‾(π))≤βvi(π)+δD^2_{\mathrm{H}}(M_i(\pi), \overline{M}(\pi)) \le \beta v_i(\pi) + \deltaDH2​(Mi​(π),M(π))≤βvi​(π)+δ. Ber(q)\mathrm{Ber}(q)Ber(q) is the Bernoulli distribution with mean qqq.

Formalization targets

Goal: Proposition 5.3

For A≥2A \ge 2A≥2 arms, M‾(π)=Ber(12)\overline{M}(\pi) = \mathrm{Ber}(\tfrac12)M(π)=Ber(21​) and every γ≥A/3\gamma \ge A/3γ≥A/3, with εγ=A12γ\varepsilon_\gamma = \frac{A}{12\gamma}εγ​=12γA​,

decγ(Mεγ∞(M‾),M‾)≥2−6⋅Aγ.\mathsf{dec}_\gamma\bigl(\mathcal{M}^\infty_{\varepsilon_\gamma}(\overline{M}), \overline{M}\bigr) \ge 2^{-6} \cdot \frac{A}{\gamma}.decγ​(Mεγ​∞​(M),M)≥2−6⋅γA​.

Milestones

  1. Lemma A.7: DH2(Ber(p),Ber(q))≤(p−q)2(1p+q+12−p−q)D^2_{\mathrm{H}}(\mathrm{Ber}(p), \mathrm{Ber}(q)) \le (p-q)^2\bigl(\frac{1}{p+q} + \frac{1}{2-p-q}\bigr)DH2​(Ber(p),Ber(q))≤(p−q)2(p+q1​+2−p−q1​), and in particular DH2(Ber(12+Δ),Ber(12))≤3Δ2D^2_{\mathrm{H}}(\mathrm{Ber}(\tfrac12 + \Delta), \mathrm{Ber}(\tfrac12)) \le 3\Delta^2DH2​(Ber(21​+Δ),Ber(21​))≤3Δ2 for Δ∈(0,12)\Delta \in (0, \tfrac12)Δ∈(0,21​).
  2. Lemma 5.1: an (α,β,δ)(\alpha, \beta, \delta)(α,β,δ)-family M′\mathcal{M}'M′ with β≥0\beta \ge 0β≥0 satisfies decγ(M′,M‾)≥α2−γ(βN+δ)\mathsf{dec}_\gamma(\mathcal{M}', \overline{M}) \ge \frac{\alpha}{2} - \gamma(\frac{\beta}{N} + \delta)decγ​(M′,M)≥2α​−γ(Nβ​+δ) for all γ>0\gamma > 0γ>0.
  3. The hard family of the proof: Mi(π)=Ber(12+Δ I{π=i})M_i(\pi) = \mathrm{Ber}(\tfrac12 + \Delta\,\mathbb{I}\{\pi = i\})Mi​(π)=Ber(21​+ΔI{π=i}) is a (Δ,3Δ2,0)(\Delta, 3\Delta^2, 0)(Δ,3Δ2,0)-family around Ber(12)\mathrm{Ber}(\tfrac12)Ber(21​).
  4. Localization: this family lies in MΔ∞(M‾)\mathcal{M}^\infty_\Delta(\overline{M})MΔ∞​(M).

Companions

The mission also states monotonicity of the DEC in the class (the step from the family to the localized class), Proposition 5.1 (posterior sampling bounds the dual square-loss DEC by A/γA/\gammaA/γ), and Proposition 6.7, the gap-dependent variant: for bandits with gap exactly Δ∈(0,18)\Delta \in (0, \tfrac18)Δ∈(0,81​), some M‾\overline{M}M in the class has decγ(MΔ∞(M‾),M‾)≥Δ4I{γ≤A48Δ}\mathsf{dec}_\gamma(\mathcal{M}^\infty_\Delta(\overline{M}), \overline{M}) \ge \frac{\Delta}{4}\mathbb{I}\{\gamma \le \frac{A}{48\Delta}\}decγ​(MΔ∞​(M),M)≥4Δ​I{γ≤48ΔA​}. Proposition 5.2 (inverse gap weighting) and the decoupling step of Proposition 5.1 are already proved on the platform and are included as references.

Significance

Proposition 5.3 is the DEC half of the Ω(AT)\Omega(\sqrt{AT})Ω(AT​) minimax lower bound for multi-armed bandits: combined with Theorem 3.2 and γ≍AT\gamma \asymp \sqrt{AT}γ≍AT​, it shows that the general DEC theory loses nothing in rate on the most basic interactive problem, and that the A/γA/\gammaA/γ upper bounds of Propositions 5.1 and 5.2 are tight up to constants. Lemma 5.1 is reusable: any lower bound on the DEC of a class can be obtained by exhibiting a hard family, and the paper uses it for linear bandits (Proposition 6.2), ReLU bandits (Proposition 6.6) and gap-dependent bounds (Propositions 6.7, 6.8).

The results are proved in the paper. As far as is known, none of them is machine-checked. A formal development provides a verified Bernoulli Hellinger bound, a reusable hard-family lemma for the DEC as published on the platform, and the first verified DEC lower bound for a concrete class.

Difficulty

The DEC is an inf-sup, so a lower bound must hold against every exploration distribution ppp. A natural first attempt picks a single alternative model and compares; this fails because ppp can concentrate on the arm where that alternative differs and then the Hellinger penalty cancels the regret. The family must be spread so that no single ppp is good for all of its members at once, and the constants must balance the regret term against the information term at the chosen Δ\DeltaΔ. In Lean, an additional difficulty is that the DEC is written with real sInf/sSup, which are 000 on empty or unbounded sets: every lower bound requires showing that the relevant sets are nonempty and bounded, and passing from the family to the localized class needs monotonicity of the DEC, which the paper uses without comment.

Formalization scope

  • Finite alphabets. The paper works over measurable spaces; this mission formalizes the finite-alphabet case, as the published FoundationsRL layer does. A model is m : S → Y → ℝ with every m π a probability vector; expectations are finite sums. The DEC is the published FoundationsRL.GeneralDM.decGf, the Hellinger distance hellingerSq, the mean reward fM.
  • Bernoulli subclass. The bandit class is the Bernoulli bandits (Y = Bool, reward 111 on true), a subclass of the paper's all-distributions-on-[0,1][0,1][0,1] class. The DEC is monotone in the class and the localized class shrinks with the class, so the Bernoulli lower bounds imply the paper's.
  • Optimal-arm selector. πM\pi_MπM​ is a parameter piStar with the hypothesis that it picks a maximizer of fMf^MfM; every statement holds for every such selector, so for whichever maximizer the paper meant (the reference Ber(12)\mathrm{Ber}(\tfrac12)Ber(21​) has every arm optimal).
  • Added or dropped hypotheses. Lemma 5.1 makes the page's nonempty decision space explicit; its optimal-arm selector also implies nonemptiness. It adds β≥0\beta \ge 0β≥0, necessary for the printed averaging argument: without it, even identical models can contradict (52). The family definition includes membership of the reference and family models in a probability-model class. The localization claim inherits A≥2A \ge 2A≥2 from Proposition 5.3. Proposition 6.7 states A≥2A \ge 2A≥2, implicit on the page since a gap needs two arms. Proposition 5.1 represents models by their mean vectors and priors as finitely supported (footnote 5).
  • Constants. All constants are the page's: 2−62^{-6}2−6, A/(12γ)A/(12\gamma)A/(12γ), γ≥A/3\gamma \ge A/3γ≥A/3, 3Δ23\Delta^23Δ2, 5Δ25\Delta^25Δ2, Δ/4\Delta/4Δ/4, A/(48Δ)A/(48\Delta)A/(48Δ). The page's proof of Proposition 6.7 names a (Δ,5Δ2,0)(\Delta, 5\Delta^2, 0)(Δ,5Δ2,0)-family but computes with 12Δ212\Delta^212Δ2; the stated proposition is the one its arithmetic supports.
  • No trivialization. A formalization in which the localized class could be empty, or the payoffs unbounded, would make decGf evaluate to the junk value 000; the goal is stated over a class containing M‾\overline{M}M with bounded rewards, so the inequality is the genuine one, and its right side 2−6A/γ2^{-6} A/\gamma2−6A/γ is strictly positive.

Contributions welcome: proofs of the milestones in attack order, sorry-free lemmas on nonemptiness and boundedness of the DEC's inner sets (reusable for every DEC lower bound in this series), and a proof of the gap-dependent companion.

Selected references

  • D. J. Foster, S. M. Kakade, J. Qian, A. Rakhlin, The Statistical Complexity of Interactive Decision Making, arXiv preprint, 2021 (v3, 2023). https://arxiv.org/abs/2112.13487
  • D. J. Foster, A. Rakhlin, Foundations of Reinforcement Learning and Interactive Decision Making, lecture notes, 2023. https://arxiv.org/abs/2312.16730
  • P. Auer, N. Cesa-Bianchi, Y. Freund, R. E. Schapire, The Nonstochastic Multiarmed Bandit Problem, SIAM Journal on Computing 32(1), 2002. https://doi.org/10.1137/S0097539701398375
  • J.-Y. Audibert, S. Bubeck, Regret Bounds and Minimax Policies under Partial Monitoring, Journal of Machine Learning Research 11, 2010. https://www.jmlr.org/papers/v11/audibert10a.html
  • D. J. Foster, A. Rakhlin, Beyond UCB: Optimal and Efficient Contextual Bandits with Regression Oracles, ICML, 2020. https://arxiv.org/abs/2002.04926
13 thms3 active usersReviewed
PreviousPage 8 of 27Next

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