Prove2Me
Navigate
DiscoverCollectionsFormalpediaBlogsUsersMomentumMy Missions+
Prove2Me
⌕
Log in
← Collections

The OR Formalization Drive

Help us formalize the operations research literature in Lean.

449 open missions

Missions

181–200 of 449
OpenCompletedAll
Bandit AlgorithmsMachine LearningOperations Research·Captain: mikedeng1

Regret Analysis of Stochastic and Nonstochastic Multi-armed Bandit Problems II: High-Probability and Expected Regret of Exp3.PTextbook

Motivation

In the adversarial (non-stochastic) multi-armed bandit problem a forecaster repeatedly chooses one of KKK actions while an opponent sets the rewards, and only the reward of the chosen action is revealed. The model was proposed as a way of playing an unknown repeated game: Baños (1968) studied the repeated game in which the player observes only its own payoff, which is exactly the bandit problem against an opponent who reacts to the player's past moves. It is the basic model of online decision making under partial feedback without statistical assumptions, and it underlies regret minimization in games, adversarial routing and online advertising. Chapter 3 of Bubeck and Cesa-Bianchi's monograph (arXiv:1204.5721v2) collects its fundamental results: the Exp3 forecaster of Auer, Cesa-Bianchi, Freund and Schapire (SIAM J. Comput. 2002), its high-probability variant Exp3.P, and the nK\sqrt{nK}nK​ minimax lower bound.

Setting

There are K≥2K \ge 2K≥2 arms and rounds t=1,2,…,nt = 1, 2, \dots, nt=1,2,…,n. At each round an adversary assigns a gain gi,t∈[0,1]g_{i,t} \in [0,1]gi,t​∈[0,1] to every arm iii; the forecaster picks an arm ItI_tIt​, possibly at random, and observes only gIt,tg_{I_t,t}gIt​,t​. The adversary may be non-oblivious (adaptive): gi,t=gi,t(I1,…,It−1)g_{i,t} = g_{i,t}(I_1,\dots,I_{t-1})gi,t​=gi,t​(I1​,…,It−1​) may depend on the forecaster's past actions. A forecaster rule maps the past actions to a probability vector ptp_tpt​ on the arms, and a run is a sequence of random arms with It∼ptI_t \sim p_tIt​∼pt​ given the past. The regret is the random variable

Rn=max⁡i=1,…,K∑t=1ngi,t−∑t=1ngIt,t,R_n = \max_{i=1,\dots,K}\sum_{t=1}^n g_{i,t} - \sum_{t=1}^n g_{I_t,t},Rn​=i=1,…,Kmax​t=1∑n​gi,t​−t=1∑n​gIt​,t​,

and, in the loss version ℓi,t∈[0,1]\ell_{i,t} \in [0,1]ℓi,t​∈[0,1], the pseudo-regret is R‾n=E∑tℓIt,t−min⁡iE∑tℓi,t\overline R_n = \mathbb E\sum_t \ell_{I_t,t} - \min_i \mathbb E\sum_t \ell_{i,t}Rn​=E∑t​ℓIt​,t​−mini​E∑t​ℓi,t​. Since the maximum sits inside the expectation, R‾n≤ERn\overline R_n \le \mathbb E R_nRn​≤ERn​ in the gain version, and against an adaptive adversary the two can differ.

Exp3 draws ItI_tIt​ from exponential weights pi,t+1∝exp⁡(−ηtL~i,t)p_{i,t+1} \propto \exp(-\eta_t \tilde L_{i,t})pi,t+1​∝exp(−ηt​L~i,t​) of importance-weighted cumulative loss estimates L~i,t=∑s≤tℓi,s1Is=i/pi,s\tilde L_{i,t} = \sum_{s \le t} \ell_{i,s}\mathbb 1_{I_s = i}/p_{i,s}L~i,t​=∑s≤t​ℓi,s​1Is​=i​/pi,s​. Exp3.P uses biased gain estimates g~i,t=(gi,t1It=i+β)/pi,t\tilde g_{i,t} = (g_{i,t}\mathbb 1_{I_t=i} + \beta)/p_{i,t}g~​i,t​=(gi,t​1It​=i​+β)/pi,t​ and mixes in the uniform distribution:

pi,t+1=(1−γ)exp⁡(ηG~i,t)∑kexp⁡(ηG~k,t)+γK,G~i,t=∑s=1tg~i,s.p_{i,t+1} = (1-\gamma)\frac{\exp(\eta\tilde G_{i,t})}{\sum_k \exp(\eta \tilde G_{k,t})} + \frac{\gamma}{K}, \qquad \tilde G_{i,t} = \sum_{s=1}^t \tilde g_{i,s}.pi,t+1​=(1−γ)∑k​exp(ηG~k,t​)exp(ηG~i,t​)​+Kγ​,G~i,t​=s=1∑t​g~​i,s​.

Formalization targets

Goal: Theorem 3.3 (expected regret of Exp3.P)

With β=ln⁡K/(nK)\beta = \sqrt{\ln K/(nK)}β=lnK/(nK)​, η=0.95ln⁡K/(nK)\eta = 0.95\sqrt{\ln K/(nK)}η=0.95lnK/(nK)​, γ=1.05Kln⁡K/n\gamma = 1.05\sqrt{K\ln K/n}γ=1.05KlnK/n​, against every adaptive adversary,

ERn≤5.15nKln⁡K+nKln⁡K.\mathbb E R_n \le 5.15\sqrt{nK\ln K} + \sqrt{\frac{nK}{\ln K}}.ERn​≤5.15nKlnK​+lnKnK​​.

Milestones

  • Lemma 3.1: for β∈(0,1]\beta \in (0,1]β∈(0,1] and a fixed arm iii, with probability at least 1−δ1-\delta1−δ, ∑tgi,t≤∑tg~i,t+ln⁡(δ−1)/β\sum_t g_{i,t} \le \sum_t \tilde g_{i,t} + \ln(\delta^{-1})/\beta∑t​gi,t​≤∑t​g~​i,t​+ln(δ−1)/β.
  • Eq. (3.12): if γ≤1/2\gamma \le 1/2γ≤1/2 and (1+β)Kη≤γ(1+\beta)K\eta \le \gamma(1+β)Kη≤γ, then with probability at least 1−δ1-\delta1−δ,
Rn≤βnK+γn+(1+β)ηKn+ln⁡(Kδ−1)β+ln⁡Kη.R_n \le \beta nK + \gamma n + (1+\beta)\eta Kn + \frac{\ln(K\delta^{-1})}{\beta} + \frac{\ln K}{\eta}.Rn​≤βnK+γn+(1+β)ηKn+βln(Kδ−1)​+ηlnK​.
  • Theorem 3.2: with β=ln⁡(Kδ−1)/(nK)\beta = \sqrt{\ln(K\delta^{-1})/(nK)}β=ln(Kδ−1)/(nK)​, Rn≤5.15nKln⁡(Kδ−1)R_n \le 5.15\sqrt{nK\ln(K\delta^{-1})}Rn​≤5.15nKln(Kδ−1)​ (3.10); with β=ln⁡K/(nK)\beta = \sqrt{\ln K/(nK)}β=lnK/(nK)​, Rn≤nK/ln⁡K ln⁡(δ−1)+5.15nKln⁡KR_n \le \sqrt{nK/\ln K}\,\ln(\delta^{-1}) + 5.15\sqrt{nK\ln K}Rn​≤nK/lnK​ln(δ−1)+5.15nKlnK​ (3.11), each with probability at least 1−δ1-\delta1−δ.
  • Theorem 3.1: Exp3 with η=2ln⁡K/(nK)\eta = \sqrt{2\ln K/(nK)}η=2lnK/(nK)​ has R‾n≤2nKln⁡K\overline R_n \le \sqrt{2nK\ln K}Rn​≤2nKlnK​ (3.2); with ηt=ln⁡K/(tK)\eta_t = \sqrt{\ln K/(tK)}ηt​=lnK/(tK)​, R‾n≤2nKln⁡K\overline R_n \le 2\sqrt{nK\ln K}Rn​≤2nKlnK​ (3.3).
  • Lemma 3.2 and Theorem 3.4: for n≥K≥2n \ge K \ge 2n≥K≥2 and every forecaster there is a Bernoulli instance with max⁡iE∑tYi,t−E∑tYIt,t≥nK/20\max_i \mathbb E\sum_t Y_{i,t} - \mathbb E\sum_t Y_{I_t,t} \ge \sqrt{nK}/20maxi​E∑t​Yi,t​−E∑t​YIt​,t​≥nK​/20.

Significance

The goal bounds the expected regret, not the pseudo-regret, against an opponent that adapts to the forecaster's randomized past choices. A pseudo-regret bound says nothing about ERn\mathbb E R_nERn​ in that setting, and the book obtains the expected-regret bound by first proving a high-probability bound valid at every confidence level, (3.11), and integrating its tail. Together with Theorem 3.4 the chapter shows that nK\sqrt{nK}nK​ is the minimax rate of adversarial bandits up to a ln⁡K\sqrt{\ln K}lnK​ factor. Lemma 3.1, the concentration of biased importance-weighted estimates, holds for any forecaster rule and is the step that turns exponential weights into a high-probability guarantee.

All results are proved in the book. On the formal side, the platform has the pseudo-regret bound of Exp3 against an oblivious adversary (a fixed reward table, Bandit Algorithms V) and an Exp3-IX high-probability bound; it has no Exp3.P, no regret bound against adaptive adversaries and no Bernoulli nK/20\sqrt{nK}/20nK​/20 lower bound. This mission adds an explicit model of adaptive adversaries and randomized forecaster runs, and the chapter's statements with the book's exact constants.

Difficulty

Against an adaptive adversary the gains are random and depend on the forecaster's own past draws, so the argument used for a fixed reward table (take expectations of an inequality that holds for every fixed sequence) does not control ERn\mathbb E R_nERn​: the maximum over arms does not commute with the expectation. Unbiased estimates do not help either, because the variance of ℓi,t/pi,t\ell_{i,t}/p_{i,t}ℓi,t​/pi,t​ is of order 1/pi,t1/p_{i,t}1/pi,t​, which can be arbitrarily large; even with uniform mixing at rate n−1/2n^{-1/2}n−1/2 the cumulative variance is of order n3/2n^{3/2}n3/2. The bias β\betaβ and the mixing γ\gammaγ have to be tuned jointly so that the estimate concentrates while the exponential-weights analysis survives, and the constants 0.950.950.95, 1.051.051.05 and 5.155.155.15 come out of that tuning. The lower bound needs an information-theoretic comparison of a forecaster's behaviour on K+1K+1K+1 Bernoulli instances, against forecasters that may be randomized.

Formalization scope

Arms are Fin K with K≥2K \ge 2K≥2; rounds are numbered 1,…,n1,\dots,n1,…,n; logarithms are natural. Action sequences are functions N→\mathbb N \toN→ Fin K whose entry 000 is ignored. An adversary is a structure holding values in [0,1][0,1][0,1] that may depend on the past actions only (gains for Exp3.P, losses for Exp3); a randomized adversary with independent external randomness reduces to this case by conditioning. A run of a forecaster rule ppp on a probability space is pinned down by the cylinder identity P(I1=h1,…,It=ht)=P(I1=h1,…,It−1=ht−1) pt(h)(ht)\mathbb P(I_1 = h_1,\dots,I_t = h_t) = \mathbb P(I_1=h_1,\dots,I_{t-1}=h_{t-1})\,p_t(h)(h_t)P(I1​=h1​,…,It​=ht​)=P(I1​=h1​,…,It−1​=ht−1​)pt​(h)(ht​), which determines the law of (I1,…,In)(I_1,\dots,I_n)(I1​,…,In​). "With probability at least 1−δ1-\delta1−δ" is P(event)≥1−δ\mathbb P(\text{event}) \ge 1-\deltaP(event)≥1−δ for δ∈(0,1)\delta \in (0,1)δ∈(0,1), and ERn\mathbb E R_nERn​ is the Bochner integral of the bounded, measurable regret. The lower bounds use a stochastic model in which the forecaster sees past actions and the rewards of the played arms, and rewards are i.i.d. product Bernoulli.

Constants and conventions:

  • Every constant is the book's exact one: 0.950.950.95, 1.051.051.05, 5.155.155.15, 1/201/201/20. No O(⋅)O(\cdot)O(⋅) is involved.
  • Exp3.P with 1.05Kln⁡K/n>11.05\sqrt{K\ln K/n} > 11.05KlnK/n​>1 is outside the box's range γ∈[0,1]\gamma \in [0,1]γ∈[0,1]; its vector can then have negative entries, and if it does on a history of positive probability no run exists. This happens only when n<1.11 Kln⁡Kn < 1.11\,K\ln Kn<1.11KlnK, where the printed bounds already follow from Rn≤nR_n \le nRn​≤n, so the statements are true there whether or not a run exists.
  • Corrected misprints: the Exp3 box's ℓ~i,s\tilde\ell_{i,s}ℓ~i,s​ is ℓ~i,t\tilde\ell_{i,t}ℓ~i,t​; the sign in (3.16) is the box's exp⁡(+ηG~)\exp(+\eta\tilde G)exp(+ηG~); the proof of (3.10) says the bound is trivial "if n≥5.15⋯n \ge 5.15\sqrt{\cdots}n≥5.15⋯​", which should be n≤n \len≤. The statements carry no lower bound on nnn.
  • Added standing hypotheses: K≥2K \ge 2K≥2 everywhere, n≥Kn \ge Kn≥K in Theorem 3.4 (from the protocol box, p. 6; Theorem 3.4 is false without it), β>0\beta > 0β>0 and pi,t>0p_{i,t} > 0pi,t​>0 in Lemma 3.1.
  • Theorem 3.4 is stated as "for every forecaster there is a Bernoulli instance with regret at least nK/20\sqrt{nK}/20nK​/20", which implies the book's inf⁡sup⁡\inf\supinfsup (3.18).

A trivializing formalization is ruled out: the forecasters are fixed rules of the observed history drawn with fresh randomness, the adversary is not restricted to a fixed sequence, and the lower bounds quantify over all forecasters and exhibit the instance.

Welcome contributions: a reusable construction of runs (existence of a probability space carrying a run for every rule), the supermartingale form of Lemma 3.1, the exponential-weights potential argument, a tail-integration lemma EW≤∫01δ−1P(W>ln⁡δ−1) dδ\mathbb E W \le \int_0^1 \delta^{-1}\mathbb P(W > \ln\delta^{-1})\,d\deltaEW≤∫01​δ−1P(W>lnδ−1)dδ, and a KL/Pinsker comparison for bandit runs.

Selected references

  • S. Bubeck, N. Cesa-Bianchi, Regret Analysis of Stochastic and Nonstochastic Multi-armed Bandit Problems, Foundations and Trends in Machine Learning 5(1), 2012. arXiv:1204.5721v2, doi:10.1561/2200000024
  • P. Auer, N. Cesa-Bianchi, Y. Freund, R. E. Schapire, The nonstochastic multiarmed bandit problem, SIAM Journal on Computing 32(1), 2002. doi:10.1137/S0097539701398375
  • J.-Y. Audibert, S. Bubeck, Regret bounds and minimax policies under partial monitoring, Journal of Machine Learning Research 11, 2010. jmlr.org/papers/v11/audibert10a
  • N. Cesa-Bianchi, G. Lugosi, Prediction, Learning, and Games, Cambridge University Press, 2006. doi:10.1017/CBO9780511546921
13 thms1 active userReviewed
Discrete GeometryLinear OptimizationOperations Research+1·Captain: mikedeng1

Sensitivity Theorems in Integer Linear Programming: Every Integral m×n Matrix Has Chvátal Rank at Most 2^(n³+1)·n^(5n)·Δ(A)^(n+1)Research Paper

Motivation

An integer linear program max⁡{wx:Ax≤b, x integral}\max\{wx : Ax \le b,\ x \text{ integral}\}max{wx:Ax≤b, x integral} is usually attacked through its linear programming relaxation max⁡{wx:Ax≤b}\max\{wx : Ax \le b\}max{wx:Ax≤b}, which drops the integrality constraint. Two questions follow at once. How far can an optimal solution of the relaxation be from an optimal integer solution? And how many rounds of rounding-based cutting planes are needed before the relaxation describes the integer points exactly? Branch-and-bound, cutting-plane methods and the parametric analysis of integer programs all depend on the answers.

W. Cook, A.M.H. Gerards, A. Schrijver and É. Tardos, Sensitivity theorems in integer linear programming (Math. Programming 34 (1986) 251–264), answer both in terms of the number of variables nnn and the largest subdeterminant Δ(A)\Delta(A)Δ(A) of the constraint matrix, independently of the right-hand side.

Timeline.

  • 1958–1963: Gomory introduces integer rounding cuts. In 1973 Chvátal (Discrete Math. 4) shows that finitely many rounds reach the integer hull of a bounded polyhedron.
  • 1977–1979: Blair and Jeroslow prove that for a fixed matrix AAA the distance between LP and IP optima, and the gap between their values, are bounded by constants depending on AAA.
  • 1980: Schrijver (Ann. Discrete Math. 9) proves that the Chvátal closure of a rational polyhedron is a polyhedron, and that every rational polyhedron, bounded or not, reaches its integer hull after finitely many rounds.
  • 1986: Cook, Gerards, Schrijver and Tardos prove the explicit bounds of this mission, nΔ(A)n\Delta(A)nΔ(A) for proximity, and show that every integral matrix has finite Chvátal rank.
  • Later work, for example Eisenbrand and Weismantel (2018), replaces the ℓ∞\ell_\inftyℓ∞​ proximity bound by ℓ1\ell_1ℓ1​ bounds for programs in standard form.

Setting

All matrices, vectors and polyhedra are rational. Let AAA be an integral m×nm\times nm×n matrix. A square submatrix of order kkk, where 1≤k≤min⁡(m,n)1\le k\le\min(m,n)1≤k≤min(m,n), keeps kkk rows and kkk columns of AAA. The quantity Δ(A)\Delta(A)Δ(A) is the largest ∣det⁡B∣|\det B|∣detB∣ over all such submatrices BBB. So Δ(0)=0\Delta(0)=0Δ(0)=0, and Δ(A)≥1\Delta(A)\ge 1Δ(A)≥1 whenever A≠0A\ne 0A=0. Norms are ∥x∥∞=max⁡i∣xi∣\|x\|_\infty=\max_i|x_i|∥x∥∞​=maxi​∣xi​∣ and ∥x∥1=∑i∣xi∣\|x\|_1=\sum_i|x_i|∥x∥1​=∑i​∣xi​∣.

For b∈Qmb\in\mathbb{Q}^mb∈Qm write P={x∈Qn:Ax≤b}P=\{x\in\mathbb{Q}^n : Ax\le b\}P={x∈Qn:Ax≤b}. An optimal solution of max⁡{wx:Ax≤b}\max\{wx : Ax\le b\}max{wx:Ax≤b} is a point of PPP maximizing wxwxwx. For max⁡{wx:Ax≤b, x integral}\max\{wx : Ax\le b,\ x\text{ integral}\}max{wx:Ax≤b, x integral} it is an integral point of PPP maximizing wxwxwx among the integral points of PPP. A rational polyhedron is a set {x:Dx≤d}\{x : Dx\le d\}{x:Dx≤d} with DDD, ddd rational. The integer hull PIP_IPI​ is the convex hull of the integral points of PPP.

If ay≤βay\le\betaay≤β for all y∈Py\in Py∈P, with aaa integral and β\betaβ rational, then every integral point of PPP satisfies the Chvátal cut ax≤⌊β⌋ax\le\lfloor\beta\rfloorax≤⌊β⌋. The Chvátal closure P′P'P′ is the set of points satisfying all Chvátal cuts. Set P(0)=PP^{(0)}=PP(0)=P and P(i)=(P(i−1))′P^{(i)}=(P^{(i-1)})'P(i)=(P(i−1))′. Then PI⊆P(i)P_I\subseteq P^{(i)}PI​⊆P(i) for all iii. The Chvátal rank of PPP is the least ttt with P(t)=PIP^{(t)}=P_IP(t)=PI​. The Chvátal rank of the matrix AAA is the supremum of the Chvátal ranks of {x:Ax≤b}\{x : Ax\le b\}{x:Ax≤b} over all integral vectors bbb.

Formalization targets

Goal: Theorem 10 (p. 260)

sup⁡b∈Zm rank⁡{x:Ax≤b} ≤ 2n3+1 n5n Δ(A)n+1.\sup_{b\in\mathbb{Z}^m}\ \operatorname{rank}\{x : Ax\le b\}\ \le\ 2^{n^3+1}\,n^{5n}\,\Delta(A)^{n+1}.b∈Zmsup​ rank{x:Ax≤b} ≤ 2n3+1n5nΔ(A)n+1.

In particular, every integral matrix has finite Chvátal rank, and the bound does not depend on mmm or on bbb.

Milestones, in attack order

  1. Theorem 1 (p. 252). Suppose Ax≤bAx\le bAx≤b has an integral solution and the LP maximum exists. Then every LP optimum has an IP optimum within ℓ∞\ell_\inftyℓ∞​-distance nΔ(A)n\Delta(A)nΔ(A), and every IP optimum has an LP optimum within the same distance.
  2. Corollary 2 (p. 253). Under the same hypotheses, max⁡{wx:Ax≤b}−max⁡{wx:Ax≤b, x integral}≤nΔ(A)∥w∥1\max\{wx: Ax\le b\}-\max\{wx : Ax\le b,\ x\text{ integral}\}\le n\Delta(A)\|w\|_1max{wx:Ax≤b}−max{wx:Ax≤b, x integral}≤nΔ(A)∥w∥1​.
  3. Theorem 5 (p. 255). Changing bbb to b′b'b′ moves LP optima by at most nΔ(A)∥b−b′∥∞n\Delta(A)\|b-b'\|_\inftynΔ(A)∥b−b′∥∞​ and IP optima by at most nΔ(A)(∥b−b′∥∞+2)n\Delta(A)(\|b-b'\|_\infty+2)nΔ(A)(∥b−b′∥∞​+2). This result is off the goal's path.
  4. Theorem 6 (p. 256). A non-optimal integral solution can be improved by an integral solution within ℓ∞\ell_\inftyℓ∞​-distance nΔ(A)n\Delta(A)nΔ(A).
  5. Theorem 7 (p. 257). A single integral matrix MMM, with entries at most n2nΔ(A)nn^{2n}\Delta(A)^nn2nΔ(A)n in absolute value, gives {x:Ax≤b}I={x:Mx≤db}\{x: Ax\le b\}_I=\{x : Mx\le d_b\}{x:Ax≤b}I​={x:Mx≤db​} for every bbb for which Ax≤bAx\le bAx≤b has an integral solution.
  6. Theorem 8, printed "Theorem 9" (p. 259). If a rational polyhedron P⊆QnP\subseteq\mathbb{Q}^nP⊆Qn has no integral point, then P(n2n2n3)=∅P^{(n^{2n}2^{n^3})}=\emptysetP(n2n2n3)=∅.
  7. Corollary 9 (p. 260). Let q=max⁡{wx:x∈PI}q=\max\{wx : x\in P_I\}q=max{wx:x∈PI​} with www integral. Then P(r)⊆{x:wx≤q}P^{(r)}\subseteq\{x : wx\le q\}P(r)⊆{x:wx≤q} for r=(n2n2n3+1)(⌊max⁡{wx:x∈P}⌋−q)+1r=(n^{2n}2^{n^3}+1)(\lfloor\max\{wx : x\in P\}\rfloor-q)+1r=(n2n2n3+1)(⌊max{wx:x∈P}⌋−q)+1.

Significance

The result. Theorem 10 shows that the number of Gomory–Chvátal rounding rounds needed for {x:Ax≤b}\{x : Ax\le b\}{x:Ax≤b} is controlled by AAA alone. It is the first general finite bound on the Chvátal rank of a matrix. Earlier, the matrices of Chvátal rank 0 had been characterized by Hoffman and Kruskal: they are the matrices whose transpose is unimodular. Some classes of rank 1 had also been characterized (Edmonds–Johnson, Gerards–Schrijver). The proximity results of §2 are used on their own. They bound the work needed to solve an integer program from an LP optimum, and they show that the optimal value of an integer program changes at most affinely with bbb. They are also the standard starting point for the later proximity literature.

Formalizing it. All results are proved in the paper. As far as is known, none of them has a machine-checked proof: the Prove2Me corpus holds no Chvátal rank bound, and its existing proximity theorems concern a different bound, the ℓ1\ell_1ℓ1​ bound with Δ\DeltaΔ the largest entry. This mission asks for Lean proofs of the paper's statements with the constants exactly as printed. It also builds a reusable layer over Q\mathbb{Q}Q: polyhedra, LP and IP optimality, integer hulls, the Chvátal closure and the Chvátal rank.

Difficulty

The proximity theorems need a conic decomposition xˉ−zˉ=∑λigi\bar x-\bar z=\sum\lambda_i g^ixˉ−zˉ=∑λi​gi into integral generators with entries bounded by Δ(A)\Delta(A)Δ(A). That requires Cramer's rule bounds on cone generators and Carathéodory's theorem, and neither is in Mathlib in this form for rational polyhedral cones.

Theorem 7 needs finite generation of integral cones with explicit coefficient bounds, together with LP duality.

The Chvátal-rank part is harder. The obvious induction on the value of a valid inequality fails, because the value gap ⌊max⁡Pwx⌋−q\lfloor\max_P wx\rfloor-q⌊maxP​wx⌋−q is not bounded independently of bbb until Theorem 7 and Corollary 2 bound it by n2n+2Δ(A)n+1n^{2n+2}\Delta(A)^{n+1}n2n+2Δ(A)n+1. Theorem 8 itself rests on a flatness theorem for lattice-free polyhedra (Lenstra; Grötschel–Lovász–Schrijver), which the paper cites without proof. It also needs Schrijver's lemma that P(k)∩F⊆F(k)P^{(k)}\cap F\subseteq F^{(k)}P(k)∩F⊆F(k) for faces FFF, and invariance under unimodular affine maps. None of these is in Mathlib.

Formalization scope

  • Rationality. Everything is over Q\mathbb{Q}Q, following the paper's standing assumption on p. 252. Points are Fin n → ℚ, AAA is Matrix (Fin m) (Fin n) ℤ cast to Q\mathbb{Q}Q, and a polyhedron is a finite system of rational inequalities.
  • Δ(A)\Delta(A)Δ(A). Only nonempty submatrices count, so Δ(0)=0\Delta(0)=0Δ(0)=0.
  • Optimality. "The maximum exists" means an optimal solution exists. Existence claims that the paper proves are part of the conclusions: the IP optimum in Theorem 1 and Corollary 2, and max⁡{wx:x∈P}\max\{wx : x\in P\}max{wx:x∈P} in Corollary 9.
  • Chvátal closure. It is defined for every subset of Qn\mathbb{Q}^nQn, using all integral aaa and rational β\betaβ. The rank is valued in N∪{∞}\mathbb{N}\cup\{\infty\}N∪{∞}, with ∞\infty∞ if no iterate equals PIP_IPI​. A version with junk value 000 would make the goal trivial and is not used. The matrix rank is a supremum over integral bbb, as printed.
  • Added hypotheses. Theorems 1 and 6 carry the hypothesis A≠0A\ne0A=0. For A=0A=0A=0 the bound nΔ(A)=0n\Delta(A)=0nΔ(A)=0 makes both statements false, and the proof on p. 257 assumes A≠0A\ne0A=0 as well. In Corollary 9 the value qqq is taken to be an integer. This loses nothing, because a maximum of an integral www over PIP_IPI​ is attained at an integral point.
  • Constants. All constants are exactly as printed, written in N\mathbb{N}N with 00=10^0=100=1.

A complete development needs the following:

  • cone generation with Cramer bounds and Carathéodory's theorem;
  • LP duality and Farkas' lemma over Q\mathbb{Q}Q;
  • the polyhedrality of P′P'P′ for rational polyhedra (Schrijver 1980);
  • Schrijver's face lemma and unimodular invariance;
  • a flatness theorem.

The LP, cone and Chvátal-closure layers are reusable beyond this mission. Proofs of any milestone are welcome, and so is groundwork such as polyhedrality of the Chvátal closure or the flatness theorem, submitted as separate theorems.

Selected references

  • W. Cook, A.M.H. Gerards, A. Schrijver, É. Tardos, Sensitivity theorems in integer linear programming, Mathematical Programming 34 (1986) 251–264. https://doi.org/10.1007/BF01582230
  • V. Chvátal, Edmonds polytopes and a hierarchy of combinatorial problems, Discrete Mathematics 4 (1973) 305–337. https://doi.org/10.1016/0012-365X(73)90167-2
  • A. Schrijver, On cutting planes, Annals of Discrete Mathematics 9 (1980) 291–296. https://doi.org/10.1016/S0167-5060(08)70085-2
  • W. Cook, C.R. Coullard, Gy. Turán, On the complexity of cutting-plane proofs, Discrete Applied Mathematics 18 (1987) 25–38. https://doi.org/10.1016/0166-218X(87)90039-4
  • F. Eisenbrand, R. Weismantel, Proximity results and faster algorithms for integer programming using the Steinitz lemma, ACM Transactions on Algorithms 16 (2020), Art. 5. https://doi.org/10.1145/3340322
11 thms1 active userReviewed
Bandit AlgorithmsConvex OptimizationMachine Learning+2·Captain: mikedeng1

Regret Analysis of Stochastic and Nonstochastic Multi-armed Bandit Problems IV: Online Stochastic Mirror Descent for Combinatorial Semi-BanditsTextbook

Motivation

Many sequential decision problems ask a learner to choose, round after round, a combination of items: a set of mmm ads out of ddd, a path in a network, a matching. After each choice the learner sees the loss of the items it used, not of those it did not. This is online combinatorial optimization with semi-bandit feedback. It contains the classical adversarial multi-armed bandit (choose one of ddd arms) and is a standard model in online advertising, routing and ranking.

Chapter 5 of Bubeck and Cesa-Bianchi's monograph arXiv:1204.5721v2 treats this problem with one algorithm, Online Stochastic Mirror Descent (OSMD). Every regret bound in the chapter comes from a single mirror-descent inequality, specialized through the choice of a convex "regularizer". The chapter's capstone, Theorem 5.7, shows that a polynomial regularizer gives pseudo-regret O(mdn)O(\sqrt{mdn})O(mdn​) with no logarithmic factor. For m=1m=1m=1 this is the minimax-optimal rate of the adversarial bandit, first attained by the INF strategy of Audibert and Bubeck (2009). The semi-bandit version is due to Audibert, Bubeck and Lugosi (2014).

Setting

Vectors live in Rd\mathbb R^dRd. The arm set is a nonempty C⊆{0,1}d\mathcal C\subseteq\{0,1\}^dC⊆{0,1}d with ∥v∥1=m\|v\|_1=m∥v∥1​=m for every v∈Cv\in\mathcal Cv∈C, and K=Conv(C)\mathcal K=\mathrm{Conv}(\mathcal C)K=Conv(C). An oblivious adversary fixes loss vectors ℓ1,…,ℓn∈[0,1]d\ell_1,\dots,\ell_n\in[0,1]^dℓ1​,…,ℓn​∈[0,1]d. In round ttt the learner plays a random arm vt∈Cv_t\in\mathcal Cvt​∈C, pays ℓt⊤vt\ell_t^\top v_tℓt⊤​vt​, and observes (ℓt(1)vt(1),…,ℓt(d)vt(d))(\ell_t(1)v_t(1),\dots,\ell_t(d)v_t(d))(ℓt​(1)vt​(1),…,ℓt​(d)vt​(d)). The pseudo-regret is

Rˉn=E∑t=1nℓt⊤vt−min⁡x∈K∑t=1nℓt⊤x.\bar R_n=\mathbb E\sum_{t=1}^n\ell_t^\top v_t-\min_{x\in\mathcal K}\sum_{t=1}^n\ell_t^\top x .Rˉn​=Et=1∑n​ℓt⊤​vt​−x∈Kmin​t=1∑n​ℓt⊤​x.

A Legendre function on Dˉ\bar DDˉ, for a nonempty open convex DDD, is a continuous F:Dˉ→RF:\bar D\to\mathbb RF:Dˉ→R that is strictly convex and C1C^1C1 on DDD and whose gradient norm tends to +∞+\infty+∞ at Dˉ∖D\bar D\setminus DDˉ∖D. Its Bregman divergence is DF(x,y)=F(x)−F(y)−(x−y)⊤∇F(y)D_F(x,y)=F(x)-F(y)-(x-y)^\top\nabla F(y)DF​(x,y)=F(x)−F(y)−(x−y)⊤∇F(y), and its Legendre–Fenchel transform is F∗(u)=sup⁡x∈Dˉ(x⊤u−F(x))F^*(u)=\sup_{x\in\bar D}(x^\top u-F(x))F∗(u)=supx∈Dˉ​(x⊤u−F(x)).

Online Mirror Descent with learning rate η>0\eta>0η>0 and vectors gtg_tgt​ starts at x1∈arg⁡min⁡KFx_1\in\arg\min_{\mathcal K}Fx1​∈argminK​F. It then sets ∇F(wt+1)=∇F(xt)−ηgt\nabla F(w_{t+1})=\nabla F(x_t)-\eta g_t∇F(wt+1​)=∇F(xt​)−ηgt​ and xt+1=arg⁡min⁡y∈KDF(y,wt+1)x_{t+1}=\arg\min_{y\in\mathcal K}D_F(y,w_{t+1})xt+1​=argminy∈K​DF​(y,wt+1​). OSMD uses a random estimate gt=ℓ~tg_t=\tilde\ell_tgt​=ℓ~t​ of the loss. In the semi-bandit case it plays vtv_tvt​ with E[vt∣xt]=xt\mathbb E[v_t\mid x_t]=x_tE[vt​∣xt​]=xt​ and uses

ℓ~t(i)=ℓt(i) vt(i)xt(i).(5.5)\tilde\ell_t(i)=\frac{\ell_t(i)\,v_t(i)}{x_t(i)}. \tag{5.5}ℓ~t​(i)=xt​(i)ℓt​(i)vt​(i)​.(5.5)

A 000-potential is a convex, C1C^1C1, increasing ψ:(−∞,a)→(0,∞)\psi:(-\infty,a)\to(0,\infty)ψ:(−∞,a)→(0,∞) with ψ(−∞)=0\psi(-\infty)=0ψ(−∞)=0, ψ(a−)=+∞\psi(a^-)=+\inftyψ(a−)=+∞ and ∫01∣ψ−1∣<∞\int_0^1|\psi^{-1}|<\infty∫01​∣ψ−1∣<∞. It defines the Legendre function Fψ(x)=∑i∫0xiψ−1(s) dsF_\psi(x)=\sum_i\int_0^{x_i}\psi^{-1}(s)\,dsFψ​(x)=∑i​∫0xi​​ψ−1(s)ds on [0,∞)d[0,\infty)^d[0,∞)d. With ψ=exp⁡\psi=\expψ=exp this is the negative entropy.

Formalization targets

Goal: Theorem 5.7 (p. 80)

For every 000-potential ψ\psiψ and non-negative unbiased estimates,

Rˉn≤sup⁡KFψ−Fψ(x1)η+η2∑t=1n∑i=1dE[ℓ~t(i)2(ψ−1)′(xt(i))].\bar R_n\le\frac{\sup_{\mathcal K}F_\psi-F_\psi(x_1)}{\eta}+\frac\eta2\sum_{t=1}^n\sum_{i=1}^d\mathbb E\left[\frac{\tilde\ell_t(i)^2}{(\psi^{-1})'(x_t(i))}\right].Rˉn​≤ηsupK​Fψ​−Fψ​(x1​)​+2η​t=1∑n​i=1∑d​E[(ψ−1)′(xt​(i))ℓ~t​(i)2​].

For ψ(x)=(−x)−q\psi(x)=(-x)^{-q}ψ(x)=(−x)−q with q>1q>1q>1, the estimate (5.5) and η=2q−1 m1−2/q/(n d1−2/q)\eta=\sqrt{\tfrac{2}{q-1}\,m^{1-2/q}/(n\,d^{1-2/q})}η=q−12​m1−2/q/(nd1−2/q)​,

Rˉn≤q2q−1 mdn,and  Rˉn≤22mdn  at q=2.\bar R_n\le q\sqrt{\tfrac{2}{q-1}\,mdn},\qquad\text{and }\ \bar R_n\le2\sqrt{2mdn}\ \text{ at }q=2.Rˉn​≤qq−12​mdn​,and  Rˉn​≤22mdn​  at q=2.

Milestones

  1. Lemma 5.1: F∗∗=FF^{**}=FF∗∗=F, ∇F∗=(∇F)−1\nabla F^*=(\nabla F)^{-1}∇F∗=(∇F)−1 on D∗D^*D∗, and DF(x,y)=DF∗(∇F(y),∇F(x))D_F(x,y)=D_{F^*}(\nabla F(y),\nabla F(x))DF​(x,y)=DF∗​(∇F(y),∇F(x)).
  2. Lemma 5.2: existence, uniqueness and the Pythagorean inequality of Bregman projections.
  3. Theorem 5.3: ∑tℓt(xt)−∑tℓt(x)≤F(x)−F(x1)η+1η∑tDF∗(∇F(xt)−η∇ℓt(xt),∇F(xt))\sum_t\ell_t(x_t)-\sum_t\ell_t(x)\le\frac{F(x)-F(x_1)}\eta+\frac1\eta\sum_tD_{F^*}(\nabla F(x_t)-\eta\nabla\ell_t(x_t),\nabla F(x_t))∑t​ℓt​(xt​)−∑t​ℓt​(x)≤ηF(x)−F(x1​)​+η1​∑t​DF∗​(∇F(xt​)−η∇ℓt​(xt​),∇F(xt​)).
  4. Theorem 5.5, linear losses, and its corrected general form.
  5. Lemma 5.3: FψF_\psiFψ​ is Legendre and DFψ∗(u,v)≤12∑iψ′(vi)(ui−vi)2D_{F_\psi^*}(u,v)\le\frac12\sum_i\psi'(v_i)(u_i-v_i)^2DFψ∗​​(u,v)≤21​∑i​ψ′(vi​)(ui​−vi​)2 for u≤vu\le vu≤v.
  6. Theorem 5.6: with the negative entropy, Rˉn≤2mdnln⁡(d/m)\bar R_n\le\sqrt{2mdn\ln(d/m)}Rˉn​≤2mdnln(d/m)​.

Significance

Theorem 5.7 is the sharpest semi-bandit bound in the monograph. It shows that removing the ln⁡(d/m)\sqrt{\ln(d/m)}ln(d/m)​ factor of the exponential-weights analysis (Theorem 5.6) is a matter of the regularizer, not of a new algorithm. The same OSMD template gives the Euclidean-ball bound of Theorem 5.8 and is reused for bandit convex optimization in Chapter 6. Lemma 5.1, Lemma 5.2 and Theorem 5.3 are the standard mirror-descent toolkit, used throughout online learning and optimization.

All results of the chapter are proved in the book. Lemmas 5.1 and 5.2 are cited from Cesa-Bianchi and Lugosi (2006). None of them is formalized on Prove2Me. The published mirror-descent bound of Bandit Algorithms XII treats linear losses with a comparator inside DDD and Euclidean-space vectors; it is not Theorem 5.3. The mission adds a machine-checked version of the whole chain, from Legendre duality to the explicit constant q2mdn/(q−1)q\sqrt{2mdn/(q-1)}q2mdn/(q−1)​, with two of the printed statements corrected (below).

Difficulty

The pathwise mirror-descent inequality is a telescoping argument, but several of its steps rest on convex analysis that Mathlib does not package. One is the existence and interior location of Bregman projections onto a set that touches the boundary of DDD. Another is the differentiability of F∗F^*F∗ on the open dual space and the identity ∇F∗=(∇F)−1\nabla F^*=(\nabla F)^{-1}∇F∗=(∇F)−1. A third is the closed form of Fψ∗F_\psi^*Fψ∗​ for a potential defined through an improper integral of ψ−1\psi^{-1}ψ−1.

The probabilistic step is not a martingale argument. Only conditioning on the current iterate xtx_txt​ is available. The estimate (5.5) divides by xt(i)x_t(i)xt​(i), so its integrability and unbiasedness have to be derived from the fact that the iterates stay in the open orthant. Finally, the explicit constant requires a Hölder step, ∑ix1(i)1−1/q≤m(q−1)/qd1/q\sum_ix_1(i)^{1-1/q}\le m^{(q-1)/q}d^{1/q}∑i​x1​(i)1−1/q≤m(q−1)/qd1/q, and the matching bound ∑ixt(i)1/q≤m1/qd1−1/q\sum_ix_t(i)^{1/q}\le m^{1/q}d^{1-1/q}∑i​xt​(i)1/q≤m1/qd1−1/q.

Formalization scope

Vectors are Fin d → ℝ. The arm set is a Set of 0/10/10/1 vectors with coordinate sum mmm, and K\mathcal KK is convexHull ℝ C. Rounds are t=1,…,nt=1,\dots,nt=1,…,n, sums run over Finset.Icc 1 n, and index 000 is unused. A randomized run is a family of measurable processes xt,vt,ℓ~t,wtx_t, v_t, \tilde\ell_t, w_txt​,vt​,ℓ~t​,wt​ on a probability space, with the deterministic OMD recursion holding on every sample path. E[⋅∣xt]\mathbb E[\cdot\mid x_t]E[⋅∣xt​] is the coordinatewise conditional expectation given σ(xt)\sigma(x_t)σ(xt​), which is exactly what the book's proofs use. Losses are oblivious, so Rˉn≤B\bar R_n\le BRˉn​≤B is stated as "for every x∈Kx\in\mathcal Kx∈K, E∑tℓt⊤vt−∑tℓt⊤x≤B\mathbb E\sum_t\ell_t^\top v_t-\sum_t\ell_t^\top x\le BE∑t​ℓt⊤​vt​−∑t​ℓt⊤​x≤B". F∗F^*F∗ is valued in EReal, and DF∗D_{F^*}DF∗​ is evaluated only on the open dual space, where F∗F^*F∗ is finite. Wherever an expectation of a possibly non-integrable quantity appears on a right-hand side, its integrability is assumed: the book's bound is then +∞+\infty+∞ and trivial, while Lean's integral would be 000.

Corrections and instantiations, each labelled in the item's Formalization Note:

  • Theorem 5.7, corrected misprint. The book prints η=2q−1m1−2/qd1−2/q\eta=\sqrt{\frac2{q-1}\frac{m^{1-2/q}}{d^{1-2/q}}}η=q−12​d1−2/qm1−2/q​​. The proof (p. 81) gives the stated bound only for η=2q−1m1−2/qn d1−2/q\eta=\sqrt{\frac2{q-1}\frac{m^{1-2/q}}{n\,d^{1-2/q}}}η=q−12​nd1−2/qm1−2/q​​, which is stated. At q=2q=2q=2 this is η=2/n\eta=\sqrt{2/n}η=2/n​.
  • Theorem 5.5, corrected misprint. In the first bound the book prints E[∥xt−x~t∥ ∥g~t∥∗]\mathbb E[\|x_t-\tilde x_t\|\,\|\tilde g_t\|_*]E[∥xt​−x~t​∥∥g~​t​∥∗​]. That statement fails for ℓt(x)=x2\ell_t(x)=x^2ℓt​(x)=x2 on [−1,1][-1,1][−1,1] with F=x2/2F=x^2/2F=x2/2 and x~t=±1\tilde x_t=\pm1x~t​=±1. The version stated uses ∥∇ℓt(x~t)∥∗\|\nabla\ell_t(\tilde x_t)\|_*∥∇ℓt​(x~t​)∥∗​, as the proof's first inequality does. The linear-loss bound is stated as printed.
  • Lemma 5.2. "For all z∈K∩Dz\in K\cap Dz∈K∩D" is read as "for the projection zzz", which lies in K∩DK\cap DK∩D.
  • Hypotheses made explicit: q>1q>1q>1; non-negativity of the estimates in Theorem 5.6 (used in its proof); unbiasedness E[ℓ~t∣xt]=ℓt\mathbb E[\tilde\ell_t\mid x_t]=\ell_tE[ℓ~t​∣xt​]=ℓt​ in the general parts of Theorems 5.6 and 5.7; K∩(0,∞)d≠∅\mathcal K\cap(0,\infty)^d\ne\emptysetK∩(0,∞)d=∅ (OMD's requirement K∩D≠∅K\cap D\ne\emptysetK∩D=∅); a subgradient selection as an explicit input.
  • Theorem 5.6's particular bound uses the book's η=2mndln⁡dm\eta=\sqrt{\frac{2m}{nd}\ln\frac dm}η=nd2m​lnmd​​ as printed. There are no O(·) constants in the chapter's statements.

A trivializing formalization would let η\etaη, xtx_txt​ or the estimate be junk values: an OSMD step at η=0\eta=0η=0, a Lean division x/0=0x/0=0x/0=0, or a regret written as a real infimum over an unbounded set. Here every run is the book's algorithm on the open orthant, and each bound is stated against every comparator in K\mathcal KK.

Reusable beyond this mission: the Legendre/Bregman layer, the OMD run predicate and the ω\omegaω-potential layer. Proofs of Lemmas 5.1 and 5.2 in this generality would be welcome additions to the library.

Selected references

  • S. Bubeck, N. Cesa-Bianchi, Regret Analysis of Stochastic and Nonstochastic Multi-armed Bandit Problems, Foundations and Trends in Machine Learning 5(1), 2012; arXiv:1204.5721v2. https://arxiv.org/abs/1204.5721
  • N. Cesa-Bianchi, G. Lugosi, Prediction, Learning, and Games, Cambridge University Press, 2006. https://doi.org/10.1017/CBO9780511546921
  • 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
  • J.-Y. Audibert, S. Bubeck, G. Lugosi, Regret in online combinatorial optimization, Mathematics of Operations Research 39(1), 2014. https://doi.org/10.1287/moor.2013.0598
12 thms1 active userReviewed
Operations ResearchOptimizationProbability+1·Captain: mikedeng1

Dimensioning Large Call Centers I: The Rationalized Staffing Function Is Asymptotically OptimalResearch Paper

Motivation

A call center with NNN agents facing Poisson arrivals at rate λ\lambdaλ and exponential service at rate μ\muμ is the M/M/N (Erlang-C) queue. Choosing NNN trades the cost of agents against the cost of customers waiting, and in practice it is done with the square-root safety-staffing rule N≈R+yRN \approx R + y\sqrt RN≈R+yR​, where R=λ/μR = \lambda/\muR=λ/μ is the offered load. Borst, Mandelbaum and Reiman (CWI Report PNA-R0015, 2000; published in Operations Research 52(1), 2004, doi:10.1287/opre.1030.0081) turned that rule of thumb into an optimization result: for a general convex staffing cost and a general waiting-cost function, they identify the safety factor yyy that makes the rule asymptotically optimal as the arrival rate grows.

Timeline of the asymptotic regime the paper builds on:

  • 1917. Erlang's delay formula π(N,ν)\pi(N,\nu)π(N,ν) for the M/M/N queue.
  • 1981. Halfin and Whitt (Oper. Res. 29(3)) show that with N=R+βRN = R + \beta\sqrt RN=R+βR​ servers the probability of waiting converges to a limit P(β)∈(0,1)P(\beta) \in (0,1)P(β)∈(0,1), the quality-and-efficiency-driven regime.
  • 2000/2004. Borst, Mandelbaum and Reiman classify cost structures into a rationalized, an efficiency-driven and a quality-driven regime, and prove asymptotic optimality of an explicit staffing rule in each.

This mission is the first of a series of four on that paper and covers the rationalized regime (Section 5), where staffing and waiting costs are of the same order.

Setting

The service rate μ>0\mu > 0μ>0 is fixed and the arrival rate λ\lambdaλ grows. A staffing cost FFF, defined on (0,∞)(0,\infty)(0,∞), is convex and strictly increasing; it does not depend on λ\lambdaλ. For each λ>0\lambda > 0λ>0 a waiting-cost function DλD_\lambdaDλ​ satisfies Dλ(0)=0D_\lambda(0)=0Dλ​(0)=0, is strictly increasing on [0,∞)[0,\infty)[0,∞), and makes

G(N,λ)=(Nμ−λ)∫0∞Dλ(t) e−(Nμ−λ)t dtG(N,\lambda) = (N\mu-\lambda)\int_0^\infty D_\lambda(t)\,e^{-(N\mu-\lambda)t}\,dtG(N,λ)=(Nμ−λ)∫0∞​Dλ​(t)e−(Nμ−λ)tdt

finite for every N>λ/μN > \lambda/\muN>λ/μ. With the Erlang-C formula

π(N,ν)=νNN!{(1−νN)∑n=0N−1νnn!+νNN!}−1,\pi(N,\nu) = \frac{\nu^N}{N!}\Big\{\big(1-\tfrac{\nu}{N}\big)\sum_{n=0}^{N-1}\frac{\nu^n}{n!}+\frac{\nu^N}{N!}\Big\}^{-1},π(N,ν)=N!νN​{(1−Nν​)n=0∑N−1​n!νn​+N!νN​}−1,

the expected total cost of staffing N>λ/μN > \lambda/\muN>λ/μ agents is C(N,λ)=F(N)+λ π(N,λ/μ) G(N,λ)C(N,\lambda) = F(N) + \lambda\,\pi(N,\lambda/\mu)\,G(N,\lambda)C(N,λ)=F(N)+λπ(N,λ/μ)G(N,λ), and Nλ∗N^*_\lambdaNλ∗​ is any integer N>λ/μN > \lambda/\muN>λ/μ minimizing it (7).

In normalized units Nλ(x)=λ/μ+xλ/μN_\lambda(x) = \lambda/\mu + x\sqrt{\lambda/\mu}Nλ​(x)=λ/μ+xλ/μ​ the paper defines Fλ(x)=F(Nλ(x))−F(λ/μ)F_\lambda(x) = F(N_\lambda(x)) - F(\lambda/\mu)Fλ​(x)=F(Nλ​(x))−F(λ/μ), Gλ(x)=λG(Nλ(x),λ)G_\lambda(x) = \lambda G(N_\lambda(x),\lambda)Gλ​(x)=λG(Nλ​(x),λ), the continuous delay probability πλ(x)=H(Nλ(x),λ/μ)\pi_\lambda(x) = H(N_\lambda(x),\lambda/\mu)πλ​(x)=H(Nλ​(x),λ/μ) with

H(M,α)={α∫0∞e−αt t (1+t)M−1 dt}−1,H(M,\alpha) = \Big\{\alpha\int_0^\infty e^{-\alpha t}\,t\,(1+t)^{M-1}\,dt\Big\}^{-1},H(M,α)={α∫0∞​e−αtt(1+t)M−1dt}−1,

and Cλ(x)=Fλ(x)+πλ(x)Gλ(x)C_\lambda(x) = F_\lambda(x) + \pi_\lambda(x)G_\lambda(x)Cλ​(x)=Fλ​(x)+πλ​(x)Gλ​(x), minimized at xλ∗x^*_\lambdaxλ∗​ (8). A surrogate C[z;F^,π^,G^]=F^(z)+π^(z)G^(z)C[z;\hat F,\hat\pi,\hat G] = \hat F(z)+\hat\pi(z)\hat G(z)C[z;F^,π^,G^]=F^(z)+π^(z)G^(z) approximates it. Rounding is measured by

Sλ(x)=min⁡{C(⌊Nλ(x)⌋,λ), C(⌈Nλ(x)⌉,λ)}.(10)S_\lambda(x) = \min\{C(\lfloor N_\lambda(x)\rfloor,\lambda),\,C(\lceil N_\lambda(x)\rceil,\lambda)\}. \tag{10}Sλ​(x)=min{C(⌊Nλ​(x)⌋,λ),C(⌈Nλ​(x)⌉,λ)}.(10)

The Halfin–Whitt delay function is P(x)=(1+x/h(−x))−1P(x) = \big(1 + x/h(-x)\big)^{-1}P(x)=(1+x/h(−x))−1, with h=ϕ/(1−Φ)h = \phi/(1-\Phi)h=ϕ/(1−Φ) the standard normal hazard rate (11). Asymptotic equality aλ≈∞bλa_\lambda \stackrel{\infty}{\approx} b_\lambdaaλ​≈∞bλ​ means aλ/bλ→1a_\lambda/b_\lambda \to 1aλ​/bλ​→1 as λ→∞\lambda\to\inftyλ→∞.

Formalization targets

Goal: Theorem 5.1

Assume the rationalized condition (18): for some κ>0\kappa > 0κ>0, Fλ(κ)/Gλ(κ)→γ∈(0,∞)F_\lambda(\kappa)/G_\lambda(\kappa) \to \gamma \in (0,\infty)Fλ​(κ)/Gλ​(κ)→γ∈(0,∞). Let yλ∗y^*_\lambdayλ∗​ minimize Fλ(y)+P(y)Gλ(y)F_\lambda(y) + P(y)G_\lambda(y)Fλ​(y)+P(y)Gλ​(y) over y>0y>0y>0 (19). Then

lim⁡λ→∞Sλ(yλ∗)−F(λ/μ)C(Nλ∗,λ)−F(λ/μ)=1.\lim_{\lambda\to\infty}\frac{S_\lambda(y^*_\lambda) - F(\lambda/\mu)}{C(N^*_\lambda,\lambda) - F(\lambda/\mu)} = 1.λ→∞lim​C(Nλ∗​,λ)−F(λ/μ)Sλ​(yλ∗​)−F(λ/μ)​=1.

The goal fixes no constant and no rate: it asserts only that the excess cost of the explicit rule is asymptotically the optimal excess cost.

Milestones

  • Lemma C.1: GλG_\lambdaGλ​ is strictly convex and strictly decreasing on (0,∞)(0,\infty)(0,∞).
  • Section 3, p. 12: H(N,ν)=π(N,ν)H(N,\nu) = \pi(N,\nu)H(N,ν)=π(N,ν) at integers N>ν>0N > \nu > 0N>ν>0.
  • Lemma 3.1, Lemma 3.2, Corollary 3.3: the approximation principle. If the surrogate approximates CλC_\lambdaCλ​ at both xλ∗x^*_\lambdaxλ∗​ and its own minimizer zλ∗z^*_\lambdazλ∗​, then rounding Nλ(zλ∗)N_\lambda(z^*_\lambda)Nλ​(zλ∗​) is asymptotically optimal.
  • Eqs. (13)–(14): FλF_\lambdaFλ​ preserves lim sup⁡\limsuplimsup-separation of ratios.
  • Lemma 4.1 (Halfin & Whitt): for bounded xλx_\lambdaxλ​, πλ(xλ)/P(xλ)→1\pi_\lambda(x_\lambda)/P(x_\lambda) \to 1πλ​(xλ​)/P(xλ​)→1.

Significance

The theorem justifies the square-root staffing rule from first principles for a broad cost class. In Example 5.3 of the paper (linear staffing cost ccc per agent, linear waiting cost aaa per unit time) it gives N∗≈R+y∗(a/c)RN^* \approx R + y^*(a/c)\sqrt RN∗≈R+y∗(a/c)R​, with y∗(r)y^*(r)y∗(r) the minimizer of y+rP(y)/yy + rP(y)/yy+rP(y)/y, a one-dimensional rule computable once for all loads. Corollary 3.3 is reused verbatim by the efficiency-driven and quality-driven theorems of the paper (missions II and III of this series), and Lemma 4.1 is the analytic input of all three.

The result has been proved since 2000; no machine-checked proof of it, or of the Halfin–Whitt limit for the continuous extension πλ\pi_\lambdaπλ​, is known to exist. The mission produces a formal proof of the regime theorem together with reusable formal statements of the Erlang-C function, its integral representation, and the Halfin–Whitt limit.

Difficulty

The reduction from discrete to continuous staffing (Lemmas 3.1–3.2) is elementary once unimodality of CλC_\lambdaCλ​ is available, but unimodality rests on convexity of πλ\pi_\lambdaπλ​, which the paper cites rather than proves, and on Lemma C.1, which needs differentiation under an improper integral. The central difficulty is Lemma 4.1: the paper derives it from Halfin and Whitt's limit theorem, which is stated for integer server counts, while πλ\pi_\lambdaπλ​ is evaluated at non-integer Nλ(xλ)N_\lambda(x_\lambda)Nλ​(xλ​); a proof needs a uniform Laplace-type asymptotic for the integral defining HHH. A further obstacle is bounding xλ∗x^*_\lambdaxλ∗​: the obvious route through continuity of the optimizer fails because nothing converges, and the paper instead argues by contradiction via (14).

Formalization scope

All objects live in DimCallCenters.Rationalized. The arrival rate is a real lam, and every limit is Filter.atTop on R\mathbb RR with μ\muμ fixed. The queue itself is not modelled; the paper's theorems are statements about the closed-form cost C(N,λ)C(N,\lambda)C(N,λ), and so are these. Committed conventions:

  1. The standing assumptions are a structure WaitModel (μ>0\mu>0μ>0; Dλ(0)=0D_\lambda(0)=0Dλ​(0)=0; DλD_\lambdaDλ​ strictly increasing on [0,∞)[0,\infty)[0,∞); t↦Dλ(t)e−θtt\mapsto D_\lambda(t)e^{-\theta t}t↦Dλ​(t)e−θt integrable on (0,∞)(0,\infty)(0,∞) for every θ>0\theta>0θ>0, which is the paper's finiteness of GGG). FFF is convex and strictly increasing on (0,∞)(0,\infty)(0,∞).
  2. Staffing levels in C(N,λ)C(N,\lambda)C(N,λ) are natural numbers; GGG and HHH take real NNN.
  3. Argmins (Nλ∗N^*_\lambdaNλ∗​, xλ∗x^*_\lambdaxλ∗​, zλ∗z^*_\lambdazλ∗​, yλ∗y^*_\lambdayλ∗​) are hypotheses that a given function is a minimizer, for every λ>0\lambda>0λ>0; ties are allowed and the theorems hold for every choice.
  4. In SλS_\lambdaSλ​ the floor term is omitted when ⌊Nλ(x)⌋≤λ/μ\lfloor N_\lambda(x)\rfloor \le \lambda/\mu⌊Nλ​(x)⌋≤λ/μ, where CCC is undefined.
  5. lim sup⁡\limsuplimsup and lim inf⁡\liminfliminf relations are written with ∃ᶠ/∀ᶠ, not Filter.limsup on R\mathbb RR.
  6. Added hypothesis. The goal assumes G(N,λ)→∞G(N,\lambda)\to\inftyG(N,λ)→∞ as N↓λ/μN\downarrow\lambda/\muN↓λ/μ. The paper asserts this limit on p. 12, but it does not follow from its assumptions (it fails for bounded DλD_\lambdaDλ​); it is equivalent to DλD_\lambdaDλ​ being unbounded and is what makes the continuous optimum exist.

The hypotheses are met by linear staffing and waiting costs (F(N)=cNF(N)=cNF(N)=cN, Dλ(t)=atD_\lambda(t)=atDλ​(t)=at), for which (18) holds with γ=cκ2/a\gamma = c\kappa^2/aγ=cκ2/a, so the goal is not vacuous. It is not trivialized by junk values either: the ratio's denominator is positive at every λ>0\lambda>0λ>0, and SλS_\lambdaSλ​ never evaluates CCC at an unstable level.

Needed infrastructure: Laplace asymptotics for ∫0∞e−αtt(1+t)M−1dt\int_0^\infty e^{-\alpha t}t(1+t)^{M-1}dt∫0∞​e−αtt(1+t)M−1dt, differentiation under the integral sign for GGG, and convexity of πλ\pi_\lambdaπλ​. All of these are reusable for missions II–IV. Proofs of the milestones in any order are welcome, as are proofs of the convexity facts the paper cites from its references [9], [10].

Selected references

  • S. Borst, A. Mandelbaum, M. I. Reiman, Dimensioning Large Call Centers, CWI Report PNA-R0015, 2000; Operations Research 52(1):17–34, 2004. https://doi.org/10.1287/opre.1030.0081
  • S. Halfin, W. Whitt, Heavy-Traffic Limits for Queues with Many Exponential Servers, Operations Research 29(3):567–588, 1981. https://doi.org/10.1287/opre.29.3.567
  • A. K. Erlang, Solution of some problems in the theory of probabilities of significance in automatic telephone exchanges, Elektroteknikeren 13, 1917.
21 thms2 active usersReviewed
Bandit AlgorithmsConvex OptimizationMachine Learning+2·Captain: mikedeng1

Regret Analysis of Stochastic and Nonstochastic Multi-armed Bandit Problems V: Bandit Convex Optimization with One-Point FeedbackTextbook

Motivation

In bandit convex optimization a forecaster repeatedly picks a point xtx_txt​ of a convex set K⊆Rd\mathcal K\subseteq\mathbb R^dK⊆Rd, and an adversary picks a convex loss ℓt\ell_tℓt​. The forecaster pays ℓt(xt)\ell_t(x_t)ℓt​(xt​) and observes only that number: it never sees the function, its gradient, or its value elsewhere. This is the model of online optimization with only function-value access, as in tuning a system online from measured costs, dynamic pricing with an unknown convex demand-cost curve, or routing with path costs observed only on the route taken. The question is how fast the forecaster can approach the best fixed point in hindsight.

Chapter 6 of Bubeck and Cesa-Bianchi's monograph (arXiv:1204.5721v2, Foundations and Trends in Machine Learning 5(1), 2012) treats the problem through spherical gradient estimates fed to projected gradient descent. The one-point method is due to Flaxman, Kalai and McMahan (SODA 2005, arXiv:cs/0408007), who obtained an O(n3/4)\mathcal O(n^{3/4})O(n3/4) regret bound. Agarwal, Dekel and Xiao (COLT 2010) showed that two function evaluations per round allow O(n)\mathcal O(\sqrt n)O(n​). Whether one-point feedback admits n\sqrt nn​ regret was open when the monograph was written (p. 94); Bubeck, Eldan and Lee (STOC 2017, arXiv:1607.03084) later obtained n\sqrt nn​ regret up to logarithmic and polynomial-in-ddd factors for convex losses, with a different and much more involved algorithm.

Setting

Let B={x∈Rd:∥x∥≤1}\mathbb B=\{x\in\mathbb R^d:\|x\|\le1\}B={x∈Rd:∥x∥≤1} be the closed Euclidean unit ball and S={x:∥x∥=1}\mathbb S=\{x:\|x\|=1\}S={x:∥x∥=1} the unit sphere, with unnormalized spherical measure σ\sigmaσ, so that σ(S)=d Vol(B)\sigma(\mathbb S)=d\,\mathrm{Vol}(\mathbb B)σ(S)=dVol(B). Fix δ>0\delta>0δ>0. For a loss ℓ\ellℓ, the smoothed loss is ℓ~(x)=E ℓ(x+δB)\widetilde\ell(x)=\mathbb E\,\ell(x+\delta B)ℓ(x)=Eℓ(x+δB) with BBB uniform on B\mathbb BB.

The set K\mathcal KK is closed and convex with rB⊆K⊆RBr\mathbb B\subseteq\mathcal K\subseteq R\mathbb BrB⊆K⊆RB. The losses ℓ1,ℓ2,⋯:Rd→R\ell_1,\ell_2,\dots:\mathbb R^d\to\mathbb Rℓ1​,ℓ2​,⋯:Rd→R are GGG-Lipschitz, differentiable and convex, and are fixed before the game (an oblivious adversary).

OSGD (Online Stochastic Gradient Descent) on a set K′\mathcal K'K′ with learning rate η\etaη starts at x1=0x_1=0x1​=0 and sets xt+1=argmin⁡y∈K′∥y−(xt−ηg~t(xt))∥x_{t+1}=\operatorname{argmin}_{y\in\mathcal K'}\|y-(x_t-\eta\widetilde g_t(x_t))\|xt+1​=argminy∈K′​∥y−(xt​−ηg​t​(xt​))∥, where g~t\widetilde g_tg​t​ is a gradient estimate. With S1,S2,…S_1,S_2,\dotsS1​,S2​,… independent and uniform on S\mathbb SS:

  • the two-point estimate (6.1) is g~t(xt)=d2δ(ℓt(Xt+)−ℓt(Xt−))St\widetilde g_t(x_t)=\frac d{2\delta}\big(\ell_t(X_t^+)-\ell_t(X_t^-)\big)S_tg​t​(xt​)=2δd​(ℓt​(Xt+​)−ℓt​(Xt−​))St​ with Xt±=xt±δStX_t^\pm=x_t\pm\delta S_tXt±​=xt​±δSt​; the played point is Xt+X_t^+Xt+​ or Xt−X_t^-Xt−​ by a fair coin;
  • the one-point estimate (6.3) is g~t(xt)=dδ ℓt(X~t)St\widetilde g_t(x_t)=\frac d\delta\,\ell_t(\widetilde X_t)S_tg​t​(xt​)=δd​ℓt​(Xt​)St​ with played point X~t=xt+δSt\widetilde X_t=x_t+\delta S_tXt​=xt​+δSt​.

OSGD runs on the shrunken set K′=(1−δ/r)K\mathcal K'=(1-\delta/r)\mathcal KK′=(1−δ/r)K, so that the perturbed points stay in K\mathcal KK. The pseudo-regret is

R‾n=E∑t=1nℓt(X~t)−min⁡x∈K∑t=1nℓt(x).\overline R_n=\mathbb E\sum_{t=1}^n\ell_t(\widetilde X_t)-\min_{x\in\mathcal K}\sum_{t=1}^n\ell_t(x).Rn​=Et=1∑n​ℓt​(Xt​)−x∈Kmin​t=1∑n​ℓt​(x).

Formalization targets

Goal: Theorem 6.2, tuned

If in addition ∣ℓt∣≤L|\ell_t|\le L∣ℓt​∣≤L on K\mathcal KK, and δ=(2n)−1/4RdL/((3+R/r)G)\delta=(2n)^{-1/4}\sqrt{RdL/((3+R/r)G)}δ=(2n)−1/4RdL/((3+R/r)G)​, η=(2n)−3/4R3/(dL(3+R/r)G)\eta=(2n)^{-3/4}\sqrt{R^3/(dL(3+R/r)G)}η=(2n)−3/4R3/(dL(3+R/r)G)​, then one-point OSGD satisfies

R‾n≤4n3/4RdL (3+R/r) G.\overline R_n\le 4n^{3/4}\sqrt{RdL\,(3+R/r)\,G}.Rn​≤4n3/4RdL(3+R/r)G​.

Milestones

  1. Lemma 6.1: ∇∫Bℓ(x+δb) db=1δ∫Sℓ(x+δs)s dσ(s)\nabla\int_{\mathbb B}\ell(x+\delta b)\,db=\frac1\delta\int_{\mathbb S}\ell(x+\delta s)s\,d\sigma(s)∇∫B​ℓ(x+δb)db=δ1​∫S​ℓ(x+δs)sdσ(s).
  2. Lemma 6.2: dδE[ℓ(x+δS)S]=∇E ℓ(x+δB)\frac d\delta\mathbb E[\ell(x+\delta S)S]=\nabla\mathbb E\,\ell(x+\delta B)δd​E[ℓ(x+δS)S]=∇Eℓ(x+δB).
  3. Eq. (6.2): ∣ℓ(x)−ℓ~(x)∣≤δG|\ell(x)-\widetilde\ell(x)|\le\delta G∣ℓ(x)−ℓ(x)∣≤δG.
  4. Lemma 6.3: the queried points' regret against xxx is at most the smoothed regret of the iterates against (1−ξ)x(1-\xi)x(1−ξ)x, plus 3δGn+ξGRn3\delta Gn+\xi GRn3δGn+ξGRn.
  5. Theorem 6.1: two-point OSGD has R‾n≤R2/η+η(Gd)2n+δ(3+R/r)Gn\overline R_n\le R^2/\eta+\eta(Gd)^2n+\delta(3+R/r)GnRn​≤R2/η+η(Gd)2n+δ(3+R/r)Gn, and R‾n≤2RGdn+δ(3+R/r)Gn\overline R_n\le 2RGd\sqrt n+\delta(3+R/r)GnRn​≤2RGdn​+δ(3+R/r)Gn for η=R/(Gdn)\eta=R/(Gd\sqrt n)η=R/(Gdn​).
  6. Theorem 6.2, first display: one-point OSGD has R‾n≤R2/η+(dL)2δ2ηn+δ(3+R/r)Gn\overline R_n\le R^2/\eta+\frac{(dL)^2}{\delta^2}\eta n+\delta(3+R/r)GnRn​≤R2/η+δ2(dL)2​ηn+δ(3+R/r)Gn for every 0<δ≤r0<\delta\le r0<δ≤r and η>0\eta>0η>0.

Significance

The n3/4n^{3/4}n3/4 bound shows that a single function value per round suffices for sublinear regret against any oblivious sequence of Lipschitz convex losses, with a forecaster whose only operations are a random perturbation and a Euclidean projection. The smoothing identity of Lemmas 6.1–6.2 is the basic tool of zeroth-order (derivative-free) optimization, used well beyond bandits, and Theorem 6.1 is the n\sqrt nn​ benchmark for two-point methods.

All results are proved in the source. To the best of current knowledge none is formalized: the related items of the Introduction to Online Convex Optimization series on Prove2Me (Hazan's Lemma 6.7 and Theorem 6.9) were formalized with missing hypotheses and are recorded as disproved. This mission produces machine-checked statements with every hypothesis explicit, and the formal infrastructure (sphere measure calculus, a projected stochastic gradient analysis) for later zeroth-order results.

Difficulty

Two steps resist a direct formal treatment. First, Lemma 6.1 is a divergence-theorem identity on the ball; Mathlib has the sphere measure and polar coordinates, but its divergence theorem covers boxes rather than balls, so differentiating the ball average in xxx requires either such a theorem or a direct argument about translates of the ball. Second, the regret analysis takes expectations of quantities that depend on the whole past: the iterate xtx_txt​ is a function of S1,…,St−1S_1,\dots,S_{t-1}S1​,…,St−1​, and unbiasedness E[g~t∣xt]=∇ℓ~t(xt)\mathbb E[\widetilde g_t\mid x_t]=\nabla\widetilde\ell_t(x_t)E[g​t​∣xt​]=∇ℓt​(xt​) holds only conditionally, via independence of StS_tSt​ from the past. A pathwise gradient-descent inequality must be combined with this conditional expectation round by round, with measurability of the projected iterates established along the way. The naive approach of treating the estimate as the true gradient of ℓt\ell_tℓt​ fails: it is a gradient of ℓ~t\widetilde\ell_tℓt​, and the gap is handled only by Eq. (6.2) and Lemma 6.3.

Formalization scope

Points are in EuclideanSpace ℝ (Fin d) with d≥1d\ge1d≥1; rounds are t=1,2,…t=1,2,\dotst=1,2,…, sums run over Finset.Icc 1 n. σ\sigmaσ is Mathlib's Measure.toSphere of Lebesgue measure; the uniform laws are normalized restrictions. Randomness lives on an arbitrary probability space; the directions StS_tSt​ are measurable, mutually independent (iIndepFun) and uniform on S\mathbb SS, and in Theorem 6.1 the pairs (St,Ct)(S_t,C_t)(St​,Ct​) are independent with CtC_tCt​ a fair sign independent of StS_tSt​. A run of OSGD is a predicate (start at 000, each iterate a Euclidean projection onto (1−δ/r)K(1-\delta/r)\mathcal K(1−δ/r)K), which determines the run uniquely, so the forecaster uses only observed values and its own randomness. The losses are Lipschitz, differentiable and convex on all of Rd\mathbb R^dRd; the bound ∣ℓt∣≤L|\ell_t|\le L∣ℓt​∣≤L is on K\mathcal KK, because a convex function bounded on Rd\mathbb R^dRd is constant. The minimum over K\mathcal KK is an infimum over the subtype K\mathcal KK, attained in every theorem.

Conventions and corrections, each stated in the item's Formalization Note:

  • Lemma 6.1 carries the factor 1/δ1/\delta1/δ that the printed statement omits and the proof contains (corrected misprint).
  • Theorem 6.1's second display prints η=R/(GDn)\eta=R/(GD\sqrt n)η=R/(GDn​) and a limit "for δ→0\delta\to0δ→0"; the item states R‾n≤2RGdn+δ(3+R/r)Gn\overline R_n\le 2RGd\sqrt n+\delta(3+R/r)GnRn​≤2RGdn​+δ(3+R/r)Gn for η=R/(Gdn)\eta=R/(Gd\sqrt n)η=R/(Gdn​) and every admissible δ\deltaδ, which implies the limit (corrected misprint).
  • Theorems 6.1 and 6.2 add 0<δ≤r0<\delta\le r0<δ≤r, which the proofs need for Xt±,X~t∈KX_t^\pm,\widetilde X_t\in\mathcal KXt±​,Xt​∈K; for the tuned δ\deltaδ of the goal it is a condition on nnn.
  • The goal adds G,L>0G,L>0G,L>0 and n≥1n\ge1n≥1, which its formulas for δ,η\delta,\etaδ,η need; the constant 444 is the book's rounding of 2⋅23/42\cdot2^{3/4}2⋅23/4 and is kept, as is the form R2/ηR^2/\etaR2/η.

The statements cannot be satisfied trivially: the run is pinned by its recursion, the losses are fixed before the randomness, the expectations are of bounded measurable functions (no zero-valued Bochner integrals), and the minimum is over the nonempty compact K\mathcal KK. Section 6.3 (Lemma 6.4, Theorem 6.3) is not included, because its algorithm box and proof use different stage lengths and its unimodality condition is stated on a smaller set than the proof uses.

Needed infrastructure: calculus of ball averages and sphere integrals, symmetry of the uniform sphere law, nonexpansiveness of projections onto closed convex sets, and conditional-expectation bookkeeping for adapted iterates. Each is reusable for zeroth-order optimization; contributions of any of them as separate lemmas are welcome.

Selected references

  • S. Bubeck, N. Cesa-Bianchi, Regret Analysis of Stochastic and Nonstochastic Multi-armed Bandit Problems, Foundations and Trends in Machine Learning 5(1), 2012. arXiv:1204.5721v2, doi:10.1561/2200000024
  • A. Flaxman, A. Kalai, H. B. McMahan, Online convex optimization in the bandit setting: gradient descent without a gradient, SODA 2005. arXiv:cs/0408007
  • A. Agarwal, O. Dekel, L. Xiao, Optimal algorithms for online convex optimization with multi-point bandit feedback, COLT 2010. link
  • S. Bubeck, R. Eldan, Y. T. Lee, Kernel-based methods for bandit convex optimization, STOC 2017. arXiv:1607.03084
10 thms2 active usersReviewed
Convex OptimizationOperations ResearchOptimization·Captain: mikedeng1

Project Scheduling with Time Windows and Scarce Resources VII: A Locally Quasiconcave Objective Always Has a Quasistable Optimal ScheduleTextbook

Motivation

Resource-constrained project scheduling asks for start times of the activities of a project that respect precedence-type time lags and the capacities of renewable resources (machines, crews, equipment). Classical project scheduling minimizes the project duration, a regular objective: delaying an activity never helps. Many objectives met in practice are not regular. The resource investment problem minimizes the cost of the resource capacities that must be procured; resource levelling problems minimize fluctuations of resource usage over time; the resource renting problem trades fixed procurement against time-dependent renting costs; net present value and earliness–tardiness objectives reward late as well as early starts. For such objectives the familiar fact that "some active schedule is optimal" fails, and algorithms need another finite set of candidate schedules that is guaranteed to contain an optimum.

Chapter 3 of Neumann, Schwindt and Zimmermann, Project Scheduling with Time Windows and Scarce Resources (2nd ed., Springer 2003, doi:10.1007/978-3-540-24800-2), organizes the objective functions of project scheduling into seven classes and pairs each class with a class of schedules that contains an optimal schedule. This mission formalizes §3.3 of that chapter. The classification goes back to Neumann, Nübel and Schwindt (2000) and Zimmermann (2001); the two locally defined classes, and the matching schedule classes of quasiactive and quasistable schedules, are the book's device for covering discontinuous resource-based objectives.

Setting

A project consists of activities V={0,1,…,n+1}V=\{0,1,\dots,n+1\}V={0,1,…,n+1}, n≥1n\ge 1n≥1, where 000 and n+1n+1n+1 are fictitious activities marking the project beginning and completion. Activity iii has an integer duration pip_ipi​ (p0=pn+1=0p_0=p_{n+1}=0p0​=pn+1​=0, pi>0p_i>0pi​>0 otherwise). The project network has an arc set EEE with integer weights δij\delta_{ij}δij​; a schedule is a vector S=(S0,…,Sn+1)S=(S_0,\dots,S_{n+1})S=(S0​,…,Sn+1​) of real start times with S0=0S_0=0S0​=0, S≥0S\ge 0S≥0, and it is time-feasible if Sj−Si≥δijS_j-S_i\ge\delta_{ij}Sj​−Si​≥δij​ for all ⟨i,j⟩∈E\langle i,j\rangle\in E⟨i,j⟩∈E. A maximum project duration dˉ∈N\bar d\in\mathbb Ndˉ∈N is prescribed through a backward arc ⟨n+1,0⟩\langle n+1,0\rangle⟨n+1,0⟩ of weight −dˉ-\bar d−dˉ, so Sn+1≤dˉS_{n+1}\le\bar dSn+1​≤dˉ. Each renewable resource kkk has capacity RkR_kRk​, activity iii uses rikr_{ik}rik​ units while in progress, and rk(S,t)r_k(S,t)rk​(S,t) is the total usage at time ttt. The feasible region S\mathcal SS consists of the time-feasible schedules with rk(S,t)≤Rkr_k(S,t)\le R_krk​(S,t)≤Rk​ for all kkk and ttt.

For an objective function f:R≥0n+2→Rf:\mathbb R^{n+2}_{\ge 0}\to\mathbb Rf:R≥0n+2​→R, problem PS∣temp,dˉ∣fPS|temp,\bar d|fPS∣temp,dˉ∣f asks for an optimal schedule: some S∈SS\in\mathcal SS∈S with f(S)≤f(S′)f(S)\le f(S')f(S)≤f(S′) for all S′∈SS'\in\mathcal SS′∈S.

A schedule induces the strict order O(S)={(i,j)∣i≠j, Sj≥Si+pi}O(S)=\{(i,j)\mid i\ne j,\ S_j\ge S_i+p_i\}O(S)={(i,j)∣i=j, Sj​≥Si​+pi​} of precedences it realizes. The equal-order set of SSS is

ST=(O(S))={S′ time-feasible∣Sj′≥Si′+pi ∀(i,j)∈O(S), O(S′)=O(S)},\mathcal S_T^{=}(O(S))=\{S'\text{ time-feasible}\mid S'_j\ge S'_i+p_i\ \forall (i,j)\in O(S),\ O(S')=O(S)\},ST=​(O(S))={S′ time-feasible∣Sj′​≥Si′​+pi​ ∀(i,j)∈O(S), O(S′)=O(S)},

a polytope with part of its boundary removed. The distinct equal-order sets partition S\mathcal SS into finitely many pieces.

Schedule classes are defined through shifts. A shift from a feasible SSS to a feasible S′≠SS'\ne SS′=S is order-preserving if O(S)⊆O(S′)O(S)\subseteq O(S')O(S)⊆O(S′); it is a left-shift if S′≤SS'\le SS′≤S. Two shifts from SSS to S′S'S′ and S′′S''S′′ are opposite if S′′−S=λ(S′−S)S''-S=\lambda(S'-S)S′′−S=λ(S′−S) with λ<0\lambda<0λ<0. A feasible schedule is active if no feasible left-shift exists, quasiactive if no order-preserving left-shift exists, stable if no pair of opposite shifts to feasible schedules exists, and quasistable if no pair of opposite order-preserving shifts exists.

Objective classes: fff is regular if S≤S′S\le S'S≤S′ implies f(S)≤f(S′)f(S)\le f(S')f(S)≤f(S′); quasiconcave on a set MMM if f(λS+(1−λ)S′)≥min⁡[f(S),f(S′)]f(\lambda S+(1-\lambda)S')\ge\min[f(S),f(S')]f(λS+(1−λ)S′)≥min[f(S),f(S′)] for S,S′∈MS,S'\in MS,S′∈M, λ∈[0,1]\lambda\in[0,1]λ∈[0,1]; lower semicontinuous if f(S)≤lim inf⁡S′→Sf(S′)f(S)\le\liminf_{S'\to S}f(S')f(S)≤liminfS′→S​f(S′) on R≥0n+2\mathbb R^{n+2}_{\ge 0}R≥0n+2​. Then fff is locally regular (class 6) if it is lower semicontinuous and regular on every equal-order set ST=(O(S))\mathcal S_T^{=}(O(S))ST=​(O(S)), S∈SS\in\mathcal SS∈S, and locally quasiconcave (class 7) if it is lower semicontinuous and quasiconcave on every such set.

Formalization targets

Goal: Theorem 3.3.13

For every locally quasiconcave fff,

S≠∅ ⟹ ∃ S quasistable with f(S)=min⁡S′∈Sf(S′).\mathcal S\ne\emptyset\ \Longrightarrow\ \exists\,S\ \text{quasistable with}\ f(S)=\min_{S'\in\mathcal S}f(S').S=∅ ⟹ ∃S quasistable with f(S)=S′∈Smin​f(S′).

Milestones

  • Class 1 (§3.3.2): every regular fff has an active optimal schedule when S≠∅\mathcal S\ne\emptysetS=∅.
  • Class 5 (§3.3.6): every quasiconcave fff has a stable optimal schedule when S≠∅\mathcal S\ne\emptysetS=∅.
  • Eq. (3.3.11): the equal-order sets form a finite partition of S\mathcal SS.
  • Propositions 3.3.5 and 3.3.6: the resource investment objective ∑kckmax⁡trk(S,t)\sum_k c_k\max_t r_k(S,t)∑k​ck​maxt​rk​(S,t) with ck≥0c_k\ge 0ck​≥0 is constant on each equal-order set and lower semicontinuous, hence locally regular.
  • Theorem 3.3.9: every locally regular fff has a quasiactive optimal schedule when S≠∅\mathcal S\ne\emptysetS=∅.

Significance

Quasiactive and quasistable schedules are finite in number: they are the minimal points and the vertices of the finitely many schedule polytopes. Theorem 3.3.13 therefore turns the minimization of any locally quasiconcave objective over a disconnected, non-convex feasible region into a finite search. Class 7 contains the resource levelling objectives ∑ck∑rkt2\sum c_k\sum r_{kt}^2∑ck​∑rkt2​ and ∑ck∑okt\sum c_k\sum o_{kt}∑ck​∑okt​, the total variation of the resource profiles, and the resource renting objective (Propositions 3.3.10 and 3.3.12, and Nübel 2001). The enumeration schemes and decision sets of §3.5–3.7 rest on this result, and Theorem 3.3.9 plays the same role for class 6 (resource investment, changeover times).

The results are proved in the book and the cited papers. As far as a search of the platform shows, none of them, and none of the schedule classes, has a machine-checked formalization; Mathlib supplies lower semicontinuity and quasiconcavity but nothing about schedules. The mission produces a checked version of the classification theorems in the book's exact generality: general time lags (cycles in the network allowed), real start times, and arbitrary objectives given only by their class.

Difficulty

The optimum need not exist a priori: objectives of classes 6 and 7 are discontinuous, and the feasible region is a finite union of polytopes that is in general disconnected. Existence of a minimizer needs compactness of S\mathcal SS (which depends on the deadline arc and the network's path structure) together with lower semicontinuity.

The main obstacle is that the objective is only controlled piecewise. Quasiconcavity holds on each equal-order set separately, and an equal-order set is not closed: a schedule polytope ST(O(S))\mathcal S_T(O(S))ST​(O(S)) also contains schedules inducing strictly larger orders, where the hypothesis on fff says nothing about its relation to the values on ST=(O(S))\mathcal S_T^{=}(O(S))ST=​(O(S)). The obvious argument, taking an optimal schedule and invoking quasiconcavity along the segment of a pair of opposite order-preserving shifts, only relates fff at points of one equal-order set, and it does not by itself produce a schedule that admits no such pair at all. The same issue arises for Theorem 3.3.9 with order-preserving left-shifts, which may cross from one equal-order set into another.

Formalization scope

Activities are Fin (n + 2), with 0 and Fin.last (n + 1) fictitious. Start times are real; objective functions are total functions (Fin (n + 2) → ℝ) → ℝ whose regularity, quasiconcavity and lower semicontinuity are required only on the nonnegative orthant (lower semicontinuity is Mathlib's LowerSemicontinuousOn on the orthant). The deadline Sn+1≤dˉS_{n+1}\le\bar dSn+1​≤dˉ is the network's backward arc, as in §3.1. The project structure records the book's standing property (p. 8) that from each node iii there is a path to n+1n+1n+1 of length at least pip_ipi​; this bounds every activity by dˉ\bar ddˉ. The resource constraints are imposed for all t≥0t\ge 0t≥0, which under that property is the book's 0≤t≤dˉ0\le t\le\bar d0≤t≤dˉ. The peak max⁡trk(S,t)\max_t r_k(S,t)maxt​rk​(S,t) in the resource investment objective is a supremum in N\mathbb NN over t≥0t\ge 0t≥0 of a nonempty finite set, hence attained.

"Optimal" always means minimizing fff over the whole feasible region S\mathcal SS, and the theorems quantify over every function in the class; a formalization with a fixed objective, or with optimality over a single polytope or a single equal-order set, would be a different and weaker statement. The schedule classes are defined through shifts, never as minimal or extreme points, so no statement is true by definition. The only hypothesis besides the class of fff is S≠∅\mathcal S\ne\emptysetS=∅.

The mission restates locally the project model, the induced orders and the shift classes also drafted by the companion missions on schedule classes of this series. Useful contributions beyond the milestones: compactness of S\mathcal SS and closedness of the schedule polytopes, the representation of S\mathcal SS as a finite union of feasible order polytopes, and the finiteness of the sets of quasiactive and quasistable schedules.

Selected references

  • K. Neumann, C. Schwindt, J. Zimmermann, Project Scheduling with Time Windows and Scarce Resources, 2nd ed., Springer, 2003, §3.3. doi:10.1007/978-3-540-24800-2
  • K. Neumann, H. Nübel, C. Schwindt, Active and stable project scheduling, Mathematical Methods of Operations Research 52 (2000), cited in the book as Neumann et al. (2000).
  • J. Zimmermann, Ablauforientiertes Projektmanagement: Modelle, Verfahren und Anwendungen, Gabler, 2001.
11 thms2 active usersReviewed
Bandit AlgorithmsMachine LearningOperations Research+1·Captain: mikedeng1

Regret Analysis of Stochastic and Nonstochastic Multi-armed Bandit Problems III: Contextual Bandits and the Banditron Mistake BoundTextbook

Motivation

In many sequential decision problems the learner sees side information before acting. A news site chooses an article for a visitor whose history and location it knows; an ad server chooses an advertisement for a query. Only the reward of the chosen action is observed. These are contextual bandit problems, and Chapter 4 of Bubeck and Cesa-Bianchi's monograph Regret Analysis of Stochastic and Nonstochastic Multi-armed Bandit Problems (arXiv:1204.5721v2) surveys several of their formal versions. In a contextual problem the learner is compared with the best policy, a map from contexts to arms, rather than with the best single arm.

This mission covers three of the chapter's models. The first marks each round with a context from a finite set. In the second, NNN experts give advice, as in prediction with expert advice. The third is the bandit multiclass problem: a linear classifier predicts one of KKK labels and then learns only whether its prediction was right. The goal is the mistake bound of the Banditron (Kakade, Shalev-Shwartz and Tewari, ICML 2008). The bound shows that one bit of feedback per round suffices to compete with every linear classifier, at regret O(n2/3)O(n^{2/3})O(n2/3).

Setting

There are K≥2K \ge 2K≥2 arms (or labels) {1,…,K}\{1,\dots,K\}{1,…,K} and rounds t=1,…,nt = 1, \dots, nt=1,…,n.

Adversarial losses. At round ttt an adversary assigns losses ℓi,t∈[0,1]\ell_{i,t} \in [0,1]ℓi,t​∈[0,1] to the arms and may adapt to the forecaster's past plays I1,…,It−1I_1, \dots, I_{t-1}I1​,…,It−1​. The forecaster draws ItI_tIt​ at random from a distribution ptp_tpt​ that depends on what it has observed, and it observes only ℓIt,t\ell_{I_t,t}ℓIt​,t​. Expectations E\mathbb EE are over the forecaster's draws.

Side information. Each round carries a context sts_tst​ from a finite set S\mathcal SS, and the sequence s1,s2,…s_1, s_2, \dotss1​,s2​,… is fixed in advance. The pseudo-regret against context-to-arm maps is

R‾nS=max⁡g:S→{1,…,K}E[∑t=1nℓIt,t−∑t=1nℓg(st),t].\overline R^{\mathcal S}_n = \max_{g:\mathcal S\to\{1,\dots,K\}} \mathbb E\Big[\sum_{t=1}^n \ell_{I_t,t} - \sum_{t=1}^n \ell_{g(s_t),t}\Big].RnS​=g:S→{1,…,K}max​E[t=1∑n​ℓIt​,t​−t=1∑n​ℓg(st​),t​].

The S-Exp3 forecaster runs one instance of Exp3 (Section 3.1 of the book) on each context.

Expert advice. At each round each of NNN experts jjj proposes a distribution ξtj\xi^j_tξtj​ over arms, which may depend on the forecaster's past plays. The contextual pseudo-regret is

R‾nctx=max⁡k=1,…,NE[∑t=1nℓIt,t−∑t=1nEi∼ξtkℓi,t].\overline R^{\mathrm{ctx}}_n = \max_{k=1,\dots,N}\mathbb E\Big[\sum_{t=1}^n \ell_{I_t,t} - \sum_{t=1}^n \mathbb E_{i\sim\xi^k_t}\ell_{i,t}\Big].Rnctx​=k=1,…,Nmax​E[t=1∑n​ℓIt​,t​−t=1∑n​Ei∼ξtk​​ℓi,t​].

Exp4 (Fig. 4.1) runs exponential weights over the experts with importance-weighted loss estimates.

Bandit multiclass. The examples (xt,yt)∈Rd×{1,…,K}(x_t, y_t) \in \mathbb R^d \times \{1,\dots,K\}(xt​,yt​)∈Rd×{1,…,K} are fixed in advance, with ∥xt∥=1\|x_t\| = 1∥xt​∥=1 (Euclidean). A K×dK\times dK×d matrix UUU classifies xxx by arg⁡max⁡i(Ux)i\arg\max_i (Ux)_iargmaxi​(Ux)i​. Its multiclass hinge loss on round ttt is ℓt(U)=[1−(Uxt)yt+max⁡i≠yt(Uxt)i]+\ell_t(U) = [1 - (Ux_t)_{y_t} + \max_{i\neq y_t}(Ux_t)_i]_+ℓt​(U)=[1−(Uxt​)yt​​+maxi=yt​​(Uxt​)i​]+​. Write Ln(U)=∑t≤nℓt(U)L_n(U) = \sum_{t\le n}\ell_t(U)Ln​(U)=∑t≤n​ℓt​(U) for the cumulative hinge loss, Lˉn(U)=Ln(U)/n\bar L_n(U) = L_n(U)/nLˉn​(U)=Ln​(U)/n for its average, and ∥U∥\|U\|∥U∥ for the Frobenius norm. The multiclass Perceptron predicts y^t=arg⁡max⁡i(Wtxt)i\hat y_t = \arg\max_i (W_tx_t)_iy^​t​=argmaxi​(Wt​xt​)i​ and, after seeing yty_tyt​, adds xtx_txt​ to row yty_tyt​ and subtracts it from row y^t\hat y_ty^​t​. The Banditron (p. 58) predicts YtY_tYt​ from pi,t=(1−γ)1y^t=i+γ/Kp_{i,t} = (1-\gamma)\mathbb 1_{\hat y_t = i} + \gamma/Kpi,t​=(1−γ)1y^​t​=i​+γ/K. It observes only 1Yt=yt\mathbb 1_{Y_t = y_t}1Yt​=yt​​ and updates Wt+1=Wt+X~tW_{t+1} = W_t + \widetilde X_tWt+1​=Wt​+Xt​, where (X~t)i,j=xt,j(1Yt=yt1Yt=i/pi,t−1y^t=i)(\widetilde X_t)_{i,j} = x_{t,j}\big(\mathbb 1_{Y_t=y_t}\mathbb 1_{Y_t=i}/p_{i,t} - \mathbb 1_{\hat y_t=i}\big)(Xt​)i,j​=xt,j​(1Yt​=yt​​1Yt​=i​/pi,t​−1y^​t​=i​). Its number of mistakes is Mn=∑t≤n1Yt≠ytM_n = \sum_{t\le n}\mathbb 1_{Y_t\neq y_t}Mn​=∑t≤n​1Yt​=yt​​.

Formalization targets

Goal: Theorem 4.7 (Banditron)

For n≥8Kn \ge 8Kn≥8K, γ=(K/n)1/3\gamma = (K/n)^{1/3}γ=(K/n)1/3, every example sequence as above and every K×dK\times dK×d matrix UUU,

E Mn≤Ln(U)+(1+∥U∥2Lˉn(U))K1/3n2/3+2∥U∥2K2/3n1/3+2 ∥U∥K1/6n1/3.\mathbb E\,M_n \le L_n(U) + \Big(1 + \|U\|\sqrt{2\bar L_n(U)}\Big)K^{1/3}n^{2/3} + 2\|U\|^2K^{2/3}n^{1/3} + \sqrt2\,\|U\|K^{1/6}n^{1/3}.EMn​≤Ln​(U)+(1+∥U∥2Lˉn​(U)​)K1/3n2/3+2∥U∥2K2/3n1/3+2​∥U∥K1/6n1/3.

Milestones

  1. Multiclass Perceptron bound (Section 4.4, p. 57). For every n≥1n \ge 1n≥1 and UUU, ∑t≤n1y^t≠yt≤Ln(U)+2∥U∥2+∥U∥2nLˉn(U)\sum_{t\le n}\mathbb 1_{\hat y_t\ne y_t} \le L_n(U) + 2\|U\|^2 + \|U\|\sqrt{2n\bar L_n(U)}∑t≤n​1y^​t​=yt​​≤Ln​(U)+2∥U∥2+∥U∥2nLˉn​(U)​.
  2. Theorem 4.1 (p. 44). S-Exp3 satisfies R‾nS≤2n∣S∣Kln⁡K\overline R^{\mathcal S}_n \le \sqrt{2n|\mathcal S|K\ln K}RnS​≤2n∣S∣KlnK​.
  3. Theorem 4.2 (p. 46), with corrected constants. Exp4 without mixing satisfies R‾nctx≤2nKln⁡N\overline R^{\mathrm{ctx}}_n \le \sqrt{2nK\ln N}Rnctx​≤2nKlnN​ for ηt=2ln⁡N/(nK)\eta_t = \sqrt{2\ln N/(nK)}ηt​=2lnN/(nK)​, and R‾nctx≤2nKln⁡N\overline R^{\mathrm{ctx}}_n \le 2\sqrt{nK\ln N}Rnctx​≤2nKlnN​ for ηt=ln⁡N/(tK)\eta_t = \sqrt{\ln N/(tK)}ηt​=lnN/(tK)​.
  4. Theorem 4.3 (p. 50), with corrected learning rate. Let the plays be drawn from distributions qtq_tqt​ with qi,t≥ε>0q_{i,t}\ge\varepsilon > 0qi,t​≥ε>0, and let Exp3 run on the estimates ℓi,t1It=i/qi,t\ell_{i,t}\mathbb 1_{I_t=i}/q_{i,t}ℓi,t​1It​=i​/qi,t​ with η=2εln⁡K/n\eta = \sqrt{2\varepsilon\ln K/n}η=2εlnK/n​. Then max⁡kE[∑tEi∼ptℓi,t−∑tℓk,t]≤(2n/ε)ln⁡K\max_k \mathbb E\big[\sum_t \mathbb E_{i\sim p_t}\ell_{i,t} - \sum_t\ell_{k,t}\big] \le \sqrt{(2n/\varepsilon)\ln K}maxk​E[∑t​Ei∼pt​​ℓi,t​−∑t​ℓk,t​]≤(2n/ε)lnK​.

Significance

Theorem 4.7 shows that, on any sequence of examples, the bandit version of online multiclass classification costs at most O(K1/3n2/3)O(K^{1/3}n^{2/3})O(K1/3n2/3) mistakes beyond the hinge loss of the best linear classifier. The full-information Perceptron, by comparison, pays O(n)O(\sqrt n)O(n​). The bound has no stochastic assumption and has explicit constants. Theorems 4.1–4.3 are the basic regret guarantees for side information and expert advice. Theorem 4.3 in particular lets learning algorithms serve as experts inside Exp4, which is the construction behind Theorem 4.5.

The mission produces machine-checked statements, and eventually proofs, of these results with fully explicit constants and an explicit model of adaptive adversaries and adaptive advice. To the curators' knowledge none of the Banditron, the multiclass Perceptron bound, S-Exp3 or Theorem 4.3 is formalized anywhere. The platform's Bandit Algorithms series has a proved Exp4 bound, but only for advice and rewards fixed in advance. The book proves all four milestones and the goal; two printed statements (4.2 and 4.3) contain misprints that this mission corrects.

Difficulty

The Banditron bound concerns a randomized process whose weight matrix depends on all earlier random predictions. The Perceptron argument tracks ⟨U,Wn+1⟩\langle U, W_{n+1}\rangle⟨U,Wn+1​⟩ and ∥Wn+1∥2\|W_{n+1}\|^2∥Wn+1​∥2. It carries over only in conditional expectation, and the second moment of the importance-weighted update is of order K/γK/\gammaK/γ on rounds where y^t≠yt\hat y_t \neq y_ty^​t​=yt​ and of order γ\gammaγ otherwise. Combining these into one inequality for ∑tP(y^t≠yt)\sum_t\mathbb P(\hat y_t\neq y_t)∑t​P(y^​t​=yt​) and then for EMn\mathbb E M_nEMn​ requires solving a quadratic inequality in the presence of expectations, and the constants must come out as printed. For the Exp3/Exp4 results, the obstacle is that losses and advice adapt to past plays. The standard potential argument has to be run conditionally on the history, and a version that fixes the losses in advance proves a weaker theorem.

Formalization scope

  • Rounds and laws. Rounds are numbered from 000 in Lean (Lean round ttt is the book's round t+1t+1t+1). Every forecaster is a sampling rule from past plays to weights on Fin K. The law of the first nnn plays is the product ∏tpt(ωt∣ω<t)\prod_t p_t(\omega_t\mid\omega_{<t})∏t​pt​(ωt​∣ω<t​) over sequences ω:Fin n→Fin K\omega : \mathrm{Fin}\,n\to\mathrm{Fin}\,Kω:Finn→FinK, and expectations are finite sums against it. The adversary and the experts are deterministic functions of past plays; an independent randomized adversary is a mixture of these. The examples of the Banditron are fixed.
  • Argmax. y^t\hat y_ty^​t​ uses any argmax selector; all tie-breaking rules are covered.
  • Norms. ∥xt∥=1\|x_t\| = 1∥xt​∥=1 is the Euclidean condition ∑jxt,j2=1\sum_j x_{t,j}^2 = 1∑j​xt,j2​=1; ∥U∥\|U\|∥U∥ is the Frobenius norm written out explicitly.
  • Infima and maxima. Each "inf⁡U\inf_UinfU​" and "max⁡k\max_kmaxk​" of the book is stated as "for every UUU" or "for every kkk", which is equivalent.
  • Explicit constants. Every bound is the one printed or, for the corrected items, the one the proof yields. No O(⋅)O(\cdot)O(⋅) appears.
  • Corrected misprints. Theorem 4.7 prints the examples in Rd×{−1,+1}\mathbb R^d\times\{-1,+1\}Rd×{−1,+1}; labels are in {1,…,K}\{1,\dots,K\}{1,…,K}. Theorem 4.2 prints 2nNln⁡K\sqrt{2nN\ln K}2nNlnK​ and 2nNln⁡K2\sqrt{nN\ln K}2nNlnK​; the proof gives 2nKln⁡N\sqrt{2nK\ln N}2nKlnN​ and 2nKln⁡N2\sqrt{nK\ln N}2nKlnN​. Theorem 4.3 prints η=2ln⁡K/(nK)\eta = \sqrt{2\ln K/(nK)}η=2lnK/(nK)​; (4.7) follows from the proof with η=2εln⁡K/n\eta = \sqrt{2\varepsilon\ln K/n}η=2εlnK/n​.
  • Parameter range. At n=8Kn = 8Kn=8K the Banditron's γ\gammaγ equals 1/21/21/2, outside the box's open interval (0,1/2)(0,1/2)(0,1/2). The proof uses only γ≤1/2\gamma\le 1/2γ≤1/2, so n=8Kn = 8Kn=8K is included.
  • Ruling out trivial forms. Theorem 4.1 is stated for the explicit S-Exp3 forecaster, not as an existence claim, so no forecaster tuned to the losses can witness it. The losses and the advice are allowed to adapt, so a proof for oblivious sequences does not suffice.
  • Left out. Theorem 4.4 (Exp4 with mixing) is proved in the book only by reference. The argument that reference suggests yields 32γn+Kln⁡N/γ\tfrac32\gamma n + K\ln N/\gamma23​γn+KlnN/γ, not the printed γn/2+Kln⁡N/γ\gamma n/2 + K\ln N/\gammaγn/2+KlnN/γ. Theorem 4.5 is stated with O(⋅)O(\cdot)O(⋅), Theorem 4.6 "for some constant ccc", and Eq. (4.8) is left to the reader.

Useful reusable infrastructure: the path-law expectation for history-dependent sampling, the exponential-weights potential argument under adaptive losses, and Perceptron-type inner-product arguments for matrices. Proofs of any milestone and of the goal are welcome.

Selected references

  • S. Bubeck, N. Cesa-Bianchi, Regret Analysis of Stochastic and Nonstochastic Multi-armed Bandit Problems, Foundations and Trends in Machine Learning 5(1), 2012. arXiv:1204.5721v2. https://arxiv.org/abs/1204.5721 ; https://doi.org/10.1561/2200000024
  • S. M. Kakade, S. Shalev-Shwartz, A. Tewari, Efficient Bandit Algorithms for Online Multiclass Prediction, ICML 2008. https://doi.org/10.1145/1390156.1390212
  • 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
  • O.-A. Maillard, R. Munos, Adaptive Bandits: Towards the Best History-Dependent Strategy, AISTATS 2011. https://proceedings.mlr.press/v15/maillard11a.html
11 thms2 active usersReviewed
Functional AnalysisOperations ResearchOptimization·Captain: mikedeng1

On the Variational Principle II: Under Linearly Independent Active Constraint Gradients, a Bounded-Below Problem Has ε²-Optimal Feasible Points Satisfying the Lagrange Multiplier Rule up to εResearch Paper

Motivation

The classical Lagrange multiplier rule and its inequality-constrained form, the Karush–Kuhn–Tucker (KKT) conditions, are necessary conditions satisfied at a minimizer of a constrained problem. In finite dimensions, a continuous function bounded below on a closed bounded feasible set attains its minimum, so the rule describes an actual point. In an infinite-dimensional Banach space this fails: closed bounded sets are not compact, minimizing sequences need not converge, and a smooth function bounded below on a smooth constraint set can have no minimizer at all. The multiplier rule then has nothing to describe.

Ekeland's 1974 paper On the Variational Principle (DOI) introduced a principle (its Theorem 1.1) that replaces "a minimizer exists" with "near every almost-minimizer there is a point that strictly minimizes a slightly perturbed function". Section 3 of the paper uses this to prove that, for a problem with finitely many smooth equality and inequality constraints satisfying a linear-independence regularity condition, nearly optimal feasible points satisfy the KKT conditions up to a small error, with no compactness and no existence of a minimizer. This is the first appearance of what is now called an approximate or asymptotic KKT condition, a notion central to the convergence theory of nonlinear programming algorithms.

Setting

Let VVV be a real Banach space with dual V∗V^*V∗ (continuous linear functionals) and dual norm ∥x∗∥∗=sup⁡∥h∥≤1⟨x∗,h⟩\|x^*\|_*=\sup_{\|h\|\le1}\langle x^*,h\rangle∥x∗∥∗​=sup∥h∥≤1​⟨x∗,h⟩. Let F:V→RF:V\to\mathbb RF:V→R be Fréchet-differentiable, with derivative F′(v)∈V∗F'(v)\in V^*F′(v)∈V∗, and let G1,…,Gm:V→RG_1,\dots,G_m:V\to\mathbb RG1​,…,Gm​:V→R be C1C^1C1 (continuously Fréchet-differentiable). Fix 0≤p≤m0\le p\le m0≤p≤m and consider

inf⁡F(v)subject toGi(v)=0 (1≤i≤p),Gi(v)≥0 (p+1≤i≤m).(3.1)\inf F(v)\quad\text{subject to}\quad G_i(v)=0\ (1\le i\le p),\qquad G_i(v)\ge0\ (p+1\le i\le m). \tag{3.1}infF(v)subject toGi​(v)=0 (1≤i≤p),Gi​(v)≥0 (p+1≤i≤m).(3.1)

The feasible set is C={v∈V:Gi(v)=0 for i≤p, Gi(v)≥0 for i>p}\mathcal C=\{v\in V : G_i(v)=0 \text{ for } i\le p,\ G_i(v)\ge0 \text{ for } i>p\}C={v∈V:Gi​(v)=0 for i≤p, Gi​(v)≥0 for i>p} (3.2). At v∈Cv\in\mathcal Cv∈C the saturated constraints are I(v)={i:Gi(v)=0}I(v)=\{i : G_i(v)=0\}I(v)={i:Gi​(v)=0} (3.3). The regularity assumption (3.4) is: for every v∈Cv\in\mathcal Cv∈C, the derivatives Gi′(v)G_i'(v)Gi′​(v), i∈I(v)i\in I(v)i∈I(v), are linearly independent in V∗V^*V∗.

In the Lean development these are EkelandVP.Constraints.feasibleSet p G and EkelandVP.Constraints.IsRegular p G, with constraints indexed by Fin m.

Formalization targets

Goal: Theorem 3.1 (p. 330)

Assume (3.4), C≠∅\mathcal C\ne\emptysetC=∅, and that FFF is bounded below on C\mathcal CC (3.5). Then for every ε>0\varepsilon>0ε>0 there are vε∈Cv_\varepsilon\in\mathcal Cvε​∈C and λ1,…,λm∈R\lambda_1,\dots,\lambda_m\in\mathbb Rλ1​,…,λm​∈R with

F(vε)≤inf⁡CF+ε2,λi≥0 (i>p),λi=0 if Gi(vε)≠0,F(v_\varepsilon)\le\inf_{\mathcal C}F+\varepsilon^2,\qquad \lambda_i\ge0\ (i>p),\qquad \lambda_i=0 \text{ if } G_i(v_\varepsilon)\ne0,F(vε​)≤Cinf​F+ε2,λi​≥0 (i>p),λi​=0 if Gi​(vε​)=0, ∥F′(vε)−∑i=1mλiGi′(vε)∥∗≤ε.\Big\|F'(v_\varepsilon)-\sum_{i=1}^m\lambda_iG_i'(v_\varepsilon)\Big\|_*\le\varepsilon.​F′(vε​)−i=1∑m​λi​Gi′​(vε​)​∗​≤ε.

Milestones, in the order of the paper's proof

  1. (3.8)–(3.12): a feasible vvv with F(v)≤inf⁡CF+ε2F(v)\le\inf_{\mathcal C}F+\varepsilon^2F(v)≤infC​F+ε2 and F(w)≥F(v)−ε∥w−v∥F(w)\ge F(v)-\varepsilon\|w-v\|F(w)≥F(v)−ε∥w−v∥ for all w∈Cw\in\mathcal Cw∈C (the variational principle applied to FFF restricted to C\mathcal CC; no regularity needed).
  2. (3.16): at a regular feasible point vvv, every hhh with ⟨Gi′(v),h⟩=0\langle G_i'(v),h\rangle=0⟨Gi′​(v),h⟩=0 (i≤pi\le pi≤p) and ⟨Gi′(v),h⟩≥0\langle G_i'(v),h\rangle\ge0⟨Gi′​(v),h⟩≥0 (i>pi>pi>p, i∈I(v)i\in I(v)i∈I(v)) is the initial velocity of a C1C^1C1 curve u:[0,τ]→Cu:[0,\tau]\to\mathcal Cu:[0,τ]→C with u(0)=vu(0)=vu(0)=v.
  3. Lemma 3.2: at a point with the property of milestone 1, ⟨F′(v),h⟩≥−ε∥h∥\langle F'(v),h\rangle\ge-\varepsilon\|h\|⟨F′(v),h⟩≥−ε∥h∥ for every such hhh.
  4. Lemma 3.3: an ε\varepsilonε-Farkas–Minkowski lemma in V∗V^*V∗: if ⟨w∗,h⟩≥−ε∥h∥\langle w^*,h\rangle\ge-\varepsilon\|h\|⟨w∗,h⟩≥−ε∥h∥ whenever ⟨ui∗,h⟩=0\langle u_i^*,h\rangle=0⟨ui∗​,h⟩=0 and ⟨vj∗,h⟩≥0\langle v_j^*,h\rangle\ge0⟨vj∗​,h⟩≥0, then ∥w∗−∑λiui∗−∑μjvj∗∥∗≤ε\|w^*-\sum\lambda_iu_i^*-\sum\mu_jv_j^*\|_*\le\varepsilon∥w∗−∑λi​ui∗​−∑μj​vj∗​∥∗​≤ε for some λi∈R\lambda_i\in\mathbb Rλi​∈R and μj≥0\mu_j\ge0μj​≥0.

An additional item states Corollary 3.4 (p. 333), the one-constraint case: if G(v)=0⇒G′(v)≠0G(v)=0\Rightarrow G'(v)\ne0G(v)=0⇒G′(v)=0, {G=0}≠∅\{G=0\}\neq\emptyset{G=0}=∅ and FFF is bounded below on {G=0}\{G=0\}{G=0}, then for every ε>0\varepsilon>0ε>0 there are vεv_\varepsilonvε​ with G(vε)=0G(v_\varepsilon)=0G(vε​)=0 and λε∈R\lambda_\varepsilon\in\mathbb Rλε​∈R with ∥F′(vε)−λεG′(vε)∥∗≤ε\|F'(v_\varepsilon)-\lambda_\varepsilon G'(v_\varepsilon)\|_*\le\varepsilon∥F′(vε​)−λε​G′(vε​)∥∗​≤ε.

Significance

The result. Theorem 3.1 is an existence theorem for approximate KKT points that needs neither compactness nor attainment of the infimum. It shows that every bounded-below problem with regular constraints has a sequence of feasible points whose objective values converge to the infimum and along which the KKT residual tends to zero. This is the property that later work calls approximate KKT or asymptotic KKT (AKKT) and uses as a stopping criterion and as a sequential optimality condition for nonlinear programming. Corollary 3.4 is the corresponding nonlinear eigenvalue statement: on a regular level set, F′F'F′ is approximately proportional to G′G'G′ at almost-minimizing points.

Formalizing it. The result is classical and proved in the paper; to our knowledge no machine-checked version exists. Mathlib has the Fréchet derivative, the implicit function theorem for strictly differentiable maps, Banach–Alaoglu and the Hahn–Banach separation theorems, but no Ekeland principle in this form, no Lyusternik-type tangent-curve theorem for mixed equality–inequality constraints, and no Farkas lemma in a dual Banach space. Each milestone is a reusable piece of nonlinear optimization theory in Banach spaces.

Difficulty

The obvious argument, "take a minimizer and apply the Lagrange multiplier rule", fails at the first step because no minimizer need exist. The variational principle supplies a point vεv_\varepsilonvε​ that minimizes F+ε∥⋅−vε∥F+\varepsilon\|\cdot-v_\varepsilon\|F+ε∥⋅−vε​∥ on C\mathcal CC, but that function is not differentiable at vεv_\varepsilonvε​, so the multiplier rule cannot be applied to it directly either. Two further gaps remain. Linearized feasible directions (those satisfying the derivative conditions on the active constraints) need not be directions along which one can actually move inside C\mathcal CC; closing this gap requires the regularity assumption and completeness of VVV, and must keep the active inequality constraints nonnegative, not just the equalities. And the resulting first-order inequality, which holds only up to ε∥h∥\varepsilon\|h\|ε∥h∥, must be turned into an approximate multiplier representation in V∗V^*V∗, an infinite-dimensional dual space in which the usual finite-dimensional Farkas lemma does not apply as stated.

Formalization scope

  • VVV is a real Banach space: [NormedAddCommGroup V] [NormedSpace ℝ V] [CompleteSpace V]. V∗V^*V∗ is V →L[ℝ] ℝ with the operator norm; F′(v)F'(v)F′(v) is fderiv ℝ F v.
  • Constraints are one family G : Fin m → V → ℝ, 0-based: the paper's constraint iii is Lean index i−1i-1i−1, an equality constraint iff its index is <p<p<p. p ≤ m is assumed in the goal.
  • C1C^1C1 is ContDiff ℝ 1; FFF is Differentiable ℝ F (Fréchet-differentiable everywhere).
  • No infimum over C\mathcal CC is formed: "bounded below" is BddBelow (F '' 𝒞) and "F(v)≤inf⁡CF+ε2F(v)\le\inf_{\mathcal C}F+\varepsilon^2F(v)≤infC​F+ε2" is "F(v)≤F(w)+ε2F(v)\le F(w)+\varepsilon^2F(v)≤F(w)+ε2 for all w∈Cw\in\mathcal Cw∈C". A real infimum over an empty or unbounded set would be a junk value.
  • Added hypotheses, disclosed in each item: C≠∅\mathcal C\ne\emptysetC=∅ (goal and milestone 1) and {G=0}≠∅\{G=0\}\ne\emptyset{G=0}=∅ (Corollary 3.4). Without them the paper's hypotheses hold vacuously (inf⁡∅=+∞\inf\emptyset=+\inftyinf∅=+∞) while the conclusion asks for a feasible point.
  • Lemma 3.3 is stated with ≤ε\le\varepsilon≤ε. The paper prints <ε<\varepsilon<ε in (3.20), which fails for V=RV=\mathbb RV=R, no constraints and w∗=ε idw^*=\varepsilon\,\mathrm{id}w∗=εid; its proof gives ≤\le≤, and Theorem 3.1 uses ≤\le≤.
  • Lemma 3.2 and milestone 2 assume the linear independence (3.4) only at the point vvv under consideration, and Lemma 3.2 is stated for any feasible vvv with property (3.12); this is exactly what the paper's proof uses.
  • Regularity is a linear independence of the indexed family (Gi′(v))i∈I(v)(G_i'(v))_{i\in I(v)}(Gi′​(v))i∈I(v)​, so a repeated saturated constraint violates it; the constraint qualification cannot be trivialized by collapsing duplicates.

Contributions welcome: a general Ekeland principle with the strict-minimizer conclusion, a Lyusternik–Graves tangent-curve theorem for C1C^1C1 maps with surjective derivative onto Rk\mathbb R^kRk, and a closedness result for finitely generated cones in V∗V^*V∗.

Selected references

  • I. Ekeland, On the Variational Principle, J. Math. Anal. Appl. 47 (1974) 324–353. https://doi.org/10.1016/0022-247X(74)90025-0
  • I. Ekeland, Nonconvex minimization problems, Bull. Amer. Math. Soc. (N.S.) 1 (1979) 443–474. https://doi.org/10.1090/S0273-0979-1979-14595-6
  • R. Andreani, G. Haeser, J. M. Martínez, On sequential optimality conditions for smooth constrained optimization, Optimization 60 (2011) 627–641. https://doi.org/10.1080/02331930903578700
7 thms1 active userReviewed
Convex OptimizationLinear algebraOperations Research+1·Captain: mikedeng1

A Nonlinear Programming Algorithm for Solving Semidefinite Programs via Low-rank Factorization: A Regular Local Minimum That Stays Locally Minimal After Adding a Zero Column Solves the SDPResearch Paper

Motivation

Semidefinite programs (SDPs) arise as convex relaxations of combinatorial problems such as maximum cut and the Lovász theta function, and in control and eigenvalue optimization. Interior-point methods solve them reliably but manipulate dense n×nn\times nn×n matrices, which limits the size of the instances they can handle. Burer and Monteiro (Math. Program. 95 (2003)) proposed replacing the matrix variable X⪰0X\succeq 0X⪰0 by a factorization X=RRTX=RR^{T}X=RRT with RRR having only rrr columns, and solving the resulting nonconvex program by a first-order augmented Lagrangian method. The approach rests on a theorem of Barvinok (1995) and Pataki (1998): an SDP with mmm linear constraints has an optimal solution of rank rrr with r(r+1)/2≤mr(r+1)/2\le mr(r+1)/2≤m, so a small number of columns suffices.

Because the factorized problem is nonconvex, a local minimum it returns is not automatically a solution of the SDP. Section 2 of the paper gives conditions under which it is. This mission formalizes those conditions, culminating in Proposition 2.5, which justifies the paper's strategy of increasing the rank one column at a time.

Setting

For real p×qp\times qp×q matrices, the trace inner product is A∙B=trace⁡(ATB)A\bullet B=\operatorname{trace}(A^{T}B)A∙B=trace(ATB). The data are symmetric matrices C,A1,…,Am∈SnC, A_1,\dots,A_m\in\mathcal S^nC,A1​,…,Am​∈Sn and a vector b∈Rmb\in\mathbb R^mb∈Rm. The primal SDP and dual SDP are

(1)min⁡{C∙X:Ai∙X=bi, i=1,…,m, X⪰0},(3)max⁡{bTy:S=C−∑i=1myiAi, S⪰0}.\text{(1)}\quad \min\{C\bullet X : A_i\bullet X=b_i,\ i=1,\dots,m,\ X\succeq0\},\qquad \text{(3)}\quad \max\Big\{b^{T}y : S=C-\sum_{i=1}^m y_iA_i,\ S\succeq0\Big\}.(1)min{C∙X:Ai​∙X=bi​, i=1,…,m, X⪰0},(3)max{bTy:S=C−i=1∑m​yi​Ai​, S⪰0}.

The standing assumptions of the paper are that A1,…,AmA_1,\dots,A_mA1​,…,Am​ are linearly independent and that there are feasible X∗X^*X∗ and (S∗,y∗)(S^*,y^*)(S∗,y∗) with C∙X∗=bTy∗C\bullet X^*=b^{T}y^*C∙X∗=bTy∗.

For a positive integer r≤nr\le nr≤n, the low-rank program is

(Nr)min⁡{C∙(RRT):Ai∙(RRT)=bi, i=1,…,m, R∈Rn×r}.(N_r)\qquad \min\{C\bullet(RR^{T}) : A_i\bullet(RR^{T})=b_i,\ i=1,\dots,m,\ R\in\mathbb R^{n\times r}\}.(Nr​)min{C∙(RRT):Ai​∙(RRT)=bi​, i=1,…,m, R∈Rn×r}.

Its Lagrangian is L(R,y)=C∙(RRT)−∑iyi(Ai∙(RRT)−bi)L(R,y)=C\bullet(RR^{T})-\sum_i y_i(A_i\bullet(RR^{T})-b_i)L(R,y)=C∙(RRT)−∑i​yi​(Ai​∙(RRT)−bi​), and S(y)=C−∑iyiAiS(y)=C-\sum_i y_iA_iS(y)=C−∑i​yi​Ai​. A feasible RRR is a local minimum if it minimizes the objective among nearby feasible points; it is a regular point if A1R,…,AmRA_1R,\dots,A_mRA1​R,…,Am​R are linearly independent; it is a stationary point with multiplier yyy if ∇RL(R,y)=0\nabla_RL(R,y)=0∇R​L(R,y)=0. The injection of R∈Rn×rR\in\mathbb R^{n\times r}R∈Rn×r is R^=[ R  0 ]∈Rn×(r+1)\hat R=[\,R\ \ 0\,]\in\mathbb R^{n\times(r+1)}R^=[R  0]∈Rn×(r+1), obtained by appending a zero column.

Formalization targets

Goal: Proposition 2.5

Let r<nr<nr<n and let R∗R^*R∗ be a regular local minimum of (Nr)(N_r)(Nr​) with multiplier y∗y^*y∗, S∗=S(y∗)S^*=S(y^*)S∗=S(y∗), S∗R∗=0S^*R^*=0S∗R∗=0. If R^\hat RR^ is a local minimum of (Nr+1)(N_{r+1})(Nr+1​), then

X∗=R∗(R∗)T solves (1)and(S∗,y∗) solves (3).X^*=R^*(R^*)^{T}\ \text{solves (1)}\quad\text{and}\quad (S^*,y^*)\ \text{solves (3)}.X∗=R∗(R∗)T solves (1)and(S∗,y∗) solves (3).

Milestones

  1. The derivative formulas (9): ∇R(Ai∙(RRT)−bi)=2AiR\nabla_R(A_i\bullet(RR^T)-b_i)=2A_iR∇R​(Ai​∙(RRT)−bi​)=2Ai​R, ∇RL(R,y)=2SR\nabla_RL(R,y)=2SR∇R​L(R,y)=2SR, and LRR′′(R,y)[D,D]=2S∙(DDT)L''_{RR}(R,y)[D,D]=2S\bullet(DD^T)LRR′′​(R,y)[D,D]=2S∙(DDT).
  2. Proposition 2.3: at a regular local minimum of (Nr)(N_r)(Nr​) there is a unique y∗y^*y∗ with S∗R∗=0S^*R^*=0S∗R∗=0, and S∗∙(DDT)≥0S^*\bullet(DD^T)\ge0S∗∙(DDT)≥0 for every DDD with AiR∗∙D=0A_iR^*\bullet D=0Ai​R∗∙D=0 for all iii.
  3. Proposition 2.1: feasible XXX and (S,y)(S,y)(S,y) are simultaneously optimal if and only if X∙S=0X\bullet S=0X∙S=0.
  4. Proposition 2.4: a stationary point of (Nr)(N_r)(Nr​) whose S∗S^*S∗ is positive semidefinite gives optimal X∗=R∗R∗TX^*=R^*R^{*T}X∗=R∗R∗T and (S∗,y∗)(S^*,y^*)(S∗,y∗).

Significance

Proposition 2.5 is a certificate of global optimality for a nonconvex problem obtained from local information alone. It is the basis of the rank-increase scheme described on p. 8 of the paper: compute a local minimum of (Nr)(N_r)(Nr​) for a small rrr; if the zero-column extension is still a local minimum of (Nr+1)(N_{r+1})(Nr+1​), the current point solves the SDP; otherwise a better point of (Nr+1)(N_{r+1})(Nr+1​) exists and rrr is increased. Proposition 2.4 gives the companion test, valid for every rrr: positive semidefiniteness of the multiplier matrix at a stationary point. These statements underlie the later convergence analysis of the method (Burer & Monteiro 2005) and the literature on benign landscapes of low-rank SDP formulations (Boumal, Voroninski & Bandeira 2016).

The results are proved in the paper. What this mission adds is a machine-checked version of the full chain from the standard-form SDP to the rank-increase certificate, including the matrix calculus (9), the first- and second-order necessary conditions for an equality-constrained program over rectangular matrices, and SDP complementary slackness in standard form. No machine-checked proof of these results is recorded in Mathlib or on the platform.

Difficulty

The SDP side (Propositions 2.1 and 2.4) is linear algebra: weak duality and the fact that the trace inner product of two positive semidefinite matrices is nonnegative. The substance lies in Proposition 2.3. The feasible set of (Nr)(N_r)(Nr​) is a variety cut out by mmm quadratic equations, and the multiplier rule and, especially, the second-order necessary condition require a constraint qualification and a curve in the feasible set realizing every tangent direction. Mathlib provides a first-order Lagrange multiplier rule, but not the second-order condition on the tangent space. A naive attempt to read Proposition 2.5 off Proposition 2.4 fails: local minimality of R∗R^*R∗ alone does not make S∗S^*S∗ positive semidefinite (when rrr is below the minimal optimal rank, it is not); the hypothesis on (Nr+1)(N_{r+1})(Nr+1​) is indispensable.

Formalization scope

Matrices are Matrix (Fin n) (Fin r) ℝ with 0-based indices. The trace inner product is frob A B = trace(Aᵀ * B), defined for rectangular matrices. The data carry explicit symmetry hypotheses C.IsSymm and (A i).IsSymm; without them the formulas (9) are false. Primal feasibility uses Mathlib's PosSemidef, which over R\mathbb RR includes symmetry. Optimality for (1) and (3) is defined relative to their entire feasible sets. The standing assumptions are a separate predicate carried as a hypothesis by Propositions 2.1, 2.3, 2.4 and 2.5, and every statement about (Nr)(N_r)(Nr​) carries 0<r0<r0<r and r≤nr\le nr≤n (or r<nr<nr<n). Gradients are Fréchet derivatives under the Frobenius norm, identified with matrices through the trace inner product; local minima use IsLocalMinOn on the feasible set of (Nr)(N_r)(Nr​) together with feasibility. The injection appends the zero column as the last column.

The statement admits several trivializing encodings, all excluded here: optimality defined relative to the factorized feasible set instead of the whole SDP, an empty or unconstrained (Nr)(N_r)(Nr​) (an unconstrained local minimum or a local minimum without feasibility), a stationarity notion that already includes S⪰0S\succeq0S⪰0, and an injection other than the zero-column extension.

A complete development needs the matrix calculus of R↦RRTR\mapsto RR^{T}R↦RRT, a second-order necessary optimality condition under linear independence of the constraint gradients, and standard-form SDP weak duality and complementary slackness; all of these are reusable well beyond this mission. Proofs of individual milestones, in particular the derivative formulas and Proposition 2.4, are welcome independently of the goal.

Selected references

  • S. Burer and R. D. C. Monteiro, A nonlinear programming algorithm for solving semidefinite programs via low-rank factorization, Mathematical Programming 95 (2003), 329–357. https://doi.org/10.1007/s10107-002-0352-8 (statements cited from the authors' manuscript of March 9, 2001)
  • A. Barvinok, Problems of distance geometry and convex properties of quadratic maps, Discrete & Computational Geometry 13 (1995), 189–202. https://doi.org/10.1007/BF02574037
  • G. Pataki, On the rank of extreme matrices in semidefinite programs and the multiplicity of optimal eigenvalues, Mathematics of Operations Research 23 (1998), 339–358. https://doi.org/10.1287/moor.23.2.339
  • R. D. C. Monteiro and M. Todd, Path-following methods for semidefinite programming, in Handbook of Semidefinite Programming, Kluwer, 2000 (source of Proposition 2.1).
  • S. Burer and R. D. C. Monteiro, Local minima and convergence in low-rank semidefinite programming, Mathematical Programming 103 (2005), 427–444. https://doi.org/10.1007/s10107-004-0564-1
  • N. Boumal, V. Voroninski and A. S. Bandeira, The non-convex Burer–Monteiro approach works on smooth semidefinite programs, NeurIPS 2016. https://arxiv.org/abs/1606.04970
10 thms2 active usersReviewed
Graph TheoryLinear OptimizationOperations Research+1·Captain: mikedeng1

Project Scheduling with Time Windows and Scarce Resources VIII: A Vertex Schedule Maximizes the Net Present Value iff Its Spanning-Tree Subprojects Have the Right SignsTextbook

Motivation

Long-running projects such as construction, plant engineering or software development involve payments to and from the contractor at many points in time: disbursements when activities are carried out, progress payments when milestones are reached. When the planning horizon is long, money received later is worth less, and the natural financial objective is the net present value of all cash flows. Scheduling a project to maximize its net present value subject to minimum and maximum time lags was studied by Russell (1970) and Grinold (1972), and the problem is the prototype of a nonregular objective: delaying an activity can be profitable, because disbursements lose value when they are postponed.

This mission follows Chapter 3 of Neumann, Schwindt and Zimmermann, Project Scheduling with Time Windows and Scarce Resources (2nd ed., Springer 2003). The book shows that the net present value objective belongs to the class of binary-monotone objective functions (§3.3.5), and it uses this in §3.9.1 to give a combinatorial optimality criterion for the resource-free problem: a vertex schedule is optimal exactly when the subprojects cut off by the arcs of a spanning tree have net present values of the right sign (Proposition 3.9.2). That criterion drives the book's parametric analysis of the net present value as a function of the discount rate and the deadline.

Setting

A project consists of activities V={0,1,…,n+1}V=\{0,1,\dots,n+1\}V={0,1,…,n+1}, n≥1n\ge1n≥1, where 000 is the project beginning and n+1n+1n+1 the project completion. Activity iii has an integer duration pip_ipi​, with p0=pn+1=0p_0=p_{n+1}=0p0​=pn+1​=0 and pi>0p_i>0pi​>0 otherwise. Temporal constraints are the arcs of a project network N=⟨V,E;δ⟩N=\langle V,E;\delta\rangleN=⟨V,E;δ⟩: an arc ⟨i,j⟩\langle i,j\rangle⟨i,j⟩ with integer weight δij\delta_{ij}δij​ requires Sj−Si≥δijS_j-S_i\ge\delta_{ij}Sj​−Si​≥δij​ for the start times SiS_iSi​. A maximum project duration dˉ\bar ddˉ is the arc ⟨n+1,0⟩\langle n+1,0\rangle⟨n+1,0⟩ with weight −dˉ-\bar d−dˉ. The time-feasible region is

ST={S∈R≥0n+2∣S0=0, Sj−Si≥δij (⟨i,j⟩∈E)}.\mathcal S_T=\{S\in\mathbb R^{n+2}_{\ge0}\mid S_0=0,\ S_j-S_i\ge\delta_{ij}\ (\langle i,j\rangle\in E)\}.ST​={S∈R≥0n+2​∣S0​=0, Sj​−Si​≥δij​ (⟨i,j⟩∈E)}.

Let 0<β≤10<\beta\le10<β≤1 be the discount rate (β=1/(1+I)\beta=1/(1+I)β=1/(1+I) for an interest rate III) and ciF∈Rc_i^F\in\mathbb RciF​∈R the cash flow of activity iii, paid at its completion time Ci=Si+piC_i=S_i+p_iCi​=Si​+pi​. The problem (3.9.1) is

minimize f(S)=−∑i∈VciFβSi+pisubject to S∈ST,\text{minimize } f(S)=-\sum_{i\in V}c_i^F\beta^{S_i+p_i}\quad\text{subject to } S\in\mathcal S_T,minimize f(S)=−i∈V∑​ciF​βSi​+pi​subject to S∈ST​,

and a minimizer is a time-optimal schedule. A vertex of ST\mathcal S_TST​ is an extreme point. A spanning tree G=⟨V,EG⟩G=\langle V,E^G\rangleG=⟨V,EG⟩ is associated with SSS if EG⊆EE^G\subseteq EEG⊆E, EGE^GEG has n+1n+1n+1 arcs and a connected underlying undirected graph, and SSS is the unique solution of S0=0S_0=0S0​=0, Sj−Si=δijS_j-S_i=\delta_{ij}Sj​−Si​=δij​ for ⟨i,j⟩∈EG\langle i,j\rangle\in E^G⟨i,j⟩∈EG. Deleting a tree arc ⟨i,j⟩\langle i,j\rangle⟨i,j⟩ splits GGG into two subtrees; VijV_{ij}Vij​ is the node set of the one not containing 000. The arc is forward if the tree path from 000 passes it from iii to jjj and backward otherwise, and

npvij(S)=∑h∈VijchFβSh+phnpv^{ij}(S)=\sum_{h\in V_{ij}}c_h^F\beta^{S_h+p_h}npvij(S)=h∈Vij​∑​chF​βSh​+ph​

is the net present value of the subproject VijV_{ij}Vij​. Finally, fff is binary-monotone if it is monotone on every line {S+λz≥0∣λ∈R}\{S+\lambda z\ge0\mid\lambda\in\mathbb R\}{S+λz≥0∣λ∈R} with direction z∈{0,1}n+2z\in\{0,1\}^{n+2}z∈{0,1}n+2 (Definition 3.3.2).

Formalization targets

Goal: Proposition 3.9.2, pinned reading

Assume every node is reached from 000 by a path of nonnegative length (the standing convention of §1.2) and let SSS be a vertex of ST\mathcal S_TST​.

(sufficiency)G associated with S,  npvij(S)≥0 on forward arcs, npvij(S)≤0 on backward arcs ⟹ S time-optimal;\text{(sufficiency)}\quad G \text{ associated with } S,\ \ npv^{ij}(S)\ge0 \text{ on forward arcs},\ npv^{ij}(S)\le0 \text{ on backward arcs}\ \Longrightarrow\ S \text{ time-optimal};(sufficiency)G associated with S,  npvij(S)≥0 on forward arcs, npvij(S)≤0 on backward arcs ⟹ S time-optimal; (necessity, β<1)S time-optimal ⟹ ∃ G associated with S satisfying the sign conditions.\text{(necessity, } \beta<1)\quad S \text{ time-optimal}\ \Longrightarrow\ \exists\, G \text{ associated with } S \text{ satisfying the sign conditions}.(necessity, β<1)S time-optimal ⟹ ∃G associated with S satisfying the sign conditions.

The book states "if and only if … for each arc of the corresponding spanning tree", where the corresponding tree is chosen using optimality. The two directions above are the reading that makes the statement well defined: sufficiency for every associated tree, necessity for some associated tree.

Milestones

  1. §3.3.5: the net present value objective is binary-monotone and sum-separable.
  2. §3.9.1: if ST\mathcal S_TST​ is nonempty and bounded, some vertex of ST\mathcal S_TST​ is time-optimal.
  3. Proposition 3.2.16: every vertex of ST\mathcal S_TST​ has an associated spanning tree, an outtree rooted at 000 if the vertex is a minimal point.
  4. Proposition 3.5.4: a directed forest with at least one node has a source with at most one successor or a sink with exactly one predecessor.

Significance

Proposition 3.9.2 turns a nonconvex continuous optimization problem into a finite check on a spanning tree. Read as an economic statement, it says that at an optimal schedule no subproject with positive net present value can be started earlier and no subproject with negative net present value can be postponed. The book builds on it the parametric procedure of §3.9.1, which tracks the optimal tree as the discount rate or the deadline varies (Propositions 3.9.3 and 3.9.4), and the steepest descent method of §3.5.2 terminates exactly when the criterion holds.

The results are proved in the book, partly by reference to network optimization (Ahuja et al., 1993) and to Schwindt and Zimmermann (2001, 2002). None of them is formalized on the platform or, as far as is known, anywhere else. A formal proof would give the first machine-checked optimality certificate for a nonregular project scheduling objective, and the spanning-tree description of vertices (Proposition 3.2.16) is shared with Mission VI of this series.

Difficulty

The objective fff is neither convex nor concave when cash flows of both signs occur, so local optimality at a vertex does not imply global optimality by a convexity argument, and a first-order check along the edges of ST\mathcal S_TST​ is not obviously enough. The criterion is also not a statement about one tree: a degenerate vertex, where more than n+1n+1n+1 temporal constraints are binding, has several associated trees, and the sign conditions may hold on some and fail on others. Necessity therefore requires producing a suitable tree, not checking a given one. Finally, the combinatorial objects (the subtree VijV_{ij}Vij​, forward and backward orientation relative to the root) have to be connected to the geometry of ST\mathcal S_TST​ through Proposition 3.2.16, whose proof in the book is a citation.

Formalization scope

Activities are Fin (n + 2) with 0 the project beginning and Fin.last (n+1) the project completion; start times are real; durations are natural numbers and arc weights integers. The deadline is a structure field together with the backward arc ⟨n+1,0⟩\langle n+1,0\rangle⟨n+1,0⟩ of weight −dˉ-\bar d−dˉ. βx\beta^xβx is Real.rpow, and every statement assumes 0<β≤10<\beta\le10<β≤1 as the book does (p. 203). Vertices are Set.extremePoints ℝ. A spanning tree is a Finset of n+1n+1n+1 arcs whose SimpleGraph.fromRel is connected; VijV_{ij}Vij​ is the set of nodes not reachable from 000 once the arc is deleted.

Three readings are committed and disclosed in the item statements. Necessity is stated only for β<1\beta<1β<1: at β=1\beta=1β=1 the objective is constant, every schedule is optimal, and the sign conditions can fail on every tree. The standing convention of §1.2 (a path of nonnegative length from 000 to every node) is a hypothesis of Proposition 3.2.16 and of the goal; without it a vertex can be fixed by Si≥0S_i\ge0Si​≥0 rather than by arcs of NNN, and necessity fails. The existence of an optimal vertex assumes ST\mathcal S_TST​ nonempty and bounded, which the book asserts in §3.1. Chapter 3's resource constraints do not occur in this mission, which concerns PS∞∣temp,dˉ∣fPS\infty|temp,\bar d|fPS∞∣temp,dˉ∣f only.

The goal cannot be discharged by choosing the tree freely: associated trees must consist of arcs of NNN that are binding at SSS and determine SSS uniquely, and sufficiency must hold for every such tree. Contributions welcome beyond the milestones: a proof of Proposition 3.2.16 reusable by Mission VI, and a general lemma relating binding spanning trees of difference constraints to extreme points.

Selected references

  • K. Neumann, C. Schwindt, J. Zimmermann, Project Scheduling with Time Windows and Scarce Resources, 2nd ed., Springer, 2003, §3.1 (p. 203), §3.3.5 (pp. 224–225), §3.5.2 (p. 252), §3.9.1 (pp. 333–334). https://doi.org/10.1007/978-3-540-24800-2
  • A. H. Russell, "Cash flows in networks", Management Science 16 (1970), 357–373. https://doi.org/10.1287/mnsc.16.5.357
  • R. C. Grinold, "The payment scheduling problem", Naval Research Logistics Quarterly 19 (1972), 123–136.
  • C. Schwindt, J. Zimmermann, "A steepest ascent approach to maximizing the net present value of projects", Mathematical Methods of Operations Research 53 (2001), 435–450.
  • C. Schwindt, J. Zimmermann, "Parametrische Optimierung als Instrument zur Bewertung von Investitionsprojekten", Zeitschrift für Betriebswirtschaft 72 (2002), 593–617.
  • R. K. Ahuja, T. L. Magnanti, J. B. Orlin, Network Flows, Prentice Hall, 1993.
  • C. Berge, Graphs and Hypergraphs, North-Holland, Amsterdam, 1976.
9 thms2 active usersReviewed
Algorithmic Game TheoryMechanism DesignOperations Research+1·Captain: mikedeng1

Multi-parameter Mechanism Design and Sequential Posted Pricing 4: A 6.75-Approximate Truthful Posted-Price Menu for Unit-Demand Buyers of Multiple ItemsResearch Paper

Motivation

A hotel sells rooms of several types, in limited numbers, to guests who each want one room. The revenue-optimal way to sell is known only in special cases: for buyers with several private values, optimal mechanisms can be randomized, involve lotteries, and lack a closed form (Manelli–Vincent 2007; Chawla, Hartline, Kleinberg 2007). In practice sellers post prices. The question is how much revenue posting prices gives up.

Chawla, Hartline, Malec and Sivan (arXiv:0907.2435v2, STOC 2010) answer it for a broad class of single- and multi-parameter problems. For unit-demand buyers of multiple copies of multiple items they show that a menu of posted prices, offered to the buyers in whatever order they arrive, earns at least 1/6.751/6.751/6.75 of the revenue of any deterministic truthful mechanism (Theorem 14). This mission formalizes that result together with the two steps it is built from: a reduction from the multi-parameter problem to a single-parameter one with "copies" of each buyer (Lemma 3, Theorem 4), and an order-oblivious pricing for the intersection of two partition matroids (Theorem 13).

Setting

Single-parameter problem (BSMD). Finitely many agents iii have independent private values vi∼Fiv_i \sim F_ivi​∼Fi​, each with a density on a bounded interval. A seller may serve any set in a downward-closed set system J\mathcal JJ. A deterministic mechanism MMM maps reported values vvv to a served set M(v)∈JM(v) \in \mathcal JM(v)∈J and payments πi(v)\pi_i(v)πi​(v); it is truthful if reporting the true value is a dominant strategy and no agent ends with negative utility. Its expected revenue is RM=Ev[∑iπi(v)]\mathcal R^M = \mathbb E_v[\sum_i \pi_i(v)]RM=Ev​[∑i​πi​(v)]. For prices ppp, agent iii desires service if pi≤vip_i \le v_ipi​≤vi​, and Sv\mathcal S_vSv​ is the class of maximal feasible sets of desiring agents. The order-oblivious revenue is

Rpobl=Ev[min⁡S∈Sv∑i∈Spi],\mathcal R^{\mathrm{obl}}_{\mathbf p} = \mathbb E_{v}\Big[\min_{S \in \mathcal S_v} \sum_{i \in S} p_i\Big],Rpobl​=Ev​[S∈Sv​min​i∈S∑​pi​],

a lower bound on the revenue of posting the prices ppp to the agents in an adversarial order.

Multi-parameter unit-demand problem (BMUMD). There are mmm buyers and a finite set JJJ of services, partitioned into the groups JiJ_iJi​ of services targeted at buyer iii. Buyer iii has value vjv_jvj​ for each j∈Jij \in J_ij∈Ji​, all values independent with vj∼Fjv_j \sim F_jvj​∼Fj​, and the set system J⊆2J\mathcal J \subseteq 2^JJ⊆2J is unit-demand: ∣S∩Ji∣≤1|S \cap J_i| \le 1∣S∩Ji​∣≤1 for feasible SSS. A mechanism A\mathcal AA is truthful if no buyer gains by misreporting its whole vector (vj)j∈Ji(v_j)_{j \in J_i}(vj​)j∈Ji​​, and individually rational if a buyer receiving jjj pays at most vjv_jvj​ and a buyer receiving nothing pays 000.

Copies. The instance Icopies\mathcal I^{\mathrm{copies}}Icopies replaces each buyer iii by ∣Ji∣|J_i|∣Ji​∣ single-parameter agents, one per service j∈Jij \in J_ij∈Ji​ with value vjv_jvj​, under the same J\mathcal JJ.

Price menus. Given prices (pj)(p_j)(pj​) and an arrival order σ\sigmaσ, the price-menu mechanism approaches the buyers in order; buyer iii is offered the services of JiJ_iJi​ that can still be feasibly allocated, at prices pjp_jpj​, and buys a utility-maximizing one if some has pj≤vjp_j \le v_jpj​≤vj​.

Multiple copies of items. With items KKK and cap(k)\mathrm{cap}(k)cap(k) copies of item kkk, services are pairs (i,k)(i,k)(i,k) and a set of services is feasible if it gives each buyer at most one item and uses at most cap(k)\mathrm{cap}(k)cap(k) copies of kkk: the intersection of two partition matroids.

Formalization targets

Goal: Theorem 14

For regular distributions there are prices ppp such that, for every arrival order σ\sigmaσ, the price-menu mechanism Pσ\mathcal P_\sigmaPσ​ is truthful and

RA≤274 RPσ\mathcal R^{\mathcal A} \le \tfrac{27}{4}\,\mathcal R^{\mathcal P_\sigma}RA≤427​RPσ​

for every individually rational, truthful deterministic mechanism A\mathcal AA.

Milestones

  • Truthful BMUMD mechanisms are weakly monotone (p. 13), and the allocation of Acopies\mathcal A^{\mathrm{copies}}Acopies is monotone in each vjv_jvj​ (p. 13).
  • Lemma 3: RA≤RA′\mathcal R^{\mathcal A} \le \mathcal R^{\mathcal A'}RA≤RA′ for some truthful A′\mathcal A'A′ on Icopies\mathcal I^{\mathrm{copies}}Icopies.
  • The price-menu mechanism allocates a maximal feasible set of services (p. 14).
  • Theorem 4: if RM′≤α Rpobl\mathcal R^{M'} \le \alpha\,\mathcal R^{\mathrm{obl}}_{\mathbf p}RM′≤αRpobl​ for every truthful M′M'M′ on Icopies\mathcal I^{\mathrm{copies}}Icopies, then RA≤α RPσ\mathcal R^{\mathcal A} \le \alpha\,\mathcal R^{\mathcal P_\sigma}RA≤αRPσ​ for every σ\sigmaσ and every truthful IR A\mathcal AA.
  • Lemma 2 (regular part): RM≤∑ipiMqiM\mathcal R^M \le \sum_i p^M_i q^M_iRM≤∑i​piM​qiM​, with qiMq^M_iqiM​ the probability that MMM serves iii and Fi(piM)=1−qiMF_i(p^M_i) = 1 - q^M_iFi​(piM​)=1−qiM​.
  • Theorem 19 (existence form): a revenue-optimal truthful mechanism exists.
  • The claim ci≥4/9c_i \ge 4/9ci​≥4/9 of App. D.4: under ∑i′∈Pqi′≤cap(P)/3\sum_{i' \in P} q_{i'} \le \mathrm{cap}(P)/3∑i′∈P​qi′​≤cap(P)/3 in every part, with probability at least 4/94/94/9 neither part of iii is full without iii.
  • Theorem 13: for two partition matroids there are prices with RM≤274 Rpobl\mathcal R^M \le \tfrac{27}{4}\,\mathcal R^{\mathrm{obl}}_{\mathbf p}RM≤427​Rpobl​ for every truthful MMM.

Significance

The result shows that for unit-demand buyers, a seller loses at most a constant factor by replacing the optimal, possibly opaque, truthful mechanism with a menu of prices that does not depend on the order in which buyers arrive. The reduction of Theorem 4 is generic: any order-oblivious pricing for the single-parameter instance with copies, under any unit-demand constraint, transfers to the multi-parameter instance with the same factor. Theorem 13 supplies one such pricing for the intersection of two partition matroids, which is exactly the shape of the multi-unit, multi-item constraint.

All results here are proved in the paper and none is formalized elsewhere; the platform has Myerson's single-unit optimal auction and weak monotonicity in an abstract quasilinear model (Börgers), but no posted-price approximation, no copies reduction, and no order-oblivious revenue. The formal development adds a machine-checked account of the reduction (in particular that the price-menu mechanism is truthful and allocates a maximal feasible set for every order), a precise version of the probabilistic claim behind the constant 6.756.756.75, and reusable definitions of order-oblivious revenue and of multi-parameter truthfulness with the paper's individual rationality.

Difficulty

Lemma 3 needs more than the observation that the copies instance has more competition: one must build a truthful single-parameter mechanism with at least the same revenue. The allocation is copied, but the payments must be threshold payments of the copies mechanism, and showing they dominate the original payments uses both weak monotonicity and the paper's individual rationality, through the taxation principle.

Theorem 13 compares order-oblivious revenue with Myerson's revenue through the bound of Lemma 2, at prices built from Myerson's service probabilities scaled by 1/31/31/3. The step that is easy to get wrong is the probability that an agent is considered: the events "part P1P_1P1​ is not full" and "part P2P_2P2​ is not full" depend on overlapping agents, so the product bound (2/3)(2/3)(2/3)(2/3)(2/3)(2/3) does not follow from Markov's inequality alone; it holds because both events are decreasing in the set of desiring agents (Harris' inequality). The comparison must also be uniform: one set of prices must serve against every truthful mechanism, which requires an optimal mechanism to exist.

Formalization scope

  • Distributions (P1): each FjF_jFj​ has a measurable density, strictly positive on a bounded interval [v‾j,v‾j]⊆[0,∞)[\underline v_j, \overline v_j] \subseteq [0, \infty)[v​j​,vj​]⊆[0,∞), with no mass outside. Values are independent (product prior).
  • Regularity (P2): the virtual value ϕ(v)=v−(1−F(v))/f(v)\phi(v) = v - (1 - F(v))/f(v)ϕ(v)=v−(1−F(v))/f(v) is non-decreasing on the support. It is assumed in Lemma 2, Theorem 19, Theorem 13 and the goal. Theorem 14 does not state it, but its proof goes through Theorem 13, which the paper proves for regular distributions; the non-regular extension (App. E, randomized prices) is out of scope, as is the second paragraph of Lemma 2.
  • Mechanisms (P3): deterministic; dominant-strategy truthful with misreports in the support (a buyer misreports all coordinates of JiJ_iJi​ at once); single-parameter IR is ex-post nonnegative utility; multi-parameter IR is the paper's (πi≤vj\pi_i \le v_jπi​≤vj​ if served jjj, πi=0\pi_i = 0πi​=0 if unserved); allocation events and payments measurable, payments integrable.
  • Benchmarks (P4): Myerson's mechanism is not constructed. "Approximates RM\mathcal R^{\mathcal M}RM" is stated against every truthful mechanism, and Lemma 3 and Theorem 19 in existence form.
  • Price menus: ties between utility-maximizing services are broken by a fixed enumeration of JJJ; a service of utility 000 is bought. Theorem 4 assumes α≥0\alpha \ge 0α≥0.
  • Dropped: the last sentence of Theorem 14 (polynomial-time computability of the prices) has no cost model here.
  • Constant: 6.756.756.75 is written 27/427/427/4 everywhere.
  • Not trivializable: the prices in Theorem 13 and the goal are chosen before the mechanism, and the benchmark includes every truthful mechanism, so a degenerate price vector cannot meet the bound; Rpobl\mathcal R^{\mathrm{obl}}_{\mathbf p}Rpobl​ is a genuine minimum over a nonempty finite class.

Needed infrastructure, reusable beyond this mission: Myerson's characterization of truthful single-parameter mechanisms and the revenue–virtual-surplus identity for densities on intervals, Harris' inequality for product measures, and the taxation principle for deterministic multi-parameter mechanisms. Contributions on any of these are welcome.

Selected references

  • S. Chawla, J. D. Hartline, D. Malec, B. Sivan, Multi-parameter Mechanism Design and Sequential Posted Pricing, STOC 2010; arXiv:0907.2435v2, 2010. https://arxiv.org/abs/0907.2435
  • R. Myerson, Optimal Auction Design, Mathematics of Operations Research 6(1), 1981. https://doi.org/10.1287/moor.6.1.58
  • S. Chawla, J. D. Hartline, R. Kleinberg, Algorithmic Pricing via Virtual Valuations, EC 2007. https://arxiv.org/abs/0711.3203
  • A. M. Manelli, D. R. Vincent, Multidimensional mechanism design: Revenue maximization and the multiple-good monopoly, Journal of Economic Theory 137(1), 2007. https://doi.org/10.1016/j.jet.2006.12.007
  • T. E. Harris, A lower bound for the critical probability in a certain percolation process, Proc. Cambridge Philos. Soc. 56, 1960. https://doi.org/10.1017/S0305004100034241
15 thms2 active usersReviewed
Operations ResearchProbabilityStochastic Systems·Captain: mikedeng1

Open Queueing Networks in Heavy Traffic: Reflected Brownian Motion Limit for the Queue Length ProcessResearch Paper

Motivation

Open networks of single-server queues with general interarrival and service distributions are the standard model of job shops, communication networks and service systems. Outside the product-form (Jackson) case their queue-length distributions are not known in closed form. When every station is close to saturation, a heavy-traffic limit replaces the network by a diffusion process. Martin I. Reiman's paper Open Queueing Networks in Heavy Traffic (Mathematics of Operations Research 9(3), 1984) proves such a limit for the vector of queue lengths of a general open network. The limit is a reflected Brownian motion on the nonnegative orthant. That process has since become the default diffusion approximation for open networks, and it is the starting point of later work on its stationary distribution and on control of networks in heavy traffic.

Timeline:

  • Iglehart and Whitt (1970a,b) proved heavy-traffic limits for a single multiple-server station and for acyclic networks, in which no customer visits a station twice.
  • Harrison (1973, 1978) treated tandem queues; the 1978 paper introduced reflected Brownian motion on the nonnegative orthant as the diffusion limit.
  • Harrison and Reiman (1981a, Ann. Probab. 9:302–308) constructed reflected Brownian motion on the orthant through a continuous reflection mapping. That paper is the source of Lemma 1 here.

(These attributions follow Reiman's own account, pp. 441–442 of the 1984 paper.)

  • Reiman (1984) proved the limit for general open networks with Markovian routing (Theorem 1). The paper also proves a limit for sojourn times along fixed routes (Theorem 2).

Setting

There are KKK single-server stations and a nonempty set J⊆{1,…,K}\mathcal J\subseteq\{1,\dots,K\}J⊆{1,…,K} of stations that receive customers from outside. The primitives are mutually independent sequences of IID random variables: interarrival times uki>0u_k^i>0uki​>0 (k∈Jk\in\mathcal Jk∈J), service times vki>0v_k^i>0vki​>0, and routing indicators ϕki∈{0,1,…,K}\phi_k^i\in\{0,1,\dots,K\}ϕki​∈{0,1,…,K}. When the iiith customer served at station kkk finishes, it moves to station ϕki\phi_k^iϕki​, or leaves if ϕki=0\phi_k^i=0ϕki​=0. The parameters are the service rates μk=(Evk1)−1\mu_k=(E v_k^1)^{-1}μk​=(Evk1​)−1, the service-time variances sk=var⁡vk1s_k=\operatorname{var} v_k^1sk​=varvk1​, the arrival rates λk=(Euk1)−1\lambda_k=(E u_k^1)^{-1}λk​=(Euk1​)−1 (with λk=0\lambda_k=0λk​=0 for k∉Jk\notin\mathcal Jk∈/J), and the interarrival variances ak=var⁡uk1a_k=\operatorname{var} u_k^1ak​=varuk1​. The routing matrix P=(pkj)P=(p_{kj})P=(pkj​), pkj=P{ϕk1=j}p_{kj}=P\{\phi_k^1=j\}pkj​=P{ϕk1​=j}, has spectral radius strictly less than one, so every customer eventually leaves.

Let Ak(t)A_k(t)Ak​(t) be the number of exogenous arrivals to station kkk by time ttt, and Sk(t)S_k(t)Sk​(t) the number of service completions at kkk in ttt units of busy time. Let S^k(t)=∑i≤Sk(t)eϕki−Sk(t)ek\hat S_k(t)=\sum_{i\le S_k(t)}e_{\phi_k^i}-S_k(t)e_kS^k​(t)=∑i≤Sk​(t)​eϕki​​−Sk​(t)ek​, with e0=0e_0=0e0​=0. The queue length Q(t)∈Z+KQ(t)\in\mathbb Z_+^KQ(t)∈Z+K​ and the busy time B(t)B(t)B(t) are the unique solution of

Q(t)=A(t)+∑k=1KS^k(Bk(t)),Bk(t)=∫0t1{Qk(s)>0} ds,B(0)=0.Q(t)=A(t)+\sum_{k=1}^K\hat S_k(B_k(t)),\qquad B_k(t)=\int_0^t1_{\{Q_k(s)>0\}}\,ds,\qquad B(0)=0 .Q(t)=A(t)+k=1∑K​S^k​(Bk​(t)),Bk​(t)=∫0t​1{Qk​(s)>0}​ds,B(0)=0.

A sequence of such networks, indexed by nnn, shares KKK, J\mathcal JJ and PPP. Its parameters μ(n),s(n),λ(n),a(n)\mu(n),s(n),\lambda(n),a(n)μ(n),s(n),λ(n),a(n) converge to finite limits μ,s,λ,a\mu,s,\lambda,aμ,s,λ,a. With ν(n)=λ(n)+μ(n)P\nu(n)=\lambda(n)+\mu(n)Pν(n)=λ(n)+μ(n)P, the heavy-traffic condition is

ck(n)=n (νk(n)−μk(n))→ck.c_k(n)=\sqrt n\,(\nu_k(n)-\mu_k(n))\to c_k .ck​(n)=n​(νk​(n)−μk​(n))→ck​.

Moments of order 2+ϵ2+\epsilon2+ϵ of the interarrival and service times are bounded uniformly in nnn. The scaled queue length is Zn(t)=n−1/2Qn(nt)Z^n(t)=n^{-1/2}Q^n(nt)Zn(t)=n−1/2Qn(nt), 0≤t≤10\le t\le10≤t≤1.

Formalization targets

Goal: Theorem 1

Let ξ\xiξ be a Brownian motion with drift ccc and covariance matrix A\mathcal AA, where

Aii=λi3ai+μi3si(1−2pii)+∑jμjpji(1−pji+pjiμj2sj),\mathcal A_{ii}=\lambda_i^3a_i+\mu_i^3s_i(1-2p_{ii})+\sum_j\mu_jp_{ji}(1-p_{ji}+p_{ji}\mu_j^2s_j),Aii​=λi3​ai​+μi3​si​(1−2pii​)+j∑​μj​pji​(1−pji​+pji​μj2​sj​), Aij=−[μi3sipij+μj3sjpji+∑kμkpkipkj(1−μk2sk)](i≠j).\mathcal A_{ij}=-\Big[\mu_i^3s_ip_{ij}+\mu_j^3s_jp_{ji}+\sum_k\mu_kp_{ki}p_{kj}(1-\mu_k^2s_k)\Big]\quad(i\ne j).Aij​=−[μi3​si​pij​+μj3​sj​pji​+k∑​μk​pki​pkj​(1−μk2​sk​)](i=j).

Let Z=ϕ(ξ)Z=\phi(\xi)Z=ϕ(ξ) be its reflection with reflection matrix I−PI-PI−P. Then

Zn⇒Zin D[0,1] (Skorohod topology).Z^n\Rightarrow Z\quad\text{in } D[0,1]\text{ (Skorohod topology)}.Zn⇒Zin D[0,1] (Skorohod topology).

The goal fixes no constants beyond the parameters' limits. It is stated for every network sequence satisfying (20)–(26).

Milestones

The milestones follow the paper's proof, in order:

  • the existence and uniqueness claim for (1)–(3);
  • the representation Q=X~+Y(I−P)Q=\tilde X+Y(I-P)Q=X~+Y(I−P) (Eq. (13));
  • the least-element map fff (Proposition 1);
  • the reflection mapping ϕ\phiϕ (Lemma 1) and f=ϕf=\phif=ϕ on continuous paths (Proposition 2);
  • the netput limit ζn⇒ζ\zeta^n\Rightarrow\zetaζn⇒ζ (Proposition 3);
  • stochastic boundedness of ZnZ^nZn (Lemma 6);
  • vanishing scaled idleness n−1Ikn(n)→0n^{-1}I^n_k(n)\to0n−1Ikn​(n)→0 (Proposition 4);
  • the centred limit ζ~n⇒ζ\tilde\zeta^n\Rightarrow\zetaζ~​n⇒ζ (Proposition 5).

Significance

Theorem 1 justifies the diffusion approximation of a heavily loaded open network. Writing Qn(t)≈n Z(t/n)Q^n(t)\approx\sqrt n\,Z(t/n)Qn(t)≈n​Z(t/n) reduces questions about the network to questions about one reflected Brownian motion, whose data are explicit functions of the first two moments of the primitives and of the routing matrix. The same limit, with Lemma 2, gives the paper's Theorem 2 on sojourn times. It is the model case for the multiclass heavy-traffic theory that followed.

The result has been proved since 1984. No machine-checked version exists. The mission's contributions would be:

  • a formal statement of the network, of its Harrison representation, and of weak convergence in DDD;
  • a formal proof of the reflection-mapping facts (Proposition 1, Lemma 1, Proposition 2), which are deterministic and reusable;
  • eventually, a formal proof of the full limit theorem.

Difficulty

The obvious route applies a functional central limit theorem to QnQ^nQn directly. That fails because QnQ^nQn is not a sum of independent terms: each station serves only while its queue is nonempty, so the service process is evaluated at the random busy time Bk(t)B_k(t)Bk​(t), which depends on the whole network. The proof therefore has to separate the netput process, which obeys a central limit theorem, from the regulator YYY. It then has to show that the random time change Bkn(nt)/nB^n_k(nt)/nBkn​(nt)/n converges to the identity, i.e. that idleness vanishes on the diffusion scale. Weak convergence must also be transported through a reflection map that is defined on all of DDD but is known to be continuous only at continuous paths.

Formalization scope

The Lean development uses the following conventions:

  • Stations are Fin K, vectors are row vectors Fin K → ℝ, and a row vector times a matrix is Matrix.vecMul.
  • A routing indicator lives in Fin (K+1), with 0 meaning "leaves" and j.succ meaning station jjj.
  • The primitives are mutually independent (iIndep of their σ-algebras), IID within each sequence, everywhere positive and square integrable.
  • "Spectral radius <1<1<1" is stated as Pm→0P^m\to0Pm→0.
  • (Qn,Bn)(Q^n,B^n)(Qn,Bn) is any pair solving (1)–(3) almost surely, with measurable paths so that (2) is a Lebesgue integral.
  • The networks are indexed by ℕ; (25)–(26) are imposed for n≥1n\ge1n≥1, (22) and (26) over k∈Jk\in\mathcal Jk∈J, and J\mathcal JJ is the same for all nnn.
  • Brownian motion with drift ccc and covariance A\mathcal AA lives on [0,∞)[0,\infty)[0,∞). It is defined by continuity, ξ(0)=0\xi(0)=0ξ(0)=0, independent increments, and the Gaussian characteristic function of increments.
  • ZZZ is the reflection of ξ\xiξ in the sense of (14)–(17).
  • Weak convergence in DDD is stated in Skorohod-representation form: a coupling with almost-sure J1_11​ convergence on [0,1][0,1][0,1]. This form accommodates a separate probability space for each nnn.

Added hypotheses, each implicit on the page:

  1. The existence item assumes Uk(l),Vk(l)→∞U_k(l),V_k(l)\to\inftyUk​(l),Vk​(l)→∞ at the sample point; without it the maxima defining Ak(t)A_k(t)Ak​(t) and Sk(t)S_k(t)Sk​(t) need not exist.
  2. Solutions of (1)–(3) have measurable paths.

No positivity hypothesis on the limits μk\mu_kμk​ is added: (25) and (26) bound the means of the service and interarrival times, so the limits are positive.

The statement is not to be weakened. Ruled out are:

  • convergence of finite-dimensional distributions only;
  • a single network without the index nnn;
  • uniform convergence used in place of the Skorohod topology without the coupling;
  • a Brownian motion that is not required to have independent Gaussian increments.

Each of these is a different theorem.

Useful contributions, all reusable beyond this mission:

  • the deterministic reflection-map results;
  • Donsker-type theorems for renewal counting processes in DDD;
  • the random time-change lemma (Billingsley);
  • the continuous mapping theorem in coupling form.

Selected references

  • M. I. Reiman, Open Queueing Networks in Heavy Traffic, Mathematics of Operations Research 9(3):441–458, 1984. https://doi.org/10.1287/moor.9.3.441
  • J. M. Harrison and M. I. Reiman, Reflected Brownian Motion on an Orthant, Annals of Probability 9:302–308, 1981 (cited in Reiman 1984 as [6]).
  • J. M. Harrison, The Diffusion Approximation for Tandem Queues in Heavy Traffic, Advances in Applied Probability 10:886–905, 1978 (Reiman 1984, [5]).
  • J. M. Harrison, The Heavy Traffic Approximation for Single Server Queues in Series, Journal of Applied Probability 10:613–629, 1973 (Reiman 1984, [4]).
  • D. L. Iglehart and W. Whitt, Multiple Channel Queues in Heavy Traffic, I and II: Sequences, Networks, and Batches, Advances in Applied Probability 2:150–177 and 355–364, 1970 (Reiman 1984, [8], [9]).
  • P. Billingsley, Convergence of Probability Measures, Wiley, New York, 1968 (Reiman 1984, [1]).
13 thms2 active usersReviewed
CombinatoricsGraph TheoryOperations Research·Captain: mikedeng1

The Strong Perfect Graph Theorem I: A Graph Is Perfect If and Only If It Is BergeResearch Paper

Motivation

A perfect graph is one whose coloring problem has a particularly sharp answer on every induced subgraph: the fewest colors needed is exactly the size of its largest clique. This makes a local obstruction, a clique, certify the optimum number of colors throughout the graph. Claude Berge proposed in 1961 that perfection could be recognized by the absence of two kinds of induced odd cycles, one in the graph and one in its complement. The equivalence became known as the strong perfect graph conjecture. Chudnovsky, Robertson, Seymour, and Thomas proved it in their 2006 paper, which also proves a structural decomposition of the graphs under study. The paper connects this question to graph coloring, Shannon capacity, and linear and integer programming. Chudnovsky et al., pp. 51–54

The earlier complement theorem was proved by Lovász in 1972 and appears as Theorem 1.1 in the paper. The strong conjecture remained unresolved for roughly four decades; the authors' Theorem 1.2 settles it. Their proof places a second result, Theorem 1.3, beside the equivalence: a graph with no forbidden odd hole or antihole must either belong to a basic class or admit one of several specified decompositions. The graph classes and decompositions are therefore part of the statement of the route to the main result, not merely vocabulary for a proof. Chudnovsky et al., pp. 52–56

Setting

All graphs here are finite and simple. The complement G‾\overline GG has the same vertices as GGG, and two distinct vertices are adjacent in G‾\overline GG exactly when they are not adjacent in GGG. For a vertex set XXX, the notation G∣XG|XG∣X means the induced subgraph on XXX. A clique is a set of pairwise adjacent vertices. Its largest possible size in a graph HHH is ω(H)\omega(H)ω(H), and χ(H)\chi(H)χ(H) is the minimum number of colors in a proper vertex coloring of HHH.

A hole is an induced cycle of length at least four. An antihole of GGG is a hole in G‾\overline GG. A graph is Berge if every hole and antihole has even length. Thus a perfect graph requires χ(G∣X)=ω(G∣X)\chi(G|X)=\omega(G|X)χ(G∣X)=ω(G∣X) for every X⊆V(G)X\subseteq V(G)X⊆V(G), while a Berge graph satisfies a restriction on induced cycles in both GGG and G‾\overline GG. “Induced” matters: a cycle with a chord is not a hole. Chudnovsky et al., pp. 51–52

For the structural milestones, a basic graph is a bipartite graph, the complement of one, a line graph of a bipartite graph, the complement of such a line graph, or a double split graph. The latter consists of paired vertices ai,bia_i,b_iai​,bi​ and cj,djc_j,d_jcj​,dj​ with the within-pair and between-pair adjacencies specified on pp. 52–53. A proper 2-join partitions the vertices into two sides with two prescribed complete cross-edge blocks, connected-component conditions on both sides, and a special odd-path condition. A proper homogeneous pair is a pair of vertex sets whose outside vertices split into four nonempty adjacency classes. A balanced skew partition has one side disconnected and the other disconnected in the complement, together with parity restrictions on induced paths and antipaths. Chudnovsky et al., pp. 52–54

Formalization targets

Theorem 1.2: perfection and the Berge property

The goal is the exact equivalence for every finite simple graph:

G is perfect⟺G is Berge.G\text{ is perfect}\quad\Longleftrightarrow\quad G\text{ is Berge}.G is perfect⟺G is Berge.

No order bound, chosen graph class, or decomposition hypothesis is attached to the goal. Chudnovsky et al., p. 52, 1.2

Structural and reduction milestones

Theorem 1.1 says that GGG perfect implies G‾\overline GG perfect. Theorem 1.5 says that a minimum imperfect graph, a Berge nonperfect graph with the smallest vertex count among all such graphs, cannot admit a balanced skew partition. Theorem 13.5 says that a recalcitrant graph—a Berge graph with the listed line-graph, double-split, 2-join, homogeneous-pair, and balanced-skew outcomes absent—has GGG or G‾\overline GG bipartite. Theorem 1.3 states the decomposition conclusion:

G Berge⟹G basic ∨ G or G‾ has a proper 2-join ∨ G has a proper homogeneous pair ∨ G has a balanced skew partition.G\text{ Berge}\Longrightarrow G\text{ basic}\ \lor\ G\text{ or }\overline G\text{ has a proper 2-join}\ \lor\ G\text{ has a proper homogeneous pair}\ \lor\ G\text{ has a balanced skew partition}.G Berge⟹G basic ∨ G or G has a proper 2-join ∨ G has a proper homogeneous pair ∨ G has a balanced skew partition.

The milestone order records the two reduction results, the later structural capstone, and the decomposition statement it yields. The paper's other section results that establish 13.5 are posed in the remaining missions of this series. Chudnovsky et al., pp. 52, 54–55, 154

Significance

Theorem 1.2 gives a forbidden-induced-subgraph characterization of perfect graphs. Its cycle condition is intrinsic to the graph and its complement; its coloring condition quantifies over every induced subgraph. Together with Theorem 1.3, it ties a numerical property of colorings to explicit graph structures and separations. The complement theorem and the exclusion of decompositions for a minimum imperfect graph explain why the structural alternatives have the strength needed for the equivalence. Chudnovsky et al., pp. 52–56

The mathematical theorem was proved in the cited paper. This mission poses its statements in Lean and seeks machine-checked proofs; the draft theorem declarations are open targets. The definition layer is useful beyond this mission: induced holes, Berge graphs, perfection, balanced skew partitions, and the decomposition predicates can support the later missions without changing what each source statement means. No machine-checked proof of these draft targets is claimed here.

Difficulty

The forward implication can be tested on induced odd cycles, but that observation does not settle the converse. Excluding odd holes and antiholes does not give an immediate coloring of an arbitrary induced subgraph. The paper instead establishes a detailed account of what a Berge graph can look like when it is not in a basic class. The delicate point in turning this account into Theorem 1.2 is that each decomposition outcome must be incompatible with a minimum imperfect graph. Ordinary skew partitions are too broad for that role; the balanced parity conditions are part of the statement. The structural conclusion 13.5 collects restrictions established across many later sections, so formalizing its prerequisites is a substantial graph-theoretic task. Chudnovsky et al., pp. 52–56, 154

Formalization scope

Lean uses SimpleGraph V with finite vertices, decidable vertex equality, and the Mathlib complement, induced subgraph, chromatic number, clique number, bipartiteness, and line graph. A hole is a list in cyclic order whose adjacency relation agrees exactly with the cycle edges; a path is likewise listed in one orientation with exactly its consecutive edges. Antiholes and antipaths use the complement graph. The empty vertex set is connected, matching p. 54. The double split partition is encoded by an equivalence from the four indexed parts to the whole vertex type, making disjointness and coverage explicit. A minimum imperfect graph is globally minimal by vertex count, among all finite Berge graphs.

No hypothesis beyond the paper's finite, simple graph convention is added to the goal or numbered milestones. In particular, “Berge” includes holes in both GGG and G‾\overline GG; “perfect” ranges over every induced subgraph; and the complement occurs only in those decomposition outcomes where the paper places it. A non-induced cycle or a missing complement condition would make the formal target different. The definitions and Lean proofs of the graph classes, their boundary cases, and the structural milestones are welcome contributions. Later missions pose the paper's intervening numbered results rather than duplicating them here.

Selected references

  • Maria Chudnovsky, Neil Robertson, Paul Seymour, and Robin Thomas, The strong perfect graph theorem, Annals of Mathematics 164 (2006), 51–229. DOI 10.4007/annals.2006.164.51.
16 thms3 active usersReviewed
Discrete GeometryLinear OptimizationOperations Research+2·Captain: mikedeng1

Smoothed Analysis of Algorithms: Why the Simplex Algorithm Usually Takes Polynomial Time 1: The Expected Shadow of a Gaussian-Perturbed Polytope Has Polynomially Many VerticesResearch Paper

Why the shadow of a perturbed polytope matters

The simplex method solves linear programs very fast in practice, yet for most pivot rules there are inputs on which it takes exponentially many steps (Klee and Minty, 1972, for Dantzig's rule; Goldfarb, 1983, for the shadow-vertex rule). Average-case analyses (Borgwardt, 1980s; Smale, 1983) explained good behaviour on random inputs, but random inputs look nothing like real ones. Spielman and Teng introduced smoothed analysis to close this gap: the input is chosen by an adversary and then perturbed by a small Gaussian, and the running time is measured in expectation over the perturbation. They proved that the shadow-vertex simplex method has smoothed complexity polynomial in the number of constraints nnn, the dimension ddd and 1/σ1/\sigma1/σ (Spielman–Teng, J. ACM 2004; this mission follows the preprint arXiv:cs/0111050v7). The work received the Gödel Prize (2008) and the Fulkerson Prize (2009).

Timeline. Borgwardt (1977–1987) bounded the expected number of shadow-vertex pivots for rotationally symmetric random data. Spielman and Teng (2001, STOC; journal 2004) proved the first smoothed bound, with a shadow bound of order nd3/σ6nd^3/\sigma^6nd3/σ6 — the theorem of this mission. Deshpande and Spielman (FOCS 2005) improved the shadow bound, Vershynin (2009) reduced the dependence on nnn to polylogarithmic, and Dadush and Huiberts (STOC 2018) obtained O(d2log⁡n σ−2)O(d^2\sqrt{\log n}\,\sigma^{-2})O(d2logn​σ−2) for small σ\sigmaσ.

Setting

Fix d≥3d\ge3d≥3 and n>dn>dn>d. The data are vectors a1,…,an∈Rda_1,\dots,a_n\in\mathbb R^da1​,…,an​∈Rd, the constraint vectors of the linear program max⁡⟨z∣x⟩\max\langle z|x\ranglemax⟨z∣x⟩ subject to ⟨ai∣x⟩≤1\langle a_i|x\rangle\le1⟨ai​∣x⟩≤1 for all iii. Each aia_iai​ is a Gaussian of standard deviation σ\sigmaσ centered at a point aˉi\bar a_iaˉi​ with ∥aˉi∥≤1\|\bar a_i\|\le1∥aˉi​∥≤1: it has density

μi(a)=(12π σ)de−∥a−aˉi∥2/2σ2,\mu_i(a)=\Big(\tfrac{1}{\sqrt{2\pi}\,\sigma}\Big)^d e^{-\|a-\bar a_i\|^2/2\sigma^2},μi​(a)=(2π​σ1​)de−∥a−aˉi​∥2/2σ2,

and the aia_iai​ are independent (joint density ∏iμi(ai)\prod_i\mu_i(a_i)∏i​μi​(ai​)).

For a direction q∈Rdq\in\mathbb R^dq∈Rd, optSimpq(a1,…,an)\mathrm{optSimp}_q(a_1,\dots,a_n)optSimpq​(a1​,…,an​) is the set of index sets I⊆{1,…,n}I\subseteq\{1,\dots,n\}I⊆{1,…,n} with ∣I∣=d|I|=d∣I∣=d such that (ai)i∈I(a_i)_{i\in I}(ai​)i∈I​ is linearly independent, the simplex △(AI)=ConvHull(ai:i∈I)\triangle(A_I)=\mathrm{ConvHull}(a_i:i\in I)△(AI​)=ConvHull(ai​:i∈I) is a facet of ConvHull(0,a1,…,an)\mathrm{ConvHull}(0,a_1,\dots,a_n)ConvHull(0,a1​,…,an​), and qqq lies in the cone {∑i∈Iαiai:αi≥0}\{\sum_{i\in I}\alpha_ia_i:\alpha_i\ge0\}{∑i∈I​αi​ai​:αi​≥0}. In polar terms, III is the set of tight constraints at the vertex of the feasible polyhedron that maximizes ⟨q∣x⟩\langle q|x\rangle⟨q∣x⟩.

For linearly independent t,zt,zt,z, the shadow Shadowt,z(a1,…,an)\mathrm{Shadow}_{t,z}(a_1,\dots,a_n)Shadowt,z​(a1​,…,an​) is the set of index sets III that belong to optSimpq\mathrm{optSimp}_qoptSimpq​ for some nonzero q∈Span(t,z)q\in\mathrm{Span}(t,z)q∈Span(t,z). Its size is the number of vertices of the projection of the feasible polyhedron onto the plane Span(t,z)\mathrm{Span}(t,z)Span(t,z); the shadow-vertex method walks along this polygon, one pivot per vertex. Finally

D(n,d,σ)=58,888,678 nd3min⁡(σ, 1/(3dln⁡n))6.\mathcal D(n,d,\sigma)=\frac{58{,}888{,}678\,nd^3}{\min\big(\sigma,\,1/(3\sqrt{d\ln n})\big)^6}.D(n,d,σ)=min(σ,1/(3dlnn​))658,888,678nd3​.

Formalization targets

Goal: Theorem 4.0.1 (Shadow Size)

Ea1,…,an[ ∣Shadowt,z(a1,…,an)∣ ]≤D(n,d,σ)\mathbb E_{a_1,\dots,a_n}\big[\,|\mathrm{Shadow}_{t,z}(a_1,\dots,a_n)|\,\big]\le\mathcal D(n,d,\sigma)Ea1​,…,an​​[∣Shadowt,z​(a1​,…,an​)∣]≤D(n,d,σ)

for every d≥3d\ge3d≥3, n>dn>dn>d, every pair of linearly independent t,zt,zt,z, every σ>0\sigma>0σ>0 and all centers of norm at most 111.

Milestones

The milestones follow the paper's proof, leaves first.

  • Probability tools: the chi-square bound (Corollary 2.4.6), the combination lemma (Lemma 2.3.5), almost polynomial densities (Lemma 2.3.7), and comparing Gaussian tails (Lemma 2.4.11).
  • Reduction: the measure of the event P={∥ai∥≤2 ∀i}P=\{\|a_i\|\le2\ \forall i\}P={∥ai​∥≤2 ∀i} (Proposition 4.0.5), and the discretization of the shadow into mmm equally spaced directions (Lemma 4.0.6).
  • Angle bound: the probability, conditioned on PPP, that the ray through a fixed unit vector qqq passes within angle ε\varepsilonε of the boundary of its optimal facet is O(nd3ε/σ6)O(nd^3\varepsilon/\sigma^6)O(nd3ε/σ6) (Lemma 4.0.7, from Lemma 4.0.11).
  • Distance and incidence: in Blaschke coordinates ai=Rωbi+sqa_i=R_\omega b_i+sqai​=Rω​bi​+sq, a deterministic split (Lemma 4.0.12), a distance bound (Lemmas 4.1.1–4.1.3) and an angle-of-incidence bound (Lemmas 4.2.1–4.2.3).

Significance

The result. Theorem 4.0.1 is the geometric heart of the smoothed analysis of the simplex method. Section 4.3 of the paper extends it to arbitrary centers, covariances and right-hand sides, and Section 5 combines these extensions with a two-phase method to show that the simplex method has polynomial smoothed complexity. The same shadow bound underlies later analyses of the simplex method, of perturbed polytopes' diameters, and of condition numbers of random linear programs.

Formalizing it. The theorem has been proved, and improved constants are known, but none of this is machine-checked. A formal proof would verify a long and delicate argument: a change of variables of integral geometry (Blaschke's formula), several conditional-density estimates, and explicit constants in the millions. The mission also produces reusable statements about Gaussian vectors and convex hulls of random points.

Difficulty

The obvious approach is to count, for each candidate facet III, the probability that III appears in the shadow; there are (nd)\binom nd(dn​) candidates, so a union bound is exponential in ddd. The paper avoids this by discretizing the angle of qqq (Lemma 4.0.6) and bounding, for each fixed direction, the probability that the optimal facet changes within a small angular step. That needs a lower bound on the angle between qqq and the boundary of its optimal facet, conditioned on the facet being optimal. The conditioning changes the distribution of a1,…,ada_1,\dots,a_da1​,…,ad​, so the bound cannot come from the Gaussian density alone. The proof changes variables to the facet's normal ω\omegaω, offset sss and in-plane coordinates bib_ibi​ (Corollary 2.5.3), whose Jacobian contributes the factors ⟨ω∣q⟩\langle\omega|q\rangle⟨ω∣q⟩ and Vol(△(b))\mathrm{Vol}(\triangle(b))Vol(△(b)). It then shows that both the distance of the origin to a face of the in-plane simplex and the angle of incidence ⟨ω∣q⟩\langle\omega|q\rangle⟨ω∣q⟩ are unlikely to be small. Measure-theoretic bookkeeping is as hard as the geometry: densities known only up to normalization, conditioning on events of positive measure, and the measure-zero degeneracies the paper sets aside.

Formalization scope

Points live in EuclideanSpace ℝ (Fin d). Constraint vectors are indexed by Fin n (0-based), so the paper's {1,…,d}\{1,\dots,d\}{1,…,d} is {i:i<d}\{i:i<d\}{i:i<d}. The Gaussian of standard deviation σ\sigmaσ centered at ccc is Lebesgue measure with the density above, and the joint law is the product measure. Lemma 4.0.6 also uses Mathlib's multivariateGaussian with a positive definite covariance. Expectations of shadow sizes are lower Lebesgue integrals of [0,∞][0,\infty][0,∞]-valued counts, and their measurability is part of each conclusion. "Density proportional to ν\nuν" and conditional probabilities are stated cross-multiplied, ∫Eν≤bound⋅∫ν\int_{E}\nu\le\text{bound}\cdot\int\nu∫E​ν≤bound⋅∫ν, so no 0/00/00/0 appears.

The shadow is the set of index sets III, and the direction q=0q=0q=0 is excluded. Including it would add every facet of ConvHull(0,a1,…,an)\mathrm{ConvHull}(0,a_1,\dots,a_n)ConvHull(0,a1​,…,an​) to the shadow, since 000 lies in every cone, and make the goal false. ang(q,∅)=∞\mathrm{ang}(q,\emptyset)=\inftyang(q,∅)=∞ is represented exactly in [0,∞][0,\infty][0,∞], never by a real infimum. Where the paper omits a hypothesis it uses, it is added and recorded in the item: the standing assumptions d≥3d\ge3d≥3, n>dn>dn>d and σ≤1/(3dln⁡n)\sigma\le1/(3\sqrt{d\ln n})σ≤1/(3dlnn​) (Lemma 4.2.3 is false without a bound on σ\sigmaσ), unit length of the reference vector qqq, s≥0s\ge0s≥0, and ε>0\varepsilon>0ε>0 for strict inequalities. Lemma 2.3.7 is stated with ≤\le≤ rather than the page's <<<, which fails in an edge case.

Infrastructure a complete development needs: Gaussian tail and chi-square estimates; faces and facets of convex hulls; the Blaschke change of variables and the latitude–longitude change of variables on the sphere (not in Mathlib); surface measure on Sd−1S^{d-1}Sd−1 (Mathlib's Measure.toSphere); and the disintegration of the joint law used in the combination lemma. The Gaussian estimates, the combination lemma and the Blaschke formula are useful beyond this mission. Proofs of any milestone, and of supporting lemmas such as the change-of-variables formulas, are welcome.

Selected references

  • D. A. Spielman, S.-H. Teng, Smoothed Analysis of Algorithms: Why the Simplex Algorithm Usually Takes Polynomial Time, arXiv:cs/0111050v7, 2003. https://arxiv.org/abs/cs/0111050v7
  • D. A. Spielman, S.-H. Teng, Smoothed analysis of algorithms: Why the simplex algorithm usually takes polynomial time, J. ACM 51(3):385–463, 2004. https://doi.org/10.1145/990308.990310
  • K. H. Borgwardt, The Simplex Method: A Probabilistic Analysis, Springer, 1987.
  • V. Klee, G. J. Minty, How good is the simplex algorithm?, in Inequalities III, Academic Press, 1972, 159–175.
  • A. Deshpande, D. A. Spielman, Improved smoothed analysis of the shadow vertex simplex method, FOCS 2005, 387–396.
  • R. Vershynin, Beyond Hirsch conjecture: walks on random polytopes and smoothed complexity of the simplex method, SIAM J. Comput. 39(2):646–678, 2009. https://doi.org/10.1137/070683386
  • D. Dadush, S. Huiberts, A friendly smoothed analysis of the simplex method, STOC 2018; arXiv:1711.05667. https://arxiv.org/abs/1711.05667
29 thms2 active usersReviewed
Operations ResearchOptimizationProbability+1·Captain: mikedeng1

Asymptotic Behavior of Statistical Estimators and of Optimal Solutions of Stochastic Optimization Problems: Optimal Solutions Under Estimated Distributions Are Strongly ConsistentResearch Paper

Motivation

Many estimation procedures in statistics, and most stochastic optimization models in operations research, have the same shape: a decision or parameter x∈Rnx\in\mathbb R^nx∈Rn is chosen to minimize an expected loss Ef(x)=∫f(x,ξ) P(dξ)Ef(x)=\int f(x,\xi)\,P(d\xi)Ef(x)=∫f(x,ξ)P(dξ) under a distribution PPP that is not known. In practice PPP is replaced by an estimate PνP^\nuPν built from the information available at stage ν\nuν (an empirical measure, a smoothed or parametric fit, a Bayesian posterior), and the minimizer of the estimated problem is used in place of the true one. The basic question is whether this is justified: do the estimated solutions converge to a true solution, and the estimated optimal values to the true optimal value, as information accumulates?

For maximum likelihood this is Wald's consistency theorem (Wald 1949); Huber extended it to M-estimators under non-standard conditions (Huber 1967). Both settings are unconstrained, or constrained to an open set, and assume finite-valued criteria. Constrained least squares, L1L^1L1 and Huber regression with inequality constraints, variance-component models with Heywood cases, and two-stage stochastic programs with recourse all lead instead to criteria that take the value +∞+\infty+∞ off a closed feasible set and are only lower semicontinuous in xxx.

J. Dupačová and R. Wets (IIASA WP-86-41, 1986; journal version Ann. Statist. 16 (1988)) proved consistency in this generality by combining epi-convergence of functions with the theory of measurable multifunctions and normal integrands. This mission formalizes their §3.

Setting

Ξ\XiΞ is a Polish space with its Borel σ\sigmaσ-field and PPP is a probability measure on it. The integrand is f:Rn×Ξ→(−∞,∞]f:\mathbb R^n\times\Xi\to(-\infty,\infty]f:Rn×Ξ→(−∞,∞], and the true problem is to minimize

Ef(x)=∫Ξf(x,ξ) P(dξ),Ef(x)=\int_\Xi f(x,\xi)\,P(d\xi),Ef(x)=∫Ξ​f(x,ξ)P(dξ),

with the convention that Ef(x)=+∞Ef(x)=+\inftyEf(x)=+∞ whenever ξ↦f(x,ξ)\xi\mapsto f(x,\xi)ξ↦f(x,ξ) is not bounded above by a summable function. The effective domain of a function h:Rn→[−∞,∞]h:\mathbb R^n\to[-\infty,\infty]h:Rn→[−∞,∞] is dom⁡h={x:h(x)<∞}\operatorname{dom}h=\{x: h(x)<\infty\}domh={x:h(x)<∞}, and argmin⁡h={x:h(x)=inf⁡h}\operatorname{argmin}h=\{x: h(x)=\inf h\}argminh={x:h(x)=infh}.

Information arrives on a probability space (Z,F,μ)(Z,\mathcal F,\mu)(Z,F,μ) with an increasing sequence of σ\sigmaσ-fields F1⊆F2⊆⋯⊆F\mathcal F^1\subseteq\mathcal F^2\subseteq\dots\subseteq\mathcal FF1⊆F2⊆⋯⊆F. Each sample ζ∈Z\zeta\in Zζ∈Z yields probability measures Pν(⋅,ζ)P^\nu(\cdot,\zeta)Pν(⋅,ζ) on Ξ\XiΞ, and ζ↦Pν(A,ζ)\zeta\mapsto P^\nu(A,\zeta)ζ↦Pν(A,ζ) is Fν\mathcal F^\nuFν-measurable for every Borel AAA: the estimate at stage ν\nuν uses only stage-ν\nuν information. The estimated problem minimizes

Eνf(x,ζ)=∫Ξf(x,ξ) Pν(dξ,ζ).E^\nu f(x,\zeta)=\int_\Xi f(x,\xi)\,P^\nu(d\xi,\zeta).Eνf(x,ζ)=∫Ξ​f(x,ξ)Pν(dξ,ζ).

A sequence gνg^\nugν epi-converges to ggg if, at every xxx, lim inf⁡gν(xν)≥g(x)\liminf g^\nu(x^\nu)\ge g(x)liminfgν(xν)≥g(x) along every sequence xν→xx^\nu\to xxν→x, and lim sup⁡gν(xν)≤g(x)\limsup g^\nu(x^\nu)\le g(x)limsupgν(xν)≤g(x) along some sequence xν→xx^\nu\to xxν→x.

The standing hypotheses are Assumption 3.4: dom⁡f=S×Ξ\operatorname{dom}f=S\times\Xidomf=S×Ξ with SSS closed and nonempty; f(x,⋅)f(x,\cdot)f(x,⋅) is continuous for x∈Sx\in Sx∈S; f(⋅,ξ)f(\cdot,\xi)f(⋅,ξ) is lower semicontinuous; and fff is locally lower Lipschitz on SSS with a bounded continuous modulus β(ξ)\beta(\xi)β(ξ). Assumption 3.5 asks that, for μ\muμ-almost every ζ\zetaζ, Pν(⋅,ζ)P^\nu(\cdot,\zeta)Pν(⋅,ζ) converge in distribution to PPP, that ∣f(x,⋅)∣|f(x,\cdot)|∣f(x,⋅)∣ be uniformly tight along P=P0,P1,…P=P^0,P^1,\dotsP=P0,P1,… for each x∈Sx\in Sx∈S, and that ∫inf⁡xf(x,ξ) Pν(dξ,ζ)>−∞\int\inf_x f(x,\xi)\,P^\nu(d\xi,\zeta)>-\infty∫infx​f(x,ξ)Pν(dξ,ζ)>−∞ for all ν\nuν.

Formalization targets

Goal: Theorem 3.9, "In particular" (pp. 21–22)

Let D⊆RnD\subseteq\mathbb R^nD⊆Rn be compact, suppose (argmin⁡Eνf)∩D≠∅(\operatorname{argmin}E^\nu f)\cap D\neq\emptyset(argminEνf)∩D=∅ μ\muμ-a.s. for every ν\nuν, and suppose {x∗}=argmin⁡Ef∩D\{x^*\}=\operatorname{argmin}Ef\cap D{x∗}=argminEf∩D. Then there are Fν\mathcal F^\nuFν-measurable selections xνx^\nuxν of argmin⁡Eνf\operatorname{argmin}E^\nu fargminEνf with

xν(ζ)→x∗andinf⁡Eνf(⋅,ζ)→inf⁡Effor μ-almost every ζ.x^\nu(\zeta)\to x^*\quad\text{and}\quad \inf E^\nu f(\cdot,\zeta)\to\inf Ef\qquad\text{for }\mu\text{-almost every }\zeta .xν(ζ)→x∗andinfEνf(⋅,ζ)→infEffor μ-almost every ζ.

The goal does not assume that EfEfEf has a unique global minimizer, and it does not assume convexity.

Milestones

In attack order:

  • Proposition 3.3: epi-convergence gives lim sup⁡(inf⁡gν)≤inf⁡g\limsup(\inf g^\nu)\le\inf glimsup(infgν)≤infg, limits of minimizers are minimizers, and the minimum is attained in the closure of a bounded DDD.
  • Lemma 3.6: almost surely, EfEfEf and every EνfE^\nu fEνf are proper and l.s.c., with domain SSS.
  • Theorem 3.7: almost surely, EνfE^\nu fEνf epi-converges and converges pointwise to EfEfEf.
  • Theorem 3.8: almost surely, the epigraphs of EνfE^\nu fEνf are closed, and they depend Fν\mathcal F^\nuFν-measurably on ζ\zetaζ.
  • Theorem 3.9:
    • (3.14) lim sup⁡(inf⁡Eνf)≤inf⁡Ef\limsup(\inf E^\nu f)\le\inf Eflimsup(infEνf)≤infEf a.s.;
    • (i) cluster points of estimated minimizers minimize EfEfEf;
    • (ii) ζ↦argmin⁡Eνf(⋅,ζ)\zeta\mapsto\operatorname{argmin}E^\nu f(\cdot,\zeta)ζ↦argminEνf(⋅,ζ) is closed-valued and Fν\mathcal F^\nuFν-measurable.
  • Proposition 3.1: the measurable selection theorem.

Significance

The result separates two things: the statistical input, which is only convergence in distribution of PνP^\nuPν plus a tightness condition, and the variational output, which is convergence of optimal values and solutions. It therefore applies to any estimator PνP^\nuPν that converges weakly almost surely: empirical measures, kernel estimates, parametric fits. It also covers constrained and nonsmooth problems: the feasible set enters through f=+∞f=+\inftyf=+∞ off SSS, and only lower semicontinuity in xxx is required. Asymptotic distribution results for constrained estimators, such as the second part of the same paper and the subsequent literature on sample average approximation, start from this consistency.

The theorem is proved on paper. To the best of current knowledge none of it is machine-checked. Mathlib has weak convergence of probability measures, lower semicontinuity and extended-real integrals. It does not have epi-convergence, Effros-measurable multifunctions, normal integrands or the Kuratowski–Ryll-Nardzewski selection theorem. A formal proof produces these as reusable components. It also has to supply the details that the paper's proof of Theorem 3.8 leaves as a sketch.

Difficulty

Pointwise convergence Eνf(x)→Ef(x)E^\nu f(x)\to Ef(x)Eνf(x)→Ef(x) is not enough to move minimizers to the limit, and uniform convergence fails because fff is +∞+\infty+∞ off SSS and need not be bounded. Epi-convergence is the right notion. Proving it needs a liminf inequality along moving points xν→xx^\nu\to xxν→x under moving measures PνP^\nuPν. That combines Fatou's lemma, the lower Lipschitz bound and the tightness condition, and the integrands are extended-real-valued, so care is needed.

The second difficulty is measurability. The exceptional null set lies in F\mathcal FF but not in Fν\mathcal F^\nuFν, so "Fν\mathcal F^\nuFν-measurable" has to be understood on a full-measure set in the trace σ\sigmaσ-field. The paper's argument for Theorem 3.8 appeals to continuity of P↦epi⁡EPfP\mapsto\operatorname{epi}E_PfP↦epiEP​f in the epi-topology, and it remarks itself that Theorem 3.7 gives this only along sequences satisfying Assumption 3.5. A solver will have to rebuild this step, for example through the normal-integrand structure of EνfE^\nu fEνf.

Formalization scope

  • Rn\mathbb R^nRn is EuclideanSpace ℝ (Fin n) with its Euclidean norm.
  • Ξ\XiΞ is a Polish space with its Borel σ\sigmaσ-algebra. This is exactly a closed subset of a Polish space with the relative Borel field.
  • fff is EReal-valued. Every expectation is the mission's expect: +∞+\infty+∞ when ∫f+=∞\int f^+=\infty∫f+=∞, and ∫f+−∫f−\int f^+-\int f^-∫f+−∫f− otherwise, both computed as Lebesgue integrals of [0,∞][0,\infty][0,∞]-valued functions. A Bochner integral, which would assign 000 to non-integrable functions, is never used for EfEfEf or EνfE^\nu fEνf.
  • The sample index is shifted: Lean's Pν k and 𝔽 k are the paper's Pk+1P^{k+1}Pk+1 and Fk+1\mathcal F^{k+1}Fk+1, and P=P0P=P^0P=P0 is a separate argument.
  • Infima, lim inf⁡\liminfliminf and lim sup⁡\limsuplimsup are taken in [−∞,∞][-\infty,\infty][−∞,∞].
  • Measurability on the full-measure set Z0Z_0Z0​ uses the trace σ\sigmaσ-field.
  • The selections in the goal are total, Fν\mathcal F^\nuFν-measurable maps Z→RnZ\to\mathbb R^nZ→Rn that select almost surely. This is equivalent to the paper's maps Z0→RnZ_0\to\mathbb R^nZ0​→Rn.
  • "Random l.s.c. function" in Theorem 3.8 is encoded by the equivalent conditions (3.4i)–(3.4ii): nonempty, closed and measurable epigraphs.
  • Lower Lipschitz (3.10) is written additively.
  • The hypothesis that Ξ\XiΞ is the support of PPP is omitted. It is unused in §3, and omitting it strengthens every statement.
  • Nothing beyond the page is assumed: no convexity, no compact SSS, no bounded fff, no unique minimizer, no i.i.d. sampling, no empirical PνP^\nuPν, no completeness of μ\muμ or Fν\mathcal F^\nuFν.

The hypotheses are not vacuous. A sorry-free check verifies all of them, including those of the goal, for f(x,ξ)=∥x∥2f(x,\xi)=\|x\|^2f(x,ξ)=∥x∥2 with Dirac measures. Defining the expectation through a Bochner integral, or dropping S≠∅S\neq\emptysetS=∅ (which makes every argmin⁡\operatorname{argmin}argmin all of Rn\mathbb R^nRn), would trivialize or change the statements; the definitions above rule both out.

Contributions are welcome on any milestone. Proposition 3.1 (Kuratowski–Ryll-Nardzewski for Rm\mathbb R^mRm-valued multifunctions) and Proposition 3.3 (deterministic epi-convergence facts) are independent of the probabilistic setting and reusable beyond this mission.

Selected references

  • J. Dupačová, R. Wets, Asymptotic Behavior of Statistical Estimators and Optimal Solutions for Stochastic Optimization Problems, IIASA Working Paper WP-86-41, 1986. https://pure.iiasa.ac.at/id/eprint/2818/ — journal version: Ann. Statist. 16(4), 1517–1549, 1988. https://doi.org/10.1214/aos/1176351052
  • A. Wald, Note on the consistency of the maximum likelihood estimate, Ann. Math. Statist. 20, 595–601, 1949. https://doi.org/10.1214/aoms/1177729938
  • P. J. Huber, The behavior of maximum likelihood estimates under nonstandard conditions, Proc. Fifth Berkeley Symp. Math. Statist. Probab. 1, 221–233, 1967. https://projecteuclid.org/euclid.bsmsp/1200512988
  • R. T. Rockafellar, R. J.-B. Wets, Variational Analysis, Springer, 1998 (Ch. 7 epi-convergence; Ch. 14 measurable multifunctions and normal integrands). https://doi.org/10.1007/978-3-642-02431-3
13 thms1 active userReviewed
Linear OptimizationOperations ResearchOptimization+2·Captain: mikedeng1

Smoothed Analysis of Algorithms: Why the Simplex Algorithm Usually Takes Polynomial Time 2: The Two-Phase Shadow-Vertex Simplex Method Has Polynomial Smoothed ComplexityResearch Paper

Motivation

The simplex method solves linear programs by moving between vertices of a feasible polyhedron. Its worst-case number of moves can grow exponentially, yet it often performs well on ordinary inputs. Worst-case examples alone therefore give an incomplete account of the method’s behavior. Spielman and Teng introduced smoothed analysis to measure expected performance after small random perturbations of an arbitrary input. Their result for a two-phase shadow-vertex simplex method gives a polynomial bound in the input dimensions and inverse perturbation scale. The pinned preprint is the source for every theorem number and constant in this mission.

The paper separates a geometric result about the expected size of a polytope’s shadow (Theorem 4.0.1) from the algorithmic result here (Theorem 5.0.1). That separation matters: a plane chosen before perturbation and a plane chosen by a running algorithm have different distributions. This mission addresses the latter. It complements the standard-form simplex theorems already formalized in the Introduction to Linear Optimization series and the worst-case Klee–Minty result in the Smale’s Ninth Problem mission; those results concern different algorithms or input models and are context rather than imported statements.

Setting

A linear program is specified by vectors a1,…,an∈Rda_1,\ldots,a_n\in\mathbb R^da1​,…,an​∈Rd, right-hand sides y1,…,yn∈Ry_1,\ldots,y_n\in\mathbb Ry1​,…,yn​∈R, and an objective vector z∈Rdz\in\mathbb R^dz∈Rd:

max⁡x⟨z,x⟩subject to⟨ai,x⟩≤yi(1≤i≤n).\max_x\langle z,x\rangle\quad\text{subject to}\quad \langle a_i,x\rangle\le y_i\qquad(1\le i\le n).xmax​⟨z,x⟩subject to⟨ai​,x⟩≤yi​(1≤i≤n).

The paper’s two-phase shadow-vertex method first draws a collection I\mathcal II of ddd-element subsets of [n][n][n] and chooses one whose constraint matrix AIA_IAI​ has the largest smallest singular value. It sets a power-of-two scale MMM from the input norm and a power-of-two scale κ\kappaκ from that singular value. These determine positive relaxed right-hand sides yi′y'_iyi′​: MMM for i∈Ii\in Ii∈I and dM2/(4κ)\sqrt d M^2/(4\kappa)d​M2/(4κ) otherwise. A coefficient vector α\alphaα is chosen uniformly from A1/d2={α:∑i∈Iαi=1, αi≥1/d2}A_{1/d^2}=\{\alpha:\sum_{i\in I}\alpha_i=1,\ \alpha_i\ge1/d^2\}A1/d2​={α:∑i∈I​αi​=1, αi​≥1/d2}. The first phase solves the relaxed program LP′ from the objective AIαA_I\alphaAI​α.

The second phase uses a lifted program LP⁺ in Rd+1\mathbb R^{d+1}Rd+1. For each original constraint it forms ai+=((yi′−yi)/2,ai)a_i^+=((y'_i-y_i)/2,a_i)ai+​=((yi′​−yi​)/2,ai​) and yi+=(yi′+yi)/2y_i^+=(y'_i+y_i)/2yi+​=(yi′​+yi​)/2, together with two artificial constraints at first coordinates 111 and −1-1−1. LP⁺ connects LP′ to the original program and makes infeasibility detectable. Its shadow is taken in the plane of (0,z)(0,z)(0,z) and z+=(1,0,…,0)z^+=(1,0,\ldots,0)z+=(1,0,…,0).

For positive right-hand sides, an optimal polar simplex is a ddd-subset of constraints whose scaled vectors ai/yia_i/y_iai​/yi​ form a facet of ConvHull⁡(0,a1/y1,…,an/yn)\operatorname{ConvHull}(0,a_1/y_1,\ldots,a_n/y_n)ConvHull(0,a1​/y1​,…,an​/yn​) and whose unscaled cone contains an objective qqq. The shadow for objectives t,zt,zt,z is the union of these simplices over all qqq in Span⁡(t,z)\operatorname{Span}(t,z)Span(t,z). Its size bounds the number of polar pivots. In Section 5 the paper writes Sz′S'_zSz′​ for the first-phase shadow size and Sz+S_z^+Sz+​ for the second-phase shadow size without the two artificial pivots.

The input is perturbed by independent Gaussians: each coordinate of aia_iai​ and each yiy_iyi​ has its prescribed center and common standard deviation σR\sigma RσR, where R=max⁡i∥(yˉi,aˉi)∥2R=\max_i\|(\bar y_i,\bar a_i)\|_2R=maxi​∥(yˉ​i​,aˉi​)∥2​. The algorithm has separate random choices of I\mathcal II and α\alphaα.

Formalization targets

The immediate targets bound the two phases: Lemma 5.2.1 gives an explicit expectation bound for Sz′S'_zSz′​ and Lemma 5.3.1 gives one for Sz+S_z^+Sz+​. Lemma 5.1.1 and its corollaries control the chance that the chosen basis has a very small singular value. Corollary 4.3.3 extends the geometric shadow bound to positive, unequal right-hand sides and general Gaussian covariance. These are the mission’s milestone targets.

The goal is the shape of Theorem 5.0.1. With C(A,y,z)=EI,α(Sz′+Sz++2)C(A,y,z)=\mathbb E_{\mathcal I,\alpha}(S'_z+S_z^++2)C(A,y,z)=EI,α​(Sz′​+Sz+​+2), there are a single polynomial P\mathcal PP and a positive constant σ0\sigma_0σ0​ such that, for all n>d≥3n>d\ge3n>d≥3 and all centers and objectives,

EA,yC(A,y,z)≤min⁡{P(d,n,1min⁡(σ,σ0)),(nd)+(nd+1)+2}.\mathbb E_{A,y}C(A,y,z)\le \min\left\{\mathcal P\left(d,n,\frac1{\min(\sigma,\sigma_0)}\right), \binom nd+\binom n{d+1}+2\right\}.EA,y​C(A,y,z)≤min{P(d,n,min(σ,σ0​)1​),(dn​)+(d+1n​)+2}.

The polynomial is uniform over the dimensions and inputs; its coefficients are not prescribed. The bound on CCC implies the corresponding result for the actual pivot count through the paper’s step-to-shadow comparison. The goal is stated with a positive center scale RRR, the case in which the paper’s Gaussian rescaling applies.

Significance

The theorem places the number of pivots of a complete simplex method under one explicit perturbation model, including the work needed to find a starting feasible basis and handle an arbitrary right-hand side. The trivial binomial bound is retained because it controls rare events in the proof and is part of the stated result. The polynomial bound says that even when the unperturbed LP is adversarial, Gaussian noise of a controlled scale makes the expected shadow-size cost polynomial.

The paper proves the mathematical result. This mission asks for machine-checked proofs of its statement and the listed milestones; the draft Lean declarations are targets with sorry, not completed proofs. The reusable formal infrastructure is the finite polar simplex and shadow construction, product Gaussian input law, smallest-singular-value events for sampled minors, and the uniform truncated-simplex coefficient law. The two shadow-size lemmas also require explicit handling of measurable finite-valued counts and their expectations.

Difficulty

The basic shadow estimate fixes its projection plane before perturbing the constraints. In LP′, the initial objective AIαA_I\alphaAI​α uses a basis selected after the perturbation, so the relevant plane depends on the random LP. The fixed-plane theorem cannot be substituted directly. For LP⁺, the normalized lifted vectors ai+/yi+a_i^+/y_i^+ai+​/yi+​ are nonlinear functions of Gaussian data; they are generally not Gaussian vectors. Thus the same shadow estimate does not apply directly to their law either. A further issue is that a poor sampled basis can make y′y'y′ very large. These are distinct obstacles, reflected in the milestone groups from Sections 5.1, 5.2, and 5.3.

Formalization scope

Vectors are EuclideanSpace ℝ (Fin d), constraints are Fin n → EuclideanSpace ℝ (Fin d), and index families are finite sets of Fin n. The paper’s [n][n][n] starts at one; Fin n starts at zero. The Gaussian constructor receives variance σ2\sigma^2σ2, not standard deviation σ\sigmaσ. The 3ndln⁡n3nd\ln n3ndlnn draws are rounded upward and are independent uniform draws with replacement. Equal singular values are resolved by the first sampled set. The uniform law on AδA_\deltaAδ​ is represented by normalized independent exponential weights followed by the affine shift that imposes αi≥δ\alpha_i\ge\deltaαi​≥δ.

The Lean definition of CCC is exactly the Section 5 shadow-size upper bound E(Sz′+Sz++2)\mathbb E(S'_z+S_z^++2)E(Sz′​+Sz+​+2), computed from the sampled LP data. It is not an arbitrary cost variable. The actual algorithmic step bound needs the paper’s polar algorithm and Lemma 3.3.5. The goal explicitly asks for inner and outer integrability so Lean’s default value for a nonintegrable Bochner integral cannot make the result vacuous. The source’s all-zero center scale is excluded because it gives zero perturbation and defeats the rescaling used in Theorem 5.0.1.

For LP⁺ the vectors live in Rd+1\mathbb R^{d+1}Rd+1, so the two LP⁺ milestone bounds use D(n,d+1,⋅)\mathcal D(n,d+1,\cdot)D(n,d+1,⋅). The preprint prints ddd in those calls even though the preceding extension theorem would be applied in dimension d+1d+1d+1. Lemma 5.2.1 is written as an inequality: its printed equality is stronger than the bound established on page 71. These corrections are visible in the theorem titles and notes. Contributions that prove the exact statements, establish the measurability and Gaussian law facts, or formalize the step-to-shadow comparison are welcome.

Selected references

  • Daniel A. Spielman and Shang-Hua Teng, Smoothed Analysis of Algorithms: Why the Simplex Algorithm Usually Takes Polynomial Time, arXiv:cs/0111050v7, 2003, preprint. The PDF used here is the 96-page version with printed and PDF page numbers aligned.
22 thms2 active usersReviewed
Dynamic ProgrammingMarkov ChainOperations Research+1·Captain: mikedeng1

An Inventory Model with Limited Production Capacity and Uncertain Demands I. The Average-Cost Criterion: With Finite Storage a Modified Base-Stock Policy Is Strongly Average-Cost OptimalResearch Paper

Motivation

A manufacturer that makes one product to stock faces random demand, can produce at most bbb units per period, and can store at most UUU units. The classical result without the production limit is that a base-stock policy is optimal: raise inventory to a fixed level yˉ\bar yyˉ​ each period. With a production limit, the natural modification is to produce up to yˉ\bar yyˉ​ when that is possible and to produce at full capacity otherwise. Federgruen and Zipkin (1986) proved that this modified base-stock (critical-number) policy is optimal under the long-run average-cost criterion, for discrete demand with a general convex cost. Production-capacity models of this type are standard in operations management texts, and the result underlies the computational and comparative-static work that followed, starting with Part II of the same paper, which treats discounted costs.

Timeline.

  • 1950s–60s: optimality of base-stock (critical-number) policies for uncapacitated periodic-review models; see Heyman and Sobel's Stochastic Models in Operations Research, Vol. II (1984).
  • 1986: Federgruen and Zipkin, Part I (average cost, MOR 11(2):193–207) and Part II (discounted cost, MOR 11(2):208–215) establish the capacitated case. Part I handles the unbounded state space with a general average-cost theory for countable-state Markov decision processes by Federgruen, Schweitzer and Tijms (1983).

Setting

Time is divided into periods t=0,1,…t = 0, 1, \dotst=0,1,…. The demands D0,D1,…D_0, D_1, \dotsD0​,D1​,… are independent copies of a random variable DDD with values in {0,1,2,… }\{0, 1, 2, \dots\}{0,1,2,…} and probability mass function p(j)p(j)p(j); write μ=E(D)\mu = E(D)μ=E(D) and P(j)=Pr⁡{D≤j}P(j) = \Pr\{D \le j\}P(j)=Pr{D≤j}. At the start of period ttt the inventory is an integer xtx_txt​ (negative values are backorders). The decision maker raises it to

yt∈Y(xt)={y∈Z:xt≤y≤xt+b, y≤U},y_t \in Y(x_t) = \{y \in \mathbb Z : x_t \le y \le x_t + b,\ y \le U\},yt​∈Y(xt​)={y∈Z:xt​≤y≤xt​+b, y≤U},

pays the expected one-period cost G(yt)G(y_t)G(yt​), and demand is subtracted: xt+1=yt−Dtx_{t+1} = y_t - D_txt+1​=yt​−Dt​. The order cost per unit is set to zero, as in the paper; this loses no generality because every policy with finite average cost has the same average order cost.

The standing assumptions are: G≥0G \ge 0G≥0 is convex and G(y)→∞G(y) \to \inftyG(y)→∞ as ∣y∣→∞|y| \to \infty∣y∣→∞ (Assumption 1); the characteristic function of DDD is analytic at the origin (Assumption 2), and 0<μ0 < \mu0<μ; G(y)≤A+B∣y∣ρG(y) \le A + B|y|^\rhoG(y)≤A+B∣y∣ρ for some positive integer ρ\rhoρ (Assumption 3); b>μb > \mub>μ and P(b)<1P(b) < 1P(b)<1 (Assumption 4). The smallest global minimizer of GGG is yˉ∞\bar y^\inftyyˉ​∞, and U≥yˉ∞U \ge \bar y^\inftyU≥yˉ​∞.

A Markov policy is a sequence π=(π0,π1,… )\pi = (\pi_0, \pi_1, \dots)π=(π0​,π1​,…) of maps with πt(x)∈Y(x)\pi_t(x) \in Y(x)πt​(x)∈Y(x). The critical-number policy with critical number yˉ\bar yyˉ​ is δ[yˉ](x)=max⁡(x,min⁡(yˉ,x+b))\delta[\bar y](x) = \max(x, \min(\bar y, x + b))δ[yˉ​](x)=max(x,min(yˉ​,x+b)). A stationary policy δ\deltaδ is strongly optimal with average cost ggg if, from every initial state x≤Ux \le Ux≤U, its average cost t−1E{∑i<tG(yi)}t^{-1}E\{\sum_{i<t} G(y_i)\}t−1E{∑i<t​G(yi​)} converges to ggg, while every Markov policy has lim-inf average cost at least ggg from every initial state.

The analysis uses the operators Rv(y)=G(y)+E v(y−D)Rv(y) = G(y) + E\,v(y - D)Rv(y)=G(y)+Ev(y−D) and Sv(x)=min⁡y∈Y(x)Rv(y)Sv(x) = \min_{y \in Y(x)} Rv(y)Sv(x)=miny∈Y(x)​Rv(y), and the optimality equation

g+v(x)=Sv(x),x≤U.(6)g + v(x) = Sv(x),\qquad x \le U. \tag{6}g+v(x)=Sv(x),x≤U.(6)

For an interval ι=[l,u]\iota = [l, u]ι=[l,u], Hιv(x)H_\iota v(x)Hι​v(x) is the largest expected sum of v(yt)v(y_t)v(yt​), over policies forced to produce at capacity below lll and to produce nothing above uuu, until the inventory first returns to ι\iotaι.

Formalization targets

Goal: Theorem 1 (p. 202)

There exist g∗g^*g∗, v∗v^*v∗ and y∗≥yˉ∞y^* \ge \bar y^\inftyy∗≥yˉ​∞ such that (g∗,v∗)(g^*, v^*)(g∗,v∗) solves (6), v∗v^*v∗ is convex with global minimizer y∗y^*y∗, and

δ∗=δ[y∗] is strongly optimal with average cost g∗.\delta^* = \delta[y^*] \text{ is strongly optimal with average cost } g^*.δ∗=δ[y∗] is strongly optimal with average cost g∗.

The y∗y^*y∗ in the optimality claim is the minimizer constructed in part (a).

Milestones

  • Lemma 2(a)–(c) (pp. 196–197): a normal-tail inequality and two series estimates.
  • Lemma 3 (p. 198): if v(x)=O(∣x∣q)v(x) = O(|x|^q)v(x)=O(∣x∣q) then Hιv(x)=O(∣x∣q+3)H_\iota v(x) = O(|x|^{q+3})Hι​v(x)=O(∣x∣q+3).
  • Corollary 1 (p. 200): Hι1=O(∣x∣3)H_\iota 1 = O(|x|^3)Hι​1=O(∣x∣3) and HιG=O(∣x∣ρ+3)H_\iota G = O(|x|^{\rho+3})Hι​G=O(∣x∣ρ+3), both finite.
  • Corollary 2 (p. 201): (t+1)−1P[δ0t]⋯P[δtt](Hι1+HιG)(x)→0(t+1)^{-1}P[\delta_{0t}]\cdots P[\delta_{tt}](H_\iota 1 + H_\iota G)(x) \to 0(t+1)−1P[δ0t​]⋯P[δtt​](Hι​1+Hι​G)(x)→0.
  • Lemma 4 (p. 201): reachability of every state in [L,U−D−][L, U - D_-][L,U−D−​] under some policy that produces at capacity below LLL.
  • Lemma 5 (p. 202): SSS and QQQ preserve the class VVV of convex functions of growth O(∣x∣ρ+3)O(|x|^{\rho+3})O(∣x∣ρ+3) that are nonincreasing below yˉ∞\bar y^\inftyyˉ​∞.

Significance

The result. Theorem 1 reduces an infinite-state average-cost control problem to a one-parameter search over critical numbers. The paper then evaluates the average cost of δ[yˉ]\delta[\bar y]δ[yˉ​] by a renewal formula, proves it convex in yˉ\bar yyˉ​ (Theorem 2), and in §5 extends optimality to unlimited storage. The strong form of optimality matters: it compares with every Markov policy from every starting state, and it compares lim-infs, not only lim-sups.

Formalizing it. The theorem has a published proof, but no machine-checked one, and its proof relies on external results that are themselves unformalized: the countable-state average-cost theory of Federgruen, Schweitzer and Tijms, a fixed-point theorem on a compact convex subset of a product space, and a large-deviation estimate quoted from Feller. A formal development produces reusable infrastructure: expected first-passage sums for integer-valued random walks with a reflecting control, polynomial moment bounds for them, and the convexity-preservation argument for capacitated value iteration.

Difficulty

The state space is unbounded below, so the finite-state theory of average-cost Markov decision processes does not apply, and the one-period cost is unbounded. The obvious approach, letting the discount factor tend to one in the discounted problem, needs uniform bounds on relative value functions. Those bounds come from the expected cost until the inventory returns to a fixed interval, and with capacity limits that expectation must be controlled with growth O(∣x∣ρ+3)O(|x|^{\rho+3})O(∣x∣ρ+3) uniformly over a class of policies. This is the content of Lemma 3, whose proof combines a large-deviation estimate for the demand sums with a renewal-type recursion. A second obstacle is strong optimality: comparing with policies whose lim-inf average cost is smaller requires that the relative value function grows sublinearly along every admissible trajectory (Corollary 2).

Formalization scope

All objects are in the namespace FedergruenZipkin.AvgCost, defined in one file. States x,yx, yx,y and the capacity UUU are integers; demands are natural numbers with a real probability mass function p; bbb is a positive natural number. Convexity on Z\mathbb ZZ is the second-difference inequality. Expectations of a real function are series ∑jp(j) v(y−j)\sum_j p(j)\,v(y-j)∑j​p(j)v(y−j); expected policy costs and hitting sums are [0,∞][0,\infty][0,∞]-valued and need no integrability side condition. Feasibility and all properties of value functions are required only on states x≤Ux \le Ux≤U, which are the only states visited. Assumption 2 is stated literally, as real-analyticity of θ↦∑jp(j)eiθj\theta \mapsto \sum_j p(j)e^{i\theta j}θ↦∑j​p(j)eiθj at 000. The order cost is zero, as in the paper. yˉ∞\bar y^\inftyyˉ​∞ is a parameter characterised as the least minimizer of GGG, not an infimum.

"Strongly optimal" has no displayed definition in the paper; it is read from eq. (7) in the proof of Theorem 1(b): convergence of the average cost of δ∗\delta^*δ∗ to g∗g^*g∗ from every state, together with a lim-inf lower bound for every Markov (memoryless, possibly nonstationary) policy from every state. The class is neither widened to history-dependent policies nor narrowed to stationary ones. The goal additionally records that E v∗(y−D)E\,v^*(y-D)Ev∗(y−D) converges, that v∗v^*v∗ has growth O(∣x∣ρ+3)O(|x|^{\rho+3})O(∣x∣ρ+3), and that g∗≥0g^* \ge 0g∗≥0; all three follow from the paper's proof.

A trivializing reading is ruled out: the existence of ggg, vvv and y∗y^*y∗ is one existential, so y∗y^*y∗ cannot be decoupled from the solution of (6), and strong optimality includes the convergence of δ∗\delta^*δ∗'s own average cost to g∗g^*g∗, so g=0g = 0g=0 does not satisfy it vacuously.

Not posed: Lemma 1 (quoted from Feller, and replaceable by a Chernoff bound); the renewal formulas (10)–(11) and Theorem 2; and §5 (unlimited storage). Useful contributions include a formal theory of expected hitting sums for skip-free-upward random walks, and a proof of Lemma 3 by any route.

Selected references

  • A. Federgruen and P. Zipkin, An Inventory Model with Limited Production Capacity and Uncertain Demands I. The Average-Cost Criterion, Mathematics of Operations Research 11(2):193–207, 1986. https://doi.org/10.1287/moor.11.2.193
  • A. Federgruen and P. Zipkin, An Inventory Model with Limited Production Capacity and Uncertain Demands II. The Discounted-Cost Criterion, Mathematics of Operations Research 11(2):208–215, 1986. https://doi.org/10.1287/moor.11.2.208
  • A. Federgruen, P. J. Schweitzer and H. C. Tijms, Denumerable Undiscounted Semi-Markov Decision Processes with Unbounded Rewards, Mathematics of Operations Research 8(2):298–314, 1983. https://doi.org/10.1287/moor.8.2.298
  • D. P. Heyman and M. J. Sobel, Stochastic Models in Operations Research, Vol. II, McGraw-Hill, 1984.
  • W. Feller, An Introduction to Probability Theory and Its Applications, Vol. II, 2nd ed., Wiley, 1971.
10 thms3 active usersReviewed
Dynamical SystemsOperations ResearchProbability+1·Captain: mikedeng1

Dynamics of Stochastic Approximation Algorithms 3: Martingale Noise with Bounded q-th Moments and Summable γ_n^(1+q/2) Satisfies Assumption A1 Almost SurelyResearch Paper

Motivation

A stochastic approximation algorithm is a recursion

xn+1−xn=γn+1(F(xn)+Un+1)x_{n+1}-x_n=\gamma_{n+1}\big(F(x_n)+U_{n+1}\big)xn+1​−xn​=γn+1​(F(xn​)+Un+1​)

in Rd\mathbb R^dRd, where FFF is a vector field, γn\gamma_nγn​ are small step sizes and Un+1U_{n+1}Un+1​ is noise. Such recursions go back to Robbins and Monro's root-finding scheme (Robbins–Monro 1951) and underlie stochastic gradient descent, temporal-difference learning, adaptive control and learning in games. The ODE method studies them by comparing the iterates with the trajectories of x˙=F(x)\dot x=F(x)x˙=F(x).

Benaïm's lecture notes (Benaïm 1999) organize the ODE method in two steps. A deterministic step, Proposition 4.1, shows that whenever the noise satisfies a condition called A1 (and the iterates are bounded, or FFF is Lipschitz and bounded on a neighbourhood of them), the interpolated process is an asymptotic pseudotrajectory of the flow of FFF. A probabilistic step then verifies A1 for concrete noise models. This mission formalizes the first such verification, Proposition 4.2: martingale difference noise with bounded qqq-th moments and step sizes with ∑nγn1+q/2<∞\sum_n\gamma_n^{1+q/2}<\infty∑n​γn1+q/2​<∞. The result is described as a particular case of a general theorem of Métivier and Priouret (1987); the same estimates reappear later in the notes.

Setting

Let {γn}n≥1\{\gamma_n\}_{n\ge1}{γn​}n≥1​ be a deterministic sequence with γn≥0\gamma_n\ge0γn​≥0, ∑nγn=∞\sum_n\gamma_n=\infty∑n​γn​=∞ and γn→0\gamma_n\to0γn​→0 (a step sequence). Put τ0=0\tau_0=0τ0​=0, τn=∑i=1nγi\tau_n=\sum_{i=1}^n\gamma_iτn​=∑i=1n​γi​, and let

m(t)=sup⁡{k≥0: t≥τk}m(t)=\sup\{k\ge0:\ t\ge\tau_k\}m(t)=sup{k≥0: t≥τk​}

be the index of the step that contains time t≥0t\ge0t≥0. For a sequence {Un}n≥1\{U_n\}_{n\ge1}{Un​}n≥1​ define the piecewise constant processes Uˉ(t)=Um(t)+1\bar U(t)=U_{m(t)+1}Uˉ(t)=Um(t)+1​ and γˉ(t)=γm(t)+1\bar\gamma(t)=\gamma_{m(t)+1}γˉ​(t)=γm(t)+1​, so that step n+1n+1n+1 occupies the time interval [τn,τn+1)[\tau_n,\tau_{n+1})[τn​,τn+1​) of length γn+1\gamma_{n+1}γn+1​.

Assumption A1 asks that for every T>0T>0T>0

lim⁡n→∞sup⁡{∥∑i=nk−1γi+1Ui+1∥: k=n+1,…,m(τn+T)}=0,\lim_{n\to\infty}\sup\Big\{\Big\|\sum_{i=n}^{k-1}\gamma_{i+1}U_{i+1}\Big\|:\ k=n+1,\dots,m(\tau_n+T)\Big\}=0,n→∞lim​sup{​i=n∑k−1​γi+1​Ui+1​​: k=n+1,…,m(τn​+T)}=0,

or, in the form the notes call equivalent, lim⁡t→∞Δ(t,T)=0\lim_{t\to\infty}\Delta(t,T)=0limt→∞​Δ(t,T)=0 for every T>0T>0T>0, where

Δ(t,T)=sup⁡0≤h≤T∥∫tt+hUˉ(s) ds∥.\Delta(t,T)=\sup_{0\le h\le T}\Big\|\int_t^{t+h}\bar U(s)\,ds\Big\|.Δ(t,T)=0≤h≤Tsup​​∫tt+h​Uˉ(s)ds​.

Let (Ω,F,P)(\Omega,\mathcal F,P)(Ω,F,P) be a probability space with a nondecreasing sequence {Fn}\{\mathcal F_n\}{Fn​} of sub-σ\sigmaσ-algebras, and F:Rd→RdF:\mathbb R^d\to\mathbb R^dF:Rd→Rd continuous. A sequence {xn}\{x_n\}{xn​} given by the recursion above is a Robbins–Monro algorithm if γ\gammaγ is deterministic, UnU_nUn​ is Fn\mathcal F_nFn​-measurable, and E(Un+1∣Fn)=0E(U_{n+1}\mid\mathcal F_n)=0E(Un+1​∣Fn​)=0.

Formalization targets

Goal: Proposition 4.2

For a Robbins–Monro algorithm and some real q≥2q\ge2q≥2, if

sup⁡nE(∥Un+1∥q)<∞and∑nγn1+q/2<∞,\sup_nE\big(\|U_{n+1}\|^q\big)<\infty\qquad\text{and}\qquad\sum_n\gamma_n^{1+q/2}<\infty,nsup​E(∥Un+1​∥q)<∞andn∑​γn1+q/2​<∞,

then with probability one the realised noise sequence satisfies A1, in both of its forms, simultaneously for all T>0T>0T>0.

Milestones

  1. Eq. (13), an instance of Burkholder's inequality with a universal constant CqC_qCq​:
E{sup⁡n<k≤m(τn+T)∥∑i=nk−1γi+1Ui+1∥q}≤Cq E{[∑i=nm(τn+T)−1γi+12∥Ui+1∥2]q/2}.E\Big\{\sup_{n<k\le m(\tau_n+T)}\Big\|\sum_{i=n}^{k-1}\gamma_{i+1}U_{i+1}\Big\|^q\Big\}\le C_q\,E\Big\{\Big[\sum_{i=n}^{m(\tau_n+T)-1}\gamma_{i+1}^2\|U_{i+1}\|^2\Big]^{q/2}\Big\}.E{n<k≤m(τn​+T)sup​​i=n∑k−1​γi+1​Ui+1​​q}≤Cq​E{[i=n∑m(τn​+T)−1​γi+12​∥Ui+1​∥2]q/2}.
  1. Inequality (14), for finite families with αi≥0\alpha_i\ge0αi​≥0, u>1u>1u>1, 0<δ<10<\delta<10<δ<1:
(∑i∣αiβi∣)u≤(∑iαiδu/(u−1))u−1∑iαi(1−δ)u∣βi∣u.\Big(\sum_i|\alpha_i\beta_i|\Big)^u\le\Big(\sum_i\alpha_i^{\delta u/(u-1)}\Big)^{u-1}\sum_i\alpha_i^{(1-\delta)u}|\beta_i|^u.(i∑​∣αi​βi​∣)u≤(i∑​αiδu/(u−1)​)u−1i∑​αi(1−δ)u​∣βi​∣u.
  1. Eq. (16): for every T>0T>0T>0 there is C(q,T)C(q,T)C(q,T) with E(Δ(t,T)q)≤C(q,T)∫tt+Tγˉq/2(s) dsE(\Delta(t,T)^q)\le C(q,T)\int_t^{t+T}\bar\gamma^{q/2}(s)\,dsE(Δ(t,T)q)≤C(q,T)∫tt+T​γˉ​q/2(s)ds for all t≥0t\ge0t≥0.
  2. Eq. (17): ∑k≥0E(Δ(kT,T)q)<∞\sum_{k\ge0}E(\Delta(kT,T)^q)<\infty∑k≥0​E(Δ(kT,T)q)<∞ for every T>0T>0T>0.
  3. Block comparison: Δ(t,T)≤2Δ(kT,T)+Δ((k+1)T,T)\Delta(t,T)\le2\Delta(kT,T)+\Delta((k+1)T,T)Δ(t,T)≤2Δ(kT,T)+Δ((k+1)T,T) for kT≤t<(k+1)TkT\le t<(k+1)TkT≤t<(k+1)T.

Significance

Proposition 4.2 is the standard sufficient condition under which the ODE method applies to stochastic gradient-type recursions with martingale noise. With q=2q=2q=2 it covers step sizes with ∑γn2<∞\sum\gamma_n^2<\infty∑γn2​<∞ (for example γn=1/n\gamma_n=1/nγn​=1/n) and noise with bounded variance; larger qqq trades stronger moment assumptions for slower decay of the steps, down to ∑γn1+q/2<∞\sum\gamma_n^{1+q/2}<\infty∑γn1+q/2​<∞. Combined with Proposition 4.1 it shows that the interpolated process of a Robbins–Monro algorithm with bounded iterates is almost surely an asymptotic pseudotrajectory of the flow of FFF; the limit set theorems of the notes then locate the limit points of the algorithm.

The result is proved in the notes and in the cited literature; it has not, to our knowledge, been machine-checked. A formal proof would supply reusable pieces that Mathlib currently lacks, most notably a Burkholder (or Burkholder–Davis–Gundy) inequality for discrete-time vector martingales in LqL^qLq, and the continuous-time bookkeeping of the step processes Uˉ\bar UUˉ, γˉ\bar\gammaγˉ​ and the noise deviation Δ\DeltaΔ, shared by the other missions of this series.

Difficulty

The obvious argument controls each window by Doob's L2L^2L2 maximal inequality and sums over windows. That works for q=2q=2q=2 only. For q>2q>2q>2 the second moment of the window sums is not summable under ∑γn1+q/2<∞\sum\gamma_n^{1+q/2}<\infty∑γn1+q/2​<∞, and one needs an LqL^qLq maximal inequality whose right-hand side is the q/2q/2q/2-th moment of the square function. That inequality, Burkholder's, is not in Mathlib. Converting the square function into the moment bound requires a Hölder-type inequality with tuned exponents, and passing from the discrete sums to Δ(t,T)\Delta(t,T)Δ(t,T) requires handling partial steps at both ends of [t,t+h][t,t+h][t,t+h]. A second subtlety is that A1 quantifies over all T>0T>0T>0: the almost-sure statement must hold on a single event of full probability for every TTT, not on an event that depends on TTT.

Formalization scope

The space is Rd\mathbb R^dRd as EuclideanSpace ℝ (Fin d); time is real; qqq is a real number with q≥2q\ge2q≥2, and all powers are real powers of nonnegative quantities. The sequences γ\gammaγ and UUU are indexed by N\mathbb NN, and their values at 000 are unused, as the paper indexes them from 111. The filtration is a Mathlib Filtration ℕ, Un+1U_{n+1}Un+1​ is required to be Fn+1\mathcal F_{n+1}Fn+1​-strongly measurable and integrable, and the martingale difference condition is E(Un+1∣Fn)=0E(U_{n+1}\mid\mathcal F_n)=0E(Un+1​∣Fn​)=0 almost surely. Expectations of nonnegative quantities, the suprema in A1 and Δ\DeltaΔ, and the moment bound are taken in [0,∞][0,\infty][0,∞], so no default value of a non-integrable expectation or of an empty supremum can make a statement hold vacuously; the supremum over an empty range of kkk is 000.

The following readings are excluded and are not acceptable formalizations: a moment hypothesis that holds vacuously, a conditional expectation hypothesis on non-integrable noise, the conclusion "for each TTT, A1 holds almost surely" in place of "almost surely, A1 holds for all TTT", and qqq fixed to 222 or restricted to integers.

All hypotheses are satisfiable: U=0U=0U=0, x=0x=0x=0, F=0F=0F=0 and γn=1/n\gamma_n=1/nγn​=1/n with q=2q=2q=2 satisfy every one of them.

Contributions welcome: a general Burkholder inequality for discrete-time martingales in finite-dimensional spaces (reusable well beyond this mission), lemmas on the step processes and Δ\DeltaΔ (measurability, local integrability, additivity), and the proofs of the milestones.

Selected references

  • M. Benaïm, Dynamics of Stochastic Approximation Algorithms, Séminaire de Probabilités XXXIII, Lecture Notes in Mathematics 1709, Springer, 1999, pp. 1–68. https://doi.org/10.1007/BFb0096509
  • M. Benaïm and M. W. Hirsch, Asymptotic pseudotrajectories and chain recurrent flows, with applications, Journal of Dynamics and Differential Equations 8 (1996), 141–176. https://doi.org/10.1007/BF02218617
  • M. Métivier and P. Priouret, Théorèmes de convergence presque sûre pour une classe d'algorithmes stochastiques à pas décroissant, Probability Theory and Related Fields 74 (1987), 403–428.
  • D. L. Burkholder, Distribution function inequalities for martingales, Annals of Probability 1 (1973), 19–42. https://doi.org/10.1214/aop/1176997023
  • D. W. Stroock, Probability Theory: An Analytic View, Cambridge University Press, 1993.
  • H. Robbins and S. Monro, A stochastic approximation method, Annals of Mathematical Statistics 22 (1951), 400–407. https://doi.org/10.1214/aoms/1177729586
  • H. J. Kushner and G. G. Yin, Stochastic Approximation Algorithms and Applications, Springer, 1997.
10 thms2 active usersReviewed
Dynamical SystemsOperations ResearchProbability+1·Captain: mikedeng1

Dynamics of Stochastic Approximation Algorithms 6: An Attractor Whose Basin Meets the Attainable Set Contains the Limit Set with Positive ProbabilityResearch Paper

Motivation

A stochastic approximation algorithm is a recursion xn+1=xn+γn+1(F(xn)+Un+1)x_{n+1}=x_n+\gamma_{n+1}\big(F(x_n)+U_{n+1}\big)xn+1​=xn​+γn+1​(F(xn​)+Un+1​) with decreasing step sizes γn\gamma_nγn​ and a noise term Un+1U_{n+1}Un+1​. Recursions of this form appear in stochastic gradient methods, adaptive control, learning in games (fictitious play, reinforcement learning) and urn models. The ODE method compares the iterates with the solutions of x˙=F(x)\dot x=F(x)x˙=F(x). In Benaïm's lecture notes (Benaïm 1999), the comparison is phrased through the continuous-time interpolated process XXX. Under standard noise conditions, XXX is almost surely an asymptotic pseudotrajectory of the flow of FFF, and its limit set is almost surely internally chain transitive.

That theorem constrains where the process may end up. It does not say which of several candidate sets the process actually reaches. When the ODE has several attractors, for example several stable equilibria of a learning dynamic or several stable compositions of an urn, an application needs to know that each attractor is reached with positive probability. Section 7 of the notes answers this question. The answer is a criterion of attainability: if the process can, with positive probability and at arbitrarily late times, enter the basin of an attractor, then it converges to that attractor with positive probability.

Timeline.

  • Kushner and Clark (1978) proved convergence statements for processes that visit a compact subset of the domain of attraction of an asymptotically stable equilibrium infinitely often.
  • Arthur, Ermoliev and Kaniovski (1983) and Pemantle (1990) studied urn processes whose limit points depend on the trajectory.
  • Benaïm (1997) and Duflo (1997, Random Iterative Models) developed the attainability argument for general stochastic approximation processes.
  • Benaïm (1999) states it for arbitrary attractors of a semiflow on a locally compact metric space, under a single conditional shadowing condition (24).

Setting

Let (M,d)(M,d)(M,d) be a metric space and let Φ=(Φt)t≥0\Phi=(\Phi_t)_{t\ge0}Φ=(Φt​)t≥0​ be a semiflow on MMM: a continuous map (t,x)↦Φt(x)(t,x)\mapsto\Phi_t(x)(t,x)↦Φt​(x) with Φ0=Id\Phi_0=\mathrm{Id}Φ0​=Id and Φt+s=Φt∘Φs\Phi_{t+s}=\Phi_t\circ\Phi_sΦt+s​=Φt​∘Φs​.

  • A set AAA is invariant if Φt(A)=A\Phi_t(A)=AΦt​(A)=A for all t≥0t\ge0t≥0, and positively invariant if Φt(A)⊂A\Phi_t(A)\subset AΦt​(A)⊂A.
  • An attractor is a nonempty compact invariant set AAA with a neighbourhood WWW on which dist⁡(Φtx,A)→0\operatorname{dist}(\Phi_tx,A)\to0dist(Φt​x,A)→0 uniformly. Its basin B(A)B(A)B(A) is the set of points xxx with dist⁡(Φtx,A)→0\operatorname{dist}(\Phi_tx,A)\to0dist(Φt​x,A)→0.
  • A continuous curve X:R+→MX:\mathbb R_+\to MX:R+​→M is an asymptotic pseudotrajectory if sup⁡0≤h≤Td(X(t+h),Φh(X(t)))→0\sup_{0\le h\le T}d(X(t+h),\Phi_h(X(t)))\to0sup0≤h≤T​d(X(t+h),Φh​(X(t)))→0 as t→∞t\to\inftyt→∞, for every T>0T>0T>0.
  • The limit set of XXX is L(X)=⋂t≥0X([t,∞))‾L(X)=\bigcap_{t\ge0}\overline{X([t,\infty))}L(X)=⋂t≥0​X([t,∞))​.
  • For T>0T>0T>0, dX(T)=sup⁡k∈Nd(ΦT(X(kT)),X(kT+T))d_X(T)=\sup_{k\in\mathbb N}d(\Phi_T(X(kT)),X(kT+T))dX​(T)=supk∈N​d(ΦT​(X(kT)),X(kT+T)).

Now let X=(X(t))t≥0X=(X(t))_{t\ge0}X=(X(t))t≥0​ be a process on a probability space (Ω,F,P)(\Omega,\mathcal F,P)(Ω,F,P) with continuous paths in MMM, adapted to a filtration (Ft)t≥0(\mathcal F_t)_{t\ge0}(Ft​)t≥0​. The standing assumption of Section 7 is that for all δ>0\delta>0δ>0, T>0T>0T>0 and t≥0t\ge0t≥0,

P(sup⁡s≥t sup⁡0≤h≤Td(X(s+h),Φh(X(s)))≥δ ∣ Ft)≤w(t,δ,T)(24)P\Big(\sup_{s\ge t}\ \sup_{0\le h\le T}d\big(X(s+h),\Phi_h(X(s))\big)\ge\delta\ \Big|\ \mathcal F_t\Big)\le w(t,\delta,T)\tag{24}P(s≥tsup​ 0≤h≤Tsup​d(X(s+h),Φh​(X(s)))≥δ ​ Ft​)≤w(t,δ,T)(24)

for a function w≥0w\ge0w≥0 with w(t,δ,T)↓0w(t,\delta,T)\downarrow0w(t,δ,T)↓0 as t→∞t\to\inftyt→∞.

A point ppp is attainable if P(∃s≥t:X(s)∈U)>0P(\exists s\ge t: X(s)\in U)>0P(∃s≥t:X(s)∈U)>0 for every t>0t>0t>0 and every open neighbourhood UUU of ppp. Att(X)\mathrm{Att}(X)Att(X) is the set of attainable points.

Formalization targets

Goal: Theorem 7.3, first statement

If MMM is locally compact, AAA is an attractor of Φ\PhiΦ, and Att(X)∩B(A)≠∅\mathrm{Att}(X)\cap B(A)\neq\emptysetAtt(X)∩B(A)=∅, then

P(L(X)⊂A)>0.P\big(L(X)\subset A\big)>0 .P(L(X)⊂A)>0.

This statement contains no constants and no rates, so it does not depend on how (24) is quantified for a particular algorithm.

Theorem 7.3, second statement

If UUU is open and relatively compact with U‾⊂B(A)\overline U\subset B(A)U⊂B(A), there are T,δ>0T,\delta>0T,δ>0, depending only on UUU (and on Φ\PhiΦ, AAA), such that for every process satisfying the standing assumption and every t≥0t\ge0t≥0

P(L(X)⊂A)≥(1−w(t,δ,T)) P(∃s≥t: X(s)∈U).P\big(L(X)\subset A\big)\ge\big(1-w(t,\delta,T)\big)\,P\big(\exists s\ge t:\ X(s)\in U\big).P(L(X)⊂A)≥(1−w(t,δ,T))P(∃s≥t: X(s)∈U).

Milestones

  • Lemma 6.8. For a nonempty compact K⊂B(A)K\subset B(A)K⊂B(A) there are T,δ>0T,\delta>0T,δ>0 such that every asymptotic pseudotrajectory with X(0)∈KX(0)\in KX(0)∈K and dX(T)<δd_X(T)<\deltadX​(T)<δ has L(X)⊂AL(X)\subset AL(X)⊂A.
  • Lemma 7.1, in three parts:
    • Att(X)\mathrm{Att}(X)Att(X) is closed;
    • it is positively invariant;
    • it contains L(X)L(X)L(X) almost surely.

Significance

The result. Theorem 7.3 turns a question about the long-run limit of a random process into a question about where the process can go. Attainability is usually checked by a controllability argument: the noise can push the iterates in every direction. For urn processes with an urn function mapping the simplex into its interior, every point is attainable (Example 7.2 of the notes). Then every attractor of the mean ODE is reached with positive probability. Combined with nonconvergence results for unstable sets (Section 9 of the notes), this characterizes the possible limits of many learning and urn processes. Theorem 7.3 is the positive half of that picture.

Formalizing it. The theorem has a published proof (p. 32 of the notes) and no machine-checked version. A formal proof needs the following:

  • a precise reading of the conditional shadowing condition (24) as a conditional expectation of an indicator;
  • the stopping-time decomposition of the event {∃s≥t:X(s)∈U}\{\exists s\ge t: X(s)\in U\}{∃s≥t:X(s)∈U} over dyadic times;
  • the deterministic Lemma 6.8, which rests on the limit set theorem for precompact asymptotic pseudotrajectories (Theorem 5.7 of the notes, the subject of mission 1 of this series).

Difficulty

The obvious argument says: once XXX enters a compact part of the basin, the flow carries it into AAA. That fails because XXX is not a trajectory of the flow. Each window of length TTT adds an error, and errors over infinitely many windows can push the process out of the basin.

Two things are needed instead:

  • A uniform version of the deterministic statement, with TTT and δ\deltaδ fixed in advance from the compact set alone. This is Lemma 6.8, which needs local compactness of MMM and the structure of limit sets of asymptotic pseudotrajectories.
  • A probabilistic step that applies (24) at the random time when XXX first enters UUU. That time is not a stopping time on a continuum, and conditioning at it needs care.

A naive union bound over all times is useless: it does not use the conditional form of (24).

Formalization scope

  • The semiflow is Mathlib's Flow ℝ≥0 M on a metric space; local compactness is LocallyCompactSpace M.
  • The process is X : ℝ≥0 → Ω → M with continuous paths. The paper's alternative of càdlàg paths is not covered.
  • Adaptedness is Borel measurability of X(t)X(t)X(t) with respect to Ft\mathcal F_tFt​, for a Mathlib Filtration ℝ≥0. PPP is a probability measure.
  • The suprema in (24) and in dX(T)d_X(T)dX​(T) are computed in [0,∞][0,\infty][0,∞] with the extended distance, so that "sup ≥δ\ge\delta≥δ" and "sup <δ<\delta<δ" are exact even when the supremum is infinite or not attained.
  • The conditional probability in (24) is the conditional expectation of the indicator of the event. The event is required to be measurable, so the condition cannot hold vacuously through a junk conditional expectation.
  • www is required to be both nonincreasing in ttt and convergent to 000.
  • The events {L(X)⊂A}\{L(X)\subset A\}{L(X)⊂A} and {∃s≥t:X(s)∈U}\{\exists s\ge t: X(s)\in U\}{∃s≥t:X(s)∈U} are measured with PPP as an outer measure, so no measurability hypothesis is added for them.
  • In the second statement of Theorem 7.3, TTT and δ\deltaδ are chosen before the probability space, the process, www and ttt.
  • Invariance in the definition of an attractor is the equality Φt(A)=A\Phi_t(A)=AΦt​(A)=A, not inclusion.
  • The almost-sure clause of Lemma 7.1 is stated for separable MMM. Without separability it cannot be proved in ordinary set theory.

These choices rule out the trivializing formalizations: a conditional-probability hypothesis that holds vacuously, invariance read as inclusion, constants T,δT,\deltaT,δ that depend on the process or on ttt, and a probability bound www without monotonicity.

A complete development needs:

  • limit sets of asymptotic pseudotrajectories and the fact that an internally chain transitive set meeting the basin of an attractor lies in the attractor (shared with missions 1 and 5 of this series);
  • measurability of path functionals of continuous processes;
  • conditioning on events of the form {τ=tn(k)}\{\tau=t_n(k)\}{τ=tn​(k)} at dyadic times.

The first and second are reusable well beyond this mission. Contributions toward either, and alternative proofs of Lemma 6.8, are welcome.

Selected references

  • M. Benaïm, Dynamics of stochastic approximation algorithms, Séminaire de Probabilités XXXIII, Lecture Notes in Mathematics 1709, Springer, 1999, pp. 1–68. https://doi.org/10.1007/BFb0096509
  • M. Benaïm, M. W. Hirsch, Asymptotic pseudotrajectories and chain recurrent flows, with applications, Journal of Dynamics and Differential Equations 8 (1996), 141–176. https://doi.org/10.1007/BF02218617
  • M. Benaïm, Vertex-reinforced random walks and a conjecture of Pemantle, Annals of Probability 25 (1997), 361–392. https://doi.org/10.1214/aop/1024404292
  • M. Duflo, Random Iterative Models, Applications of Mathematics 34, Springer, 1997. https://doi.org/10.1007/978-3-662-12880-0
  • H. J. Kushner, D. S. Clark, Stochastic Approximation Methods for Constrained and Unconstrained Systems, Springer, 1978. https://doi.org/10.1007/978-1-4684-9352-8
  • C. Conley, Isolated Invariant Sets and the Morse Index, CBMS Regional Conference Series in Mathematics 38, AMS, 1978. https://doi.org/10.1090/cbms/038
11 thms2 active usersReviewed
Bandit AlgorithmsMachine LearningOperations Research+1·Captain: mikedeng1

Analysis of Thompson Sampling for the Multi-armed Bandit Problem 2: Logarithmic Regret for N ArmsResearch Paper

Motivation

Thompson Sampling is the oldest heuristic for the multi-armed bandit problem: proposed by Thompson in 1933, it plays each arm with the posterior probability that the arm is the best one. It is simple to implement, performs well empirically (Chapelle and Li, NIPS 2011), and has been used in production systems such as click-through-rate prediction for search advertising. For a long time, however, no finite-time regret guarantee was known for it: the analyses available before 2012 gave only o(T)o(T)o(T) regret.

Agrawal and Goyal (arXiv:1111.1797, COLT 2012) gave the first logarithmic bounds on the expected regret of Thompson Sampling. This mission formalizes their bound for the general case of NNN arms (their Theorem 2). A companion mission of the same series formalizes their two-armed bound (Theorem 1), whose proof is independent.

Timeline. Lai and Robbins (1985) proved that every consistent algorithm has regret at least of order ∑iΔiD(μi∥μ1)ln⁡T\sum_i \frac{\Delta_i}{D(\mu_i\|\mu_1)}\ln T∑i​D(μi​∥μ1​)Δi​​lnT. Auer, Cesa-Bianchi and Fischer (2002) showed that UCB1 achieves O(∑iln⁡T/Δi)O(\sum_i \ln T/\Delta_i)O(∑i​lnT/Δi​) in finite time. Agrawal and Goyal (2012) proved O((∑a1/Δa2)2ln⁡T)O((\sum_a 1/\Delta_a^2)^2\ln T)O((∑a​1/Δa2​)2lnT) for Thompson Sampling with NNN arms; Kaufmann, Korda and Munos (2012) and Agrawal and Goyal (2013) later proved the asymptotically optimal constant for Bernoulli rewards.

Setting

A stochastic NNN-armed bandit has arms 1,…,N1,\dots,N1,…,N. Arm iii, when played, yields a random reward drawn from a fixed distribution νi\nu_iνi​ supported in [0,1][0,1][0,1], with mean μi\mu_iμi​; rewards of an arm are i.i.d. and independent of the other arms. Arm 111 is assumed to be the unique optimal arm, μ1>μi\mu_1>\mu_iμ1​>μi​ for i≠1i\ne1i=1, and Δi=μ1−μi>0\Delta_i=\mu_1-\mu_i>0Δi​=μ1​−μi​>0 is the gap of arm iii.

Thompson Sampling for general stochastic bandits (Algorithm 2 of the paper) keeps, for each arm iii, a count SiS_iSi​ of successes and FiF_iFi​ of failures, both starting at 000. In each round ttt it draws θi(t)∼Beta(Si+1,Fi+1)\theta_i(t)\sim\mathrm{Beta}(S_i+1,F_i+1)θi​(t)∼Beta(Si​+1,Fi​+1) independently for every arm, plays i(t)=arg⁡max⁡iθi(t)i(t)=\arg\max_i\theta_i(t)i(t)=argmaxi​θi​(t), observes a reward r~t∼νi(t)\tilde r_t\sim\nu_{i(t)}r~t​∼νi(t)​, performs a Bernoulli trial with success probability r~t\tilde r_tr~t​, and increments Si(t)S_{i(t)}Si(t)​ on success and Fi(t)F_{i(t)}Fi(t)​ on failure.

The expected regret in time TTT is

E[R(T)]=E[∑t=1T(μ∗−μi(t))],μ∗=max⁡iμi,\mathbb E[\mathcal R(T)]=\mathbb E\Big[\sum_{t=1}^T(\mu^*-\mu_{i(t)})\Big],\qquad \mu^*=\max_i\mu_i,E[R(T)]=E[t=1∑T​(μ∗−μi(t)​)],μ∗=imax​μi​,

the expectation being over the rewards, the Bernoulli trials and the posterior samples.

The proof works with the reward stacks Zi,mZ_{i,m}Zi,m​: the outcome of the mmm-th Bernoulli trial of arm iii, all independent. Then s(j)=∑m≤jZ1,ms(j)=\sum_{m\le j}Z_{1,m}s(j)=∑m≤j​Z1,m​, the number of successes in the first jjj plays of arm 111, is a Binomial(j,μ1)\mathrm{Binomial}(j,\mu_1)Binomial(j,μ1​) random variable. The other objects of the proof are the threshold Li=24ln⁡T/Δi2L_i=24\ln T/\Delta_i^2Li​=24lnT/Δi2​, the saturated set C(t)C(t)C(t) of suboptimal arms with at least LiL_iLi​ plays before round ttt, the intervals IjI_jIj​ between the jjj-th and (j+1)(j+1)(j+1)-th plays of arm 111, and the counts γj\gamma_jγj​ and Vjℓ,aV_j^{\ell,a}Vjℓ,a​ defined in §4.

Formalization targets

Goal: Theorem 2

There is an absolute constant C>0C>0C>0 such that for every N≥2N\ge2N≥2, every instance as above and every horizon T≥2T\ge2T≥2,

E[R(T)]≤C(∑a=2N1Δa2)2ln⁡T.\mathbb E[\mathcal R(T)]\le C\Big(\sum_{a=2}^N\frac{1}{\Delta_a^2}\Big)^2\ln T .E[R(T)]≤C(a=2∑N​Δa2​1​)2lnT.

CCC does not depend on NNN, on the reward distributions or on TTT.

Milestones

  1. Lemma 4: with E(t)E(t)E(t) the event that every saturated arm's sample lies within Δi/2\Delta_i/2Δi​/2 of its mean, Pr⁡(E(t))≥1−4(N−1)/T2\Pr(E(t))\ge1-4(N-1)/T^2Pr(E(t))≥1−4(N−1)/T2, also conditionally on s(j)=ss(j)=ss(j)=s.
  2. Lemma 5 (Eq. (7)): the expected regret from saturated arms inside IjI_jIj​ is at most E[E[γj+1∣s(j)]∑aΔaE[min⁡{X(j,s(j),μa+Δa/2),T}∣s(j)]]\mathbb E\big[\mathbb E[\gamma_j+1\mid s(j)]\sum_a\Delta_a\mathbb E[\min\{X(j,s(j),\mu_a+\Delta_a/2),T\}\mid s(j)]\big]E[E[γj​+1∣s(j)]∑a​Δa​E[min{X(j,s(j),μa​+Δa​/2),T}∣s(j)]].
  3. Lemma 1: E[X(j,s,y)]=1/Fj+1,yB(s)−1\mathbb E[X(j,s,y)]=1/F^B_{j+1,y}(s)-1E[X(j,s,y)]=1/Fj+1,yB​(s)−1, where X(j,s,y)X(j,s,y)X(j,s,y) counts the trials before an independent Beta(s+1,j−s+1)\mathrm{Beta}(s+1,j-s+1)Beta(s+1,j−s+1) sample exceeds yyy.
  4. Lemma 3: a three-case bound on E[E[min⁡{X(j,s(j),y),T}∣s(j)]]\mathbb E[\mathbb E[\min\{X(j,s(j),y),T\}\mid s(j)]]E[E[min{X(j,s(j),y),T}∣s(j)]] in terms of the Bernoulli KL divergence DDD between yyy and μ1\mu_1μ1​.

Significance

The result. Theorem 2 shows that Thompson Sampling, a randomized Bayesian heuristic, achieves regret logarithmic in the horizon for any number of arms with bounded rewards, matching the order in TTT of the Lai–Robbins lower bound. Its dependence on the gaps, (∑aΔa−2)2(\sum_a\Delta_a^{-2})^2(∑a​Δa−2​)2, is worse than UCB1's; the paper's own Remark 1 and later work improve it. The proof introduced the device of bounding the waiting time between plays of the optimal arm through geometric variables with Beta-cdf parameters (Lemmas 1 and 3), which reappears in later analyses of Thompson Sampling.

Formalizing it. The theorem is proved on paper; it has not been machine-checked. Bandit theory in Lean (bandit environments, regret, UCB-type analyses) is still young, and no Beta–Bernoulli Thompson Sampling result is formalized. The mission produces a Lean model of Algorithm 2 for general [0,1][0,1][0,1] rewards with the paper's stack coupling, the §4 bookkeeping of saturated arms and intervals, and the paper's lemmas as separate targets.

Difficulty

Two difficulties are specific to the NNN-armed analysis. First, the arm that competes with arm 111 changes over time: the set of saturated arms grows, and which saturated arm is "best" depends on the history, so the waiting time between plays of arm 111 cannot be compared with a single geometric variable as in the two-armed case. Second, the number γj\gamma_jγj​ of rounds at which arm 111's sample is large but arm 111 is not played is not independent of the counts Vjℓ,aV_j^{\ell,a}Vjℓ,a​: both depend on the same posterior samples, and Lemma 5 needs a careful conditioning on the history to separate them. The obvious union bound over arms, treating each suboptimal arm as in the two-armed proof, fails because it ignores the interruptions by unsaturated arms, whose number is the source of the squared sum in the bound.

Formalization scope

  • Probability space. Algorithm 2 is realized on a product of three independent i.i.d. tables: Beta draws indexed by (arm, round, successes, failures), rewards indexed by (arm, round) and uniform variables indexed by (arm, round); the Bernoulli trial of a round succeeds when the played arm's uniform variable is below its reward. The law of the run is that of Algorithm 2, which runs for every round t=1,2,…t=1,2,\dotst=1,2,…. s(j)s(j)s(j) is the number of successful trials among the first jjj plays of arm 111 in this infinite run (possibly after the horizon TTT), so it is a Binomial(j,μ1)\mathrm{Binomial}(j,\mu_1)Binomial(j,μ1​) random variable for every jjj, as the paper's independent Z1,mZ_{1,m}Z1,m​ make it. Ties in the arg max (probability 000) go to the smallest index.
  • Indexing. Arms are Fin N, and Lean arm 0 is the paper's arm 111. Rounds are 0,…,T−10,\dots,T-10,…,T−1; Lean round ttt is the paper's round t+1t+1t+1. Sums over a=2,…,Na=2,\dots,Na=2,…,N are sums over a≠0a\ne0a=0.
  • Expectations are lower Lebesgue integrals in [0,∞][0,\infty][0,∞], which has no junk value for non-integrable functions. Conditional expectations given s(j)s(j)s(j) are written as finite sums over the values of s(j)s(j)s(j).
  • The O(⋅)O(\cdot)O(⋅). The paper writes O(⋅)O(\cdot)O(⋅) in the sense of its footnote 1 (f(n)≤c g(n)f(n)\le c\,g(n)f(n)≤cg(n) for n≥n0n\ge n_0n≥n0​). The goal states it with one universal constant CCC, quantified before NNN, the instance and TTT, for every T≥2T\ge2T≥2. The explicit constants printed in App. D are not formalized: expanding the paper's Eq. (21) gives terms 288(N−1)(ln⁡T)∑aΔa−2288(N-1)(\ln T)\sum_a\Delta_a^{-2}288(N−1)(lnT)∑a​Δa−2​ and 48(N−1)248(N-1)^248(N−1)2 where the paper prints 288(ln⁡T)∑iΔi−2288(\ln T)\sum_i\Delta_i^{-2}288(lnT)∑i​Δi−2​, and Eq. (22) drops a factor ln⁡T\ln TlnT in its 192/Δa2192/\Delta_a^2192/Δa2​ term. The O(⋅)O(\cdot)O(⋅) claim does not depend on these slips; a statement pinned to the printed numerals might be false.
  • Ruled out. A constant depending on NNN, on the means or on TTT; a fixed number of arms; Bernoulli rewards only; or any algorithm other than Algorithm 2 would each make the goal a different and weaker theorem. The statement quantifies over all N≥2N\ge2N≥2 and all reward distributions on [0,1][0,1][0,1].
  • Not included. Eq. (8), the bound ∑jE[γj∣s(j)]≤∑uLu+4(N−1)\sum_{j}\mathbb E[\gamma_j\mid s(j)]\le\sum_uL_u+4(N-1)∑j​E[γj​∣s(j)]≤∑u​Lu​+4(N−1) "for all instantiations", is not a milestone: each term is conditioned on a different s(j)s(j)s(j), and the pointwise reading does not follow from the argument given. Remark 1 (an alternate bound) and App. A (several optimal arms) are not part of this mission.
  • Contributions welcome: Beta–Binomial identities, geometric waiting times, Hoeffding bounds for binomial cdfs, and the stopping-time arguments behind Lemma 5. Lemma 1 and Lemma 3 are shared with the two-armed mission of this series.

Selected references

  • S. Agrawal and N. Goyal, Analysis of Thompson Sampling for the Multi-armed Bandit Problem, COLT 2012; arXiv:1111.1797v3, 2012. https://arxiv.org/abs/1111.1797
  • W. R. Thompson, On the likelihood that one unknown probability exceeds another in view of the evidence of two samples, Biometrika 25, 1933. https://doi.org/10.1093/biomet/25.3-4.285
  • T. L. Lai and H. Robbins, Asymptotically efficient adaptive allocation rules, Advances in Applied Mathematics 6, 1985. https://doi.org/10.1016/0196-8858(85)90002-8
  • P. Auer, N. Cesa-Bianchi and P. Fischer, Finite-time analysis of the multiarmed bandit problem, Machine Learning 47, 2002. https://doi.org/10.1023/A:1013689704352
  • O. Chapelle and L. Li, An empirical evaluation of Thompson Sampling, NIPS 2011. https://papers.nips.cc/paper/4321-an-empirical-evaluation-of-thompson-sampling
  • E. Kaufmann, N. Korda and R. Munos, Thompson Sampling: an asymptotically optimal finite-time analysis, ALT 2012. https://arxiv.org/abs/1205.4217
  • S. Agrawal and N. Goyal, Further optimal regret bounds for Thompson Sampling, AISTATS 2013. https://arxiv.org/abs/1209.3353
11 thms3 active usersReviewed
PreviousNext

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