Prove2Me
Navigate
DiscoverCollectionsFormalpediaBlogsUsersMomentumMy Missions+
Prove2Me
⌕
Log in

Get started

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

Probability

1,628 missions · 697 completed

Missions

Open931Completed697All1628
Algorithmic Game TheoryMechanism DesignOperations Research·Captain: mikedeng1

Approximate Revenue Maximization with Multiple Items 6: With n Buyers and Dominant-Strategy Implementation, Selling Two Independent Goods Separately Earns at Least Half the Optimal RevenueResearch Paper

Motivation

A seller who owns several goods and faces buyers with private values must decide how to sell them. The revenue-maximizing mechanism is known in closed form only for one good: for a single buyer it is a posted price, and for several buyers it is Myerson's optimal auction (Myerson 1981). With two or more goods the optimal mechanism may use lotteries, may be impossible to describe compactly, and may behave non-monotonically in the value distribution. This makes simple mechanisms with provable revenue guarantees the practical question: how much of the optimal revenue does a seller lose by selling each good separately with its own optimal one-good mechanism?

Hart and Nisan (arXiv:1204.1846v3) answer this for one buyer and two independent goods (their Theorem A): selling separately earns at least half of the optimal revenue. Their Theorem 18 (p. 29), proved in Appendix A.7 as Theorems 33 and 34, extends the factor 12\tfrac1221​ to any number nnn of buyers. This mission formalizes the dominant-strategy half, Theorem 33, in which the buyers may be arbitrarily correlated with one another, as long as the two goods are independent.

Setting

One seller sells kkk goods to n≥1n \ge 1n≥1 buyers. Buyer jjj's value for good iii is a number xij≥0x^j_i \ge 0xij​≥0; a valuation profile is x=(xij)∈R+knx = (x^j_i) \in \mathbb R^{kn}_+x=(xij​)∈R+kn​, and xj=(xij)ix^j = (x^j_i)_ixj=(xij​)i​ is buyer jjj's valuation vector. The profile is random, XXX, with a commonly known joint law.

A mechanism μ=(qj,sj)j=1,…,n\mu = (q^j, s^j)_{j=1,\dots,n}μ=(qj,sj)j=1,…,n​ gives each buyer jjj an allocation qj(x)∈[0,1]kq^j(x) \in [0,1]^kqj(x)∈[0,1]k (the probability that jjj receives each good) and a payment sj(x)∈Rs^j(x) \in \mathbb Rsj(x)∈R, as functions of the reported profile, subject to feasibility: ∑jqij(x)≤1\sum_j q^j_i(x) \le 1∑j​qij​(x)≤1 for every good iii. Buyer jjj's payoff is bj(x)=qj(x)⋅xj−sj(x)b^j(x) = q^j(x)\cdot x^j - s^j(x)bj(x)=qj(x)⋅xj−sj(x), and the seller's revenue is S(x)=∑jsj(x)S(x) = \sum_j s^j(x)S(x)=∑j​sj(x).

  • IC-DS (dominant-strategy incentive compatibility): bj(x)≥qj(x~j,x−j)⋅xj−sj(x~j,x−j)b^j(x) \ge q^j(\tilde x^j, x^{-j})\cdot x^j - s^j(\tilde x^j, x^{-j})bj(x)≥qj(x~j,x−j)⋅xj−sj(x~j,x−j) for every buyer jjj, every profile xxx and every report x~j\tilde x^jx~j, where (x~j,x−j)(\tilde x^j, x^{-j})(x~j,x−j) replaces jjj's vector in xxx by x~j\tilde x^jx~j.
  • IR-DS: bj(x)≥0b^j(x) \ge 0bj(x)≥0 for every jjj and xxx.
  • NPT (no positive transfer): sj≥0s^j \ge 0sj≥0.

The expected revenue is R(μ;X)=E[S(X)]R(\mu; X) = \mathbb E[S(X)]R(μ;X)=E[S(X)], and

RevDS(X)=sup⁡{R(μ;X):μ feasible, IC-DS and IR-DS}∈[0,∞].\mathrm{Rev}^{DS}(X) = \sup\{R(\mu;X) : \mu \text{ feasible, IC-DS and IR-DS}\} \in [0, \infty].RevDS(X)=sup{R(μ;X):μ feasible, IC-DS and IR-DS}∈[0,∞].

For two goods write Y=X1=(X1j)jY = X_1 = (X^j_1)_jY=X1​=(X1j​)j​ and Z=X2=(X2j)jZ = X_2 = (X^j_2)_jZ=X2​=(X2j​)j​ for the vectors of the buyers' values for good 1 and good 2, so X=(Y,Z)X = (Y, Z)X=(Y,Z). RevDS(Y)\mathrm{Rev}^{DS}(Y)RevDS(Y) is the optimal dominant-strategy revenue from good 1 alone, sold to the same nnn buyers. Finally a(1)=max⁡jaja^{(1)} = \max_j a^ja(1)=maxj​aj.

Formalization targets

Goal: Theorem 33 (p. 49)

If YYY and ZZZ are independent, then

RevDS(X1)+RevDS(X2)≥12 RevDS(X1,X2).\mathrm{Rev}^{DS}(X_1) + \mathrm{Rev}^{DS}(X_2) \ge \tfrac12\,\mathrm{Rev}^{DS}(X_1, X_2).RevDS(X1​)+RevDS(X2​)≥21​RevDS(X1​,X2​).

No independence among buyers is assumed: the values Xi1,…,XinX^1_i, \dots, X^n_iXi1​,…,Xin​ for the same good may have any joint law.

Milestones (the steps of the paper's proof)

  1. Remark (b) of A.7 (p. 48): an IC-DS, IR-DS mechanism is NPT iff sj(0,x−j)=0s^j(0, x^{-j}) = 0sj(0,x−j)=0 for all jjj, x−jx^{-j}x−j, and NPT is without loss of generality for RevDS\mathrm{Rev}^{DS}RevDS.
  2. Remark (b), subdomain property (p. 48): E[∑jsj(X)1X∈A]≤RevDS(X1X∈A)≤RevDS(X)\mathbb E[\sum_j s^j(X)\mathbf 1_{X\in A}] \le \mathrm{Rev}^{DS}(X\mathbf 1_{X\in A}) \le \mathrm{Rev}^{DS}(X)E[∑j​sj(X)1X∈A​]≤RevDS(X1X∈A​)≤RevDS(X).
  3. Display (21) (p. 50): z(1) P[Y(1)≥z(1)]≤RevDS(Y)z^{(1)}\,\mathbb P[Y^{(1)} \ge z^{(1)}] \le \mathrm{Rev}^{DS}(Y)z(1)P[Y(1)≥z(1)]≤RevDS(Y).
  4. The marginal mechanism (p. 50): for fixed second-good values zzz, q^j(y)=q1j(y,z)\hat q^j(y) = q^j_1(y,z)q^​j(y)=q1j​(y,z), s^j(y)=sj(y,z)−q2j(y,z)zj\hat s^j(y) = s^j(y,z) - q^j_2(y,z)z^js^j(y)=sj(y,z)−q2j​(y,z)zj is feasible, IC-DS and IR-DS.
  5. Display (19) (p. 49): E[S(Y,Z)1Y(1)≥Z(1)]≤2 RevDS(Y)\mathbb E[S(Y,Z)\mathbf 1_{Y^{(1)} \ge Z^{(1)}}] \le 2\,\mathrm{Rev}^{DS}(Y)E[S(Y,Z)1Y(1)≥Z(1)​]≤2RevDS(Y) for an NPT, IC-DS, IR-DS mechanism.

Companion: Theorem 34 (p. 49)

With nnn independent buyers and two independent goods, so that all 2n2n2n values are independent (footnote 27), the same inequality holds for Bayesian Nash implementation: RevBN(X1)+RevBN(X2)≥12RevBN(X1,X2)\mathrm{Rev}^{BN}(X_1) + \mathrm{Rev}^{BN}(X_2) \ge \tfrac12\mathrm{Rev}^{BN}(X_1, X_2)RevBN(X1​)+RevBN(X2​)≥21​RevBN(X1​,X2​).

Significance

The result. Theorem 33 shows that the factor-12\tfrac1221​ guarantee of separate selling does not depend on there being a single buyer, nor on the buyers being independent, provided incentive compatibility is required in dominant strategies. It reduces a two-good, nnn-buyer design problem, whose optimum is not characterized, to two one-good problems, each of which can be handled by known single-good auction theory. Remark (a) on p. 29 records that the dominant-strategy statement holds for correlated buyers, while Remark (b) there explains, via Crémer and McLean (1988), why the Bayesian Nash analogue cannot be extended to correlated buyers.

Formalizing it. The result is proved in the paper; to our knowledge no machine-checked version exists. A formal proof requires multi-agent mechanisms on R+kn\mathbb R^{kn}_+R+kn​, the Lebesgue integral of a possibly non-integrable revenue, conditioning on an independent coordinate of a product law, and the subdomain property of Remark (b). The definitions here (feasible multi-buyer mechanisms, IC-DS, IR-DS, RevDS\mathrm{Rev}^{DS}RevDS, and the Bayesian Nash counterparts) are reusable for any multi-buyer revenue-maximization result.

Difficulty

The obvious first idea is to apply the one-buyer Theorem A buyer by buyer. It fails: the feasibility constraint ∑jqij≤1\sum_j q^j_i \le 1∑j​qij​≤1 couples the buyers, so a multi-buyer mechanism does not decompose into single-buyer mechanisms, and the revenue extracted from one buyer depends on the reports of all the others. Any reduction to one-good problems must produce, from a two-good mechanism, one-good mechanisms that remain feasible and dominant-strategy incentive compatible for every buyer at once. On the measure-theoretic side, payments of an incentive compatible mechanism may be negative and non-integrable, so expectations must be handled in the extended reals, and the passage from bounds that hold for each fixed value of one good to a bound in expectation needs the product structure of the law of (Y,Z)(Y, Z)(Y,Z).

Formalization scope

  • A profile is x : Fin n → ι → ℝ≥0 (x j i =xij= x^j_i=xij​); a random profile is given by its law, a probability measure on that space. One good is ι = Unit, two goods are ι = Fin 2 (good 1 is index 0). The law of two independent goods is twoGoodsN μY μZ, the image of the product law μY⊗μZ\mu_Y\otimes\mu_ZμY​⊗μZ​; μY and μZ are arbitrary probability laws on R+n\mathbb R^n_+R+n​, so correlation among buyers is allowed.
  • (x~j,x−j)(\tilde x^j, x^{-j})(x~j,x−j) is Function.update x j x'. Feasibility ∑jqij≤1\sum_j q^j_i \le 1∑j​qij​≤1 is part of admissibility; without it a mechanism could sell each good to every buyer.
  • Payments are required measurable (the paper's footnote 12). R(μ;X)R(\mu;X)R(μ;X) is the EReal difference of the integrals of S+S^+S+ and S−S^-S−, never a Bochner integral, which would vanish on non-integrable revenues. RevDS\mathrm{Rev}^{DS}RevDS takes values in [0,∞][0,\infty][0,∞].
  • The subdomain property requires the set AAA to be measurable, so that the law of X1X∈AX\mathbf 1_{X\in A}X1X∈A​ is defined. It also requires the signed expectation on AAA to be defined: at least one of its positive and negative part integrals must be finite. The hypothesis n≥1n \ge 1n≥1 is the paper's.
  • For Theorem 34, buyers are independent with laws ρj\rho_jρj​; interim payoffs integrate over the others' product law, and allocations as well as payments are required measurable.
  • The constant 12\tfrac1221​ is stated exactly as printed. A statement over an arbitrary joint law with given marginals, or with the feasibility constraint dropped, would be a different (and false, or trivial) theorem.

Contributions of every kind are welcome: proofs of the milestones, of Theorem 34, and general lemmas on multi-buyer mechanisms.

Selected references

  • S. Hart and N. Nisan, Approximate Revenue Maximization with Multiple Items, arXiv:1204.1846v3, 2017; J. Economic Theory 172 (2017). https://arxiv.org/abs/1204.1846
  • R. B. Myerson, Optimal Auction Design, Mathematics of Operations Research 6 (1981), 58–73. https://doi.org/10.1287/moor.6.1.58
  • J. Crémer and R. P. McLean, Full Extraction of the Surplus in Bayesian and Dominant Strategy Auctions, Econometrica 56 (1988), 1247–1257. https://doi.org/10.2307/1913096
  • S. Hart and P. J. Reny, Maximal Revenue with Multiple Goods: Nonmonotonicity and Other Observations, Theoretical Economics 10 (2015), 893–922. https://doi.org/10.3982/TE1517
8 thms1 active userReviewed
Algorithmic Game TheoryMechanism DesignOperations Research·Captain: mikedeng1

Approximate Revenue Maximization with Multiple Items 3: Selling k ≥ 2 Independent Goods Separately Guarantees a Fraction c/log² k of the Optimal RevenueResearch Paper

Motivation

A seller who knows the distributions of a buyer's values for several goods can, in principle, design a mechanism that offers lotteries and prices conditional on the buyer's reported values. Finding the best such mechanism can be difficult even with two goods. Separate selling is much simpler: choose a one-good price for each item and let the buyer decide item by item. The question is how much revenue this simplicity can cost when the goods' values are independent. Hart and Nisan establish a guarantee that declines only with the square of the logarithm of the number of goods, even though the distributions need not be identical or bounded. Their result is Theorem C of Approximate Revenue Maximization with Multiple Items, stated on p. 7 and proved on p. 26.

The paper grew from the problem of comparing easy selling formats with fully optimal mechanisms. For one good, a posted price characterizes optimal revenue. With multiple goods, the choice of a mechanism becomes genuinely multidimensional: its allocation rule can correlate the disposition of the goods and use lotteries. The guarantee in this mission applies uniformly over all independent distributions on nonnegative valuations, so a seller does not need a regularity, density, or finite-moment assumption to use it.

Setting

There is one seller, one risk-neutral buyer, and k≥2k\ge2k≥2 goods. The buyer's valuation vector is x=(xi)i=1k∈R+kx=(x_i)_{i=1}^k\in\mathbb R_+^kx=(xi​)i=1k​∈R+k​; receiving a set of goods gives the sum of the values of its members. The seller knows the joint law of the random vector XXX, but not its realization. In this mission the coordinates XiX_iXi​ are independent, with possibly different probability laws νi\nu_iνi​.

A direct mechanism consists of an allocation rule qi(x)∈[0,1]q_i(x)\in[0,1]qi​(x)∈[0,1] for each good and a real payment s(x)s(x)s(x). At a truthful report xxx, the buyer's payoff is b(x)=∑iqi(x)xi−s(x)b(x)=\sum_i q_i(x)x_i-s(x)b(x)=∑i​qi​(x)xi​−s(x). The mechanism is incentive compatible if truthful reporting is at least as good as every alternative report, and individually rational if b(x)≥0b(x)\ge0b(x)≥0 for every xxx. The seller's expected payment is its revenue. The optimal revenue Rev⁡(X)\operatorname{Rev}(X)Rev(X) is the supremum of revenue over such mechanisms. Payments are taken to be measurable, following footnote 12 on p. 12 of Hart and Nisan.

For one good, Rev⁡(Xi)\operatorname{Rev}(X_i)Rev(Xi​) equals the supremum of pPr⁡[Xi≥p]p\Pr[X_i\ge p]pPr[Xi​≥p] over prices p≥0p\ge0p≥0, as in display (1) on p. 12. The separate revenue is SRev⁡(X)=∑iRev⁡(Xi)\operatorname{SRev}(X)=\sum_i\operatorname{Rev}(X_i)SRev(X)=∑i​Rev(Xi​). The bundled revenue is BRev⁡(X)=Rev⁡(∑iXi)\operatorname{BRev}(X)=\operatorname{Rev}(\sum_i X_i)BRev(X)=Rev(∑i​Xi​), where the sum is treated as one good. Both formats are defined through one-good revenue, even when the full optimum can use multidimensional allocations.

Formalization targets

Theorem C: a uniform guarantee for separate selling

The goal is the inequality form of Theorem C, with one constant that works for every number and every family of independent goods:

∃c>0  ∀k≥2  ∀(νi)i=1k,SRev⁡(X1,…,Xk)≥c(log⁡k)2Rev⁡(X1,…,Xk).\exists c>0\;\forall k\ge2\;\forall (\nu_i)_{i=1}^k,\qquad \operatorname{SRev}(X_1,\ldots,X_k) \ge \frac{c}{(\log k)^2}\operatorname{Rev}(X_1,\ldots,X_k).∃c>0∀k≥2∀(νi​)i=1k​,SRev(X1​,…,Xk​)≥(logk)2c​Rev(X1​,…,Xk​).

Here XiX_iXi​ has law νi\nu_iνi​, and the joint law is their product. The paper phrases the headline as a guarantee for an approximation ratio; the displayed inequality is the form used in its proof on p. 26. It handles zero and infinite optimal revenues without introducing a ratio with an undefined denominator.

Supporting targets

The milestone list records the one-good posted-price formula, monotonicity under first-order stochastic domination, the equal-revenue comparison of Lemma 20, Proposition 13(ii)'s comparison of separate and bundled revenue, the two inequalities of Theorem 7, and the two-good base case of Theorem A. It also records the power-of-two estimate and the zero-padding assertion stated in the proof of Theorem C. Each milestone is attached to its source page in the pinned preprint.

Significance

The theorem gives a distribution-independent performance floor for a mechanism that can be run using one price per good. It applies to unequal and unbounded independent distributions, for which expected values can be infinite while one-good optimal revenues remain finite. The comparison is with the best incentive-compatible and individually rational multi-good mechanism, including randomized allocations, rather than with another simple selling format. The result is proved in the 2017 preprint; the open work here is a machine-checked formalization of that known result and its selected supporting statements.

The formal development also supplies reusable objects: a law-based direct-mechanism model, extended-real expected revenue, product laws for independent groups of goods, and first-order stochastic domination of one-good laws. These are useful for other comparisons among separate selling, bundling, and optimal mechanisms in the same paper. The mission does not pose the paper's distinct i.i.d. bundling theorem, multi-buyer results, or tightness examples.

Difficulty

The main obstacle is that optimal multi-good revenue is not monotone in the values of the goods. Increasing each coordinate of a valuation vector does not, by itself, give a valid comparison between the optimal revenues of the two vectors' laws; Hart and Nisan, p. 23 explicitly warn against that inference. One-good revenue does have the needed monotonicity, but the full optimum can depend on how the goods interact in an incentive-compatible mechanism. Theorem 7 is a substantive bridge between these two regimes. Infinite expected values also make a real-valued integral or an ordinary ratio unsuitable for the general statement.

Formalization scope

Goods are indexed by finite types, and valuations lie in R≥0k\mathbb R_{\ge0}^kR≥0k​. Each good's law is a probability measure; independence means the product law, while Theorem 7 allows dependence among coordinates within each of its two independent groups. A one-good law is represented by a singleton coordinate type. Revenue takes values in R≥0∪{∞}\mathbb R_{\ge0}\cup\{\infty\}R≥0​∪{∞}, and the expectation of a mechanism's possibly negative payment is represented in extended reals. The supremum defining optimal revenue ranges only over feasible, incentive-compatible, individually rational mechanisms with measurable payments. Theorem C quantifies ccc before kkk and the laws, and the restriction k≥2k\ge2k≥2 keeps the logarithmic coefficient in its stated domain.

The equal-revenue law is represented by the density x−2x^{-2}x−2 on [1,∞)[1,\infty)[1,∞), mapped to nonnegative reals. First-order stochastic domination compares the upper-tail probabilities at nonnegative thresholds, which is equivalent to testing every real threshold for nonnegative valuations. A complete proof will need measure-theoretic facts about that law, the posted-price characterization, product measures, and operations on extended nonnegative revenue. Contributions to those reusable facts and to the selected milestone theorems are within scope.

Selected references

  • Sergiu Hart and Noam Nisan, Approximate Revenue Maximization with Multiple Items, arXiv:1204.1846v3, 2017. Preprint.
10 thms1 active userReviewed
Algorithmic Game TheoryMechanism DesignOperations Research·Captain: mikedeng1

Approximate Revenue Maximization with Multiple Items 5: Bundling Is Optimal for Two I.I.D. Goods Whose Density on [a, ∞) Satisfies x f′(x) + (3/2) f(x) ≤ 0Research Paper

Motivation

A seller with several goods and a single buyer of unknown, additive valuation faces the multi-dimensional revenue maximization problem. For one good the answer is classical: by Myerson's theorem (Myerson 1981) a single posted price is optimal. For two or more goods no such simplification exists in general: optimal mechanisms may randomize, may offer menus of lotteries, and are rarely known in closed form. Selling the grand bundle at one price is the simplest multi-good mechanism, and it is used widely in practice.

Hart and Nisan (arXiv:1204.1846, J. Econ. Theory 2017) quantify how much revenue simple mechanisms lose. Most of their results are approximation bounds. Section 7.2 of the paper also contains an exact result: for two independent and identically distributed goods whose density decays fast enough, in the sense of a differential inequality, bundling is exactly optimal. This mission formalizes that theorem, Theorem 16, and the steps of its proof in Appendix A.6.

Timeline. Manelli and Vincent (2006) gave conditions under which a bundled mechanism is optimal for a multiple-good monopolist. Hart and Nisan's preprint (2012, revised 2017) proved Theorem 16 with the condition xf′(x)+32f(x)≤0x f'(x) + \tfrac32 f(x) \le 0xf′(x)+23​f(x)≤0. Daskalakis, Deckelbaum and Tzamos (2017) later gave a strong duality theory for the multiple-good monopolist, based on optimal transport.

Setting

There are two goods and one buyer. A valuation is a vector x=(y,z)∈R+2x = (y, z) \in \mathbb{R}^2_+x=(y,z)∈R+2​. The buyer's values Y,ZY, ZY,Z are i.i.d. with a one-good law FFF on R+\mathbb{R}_+R+​.

A mechanism μ=(q,s)\mu = (q, s)μ=(q,s) assigns to each reported valuation xxx an allocation vector q(x)∈[0,1]2q(x) \in [0,1]^2q(x)∈[0,1]2 (qi(x)q_i(x)qi​(x) is the probability that good iii is sold) and a payment s(x)∈Rs(x) \in \mathbb{R}s(x)∈R. The buyer's payoff is b(x)=q(x)⋅x−s(x)b(x) = q(x)\cdot x - s(x)b(x)=q(x)⋅x−s(x). The mechanism is

  • incentive compatible (IC) if b(x)≥q(x~)⋅x−s(x~)b(x) \ge q(\tilde x)\cdot x - s(\tilde x)b(x)≥q(x~)⋅x−s(x~) for all x,x~∈R+2x, \tilde x \in \mathbb{R}^2_+x,x~∈R+2​;
  • individually rational (IR) if b(x)≥0b(x) \ge 0b(x)≥0 for all xxx;
  • NPT (no positive transfer) if s≥0s \ge 0s≥0.

The revenue of μ\muμ is R(μ;X)=E[s(X)]R(\mu; X) = \mathbb{E}[s(X)]R(μ;X)=E[s(X)], and the optimal revenue is Rev(X)=sup⁡R(μ;X)\mathrm{Rev}(X) = \sup R(\mu; X)Rev(X)=supR(μ;X) over all IC and IR mechanisms. The bundled revenue is BRev(X)=Rev(Y+Z)\mathrm{BRev}(X) = \mathrm{Rev}(Y + Z)BRev(X)=Rev(Y+Z), the optimal revenue from selling the bundle as a single good.

A two-good mechanism is symmetric if q1(y,z)=q2(z,y)q_1(y,z) = q_2(z,y)q1​(y,z)=q2​(z,y) and s(y,z)=s(z,y)s(y,z) = s(z,y)s(y,z)=s(z,y). For a level a≥0a \ge 0a≥0, the bundled majorant μ^=(q^,s^)\hat\mu = (\hat q, \hat s)μ^​=(q^​,s^) of μ\muμ is

q^(y,z)=(q1(y+z−a,a), q1(y+z−a,a)),s^(y,z)=q^(y,z)⋅(y,z)−b(y+z−a,a).\hat q(y,z) = \bigl(q_1(y+z-a, a),\ q_1(y+z-a, a)\bigr), \qquad \hat s(y,z) = \hat q(y,z)\cdot(y,z) - b(y+z-a, a).q^​(y,z)=(q1​(y+z−a,a), q1​(y+z−a,a)),s^(y,z)=q^​(y,z)⋅(y,z)−b(y+z−a,a).

It depends on (y,z)(y,z)(y,z) only through y+zy + zy+z, so it is a bundled mechanism.

Formalization targets

Goal: Theorem 16

Let a>0a > 0a>0 and let FFF have a density fff supported in [a,∞)[a, \infty)[a,∞), differentiable on (a,∞)(a,\infty)(a,∞), with

xf′(x)+32f(x)≤0(x>a).(9)x f'(x) + \tfrac32 f(x) \le 0 \qquad (x > a). \tag{9}xf′(x)+23​f(x)≤0(x>a).(9)

Then for two i.i.d.-FFF goods,

Rev(X1,X2)=BRev(X1,X2)=Rev(X1+X2).\mathrm{Rev}(X_1, X_2) = \mathrm{BRev}(X_1, X_2) = \mathrm{Rev}(X_1 + X_2).Rev(X1​,X2​)=BRev(X1​,X2​)=Rev(X1​+X2​).

The substance is Rev≤BRev\mathrm{Rev} \le \mathrm{BRev}Rev≤BRev. The reverse inequality holds for every law, and the second equality is the definition of BRev\mathrm{BRev}BRev.

Milestones

  1. Proposition 5 (at least one good): IC holds if and only if bbb is convex and q(x)q(x)q(x) is a subgradient of bbb at every xxx.
  2. Proposition 6 (iii) (at least one good): the supremum defining Rev\mathrm{Rev}Rev may be restricted to NPT mechanisms.
  3. Symmetrization (p. 32): for i.i.d. goods, the symmetrization μˉ\bar\muμˉ​ is IC, IR, NPT and symmetric, and it has the same revenue.
  4. Display (18): for symmetric IC μ\muμ and y,z≥ay, z \ge ay,z≥a, b(y,z)≤b(y+z−a,a)b(y,z) \le b(y+z-a, a)b(y,z)≤b(y+z−a,a).
  5. μ^\hat\muμ^​ is feasible, IC and IR on the quadrant [a,∞)2[a,\infty)^2[a,∞)2, and bundled.
  6. Under (9), R(μ;X)≤R(μ^;X)R(\mu; X) \le R(\hat\mu; X)R(μ;X)≤R(μ^​;X) for symmetric IC, IR and NPT μ\muμ.

Companion: Corollary 17

For two i.i.d. goods with the equal-revenue law P[V≥p]=1/p\mathbb{P}[V \ge p] = 1/pP[V≥p]=1/p (p≥1p \ge 1p≥1),

Rev(V1,V2)=BRev(V1,V2)=2(w+1)≈2.56,\mathrm{Rev}(V_1,V_2) = \mathrm{BRev}(V_1,V_2) = 2(w+1) \approx 2.56,Rev(V1​,V2​)=BRev(V1​,V2​)=2(w+1)≈2.56,

where wew+1=1w e^{w+1} = 1wew+1=1.

Significance

Theorem 16 is one of the few closed-form optimality results for multi-good monopoly. It gives an explicit, checkable condition: x3/2f(x)x^{3/2} f(x)x3/2f(x) is nonincreasing. The condition covers the equal-revenue distribution and Pareto distributions with index α≥1/2\alpha \ge 1/2α≥1/2. On these distributions the optimal revenue reduces to a one-dimensional pricing problem for Y+ZY + ZY+Z. Corollary 17 is the resulting exact value for two equal-revenue goods. Elsewhere in the paper, Theorems C and D use the equal-revenue law as the extremal example for the approximation ratios.

The theorem is proved in the paper. As far as is known it has not been machine-checked. The mission formalizes the model and every step of the proof as a separate milestone. These steps are Rochet's characterization (Proposition 5), the NPT reduction, symmetrization, the convexity majorization (18), and the revenue comparison under (9). Propositions 5 and 6 (iii) and symmetrization are shared with other missions of this series, so they are reusable beyond Theorem 16.

Difficulty

The obvious approach is to show that the optimal mechanism is bundled by exhibiting it. This fails because nothing is known about the optimal mechanism, and it need not exist. The proof instead compares an arbitrary symmetric mechanism with its bundled majorant μ^\hat\muμ^​. A pointwise comparison fails: b^≥b\hat b \ge bb^≥b does not imply s^≥s\hat s \ge ss^≥s. The comparison holds only in expectation and only for densities satisfying (9); without (9) bundling is not optimal in general.

Making this rigorous requires several things:

  • the a.e. differentiability of the convex function bbb, together with q1=∂b/∂yq_1 = \partial b/\partial yq1​=∂b/∂y a.e.;
  • integration by parts on truncated squares [a,M]2[a, M]^2[a,M]2 with possibly infinite total revenue;
  • control of the boundary term at y=ay = ay=a, where fff need not be bounded or continuous;
  • the extension of μ^\hat\muμ^​, defined on [a,∞)2[a,\infty)^2[a,∞)2, to an IC and IR mechanism on all of R+2\mathbb{R}^2_+R+2​.

Formalization scope

Valuations are functions Fin 2→R≥0\mathrm{Fin}\ 2 \to \mathbb{R}_{\ge 0}Fin 2→R≥0​, and random valuations are given by their laws. "Two i.i.d. goods" is the product measure ν⊗ν\nu \otimes \nuν⊗ν. The law ν\nuν is the image on R≥0\mathbb{R}_{\ge0}R≥0​ of f(x) dxf(x)\,dxf(x)dx, so no mass is lost because f=0f = 0f=0 below a>0a > 0a>0. "Density function" means: fff measurable, f≥0f \ge 0f≥0, ∫f=1\int f = 1∫f=1, and f=0f = 0f=0 on (−∞,a)(-\infty, a)(−∞,a). Continuity of FFF is then automatic. Differentiability of fff is required at every x>ax > ax>a, the range of (9). Nothing is assumed about fff at aaa; in particular fff may blow up there.

Admissible mechanisms are feasible, IC, IR and have a measurable payment function. The paper allows measurability without loss of generality (footnote 12). Revenue R(μ;X)R(\mu;X)R(μ;X) is an extended real, so a non-integrable payment never collapses to 000. Rev\mathrm{Rev}Rev takes values in [0,∞][0,\infty][0,∞] and may be infinite under (9), for instance for f(x)=c x−3/2f(x) = c\,x^{-3/2}f(x)=cx−3/2 on [a,∞)[a,\infty)[a,∞). Milestone 5 states IC and IR of μ^\hat\muμ^​ on the quadrant only, as the paper defines μ^\hat\muμ^​ there. Off the quadrant, y+z−ay + z - ay+z−a is truncated at 000.

The goal must not be stated as the second equality alone, which is definitional. Nor may it assume the law is bounded or that f(a)f(a)f(a) is finite: these would replace the theorem with a special case. The display (17) and the limit formula (16) are not posed, because their boundary term 2af(a)∫b(a,z)f(z) dz2a f(a)\int b(a,z) f(z)\,dz2af(a)∫b(a,z)f(z)dz refers to the value of fff at aaa, which the hypotheses do not control.

A complete development needs Rochet's characterization, a.e. gradients of convex functions on R+2\mathbb{R}^2_+R+2​, Fubini–Tonelli on the product law, one-dimensional integration by parts with an absolutely continuous weight, and Myerson's one-good formula for Corollary 17. Contributions of any of these as standalone lemmas are welcome.

Selected references

  • S. Hart and N. Nisan, Approximate Revenue Maximization with Multiple Items, J. Econ. Theory 172 (2017); arXiv:1204.1846v3. https://arxiv.org/abs/1204.1846
  • R. Myerson, Optimal Auction Design, Math. Oper. Res. 6 (1981) 58–73. https://doi.org/10.1287/moor.6.1.58
  • J.-C. Rochet, The taxation principle and multi-time Hamilton–Jacobi equations, J. Math. Econ. 14 (1985) 113–128. https://doi.org/10.1016/0304-4068(85)90015-1
  • A. Manelli and D. Vincent, Bundling as an optimal selling mechanism for a multiple-good monopolist, J. Econ. Theory 127 (2006) 1–35. https://doi.org/10.1016/j.jet.2005.08.007
  • S. Hart and P. Reny, Maximal revenue with multiple goods: nonmonotonicity and other observations, Theoretical Economics 10 (2015) 893–922. https://doi.org/10.3982/TE1517
  • C. Daskalakis, A. Deckelbaum and C. Tzamos, Strong duality for a multiple-good monopolist, Econometrica 85 (2017) 735–767. https://doi.org/10.3982/ECTA12618
9 thms1 active userReviewed
Harmonic AnalysisMachine Learning·Captain: mikedeng1

Approximation Bounds for Random Neural Networks and Reservoir Systems I: Random ReLU Networks with Sampled Inner Weights Approximate Fourier-Represented Functions on a Hilbert Space (Theorem 1)Research Paper

Motivation

In random-feature models, also called extreme learning machines or random neural networks, the inner weights of a one-hidden-layer network are drawn at random and only the linear readout is trained. Training is then a least-squares problem, which is why such models, and their dynamic counterparts reservoir computing and echo state networks, are used in practice. The question is how well such a network can approximate a target when it is not allowed to choose its inner weights. Gonon, Grigoryeva and Ortega (Ann. Appl. Probab. 33 (2023) 28–69) answer it with explicit, non-asymptotic L2L^2L2 bounds. Their first result, Theorem 1, covers static ReLU networks on an arbitrary separable Hilbert space. The later results of the paper, on finite-dimensional inputs and on echo state networks, are built on it.

Timeline:

  • 1993: Barron (IEEE Trans. Inf. Theory 39) proves O(N−1/2)O(N^{-1/2})O(N−1/2) approximation rates for functions with a finite first Fourier moment, using networks whose inner weights are optimized.
  • 2008: Rahimi and Recht (NeurIPS 2008) prove random kitchen-sink bounds for random features sampled from a fixed law.
  • 2023: Gonon, Grigoryeva and Ortega prove bounds with explicit constants for ReLU units with importance-sampled inner weights, on infinite-dimensional input spaces, under weak conditions on the sampling law.

Setting

X\mathcal XX is a separable real Hilbert space with inner product ⟨⋅,⋅⟩\langle\cdot,\cdot\rangle⟨⋅,⋅⟩, and σ(x)=max⁡(x,0)\sigma(x)=\max(x,0)σ(x)=max(x,0) is the rectifier. The inner weights θi=(Ai,ζi)∈X×R\theta_i=(A_i,\zeta_i)\in\mathcal X\times\mathbb Rθi​=(Ai​,ζi​)∈X×R, i=1,…,Ni=1,\dots,Ni=1,…,N, are i.i.d. with law π\piπ. For a readout matrix W∈Mm,NW\in\mathbb M_{m,N}W∈Mm,N​ the random ReLU network is

HWA,ζ(z)=Wσ(Az+ζ),(Az+ζ)i=⟨Ai,z⟩+ζi.H_W^{A,\zeta}(z)=W\sigma(Az+\zeta),\qquad (Az+\zeta)_i=\langle A_i,z\rangle+\zeta_i .HWA,ζ​(z)=Wσ(Az+ζ),(Az+ζ)i​=⟨Ai​,z⟩+ζi​.

The target H∗:X→RmH^*:\mathcal X\to\mathbb R^mH∗:X→Rm is assumed to have a Fourier representation on the ball of radius MMM: Hj∗(z)=∫ei⟨w,z⟩μ^j(dw)H^*_j(z)=\int e^{i\langle w,z\rangle}\hat\mu_j(dw)Hj∗​(z)=∫ei⟨w,z⟩μ^​j​(dw) for ∥z∥≤M\|z\|\le M∥z∥≤M, with complex measures μ^j\hat\mu_jμ^​j​ whose total variations ∣μ^j∣|\hat\mu_j|∣μ^​j​∣ have finite second moments (condition (15)). The sampling law has the product form π=πX⊗(πR(x)dx)\pi=\pi_{\mathcal X}\otimes(\pi_{\mathbb R}(x)dx)π=πX​⊗(πR​(x)dx). Two further quantities enter: Fπ(x)=2∫−x0πR(u)−1duF_\pi(x)=2\int_{-x}^0\pi_{\mathbb R}(u)^{-1}duFπ​(x)=2∫−x0​πR​(u)−1du, and the density gj=d(∣μ^j∣+∣μ^j∣−)/dπXg_j=d(|\hat\mu_j|+|\hat\mu_j|^-)/d\pi_{\mathcal X}gj​=d(∣μ^​j​∣+∣μ^​j​∣−)/dπX​, where ν−(⋅)=ν(−⋅)\nu^-(\cdot)=\nu(-\cdot)ν−(⋅)=ν(−⋅). The input ZZZ is an X\mathcal XX-valued random variable with law μZ\mu_ZμZ​ and ∥Z∥≤M\|Z\|\le M∥Z∥≤M, independent of the inner weights.

Formalization targets

Goal: Theorem 1 (pp. 37–38)

If ∣μ^j∣+∣μ^j∣−≪πX|\hat\mu_j|+|\hat\mu_j|^-\ll\pi_{\mathcal X}∣μ^​j​∣+∣μ^​j​∣−≪πX​, condition (i) or (ii) on πR\pi_{\mathbb R}πR​ holds, and the integrals (16) of Fπ(M∥w∥)∥w∥2gj2F_\pi(M\|w\|)\|w\|^2g_j^2Fπ​(M∥w∥)∥w∥2gj2​ and max⁡(∥w∥2,1)gj2\max(\|w\|^2,1)g_j^2max(∥w∥2,1)gj2​ are finite, then there is a readout WWW, a measurable function of (A,ζ)(A,\zeta)(A,ζ), with

E[∥HWA,ζ(Z)−H∗(Z)∥2]≤C∗N,\mathbb E\big[\|H_W^{A,\zeta}(Z)-H^*(Z)\|^2\big]\le\frac{C^*}{N},E[∥HWA,ζ​(Z)−H∗(Z)∥2]≤NC∗​,

and for every δ∈(0,1)\delta\in(0,1)δ∈(0,1), with probability at least 1−δ1-\delta1−δ, ∥HWA,ζ−H∗∥L2(μZ)≤C∗/(δN)\|H_W^{A,\zeta}-H^*\|_{L^2(\mu_Z)}\le\sqrt{C^*}/(\delta\sqrt N)∥HWA,ζ​−H∗∥L2(μZ​)​≤C∗​/(δN​). Here C∗=∑jCj∗C^*=\sum_jC^*_jC∗=∑j​Cj∗​ and

Cj∗=M2 ⁣∫Fπ(M∥w∥)∥w∥2gj2 dπX+32max⁡(M2,1)(Fπ(1)−Fπ(−1)) ⁣∫max⁡(∥w∥2,1)gj2 dπX.C^*_j=M^2\!\int F_\pi(M\|w\|)\|w\|^2g_j^2\,d\pi_{\mathcal X}+32\max(M^2,1)\big(F_\pi(1)-F_\pi(-1)\big)\!\int\max(\|w\|^2,1)g_j^2\,d\pi_{\mathcal X}.Cj∗​=M2∫Fπ​(M∥w∥)∥w∥2gj2​dπX​+32max(M2,1)(Fπ​(1)−Fπ​(−1))∫max(∥w∥2,1)gj2​dπX​.

Milestones

  1. Lemma 1 (p. 37): (ω,z)↦HW(ω)A(ω),ζ(ω)(z)(\omega,z)\mapsto H_{W(\omega)}^{A(\omega),\zeta(\omega)}(z)(ω,z)↦HW(ω)A(ω),ζ(ω)​(z) is product-measurable.
  2. Step 1 (p. 39): ∫i⟨v,w⟩μ^(dw)\int i\langle v,w\rangle\hat\mu(dw)∫i⟨v,w⟩μ^​(dw) is real, and the identity (19).
  3. (20) (p. 40): H∗(v)=∫σ(⟨v,w⟩+u) α(dw,du)H^*(v)=\int\sigma(\langle v,w\rangle+u)\,\alpha(dw,du)H∗(v)=∫σ(⟨v,w⟩+u)α(dw,du) for ∥v∥≤M\|v\|\le M∥v∥≤M, for an explicit signed measure α\alphaα.
  4. (21) (p. 40): α≪π\alpha\ll\piα≪π, with ∣dα/dπ(w,u)∣≤(1(−M∥w∥,0](u)+22 1[−1,1](u)) g(w)/πR(u)|d\alpha/d\pi(w,u)|\le(\mathbb 1_{(-M\|w\|,0]}(u)+2\sqrt2\,\mathbb 1_{[-1,1]}(u))\,g(w)/\pi_{\mathbb R}(u)∣dα/dπ(w,u)∣≤(1(−M∥w∥,0]​(u)+22​1[−1,1]​(u))g(w)/πR​(u).
  5. (25) (p. 40): the mean-square error of the i.i.d. average is 1NVar⁡\frac1N\operatorname{Var}N1​Var, which is at most 1N\frac1NN1​ times the second moment.
  6. (26) (p. 41): that second moment is at most C∗C^*C∗.

Milestones 2–6 are stated for a single output component (m=1m=1m=1), to which the proof reduces.

Significance

Theorem 1 says that only the readout has to be trained. With inner weights drawn from an essentially arbitrary product law satisfying (i) or (ii), a least-squares readout achieves L2L^2L2 error O(N−1/2)O(N^{-1/2})O(N−1/2) with a constant computable from the Fourier data of H∗H^*H∗ and the sampling law (Remark 3, p. 38). The finite-dimensional rates of Section 4.2 (Proposition 3, Corollary 1) follow from it, and so do, through the reservoir functional's Fourier representation, the echo state network bounds of Sections 5–6.

The result is proved on paper but has no machine-checked proof. A formal proof needs integral representations of Fourier-type functions by ReLU units, Radon–Nikodym arguments for signed measures on product spaces, and the variance identity for i.i.d. averages. A side product of formalizing it is a correction: the printed constant 32M232M^232M2 is false for small MMM and has to be replaced by 32max⁡(M2,1)32\max(M^2,1)32max(M2,1), which agrees with it for M≥1M\ge1M≥1.

Difficulty

There is an obvious argument: sample features from ∣μ^∣|\hat\mu|∣μ^​∣ and use Maurey's empirical method. It does not apply, because the features here are ReLU ridge functions of a real-valued bias drawn from πR\pi_{\mathbb R}πR​, not complex exponentials. The proof must first rewrite the oscillatory representation as an integral of rectifiers, (18)–(20). That step uses the realness of ∫i⟨v,w⟩μ^(dw)\int i\langle v,w\rangle\hat\mu(dw)∫i⟨v,w⟩μ^​(dw), which is not given and follows only from H∗H^*H∗ being real near 000. The signed measure α\alphaα obtained this way must then be absolutely continuous with respect to the sampling law π\piπ, with a density whose second moment is controlled by (16). In case (ii), πR\pi_{\mathbb R}πR​ may vanish outside a bounded interval, and the support of α\alphaα has to be matched to it.

Formalization scope

  • Space. X\mathcal XX is a complete, separable real inner product space with its Borel σ\sigmaσ-algebra.
  • Complex measures. Each μ^j\hat\mu_jμ^​j​ is given in polar form hj νjh_j\,\nu_jhj​νj​, with νj=∣μ^j∣\nu_j=|\hat\mu_j|νj​=∣μ^​j​∣ a finite measure and hjh_jhj​ measurable with ∣hj∣=1|h_j|=1∣hj​∣=1. This loses no generality (Rudin, Thm 6.12).
  • Randomness. The inner weights live on (X×R)N(\mathcal X\times\mathbb R)^N(X×R)N with law π⊗N\pi^{\otimes N}π⊗N, and ZZZ has its own law μZ\mu_ZμZ​. Expectations are iterated lower integrals; this is the standing independence of §4.1. "With probability 1−δ1-\delta1−δ" bounds the π⊗N\pi^{\otimes N}π⊗N-measure of the bad set.
  • Readout. The readout is a measurable function of the inner weights only. A readout allowed to depend on the input ZZZ would make the statement trivial, and it is ruled out.
  • Norms and normalisation. Norms on Rm\mathbb R^mRm are Euclidean, as a sum of squares. πX\pi_{\mathcal X}πX​ is a probability measure and ∫πR=1\int\pi_{\mathbb R}=1∫πR​=1; rescaling πX\pi_{\mathcal X}πX​ by ccc and πR\pi_{\mathbb R}πR​ by 1/c1/c1/c changes neither π\piπ nor C∗C^*C∗.
  • Encodings. "∣Fπ(x)∣<∞|F_\pi(x)|<\infty∣Fπ​(x)∣<∞" is integrability of 1/πR1/\pi_{\mathbb R}1/πR​ on the interval between −x-x−x and 000, and FπF_\piFπ​ is an oriented interval integral. (15) is a finite lower integral and (16) is integrability, so none of them holds vacuously.
  • Departures from the print. The page's "C∗>0C^*>0C∗>0" is dropped in favour of the explicit value. The constant uses 32max⁡(M2,1)32\max(M^2,1)32max(M2,1).

A complete development needs the following:

  • the rectifier integral identity (18), which is identity (8) of Proposition 2 and is posed in the companion mission on echo state networks;
  • the Radon–Nikodym and Fubini arguments of Step 1;
  • the variance of i.i.d. averages under product measures;
  • Markov's inequality.

The variance identity (25) and the bound (26) are independent of the Fourier setting and reusable. Contributions to any milestone are welcome, and so are alternative proofs of (19) and (20).

Selected references

  • L. Gonon, L. Grigoryeva, J.-P. Ortega, Approximation bounds for random neural networks and reservoir systems, Ann. Appl. Probab. 33(1), 2023, 28–69. https://doi.org/10.1214/22-AAP1806
  • A. R. Barron, Universal approximation bounds for superpositions of a sigmoidal function, IEEE Trans. Inf. Theory 39(3), 1993, 930–945. https://doi.org/10.1109/18.256500
  • A. Rahimi, B. Recht, Weighted sums of random kitchen sinks: replacing minimization with randomization in learning, NeurIPS 2008. https://papers.nips.cc/paper/2008/hash/0efe32849d230d7f53049ddc4a4b0c60-Abstract.html
  • W. Rudin, Real and Complex Analysis, 3rd ed., McGraw-Hill, 1987 (Chapter 6, Theorem 6.12).
8 thms1 active userReviewed
Algorithmic Game TheoryMechanism DesignOperations Research·Captain: mikedeng1

Approximate Revenue Maximization with Multiple Items 1: For Independent Groups of Goods Y and Z, Rev(Y, Z) ≤ Rev(Y) + Rev(Z) + BRev(Y) + BRev(Z)Research Paper

Motivation

A seller who offers several goods to a buyer with private values has to choose a selling mechanism. For one good the optimal mechanism is known: post the best take-it-or-leave-it price (Myerson 1981). For two or more goods the optimal mechanism can be complicated: it may use lotteries, its menu may be infinite, and it may not be computable in closed form even for simple value distributions. Sellers therefore use simple mechanisms, such as selling each good separately or selling all goods as one bundle, and the question is how much revenue this costs.

Hart and Nisan (arXiv:1204.1846v3, J. Econ. Theory 2017) answer this with guaranteed fractions of the optimal revenue. Their basic tool is a decomposition theorem: when the goods split into two groups whose values are independent, the optimal revenue from all goods together is bounded by revenues computed from each group separately. This mission formalizes that theorem (Theorem 7) and the paper's first headline result, Theorem A, which is its two-good case.

Setting

There are k≥1k \ge 1k≥1 goods and one buyer with additive values. A valuation is a vector x∈R+kx \in \mathbb{R}^k_+x∈R+k​, and a kkk-good random valuation XXX is a random vector in R+k\mathbb{R}^k_+R+k​ whose law the seller knows.

A mechanism μ=(q,s)\mu = (q, s)μ=(q,s) consists of an allocation q:R+k→[0,1]kq : \mathbb{R}^k_+ \to [0,1]^kq:R+k​→[0,1]k, where qi(x)q_i(x)qi​(x) is the probability that the buyer receives good iii when reporting xxx, and a payment s:R+k→Rs : \mathbb{R}^k_+ \to \mathbb{R}s:R+k​→R. The buyer's payoff from truthful reporting is b(x)=q(x)⋅x−s(x)b(x) = q(x)\cdot x - s(x)b(x)=q(x)⋅x−s(x). The mechanism is

  • incentive compatible (IC) if b(x)≥q(x~)⋅x−s(x~)b(x) \ge q(\tilde x)\cdot x - s(\tilde x)b(x)≥q(x~)⋅x−s(x~) for all x,x~∈R+kx, \tilde x \in \mathbb{R}^k_+x,x~∈R+k​;
  • individually rational (IR) if b(x)≥0b(x) \ge 0b(x)≥0 for all xxx;
  • no positive transfer (NPT) if s(x)≥0s(x) \ge 0s(x)≥0 for all xxx.

Let M\mathcal MM be the class of IC and IR mechanisms. The expected revenue of μ∈M\mu \in \mathcal Mμ∈M is R(μ;X)=E[s(X)]R(\mu; X) = \mathbb{E}[s(X)]R(μ;X)=E[s(X)], which lies in (−∞,+∞](-\infty, +\infty](−∞,+∞], and the optimal revenue is

Rev(X)=sup⁡μ∈MR(μ;X)∈[0,∞].\mathrm{Rev}(X) = \sup_{\mu \in \mathcal M} R(\mu; X) \in [0, \infty].Rev(X)=μ∈Msup​R(μ;X)∈[0,∞].

Two benchmarks reduce to one-good problems: the separate revenue SRev(X)=Rev(X1)+⋯+Rev(Xk)\mathrm{SRev}(X) = \mathrm{Rev}(X_1) + \dots + \mathrm{Rev}(X_k)SRev(X)=Rev(X1​)+⋯+Rev(Xk​) and the bundled revenue BRev(X)=Rev(X1+⋯+Xk)\mathrm{BRev}(X) = \mathrm{Rev}(X_1 + \dots + X_k)BRev(X)=Rev(X1​+⋯+Xk​). Finally Val(X)=E[∑iXi]\mathrm{Val}(X) = \mathbb{E}[\sum_i X_i]Val(X)=E[∑i​Xi​] is the expected total value.

For the decomposition, YYY is a random valuation for k1≥1k_1 \ge 1k1​≥1 goods and ZZZ one for k2≥1k_2 \ge 1k2​≥1 goods. The coordinates of YYY may be arbitrarily dependent, and so may those of ZZZ, but the vectors YYY and ZZZ are independent. Rev(Y,Z)\mathrm{Rev}(Y, Z)Rev(Y,Z) is the optimal revenue from selling all k1+k2k_1 + k_2k1​+k2​ goods together, and X1X∈AX\mathbf 1_{X \in A}X1X∈A​ is the random valuation equal to XXX on the event X∈AX \in AX∈A and to 000 otherwise.

Formalization targets

Goal: Theorem 7 (p. 18)

If YYY and ZZZ are independent, then

Rev(Y,Z)≤Rev(Y)+Rev(Z)+BRev(Y)+BRev(Z)≤2(Rev(Y)+Rev(Z)).\mathrm{Rev}(Y, Z) \le \mathrm{Rev}(Y) + \mathrm{Rev}(Z) + \mathrm{BRev}(Y) + \mathrm{BRev}(Z) \le 2\big(\mathrm{Rev}(Y) + \mathrm{Rev}(Z)\big).Rev(Y,Z)≤Rev(Y)+Rev(Z)+BRev(Y)+BRev(Z)≤2(Rev(Y)+Rev(Z)).

Both inequalities are part of the goal, with the same middle term.

Milestones, in the order of the paper's proof

  1. Proposition 6 (iii), p. 14: Rev(X)\mathrm{Rev}(X)Rev(X) is the supremum of R(μ;X)R(\mu; X)R(μ;X) over IC, IR and NPT mechanisms.
  2. The NPT split, p. 20: Rev(Y,Z)≤Rev((Y,Z)1∑iYi≥∑jZj)+Rev((Y,Z)1∑jZj≥∑iYi)\mathrm{Rev}(Y,Z) \le \mathrm{Rev}\big((Y,Z)\mathbf 1_{\sum_i Y_i \ge \sum_j Z_j}\big) + \mathrm{Rev}\big((Y,Z)\mathbf 1_{\sum_j Z_j \ge \sum_i Y_i}\big)Rev(Y,Z)≤Rev((Y,Z)1∑i​Yi​≥∑j​Zj​​)+Rev((Y,Z)1∑j​Zj​≥∑i​Yi​​).
  3. Proposition 6 (iv), p. 14: for μ∈M\mu \in \mathcal Mμ∈M and a set AAA of values, E[s(X)1X∈A]≤Rev(X1X∈A)≤Rev(X)\mathbb{E}[s(X)\mathbf 1_{X\in A}] \le \mathrm{Rev}(X\mathbf 1_{X\in A}) \le \mathrm{Rev}(X)E[s(X)1X∈A​]≤Rev(X1X∈A​)≤Rev(X).
  4. Lemma 8, p. 19: Rev((Y,Z)1(Y,Z)∈A)≤Rev(Y)+Val(Z1(Y,Z)∈A)\mathrm{Rev}\big((Y,Z)\mathbf 1_{(Y,Z)\in A}\big) \le \mathrm{Rev}(Y) + \mathrm{Val}\big(Z\mathbf 1_{(Y,Z)\in A}\big)Rev((Y,Z)1(Y,Z)∈A​)≤Rev(Y)+Val(Z1(Y,Z)∈A​).
  5. Lemma 9, p. 20: for one-dimensional independent Y,ZY, ZY,Z, Val(Z1Y≥Z)≤Rev(Y)\mathrm{Val}(Z\mathbf 1_{Y \ge Z}) \le \mathrm{Rev}(Y)Val(Z1Y≥Z​)≤Rev(Y).
  6. Lemma 10, p. 20: Val(Z1∑iYi≥∑jZj)≤BRev(Y)\mathrm{Val}\big(Z\mathbf 1_{\sum_i Y_i \ge \sum_j Z_j}\big) \le \mathrm{BRev}(Y)Val(Z1∑i​Yi​≥∑j​Zj​​)≤BRev(Y).
  7. BRev ≤ Rev, p. 18.

Companions

Theorem A in the form (2), p. 16: for two independent goods Y,ZY, ZY,Z, Rev(Y,Z)≤2(Rev(Y)+Rev(Z))\mathrm{Rev}(Y, Z) \le 2(\mathrm{Rev}(Y) + \mathrm{Rev}(Z))Rev(Y,Z)≤2(Rev(Y)+Rev(Z)). Also posed is display (5), p. 17: z⋅P[Y≥z]≤Rev(Y)z\cdot\mathbb{P}[Y \ge z] \le \mathrm{Rev}(Y)z⋅P[Y≥z]≤Rev(Y) for one good and every z≥0z \ge 0z≥0.

Significance

Theorem A says that for two independent goods, selling them separately guarantees half of the optimal revenue, whatever the distributions. Theorem 7 extends this to two independent groups of goods. It is the induction step behind the paper's bounds for many goods: separate selling guarantees a fraction c/log⁡2kc/\log^2 kc/log2k of the optimal revenue for kkk independent goods (Theorem C), and bundling guarantees c/log⁡kc/\log kc/logk for kkk i.i.d. goods (Theorem D). These are other missions of this series, and both use Theorem 7 as a black box. The same approach underlies later work on simple versus optimal mechanisms, for example Babaioff, Immorlica, Lucier and Weinberg (FOCS 2014) and Li and Yao (PNAS 2013).

All results posed here are proved in the paper. As far as we know, none has a machine-checked proof. A formal proof would give a reusable development of mechanism design with multi-dimensional types: IC and IR mechanisms on R+k\mathbb{R}^k_+R+k​, revenue as an extended-real expectation, restriction of mechanisms to subdomains, and marginal mechanisms obtained by fixing some coordinates of an independent product.

Difficulty

The proofs are short on paper but rest on measure theory the prose passes over. Lemma 8 conditions on Z=zZ = zZ=z and applies Proposition 6 (iv) to a one-group mechanism for each fixed zzz. Formally this is a Tonelli argument on a product measure, with an integrand that may be negative on part of the space and +∞+\infty+∞ in expectation. The induced mechanism sz(y)=s(y,z)−∑jqj(y,z)zjs^z(y) = s(y, z) - \sum_j q_j(y,z) z_jsz(y)=s(y,z)−∑j​qj​(y,z)zj​ involves the allocation qqq, which is not assumed measurable. Proposition 6 (iv) bounds a mechanism with possibly negative payments by the optimal revenue on a modified valuation X1X∈AX\mathbf 1_{X \in A}X1X∈A​, and this requires shifting the mechanism to make it NPT. Finally, every revenue may be infinite, so the inequalities must be handled in [0,∞][0,\infty][0,∞] and in the extended reals, without cancelling.

A natural first idea is to apply Theorem A coordinate by coordinate. This fails, because the coordinates inside YYY may be dependent. Only the independence of YYY from ZZZ is available, and the bound must go through the bundled revenue BRev\mathrm{BRev}BRev.

Formalization scope

  • Valuations as laws. A random valuation is represented by its law, a probability measure on ι → ℝ≥0 for a finite type ι of goods. Every quantity of the paper depends only on the law. One good is ι = Unit.
  • Independence as a product law. The law of (Y,Z)(Y,Z)(Y,Z) is the product of the laws of YYY and ZZZ, transported to ι₁ ⊕ ι₂ → ℝ≥0. Theorems that need independence are never stated over an arbitrary joint law with given marginals, since they are false for correlated groups. The NPT split in the proof of Theorem 7 is stated for independent groups, as on the page. Proposition 6 applies to an arbitrary law.
  • The class M\mathcal MM. M\mathcal MM requires q∈[0,1]kq \in [0,1]^kq∈[0,1]k, IC, IR and a measurable payment sss, the last as allowed by the paper's footnote 12. qqq is not required to be measurable.
  • Revenue. R(μ;X)R(\mu;X)R(μ;X) is ∫s+−∫s−\int s^+ - \int s^-∫s+−∫s− in the extended reals, so that E[s(X)]=+∞\mathbb{E}[s(X)] = +\inftyE[s(X)]=+∞ is kept. Since s≥s(0)s \ge s(0)s≥s(0) under IC, the negative part is finite. The Bochner integral is not used, because it would assign 000 to a non-integrable payment. Rev\mathrm{Rev}Rev, SRev\mathrm{SRev}SRev, BRev\mathrm{BRev}BRev and Val\mathrm{Val}Val take values in [0,∞][0,\infty][0,∞].
  • Subdomains. X1X∈AX\mathbf 1_{X\in A}X1X∈A​ is the image law under x↦x1A(x)x \mapsto x\mathbf 1_A(x)x↦x1A​(x), and AAA is assumed measurable. The regions of the proof of Theorem 7 are both closed, so the diagonal is counted twice, as in the paper.
  • Not trivializable. A formalization in which Rev\mathrm{Rev}Rev is computed with a Bochner integral, the expected revenue is clipped before comparing a single mechanism with Rev\mathrm{Rev}Rev, or independence is dropped does not state the paper's theorem, and none of these is used.
  • Not posed. The ratio form of Theorem A (GFOR≥1/2\mathrm{GFOR} \ge 1/2GFOR≥1/2) is not posed, because of its conventions 0/0=10/0 = 10/0=1 and ∞/∞\infty/\infty∞/∞. Theorems B, C, D, 16 and 33 are separate missions of this series.

Contributions are welcome at every level: Proposition 6 (the NPT shift and subdomain restriction), the marginal-mechanism construction behind Lemma 8, and the one-good facts (5) and Lemma 9 are all reusable in the sibling missions on Theorems C and D.

Selected references

  • S. Hart and N. Nisan, Approximate Revenue Maximization with Multiple Items, Journal of Economic Theory 172 (2017); preprint arXiv:1204.1846v3, 2017. https://arxiv.org/abs/1204.1846
  • R. B. Myerson, Optimal Auction Design, Mathematics of Operations Research 6 (1981) 58–73. https://doi.org/10.1287/moor.6.1.58
  • S. Hart and P. J. Reny, Maximal Revenue with Multiple Goods: Nonmonotonicity and Other Observations, Theoretical Economics 10 (2015) 893–922. https://doi.org/10.3982/TE1517
  • M. Babaioff, N. Immorlica, B. Lucier and S. M. Weinberg, A Simple and Approximately Optimal Mechanism for an Additive Buyer, FOCS 2014. https://arxiv.org/abs/1405.6146
  • X. Li and A. C.-C. Yao, On Revenue Maximization for Selling Multiple Independently Distributed Items, PNAS 110 (2013) 11232–11237. https://doi.org/10.1073/pnas.1309533110
9 thms1 active userReviewed
CombinatoricsOperations ResearchOptimization·Captain: mikedeng1

Online Contention Resolution Schemes with Applications to Bayesian Selection Problems 2: The Degree Relaxation of Matching Has a (b, e^(−2b))-Selectable Randomized Greedy OCRSResearch Paper

Why online matching needs a stronger rounding guarantee

In online selection, candidate items reveal themselves one at a time and an algorithm must decide immediately whether to keep each one. When the kept items must form a matching, accepting one edge can prevent a later edge sharing an endpoint from being accepted. A fractional solution supplies useful marginal probabilities, but it does not by itself specify which active edges can be kept under an unfavorable arrival order. Online contention resolution schemes provide that missing rounding step: they turn independent edge activity into a feasible online selection while guaranteeing each edge a chance of selection. Feldman, Svensson, and Zenklusen introduced the strong selectability guarantee used here, which protects an edge even against an arrival order chosen with knowledge of all random outcomes (Feldman–Svensson–Zenklusen, §1.1).

This mission isolates the matching result of their paper. It concerns a degree relaxation that is weaker than the full matching polytope, so the result applies even when the input fractional point satisfies only local vertex constraints. The source states the existence theorem as Theorem 2.7 and gives its numerical guarantee explicitly (Theorem 2.7, p. 13).

Finite graphs, fractional edges, and greedy families

Let G=(V,E)G=(V,E)G=(V,E) be a finite graph. An edge ggg has two distinct endpoints; multiple edges may have the same endpoints. A matching is a set of edges with no shared vertex. For a vertex uuu, write δ(u)\delta(u)δ(u) for its incident edges. The degree relaxation is

PG={x∈R≥0E:∑g∈δ(u)xg≤1 for every u∈V}.P_G=\left\{x\in\mathbb R_{\ge0}^{E}:\sum_{g\in\delta(u)}x_g\le1\text{ for every }u\in V\right\}.PG​=⎩⎨⎧​x∈R≥0E​:g∈δ(u)∑​xg​≤1 for every u∈V⎭⎬⎫​.

Since each edge meets a vertex, these constraints also give xg≤1x_g\le1xg​≤1. Every matching has a characteristic vector in PGP_GPG​. The degree relaxation can contain points outside the convex hull of matchings, which makes the theorem stronger than a statement restricted to the matching polytope (§2.2, p. 13).

For x∈[0,1]Ex\in[0,1]^Ex∈[0,1]E, the active set R(x)R(x)R(x) includes each edge ggg independently with probability xgx_gxg​. A greedy family Fx\mathcal F_xFx​ is a downward closed family of matchings containing the empty set. A greedy online scheme accepts an active edge if adding it to the already accepted set remains in Fx\mathcal F_xFx​. The family may itself be random, sampled independently of R(x)R(x)R(x).

An edge ggg is selectable for a realized active set and greedy family when every member I∈FxI\in\mathcal F_xI∈Fx​ with I⊆R(x)I\subseteq R(x)I⊆R(x) can be extended by ggg while staying in Fx\mathcal F_xFx​. Thus selectability promises acceptance regardless of which eligible active edges preceded ggg. A scheme is (b,c)(b,c)(b,c)-selectable if this event has probability at least ccc for every x∈bPGx\in bP_Gx∈bPG​ and every edge. The probability includes both the active-set draw and any random family draw (Definitions 1.3–1.5, pp. 3–4).

Formalization targets

Theorem 2.7: a randomized matching scheme

For every b∈[0,1]b\in[0,1]b∈[0,1] and every finite loopless graph GGG,

∃ a randomized greedy scheme for PG that is (b,e−2b)-selectable.\exists\text{ a randomized greedy scheme for }P_G\text{ that is }(b,e^{-2b})\text{-selectable}.∃ a randomized greedy scheme for PG​ that is (b,e−2b)-selectable.

Equivalently, for each x∈bPGx\in bP_Gx∈bPG​ and each edge g∈Eg\in Eg∈E, the probability that ggg is selectable is at least e−2be^{-2b}e−2b. The statement includes edges with xg=0x_g=0xg​=0; the selectability event concerns the ability to add an edge, rather than the probability that it becomes active. The milestones are the steps the source states in its proof, in attack order (proof of Theorem 2.7, pp. 13–14):

  1. for t>0t>0t>0, 0≤(1−e−t)/t≤10\le(1-e^{-t})/t\le10≤(1−e−t)/t≤1, so that an edge can be put into a random set KKK of potential edges with probability (1−e−xg)/xg(1-e^{-x_g})/x_g(1−e−xg​)/xg​;
  2. drawing KKK this way and using the matchings inside KKK as the family is a randomized greedy scheme;
  3. an edge g′g'g′ with ends u,vu,vu,v is selectable exactly when g′∈Kg'\in Kg′∈K and A∩K∩(δ(u)∪δ(v)∖{g′})=∅A\cap K\cap(\delta(u)\cup\delta(v)\setminus\{g'\})=\varnothingA∩K∩(δ(u)∪δ(v)∖{g′})=∅;
  4. the probability of that event equals 1−e−xg′xg′ e−∑g∈δ(u)∪δ(v)∖{g′}xg\frac{1-e^{-x_{g'}}}{x_{g'}}\,e^{-\sum_{g\in\delta(u)\cup\delta(v)\setminus\{g'\}}x_g}xg′​1−e−xg′​​e−∑g∈δ(u)∪δ(v)∖{g′}​xg​;
  5. for x∈bPGx\in bP_Gx∈bPG​, ∑g∈δ(u)∪δ(v)∖{g′}xg≤2(b−xg′)\sum_{g\in\delta(u)\cup\delta(v)\setminus\{g'\}}x_g\le2(b-x_{g'})∑g∈δ(u)∪δ(v)∖{g′}​xg​≤2(b−xg′​);
  6. for t>0t>0t>0, et(et−1)/t≥1e^{t}(e^{t}-1)/t\ge1et(et−1)/t≥1.

Remark (p. 14): the deterministic variant

Taking K=EK=EK=E, the deterministic greedy scheme whose family is the set of all matchings is (b,(1−b)2)(b,(1-b)^2)(b,(1−b)2)-selectable for PGP_GPG​, b∈[0,1]b\in[0,1]b∈[0,1]:

Pr⁡A∼R(x)[g is selectable for A] ≥ (1−b)2∀x∈bPG, ∀g∈E.\Pr_{A\sim R(x)}\big[g\text{ is selectable for }A\big]\ \ge\ (1-b)^2\qquad\forall x\in bP_G,\ \forall g\in E.A∼R(x)Pr​[g is selectable for A] ≥ (1−b)2∀x∈bPG​, ∀g∈E.

This companion is weaker than Theorem 2.7 for every b∈(0,1]b\in(0,1]b∈(0,1] and is not used by it (Remark, p. 14).

What the result contributes

The theorem gives a constant probability bound against the paper's strongest, “almighty” arrival-order adversary. The guarantee is available from vertex degree inequalities alone. Consequently it applies to points of the exact matching polytope as well, without asking the online rounding procedure to exploit its additional odd-set inequalities. The paper places this result beside analogous existence claims for matroids and knapsacks in Theorem 1.8, and uses greedy schemes as components in more general selection results (Theorem 1.8, p. 4).

The existence theorem is proved in the source; it is an open formalization target here. A completed Lean development would supply a reusable finite model of independent active sets, random greedy families, edge selectability, and matching degree relaxations. It would also put the exact e−2be^{-2b}e−2b constant and its endpoint cases under machine checking. The definitions can support future formalizations of the paper's combination and application results, although this mission claims only Theorem 2.7 and its local proof milestones.

Where the argument is delicate

The local degree constraints bound the fractional weight at each endpoint separately, while selectability must survive conflicts at either endpoint under any arrival order. Merely ensuring that the algorithm outputs a matching gives no lower bound for a particular edge: the adversary can expose competing active edges first. The theorem therefore asks for a stronger event that certifies an edge could be added after every eligible previously selected matching. The numerical guarantee must hold uniformly over all edges, even when a coordinate is zero (Definitions 1.4–1.5, pp. 3–4; proof of Theorem 2.7, p. 13).

Formalization scope

The ground set is the finite edge type EEE. The published EdmondsMatching65.Polyhedron.Graph structure represents a graph by a map from edges to unordered pairs of vertices and includes a proof excluding self-loops; its matching predicate is reused directly. This admits parallel edges, a harmless generalization of the page's graph language. Incident edge sets and matchings are finite sets. The polytope is the set of real edge vectors satisfying precisely the nonnegativity and vertex degree inequalities above. Scaling is the pointwise set operation bPGbP_GbPG​, so b=0b=0b=0 remains meaningful. The parameter range 0≤b≤10\le b\le10≤b≤1 is explicit.

Probability is represented by finite sums of product weights. A randomized scheme is a weight function on finite greedy families, independent of the active-set draw. Each family contains the empty set and lies inside the matchings; allowing an empty family would make selectability vacuous. The degree inequalities keep relevant coordinates in [0,1][0,1][0,1], so active-set weights are genuine probabilities. The source's ratio (1−e−xg)/xg(1-e^{-x_g})/x_g(1−e−xg​)/xg​ is given its continuous value 111 at xg=0x_g=0xg​=0. For inputs outside the nonnegative domain, which the theorem never uses, that inclusion probability is also assigned 111 to keep the distribution well formed. These choices and the loopless graph convention are recorded with the draft statements.

The formalization needs finite product probability identities, elementary exponential inequalities, and set facts about matchings and incident edges. Contributions proving the six listed milestones, the deterministic Remark, or the full existential theorem are welcome. The finite model of active sets, greedy families and selectability is shared with the other missions of this series (matroids, knapsack, combinations of schemes, submodular rounding).

Selected references

  • Moran Feldman, Ola Svensson, and Rico Zenklusen, Online Contention Resolution Schemes, arXiv:1508.00142v2, 2015, preprint. Published as Online Contention Resolution Schemes with Applications to Bayesian Selection Problems, SIAM Journal on Computing, 2021.
11 thms1 active userReviewed
Algorithmic Game TheoryCombinatoricsOperations Research·Captain: mikedeng1

Online Contention Resolution Schemes with Applications to Bayesian Selection Problems 5: A (b, c)-Selectable Greedy OCRS Rounds x ∈ bP to E[f(S)] ≥ c·F(x) for Monotone Submodular fResearch Paper

Motivation

Many online selection problems ask an algorithm to accept or reject items one at a time, irrevocably, while keeping the accepted set feasible: independent in a matroid, a matching in a graph, or within a knapsack budget. In Bayesian versions, such as prophet inequalities and posted-price mechanisms, a fractional point xxx of a relaxation of the feasible sets is available in advance, and each item is active independently with probability xex_exe​. An online contention resolution scheme (OCRS), introduced by Feldman, Svensson and Zenklusen, turns this random active set into a feasible selection as items arrive.

The quality of an OCRS is measured by its selectability: the probability that every element is accepted when it is active, whatever the arrival order. Applications, however, care about the value of the selected set, often a submodular function such as a coverage or a utility with diminishing returns. Theorem 1.10 of the paper is the bridge between the two: a selectability guarantee ccc yields an approximation guarantee ccc for every monotone submodular objective, and c/4c/4c/4 for non-monotone ones, even against an adversary who sees every random outcome before fixing the order. The other four missions of this series construct selectable OCRSs for matroids, matchings and knapsacks and combine them; this mission formalizes what such a guarantee is worth.

The offline counterpart is the theory of contention resolution schemes (CRSs) of Chekuri, Vondrák and Zenklusen (SIAM J. Comput. 2014), whose Theorem 1.3 shows that a monotone balanced CRS rounds the multilinear extension of a submodular function with the balancedness factor. The proof of Theorem 1.10 reduces the online statement to that offline result.

Setting

Let NNN be a finite ground set and F⊆2N\mathcal F\subseteq 2^NF⊆2N a down-closed family of feasible sets. A polytope P⊆[0,1]NP\subseteq[0,1]^NP⊆[0,1]N is a relaxation of F\mathcal FF if it contains exactly the characteristic vectors 1I\mathbf 1_I1I​ with I∈FI\in\mathcal FI∈F among the {0,1}\{0,1\}{0,1}-points. For x∈[0,1]Nx\in[0,1]^Nx∈[0,1]N, the random set R(x)R(x)R(x) contains every eee independently with probability xex_exe​.

A greedy OCRS fixes, for each input xxx, a down-closed family Fx⊆F\mathcal F_x\subseteq\mathcal FFx​⊆F containing ∅\emptyset∅ (deterministically, or drawn from a distribution independent of R(x)R(x)R(x)). The elements then arrive one by one; an arriving active element is selected if, together with the elements already selected, it forms a set in Fx\mathcal F_xFx​. An element eee is selectable for the active set AAA if I∪{e}∈FxI\cup\{e\}\in\mathcal F_xI∪{e}∈Fx​ for every I⊆AI\subseteq AI⊆A with I∈FxI\in\mathcal F_xI∈Fx​. For b,c∈[0,1]b,c\in[0,1]b,c∈[0,1] the OCRS is (b,c)(b,c)(b,c)-selectable if for every x∈bPx\in bPx∈bP and every e∈Ne\in Ne∈N

Pr⁡[I∪{e}∈Fx  ∀I⊆R(x), I∈Fx]≥c.\Pr\big[I\cup\{e\}\in\mathcal F_x\ \ \forall I\subseteq R(x),\ I\in\mathcal F_x\big]\ge c .Pr[I∪{e}∈Fx​  ∀I⊆R(x), I∈Fx​]≥c.

The arrival order is chosen by the almighty adversary, which knows R(x)R(x)R(x), the realized family Fx\mathcal F_xFx​ and every random bit of the algorithm before choosing it.

For f:2N→Rf:2^N\to\mathbb Rf:2N→R, the multilinear extension is F(x)=E[f(R(x))]F(x)=\mathbb E[f(R(x))]F(x)=E[f(R(x))]. The function fff is submodular if f(A∪B)+f(A∩B)≤f(A)+f(B)f(A\cup B)+f(A\cap B)\le f(A)+f(B)f(A∪B)+f(A∩B)≤f(A)+f(B) for all A,BA,BA,B, and monotone if A⊆BA\subseteq BA⊆B implies f(A)≤f(B)f(A)\le f(B)f(A)≤f(B). For a set SSS, S(1/2)S(1/2)S(1/2) keeps every element of SSS independently with probability 1/21/21/2.

An (offline) CRS π\piπ for PPP maps each S⊆NS\subseteq NS⊆N, possibly at random, to a subset π(S)⊆S\pi(S)\subseteq Sπ(S)⊆S with 1π(S)∈P\mathbf 1_{\pi(S)}\in P1π(S)​∈P. It is (b,c)(b,c)(b,c)-balanced if Pr⁡[e∈π(R(x))∣e∈R(x)]≥c\Pr[e\in\pi(R(x))\mid e\in R(x)]\ge cPr[e∈π(R(x))∣e∈R(x)]≥c for x∈bPx\in bPx∈bP and xe>0x_e>0xe​>0, and monotone if Pr⁡[e∈π(S1)]≥Pr⁡[e∈π(S2)]\Pr[e\in\pi(S_1)]\ge\Pr[e\in\pi(S_2)]Pr[e∈π(S1​)]≥Pr[e∈π(S2​)] whenever e∈S1⊆S2e\in S_1\subseteq S_2e∈S1​⊆S2​. The characteristic CRS of a greedy OCRS is πˉ(A)={e∈A:e selectable for A}\bar\pi(A)=\{e\in A : e \text{ selectable for } A\}πˉ(A)={e∈A:e selectable for A}.

Formalization targets

Goal: Theorem 1.10

Let PPP be a relaxation of F\mathcal FF, let π\piπ be a (b,c)(b,c)(b,c)-selectable greedy OCRS for PPP, let x∈bPx\in bPx∈bP, and let f≥0f\ge0f≥0 be submodular. Against every almighty adversary, the selected set SSS satisfies

E[f(S)]≥c⋅F(x)if f is monotone,E[f(S(1/2))]≥c4⋅F(x)in general,\mathbb E[f(S)]\ge c\cdot F(x)\quad\text{if } f \text{ is monotone},\qquad \mathbb E[f(S(1/2))]\ge \frac c4\cdot F(x)\quad\text{in general},E[f(S)]≥c⋅F(x)if f is monotone,E[f(S(1/2))]≥4c​⋅F(x)in general,

where in the second case the coins defining S(1/2)S(1/2)S(1/2) are known to the adversary. The constants are the paper's and are not parameters to be improved here.

Milestones

  1. Observation 3.3: for every greedy family, active set and order, πˉ(A)\bar\pi(A)πˉ(A) is contained in the selected set.
  2. The remark after Definition 3.2: πˉ\bar\piπˉ is a CRS for PPP.
  3. Lemma 3.4: the characteristic CRS of a (b,c)(b,c)(b,c)-selectable greedy OCRS is (b,c)(b,c)(b,c)-balanced and monotone.
  4. Lemma 3.5 (from Chekuri–Vondrák–Zenklusen): for every non-negative submodular fff there is a pruning map ηf\eta_fηf​, with ηf(S)⊆S\eta_f(S)\subseteq Sηf​(S)⊆S, such that E[f(ηf(π(R(x))))]≥c⋅F(x)\mathbb E[f(\eta_f(\pi(R(x))))]\ge c\cdot F(x)E[f(ηf​(π(R(x))))]≥c⋅F(x) for every monotone (b,c)(b,c)(b,c)-balanced CRS π\piπ and every x∈bPx\in bPx∈bP.
  5. Lemma 3.6 (Feige–Mirrokni–Vondrák, Lemma 2.2): if Tp⊆TT_p\subseteq TTp​⊆T contains every element of TTT with probability ppp, not necessarily independently, then E[g(Tp)]≥(1−p)g(∅)+p g(T)\mathbb E[g(T_p)]\ge(1-p)g(\emptyset)+p\,g(T)E[g(Tp​)]≥(1−p)g(∅)+pg(T) for submodular ggg.
  6. Lemma 3.7 (Buchbinder–Feldman–Naor–Schwartz, Lemma 2.2): if NpN_pNp​ contains every element with probability at most ppp, then E[g(Np)]≥(1−p)g(∅)\mathbb E[g(N_p)]\ge(1-p)g(\emptyset)E[g(Np​)]≥(1−p)g(∅) for non-negative submodular ggg.

Significance

Theorem 1.10 is what makes selectability useful. Combined with the paper's OCRSs for matroids, matchings and knapsacks (Theorem 1.8) and the combination rule (Theorem 1.9), it gives constant-factor online algorithms for submodular objectives under intersections of these constraints, and through them the paper's applications to prophet inequalities, posted-price mechanisms and stochastic probing. The guarantee holds against the strongest adversary considered, so it applies whatever arrival order the application imposes.

The theorem and its proof are published; none of the statements of this mission has a machine-checked proof. Formalizing it produces a reusable interface between online and offline contention resolution, and machine-checked versions of two standard lemmas on correlated random subsets of submodular functions (Lemmas 3.6 and 3.7) used throughout submodular maximization. Lemma 3.5 is the substantial imported ingredient: its formal proof amounts to formalizing the rounding theorem of Chekuri, Vondrák and Zenklusen for non-monotone objectives.

Difficulty

The first idea is to compare the online selection directly with F(x)F(x)F(x) element by element: each element is selected with probability at least c xec\,x_ecxe​. For a linear objective this suffices. For a submodular objective it does not, because the marginal value of an element depends on which other elements are selected, and the adversary controls this through the order. The reduction to the characteristic CRS removes the order, but the offline rounding result (Lemma 3.5) then needs the full CRS machinery, including a pruning step that depends on fff alone.

For non-monotone fff, selecting extra elements can lose value, and the adversary decides which elements beyond πˉ(A)\bar\pi(A)πˉ(A) are selected, after seeing the coins. Controlling this requires the correlated-sampling bounds of Lemmas 3.6 and 3.7 rather than independence.

Formalization scope

The ground set is a finite type α : Type with decidable equality; sets are Finset α; vectors are α → ℝ. All probabilities and expectations are exact finite sums of product weights, without measure theory. bPbPbP is the pointwise scaling b • P, which handles b=0b=0b=0. The multilinear extension and submodularity are the published definitions NonmonotoneSubmod.Shared.F and NonmonotoneSubmod.Shared.Submodular.

A randomized greedy OCRS is a weight function on families Finset (Finset α), a CRS a weight function on maps Finset α → Finset α; a deterministic OCRS is a point mass. The almighty adversary is any function from the realized family and active set (and, in part 2, the coin set CCC, with S(1/2)=S∩CS(1/2)=S\cap CS(1/2)=S∩C and CCC uniform and independent of everything else) to an order listing every element once. The standing assumptions of Definition 1.3 are explicit hypotheses: F\mathcal FF down-closed, PPP a polytope contained in [0,1]N[0,1]^N[0,1]N and a relaxation of F\mathcal FF, and b,c∈[0,1]b,c\in[0,1]b,c∈[0,1]. A polytope is encoded as the convex hull of finitely many vectors. Every greedy family contains ∅\emptyset∅, so the empty family cannot make every element vacuously selectable, and the inputs x∈bPx\in bPx∈bP lie in [0,1]N[0,1]^N[0,1]N, so the product weights are probabilities.

Definition 3.1 is printed with "≥c⋅xe\ge c\cdot x_e≥c⋅xe​"; the formalization uses "≥c\ge c≥c", as the proof of Lemma 3.4 and Definition 4.5 do. Lemma 3.5 keeps the paper's quantifier order: ηf\eta_fηf​ is fixed before PPP, bbb, ccc, π\piπ and xxx.

Welcome contributions: proofs of the milestones in any order; Lemmas 3.6 and 3.7 are self-contained and reusable; Lemma 3.5 may be split further along the proof of Chekuri–Vondrák–Zenklusen.

Selected references

  • M. Feldman, O. Svensson, R. Zenklusen, Online Contention Resolution Schemes, arXiv:1508.00142v2, 2015; Online Contention Resolution Schemes with Applications to Bayesian Selection Problems, SIAM J. Comput. 50(2), 2021. https://arxiv.org/abs/1508.00142v2
  • C. Chekuri, J. Vondrák, R. Zenklusen, Submodular Function Maximization via the Multilinear Relaxation and Contention Resolution Schemes, SIAM J. Comput. 43(6), 2014. https://doi.org/10.1137/110839655
  • U. Feige, V. Mirrokni, J. Vondrák, Maximizing Non-Monotone Submodular Functions, SIAM J. Comput. 40(4), 2011. https://doi.org/10.1137/090779346
  • N. Buchbinder, M. Feldman, J. Naor, R. Schwartz, Submodular Maximization with Cardinality Constraints, SODA 2014. https://doi.org/10.1137/1.9781611973402.106
12 thms1 active userReviewed
Algorithmic Game TheoryCombinatoricsOperations Research·Captain: mikedeng1

Online Contention Resolution Schemes with Applications to Bayesian Selection Problems 4: Intersecting the Families of (b, c₁)- and (b, c₂)-Selectable Greedy OCRSs Is (b, c₁c₂)-SelectableResearch Paper

Motivation

Online selection problems ask an algorithm to accept or reject elements as they arrive, while keeping the accepted set feasible. The decisions are irrevocable, and the arrival order can be chosen after the random activity of the elements is known. A guarantee for one constraint is therefore useful only if it remains meaningful when several constraints apply at once. In Feldman, Svensson and Zenklusen's paper, the constraints include matroid independence, matching, and knapsack capacity. A set satisfying two such requirements lies in the intersection of their feasible families. Theorem 1.9 establishes how their online guarantees combine: the two selectability factors multiply.

The result belongs to a sequence of statements in the paper. The authors first define greedy online contention resolution schemes and their selectability guarantee in Definitions 1.3–1.5. They then state the combination theorem in the introduction and identify the intersection construction in Definition 2.10. Lemma 2.11 supplies the precise selectability assertion for that construction. This mission isolates those claims so that later formalizations of specific constraints can reuse them without repeating the probability argument.

Setting

Let NNN be a finite ground set and F⊆2N\mathcal F\subseteq 2^NF⊆2N a family of feasible sets. A fractional point x∈[0,1]Nx\in[0,1]^Nx∈[0,1]N determines a random active set R(x)R(x)R(x): each e∈Ne\in Ne∈N belongs to R(x)R(x)R(x) independently with probability xex_exe​. Thus a particular set A⊆NA\subseteq NA⊆N has weight

Pr⁡[R(x)=A]=∏e∈Axe∏e∉A(1−xe).\Pr[R(x)=A]=\prod_{e\in A}x_e\prod_{e\notin A}(1-x_e).Pr[R(x)=A]=e∈A∏​xe​e∈/A∏​(1−xe​).

A greedy online contention resolution scheme, or greedy OCRS, chooses a down-closed family Gx⊆F\mathcal G_x\subseteq\mathcal FGx​⊆F. As active elements arrive, it accepts an element when adding that element to those already accepted remains in Gx\mathcal G_xGx​. In a randomized scheme, the choice of Gx\mathcal G_xGx​ is random and independent of R(x)R(x)R(x). The family contains the empty set, as the procedure begins with no accepted elements.

An element eee is selectable for a realized active set AAA and family Gx\mathcal G_xGx​ if every I⊆AI\subseteq AI⊆A with I∈GxI\in\mathcal G_xI∈Gx​ still permits eee: I∪{e}∈GxI\cup\{e\}\in\mathcal G_xI∪{e}∈Gx​. This condition covers every possible set that might have been accepted before eee arrives. A greedy OCRS is (b,c)(b,c)(b,c)-selectable for a relaxation PPP when the probability of this event is at least ccc for every x∈bPx\in bPx∈bP and every e∈Ne\in Ne∈N. The probability includes both draws when the family is randomized. Here b,c∈[0,1]b,c\in[0,1]b,c∈[0,1] and bP={by:y∈P}bP=\{by:y\in P\}bP={by:y∈P}.

For two feasible families F1,F2\mathcal F_1,\mathcal F_2F1​,F2​, the combination of their greedy OCRSs independently draws G1,x\mathcal G_{1,x}G1,x​ and G2,x\mathcal G_{2,x}G2,x​, then uses Gx=G1,x∩G2,x\mathcal G_x=\mathcal G_{1,x}\cap\mathcal G_{2,x}Gx​=G1,x​∩G2,x​. This is a family of sets satisfying both constraints. The proof also uses χe(A,F,F′)\chi_e(A,F,F')χe​(A,F,F′), the indicator that I∪{e}∈FI\cup\{e\}\in FI∪{e}∈F for every I⊆AI\subseteq AI⊆A lying in F′F'F′. For selectability itself, F′=FF'=FF′=F.

Formalization targets

Goal: Theorem 1.9

For P1,P2⊆[0,1]NP_1,P_2\subseteq[0,1]^NP1​,P2​⊆[0,1]N, let the first scheme be (b,c1)(b,c_1)(b,c1​)-selectable and the second (b,c2)(b,c_2)(b,c2​)-selectable. The mission asks for a greedy OCRS on the intersection with

Pr⁡[e selectable at x]≥c1c2for all x∈b(P1∩P2), e∈N.\Pr[e\text{ selectable at }x]\ge c_1c_2 \quad\text{for all }x\in b(P_1\cap P_2),\ e\in N.Pr[e selectable at x]≥c1​c2​for all x∈b(P1​∩P2​), e∈N.

The target records existence of the combined scheme, independent of a particular proof of its selectability. The paper also claims preservation of polynomial-time efficiency; the finite mathematical model here does not assign a complexity measure to a family sampler or membership oracle, so that clause is outside this target.

Construction: Lemma 2.11

The lemma names the combined scheme: independently draw both greedy families and use their intersection. It asserts the same (b,c1c2)(b,c_1c_2)(b,c1​c2​) bound for that specific scheme. The earlier milestones record the two order properties of χe\chi_eχe​ and the finite product-measure correlation inequality appearing in the lemma's proof. Their statements keep the original quantifiers over active sets and families visible.

Significance

The theorem turns guarantees for individual constraints into a guarantee for simultaneous constraints. If a matroid scheme and a matching scheme meet the same scale parameter bbb, their intersection is governed by the product of their selectability factors. The factor loss is explicit, which lets subsequent algorithmic results track it without redefining the combined scheme. The lemma additionally fixes the actual distribution over intersected families, avoiding ambiguity about how two randomized procedures are coupled.

Mathlib already provides a general finite FKG inequality, but it does not state the paper's OCRS objects, the selectable event, or the intersection theorem. This mission adds those interfaces and the specialized correlation statement. The source results are proved in the cited paper; the local theorem files in this proposal are statements for later machine-checked proofs. A complete formalization would connect the finite active-set weights to Mathlib's lattice correlation theorem and then establish the randomized family bookkeeping. The finite probability definitions can also serve other OCRS results in the same paper.

Difficulty

The obvious argument would multiply the two probabilities that eee is selectable for the separate families. Both events depend on the same active set R(x)R(x)R(x), so they are not independent even when the families are drawn independently. Moreover, selectability for the intersection checks only sets belonging to both families, whereas the original bounds check each family separately. The formal statements must keep these two dependencies distinct. The paper's two monotonicity claims and its FKG display identify the needed intermediate assertions, but their use with an arbitrary random choice of families still requires careful finite-sum reasoning.

Formalization scope

The ground set is a finite Lean type and subsets are Finsets. Feasibility is a predicate on those finite sets. Fractional points are real-valued functions, and each input relaxation is represented as a set of such functions with an explicit coordinate bound 0≤ye≤10\le y_e\le10≤ye​≤1. The scale and selectability parameters belong to [0,1][0,1][0,1]. The point xxx ranges over b(P1∩P2)b(P_1\cap P_2)b(P1​∩P2​), not over an independently chosen intersection of scaled sets, so the case b=0b=0b=0 retains the paper's meaning. The finite products define the law of R(x)R(x)R(x); no measure-theoretic integral is used.

A randomized scheme is a probability weight on possible greedy families for every input point. The combined weight is the product of the two original weights summed over pairs with the required intersection, expressing independent draws. The family distribution is independent of the active-set draw. This sampling convention is necessary for the stated product guarantee and is made explicit here although Definition 2.10 does not spell it out. Each greedy family contains ∅\varnothing∅ and lies in its own feasible family. These requirements exclude the empty-family and unconstrained-family encodings that would make selectability trivial. The unit-cube bounds prevent negative active-set weights.

The Lean declarations focus on the finite selectability guarantee. They abstract a polytope by its set of fractional points and use the paper's unit-cube condition; they do not formalize a polyhedral representation or test that its binary points equal the feasible sets. A full development of the paper's algorithmic setting would additionally represent online arrival orders, family sampling algorithms, membership oracles, and polynomial running time. Contributions establishing the finite probability lemmas, the specialized FKG statement, and Lemma 2.11 are directly useful here. The general active-set and selectability definitions are reusable across the paper's other missions.

Selected references

  • Moran Feldman, Ola Svensson and Rico Zenklusen, Online Contention Resolution Schemes, arXiv preprint, 2015, arXiv:1508.00142v2. Definitions 1.3–1.5, Theorem 1.9, Definition 2.10 and Lemma 2.11.
  • Moran Feldman, Ola Svensson and Rico Zenklusen, Online Contention Resolution Schemes with Applications to Bayesian Selection Problems, SIAM Journal on Computing 50(2), 2021, DOI:10.1137/18M1226130. The mission's statement indices follow the arXiv version above.
7 thms1 active userReviewed
CombinatoricsOperations ResearchOptimization·Captain: mikedeng1

Online Contention Resolution Schemes with Applications to Bayesian Selection Problems 3: A Knapsack Polytope Has a (b, (1−2b)/(2−2b))-Selectable Randomized Greedy OCRS for b ≤ 1/2Research Paper

Motivation

In online selection, useful items arrive one at a time and must be accepted or rejected before the next arrival. A knapsack constraint limits the total size of accepted items to one unit. When the order may depend on all of the random outcomes, an algorithm needs a guarantee that survives every such order. Feldman, Svensson, and Zenklusen introduced online contention resolution schemes (OCRSs) to provide this guarantee for fractional relaxations of selection problems. Their knapsack theorem gives an explicit guarantee for a randomized greedy scheme and contrasts with a size-dependent obstruction for deterministic greedy schemes (Feldman–Svensson–Zenklusen, §2.3).

The knapsack result is one of the three concrete existence results summarized in Theorem 1.8 of the paper. The surrounding framework is designed to convert fractional solutions into feasible online selections under an almighty adversary, who knows the active items and the scheme's random choices before choosing their arrival order. This makes the selectability requirement stronger than a performance claim for one fixed order (Feldman–Svensson–Zenklusen, Definitions 1.3–1.5).

Setting

Let NNN be a finite ground set. Each element e∈Ne\in Ne∈N has a size se∈[0,1]s_e\in[0,1]se​∈[0,1]. A subset I⊆NI\subseteq NI⊆N is feasible when its total size is at most one. The natural fractional polytope is

P={x∈[0,1]N:∑e∈Nsexe≤1}.P=\left\{x\in[0,1]^N:\sum_{e\in N}s_ex_e\le1\right\}.P={x∈[0,1]N:e∈N∑​se​xe​≤1}.

For a vector x∈[0,1]Nx\in[0,1]^Nx∈[0,1]N, the random active set R(x)R(x)R(x) includes each element eee independently with probability xex_exe​. A greedy family FxF_xFx​ is a downward-closed family of feasible subsets. The online rule accepts an active arriving element when adding it to the accepted set keeps that set inside FxF_xFx​. A randomized greedy OCRS chooses the family at random, independently of R(x)R(x)R(x). The paper studies the probability that an element eee is selectable: for every possible accepted set I⊆R(x)I\subseteq R(x)I⊆R(x) with I∈FxI\in F_xI∈Fx​, the enlarged set I∪{e}I\cup\{e\}I∪{e} still belongs to FxF_xFx​. Such an event makes eee acceptable regardless of the arrival order (Feldman–Svensson–Zenklusen, §1.1).

Write bP={bp:p∈P}bP=\{bp:p\in P\}bP={bp:p∈P}. An OCRS is (b,c)(b,c)(b,c)-selectable when, for every x∈bPx\in bPx∈bP and every e∈Ne\in Ne∈N, the probability that eee is selectable is at least ccc. The paper divides elements at size 1/21/21/2: Nbig={e:se>1/2}N_{\mathrm{big}}=\{e:s_e>1/2\}Nbig​={e:se​>1/2}, while equality belongs to the small elements. Its construction randomly chooses between the feasible subsets using only big elements and those using only small elements (Feldman–Svensson–Zenklusen, proof of Theorem 2.9).

Formalization targets

Main theorem

Theorem 2.9 asserts that for every b∈[0,1/2]b\in[0,1/2]b∈[0,1/2] and every such knapsack polytope, a randomized greedy OCRS exists with

c=1−2b2−2b.c=\frac{1-2b}{2-2b}.c=2−2b1−2b​.

The goal quantifies over every finite ground set and every allowed size vector. At the endpoint b=1/2b=1/2b=1/2, the claimed lower bound is zero; at b=0b=0b=0, it is 1/21/21/2. The denominator remains positive throughout the stated interval. The goal is an existence statement; its witnesses are arbitrary randomized greedy families meeting the definition, so the paper's particular construction is named separately in the supporting definitions (Feldman–Svensson–Zenklusen, Theorem 2.9, p. 14).

Supporting targets

The milestones include the paper's bound 0≤bbig≤b0\le b_{\mathrm{big}}\le b0≤bbig​≤b for bbig=∑e∈Nbigsexeb_{\mathrm{big}}=\sum_{e\in N_{\mathrm{big}}}s_ex_ebbig​=∑e∈Nbig​​se​xe​, the big and small element selectability bounds for the paper's construction, and Proposition 2.8. That proposition says that, for every ground-set size n≥1n\ge1n≥1, some knapsack constraint defeats every deterministic greedy OCRS whenever c>(1−b)n−1c>(1-b)^{n-1}c>(1−b)n−1 for b∈[0,1]b\in[0,1]b∈[0,1] (Feldman–Svensson–Zenklusen, §2.3, pp. 14–15).

Significance

The theorem supplies a concrete selectable OCRS for the natural fractional knapsack relaxation. Its guarantee applies to every element and every fractional point in bPbPbP, under the paper's strongest arrival-order adversary. Proposition 2.8 explains why randomization matters here: a deterministic greedy scheme cannot retain an analogous positive constant independent of the number of elements on every knapsack instance. These results serve as the knapsack component in the paper's broader OCRS framework and can be combined with schemes for other constraints in the paper's later combination result (Feldman–Svensson–Zenklusen, Theorem 1.9).

The mathematical theorem and proposition have proofs in the cited paper. This mission asks for machine-checked proofs of their formal statements and of the stated intermediate bounds. The current draft declarations contain proof placeholders; compilation establishes only that their statements and definitions are well formed. A completed development would also leave reusable finite active-set probabilities, greedy-family predicates, and explicit knapsack constructions for later online-rounding work.

Difficulty

A family that contains every feasible subset is not automatically useful for online selection: one active set may contain several subsets whose accepted elements block the tested element. The requirement quantifies over every feasible accepted subset of the realized active set, reflecting the almighty adversary rather than a fixed order. In a knapsack, elements larger than one half behave differently from smaller ones, and a single deterministic family cannot give a ground-set-independent selectability guarantee on all instances. The randomized construction must assign valid probabilities to its two families and obtain the same lower bound in both size regimes (Feldman–Svensson–Zenklusen, Proposition 2.8 and Theorem 2.9).

Formalization scope

The Lean ground set is a finite type, and subsets are Finsets. The polytope is the set of real vectors satisfying both the box inequalities and the unit-capacity inequality. The feasible family uses the same sizes and unit capacity. Active-set and family-choice probabilities are finite sums, with independent Bernoulli factors for R(x)R(x)R(x). The scaling bPbPbP is pointwise set scalar multiplication, so the case b=0b=0b=0 retains its intended meaning. The size assumptions se∈[0,1]s_e\in[0,1]se​∈[0,1] and the ranges of bbb are explicit theorem hypotheses. No nonemptiness assumption is added to the ground set; the statements include the empty case.

The greedy-family definition requires the empty set, making it impossible to satisfy selectability by using an empty family. Every family chosen with positive probability is downward closed and consists of feasible sets. The polytope's box constraints keep every active-set weight a probability. The construction's mixing probability is used exactly on x∈bPx\in bPx∈bP; on other vectors its totalized value chooses the small family, where the paper makes no performance claim. The two branch weights are added if their families coincide, as they can on an empty ground set.

One display on page 15 includes the tested element in a load sum, while the next displayed equality treats that element as excluded. The pointwise supporting item uses the excluded-element form needed for the stated argument; the printed display is preserved in the milestone quotation and the discrepancy is detailed in MODERATION_NOTES.md. Contributions toward either size regime, the finite probability bounds, or the deterministic obstruction are welcome. The finite probability and greedy-family definitions are the parts most directly reusable beyond knapsack.

Selected references

  • Moran Feldman, Ola Svensson, and Rico Zenklusen, Online Contention Resolution Schemes, arXiv:1508.00142v2, 2015, §1.1 and §2.3, pp. 2–4 and 14–15. arXiv.
8 thms1 active userReviewed
Algorithmic Game TheoryCombinatoricsOperations Research·Captain: mikedeng1

Online Contention Resolution Schemes with Applications to Bayesian Selection Problems 1: Every Loopless Matroid Polytope Has a (b, 1−b)-Selectable Deterministic Greedy OCRSResearch Paper

Motivation

In online selection, items reveal information one at a time, and each decision to accept an item is irrevocable. A matroid constraint permits some subsets and excludes others; it covers choices such as selecting linearly independent vectors or acyclic edges. A fractional solution can assign an acceptance probability to each item, but independent activations can produce a set that violates the constraint. An online contention resolution scheme (OCRS) turns these activations into a feasible selection while respecting the order in which items arrive. The matroid result of Feldman, Svensson and Zenklusen gives a guarantee that survives even when an adversary knows all random outcomes before choosing the arrival order.

This mission concerns the paper's matroid theorem, Theorem 2.1. The same paper develops schemes for matching and knapsack constraints and a rule for combining schemes for intersecting constraints. Those results are separate missions in this series. Here the central question is how much of a fractional matroid point can be retained by a deterministic greedy rule as each active item is revealed.

Setting

Let NNN be a finite ground set and let MMM be a matroid on NNN. Its independent subsets form a down-closed family F\mathcal FF: if III is independent, every subset of III is independent. The rank r(S)r(S)r(S) of S⊆NS\subseteq NS⊆N is the maximum size of an independent subset of SSS. The matroid independence polytope is

PF={x∈R≥0N:x(S):=∑e∈Sxe≤r(S) for all S⊆N}.P_{\mathcal F}=\left\{x\in\mathbb R^N_{\ge0}:x(S):=\sum_{e\in S}x_e\le r(S)\text{ for all }S\subseteq N\right\}.PF​={x∈R≥0N​:x(S):=e∈S∑​xe​≤r(S) for all S⊆N}.

Its points have 0≤xe≤10\le x_e\le10≤xe​≤1. For a point xxx, the active set R(x)R(x)R(x) contains each eee independently with probability xex_exe​. A deterministic greedy OCRS assigns a down-closed family Fx⊆F\mathcal F_x\subseteq\mathcal FFx​⊆F to xxx and accepts an active element when adding it to the already accepted set stays in Fx\mathcal F_xFx​. The empty set belongs to every Fx\mathcal F_xFx​. The rule depends on xxx before the active set and arrival order are revealed.

An element eee is selectable for an active set AAA if I∪{e}∈FxI\cup\{e\}\in\mathcal F_xI∪{e}∈Fx​ for every I⊆AI\subseteq AI⊆A with I∈FxI\in\mathcal F_xI∈Fx​. A greedy OCRS is (b,c)(b,c)(b,c)-selectable if this event has probability at least ccc for every x∈bPFx\in bP_{\mathcal F}x∈bPF​ and every e∈Ne\in Ne∈N. The event implies that eee can be added regardless of which eligible active elements arrived first. These definitions are from §1.1 of the paper.

Formalization targets

Matroid OCRS

For every b∈[0,1]b\in[0,1]b∈[0,1] and every finite loopless matroid MMM, the target is a deterministic greedy OCRS with

Pr⁡[I∪{e}∈Fx for all I⊆R(x), I∈Fx]≥1−bfor every x∈bPF, e∈N.\Pr\left[I\cup\{e\}\in\mathcal F_x\text{ for all }I\subseteq R(x),\ I\in\mathcal F_x\right]\ge 1-b \quad\text{for every }x\in bP_{\mathcal F},\ e\in N.Pr[I∪{e}∈Fx​ for all I⊆R(x), I∈Fx​]≥1−bfor every x∈bPF​, e∈N.

The loopless condition is a disclosed correction to the printed Theorem 2.1: a loop cannot belong to any independent singleton, so its selectability is zero. The printed bound would demand positive selectability for that loop whenever b<1b<1b<1. The theorem keeps the paper's constant 1−b1-b1−b exactly.

Supporting targets

The milestones reproduce six statements from §2.1–2.1.1, in attack order: validity of a family defined by a chain ∅=Nℓ⊊⋯⊊N0=N\varnothing=N_\ell\subsetneq\cdots\subsetneq N_0=N∅=Nℓ​⊊⋯⊊N0​=N of matroid minors (independence of the union of layer-wise independent sets is Soto, Theorem 5.1); a characterization of selectability through span in one layer; monotonicity of the first chain set as coordinates of xxx increase; the description of the base polytope PB={x∈PF:x(N)=r(N)}P_{\mathcal B}=\{x\in P_{\mathcal F}:x(N)=r(N)\}PB​={x∈PF​:x(N)=r(N)} as the set of maximal vectors of PFP_{\mathcal F}PF​; the weighted span bound of Lemma 2.4,

∑e∈NxePr⁡[e∈span⁡(R(x)∪S)]≤b⋅(x(N)+(1−b) r(S)),\sum_{e\in N}x_e\Pr[e\in\operatorname{span}(R(x)\cup S)]\le b\cdot\bigl(x(N)+(1-b)\,r(S)\bigr),e∈N∑​xe​Pr[e∈span(R(x)∪S)]≤b⋅(x(N)+(1−b)r(S)),

strict when S≠∅S\neq\varnothingS=∅; and the proper-inclusion conclusion S⊊NS\subsetneq NS⊊N of Corollary 2.5. They expose the distinct matroid and probability claims on which the main theorem rests.

Significance

The theorem gives a fixed greedy family for each fractional point, with a per-element probability guarantee that is independent of the arrival order. Offline contention resolution schemes for matroids were introduced by Chekuri, Vondrák and Zenklusen; the online, almighty-adversary version here is stronger. A solver can use that guarantee as the matroid component in larger online selection problems. In the source paper it feeds the intersection and application results; those applications require further models and statements beyond this mission. The factor 1−b1-b1−b makes the tradeoff explicit: using a more heavily scaled fractional point raises the risk of contention.

The result is proved in the cited paper. This mission asks for a machine-checked Lean proof of its corrected statement and the source's supporting claims. A successful development would add reusable finite representations of independent activation, matroid polytope inequalities, span probabilities, and chain families. It would also make the loop boundary explicit, preventing a formally false version of the printed theorem from being published as a target.

Difficulty

Independent activation does not itself respect matroid independence. A greedy family that simply contains every independent set can make an element unavailable after earlier active elements are accepted, and the arrival order can expose the worst such history. The central issue is therefore a uniform selectability event: the same element must remain addable after every feasible subset of the active set. The paper measures that obstruction through spans in successive matroid minors. The first chain set must be proper so that the construction can continue, which is why the strict part of the weighted span inequality matters.

Formalization scope

The ground set is a finite Lean type, subsets are Finsets, and Mathlib's Matroid supplies independence, closure, contraction and restriction. Matroid rank is converted from Mathlib's extended natural rank to a real number; finiteness makes this conversion finite on every subset. The independence and base polytopes are defined by the paper's inequalities. The scalar multiple bPbPbP is pointwise set scalar multiplication, including b=0b=0b=0.

Probabilities are finite sums of product weights over all active sets. In particular, the event in Lemma 2.4 uses span⁡(R(x)∪S)\operatorname{span}(R(x)\cup S)span(R(x)∪S), whereas the recurrence defining SSS uses span⁡((R(x)∪Si)∖{e})\operatorname{span}((R(x)\cup S_i)\setminus\{e\})span((R(x)∪Si​)∖{e}). These are separate definitions. Every point used for activation lies in [0,1]N[0,1]^N[0,1]N, as enforced by the polytope and b∈[0,1]b\in[0,1]b∈[0,1]. The base-polytope condition on Corollary 2.5 is the standing assumption stated before the lemma in the source.

The loopless condition applies to the goal, to Corollary 2.5, and to the strict clause of Lemma 2.4 only; its non-strict clause, the chain statements and the monotonicity statement hold for every matroid and are stated without it. Three trivializing formalizations are ruled out: every greedy family must contain the empty set, so an empty family cannot satisfy selectability vacuously; every member of a greedy family must be independent in MMM, so Fx=2N\mathcal F_x=2^NFx​=2N is not admissible; and every point used for activation lies in [0,1]N[0,1]^N[0,1]N, so Pr⁡[R(x)=A]\Pr[R(x)=A]Pr[R(x)=A] is a genuine probability. The quantified scheme depends only on xxx, and its success condition ranges over every feasible subset of the realized active set. Contributions that establish the chain family's independence, finite probability identities, and rank or closure facts are within scope. Polynomial-time implementation and the Monte Carlo estimates of Lemma 2.6 are outside this mission's mathematical target.

Selected references

  • Moran Feldman, Ola Svensson and Rico Zenklusen, Online Contention Resolution Schemes, arXiv:1508.00142v2, 2015; published as Online Contention Resolution Schemes with Applications to Bayesian Selection Problems, SIAM Journal on Computing, 2021. Preprint.
  • Chandra Chekuri, Jan Vondrák and Rico Zenklusen, Submodular Function Maximization via the Multilinear Relaxation and Contention Resolution Schemes, SIAM Journal on Computing, 2014. arXiv:1105.4593.
  • José A. Soto, Matroid Secretary Problem in the Random-Assignment Model, SIAM Journal on Computing, 2013. arXiv:1007.2140.
10 thms1 active userReviewed
AnalysisDynamical SystemsStochastic Systems·Captain: mikedeng1

Averaging Dynamics Driven by Fractional Brownian Motion: Young Solutions of dxⁿ = fₙ(t, xⁿ)dB Converge in Probability in C^α to the Solution of dx = f̄(x)dB When |fₙ − f̄|_{−κ,γ} → 0Research Paper

Why oscillating coefficients matter

An equation driven by fractional Brownian motion can respond to a rapidly varying coefficient even when that coefficient has no useful pointwise limit. Hairer and Li study this issue for fractional Brownian motion with Hurst index above one half. Their Section 3 isolates a deterministic averaging statement: if the time antiderivatives of the coefficients approach those of a limiting coefficient in a suitable negative Hölder norm, the solutions approach the limiting solution in probability. This statement supplies a common target for averaging models whose fast variables have very different origins. The authors use it later with mixing fast processes, and also present a quasi-periodic example. The result and its role in the paper are stated in Hairer and Li, Section 3.

The problem sits at the boundary of two familiar views of differential equations. A pathwise Young equation is defined from individual sufficiently regular paths, whereas the available estimates for the oscillating coefficient and the fractional Brownian driver involve probability and conditional expectation. The mission records the precise point where these descriptions meet. It treats the published article, rather than a preprint, and targets its Theorem 3.13.

Driver, coefficients, and solutions

An mmm-dimensional fractional Brownian motion B=(B1,…,Bm)B=(B^1,\ldots,B^m)B=(B1,…,Bm) with index HHH is a centered Gaussian process with independent coordinates and covariance

E[BsiBtj]=1i=j∣s∣2H+∣t∣2H−∣t−s∣2H2.\mathbb E[B_s^iB_t^j]=\mathbf1_{i=j}\frac{|s|^{2H}+|t|^{2H}-|t-s|^{2H}}2.E[Bsi​Btj​]=1i=j​2∣s∣2H+∣t∣2H−∣t−s∣2H​.

The paper works with a two-sided process and 1/2<H<11/2<H<11/2<H<1. It uses a version whose sample paths are locally β\betaβ-Hölder for every β<H\beta<Hβ<H. The mission names that version condition separately, because the covariance law alone does not select particular sample paths. The relevant past sigma algebra at time sss is generated by increments Bu−BvB_u-B_vBu​−Bv​ with u,v≤su,v\le su,v≤s. These conventions come from pp. 1826 and 1831–1832 of the article.

For a path xxx on [0,T][0,T][0,T], write ∣x∣α|x|_\alpha∣x∣α​ for its homogeneous Hölder seminorm and ∣x∣∞|x|_\infty∣x∣∞​ for its supremum norm. The norm on Cα[0,T]C^\alpha[0,T]Cα[0,T] is ∣x∣∞+∣x∣α|x|_\infty+|x|_\alpha∣x∣∞​+∣x∣α​. A Young integral ∫stfr dbr\int_s^t f_r\,db_r∫st​fr​dbr​ is the limit of left-point sums over partitions of [s,t][s,t][s,t] as the largest cell shrinks to zero. A Young solution from x0x_0x0​ is an α\alphaα-Hölder path satisfying its integral equation at every time, with value x0x_0x0​ at time zero. Here BC2BC^2BC2 means bounded, twice continuously differentiable maps with bounded derivatives through order two; Cζ,2C^{\zeta,2}Cζ,2 means ζ\zetaζ-Hölder in time with values in BC2BC^2BC2. These are the paper's conventions on pp. 1829 and 1835–1838.

The negative Hölder bound ∣f∣−κ,γ≤K|f|_{-\kappa,\gamma}\le K∣f∣−κ,γ​≤K means that the time antiderivative of fff has increments bounded by K∣t−s∣1−κK|t-s|^{1-\kappa}K∣t−s∣1−κ, and that the difference of those increments at spatial points x,yx,yx,y is bounded by K∣t−s∣1−κ∣x−y∣γK|t-s|^{1-\kappa}|x-y|^\gammaK∣t−s∣1−κ∣x−y∣γ. Both bounds apply on the time window in question when ∣t−s∣≤1|t-s|\le1∣t−s∣≤1. This is equation (3.12) in the article. It measures cancellation over time; it does not require f(t,x)f(t,x)f(t,x) itself to be small.

Formalization targets

The goal is the dBdBdB part of Theorem 3.13. Let fnf_nfn​ be Cζ,2C^{\zeta,2}Cζ,2, let the limiting fˉ\bar ffˉ​ be autonomous and in BC2BC^2BC2, and suppose ∣fn−fˉ∣−κ,γ→0|f_n-\bar f|_{-\kappa,\gamma}\to0∣fn​−fˉ​∣−κ,γ​→0. The exponents satisfy 1/2<α<H−κ1/2<\alpha<H-\kappa1/2<α<H−κ, H−κ>1/2H-\kappa>1/2H−κ>1/2, H−κ+γα>1H-\kappa+\gamma\alpha>1H−κ+γα>1, and α<ζ≤1\alpha<\zeta\le1α<ζ≤1. For Young solutions with the same starting point,

dxtn=fn(t,xtn) dBt,dxt=fˉ(xt) dBt,dx_t^n=f_n(t,x_t^n)\,dB_t,\qquad dx_t=\bar f(x_t)\,dB_t,dxtn​=fn​(t,xtn​)dBt​,dxt​=fˉ​(xt​)dBt​,

the target is

∀ε>0,P(∥xn−x∥Cα[0,T]>ε)⟶0.\forall\varepsilon>0,\qquad P\bigl(\|x^n-x\|_{C^\alpha[0,T]}>\varepsilon\bigr)\longrightarrow0.∀ε>0,P(∥xn−x∥Cα[0,T]​>ε)⟶0.

The milestone list follows the statements used around this target: the Hölder composition estimate of Lemma 2.1, the RKHS estimate of Lemma 3.3, the fractional Brownian integral estimate of Lemma 3.4, the stochastic sewing lemma of Lemma 3.5, the Young case of (3.18) from Lemmas 3.10 and 3.12, the Kolmogorov inclusions (3.24), and the deterministic stability estimate of Lemma 2.2. Each is a separate formal theorem. The authors' proofs follow their statements in Sections 2 and 3 and the Appendix; Lemma 3.5 is quoted from Lê, Theorem 2.1 without proof in the article.

What the result provides

The target gives convergence of entire solution paths in a Hölder topology, in probability. In the paper, this lets a model-specific estimate of the negative Hölder distance be used to conclude convergence of the corresponding dynamical system. Theorem A is obtained by adding such an estimate for a stationary mixing fast process; the deterministic statement here isolates the part independent of that model. These uses are described in Sections 1 and 3.

The published paper proves the result. This mission asks for machine-checked proofs of the stated Lean theorems. Its reusable work would include a pathwise Young integral over arbitrary partitions, Hölder norms that retain the value +∞+\infty+∞, a two-sided vector fractional Brownian motion interface, and stochastic estimates for time-dependent integrands. The goal and milestones in this proposal are statements with proof placeholders; the local build checks their types, not their proofs.

Main obstruction

Pointwise control of fn−fˉf_n-\bar ffn​−fˉ​ is unavailable: a rapidly oscillating function can be large at most times and still have small time antiderivatives. Applying a pathwise Young stability estimate directly to fn−fˉf_n-\bar ffn​−fˉ​ therefore leaves a coefficient norm that need not tend to zero. The integral estimate must connect cancellation measured by ∣⋅∣−κ,γ|\cdot|_{-\kappa,\gamma}∣⋅∣−κ,γ​ with LpL^pLp bounds on the driven equation. The dependence on the past of the fractional Brownian motion makes this step more delicate than an estimate for a deterministic integrand. The required estimates and their domains appear in Lemmas 3.4–3.12.

Formalization scope

Time is R\mathbb RR for the two-sided driver; coefficient and solution assertions use [0,T][0,T][0,T] with T>0T>0T>0. State spaces are finite-dimensional real Euclidean spaces. The goal takes H∈(1/2,1)H\in(1/2,1)H∈(1/2,1) from the standing setting of Section 3.1, and it assumes the locally Hölder version of BBB described on p. 1826. The solutions are required almost surely to satisfy the Young equation from the same x0x_0x0​. This requirement fixes the processes to the equations; the convergence claim cannot be obtained by choosing arbitrary paths. Probability of the norm event uses the measure of that set, and the conclusion is convergence in probability in the full CαC^\alphaCα norm.

The Lean definitions use [0,∞][0,\infty][0,∞] for path and process seminorms, so unbounded suprema are not replaced by a default real value. The negative Hölder hypothesis uses real bounds Kn→0K_n\to0Kn​→0 instead of an infimum that could be meaningless on an empty class. The Young integral quantifies over every monotone partition and the mesh limit of its left-point sums. The fBm definition includes Gaussianity, means, and every cross-coordinate covariance; its path version is explicit. Interval integrals in the coefficient bound are applied to time-continuous coefficients in the goal. Lemmas 3.4 and 3.10 use the increment filtration, the case needed for this theorem.

The paper's Lemma 2.2 needs a common initial value for its homogeneous-seminorm estimate; the formal statement includes it and records the printed issue. The paper writes the limiting coefficient with a time argument but uses an autonomous coefficient in (3.23), which the goal adopts. This proposal covers the dBdBdB equation. The final dtdtdt extension in Theorem 3.13 and the conditional-expectation estimate (3.19) remain outside its current items. Contributions that close these gaps, or supply a general conditional Wiener integral matching the article's notation, would extend the development.

Selected references

  • M. Hairer and X.-M. Li, Averaging dynamics driven by fractional Brownian motion, Annals of Probability 48(4), 2020, 1826–1860. DOI.
  • K. Lê, A stochastic sewing lemma and applications, 2018, arXiv:1810.10500, cited as [23] in Hairer and Li. arXiv.
9 thms1 active userReviewed
Machine LearningStatistics·Captain: mikedeng1

Certifying Some Distributional Robustness with Principled Adversarial Training 2: The Empirical Robustness Level E_P̂ₙ[c(T_γ(θ;Z),Z)] Concentrates Uniformly over ΘResearch Paper

Motivation

Neural networks trained by empirical risk minimization are sensitive to small adversarial perturbations of their inputs. Distributionally robust optimization addresses this by training against the worst distribution in a neighbourhood of the data distribution. Sinha, Namkoong, Volpi and Duchi (arXiv:1710.10571v5, ICLR 2018) take the neighbourhood to be a Wasserstein ball and relax the hard constraint to a Lagrangian penalty with multiplier γ\gammaγ. For large γ\gammaγ the inner maximization becomes strongly concave and is solved per data point by a transportation map TγT_\gammaTγ​, which turns each training example into its worst-case perturbation.

The average transport cost ρ^n(θ)=EP^n[c(Tγ(θ;Z),Z)]\widehat\rho_n(\theta)=\mathbb E_{\widehat P_n}[c(T_\gamma(\theta;Z),Z)]ρ​n​(θ)=EPn​​[c(Tγ​(θ;Z),Z)] is the amount of robustness the trained model actually certifies on the sample: by the paper's duality result, the empirical worst-case loss at level ρ^n(θ)\widehat\rho_n(\theta)ρ​n​(θ) is computable, and it is the certificate of robustness that the paper advertises. A certificate measured on the training sample is only useful if it transfers to the population. This mission formalizes the result that it does: ρ^n(θ)\widehat\rho_n(\theta)ρ​n​(θ) concentrates around its population counterpart, uniformly over the parameter set.

Setting

Parameters θ\thetaθ range over a set Θ⊆Rd\Theta\subseteq\mathbb R^dΘ⊆Rd and data over a set Z⊆{z∈Rm:∥z∥2≤Mz}Z\subseteq\{z\in\mathbb R^m:\|z\|_2\le M_z\}Z⊆{z∈Rm:∥z∥2​≤Mz​}, both with the Euclidean norm. A transportation cost c:Z×Z→[0,∞)c:Z\times Z\to[0,\infty)c:Z×Z→[0,∞) satisfies c(z,z)=0c(z,z)=0c(z,z)=0 and Assumption A: ccc is continuous and z↦c(z,z0)z\mapsto c(z,z_0)z↦c(z,z0​) is 111-strongly convex for each z0∈Zz_0\in Zz0​∈Z. A loss ℓ(θ;z)\ell(\theta;z)ℓ(θ;z) satisfies Assumption B: its partial gradients ∇θℓ\nabla_\theta\ell∇θ​ℓ, ∇zℓ\nabla_z\ell∇z​ℓ are Lipschitz in each variable, with constants Lθθ,Lzz,Lθz,LzθL_{\theta\theta},L_{zz},L_{\theta z},L_{z\theta}Lθθ​,Lzz​,Lθz​,Lzθ​.

For γ>Lzz\gamma>L_{zz}γ>Lzz​ the problem z↦ℓ(θ;z)−γc(z,z0)z\mapsto\ell(\theta;z)-\gamma c(z,z_0)z↦ℓ(θ;z)−γc(z,z0​) is (γ−Lzz)(\gamma-L_{zz})(γ−Lzz​)-strongly concave on ZZZ, and the transportation map is its maximizer,

Tγ(θ;z0)∈argmax⁡z∈Z{ℓ(θ;z)−γc(z,z0)}.T_\gamma(\theta;z_0)\in\operatorname*{argmax}_{z\in Z}\{\ell(\theta;z)-\gamma c(z,z_0)\}.Tγ​(θ;z0​)∈z∈Zargmax​{ℓ(θ;z)−γc(z,z0​)}.

Given i.i.d. samples Z1,…,Zn∼P0Z_1,\dots,Z_n\sim P_0Z1​,…,Zn​∼P0​ with P0(Z)=1P_0(Z)=1P0​(Z)=1 and empirical distribution P^n\widehat P_nPn​, the empirical robustness level and its population counterpart are

ρ^n(θ)=1n∑i=1nc(Tγ(θ;Zi),Zi),ρ(θ)=EP0[c(Tγ(θ;Z),Z)].\widehat\rho_n(\theta)=\frac1n\sum_{i=1}^n c(T_\gamma(\theta;Z_i),Z_i),\qquad \rho(\theta)=\mathbb E_{P_0}[c(T_\gamma(\theta;Z),Z)].ρ​n​(θ)=n1​i=1∑n​c(Tγ​(θ;Zi​),Zi​),ρ(θ)=EP0​​[c(Tγ​(θ;Z),Z)].

The complexity of Θ\ThetaΘ enters through its covering number N(Θ,ϵ,∥⋅∥)N(\Theta,\epsilon,\|\cdot\|)N(Θ,ϵ,∥⋅∥), the least number of points θ1,…,θN∈Θ\theta_1,\dots,\theta_N\in\Thetaθ1​,…,θN​∈Θ such that every θ∈Θ\theta\in\Thetaθ∈Θ is within ϵ\epsilonϵ of some θi\theta_iθi​ (+∞+\infty+∞ if there is no finite such family).

Two alternative regularity conditions close the setting, with a constant Lc>0L_c>0Lc​>0: (i) ccc is LcL_cLc​-Lipschitz on ZZZ in each argument; or (ii) ℓ(θ;z)∈[0,Mℓ]\ell(\theta;z)\in[0,M_\ell]ℓ(θ;z)∈[0,Mℓ​] and z↦ℓ(θ;z)z\mapsto\ell(\theta;z)z↦ℓ(θ;z) is γLc\gamma L_cγLc​-Lipschitz. Put K=2LcMzK=2L_cM_zK=2Lc​Mz​ in case (i) and K=Mℓ/γK=M_\ell/\gammaK=Mℓ​/γ in case (ii).

Formalization targets

Goal: the uniform tail bound of §B.7

For every t>0t>0t>0,

P(sup⁡θ∈Θ∣ρ^n(θ)−ρ(θ)∣≥t)≤2 N(Θ,(γ−Lzz)t4LcLzθ,∥⋅∥2)exp⁡(−nt28K2).\mathbb P\Bigl(\sup_{\theta\in\Theta}|\widehat\rho_n(\theta)-\rho(\theta)|\ge t\Bigr)\le 2\,N\Bigl(\Theta,\frac{(\gamma-L_{zz})t}{4L_cL_{z\theta}},\|\cdot\|_2\Bigr)\exp\Bigl(-\frac{nt^2}{8K^2}\Bigr).P(θ∈Θsup​∣ρ​n​(θ)−ρ(θ)∣≥t)≤2N(Θ,4Lc​Lzθ​(γ−Lzz​)t​,∥⋅∥2​)exp(−8K2nt2​).

In case (i), 8K2=32Lc2Mz28K^2=32L_c^2M_z^28K2=32Lc2​Mz2​ and this is exactly the last display of the proof of Theorem 4 (p. 44); case (ii) is the "identical reasoning" the paper invokes with K=Mℓ/γK=M_\ell/\gammaK=Mℓ​/γ. Θ\ThetaΘ is an arbitrary subset of Rd\mathbb R^dRd: no compactness is assumed, and an infinite covering number makes the bound vacuous rather than false.

Milestones

In the order the proof uses them:

  1. (4): strong concavity of z↦ℓ(θ;z)−γc(z,z0)z\mapsto\ell(\theta;z)-\gamma c(z,z_0)z↦ℓ(θ;z)−γc(z,z0​) with modulus γ−Lzz\gamma-L_{zz}γ−Lzz​;
  2. (16): ∥Tγ(θ1;z)−Tγ(θ2;z)∥≤Lzθγ−Lzz∥θ1−θ2∥\|T_\gamma(\theta_1;z)-T_\gamma(\theta_2;z)\|\le\frac{L_{z\theta}}{\gamma-L_{zz}}\|\theta_1-\theta_2\|∥Tγ​(θ1​;z)−Tγ​(θ2​;z)∥≤γ−Lzz​Lzθ​​∥θ1​−θ2​∥;
  3. Lemma 5: the transported cost is LcLzθγ−Lzz\frac{L_cL_{z\theta}}{\gamma-L_{zz}}γ−Lzz​Lc​Lzθ​​-Lipschitz in θ\thetaθ;
  4. the deviation ∣ρ^n(θ)−ρ(θ)∣|\widehat\rho_n(\theta)-\rho(\theta)|∣ρ​n​(θ)−ρ(θ)∣ is 2LcLzθγ−Lzz\frac{2L_cL_{z\theta}}{\gamma-L_{zz}}γ−Lzz​2Lc​Lzθ​​-Lipschitz on Θ\ThetaΘ;
  5. a (γ−Lzz)t4LcLzθ\frac{(\gamma-L_{zz})t}{4L_cL_{z\theta}}4Lc​Lzθ​(γ−Lzz​)t​-cover transfers a bound t/2t/2t/2 at its centres to a bound ttt on the supremum;
  6. 0≤c(Tγ(θ;z),z)≤K0\le c(T_\gamma(\theta;z),z)\le K0≤c(Tγ​(θ;z),z)≤K;
  7. Hoeffding's inequality at a fixed θ\thetaθ: P(∣ρ^n(θ)−ρ(θ)∣≥t/2)≤2e−nt2/(8K2)\mathbb P(|\widehat\rho_n(\theta)-\rho(\theta)|\ge t/2)\le 2e^{-nt^2/(8K^2)}P(∣ρ​n​(θ)−ρ(θ)∣≥t/2)≤2e−nt2/(8K2).

Significance

The tail bound gives a d/n\sqrt{d/n}d/n​-type generalization rate for the gap between the robustness certified on the sample and the robustness achieved on the population: with log⁡N(Θ,ϵ)≤dlog⁡(1+diam⁡(Θ)/ϵ)\log N(\Theta,\epsilon)\le d\log(1+\operatorname{diam}(\Theta)/\epsilon)logN(Θ,ϵ)≤dlog(1+diam(Θ)/ϵ) for bounded Θ\ThetaΘ, the deviation is of order K(dlog⁡(1/t)+log⁡(1/δ))/nK\sqrt{(d\log(1/t)+\log(1/\delta))/n}K(dlog(1/t)+log(1/δ))/n​. It also implies, as the paper notes, that the certified level ρ^n\widehat\rho_nρ​n​ stays positive for large nnn whenever the loss is not invariant to the data.

The paper states the result as Theorem 4, display (17), in "with probability 1−e−t1-e^{-t}1−e−t" form with the constant 4B4B4B. That form does not follow from the proved tail bound as printed: inverting the tail bound gives a deviation 42 B(t+log⁡2N)/n4\sqrt2\,B\sqrt{(t+\log 2N)/n}42​B(t+log2N)/n​ with a covering radius that depends on the deviation rather than on ttt, whereas (17) uses ttt in both places and drops 2\sqrt22​ and log⁡2\log 2log2. The mission therefore poses what the proof establishes, the tail bound above, and does not claim (17). None of this paper's results has a machine-checked proof; Hoeffding's inequality and covering numbers are in Mathlib, the argmax stability and the cover argument are not.

Difficulty

Concentration at a single fixed θ\thetaθ is Hoeffding's inequality and carries no information about the trained model, whose parameter depends on the sample. The supremum over an arbitrary, possibly uncountable Θ\ThetaΘ is the content. A union bound needs a finite family, and transferring control from a finite family to all of Θ\ThetaΘ needs the deviation to be Lipschitz in θ\thetaθ. The deviation is built from the transportation map, which is an argmax and has no closed form; its regularity in θ\thetaθ comes only from strong concavity, which holds only for γ>Lzz\gamma>L_{zz}γ>Lzz​. In case (ii) the cost ccc itself has no Lipschitz assumption, and the Lipschitz property of c(Tγ(θ;z),z)c(T_\gamma(\theta;z),z)c(Tγ​(θ;z),z) has to be extracted from the optimality of TγT_\gammaTγ​ at two different parameters. Finally, the supremum event is not obviously measurable.

Formalization scope

Spaces are EuclideanSpace ℝ (Fin d) and EuclideanSpace ℝ (Fin m); the paper's general norm is specialised to ℓ2\ell_2ℓ2​, in Assumptions A and B, in cases (i) and (ii), and in the covering number. The cost is real valued and constrained on Z×ZZ\times ZZ×Z only; the loss is typed on Rd×Rm\mathbb R^d\times\mathbb R^mRd×Rm, but Assumption B and the argmax property of the transportation map are required only for θ∈Θ\theta\in\Thetaθ∈Θ and z∈Zz\in Zz∈Z, matching the paper's ℓ:Θ×Z→R\ell:\Theta\times Z\to\mathbb Rℓ:Θ×Z→R. The transportation map is a function together with its argmax property. Samples are the coordinates of the product measure P0⊗nP_0^{\otimes n}P0⊗n​; P\mathbb PP of the supremum event is its outer measure, so no measurability of the event is assumed. The supremum is a real supremum of quantities bounded by 2K2K2K. The covering number is Mathlib's Metric.coveringNumber (closed balls, centres in Θ\ThetaΘ, values in N∪{∞}\mathbb N\cup\{\infty\}N∪{∞}). The definition on p. 9 does not say where the centres lie; the proof applies Hoeffding's inequality "for any fixed θ∈Θ\theta\in\Thetaθ∈Θ" at the centres, and the loss is only defined on Θ\ThetaΘ, so the centres are taken in Θ\ThetaΘ.

Hypotheses made explicit that the page leaves implicit:

  • γ>Lzz\gamma>L_{zz}γ>Lzz​ (p. 43, "our strong concavity assumption"), so [γ−Lzz]+=γ−Lzz[\gamma-L_{zz}]_+=\gamma-L_{zz}[γ−Lzz​]+​=γ−Lzz​;
  • Lc>0L_c>0Lc​>0 and Lzθ>0L_{z\theta}>0Lzθ​>0, since the covering radius divides by 4LcLzθ4L_cL_{z\theta}4Lc​Lzθ​;
  • ZZZ convex (needed for the transportation map to be stable);
  • z↦c(Tγ(θ;z),z)z\mapsto c(T_\gamma(\theta;z),z)z↦c(Tγ​(θ;z),z) is P0P_0P0​-measurable for every θ∈Θ\theta\in\Thetaθ∈Θ (the standing convention).

Case (ii) is assumed for θ∈Θ\theta\in\Thetaθ∈Θ only, as printed in Theorem 4.

A bound for one fixed θ\thetaθ, a supremum taken outside the probability, a covering radius not tied to ttt, or a covering number of a set other than Θ\ThetaΘ would be a different and much weaker statement; the goal has the supremum over Θ\ThetaΘ inside the event and the radius (γ−Lzz)t/(4LcLzθ)(\gamma-L_{zz})t/(4L_cL_{z\theta})(γ−Lzz​)t/(4Lc​Lzθ​).

A complete development needs the stability of maximizers of strongly concave functions under Lipschitz perturbation of the gradient, a finite-cover union bound for outer measures, and Hoeffding's inequality for bounded i.i.d. coordinates of a product measure. The first two are reusable well beyond this mission. Proofs of any milestone are welcome, as are alternative arguments for the goal.

Selected references

  • A. Sinha, H. Namkoong, R. Volpi, J. Duchi, Certifying Some Distributional Robustness with Principled Adversarial Training, ICLR 2018; arXiv:1710.10571v5, 2020. https://arxiv.org/abs/1710.10571v5
  • W. Hoeffding, Probability inequalities for sums of bounded random variables, Journal of the American Statistical Association 58 (1963), 13–30. https://doi.org/10.1080/01621459.1963.10500830
  • A. W. van der Vaart, J. A. Wellner, Weak Convergence and Empirical Processes, Springer, 1996. https://doi.org/10.1007/978-1-4757-2545-2
  • C. Villani, Optimal Transport: Old and New, Springer, 2009. https://doi.org/10.1007/978-3-540-71050-9
10 thms1 active userReviewed
Machine LearningOptimization·Captain: mikedeng1

Certifying Some Distributional Robustness with Principled Adversarial Training 1: SGD on the Wasserstein Penalty Surrogate Finds an Approximate Stationary Point at Rate 4σ√(L_φΔ_F/T)Research Paper

Motivation

Neural networks trained by empirical risk minimization are sensitive to small, deliberately chosen perturbations of their inputs ("adversarial examples"). A standard defence is adversarial training: at each step, perturb the training example to increase the loss, then take a gradient step at the perturbed point. Sinha, Namkoong, Volpi and Duchi (arXiv:1710.10571, ICLR 2018) recast adversarial training as distributionally robust optimization over a Wasserstein ball around the data distribution, and replaced the hard ball constraint by a Lagrangian penalty. The resulting procedure perturbs each example by solving a penalized inner maximization, and the paper shows that, when the penalty is large enough relative to the smoothness of the loss, this inner problem is strongly concave and therefore tractable.

This mission formalizes the paper's computational guarantee, Theorem 2: stochastic gradient descent on the penalized robust objective, with an inexact adversary, reaches an approximate stationary point at the rate of smooth nonconvex stochastic optimization. The analysis follows Ghadimi and Lan's argument for nonconvex SGD (arXiv:1309.5549), adapted to a gradient oracle that is biased by the inexact inner maximization.

Setting

Parameters θ\thetaθ live in Θ=Rd\Theta=\mathbb R^dΘ=Rd and data in a set Z⊆RmZ\subseteq\mathbb R^mZ⊆Rm; both spaces carry the Euclidean norm ∥⋅∥2\|\cdot\|_2∥⋅∥2​. Data are drawn from a probability distribution P0P_0P0​ on ZZZ.

A transportation cost c:Z×Z→[0,∞)c:Z\times Z\to[0,\infty)c:Z×Z→[0,∞) measures how expensive it is to move a data point z0z_0z0​ to zzz. Assumption A: ccc is continuous, c(z,z)=0c(z,z)=0c(z,z)=0, and for each z0∈Zz_0\in Zz0​∈Z the map z↦c(z,z0)z\mapsto c(z,z_0)z↦c(z,z0​) is 111-strongly convex. The canonical example is c(z,z0)=12∥z−z0∥22c(z,z_0)=\tfrac12\|z-z_0\|_2^2c(z,z0​)=21​∥z−z0​∥22​.

A loss ℓ(θ;z)\ell(\theta;z)ℓ(θ;z) has partial gradients ∇θℓ\nabla_\theta\ell∇θ​ℓ and ∇zℓ\nabla_z\ell∇z​ℓ. Assumption B: these gradients are Lipschitz in each argument,

∥∇θℓ(θ;z)−∇θℓ(θ′;z)∥≤Lθθ∥θ−θ′∥,∥∇zℓ(θ;z)−∇zℓ(θ;z′)∥≤Lzz∥z−z′∥,\|\nabla_\theta\ell(\theta;z)-\nabla_\theta\ell(\theta';z)\|\le L_{\theta\theta}\|\theta-\theta'\|,\quad \|\nabla_z\ell(\theta;z)-\nabla_z\ell(\theta;z')\|\le L_{zz}\|z-z'\|,∥∇θ​ℓ(θ;z)−∇θ​ℓ(θ′;z)∥≤Lθθ​∥θ−θ′∥,∥∇z​ℓ(θ;z)−∇z​ℓ(θ;z′)∥≤Lzz​∥z−z′∥, ∥∇θℓ(θ;z)−∇θℓ(θ;z′)∥≤Lθz∥z−z′∥,∥∇zℓ(θ;z)−∇zℓ(θ′;z)∥≤Lzθ∥θ−θ′∥.\|\nabla_\theta\ell(\theta;z)-\nabla_\theta\ell(\theta;z')\|\le L_{\theta z}\|z-z'\|,\quad \|\nabla_z\ell(\theta;z)-\nabla_z\ell(\theta';z)\|\le L_{z\theta}\|\theta-\theta'\|.∥∇θ​ℓ(θ;z)−∇θ​ℓ(θ;z′)∥≤Lθz​∥z−z′∥,∥∇z​ℓ(θ;z)−∇z​ℓ(θ′;z)∥≤Lzθ​∥θ−θ′∥.

For a penalty γ>0\gamma>0γ>0, the robust surrogate and the robust objective are

φγ(θ;z0)=sup⁡z∈Z{ℓ(θ;z)−γc(z,z0)},F(θ)=EP0[φγ(θ;Z)].\varphi_\gamma(\theta;z_0)=\sup_{z\in Z}\{\ell(\theta;z)-\gamma c(z,z_0)\},\qquad F(\theta)=\mathbb E_{P_0}[\varphi_\gamma(\theta;Z)].φγ​(θ;z0​)=z∈Zsup​{ℓ(θ;z)−γc(z,z0​)},F(θ)=EP0​​[φγ​(θ;Z)].

By Wasserstein duality, F(θ)F(\theta)F(θ) equals the worst-case expected loss over all distributions PPP, penalized by γ\gammaγ times the transport distance from PPP to P0P_0P0​.

Algorithm 1 (with Θ=Rd\Theta=\mathbb R^dΘ=Rd and constant stepsize α\alphaα): from θ0\theta^0θ0, for t=0,…,T−1t=0,\dots,T-1t=0,…,T−1, draw zt∼P0z^t\sim P_0zt∼P0​ independently, find an ϵ\epsilonϵ-approximate maximizer z^t∈Z\hat z^t\in Zz^t∈Z, meaning ℓ(θt;z^t)−γc(z^t,zt)≥φγ(θt;zt)−ϵ\ell(\theta^t;\hat z^t)-\gamma c(\hat z^t,z^t)\ge\varphi_\gamma(\theta^t;z^t)-\epsilonℓ(θt;z^t)−γc(z^t,zt)≥φγ​(θt;zt)−ϵ, and set θt+1=θt−α∇θℓ(θt;z^t)\theta^{t+1}=\theta^t-\alpha\nabla_\theta\ell(\theta^t;\hat z^t)θt+1=θt−α∇θ​ℓ(θt;z^t).

Formalization targets

Goal: Theorem 2 (p. 6)

Let γ>Lzz\gamma>L_{zz}γ>Lzz​, Lφ=Lθθ+LθzLzθ/(γ−Lzz)L_\varphi=L_{\theta\theta}+L_{\theta z}L_{z\theta}/(\gamma-L_{zz})Lφ​=Lθθ​+Lθz​Lzθ​/(γ−Lzz​), ΔF≥F(θ0)−inf⁡θF(θ)\Delta_F\ge F(\theta^0)-\inf_\theta F(\theta)ΔF​≥F(θ0)−infθ​F(θ), and assume E∥∇F(θ)−∇θφγ(θ;Z)∥22≤σ2\mathbb E\|\nabla F(\theta)-\nabla_\theta\varphi_\gamma(\theta;Z)\|_2^2\le\sigma^2E∥∇F(θ)−∇θ​φγ​(θ;Z)∥22​≤σ2 for all θ\thetaθ. With α=ΔF/(LφTσ2)\alpha=\sqrt{\Delta_F/(L_\varphi T\sigma^2)}α=ΔF​/(Lφ​Tσ2)​ and T≥LφΔF/σ2T\ge L_\varphi\Delta_F/\sigma^2T≥Lφ​ΔF​/σ2,

1T∑t=0T−1E[∥∇F(θt)∥22]−4Lθz2γ−Lzz ϵ≤4σLφΔFT.\frac1T\sum_{t=0}^{T-1}\mathbb E\big[\|\nabla F(\theta^t)\|_2^2\big]-\frac{4L_{\theta z}^2}{\gamma-L_{zz}}\,\epsilon\le4\sigma\sqrt{\frac{L_\varphi\Delta_F}{T}}.T1​t=0∑T−1​E[∥∇F(θt)∥22​]−γ−Lzz​4Lθz2​​ϵ≤4σTLφ​ΔF​​​.

Milestones

  1. (4), p. 4: for γ≥Lzz\gamma\ge L_{zz}γ≥Lzz​, z↦ℓ(θ;z)−γc(z,z0)z\mapsto\ell(\theta;z)-\gamma c(z,z_0)z↦ℓ(θ;z)−γc(z,z0​) is (γ−Lzz)(\gamma-L_{zz})(γ−Lzz​)-strongly concave.
  2. Lemma 1, p. 5: the supremum fˉ(θ)=sup⁡z∈Zf(θ,z)\bar f(\theta)=\sup_{z\in Z}f(\theta,z)fˉ​(θ)=supz∈Z​f(θ,z) of a λ\lambdaλ-strongly concave family satisfying Assumption B is differentiable, with ∇fˉ(θ)=∇θf(θ,z⋆(θ))\nabla\bar f(\theta)=\nabla_\theta f(\theta,z_\star(\theta))∇fˉ​(θ)=∇θ​f(θ,z⋆​(θ)), a (Lzθ/λ)(L_{z\theta}/\lambda)(Lzθ​/λ)-Lipschitz maximizer and a (Lθθ+LθzLzθ/λ)(L_{\theta\theta}+L_{\theta z}L_{z\theta}/\lambda)(Lθθ​+Lθz​Lzθ​/λ)-Lipschitz gradient.
  3. §2.1, p. 6: φγ(⋅;z0)\varphi_\gamma(\cdot;z_0)φγ​(⋅;z0​) has LφL_\varphiLφ​-Lipschitz gradient ∇θφγ(θ;z0)=∇θℓ(θ;z⋆(z0,θ))\nabla_\theta\varphi_\gamma(\theta;z_0)=\nabla_\theta\ell(\theta;z_\star(z_0,\theta))∇θ​φγ​(θ;z0​)=∇θ​ℓ(θ;z⋆​(z0​,θ)).
  4. (31) and (32), p. 40: the descent inequality for an LφL_\varphiLφ​-smooth objective, and its form with the bias δ=g−∇θφγ\delta=g-\nabla_\theta\varphi_\gammaδ=g−∇θ​φγ​ separated, for α≤1/Lφ\alpha\le1/L_\varphiα≤1/Lφ​.
  5. The bias bound, p. 41: ∥δ∥22≤2Lθz2ϵ/(γ−Lzz)\|\delta\|_2^2\le 2L_{\theta z}^2\epsilon/(\gamma-L_{zz})∥δ∥22​≤2Lθz2​ϵ/(γ−Lzz​).
  6. The one-step bound in conditional expectation and the averaged bound, p. 41:
1T∑t<TE∥∇F(θt)∥22≤2ϵ^+2ΔFαT+2Lφασ2,ϵ^=2Lθz2γ−Lzzϵ,\frac1T\sum_{t<T}\mathbb E\|\nabla F(\theta^t)\|_2^2\le2\hat\epsilon+\frac{2\Delta_F}{\alpha T}+2L_\varphi\alpha\sigma^2,\qquad\hat\epsilon=\frac{2L_{\theta z}^2}{\gamma-L_{zz}}\epsilon,T1​t<T∑​E∥∇F(θt)∥22​≤2ϵ^+αT2ΔF​​+2Lφ​ασ2,ϵ^=γ−Lzz​2Lθz2​​ϵ,

for every constant stepsize 0<α≤1/Lφ0<\alpha\le1/L_\varphi0<α≤1/Lφ​. The goal is this bound at the balancing stepsize.

Significance

The result. Theorem 2 certifies that the robust training procedure is no harder, in iteration count, than ordinary nonconvex SGD: the stationarity measure decays as T−1/2T^{-1/2}T−1/2, and the inexactness of the adversary enters only as an additive term proportional to ϵ\epsilonϵ, independent of TTT. This is what licenses running a few steps of gradient ascent for the inner problem rather than solving it exactly. It also explains the role of the penalty: the guarantee requires γ>Lzz\gamma>L_{zz}γ>Lzz​, so smoothness of the loss in the data (smooth activations rather than ReLUs) is what makes the robust problem tractable.

Formalizing it. The result and its proof are published; nothing here is open. To our knowledge none of it is machine-checked. A formal proof would add a verified Danskin-type theorem for strongly concave suprema with explicit Lipschitz constants (Lemma 1), which is reusable well beyond this paper, and a verified nonconvex SGD analysis with a biased oracle, including the measure-theoretic bookkeeping (independence of the fresh sample from the current iterate, differentiation under the integral) that the paper's proof compresses into one line.

Difficulty

The deterministic inequalities (31), (32) and the bias bound are short. The weight sits elsewhere. First, Lemma 1 requires differentiability of a supremum over a possibly unbounded set, with the maximizer depending on θ\thetaθ; the paper invokes a Danskin-type theorem that is not in Mathlib, and the cost ccc is not assumed differentiable, so the first-order optimality conditions for the inner problem must be derived from strong convexity alone. Second, the expectation step uses E[∇θφγ(θt;zt)∣θt]=∇F(θt)\mathbb E[\nabla_\theta\varphi_\gamma(\theta^t;z^t)\mid\theta^t]=\nabla F(\theta^t)E[∇θ​φγ​(θt;zt)∣θt]=∇F(θt), which needs both differentiation under the integral sign for FFF and the independence of ztz^tzt from θt\theta^tθt along the run. Third, finiteness of every expectation along the run is not assumed and must be propagated from step to step.

Formalization scope

Lean representation: Θ\ThetaΘ and Rm\mathbb R^mRm are EuclideanSpace ℝ (Fin d) and EuclideanSpace ℝ (Fin m); ZZZ is a set in the latter. The loss is a function on all of Rd×Rm\mathbb R^d\times\mathbb R^mRd×Rm; its partial gradients are required at points of ZZZ and the Lipschitz conditions on ZZZ only. Strong convexity is Mathlib's StrongConvexOn, whose modulus convention matches the paper. Suprema are sSup over the image of ZZZ. FFF is a Bochner integral, and φγ(θ;⋅)\varphi_\gamma(\theta;\cdot)φγ​(θ;⋅) is assumed P0P_0P0​-integrable for every θ\thetaθ. The run is driven by (z0,…,zT−1)∼P0⊗T(z^0,\dots,z^{T-1})\sim P_0^{\otimes T}(z0,…,zT−1)∼P0⊗T​ (Measure.pi). The expectations in the conclusion and in the variance hypothesis are lower Lebesgue integrals in [0,∞][0,\infty][0,∞]; the ϵ\epsilonϵ-term of Theorem 2 is moved to the right-hand side.

Hypotheses made explicit that the paper leaves implicit: γ>Lzz\gamma>L_{zz}γ>Lzz​ (the paper writes γ≥Lzz\gamma\ge L_{zz}γ≥Lzz​ but divides by γ−Lzz\gamma-L_{zz}γ−Lzz​); ZZZ nonempty, convex and closed; Lipschitz constants nonnegative; σ>0\sigma>0σ>0, ΔF>0\Delta_F>0ΔF​>0, T≥1T\ge1T≥1, ϵ≥0\epsilon\ge0ϵ≥0; FFF bounded below (so inf⁡θF\inf_\theta Finfθ​F is a real number); integrability of φγ(θ;⋅)\varphi_\gamma(\theta;\cdot)φγ​(θ;⋅); P0(Z)=1P_0(Z)=1P0​(Z)=1; joint measurability of the oracle (θ,z0)↦z^(θ,z0)(\theta,z_0)\mapsto\hat z(\theta,z_0)(θ,z0​)↦z^(θ,z0​); and, in (32) and the later milestones, 0<α≤1/Lφ0<\alpha\le1/L_\varphi0<α≤1/Lφ​, which the paper assumes "in the rest of the proof". Lemma 1 is specialized to the Euclidean norm, the case Theorem 2 uses. The cost ccc is not assumed differentiable.

A trivializing formalization is ruled out: the run uses an oracle tied to φγ\varphi_\gammaφγ​ by the ϵ\epsilonϵ-gap and to ZZZ by membership, φγ\varphi_\gammaφγ​ is the supremum of the penalized loss over ZZZ, FFF is its P0P_0P0​-integral, and the hypotheses are satisfiable (for instance ℓ(θ;z)=12∥θ∥2\ell(\theta;z)=\tfrac12\|\theta\|^2ℓ(θ;z)=21​∥θ∥2, c(z,z0)=12∥z−z0∥2c(z,z_0)=\tfrac12\|z-z_0\|^2c(z,z0​)=21​∥z−z0​∥2, Z=RmZ=\mathbb R^mZ=Rm, P0P_0P0​ a point mass, ϵ=0\epsilon=0ϵ=0). Expectations as lower integrals prevent a non-integrable ∥∇F(θt)∥2\|\nabla F(\theta^t)\|^2∥∇F(θt)∥2 from being read as 000.

Infrastructure a complete development needs: the descent lemma for functions with Lipschitz gradient; first-order characterizations of strong concavity on a convex set; a Danskin/envelope theorem for strongly concave suprema; differentiation under the integral for FFF; and conditional-expectation bookkeeping for Markov-type recursions driven by i.i.d. samples on a finite product space. The first three are reusable in any minimax or bilevel analysis. Contributions to any milestone, or to these supporting facts as separate lemmas, are welcome.

Selected references

  • A. Sinha, H. Namkoong, R. Volpi, J. Duchi, Certifying Some Distributional Robustness with Principled Adversarial Training, ICLR 2018; arXiv:1710.10571v5, 2020. https://arxiv.org/abs/1710.10571
  • S. Ghadimi, G. Lan, Stochastic First- and Zeroth-Order Methods for Nonconvex Stochastic Programming, SIAM Journal on Optimization 23(4), 2013. https://doi.org/10.1137/120880811
  • J. F. Bonnans, A. Shapiro, Perturbation Analysis of Optimization Problems, Springer, 2000. https://doi.org/10.1007/978-1-4612-1394-9
  • R. Gao, A. Kleywegt, Distributionally Robust Stochastic Optimization with Wasserstein Distance, Mathematics of Operations Research 48(2), 2023; arXiv:1604.02199. https://arxiv.org/abs/1604.02199
10 thms1 active userReviewed
Dynamic ProgrammingMarkov ChainOperations Research·Captain: mikedeng1

Partially Observable Total-Cost Markov Decision Processes with Weakly Continuous Transition Probabilities 2: Setwise-Continuous Observations Do Not Make the Belief Transition Weakly ContinuousResearch Paper

Motivation

A partially observable Markov decision process describes a controlled system whose state is hidden but whose observations are available to a decision maker. The distribution of the hidden state conditional on the observations is the belief state. Replacing the hidden state by this distribution turns the model into a fully observed decision process, provided the resulting belief transition has the regularity needed by the usual existence and dynamic programming results. Feinberg, Kasyanov, and Zgurovsky give sufficient conditions for that regularity in their study of total-cost POMDPs, Theorems 3.6 and 3.7.

Their observation-kernel condition is continuity in total variation. Example 4.1 of the same paper tests its sharpness: the authors construct a two-state model with setwise-continuous observations for which the belief transition fails to be weakly continuous. This matters when selecting assumptions for a POMDP theorem. Setwise convergence already requires convergence of the probabilities of every fixed observation event, so it can look adequate for belief updates. The example shows that it is still too weak for the stated conclusion.

Setting

The hidden state space is X={1,2}X=\{1,2\}X={1,2} and the observation space is Y=[0,1]Y=[0,1]Y=[0,1] with Lebesgue probability measure mmm. The action space is A={0}∪{1/n:n=1,2,…}A=\{0\}\cup\{1/n:n=1,2,\ldots\}A={0}∪{1/n:n=1,2,…}, equipped with the topology inherited from R\mathbb RR. The state does not move: its transition law P(⋅∣x,a)P(\cdot\mid x,a)P(⋅∣x,a) is the point mass at xxx for every action. Its continuity in total variation is immediate from this fixed-state behavior.

At hidden state 1 the observation law is always mmm. At hidden state 2, action zero also produces mmm. Under action 1/n1/n1/n, state 2 produces a law m(n)=f(n)mm^{(n)}=f^{(n)}mm(n)=f(n)m. The density f(n)f^{(n)}f(n) equals zero on the open left half of each cell in the 2n2^n2n-part dyadic partition and equals two elsewhere, including the cell endpoints. Equation (4.1) gives the exact intervals. These laws oscillate more rapidly as nnn grows, and the paper establishes m(n)(C)→m(C)m^{(n)}(C)\to m(C)m(n)(C)→m(C) for every fixed Borel set CCC. Thus the observation kernel Q(⋅∣a,x)Q(\cdot\mid a,x)Q(⋅∣a,x) is setwise continuous even at the accumulation action zero. Example 4.1, pp. 13–14.

A belief z∈P(X)z\in\mathcal P(X)z∈P(X) assigns a probability to each hidden state. For a current belief zzz and an action aaa, the joint law R(⋅∣z,a)R(\cdot\mid z,a)R(⋅∣z,a) covers the next state and its observation, and R′(⋅∣z,a)R'(\cdot\mid z,a)R′(⋅∣z,a) is the observation marginal. A filter H(z,a,y)H(z,a,y)H(z,a,y) is a conditional law of the next state given observation yyy: it satisfies the disintegration identity (3.3). The belief transition qH(⋅∣z,a)q_H(\cdot\mid z,a)qH​(⋅∣z,a) is the law of the random posterior H(z,a,y)H(z,a,y)H(z,a,y) when yyy has law R′(⋅∣z,a)R'(\cdot\mid z,a)R′(⋅∣z,a). These objects are defined in (3.1)–(3.5) of the paper. Different versions of HHH agree only almost everywhere for each current belief and action, yet they yield the same belief transition.

Formalization targets

Example 4.1: failure of weak belief continuity

At the uniform belief z=(1/2,1/2)z=(1/2,1/2)z=(1/2,1/2), action zero reveals nothing, while every positive-index action produces the same two-point posterior law:

qH(⋅∣z,0)=δz,qH(⋅∣z,1/n)=34δ(1/3,2/3)+14δ(1,0)(n≥1).q_H(\cdot\mid z,0)=\delta_z,\qquad q_H(\cdot\mid z,1/n)=\tfrac34\delta_{(1/3,2/3)}+\tfrac14\delta_{(1,0)} \quad(n\ge1).qH​(⋅∣z,0)=δz​,qH​(⋅∣z,1/n)=43​δ(1/3,2/3)​+41​δ(1,0)​(n≥1).

The mission's goal asserts that PPP is total-variation continuous, QQQ is a setwise-continuous probability kernel, a valid filter exists, and the positive-index laws fail to converge weakly to the zero-action law as 1/n→01/n\to01/n→0. It states the failure for every filter satisfying (3.3), so it is a property of the POMDP rather than of a selected version of its conditional distribution. The milestone list follows the paper's assertions: setwise convergence of the dyadic laws, setwise continuity of QQQ, the two posterior values almost everywhere, and the explicit laws of qHq_HqH​.

Significance

The example identifies a precise boundary of the sufficient condition in Theorems 3.6 and 3.7. Those theorems use total-variation continuity of the observation law to obtain a weakly continuous belief transition; replacing that premise by setwise continuity makes their continuity conclusion false in this model. The explicit law also shows the size of the discrepancy: the positive-index belief transition has two fixed atoms, while the zero-action transition is concentrated at the original belief. Example 4.1, pp. 13–14.

A machine-checked development would add reusable definitions of sequential kernel continuity, the joint state-observation law, conditional filters, and the belief transition on the Borel space of probability measures. It would also formalize a concrete setwise-convergent sequence of densities whose induced posterior laws do not converge weakly. The paper proves the example; the Lean statements in this proposal are open proof targets, with no completed machine-checked proof claimed here.

Difficulty

Setwise continuity controls Q(C∣a,x)Q(C\mid a,x)Q(C∣a,x) for each fixed Borel set CCC. A posterior is a probability measure determined by the observation, and the belief transition tests the distribution of those posterior measures. The observation regions distinguishing the states change with nnn, so convergence on fixed observation sets does not imply convergence of their posterior distributions. This is the gap in the immediate argument that tries to pass from setwise convergence of QQQ to weak convergence of qHq_HqH​.

The conditional filter adds a second precision issue. Equation (3.3) determines its values only almost surely under the observation marginal, and different versions may disagree on null observations. Any claim about a posterior value at every yyy would exceed the source; the posterior milestone therefore has an almost-everywhere qualifier, while the goal quantifies over every valid version of HHH.

Formalization scope

Lean represents XXX by Fin 2, with 0 for the paper's state 1 and 1 for state 2; YYY by unitInterval; and AAA by the specified subtype of R\mathbb RR. Positive integers are a subtype of N\mathbb NN, so the density and the action 1/n1/n1/n are never evaluated at n=0n=0n=0. The measure mmm is volume on unitInterval, which has total mass one. The density follows (4.1) literally, with open zero intervals and density two at their endpoints. The general transition PPP and observation law QQQ use Mathlib kernels; the example instantiates them directly. The goal includes the probability-kernel property of the concrete QQQ and existence of a valid HHH, guarding the universal filter clause against vacuity.

Weak, setwise, and total-variation convergence are sequential, as on page 4. The weak test class consists of bounded continuous real functions; the setwise test class consists of all Borel sets. The target space of qHq_HqH​ carries the Borel σ-algebra of the weak topology on P(X)\mathcal P(X)P(X), stated explicitly because Mathlib's default measurable structure on that type is the Giry σ-algebra. The filter predicate requires measurability in yyy for fixed (z,a)(z,a)(z,a) in both structures. Joint measurability in (z,a,y)(z,a,y)(z,a,y) is described in the paper but is not encoded in this local predicate; only the fixed-parameter property is used by (3.3) and (3.5). This choice and its quantifier effect are recorded in the moderation notes.

No cost, initial observation kernel, initial prior, or discount factor enters Example 4.1. Contributions on the dyadic law, kernel continuity, conditional posterior values, and the resulting belief-law identity are all within scope. A filter hard-coded to return the prior belief, a discrete topology on AAA, or a claim for only one version of HHH would change the counterexample and is excluded by the definitions and goal.

Selected references

  • E. A. Feinberg, P. O. Kasyanov, and M. Z. Zgurovsky, Partially Observable Total-Cost Markov Decision Processes with Weakly Continuous Transition Probabilities, Mathematics of Operations Research 41(2), 2016; arXiv:1401.2168v2. This mission follows the pinned preprint's numbering and pages.
6 thms1 active userReviewed
Dynamic ProgrammingMarkov ChainOperations Research·Captain: mikedeng1

Partially Observable Total-Cost Markov Decision Processes with Weakly Continuous Transition Probabilities 1: The Belief MDP Satisfies Assumption (W*) Under Total-Variation-Continuous ObservationsResearch Paper

Motivation

A partially observable Markov decision process (POMDP) models a controller that cannot see the state of the system it steers, only noisy observations of it. Inventory control with unrecorded demand, machine maintenance with imperfect inspection, and target tracking from sensor readings all have this form. The standard way to solve a POMDP is to replace it by its belief MDP, also called the completely observable MDP (COMDP): the state becomes the posterior distribution of the hidden state given everything observed so far. The belief MDP is fully observed, so the theory of Markov decision processes applies to it, but only if the belief MDP satisfies that theory's hypotheses. For MDPs with Borel state and action spaces, unbounded costs and possibly noncompact action sets, the hypothesis that yields optimality equations, convergence of value iterations and stationary optimal policies is Assumption (W*) of Feinberg, Kasyanov and Zadoianchuk (Math. Oper. Res. 2012): the cost is K\mathbb KK-inf-compact and the transition probability is weakly continuous.

Before Feinberg, Kasyanov and Zgurovsky (arXiv:1401.2168v2, published in Math. Oper. Res. 41(2), 2016), sufficient conditions for the belief MDP to be weakly continuous required a weakly continuous filter (Hernández-Lerma, 1989, as recalled on p. 11 of arXiv:1401.2168v2), a condition that fails in many natural models. The paper replaces it by a condition on the primitive data of the POMDP.

Setting

Let X\mathbb XX (states), Y\mathbb YY (observations) and A\mathbb AA (actions) be Borel subsets of Polish spaces. Let P(dx′∣x,a)P(dx' \mid x, a)P(dx′∣x,a) be the transition kernel on X\mathbb XX given X×A\mathbb X \times \mathbb AX×A, Q(dy∣a,x)Q(dy \mid a, x)Q(dy∣a,x) the observation kernel on Y\mathbb YY given A×X\mathbb A \times \mathbb XA×X, and c:X×A→R∪{+∞}c : \mathbb X \times \mathbb A \to \mathbb R \cup \{+\infty\}c:X×A→R∪{+∞} a Borel one-step cost bounded below. P(X)\mathbb P(\mathbb X)P(X) is the set of probability measures on X\mathbb XX with the topology of weak convergence (∫f dμn→∫f dμ\int f\,d\mu_n \to \int f\,d\mu∫fdμn​→∫fdμ for bounded continuous fff). A sequence μn\mu_nμn​ converges setwise if μn(C)→μ(C)\mu_n(C) \to \mu(C)μn​(C)→μ(C) for every Borel CCC, and in the total variation if sup⁡f∣∫f dμn−∫f dμ∣→0\sup_f |\int f\,d\mu_n - \int f\,d\mu| \to 0supf​∣∫fdμn​−∫fdμ∣→0 over Borel fff with values in [−1,1][-1, 1][−1,1]. A kernel is weakly (setwise, total-variation) continuous if xn→xx_n \to xxn​→x implies the corresponding convergence of the measures.

For a belief z∈P(X)z \in \mathbb P(\mathbb X)z∈P(X) and an action aaa, let R(B×C∣z,a)=∫X∫BQ(C∣a,x′) P(dx′∣x,a) z(dx)R(B \times C \mid z, a) = \int_{\mathbb X} \int_B Q(C \mid a, x')\,P(dx' \mid x, a)\,z(dx)R(B×C∣z,a)=∫X​∫B​Q(C∣a,x′)P(dx′∣x,a)z(dx) be the joint law of the next state and observation (3.1), and R′(C∣z,a)=R(X×C∣z,a)R'(C \mid z, a) = R(\mathbb X \times C \mid z, a)R′(C∣z,a)=R(X×C∣z,a) its observation marginal (3.2). A filter is a kernel H(dx∣z,a,y)H(dx \mid z, a, y)H(dx∣z,a,y) with R(B×C∣z,a)=∫CH(B∣z,a,y) R′(dy∣z,a)R(B \times C \mid z, a) = \int_C H(B \mid z, a, y)\,R'(dy \mid z, a)R(B×C∣z,a)=∫C​H(B∣z,a,y)R′(dy∣z,a) (3.3): H(z,a,y)H(z, a, y)H(z,a,y) is the posterior after observing yyy. The belief MDP has state space P(X)\mathbb P(\mathbb X)P(X), action set A\mathbb AA, cost cˉ(z,a)=∫c(x,a) z(dx)\bar c(z, a) = \int c(x, a)\,z(dx)cˉ(z,a)=∫c(x,a)z(dx) (3.8), and transition q(⋅∣z,a)q(\cdot \mid z, a)q(⋅∣z,a), the law of H(z,a,y)H(z, a, y)H(z,a,y) when y∼R′(⋅∣z,a)y \sim R'(\cdot \mid z, a)y∼R′(⋅∣z,a) (3.5). A cost is K\mathbb KK-inf-compact on S1×S2\mathbb S_1 \times \mathbb S_2S1​×S2​ if its restriction to K×S2K \times \mathbb S_2K×S2​ has compact level sets for every compact K⊆S1K \subseteq \mathbb S_1K⊆S1​.

Formalization targets

Goal: Theorem 3.6

If either Assumption (D) (0≤α<10 \le \alpha < 10≤α<1) or Assumption (P) (c≥0c \ge 0c≥0 and 0≤α≤10 \le \alpha \le 10≤α≤1) holds, ccc is bounded below and K\mathbb KK-inf-compact on X×A\mathbb X \times \mathbb AX×A, PPP is weakly continuous and QQQ is continuous in the total variation, then the belief MDP satisfies Assumption (W*):

cˉ is bounded below and K-inf-compact on P(X)×A,(zn,an)→(z,a)  ⟹  q(⋅∣zn,an)→q(⋅∣z,a) weakly.\bar c \text{ is bounded below and } \mathbb K\text{-inf-compact on } \mathbb P(\mathbb X) \times \mathbb A, \qquad (z_n, a_n) \to (z, a) \implies q(\cdot \mid z_n, a_n) \to q(\cdot \mid z, a) \text{ weakly}.cˉ is bounded below and K-inf-compact on P(X)×A,(zn​,an​)→(z,a)⟹q(⋅∣zn​,an​)→q(⋅∣z,a) weakly.

The two halves

  • Cost half (Theorem 3.4, through Lemma 6.1 and Lemma 5.1(ii)): cˉ\bar ccˉ is bounded below by the same constant as ccc and K\mathbb KK-inf-compact.
  • Transition half (Theorem 3.7): under the hypotheses on PPP and QQQ, R′R'R′ is setwise continuous and Assumption (H) holds: whenever zn→zz_n \to zzn​→z and an→aa_n \to aan​→a, along a subsequence H(znk,ank,y)→H(z,a,y)H(z_{n_k}, a_{n_k}, y) \to H(z, a, y)H(znk​​,ank​​,y)→H(z,a,y) weakly for R′(⋅∣z,a)R'(\cdot \mid z, a)R′(⋅∣z,a)-almost every yyy. Theorem 3.5 then gives weak continuity of qqq. The route runs through Theorem 5.2, Lemma 5.3, Corollary 5.4, Theorem 5.5 and Lemma 5.6.

Significance

Theorem 3.6 supplies the hypothesis under which the cited general MDP theory gives stationary optimal policies, optimality equations and convergence of value iterations for total-cost POMDPs in terms of PPP, QQQ and ccc alone, without compactness of the action set, boundedness of the cost or continuity of the filter. It covers, for example, observation kernels with densities continuous in (a,x)(a, x)(a,x), and the paper applies it to inventory control with incomplete records and to the Kalman filter. Example 4.1 of the same paper shows that total-variation continuity of QQQ cannot be weakened to setwise continuity.

The result is proved in the paper; this proposal poses its formal statement and the supporting results as proof targets. Completing them would provide a checked bridge from POMDPs to the general MDP theory under weak continuity, together with reusable facts about weak, setwise and total-variation convergence: equicontinuity of integrals against weakly continuous kernels (Theorem 5.2), the convergence-in-probability criterion (Theorem 5.5), and the generalized Fatou lemmas.

Difficulty

The posterior H(z,a,y)H(z, a, y)H(z,a,y) is defined only for R′(⋅∣z,a)R'(\cdot \mid z, a)R′(⋅∣z,a)-almost every yyy and has no continuous version in general, so weak continuity of qqq cannot be read off from continuity of HHH. The tempting argument, that continuous PPP and QQQ make the Bayes posterior move continuously, fails: Y\mathbb YY need not be countable and QQQ need not have densities. What has to be controlled is the joint law RRR on sets (O1∖O2)×C(\mathcal O_1 \setminus \mathcal O_2) \times C(O1​∖O2​)×C uniformly in the observation set CCC; convergence for each fixed CCC is not enough. Turning this uniform control into almost-sure convergence of posteriors along a subsequence requires a diagonal argument over a countable base of X\mathbb XX. On the cost side, the difficulty is that ccc may take the value +∞+\infty+∞ and is not continuous; K\mathbb KK-inf-compactness has to be transported through integration against a weakly converging sequence of measures.

Formalization scope

All definitions live in FeinbergPOMDP.WStar.Model. X,Y,A\mathbb X, \mathbb Y, \mathbb AX,Y,A are separable metrizable spaces with their Borel σ-algebras (PolishSpace is not assumed, since Borel subsets of Polish spaces need not be Polish); X\mathbb XX is additionally nonempty and standard Borel in Theorems 3.6 and 3.7, where a filter must exist. P(X)\mathbb P(\mathbb X)P(X) is Mathlib's ProbabilityMeasure with the weak topology. All continuity notions are the paper's sequential ones; total variation is the supremum over Borel functions with values in [−1,1][-1, 1][−1,1], written in ε\varepsilonε–NNN form, while the set-indexed suprema of (5.15) and (5.17) keep their set form. The cost is written c=ℓ+fc = \ell + fc=ℓ+f with ℓ∈R\ell \in \mathbb Rℓ∈R and fff taking values in [0,∞][0, \infty][0,∞], so cˉ=ℓ+∫f dz\bar c = \ell + \int f\,dzcˉ=ℓ+∫fdz, exact even when c=+∞c = +\inftyc=+∞ on a set of positive measure. K\mathbb KK-inf-compactness is the published FeinbergLiang.ACOE.KInfCompact. A filter is a map H:P(X)→A→Y→P(X)H : \mathbb P(\mathbb X) \to \mathbb A \to \mathbb Y \to \mathbb P(\mathbb X)H:P(X)→A→Y→P(X) that satisfies (3.3) and is jointly Borel measurable in (z,a,y)(z,a,y)(z,a,y), as the paper's stochastic kernel is; qqq is the image of R′R'R′ under H(z,a,⋅)H(z, a, \cdot)H(z,a,⋅) for the Borel σ-algebra of the weak topology on P(X)\mathbb P(\mathbb X)P(X).

Since qqq does not depend on the choice of filter, every statement about qqq is made for every filter satisfying (3.3). A fixed filter that ignores the observation, such as H(z,a,y)=zH(z, a, y) = zH(z,a,y)=z, would make qqq a point mass and the weak continuity trivial; it satisfies (3.3) only when observations carry no information, and is excluded by quantifying over all filters.

The goal states the Assumption (W*) conclusion of Theorem 3.6. Its tail, "statements (i)–(vi) of Theorem 3.1 hold", invokes Theorem 2.1 of Feinberg et al. (2012) applied to the belief MDP; it is cited rather than proved in this paper. Generalized Fatou's lemmas 5.1(i) and 5.1(ii) are referenced from the published Schäl (1993) items. A complete development needs weak and setwise convergence on ProbabilityMeasure, portmanteau-type arguments, disintegration (Measure.condKernel) and the identification of the Borel and Giry σ-algebras on P(X)\mathbb P(\mathbb X)P(X); the last is missing from Mathlib and is a welcome contribution in its own right. Proofs of any milestone, and alternative arguments for Theorem 3.7, are welcome.

Selected references

  • E. A. Feinberg, P. O. Kasyanov, M. Z. Zgurovsky, Partially Observable Total-Cost Markov Decision Processes with Weakly Continuous Transition Probabilities, arXiv:1401.2168v2 (2014); Mathematics of Operations Research 41(2), 2016. https://arxiv.org/abs/1401.2168v2
  • E. A. Feinberg, P. O. Kasyanov, N. V. Zadoianchuk, Average-cost Markov decision processes with weakly continuous transition probabilities, Mathematics of Operations Research 37(4), 2012. https://doi.org/10.1287/moor.1120.0555
  • D. P. Bertsekas, S. E. Shreve, Stochastic Optimal Control: The Discrete-Time Case, Academic Press, 1978. https://web.mit.edu/dimitrib/www/soc.html
  • M. Schäl, Average optimality in dynamic programming with general state space, Mathematics of Operations Research 18(1), 1993. https://doi.org/10.1287/moor.18.1.163
14 thms1 active userReviewed
Bandit AlgorithmsMachine LearningStatistics·Captain: mikedeng1

The Multi-Armed Bandit Problem with Covariates II: Adaptively Binned Successive Elimination Has Regret at Most Cn(K log K/n)^{β(α+1)/(2β+d)} Under Smoothness and Margin ConditionsResearch Paper

Motivation

An online decision maker often sees information about a customer, patient, or request before choosing an action. A bandit with covariates models this setting when the reward of every available action is initially unknown and only the reward of the chosen action is observed. A method must learn which action works best for which covariates while paying for exploratory choices. Perchet and Rigollet studied this problem for smooth reward functions and measured how quickly a policy can approach the performance of an oracle that knows every arm's mean reward Perchet and Rigollet, 2013.

A fixed grid divides the covariate space at one resolution everywhere. It can spend pulls distinguishing arms in regions where their difference is already large, while leaving too few observations in regions where they are close. Adaptively binned successive elimination (ABSE) instead refines a region only while several arms remain plausible there. The paper proves a regret rate controlled jointly by the smoothness of the reward functions and the probability of encountering near ties. Its Theorem 5.1 is the target of this mission.

Setting

There are K≥2K\ge2K≥2 arms and covariates XtX_tXt​ in the unit cube X=[0,1]d⊂Rd\mathcal X=[0,1]^d\subset\mathbb R^dX=[0,1]d⊂Rd, with d≥1d\ge1d≥1. At time ttt, the policy sees XtX_tXt​, chooses an arm, and sees that arm's reward. The random vectors consisting of the covariate and all potential arm rewards are independent across times. Covariates have a common probability law PXP_XPX​ with a density between constants 0<c‾≤c‾0<\underline c\le\overline c0<c​≤c. Every reward lies in [0,1][0,1][0,1]; for arm iii its conditional mean given Xt=xX_t=xXt​=x is f(i)(x)∈[0,1]f^{(i)}(x)\in[0,1]f(i)(x)∈[0,1]. The paper does not require the complete reward vectors to be identically distributed across time. These conventions come from §3.1, pp. 10–11.

At xxx, write f∗(x)=max⁡if(i)(x)f^*(x)=\max_i f^{(i)}(x)f∗(x)=maxi​f(i)(x). Write f♯(x)f^\sharp(x)f♯(x) for the largest mean strictly below f∗(x)f^*(x)f∗(x), and set f♯(x)=f∗(x)f^\sharp(x)=f^*(x)f♯(x)=f∗(x) when every arm ties. The difference f∗(x)−f♯(x)f^*(x)-f^\sharp(x)f∗(x)−f♯(x) is the positive gap wherever a genuinely suboptimal arm exists. The smoothness condition requires ∣f(i)(x)−f(i)(x′)∣≤L∥x−x′∥β|f^{(i)}(x)-f^{(i)}(x')|\le L\|x-x'\|^\beta∣f(i)(x)−f(i)(x′)∣≤L∥x−x′∥β for every arm, with the Euclidean norm, 0<β≤10<\beta\le10<β≤1, and L>0L>0L>0. The margin condition, with α>0\alpha>0α>0, 0<δ0<10<\delta_0<10<δ0​<1, and C0>0C_0>0C0​>0, says

PX{0<f∗(X)−f♯(X)≤δ}≤C0δα(0≤δ≤δ0).P_X\{0<f^*(X)-f^\sharp(X)\le\delta\}\le C_0\delta^\alpha\qquad(0\le\delta\le\delta_0).PX​{0<f∗(X)−f♯(X)≤δ}≤C0​δα(0≤δ≤δ0​).

These are the assumptions defining the paper's class MXK(α,β,L)\mathcal M^K_{\mathcal X}(\alpha,\beta,L)MXK​(α,β,L) §3.2, pp. 12–13.

ABSE starts with the whole cube. Each live dyadic cell runs a copy of successive elimination on its inherited arms. One round pulls each remaining arm once. After a complete round, an arm is eliminated if its empirical mean falls more than twice the radius U(s,n∣B∣d)U(s,n|B|^d)U(s,n∣B∣d) below the largest empirical mean, where ∣B∣|B|∣B∣ is the cell's side length and UUU uses log⁡‾(x)=max⁡{log⁡x,1}\overline{\log}(x)=\max\{\log x,1\}log​(x)=max{logx,1}. The cell splits after its specified stopping round if it has at least two active arms and is below the maximum depth k0k_0k0​. Each child receives the survivors. The policy takes d,β,L,K,nd,\beta,L,K,nd,β,L,K,n as parameters and uses the current covariate and rewards it actually observed Policy 3, p. 22.

Formalization targets

Goal: Theorem 5.1, finite margin exponent

The goal is the finite-α\alphaα part of Theorem 5.1, p. 22. Fix d≥1d\ge1d≥1, 0<β≤10<\beta\le10<β≤1, L>0L>0L>0, α>0\alpha>0α>0 with αβ≤d\alpha\beta\le dαβ≤d, 0<δ0<10<\delta_0<10<δ0​<1, C0>0C_0>0C0​>0 and 0<c‾≤c‾0<\underline c\le\overline c0<c​≤c. There is one constant C>0C>0C>0, depending on these parameters but not on KKK, nnn or the machine, such that for every K≥2K\ge2K≥2, every n≥Klog⁡Kn\ge K\log Kn≥KlogK and every machine of the class,

ERn(π~)≤Cn(Klog⁡Kn)β(α+1)2β+d.\mathbb E R_n(\widetilde\pi)\le Cn\left(\frac{K\log K}{n}\right)^{\frac{\beta(\alpha+1)}{2\beta+d}}.ERn​(π)≤Cn(nKlogK​)2β+dβ(α+1)​.

Here π~\widetilde\piπ is ABSE and RnR_nRn​ is its cumulative regret against an arm of maximum conditional mean at each covariate. The goal asserts the shape of the rate with an unspecified constant, which is the paper's form.

Milestones

In attack order:

  1. Lemma A.1 (p. 28): a maximal Hoeffding–Azuma inequality, P{∃t≤T:Zˉt≥2(b−a)2t−1log⁡(4T/(δt))}≤δ\mathbb P\{\exists t\le T:\bar Z_t\ge\sqrt{2(b-a)^2t^{-1}\log(4T/(\delta t))}\}\le\deltaP{∃t≤T:Zˉt​≥2(b−a)2t−1log(4T/(δt))​}≤δ for bounded martingale differences.
  2. (5.3) (p. 21): ℓB≤Cℓ∣B∣−2β log⁡‾(n∣B∣2β+d)\ell_B\le C_\ell|B|^{-2\beta}\,\overline{\log}(n|B|^{2\beta+d})ℓB​≤Cℓ​∣B∣−2βlog​(n∣B∣2β+d), the number of rounds a cell spends before it may split.
  3. (5.10) (p. 26): if arm iii has gap above 8c0∣B∣β8c_0|B|^\beta8c0​∣B∣β at a point of BBB, then the cell mean of the arm ıˇ\check\imathˇ that is best at that point exceeds the cell mean of iii by at least 32εB,ℓB\tfrac32\varepsilon_{B,\ell_B}23​εB,ℓB​​.
  4. The display before (5.6) (p. 25): ∑∣B∣=2−kqB≤c22k(d−βα)\sum_{|B|=2^{-k}}q_B\le c_2 2^{k(d-\beta\alpha)}∑∣B∣=2−k​qB​≤c2​2k(d−βα), where qBq_BqB​ is the conditional probability of a near tie inside BBB.
  5. (5.5) (p. 24): the pathwise decomposition rnborn(B)=rnlive(B)+∑B′∈burst(B)rnborn(B′)r_n^{\mathrm{born}}(B)=r_n^{\mathrm{live}}(B)+\sum_{B'\in\mathrm{burst}(B)}r_n^{\mathrm{born}}(B')rnborn​(B)=rnlive​(B)+∑B′∈burst(B)​rnborn​(B′) of regret along the tree.
  6. Theorem 5.1, α=∞\alpha=\inftyα=∞ (p. 22): when the margin region {0<f⋆−f♯≤δ0}\{0<f^\star-f^\sharp\le\delta_0\}{0<f⋆−f♯≤δ0​} has probability zero, ERn(π~)≤CKlog⁡n\mathbb E R_n(\widetilde\pi)\le CK\log nERn​(π)≤CKlogn.

Significance

The result provides a performance guarantee for a policy that adapts its spatial resolution to observed evidence, under a class of nonparametric reward functions. It permits a direct comparison with a fixed-bin policy under the same smoothness and margin conditions. The paper proves the rate and discusses optimality in a minimax sense when K=2K=2K=2 §5, pp. 22–27. This mission asks for a machine-checked proof of the stated finite-α\alphaα regime; the Lean goal and milestones are currently open statements, so the paper's proof is the mathematical source rather than a completed formal proof.

A completed development would make the stochastic model, adaptive policy, cell geometry, and concentration estimates reusable for other contextual bandit analyses.

Difficulty

The main difficulty is the feedback loop between decisions and observations. A cell receives samples only when the covariate lands there, and its available arms depend on eliminations in its ancestors. Applying an ordinary concentration inequality to a fixed sequence of arm rewards does not by itself justify the bound at these random birth and stopping times. The formal proof must track the inherited active sets and the samples collected after each birth, then combine local error bounds across the tree. The appendix's maximal inequality controls a range of round counts; the regret argument also uses the small probability of covariates with nearly tied arm means §5, pp. 23–27 and Lemma A.1, p. 28.

Formalization scope

Lean uses EuclideanSpace ℝ (Fin d) and a probability measure on its unit cube. The lower density bound ensures every valid dyadic cell has positive probability, so its conditional mean is defined by a real quotient. A consistent half-open boundary assignment turns the paper's cells, which are specified up to Lebesgue-null boundaries, into a partition on every path. The policy is a recursive state machine: at time ttt it receives XtX_tXt​, chooses an arm from its prior state, and updates with only the reward of that arm. It does not receive the unknown means or any unobserved potential reward. The existential constant in the goal is quantified before KKK, nnn, the sample space, and the machine.

The finite-α\alphaα goal includes αβ≤d\alpha\beta\le dαβ≤d. The printed theorem omits it, although the proof uses it. Without it, a two-arm machine with constant means 111 and 000 has regret from the policy's first exploratory round, while the printed bound tends to zero when α\alphaα is large. The α=∞\alpha=\inftyα=∞ clause reads the margin condition as PX(0<f⋆−f♯≤δ0)=0P_X(0<f^\star-f^\sharp\le\delta_0)=0PX​(0<f⋆−f♯≤δ0​)=0, the limit of C0δαC_0\delta^\alphaC0​δα for δ<1\delta<1δ<1, because the page does not define it. The (5.3) milestone uses log⁡‾\overline{\log}log​: the printed ordinary logarithm can be nonpositive at a finest cell, while the left side is a positive integer. The cell parameter TB=n∣B∣dT_B=n|B|^dTB​=n∣B∣d and bursting after a full round follow Policy 3 and the proof where nearby prose or pseudocode differs. Lemma A.1 adds EZ1=0\mathbb EZ_1=0EZ1​=0 and a<ba<ba<b, both implicit on the page. The constants c2c_2c2​ and CℓC_\ellCℓ​ are existential and placed after the class parameters, because the page's explicit c2=c1α/c‾c_2=c_1^\alpha/\underline cc2​=c1α​/c​ omits the margin constant C0C_0C0​.

The probabilistic error displays (5.6)–(5.9) and (5.11)–(5.14) are not milestones. They are statements about the events AB={I‾B⊆IB⊆IˉB}\mathcal A_B=\{\underline{\mathcal I}_B\subseteq\mathcal I_B\subseteq\bar{\mathcal I}_B\}AB​={I​B​⊆IB​⊆IˉB​} and GB\mathcal G_BGB​, which the page leaves undefined for cells that are never born or never finish ℓB\ell_BℓB​ rounds. The decomposition (5.8) also leaves out arms of I‾B\underline{\mathcal I}_BI​B​ that were eliminated at an ancestor of BBB. A faithful encoding of these events is welcome as a contribution.

A complete proof needs measurable recursive policies, conditional expectations and independence, dyadic-cell volume and nesting facts, concentration for arm observations at adaptive times, and sums over tree depths. Contributions establishing these interfaces and the paper's intermediate bounds are welcome. The result must apply to every machine in the stated class. Each of the following would change the theorem: a policy that reads the unknown f(i)f^{(i)}f(i) or unobserved rewards (the oracle has zero regret), a constant chosen after the machine, a covariate law that is not a probability measure on the cube, Hölder continuity in the sup norm instead of the Euclidean norm, or dropping the density lower bound (cells of zero mass would make fˉB\bar f_Bfˉ​B​ a junk value). The paper's proofs are in §5 (pp. 23–27) and its Appendix (p. 28).

Selected references

  • Jean-Yves Audibert and Alexandre B. Tsybakov, Fast learning rates for plug-in classifiers, Annals of Statistics 35(2):608–633, 2007, DOI:10.1214/009053606000001217.
  • Eyal Even-Dar, Shie Mannor and Yishay Mansour, Action elimination and stopping conditions for the multi-armed bandit and reinforcement learning problems, Journal of Machine Learning Research 7:1079–1105, 2006, jmlr.org/papers/v7/evendar06a.
  • Philippe Rigollet and Assaf Zeevi, Nonparametric bandits with covariates, COLT 2010, pp. 54–66, arXiv:1003.1630.
  • Vianney Perchet and Philippe Rigollet, The multi-armed bandit problem with covariates, Annals of Statistics 41(2):693–721, 2013, pinned arXiv:1110.6084v3, DOI:10.1214/13-AOS1101.
13 thms1 active userReviewed
Operations ResearchOptimizationStatistics·Captain: mikedeng1

Recovering Best Statistical Guarantees via the Empirical Divergence-Based Distributionally Robust Optimization 3: On a Finite Support the Empirical DRO Equals the Burg-Divergence DRO with χ²₁ RadiusResearch Paper

Motivation

Distributionally robust optimization (DRO) evaluates a decision against a collection of probability distributions near an empirical one. In a finite scenario model, the familiar construction puts a divergence ball around the histogram of observed scenarios. Lam's empirical construction instead gives one weight to each observation, even when several observations have the same value. The distinction matters because a sample can contain repeated values and can omit some points of the stated support. The paper uses a Burg-entropy ball on those observation weights to recover statistical guarantees for an expected loss, despite the empirical ball having low probability of containing the true distribution. Lam, §§2.1–2.3, pp. 5–9.

Proposition 1 identifies what this observation-level construction becomes when the random quantity has finite support: its lower and upper DRO values coincide with values from a Burg-divergence ball on the support histogram. This links a method defined on individual data points to the standard finite-scenario formulation. The radius uses a one-degree-of-freedom chi-square quantile because the paper calibrates an expected value at a fixed decision; it does not use the k−1k-1k−1 degrees of freedom associated with estimating a full kkk-point distribution. Lam, Proposition 1, p. 11, and the discussion on pp. 8–9.

Setting

Fix a decision xxx and a real loss h(x;si)h(x;s_i)h(x;si​) at each of kkk support points s1,…,sks_1,\ldots,s_ks1​,…,sk​. A sample of size nnn has an assignment c(j)c(j)c(j) telling which support point produced observation jjj. Write nin_ini​ for the number of assignments to point iii and p^i=ni/n\hat p_i=n_i/np^​i​=ni​/n for the empirical histogram. The support points need not all appear in the sample. The decision and loss only determine a finite vector of loss values; no probability space is needed for this deterministic comparison. Lam, §2.1, p. 6, and Proposition 1, p. 11.

The empirical Burg ball Un(η)\mathcal U_n(\eta)Un​(η) consists of observation weights w1,…,wnw_1,\ldots,w_nw1​,…,wn​ that are nonnegative, sum to one, and obey −n−1∑jlog⁡(nwj)≤η-n^{-1}\sum_j\log(nw_j)\le\eta−n−1∑j​log(nwj​)≤η. The support Burg ball UBurg′(η)\mathcal U'_{\mathrm{Burg}}(\eta)UBurg′​(η) consists of support weights p1,…,pkp_1,\ldots,p_kp1​,…,pk​ that are nonnegative, sum to one, are absolutely continuous with respect to p^\hat pp^​, and obey −∑ip^ilog⁡(pi/p^i)≤η-\sum_i\hat p_i\log(p_i/\hat p_i)\le\eta−∑i​p^​i​log(pi​/p^​i​)≤η. At an observed point, finite Burg divergence requires pi>0p_i>0pi​>0; at an unobserved point, absolute continuity requires pi=0p_i=0pi​=0. Lam, (13)–(15), pp. 5–6; (19), p. 8; (28), p. 11.

For a support weight vector ppp, the expected loss is ∑ipih(x;si)\sum_i p_i h(x;s_i)∑i​pi​h(x;si​). For an observation weight vector www, it is ∑jwjh(x;sc(j))\sum_j w_j h(x;s_{c(j)})∑j​wj​h(x;sc(j)​). The empirical lower and upper values are the minimum and maximum over Un(η)\mathcal U_n(\eta)Un​(η); the support values use UBurg′(η)\mathcal U'_{\mathrm{Burg}}(\eta)UBurg′​(η). The same empirical ball is used for both bounds. Lam, (26)–(28), p. 11.

Formalization targets

The goal is Proposition 1. For every n>0n>0n>0, every assignment of the sample to finite support points, every fixed decision xxx, and every κ≥0\kappa\ge0κ≥0, set η=κ/(2n)\eta=\kappa/(2n)η=κ/(2n). Then

min⁡w∈Un(η)∑j=1nwjh(x;sc(j))=min⁡p∈UBurg′(η)∑i=1kpih(x;si),\min_{w\in\mathcal U_n(\eta)}\sum_{j=1}^n w_jh(x;s_{c(j)}) =\min_{p\in\mathcal U'_{\mathrm{Burg}}(\eta)}\sum_{i=1}^k p_ih(x;s_i),w∈Un​(η)min​j=1∑n​wj​h(x;sc(j)​)=p∈UBurg′​(η)min​i=1∑k​pi​h(x;si​), max⁡w∈Un(η)∑j=1nwjh(x;sc(j))=max⁡p∈UBurg′(η)∑i=1kpih(x;si).\max_{w\in\mathcal U_n(\eta)}\sum_{j=1}^n w_jh(x;s_{c(j)}) =\max_{p\in\mathcal U'_{\mathrm{Burg}}(\eta)}\sum_{i=1}^k p_ih(x;s_i).w∈Un​(η)max​j=1∑n​wj​h(x;sc(j)​)=p∈UBurg′​(η)max​i=1∑k​pi​h(x;si​).

The printed proposition takes κ=χ1,1−α2\kappa=\chi^2_{1,1-\alpha}κ=χ1,1−α2​. The deterministic identity holds at every nonnegative κ\kappaκ, including zero. The two milestones record the source proof's statements that a feasible support vector gives a feasible observation vector with the same objective value, and conversely that a feasible observation vector gives a feasible support vector with the same objective value. Lam, Proposition 1, p. 11, proof p. 24.

Significance

The identity lets a finite-support user compute the empirical DRO values with one variable per distinct support point instead of one per observation, without changing either bound. It also shows precisely how Lam's empirical calibration relates to a conventional histogram-based divergence model: the ambiguity sets have different coordinates, but the resulting extrema of any loss constant on each class agree. This is the finite-support interpretation stated immediately after Theorem 2, and it explains the radius comparison on pp. 8–9. Lam, pp. 8–9 and Proposition 1, p. 11.

The result is proved in the paper; the remaining task is a machine-checked proof of the two transformations and the equality of extrema. The published formal library already contains the observation-level ϕ\phiϕ-divergence ball and its robust upper mean, and a shared setting supplies the Burg generator and lower empirical value. This mission adds the support histogram, its Burg ball, and the bridge between the two representations. Those definitions and transformations can also be reused for other finite-support losses. The paper's separate pointwise coverage and consistency results are goals of other missions in this series. Lam, Theorems 2–3 and Proposition 1, pp. 11–12.

Difficulty

Repeated observations make the two weight spaces different: the empirical problem has nnn coordinates and the support problem has kkk. Simply identifying coordinates loses multiplicities in both the normalization and the logarithmic divergence. Empty classes create a second boundary case, since no observation weight can be assigned to an unobserved support point. A statement that silently permits such mass changes the support optimization value. The full reduction must preserve feasible objective values in both directions, for the lower and upper extrema. Lam, proof of Proposition 1, p. 24.

Formalization scope

Lean uses Fin n for observations and Fin k for support points, with zero-based indices, a class map ccc, real loss values, and the paper's Burg generator valued in EReal. Nonpositive likelihood ratios receive +∞+\infty+∞. The empirical ball reuses the published PhiDivRobust.Counterpart.probUncertaintySet; its upper value is the published GenEmpLik.Expansion.robustMean. The lower value is a real infimum over the same ball. The condition n>0n>0n>0 and the nonnegative radius make both feasible value sets nonempty and bounded, so the real infimum and supremum mean the paper's minimum and maximum. The histogram can have zero entries. Lam, (19), (26)–(28), pp. 8, 11.

Display (28) omits a condition from the preceding definition of a divergence ball: ppp must be absolutely continuous with respect to p^\hat pp^​. Without it, an unobserved support point can receive positive mass and Proposition 1 is false. The formal support ball includes p^i=0⇒pi=0\hat p_i=0\Rightarrow p_i=0p^​i​=0⇒pi​=0 from (13)–(14). It also requires pi>0p_i>0pi​>0 where p^i>0\hat p_i>0p^​i​>0, because Lean's real logarithm returns zero at zero whereas the paper's Burg divergence is infinite there. These conditions, together with n>0n>0n>0, exclude trivial values from empty balls, division by zero, and logarithms at invalid weights. The decision set Θ\ThetaΘ is omitted because the paper fixes one xxx and the result holds for every such xxx. Lam, (13)–(14), pp. 5–6, and (28), p. 11.

The development needs finite-sum identities, elementary properties of the probability simplex, the logarithm's concavity, and bounds on the extrema of bounded nonempty finite-dimensional value sets. Contributions that establish the two direction statements, or reusable facts about aggregation of finite probability weights, directly advance this mission. The paper's uniform calibration theorem, Theorem 4, is outside this scope: it needs a Donsker class, a data-dependent chi-square process, and its supremum quantile, none of which enters Proposition 1. Lam, Theorem 4, pp. 13–14.

Selected references

  • Henry Lam, Recovering Best Statistical Guarantees via the Empirical Divergence-Based Distributionally Robust Optimization, arXiv:1605.09349v1, 2016; published in Operations Research 67(4), 2019. arXiv, DOI.
8 thms1 active userReviewed
Machine LearningOptimal TransportStatistics·Captain: mikedeng1

Wasserstein Distributionally Robust Optimization and Variation Regularization III: With Probability 1 − e^(−t), E_Pn[(ρ − d(z, D_f))₊] ≤ Cρ² + 2ρ·E[Rademacher] + ρ√(t/2n) for Every fResearch Paper

Motivation

Wasserstein distributionally robust optimization (DRO) replaces the empirical risk EPn[f]\mathbb E_{P_n}[f]EPn​​[f] of a loss fff by its worst case over all distributions within Wasserstein distance ρ\rhoρ of the empirical law PnP_nPn​. Gao, Chen and Kleywegt (arXiv:1712.06050; Operations Research, 2024, doi:10.1287/opre.2022.2383) show that, for small radii, this worst-case risk is the empirical risk plus ρ\rhoρ times a variation of the loss (a norm of its gradient or local slope), up to a remainder. This links DRO to gradient-norm regularization, as used in adversarial training and in data-driven newsvendor and pricing models.

For smooth losses the remainder is of order ρ2\rho^2ρ2. Many losses used in practice are only piecewise smooth: hinge, ReLU networks, newsvendor costs, θz∧1\theta z \wedge 1θz∧1. Near the set Df\mathcal D_fDf​ where the pieces meet, the first-order expansion fails, and the paper's Theorem 1(I)–(II) carries an extra boundary term M EPn[(ρ−d(z,Df))+]M\,\mathbb E_{P_n}[(\rho - d(z,\mathcal D_f))_+]MEPn​​[(ρ−d(z,Df​))+​]. Theorem 1(III) shows that this term is small with high probability, uniformly over the loss class, provided the data distribution does not concentrate near the non-smooth points. This mission formalizes that probabilistic bound.

Setting

Let (Z,d)(\mathcal Z, d)(Z,d) be a metric space with its Borel σ\sigmaσ-algebra, PtrueP_{\rm true}Ptrue​ a probability measure on Z\mathcal ZZ, and z1,…,znz_1, \dots, z_nz1​,…,zn​ an i.i.d. sample from PtrueP_{\rm true}Ptrue​, with empirical law Pn=1n∑iδziP_n = \frac1n\sum_i \delta_{z_i}Pn​=n1​∑i​δzi​​. The sampling law is the product P⊗=Ptrue⊗nP_\otimes = P_{\rm true}^{\otimes n}P⊗​=Ptrue⊗n​ on Zn\mathcal Z^nZn, and E⊗\mathbb E_\otimesE⊗​ is its expectation. The formal statements also assume the metric space is second countable.

Let F\mathcal FF be a set of losses f:Z→Rf : \mathcal Z \to \mathbb Rf:Z→R, and for each fff let Df⊆Z\mathcal D_f \subseteq \mathcal ZDf​⊆Z be its non-smooth set: for a piecewise smooth fff with pieces Zf,1,…,Zf,Kf\mathcal Z_{f,1}, \dots, \mathcal Z_{f,K_f}Zf,1​,…,Zf,Kf​​, the union of the pairwise intersections of their closures (Assumption 1). The paper assumes Df\mathcal D_fDf​ is PtrueP_{\rm true}Ptrue​-null. The distance to a set is d(z,D)=inf⁡z~∈Dd(z,z~)d(z, D) = \inf_{\tilde z\in D} d(z,\tilde z)d(z,D)=infz~∈D​d(z,z~), with d(z,∅)=+∞d(z,\varnothing) = +\inftyd(z,∅)=+∞.

The objects of the mission are:

  • the boundary gap (ρ−d(z,Df))+(\rho - d(z, \mathcal D_f))_+(ρ−d(z,Df​))+​ and the near indicator 1{d(z,Df)<ρ}\mathbf 1\{d(z,\mathcal D_f) < \rho\}1{d(z,Df​)<ρ};
  • the indicator class Iρ={z↦1{d(z,Df)<ρ}:f∈F, Df≠∅}\mathcal I_\rho = \{z\mapsto \mathbf 1\{d(z,\mathcal D_f)<\rho\} : f\in\mathcal F,\ \mathcal D_f\neq\varnothing\}Iρ​={z↦1{d(z,Df​)<ρ}:f∈F, Df​=∅};
  • the empirical Rademacher complexity Rn(H)=Eσ[sup⁡h∈H1n∑iσih(zi)]\mathfrak R_n(\mathcal H) = \mathbb E_\sigma\bigl[\sup_{h\in\mathcal H}\frac1n\sum_i\sigma_i h(z_i)\bigr]Rn​(H)=Eσ​[suph∈H​n1​∑i​σi​h(zi​)] of a class H\mathcal HH, with i.i.d. uniform signs σi\sigma_iσi​, and the Rademacher complexity E⊗[Rn(H)]\mathbb E_\otimes[\mathfrak R_n(\mathcal H)]E⊗​[Rn​(H)];
  • Assumption 3 (bounded density):
lim sup⁡δ↓0 sup⁡f∈F: Df≠∅Ptrue{z:0<d(z,Df)<δ}δ<∞.\limsup_{\delta\downarrow 0}\ \sup_{f\in\mathcal F:\,\mathcal D_f\ne\varnothing}\frac{P_{\rm true}\{z : 0 < d(z,\mathcal D_f) < \delta\}}{\delta} < \infty.δ↓0limsup​ f∈F:Df​=∅sup​δPtrue​{z:0<d(z,Df​)<δ}​<∞.

In Lean the boundary gap is the shared WassVarReg.PInf.boundaryGap; nearIndicator, indicatorClass, empRademacher, rademacher and Assumption3 are in WassVarReg.Boundary. The sample law is the published MinimaxWass.DataDep.sampleLaw.

Formalization targets

Goal: Theorem 1(III), first display (p. 8)

Under Assumption 3 there exist ρˉ,C>0\bar\rho, C > 0ρˉ​,C>0 such that for all n≥1n\ge1n≥1, t>0t>0t>0 and 0≤ρ<ρˉ0\le\rho<\bar\rho0≤ρ<ρˉ​, with probability at least 1−e−t1-e^{-t}1−e−t, for every f∈Ff\in\mathcal Ff∈F,

EPn[(ρ−d(z,Df))+]≤Cρ2+2ρ E⊗[Rn(Iρ)]+ρt2n.\mathbb E_{P_n}\bigl[(\rho - d(z,\mathcal D_f))_+\bigr] \le C\rho^2 + 2\rho\,\mathbb E_\otimes[\mathfrak R_n(\mathcal I_\rho)] + \rho\sqrt{\frac{t}{2n}}.EPn​​[(ρ−d(z,Df​))+​]≤Cρ2+2ρE⊗​[Rn​(Iρ​)]+ρ2nt​​.

The constants CCC and ρˉ\bar\rhoρˉ​ are left unspecified; they depend only on Assumption 3.

Milestones

  1. Lemma EC.5 (p. ec4): for a class H\mathcal HH of measurable functions with values in [0,M][0,M][0,M], with probability at least 1−e−t1-e^{-t}1−e−t, EPn[h]≤EPtrue[h]+2 E⊗[Rn(H)]+Mt/(2n)\mathbb E_{P_n}[h] \le \mathbb E_{P_{\rm true}}[h] + 2\,\mathbb E_\otimes[\mathfrak R_n(\mathcal H)] + M\sqrt{t/(2n)}EPn​​[h]≤EPtrue​​[h]+2E⊗​[Rn​(H)]+Mt/(2n)​ for all h∈Hh\in\mathcal Hh∈H.
  2. Assumption 3 bounds the true mass (proof of Lemma EC.11, p. ec11): EPtrue[1{d(z,Df)<ρ}]≤Cρ\mathbb E_{P_{\rm true}}[\mathbf 1\{d(z,\mathcal D_f)<\rho\}] \le C\rhoEPtrue​​[1{d(z,Df​)<ρ}]≤Cρ for 0≤ρ<ρˉ0\le\rho<\bar\rho0≤ρ<ρˉ​, uniformly in fff.
  3. Lemma EC.11 (p. ec11): with probability at least 1−e−t1-e^{-t}1−e−t, for every fff, EPn[1{d(z,Df)<ρ}]≤Cρ+2 E⊗[Rn(Iρ)]+t/(2n)\mathbb E_{P_n}[\mathbf 1\{d(z,\mathcal D_f)<\rho\}] \le C\rho + 2\,\mathbb E_\otimes[\mathfrak R_n(\mathcal I_\rho)] + \sqrt{t/(2n)}EPn​​[1{d(z,Df​)<ρ}]≤Cρ+2E⊗​[Rn​(Iρ​)]+t/(2n)​.
  4. Gap versus mass (p. ec11): EPn[(ρ−d(z,Df))+]≤ρ EPn[1{d(z,Df)<ρ}]\mathbb E_{P_n}[(\rho-d(z,\mathcal D_f))_+] \le \rho\,\mathbb E_{P_n}[\mathbf 1\{d(z,\mathcal D_f)<\rho\}]EPn​​[(ρ−d(z,Df​))+​]≤ρEPn​​[1{d(z,Df​)<ρ}].

Significance

The result. With ρn=ρ0/n\rho_n = \rho_0/\sqrt nρn​=ρ0​/n​ and a class Iρ\mathcal I_\rhoIρ​ of bounded complexity, the right side of the goal is Op(1/n)O_p(1/n)Op​(1/n). Combined with Theorem 1(I)–(II), it gives the paper's conclusion that ppp-Wasserstein DRO with p∈(1,∞]p\in(1,\infty]p∈(1,∞] is asymptotically equivalent to empirical variation regularization, min⁡f{EPn[f]+ρnVPn,q(f)}+Op(1/n)\min_f\{\mathbb E_{P_n}[f] + \rho_n\mathcal V_{P_n,q}(f)\} + O_p(1/n)minf​{EPn​​[f]+ρn​VPn​,q​(f)}+Op​(1/n), for non-smooth losses, uniformly over the class. Without it, the boundary term is only bounded by MρM\rhoMρ, and the expansion has no better accuracy than first order. Assumption 3 is the condition under which the non-smooth set does not matter statistically.

Formalizing it. The result is proved in the paper; it has no machine-checked proof. Its proof composes a uniform deviation bound (McDiarmid's inequality plus symmetrization) with a bounded-density estimate. The mission produces a general Rademacher deviation bound for bounded classes on an arbitrary metric space with a measurability guard, reusable by any statistical-learning mission, and fixes the printed constant of Lemma EC.5 (M2t/nM\sqrt{2t/n}M2t/n​ on the page, Mt/(2n)M\sqrt{t/(2n)}Mt/(2n)​ from its proof).

Difficulty

The deterministic parts are short. The difficulty is in Lemma EC.5. McDiarmid's inequality must be applied to the uniform deviation over an arbitrary, possibly uncountable, class. Symmetrization must then be carried out with Bochner integrals over Zn×Zn\mathcal Z^n\times\mathcal Z^nZn×Zn and a finite sign space. Neither tool is in Mathlib in this form. A tempting shortcut, bounding each fixed hhh by Hoeffding's inequality and taking a union bound, fails: the class is infinite, and the bound must hold simultaneously for all fff inside a single event. A second point is Assumption 3: it controls only the shell 0<d(z,Df)<δ0 < d(z,\mathcal D_f) < \delta0<d(z,Df​)<δ, so the mass of {d=0}\{d = 0\}{d=0} has to come from the null closure of Df\mathcal D_fDf​.

Formalization scope

  • Space and sample. Z\mathcal ZZ is a metric space with [MeasurableSpace Z] [BorelSpace Z] [SecondCountableTopology Z] in the goal-level statements. Second countability is an added technical hypothesis used throughout this paper's mission series. A sample is ω : Fin n → Z with n≥1n\ge1n≥1, and its law is Measure.pi. Empirical means are 1n∑i:Fin n\frac1n\sum_{i : \mathrm{Fin}\,n}n1​∑i:Finn​.
  • Distance to a set. Metric.infEDist, which is +∞+\infty+∞ on ∅\varnothing∅. A loss with Df=∅\mathcal D_f=\varnothingDf​=∅ thus has boundary gap 000 and is excluded from Iρ\mathcal I_\rhoIρ​, as on the page. The real-valued Metric.infDist (which is 000 on ∅\varnothing∅) would give such a loss gap ρ\rhoρ and is not used.
  • The non-smooth sets are an arbitrary family D : (Z → ℝ) → Set Z. The goal-level statements assume each Df\mathcal D_fDf​ is closed and PtrueP_{\rm true}Ptrue​-null, the standing consequence of its construction in Assumption 1. The details of the piecewise smooth loss model are not needed for this probability estimate.
  • Probabilities are written as P⊗{ω:¬ ∀f∈F,… }≤e−tP_\otimes\{\omega : \neg\,\forall f\in\mathcal F, \dots\}\le e^{-t}P⊗​{ω:¬∀f∈F,…}≤e−t, with the quantifier over F\mathcal FF inside the event (an outer measure for uncountable F\mathcal FF).
  • Rademacher complexity. The sign expectation is an average over {±1}n\{\pm1\}^n{±1}n. Suprema are real suprema of classes with values in [0,M][0,M][0,M]; an empty class gives 000. E⊗\mathbb E_\otimesE⊗​ is a Bochner integral.
  • Measurability guard (added): the uniform deviation, the empirical Rademacher complexity and the double-sample deviation of the class are measurable. This holds, for instance, for countable F\mathcal FF, and the paper leaves it implicit.
  • Quantifier order. ρˉ\bar\rhoρˉ​ and CCC are chosen before nnn and ttt, which is stronger than the page's order. The radius satisfies ρ≥0\rho\ge0ρ≥0.
  • Not stated. The second display of Theorem 1(III) and the "Moreover" part of Lemma EC.11, which use Dudley's entropy integral: their proof bounds the empirical Rademacher complexity, while Lemma EC.5 needs its expectation.
  • Trivialization ruled out. The Rademacher term is that of the specific class Iρ\mathcal I_\rhoIρ​ built from the sets Df\mathcal D_fDf​. It is not a free constant, and no hypothesis assumes the conclusion of Lemma EC.5.

Contributions welcome: a general McDiarmid inequality for product measures, a symmetrization lemma for measurable classes, and proofs of the milestones in any order.

Selected references

  • R. Gao, X. Chen, A. J. Kleywegt, Wasserstein Distributionally Robust Optimization and Variation Regularization, Operations Research, 2024. arXiv:1712.06050, doi:10.1287/opre.2022.2383
  • S. Shalev-Shwartz, S. Ben-David, Understanding Machine Learning: From Theory to Algorithms, Cambridge University Press, 2014, Theorem 26.5. doi:10.1017/CBO9781107298019
  • C. McDiarmid, On the method of bounded differences, Surveys in Combinatorics, LMS Lecture Note Series 141, 1989. doi:10.1017/CBO9781107359949.008
  • M. J. Wainwright, High-Dimensional Statistics: A Non-Asymptotic Viewpoint, Cambridge University Press, 2019, Ch. 5. doi:10.1017/9781108627771
8 thms1 active userReviewed
Algorithmic Game TheoryOperations Research·Captain: mikedeng1

Prophet Inequalities Made Easy: Stochastic Optimization by Pricing Nonstochastic Inputs II: Posted Expected Weakly (α, β₁, β₂)-Balanced Prices Earn 1/(α(2β₁ + 4β₂)) of the Expected Welfare of ALGResearch Paper

Motivation

A prophet inequality compares an online decision maker, who sees random inputs one at a time and must commit irrevocably, with a "prophet" who sees all inputs in advance. The classic single-item version (Krengel–Sucheston, Samuel-Cahn) says that a well-chosen fixed threshold earns at least half of the expected maximum. In mechanism design the same statement reads: a posted-price mechanism, which offers each arriving buyer a take-it-or-leave-it menu of prices, extracts a constant fraction of the optimal expected welfare, and does so in a way that is truthful, simple and order oblivious.

Dütting, Feldman, Kesselheim and Lucier (SIAM J. Comput. 49 (2020)) reduce the construction of such prices to a deterministic question. If, for every fixed valuation profile vvv, one can find prices pvp^vpv that are balanced — high enough to pay for the welfare an early allocation destroys, low enough that whatever is still feasible stays affordable — then the expected prices Ev~[pv~]\mathbb E_{\tilde v}[p^{\tilde v}]Ev~​[pv~], suitably scaled, give a prophet inequality. This is their Theorem 3.2 for (α,β)(\alpha,\beta)(α,β)-balanced prices. For combinatorial auctions with bundles of size at most ddd, MPH-kkk valuations and related settings, the natural prices only satisfy a weaker condition, in which the affordability bound may also spend a multiple of the benchmark welfare. Theorem 3.5, the subject of this mission, is the extension theorem for that weak balancedness; it is what turns the paper's Theorem 4.1 (bundle size ddd) and Theorem 1.6 (MPH-kkk) into O(d)O(d)O(d)- and O(k)O(k)O(k)-approximate posted-price mechanisms.

Setting

There are nnn agents N={1,…,n}N=\{1,\dots,n\}N={1,…,n}. Agent iii has an outcome space XiX_iXi​ containing a null outcome ∅\emptyset∅; outcome profiles are x∈X=X1×⋯×Xnx\in X=X_1\times\cdots\times X_nx∈X=X1​×⋯×Xn​. For S⊆NS\subseteq NS⊆N, xSx_SxS​ gives xix_ixi​ to the agents in SSS and ∅\emptyset∅ to the others, and x[i−1]x_{[i-1]}x[i−1]​ keeps only the agents before iii. The feasible profiles form a set F⊆X\mathcal F\subseteq XF⊆X that is downward closed (x∈F⇒xS∈Fx\in\mathcal F\Rightarrow x_S\in\mathcal Fx∈F⇒xS​∈F).

Agent iii has a valuation vi:Xi→[0,1]v_i:X_i\to[0,1]vi​:Xi​→[0,1], drawn independently from a known distribution Di\mathcal D_iDi​; D=∏iDi\mathcal D=\prod_i\mathcal D_iD=∏i​Di​. The welfare of xxx is v(x)=∑ivi(xi)v(x)=\sum_i v_i(x_i)v(x)=∑i​vi​(xi​), and v(OPT(v,S))=sup⁡x∈Sv(x)v(\mathrm{OPT}(v,S))=\sup_{x\in S}v(x)v(OPT(v,S))=supx∈S​v(x) (zero on S=∅S=\emptysetS=∅). An allocation rule ALG\mathrm{ALG}ALG maps every valuation profile to a feasible outcome.

A pricing rule assigns prices pi(xi∣y)∈[0,∞]p_i(x_i\mid y)\in[0,\infty]pi​(xi​∣y)∈[0,∞] to outcome xix_ixi​ offered to agent iii given the partial allocation yyy, with pi(xi∣y)=∞p_i(x_i\mid y)=\inftypi​(xi​∣y)=∞ when (xi,y−i)∉F(x_i,y_{-i})\notin\mathcal F(xi​,y−i​)∈/F. In the posted-price mechanism the agents are approached in index order; agent iii sees the prices pi(⋅∣x[i−1])p_i(\cdot\mid x_{[i-1]})pi​(⋅∣x[i−1]​) and buys a utility-maximizing outcome, utility being vi(xi)−pi(xi∣x[i−1])v_i(x_i)-p_i(x_i\mid x_{[i-1]})vi​(xi​)−pi​(xi​∣x[i−1]​).

A family (Fx)x∈X(\mathcal F_x)_{x\in X}(Fx​)x∈X​ of sets of profiles is exchange compatible if (yi,x−i)∈F(y_i,x_{-i})\in\mathcal F(yi​,x−i​)∈F for every x∈Xx\in Xx∈X, y∈Fxy\in\mathcal F_xy∈Fx​ and iii. Given α>0\alpha>0α>0 and β1,β2≥0\beta_1,\beta_2\ge0β1​,β2​≥0, a pricing rule pvp^vpv is weakly (α,β1,β2)(\alpha,\beta_1,\beta_2)(α,β1​,β2​)-balanced (Definition 3.4) with respect to ALG\mathrm{ALG}ALG and (Fx)(\mathcal F_x)(Fx​) if for every x∈Fx\in\mathcal Fx∈F

(a)∑ipiv(xi∣x[i−1]) ≥ 1α(v(ALG(v))−v(OPT(v,Fx))),\text{(a)}\quad \sum_i p^v_i(x_i\mid x_{[i-1]})\ \ge\ \tfrac1\alpha\bigl(v(\mathrm{ALG}(v))-v(\mathrm{OPT}(v,\mathcal F_x))\bigr),(a)i∑​piv​(xi​∣x[i−1]​) ≥ α1​(v(ALG(v))−v(OPT(v,Fx​))), (b)∑ipiv(xi′∣x[i−1]) ≤ β1 v(OPT(v,Fx))+β2 v(ALG(v))for all x′∈Fx.\text{(b)}\quad \sum_i p^v_i(x'_i\mid x_{[i-1]})\ \le\ \beta_1\,v(\mathrm{OPT}(v,\mathcal F_x))+\beta_2\,v(\mathrm{ALG}(v))\quad\text{for all }x'\in\mathcal F_x.(b)i∑​piv​(xi′​∣x[i−1]​) ≤ β1​v(OPT(v,Fx​))+β2​v(ALG(v))for all x′∈Fx​.

A collection (pv)v(p^v)_{v}(pv)v​ is weakly balanced if one exchange-compatible family works for every vvv.

Formalization targets

Goal: Theorem 3.5 (p. 551)

If (pv)v(p^v)_v(pv)v​ is weakly (α,β1,β2)(\alpha,\beta_1,\beta_2)(α,β1​,β2​)-balanced with β1+β2≥1/α\beta_1+\beta_2\ge1/\alphaβ1​+β2​≥1/α, then with

δ=1β1+max⁡{2β2,1/α},pi(xi∣y)=Ev~[piv~(xi∣y)],\delta=\frac1{\beta_1+\max\{2\beta_2,1/\alpha\}},\qquad p_i(x_i\mid y)=\mathbb E_{\tilde v}\bigl[p^{\tilde v}_i(x_i\mid y)\bigr],δ=β1​+max{2β2​,1/α}1​,pi​(xi​∣y)=Ev~​[piv~​(xi​∣y)],

the posted-price mechanism with prices δp\delta pδp, approaching agents in index order, satisfies

Ev[∑ivi(xi(v))] ≥ 1α(2β1+4β2) Ev[v(ALG(v))].\mathbb E_v\Bigl[\sum_i v_i(x_i(v))\Bigr]\ \ge\ \frac1{\alpha(2\beta_1+4\beta_2)}\,\mathbb E_v\bigl[v(\mathrm{ALG}(v))\bigr].Ev​[i∑​vi​(xi​(v))] ≥ α(2β1​+4β2​)1​Ev​[v(ALG(v))].

Milestones (Appendix A, pp. 559–560)

  1. Each agent's utility in the run is nonnegative.
  2. (A.1): the total expected utility is at least the value of a ghost-sample deviation x′(v,v′)∈Fx(v)x'(v,v')\in\mathcal F_{x(v)}x′(v,v′)∈Fx(v)​ minus its expected price.
  3. (A.2): property (b) bounds that price by δβ1E[v~(OPT(v~,Fx(v)))]+δβ2E[v~(ALG(v~))]\delta\beta_1\mathbb E[\tilde v(\mathrm{OPT}(\tilde v,\mathcal F_{x(v)}))]+\delta\beta_2\mathbb E[\tilde v(\mathrm{ALG}(\tilde v))]δβ1​E[v~(OPT(v~,Fx(v)​))]+δβ2​E[v~(ALG(v~))].
  4. (A.3): E[∑iui]≥(1−δβ1) E[v~(OPT(v~,Fx(v)))]−δβ2 E[v~(ALG(v~))]\mathbb E[\sum_iu_i]\ge(1-\delta\beta_1)\,\mathbb E[\tilde v(\mathrm{OPT}(\tilde v,\mathcal F_{x(v)}))]-\delta\beta_2\,\mathbb E[\tilde v(\mathrm{ALG}(\tilde v))]E[∑i​ui​]≥(1−δβ1​)E[v~(OPT(v~,Fx(v)​))]−δβ2​E[v~(ALG(v~))].
  5. (A.4): property (a) bounds the expected revenue from below.
  6. Case 1, β2≥1/(2α)\beta_2\ge1/(2\alpha)β2​≥1/(2α), δ=1/(β1+2β2)\delta=1/(\beta_1+2\beta_2)δ=1/(β1​+2β2​): the goal's bound.
  7. Case 2, β2<1/(2α)\beta_2<1/(2\alpha)β2​<1/(2α), δ=1/(β1+1/α)\delta=1/(\beta_1+1/\alpha)δ=1/(β1​+1/α): the bound with constant (1−αβ2)/(1+αβ1)(1-\alpha\beta_2)/(1+\alpha\beta_1)(1−αβ2​)/(1+αβ1​).
  8. Case 2's last step: (1−αβ2)/(1+αβ1)≥1/(α(2β1+4β2))(1-\alpha\beta_2)/(1+\alpha\beta_1)\ge1/(\alpha(2\beta_1+4\beta_2))(1−αβ2​)/(1+αβ1​)≥1/(α(2β1​+4β2​)) under the theorem's constraints.

Significance

The result. Theorem 3.5 is the second of the paper's two extension theorems. Its hypothesis is a deterministic, full-information property of prices, which is often checkable by a direct combinatorial argument; its conclusion is a Bayesian guarantee for an online, truthful mechanism. The paper instantiates it with weakly (1,1,d−1)(1,1,d-1)(1,1,d−1)-balanced fractional prices for combinatorial auctions with bundle size at most ddd (Theorem 4.1) and for MPH-kkk valuations (Theorem 1.6), obtaining approximation ratios linear in ddd and kkk, where (α,β)(\alpha,\beta)(α,β)-balancedness alone would not apply.

Formalizing it. The theorem is proved on paper; no machine-checked version exists. The mission produces a formal account of the sequential posted-price mechanism with history-dependent prices (an online run defined by recursion over the agents), of the ghost-sample argument, and of the case analysis in δ\deltaδ. The definitions are general enough that the paper's applications — and other balanced-price arguments — can be stated on top of them. The printed proof has a slip: (A.1) is printed as an equality but only the inequality "≥\ge≥" holds; the formalization states the inequality.

Difficulty

Two steps carry the content. First, the utility bound (A.1): agent iii could have bought its part of the deviation x′(v,v′)x'(v,v')x′(v,v′), but that deviation is defined from the run on a different profile. The argument needs that the partial allocation x[i−1](v)x_{[i-1]}(v)x[i−1]​(v) does not depend on viv_ivi​ — the mechanism is online — and that exchanging viv_ivi​ with an independent copy vi′v'_ivi′​ preserves the joint law. Making this exchange rigorous on a product measure, and handling an empty Fx(v)\mathcal F_{x(v)}Fx(v)​ (where no deviation exists and the bound must come from ui≥0u_i\ge0ui​≥0), is where the work lies. Second, v(OPT(v,Fx))v(\mathrm{OPT}(v,\mathcal F_x))v(OPT(v,Fx​)) is a supremum that need not be attained, so "buy your part of OPT" must be replaced by near-optimal selections and a limit.

The tempting shortcut of treating weak balancedness as (α,β1+β2)(\alpha,\beta_1+\beta_2)(α,β1​+β2​)-balancedness and applying Theorem 3.2 does not work: property (b) of Definition 3.1 bounds prices by a multiple of v(OPT(v,Fx))v(\mathrm{OPT}(v,\mathcal F_x))v(OPT(v,Fx​)) alone, which can be zero while v(ALG(v))v(\mathrm{ALG}(v))v(ALG(v)) is not.

Formalization scope

Agents are Fin n, 0-based: the paper's agent iii is index i−1i-1i−1, and "the order they are indexed" is the order of Fin n. Prices take values in ENNReal with ∞=⊤\infty=\top∞=⊤; price sums are never converted to reals, and property (a) passes its right side through ENNReal.ofReal (a negative bound becomes 000, which is the page's meaning). v(OPT(v,S))v(\mathrm{OPT}(v,S))v(OPT(v,S)) is the real supremum sSup, zero on the empty set, meaningful because values lie in [0,1][0,1][0,1].

Valuations are given by type spaces ViV_iVi​ and an evaluation vali(w,xi)∈[0,1]\mathrm{val}_i(w,x_i)\in[0,1]vali​(w,xi​)∈[0,1]. The type spaces are countable with discrete σ-algebras — a disclosed restriction of the paper's general distributions, which makes every function measurable and every bounded one integrable. The null outcome is assumed free at feasible partial allocations (piv~(∅∣y)=0p^{\tilde v}_i(\emptyset\mid y)=0piv~​(∅∣y)=0 for y∈Fy\in\mathcal Fy∈F): every pricing rule in the paper has this property, and the proof needs its consequence ui(v)≥0u_i(v)\ge0ui​(v)≥0. Exchange compatibility is required for all x∈Xx\in Xx∈X, as printed. The agents' behaviour is a choice rule that sees only the agent's own type and the earlier agents' outcomes and maximizes utility at every feasible partial allocation; at infeasible ones no requirement is made. The constant δ\deltaδ is written literally inside the posted prices, never as a free parameter.

The milestones (A.1)–(A.4) are stated for an arbitrary δ>0\delta>0δ>0, and (A.1)–(A.2) for an arbitrary selector into Fx(v)\mathcal F_{x(v)}Fx(v)​ in place of the maximizer, with an explicit zero when Fx(v)=∅\mathcal F_{x(v)}=\emptysetFx(v)​=∅; both are generalizations of the printed steps. A statement in which the agents' optimality hypothesis can never be met, or in which an infinite price sum is read as 000, would make the goal vacuous; the encoding above excludes both, and a sorry-free one-agent instance in the workspace confirms that all hypotheses of the goal can hold together.

Needed infrastructure: product measures and the swap of one coordinate with an independent copy, lower integrals of ENNReal-valued prices, the recursion defining the run, and near-maximizers of bounded suprema. The model and mechanism definitions are reusable for the paper's other extension theorem (Theorem 3.2) and its applications. Contributions are welcome on any milestone; the pure inequality of Case 2 is self-contained.

Selected references

  • P. Dütting, M. Feldman, T. Kesselheim, B. Lucier, Prophet inequalities made easy: Stochastic optimization by pricing nonstochastic inputs, SIAM Journal on Computing 49(3):540–582, 2020. https://doi.org/10.1137/20M1323850
  • U. Krengel, L. Sucheston, Semiamarts and finite values, Bulletin of the AMS 83(4):745–747, 1977. https://doi.org/10.1090/S0002-9904-1977-14378-4
  • E. Samuel-Cahn, Comparison of threshold stop rules and maximum for independent nonnegative random variables, Annals of Probability 12(4):1213–1216, 1984. https://doi.org/10.1214/aop/1176993150
  • R. Kleinberg, S. M. Weinberg, Matroid prophet inequalities, STOC 2012, 123–136. https://doi.org/10.1145/2213977.2213991
12 thms1 active userReviewed
Bandit AlgorithmsInformation TheoryOperations Research·Captain: mikedeng1

Data-Based Dynamic Pricing and Inventory Control with Censored Demand and Limited Price Changes 2: Any Policy with at Most m Price Changes Has Regret at Least KT^(1/(m+1)) on a Bernoulli InstanceResearch Paper

Why limited price changes matter

A seller learning a demand curve normally adjusts prices as new observations arrive. A price-change limit prevents that continual adjustment: a price may have to be chosen before the seller knows which demand curve it faces, then held through a block of sales. This matters when repricing is operationally costly, constrained by contracts, or visible to customers. Chen, Chao and Wang study pricing jointly with inventory replenishment under censored demand. Their Theorem 2 asks how much regret is unavoidable when the seller can change price at most mmm times, regardless of the learning policy it uses (SSRN revision, pp. 10–11, 35–42).

The lower bound is paired with their Theorem 1: in the well-separated setting an algorithm attains regret proportional to T1/(m+1)T^{1/(m+1)}T1/(m+1) for large horizons, while Theorem 2 asserts the same exponent cannot be improved uniformly over admissible policies (SSRN revision, p. 10). The paper appeared as an Operations Research technical note in 2020 (DOI 10.1287/opre.2020.1993); this mission cites the authors' SSRN revision dated 2020-02-10, whose page and result numbers are used throughout.

Bernoulli pricing setting

The horizon has TTT periods. At each period the seller posts a price pt∈[1,6]p_t\in[1,6]pt​∈[1,6] and replenishes inventory to one unit. Demand dtd_tdt​ is either zero or one, so all demand can be served and the observed sale equals demand. The unknown scalar parameter zzz belongs to [1/6,5/6][1/6,5/6][1/6,5/6]. Conditional on the posted price, the probability of one unit of demand is

ν(p,z)=max⁡{0,1−pz/2}.\nu(p,z)=\max\{0,1-pz/2\}.ν(p,z)=max{0,1−pz/2}.

Holding and shortage costs are zero. Hence expected one-period revenue is rz(p)=pν(p,z)r_z(p)=p\nu(p,z)rz​(p)=pν(p,z). The clairvoyant benchmark is G∗(z)=1/(2z)G^*(z)=1/(2z)G∗(z)=1/(2z): it is attained by the feasible price p∗=1/zp^*=1/zp∗=1/z. The local Lean sanity check proves this maximum property, so the benchmark does not rely on a default value of a supremum. The model is the explicit Bernoulli construction in the proof of Theorem 2 (SSRN revision, p. 35).

An admissible randomized policy may draw random bits and use past posted prices and observed sales when choosing its next price. Its action cannot depend on future demand or on the unknown parameter. Every realized seed and demand path must satisfy the price-change constraint

∑t=1T−11{pt≠pt+1}≤m,m≥1.\sum_{t=1}^{T-1}\mathbf 1\{p_t\ne p_{t+1}\}\le m,\qquad m\ge1.t=1∑T−1​1{pt​=pt+1​}≤m,m≥1.

The regret Rπ(T,z)R_\pi(T,z)Rπ​(T,z) is the expected sum of G∗(z)−rz(pt)G^*(z)-r_z(p_t)G∗(z)−rz​(pt​) under the policy and the demand law at parameter zzz. The paper's admissibility and change constraint are stated on pp. 5–6; its randomized policy factorization is displayed in (49) on p. 35 (SSRN revision).

Formalization targets

The goal is the minimax reading of Theorem 2. For each fixed m≥1m\ge1m≥1, there is a positive constant K7K_7K7​ and a horizon threshold T0T_0T0​, independent of the policy and the horizon, such that

∀T≥T0  ∀π admissible with at most m changes,∃z∈[1/6,5/6]: Rπ(T,z)≥K7T1/(m+1).\forall T\ge T_0\;\forall\pi\text{ admissible with at most }m\text{ changes},\quad \exists z\in[1/6,5/6]:\ R_\pi(T,z)\ge K_7T^{1/(m+1)}.∀T≥T0​∀π admissible with at most m changes,∃z∈[1/6,5/6]: Rπ​(T,z)≥K7​T1/(m+1).

The parameter is selected after the policy. The proof's objective (56) is a worst-case parameter choice for each policy (SSRN revision, p. 37). The printed theorem sentence can also be read as one fixed instance defeating every algorithm; that order is false because a policy supplied with that instance's zzz can post 1/z1/z1/z in every period.

Five milestones expose the source's quantitative ingredients. The hierarchy (50) constructs parameters zζz_\zetazζ​ from signs ζ∈{±1}m+1\zeta\in\{\pm1\}^{m+1}ζ∈{±1}m+1 and scales εℓ=T−ℓ/[2(m+1)]\varepsilon_\ell=T^{-\ell/[2(m+1)]}εℓ​=T−ℓ/[2(m+1)]. Its parameters lie in [1/6,5/6][1/6,5/6][1/6,5/6] for large TTT. Different signs are separated at scale εℓ/8\varepsilon_\ell/8εℓ​/8. Lemma 1 bounds one-period regret by the squared error in a nearest-instance selection. The pointwise inequality underlying (57) gives a cost of at least 10−4εℓ210^{-4}\varepsilon_\ell^210−4εℓ2​ for a wrong sign. A corrected Lemma 2 supplies a quadratic Bernoulli relative-entropy bound under an explicit interior condition (SSRN revision, pp. 36–37).

What the result establishes

The exponent 1/(m+1)1/(m+1)1/(m+1) quantifies the cost of the change budget: fixing mmm leaves a polynomial regret floor even though the Bernoulli demand law is one-dimensional and observed without censoring. Since the constructed instance fixes inventory at a level that serves every demand, a lower bound here also applies to policies in the joint pricing and inventory setting whenever that instance is admitted. It identifies the information loss caused by limited repricing, independently of lost sales and censoring.

The paper states the lower bound and gives a proof, but its printed Lemma 2 is false at prices where ν(p,z)=0\nu(p,z)=0ν(p,z)=0. The subsequent Lemma 3 applies that unrestricted bound. Therefore the source proof does not, as printed, establish the minimax claim for every admissible price policy. The Lean goal is an open target; the four elementary milestones and corrected information bound are also open Lean statements. Formalizing this mission requires an argument that handles prices at the zero-demand kink, or a different route to the lower bound. No machine-checked proof of this specific result is claimed here.

Where the argument is difficult

The natural comparison of two nearby parameters uses the relative entropy between their observed demand laws. Away from the kink, this information is quadratic in the parameter difference. At a price that makes demand impossible under one parameter but possible under the other, the divergence in the printed direction can be linear in that difference. The unrestricted 108∣z−z′∣2108|z-z'|^2108∣z−z′∣2 bound of Lemma 2 therefore fails at exactly the prices an arbitrary policy may choose. Restricting the policy's prices would change the target. A faithful proof must retain those prices while still controlling the policy's information and regret (SSRN revision, Lemmas 2–3, pp. 36–40).

Formalization scope

Lean represents the horizon by Fin T; its index zero is period one in the paper. Bernoulli demand paths have type Fin T → Bool. A policy's price at a period depends on a random seed and only the earlier demand bits; with inventory level one, these bits are exactly the observed sales. A probability measure on the seed models adaptive randomization, including a pre-sampled stream of random choices. Prices are measurable functions of the seed for each finite history. The budget holds on every seed and path, not only on average.

Path probabilities are finite products of conditional Bernoulli probabilities. Regret is a finite sum of nonnegative gaps integrated over the seed with lintegral, so it cannot become zero through non-integrability of a real-valued expectation. The exponent of TTT and the scales εℓ\varepsilon_\ellεℓ​ use real division. The arg-min in (51) is any genuine nearest-instance selection; a local proof checks that one exists. The hierarchy's validity and separation statements require a large-horizon threshold, which the unqualified printed prose omits. The corrected Lemma 2 adds pz≤11/6pz\le11/6pz≤11/6 to keep the Bernoulli probabilities in the interior. Its verbatim milestone retains the printed statement, and its formal statement and note disclose the correction. Lemma 3 is not stated because its printed proof relies on Lemma 2 outside that interior region.

The named Bernoulli family strengthens the theorem's existential instance wording, while also exposing a limitation: at p≥2/zp\ge2/zp≥2/z it violates the positive-mean and identifiability conditions stated earlier in the paper. That discrepancy is recorded rather than hidden in an assumption. Definitions for randomized adapted policies, finite path laws, the parameter hierarchy, and Bernoulli relative entropy are reusable; contributions that establish the open goal or repair the information comparison are in scope.

Selected references

  • B. Chen, X. Chao, and Y. Wang, Data-Based Dynamic Pricing and Inventory Control with Censored Demand and Limited Price Changes, SSRN 2700747, revision of 2020-02-10. Preprint.
  • B. Chen, X. Chao, and Y. Wang, Technical Note—Data-Based Dynamic Pricing and Inventory Control with Censored Demand and Limited Price Changes, Operations Research, 2020. DOI 10.1287/opre.2020.1993.
10 thms1 active userReviewed
Bandit AlgorithmsOperations ResearchStatistics·Captain: mikedeng1

Data-Based Dynamic Pricing and Inventory Control with Censored Demand and Limited Price Changes 1: In the Well-Separated Case, Algorithm-I with m Price Changes Has Regret at Most KT^(1/(m+1))Research Paper

Motivation

Retailers set prices and stock levels without knowing how demand responds to price. They must learn the demand curve from their own sales while still earning money, which is the exploration–exploitation trade-off of dynamic pricing with learning (den Boer 2015). Two features of practice make the problem harder than the textbook bandit. First, demand is censored: when the shelf runs empty, the firm sees only its sales min⁡{dt,yt}\min\{d_t,y_t\}min{dt​,yt​}, not the demand dtd_tdt​. Second, firms change prices rarely, because of menu costs and customer reaction, so a policy may change its price at most mmm times over the whole horizon (Cheung, Simchi-Levi and Wang 2017).

Chen, Chao and Wang study a firm facing both features at once, in a parametric model where the demand distribution is known up to a parameter. For the well-separated case they give an algorithm, Algorithm-I, whose regret grows as T1/(m+1)T^{1/(m+1)}T1/(m+1), and they show that no admissible policy does better in general. This mission formalizes the upper bound. A companion mission formalizes the lower bound.

Setting

A firm sells one product over periods t=1,…,Tt=1,\dots,Tt=1,…,T. At the start of period ttt it holds inventory xtx_txt​ (with x1=0x_1=0x1​=0). It chooses a price pt∈P=[pl,ph]p_t\in\mathcal P=[p^l,p^h]pt​∈P=[pl,ph] and an order-up-to level yt∈Y={yl,…,yh}y_t\in\mathcal Y=\{y^l,\dots,y^h\}yt​∈Y={yl,…,yh} with yt≥xty_t\ge x_tyt​≥xt​, and the order arrives at once. Demand dtd_tdt​ then has probability mass function f(⋅;pt,z)f(\cdot;p_t,z)f(⋅;pt​,z) on {dl,dl+1,…,dh}\{d^l,d^l+1,\dots,d^h\}{dl,dl+1,…,dh}, where dh≤+∞d^h\le+\inftydh≤+∞ and the parameter zzz lies in Z=[zl,zh]\mathcal Z=[z^l,z^h]Z=[zl,zh], 0≤zl≤zh0\le z^l\le z^h0≤zl≤zh. Unmet demand is lost, and the firm observes only the sales min⁡{dt,yt}\min\{d_t,y_t\}min{dt​,yt​}. Leftover stock carries over, xt+1=(yt−dt)+x_{t+1}=(y_t-d_t)^+xt+1​=(yt​−dt​)+. Each unit left over costs h≥0h\ge0h≥0 and each lost sale costs b≥0b\ge0b≥0. The single-period expected profit is

G(p,y,z)=p E[D(p,z)]−h E[y−D(p,z)]+−(b+p) E[D(p,z)−y]+.G(p,y,z)=p\,\mathbb E[D(p,z)]-h\,\mathbb E[y-D(p,z)]^+-(b+p)\,\mathbb E[D(p,z)-y]^+ .G(p,y,z)=pE[D(p,z)]−hE[y−D(p,z)]+−(b+p)E[D(p,z)−y]+.

For each yyy, py∗(z)p^*_y(z)py∗​(z) is an optimal price, max⁡p∈PG(p,y,z)=G(py∗(z),y,z)\max_{p\in\mathcal P}G(p,y,z)=G(p^*_y(z),y,z)maxp∈P​G(p,y,z)=G(py∗​(z),y,z), and the clairvoyant value is G∗(z)=max⁡(p,y)∈P×YG(p,y,z)G^*(z)=\max_{(p,y)\in\mathcal P\times\mathcal Y}G(p,y,z)G∗(z)=max(p,y)∈P×Y​G(p,y,z).

A policy's decision in period ttt may depend only on past decisions and past sales. The regret of a policy is

R(T)=T G∗(z)−E[Vϕ(T,z)]=∑t=1TE[G∗(z)−G(pt,yt,z)].R(T)=T\,G^*(z)-\mathbb E[V^\phi(T,z)]=\sum_{t=1}^T\mathbb E\big[G^*(z)-G(p_t,y_t,z)\big].R(T)=TG∗(z)−E[Vϕ(T,z)]=t=1∑T​E[G∗(z)−G(pt​,yt​,z)].

The family is well-separated (Definition 1) if for every price the pmfs f(⋅;p,z)f(\cdot;p,z)f(⋅;p,z), z∈Zz\in\mathcal Zz∈Z, are pairwise distinct.

Algorithm-I splits the horizon into m+1m+1m+1 stages of lengths Ii=⌈Ti/(m+1)⌉I_i=\lceil T^{i/(m+1)}\rceilIi​=⌈Ti/(m+1)⌉ for i≤mi\le mi≤m, with Im+1=T−∑i≤mIiI_{m+1}=T-\sum_{i\le m}I_iIm+1​=T−∑i≤m​Ii​. Stage iii uses one price p^i\hat p_ip^​i​ and the order-up-to target y~i\tilde y_iy~​i​. The target equals the planned level y^i\hat y_iy^​i​, raised by Δ≥1\Delta\ge1Δ≥1 when y^i=dl\hat y_i=d^ly^​i​=dl so that the sales carry information. At the end of the stage the algorithm computes the censored-data maximum-likelihood estimate z^i\hat z_iz^i​ (5). It then sets the next stage's decisions to a maximizer of G(⋅,⋅,z^i)G(\cdot,\cdot,\hat z_i)G(⋅,⋅,z^i​) (6). It changes its price at most mmm times.

Formalization targets

Goal: Theorem 1

Under Assumption A (regularity of GGG and of py∗p^*_ypy∗​ near the true zzz), Assumption 1 (Fisher-information bounds, a positive lower bound on the pmf, strict monotonicity of f(dl;p,⋅)f(d^l;p,\cdot)f(dl;p,⋅)), Definition 1, and the condition that every optimal order-up-to level exceeds dld^ldl, there is K6>0K_6>0K6​>0 such that for all large TTT

R(T)≤K6 T1m+1.R(T)\le K_6\,T^{\frac1{m+1}} .R(T)≤K6​Tm+11​.

The constant may depend on the model, the true parameter and the algorithm's inputs. The rate T1/(m+1)T^{1/(m+1)}T1/(m+1) is the content.

Milestones

  1. Theorem B1: if an estimator satisfies P{∥z^−z∥≥ϵ}≤K14e−nK15ϵ2+K16/n\mathbb P\{\|\hat z-z\|\ge\epsilon\}\le K_{14}e^{-nK_{15}\epsilon^2}+K_{16}/nP{∥z^−z∥≥ϵ}≤K14​e−nK15​ϵ2+K16​/n, then the plug-in decision has optimality gap G∗(z)−E[G(p^∗,y^∗,z)]≤K17/nG^*(z)-\mathbb E[G(\hat p^*,\hat y^*,z)]\le K_{17}/nG∗(z)−E[G(p^​∗,y^​∗,z)]≤K17​/n.
  2. Proposition 1: the censored MLE satisfies P{∣z−z^i∣≥ϵ}≤K3e−K4Iiϵ2+K5/Ii\mathbb P\{|z-\hat z_i|\ge\epsilon\}\le K_3e^{-K_4I_i\epsilon^2}+K_5/I_iP{∣z−z^i​∣≥ϵ}≤K3​e−K4​Ii​ϵ2+K5​/Ii​.
  3. Per-stage estimation regret (p. 25): E[G∗(z)−G(p^i,y^i,z)]≤K17/Ii−1\mathbb E[G^*(z)-G(\hat p_i,\hat y_i,z)]\le K_{17}/I_{i-1}E[G∗(z)−G(p^​i​,y^​i​,z)]≤K17​/Ii−1​.
  4. (25): the summed estimation regret is at most a constant times mT1/(m+1)mT^{1/(m+1)}mT1/(m+1).
  5. (26): P(y^i≠y~i)≤K/Ii−1\mathbb P(\hat y_i\ne\tilde y_i)\le K/I_{i-1}P(y^​i​=y~​i​)≤K/Ii−1​, so exploration after stage 1 is rare.

Significance

The upper bound shows that a budget of mmm price changes costs only a T1/(m+1)T^{1/(m+1)}T1/(m+1) regret even under censoring. With mmm of order log⁡T\log TlogT the rate drops to logarithmic (Theorem 3 of the paper). Together with the matching lower bound (Theorem 2), it identifies the exact price of limiting price changes in this model. Theorem B1 is a general transfer principle, from a large-deviation bound on an estimate to a 1/n1/n1/n optimality gap of the plug-in decision. It allows continuous and discrete decisions together and needs no unique optimizer. It is reusable for any parametric learning problem with that structure. Proposition 1 is a large-deviation bound for maximum likelihood from censored, dependent, non-identically distributed data. The paper notes that such results were not previously available in the statistics literature.

The results are proved in the paper's appendices. None of them has a machine-checked proof. A formalization would check the regret decomposition (23), the stagewise bookkeeping with ceilings, and the large-deviation argument for the censored MLE. It would also make the reading of several loosely stated steps exact (see Formalization scope).

Difficulty

The obvious argument plugs a consistent estimator into the optimizer and bounds the regret by the estimation error. It fails twice. The error of z^i\hat z_iz^i​ is of order Ii−1−1/2I_{i-1}^{-1/2}Ii−1−1/2​, so a Lipschitz bound on the profit would give regret of order Ii/Ii−1I_i/\sqrt{I_{i-1}}Ii​/Ii−1​​ per stage, which is too large. The 1/n1/n1/n rate of Theorem B1 needs the first-order condition at an interior optimum and the exact gap δ\deltaδ between optimal and suboptimal discrete levels. Second, the data are censored and dependent. The censoring level yty_tyt​ depends on carried inventory, which depends on past demand and on past prices. The classical large-deviation results for maximum likelihood assume i.i.d. uncensored samples (Borovkov 1998) and do not apply. Proposition 1 needs a separate argument that the inventory target is reached in most periods with high probability.

Formalization scope

The model is encoded with a scalar parameter. P\mathcal PP, Y\mathcal YY and Z\mathcal ZZ are as above. The pmf is positive exactly on {dl,…,dh}\{d^l,\dots,d^h\}{dl,…,dh} and sums to one, and dh=+∞d^h=+\inftydh=+∞ is allowed. The expected profit GGG is built from these sums. G∗(z)G^*(z)G∗(z) is the finite maximum max⁡yG(py∗(z),y,z)\max_{y}G(p^*_y(z),y,z)maxy​G(py∗​(z),y,z) for a price selection p∗p^*p∗, not a real supremum. Periods are 0-based in Lean, while stages keep the paper's indices. Expectations and probabilities are sums over demand paths d∈NTd\in\mathbb N^Td∈NT weighted by ∏tf(dt;pt,z)\prod_tf(d_t;p_t,z)∏t​f(dt​;pt​,z), the law of the path under the algorithm, taken in [0,∞][0,\infty][0,∞]. Every regret summand is nonnegative.

Algorithm-I reads demand only through sales. Step 2's arg max is any maximizer of the censored product likelihood, which exists by compactness. Step 3's arg max is y^\hat yy^​, any maximizer of y↦G(py∗(z^),y,z^)y\mapsto G(p^*_y(\hat z),y,\hat z)y↦G(py∗​(z^),y,z^), together with the price py^∗(z^)p^*_{\hat y}(\hat z)py^​∗​(z^) of Assumption A(iii). This is the decision the proof of Theorem B1 analyses. Proposition 1's "large enough iii" is read as "large enough TTT, uniformly in the stage and in ϵ\epsilonϵ". Theorem B1 is stated for a scalar price and level and a parameter in Rk\mathbb R^kRk.

Three hypotheses are added and disclosed. Differentiability of fff in zzz is presupposed by Assumption 1. Finiteness of the mean demand is needed for (4) when dh=∞d^h=\inftydh=∞. The condition dl≤yld^l\le y^ldl≤yl is needed because Step 1 defines the target only for y^≥dl\hat y\ge d^ly^​≥dl. "The optimal order-up-to level exceeds dld^ldl" is read for every optimal level.

A trivializing formalization is ruled out. Regret and probabilities are path sums of nonnegative terms, not Bochner integrals that vanish on non-integrable inputs. T1/(m+1)T^{1/(m+1)}T1/(m+1) is a real power. The hypotheses are jointly satisfiable: a Bernoulli instance with Y={1}\mathcal Y=\{1\}Y={1}, P=[1,2]\mathcal P=[1,2]P=[1,2], Z=[2/3,4/5]\mathcal Z=[2/3,4/5]Z=[2/3,4/5], f(0;p,z)=zp/2f(0;p,z)=zp/2f(0;p,z)=zp/2 satisfies all of them, together with explicit selections, and this has been verified in Lean.

A complete development needs the following:

  • Hoeffding-type concentration for dependent bounded sums;
  • a large-deviation bound for the censored log-likelihood;
  • Taylor's theorem with integral remainder on an interval;
  • bookkeeping lemmas for the stage lengths ⌈Ti/(m+1)⌉\lceil T^{i/(m+1)}\rceil⌈Ti/(m+1)⌉.

Theorem B1 and the model definitions are reusable for the paper's later results (Theorems 3, 5, 6), which repeat these arguments. Contributions to any milestone are welcome, including alternative proofs of Proposition 1.

Selected references

  • B. Chen, X. Chao, Y. Wang, Data-Based Dynamic Pricing and Inventory Control with Censored Demand and Limited Price Changes, Operations Research (technical note), 2020; read in the SSRN revision of 2020-02-10. https://doi.org/10.1287/opre.2020.1993, https://ssrn.com/abstract=2700747
  • W. C. Cheung, D. Simchi-Levi, H. Wang, Dynamic Pricing and Demand Learning with Limited Price Experimentation, Operations Research 65(6), 2017. https://doi.org/10.1287/opre.2017.1629
  • J. Broder, P. Rusmevichientong, Dynamic Pricing under a General Parametric Choice Model, Operations Research 60(4), 2012. https://doi.org/10.1287/opre.1120.1057
  • A. A. Borovkov, Mathematical Statistics, Gordon and Breach, 1998.
  • A. V. den Boer, Dynamic pricing and learning: Historical origins, current research, and new directions, Surveys in Operations Research and Management Science 20(1), 2015. https://doi.org/10.1016/j.sorms.2015.03.001
9 thms1 active userReviewed
Operations ResearchOptimization·Captain: mikedeng1

Optimization-Based Scenario Reduction for Data-Driven Two-Stage Stochastic Optimization II: Under a Symmetric Uncertainty Law the Mean Solves the Convex Upper-Bound Scenario ProblemResearch Paper

Why reduce scenarios with the cost in view

A two-stage stochastic programme chooses a decision zzz before an uncertain parameter ξ\xiξ is revealed and pays a cost c(z;ξ)c(z;\xi)c(z;ξ). When ξ\xiξ is known only through nnn data points, or through a law with many atoms, the expected cost is expensive to optimize, and practitioners replace the law by a small set of mmm reduced scenarios. Classical scenario reduction chooses these scenarios to be close to the data in a Wasserstein distance, without looking at the cost (Dupačová, Gröwe-Kuska and Römisch, 2003; Rujeerapaiboon, Schindler, Kuhn and Wiesemann, 2022). Bertsimas and Mundru (Operations Research, 2023) propose problem-dependent scenario reduction: the distance between two scenarios measures how much the optimal decision for one costs when evaluated on the other. Their method needs, for each cluster of data, the scenario that minimizes a cost-dependent dissimilarity, which they replace by a convex upper bound, problem (16).

Theorem 3 of the paper is the sanity check for that upper bound. For the simplest cost to which it applies and a single reduced scenario (m=1m=1m=1), under a symmetry assumption on the law, the minimizer of the upper bound is the mean of ξ\xiξ, which is also what Wasserstein scenario reduction returns for m=1m=1m=1. This mission formalizes that statement.

Setting

Decisions and parameters both live in Rd\mathbb R^dRd, with the Euclidean inner product z′ξz'\xiz′ξ. The feasible set is a polytope

Z={z∈R+d:Pz≤q},\mathcal Z = \{ z\in\mathbb R^d_+ : Pz\le q\},Z={z∈R+d​:Pz≤q},

for a matrix P∈Rr×dP\in\mathbb R^{r\times d}P∈Rr×d and a vector q∈Rrq\in\mathbb R^rq∈Rr, assumed nonempty and bounded. The cost is c(z;ξ)=max⁡{z′ξ,0}c(z;\xi)=\max\{z'\xi,0\}c(z;ξ)=max{z′ξ,0}, a special case of the piecewise bilinear costs max⁡1≤t≤kz′Atξ\max_{1\le t\le k} z'A_t\ximax1≤t≤k​z′At​ξ of (15) (take k=2k=2k=2, A1=IA_1=IA1​=I, A2=0A_2=0A2​=0).

An optimal-decision map is any function z∗:Rd→Rdz^*:\mathbb R^d\to\mathbb R^dz∗:Rd→Rd with z∗(ξ)∈arg⁡min⁡z∈Zz′ξz^*(\xi)\in\arg\min_{z\in\mathcal Z} z'\xiz∗(ξ)∈argminz∈Z​z′ξ for every ξ\xiξ; for a vector vvv, Z∗(v)=arg⁡min⁡z∈Zz′v\mathcal Z^*(v)=\arg\min_{z\in\mathcal Z} z'vZ∗(v)=argminz∈Z​z′v is the whole set of minimizers. The set of cost-reducing parameters is

U={ξ:min⁡z∈Zz′ξ<0}.\mathcal U=\Big\{\xi : \min_{z\in\mathcal Z} z'\xi<0\Big\}.U={ξ:z∈Zmin​z′ξ<0}.

For a candidate reduced scenario ζ∈Rd\zeta\in\mathbb R^dζ∈Rd, the per-scenario upper-bound loss is

L(ξ,ζ)=max⁡{max⁡z∈Zz′(ξ−2ζ), 0}+2max⁡{z∗(ξ)′ζ, 0}.L(\xi,\zeta)=\max\Big\{\max_{z\in\mathcal Z} z'(\xi-2\zeta),\,0\Big\}+2\max\{z^*(\xi)'\zeta,\,0\}.L(ξ,ζ)=max{z∈Zmax​z′(ξ−2ζ),0}+2max{z∗(ξ)′ζ,0}.

For this cost it is the optimal value, in the auxiliary variables (λ,θ,γ)(\lambda,\theta,\gamma)(λ,θ,γ), of the per-point linear programme in (16). With a random parameter ξ\xiξ of law μ\muμ, the upper-bound objective for m=1m=1m=1 is E[L(ξ,ζ)]\mathbb E[L(\xi,\zeta)]E[L(ξ,ζ)]. The squared cost for Wasserstein reduction to one atom is W(ζ)=∫∥ξ−ζ∥22 dμW(\zeta)=\int\|\xi-\zeta\|_2^2\,d\muW(ζ)=∫∥ξ−ζ∥22​dμ; (4) takes its square root, which preserves minimizers.

Assumption 5 (p. 24) asks:

  • (a) ξ\xiξ has a continuous distribution on U\mathcal UU and density 000 elsewhere;
  • (b) the mean ξˉ=E[ξ]\bar\xi=\mathbb E[\xi]ξˉ​=E[ξ] is finite and z∗(ξˉ)′ξˉ>0z^*(\bar\xi)'\bar\xi>0z∗(ξˉ​)′ξˉ​>0;
  • (c) the law of ξ\xiξ is symmetric about ξˉ\bar\xiξˉ​, that is, ξ\xiξ and 2ξˉ−ξ2\bar\xi-\xi2ξˉ​−ξ are equal in distribution.

Formalization targets

Goal: Theorem 3 (p. 24)

Under Assumption 5, for every measurable optimal-decision map z∗z^*z∗,

E[L(ξ,ξˉ)]  ≤  E[L(ξ,ζ)]andW(ξˉ)≤W(ζ)for all ζ∈Rd.\mathbb E\big[L(\xi,\bar\xi)\big]\;\le\;\mathbb E\big[L(\xi,\zeta)\big]\quad\text{and}\quad W(\bar\xi)\le W(\zeta)\qquad\text{for all }\zeta\in\mathbb R^d.E[L(ξ,ξˉ​)]≤E[L(ξ,ζ)]andW(ξˉ​)≤W(ζ)for all ζ∈Rd.

The page's "the solution … is simply the mean" is read as "the mean is a minimizer", which is what the proof concludes; uniqueness is not claimed.

Milestones (steps of the proof, pp. 24–25)

  1. For each ξ\xiξ, ζ↦L(ξ,ζ)\zeta\mapsto L(\xi,\zeta)ζ↦L(ξ,ζ) is convex.
  2. L(ξ,ζ)≤∣max⁡z∈Zz′(ξ−2ζ)∣+2∣ζ′z∗(ξ)∣≤(∥ξ−2ζ∥1+2∥ζ∥1)max⁡z∈Z∥z∥∞≤(∥ξ∥1+4∥ζ∥1)max⁡z∈Z∥z∥∞L(\xi,\zeta)\le\big|\max_{z\in\mathcal Z} z'(\xi-2\zeta)\big|+2|\zeta'z^*(\xi)|\le(\|\xi-2\zeta\|_1+2\|\zeta\|_1)\max_{z\in\mathcal Z}\|z\|_\infty\le(\|\xi\|_1+4\|\zeta\|_1)\max_{z\in\mathcal Z}\|z\|_\inftyL(ξ,ζ)≤​maxz∈Z​z′(ξ−2ζ)​+2∣ζ′z∗(ξ)∣≤(∥ξ−2ζ∥1​+2∥ζ∥1​)maxz∈Z​∥z∥∞​≤(∥ξ∥1​+4∥ζ∥1​)maxz∈Z​∥z∥∞​.
  3. With probability one, 2ξˉ−ξ∈U2\bar\xi-\xi\in\mathcal U2ξˉ​−ξ∈U and Z∗(2ξˉ−ξ)\mathcal Z^*(2\bar\xi-\xi)Z∗(2ξˉ​−ξ) is a singleton.
  4. Under symmetry, E[z∗(2ξˉ−ξ)]=E[z∗(ξ)]\mathbb E[z^*(2\bar\xi-\xi)]=\mathbb E[z^*(\xi)]E[z∗(2ξˉ​−ξ)]=E[z∗(ξ)].
  5. If z∗(ξˉ)′ξˉ>0z^*(\bar\xi)'\bar\xi>0z∗(ξˉ​)′ξˉ​>0 then z′ξˉ≥z∗(ξˉ)′ξˉ>0z'\bar\xi\ge z^*(\bar\xi)'\bar\xi>0z′ξˉ​≥z∗(ξˉ​)′ξˉ​>0 for every z∈Zz\in\mathcal Zz∈Z.

Significance

The upper bound (16) is the computational core of the paper's Algorithm 1: each alternating-minimization step solves one instance of it per cluster. Theorem 3 is the only result in the paper that identifies its minimizer in closed form, and it shows that, in a symmetric case, the cost-aware reduction agrees with the classical one at m=1m=1m=1. It is the paper's justification that the upper bound is a sensible surrogate and not an arbitrary convex relaxation.

The argument is a variant of the mean-minimizes-the-surrogate result of Elmachtoub and Grigas for the SPO+ loss (Management Science, 2022, Proposition 6(a)), with both terms clipped at zero and with a law living only on U\mathcal UU instead of all of Rd\mathbb R^dRd. Theorem 3 is proved on paper; it has not been formalized, and on Prove2Me the Elmachtoub–Grigas statement is posed but open. The work here is to formalize the known proof, including the measure-theoretic steps the paper treats in one sentence each: integrability of the loss, exchange of the subdifferential and the expectation, and the almost-sure uniqueness of linear minimizers over a polytope.

Difficulty

Convexity of the objective and a vanishing subgradient at ξˉ\bar\xiξˉ​ would settle the theorem, but two steps of that argument are not routine. First, the objective is not differentiable everywhere, and the subdifferential of ζ↦E[L(ξ,ζ)]\zeta\mapsto\mathbb E[L(\xi,\zeta)]ζ↦E[L(ξ,ζ)] is not in general the expectation of the pointwise subdifferentials: the exchange needs the integrable bound of milestone 2 and the almost-sure uniqueness of milestone 3. Second, that uniqueness is a statement about the geometry of polytopes and about the law of the reflected parameter 2ξˉ−ξ2\bar\xi-\xi2ξˉ​−ξ together: neither convexity nor absolute continuity alone gives it. The clipping at zero adds a third point: one has to check that it is inactive at ζ=ξˉ\zeta=\bar\xiζ=ξˉ​ almost surely, which is where Assumption 5a's support condition on U\mathcal UU and Assumption 5b enter. Dropping either half of Assumption 5a breaks the argument: a two-atom symmetric law can put both atoms where Z∗(2ξˉ−ξ)\mathcal Z^*(2\bar\xi-\xi)Z∗(2ξˉ​−ξ) is a whole edge, and a law charging the complement of U\mathcal UU can make the first clipping active at ξˉ\bar\xiξˉ​.

Formalization scope

  • Rd\mathbb R^dRd is EuclideanSpace ℝ (Fin d) and z′ξz'\xiz′ξ is ⟪z, ξ⟫_ℝ. Decision and parameter spaces coincide (nz=dn_z=dnz​=d), as the cost z′ξz'\xiz′ξ requires.
  • Z\mathcal ZZ is polytope P qv, with the hypotheses that it is nonempty and bounded (§1 of the paper assumes Z\mathcal ZZ nonempty and compact; the polytope is closed).
  • The maximum and minimum of z′vz'vz′v over Z\mathcal ZZ and the argmin set are the published definitions SmartPTO.Fisher.xi, SmartPTO.Fisher.zstar and SmartPTO.Fisher.Wstar; "symmetric about its mean" is the published SmartPTO.Fisher.CentrallySymmetric.
  • z∗z^*z∗ is a function zsel with the hypothesis that zsel ξ is a minimizer for every ξ\xiξ, and it is assumed measurable (needed for E[z∗(ξ)]\mathbb E[z^*(\xi)]E[z∗(ξ)]; not on the page). The theorem holds for every such selection; uniqueness of minimizers is not assumed.
  • The law is a probability measure μ\muμ. Assumption 5a is two hypotheses: μ≪\mu\llμ≪ Lebesgue and μ(Rd∖U)=0\mu(\mathbb R^d\setminus\mathcal U)=0μ(Rd∖U)=0. Full support is not assumed and would contradict the second. Assumption 5b is integrability of ξ\xiξ plus z∗(ξˉ)′ξˉ>0z^*(\bar\xi)'\bar\xi>0z∗(ξˉ​)′ξˉ​>0 with ξˉ=∫ξ dμ\bar\xi=\int\xi\,d\muξˉ​=∫ξdμ.
  • The goal is stated on E[L(ξ,ζ)]\mathbb E[L(\xi,\zeta)]E[L(ξ,ζ)], the first display of the proof, not on the linear programme (16) with its auxiliary variables; for this cost the two coincide by LP duality (Proposition 3's proof), and Theorem 3 is the population, m=1m=1m=1 version.
  • Expectations are Bochner integrals, which return 000 for non-integrable integrands. Under the stated hypotheses ξ↦L(ξ,ζ)\xi\mapsto L(\xi,\zeta)ξ↦L(ξ,ζ) is dominated by (∥ξ∥1+4∥ζ∥1)max⁡z∈Z∥z∥∞(\|\xi\|_1+4\|\zeta\|_1)\max_{z\in\mathcal Z}\|z\|_\infty(∥ξ∥1​+4∥ζ∥1​)maxz∈Z​∥z∥∞​ and is integrable, so the goal is not satisfied by junk values; integrability and measurability hypotheses must not be weakened.
  • The Wasserstein clause uses an extended nonnegative integral for WWW. This expresses the theorem without adding a second-moment hypothesis absent from Assumption 5. With a finite second moment it is the usual squared Wasserstein objective; with an infinite second moment every one-atom value is infinite, so the clause has only the extended-value meaning.

Useful infrastructure, reusable beyond this mission: support functions of polytopes and their directional derivatives; the fact that a linear function has a unique minimizer over a polytope for Lebesgue-almost every cost vector; differentiation of convex integral functionals under the integral sign. Contributions of any of these as stand-alone lemmas are welcome.

The source is the authors' accepted manuscript (MIT DSpace, "Submitted to Operations Research"), not the published article; all page and theorem numbers refer to the manuscript, whose printed page numbers equal its PDF page numbers.

Selected references

  • D. Bertsimas and N. Mundru, Optimization-based Scenario Reduction for Data-Driven Two-stage Stochastic Optimization, Operations Research, 2023 (author manuscript, MIT DSpace). https://doi.org/10.1287/opre.2022.2265
  • A. N. Elmachtoub and P. Grigas, Smart "Predict, then Optimize", Management Science 68(1), 2022. https://doi.org/10.1287/mnsc.2020.3922
  • J. Dupačová, N. Gröwe-Kuska and W. Römisch, Scenario reduction in stochastic programming: an approach using probability metrics, Mathematical Programming 95, 2003. https://doi.org/10.1007/s10107-002-0331-0
  • N. Rujeerapaiboon, K. Schindler, D. Kuhn and W. Wiesemann, Scenario Reduction Revisited: Fundamental Limits and Guarantees, Mathematical Programming 191, 2022. https://doi.org/10.1007/s10107-018-1269-1
9 thms1 active userReviewed
Operations Research·Captain: mikedeng1

Robust Dynamic Pricing with Strategic Customers: Over Horizon T the Simple Robust Pricing Policy Earns at Least V*_β(x₀)/(1 + e^{−βT}/(βT))Research Paper

Motivation

A seller with a fixed stock must choose prices before knowing when customers will arrive or what they will be willing to pay. If customers can wait for a later price, the seller must also account for how today's price changes their purchasing decisions. Chen and Farias study this problem with strategic customers and identify a simple inventory-based policy that can be compared with the best attainable revenue. Their Theorem 1 shows that this policy, with discount rate β=1/(1.42T)\beta=1/(1.42T)β=1/(1.42T) for a season of length TTT, earns at least 0.290.290.29 of the optimal expected revenue. The policy posts the price that is optimal in the canonical discounted revenue-management problem with myopic customers (Gallego and van Ryzin, 1994; Farias and Van Roy, 2010) at the current inventory level. This mission targets the lower-bound half of its argument, stated as Lemma 9: a comparison between the policy's revenue over a fixed selling season and a discounted, infinite-horizon value. The latter comparison is meaningful independently of the paper's strategic-customer benchmark.

Setting

The seller begins with inventory x0∈Nx_0\in\mathbb Nx0​∈N and sells at most one item to each customer. Customers arrive at rate λ>0\lambda>0λ>0. Their nonnegative valuations have density fff with total mass one. Write Fˉ(p)=∫p∞f(v) dv\bar F(p)=\int_p^\infty f(v)\,dvFˉ(p)=∫p∞​f(v)dv for the probability that a valuation exceeds price ppp. A customer who behaves myopically buys at a posted price ppp exactly when the valuation is high enough, so sales occur at rate λFˉ(p)\lambda\bar F(p)λFˉ(p) while that price is posted.

The paper's Assumption 1 concerns the virtual value ψ(p)=p−Fˉ(p)/f(p)\psi(p)=p-\bar F(p)/f(p)ψ(p)=p−Fˉ(p)/f(p): it is nondecreasing on nonnegative prices and has a nonnegative root. The simple policy is built from a discount rate β>0\beta>0β>0 and an infinite-horizon Bellman value V:N→RV:\mathbb N\to\mathbb RV:N→R. It has V(0)=0V(0)=0V(0)=0 and, for x≥1x\ge1x≥1, satisfies the sign-corrected recursion

βV(x)=sup⁡p≥0λFˉ(p)(p+V(x−1)−V(x)).\beta V(x)=\sup_{p\ge0}\lambda\bar F(p)\bigl(p+V(x-1)-V(x)\bigr).βV(x)=p≥0sup​λFˉ(p)(p+V(x−1)−V(x)).

At inventory x≥1x\ge1x≥1 the root price π(x)\pi(x)π(x) is the unique nonnegative solution of ψ(π(x))=V(x)−V(x−1)\psi(\pi(x))=V(x)-V(x-1)ψ(π(x))=V(x)−V(x−1). The price posted just before a sale is π(x)\pi(x)π(x), and the inventory then becomes x−1x-1x−1. At zero inventory the posted price is infinity and there are no more sales. The expected sum of prices paid at sales by time TTT is denoted J(x0,T)J(x_0,T)J(x0​,T). These objects come from §§2–3.1 of the paper.

Formalization targets

The goal is Lemma 9, for every λ,β,T>0\lambda,\beta,T>0λ,β,T>0 and every initial inventory x0x_0x0​:

V(x0)≤(1+e−βTβT)J(x0,T).V(x_0)\le\left(1+\frac{e^{-\beta T}}{\beta T}\right)J(x_0,T).V(x0​)≤(1+βTe−βT​)J(x0​,T).

The coefficient is part of the target. No choice of β\betaβ is fixed. The milestones follow the results stated in the paper: inventory monotonicity of VVV, price monotonicity, pathwise monotonicity of the price and its revenue rate (Lemma 7), the expected-rate identity, revenue-per-time monotonicity (Lemma 8), the identification of VVV with discounted policy revenue, and the comparison through an independent exponential horizon. The last milestone evaluates E[max⁡{1,X/T}]\mathbb E[\max\{1,X/T\}]E[max{1,X/T}] exactly for X∼Exp⁡(β)X\sim\operatorname{Exp}(\beta)X∼Exp(β).

Significance

The bound quantifies how much of the discounted value this explicit policy earns within a finite selling season. Its factor depends only on the dimensionless product βT\beta TβT; it does not depend on inventory, arrival rate or valuation density. Together with the paper's separate upper bound on its strategic benchmark, Lemma 9 contributes to the advertised 0.290.290.29 comparison. The full comparison requires definitions and equilibrium arguments beyond this mission, so the goal here is precisely the fixed-horizon lower bound.

Chen and Farias proved these statements in the published paper. The Lean items in this proposal are statements awaiting machine-checked proofs. A completed development would establish the discounted-value verification step and the fixed-horizon comparison for an explicit pure-death sales process. Its constructions of inventory-dependent sale times and finite-horizon revenue could also be reused in other pricing models with Poisson arrivals.

Difficulty

The Bellman value is specified by an optimization equation, while J(x0,T)J(x_0,T)J(x0​,T) is defined from the actual sales process. Their equality at an exponential horizon does not follow by unfolding a definition. The price chosen through the virtual-value equation must be connected to the Bellman optimum, and the inventory-dependent process must be related to the same value. The paper also relies on the fact that the root price falls as inventory rises, citing Lemma 1 of Farias and Van Roy; that fact does not follow merely from monotonicity of ψ\psiψ without control of the marginal values V(x)−V(x−1)V(x)-V(x-1)V(x)−V(x−1). Lemma 8 then concerns expected revenue over two different horizons, including paths that sell out early.

Formalization scope

The Lean model uses a nonnegative integrable density on R\mathbb RR, zero on negative valuations and normalized to mass one. It adds f(p)>0f(p)>0f(p)>0 for every p≥0p\ge0p≥0. This disclosed condition excludes bounded-support densities: without it the paper's quotient Fˉ(p)/f(p)\bar F(p)/f(p)Fˉ(p)/f(p) can divide by zero and Lean would assign a spurious value. The paper says the root equation has a unique solution under Assumption 1, but nondecreasing ψ\psiψ alone does not guarantee existence or uniqueness. Accordingly, the root policy is characterized by its equation and an explicit uniqueness condition. The Bellman recursion uses a least-upper-bound predicate rather than a real supremum with a default value on an empty or unbounded set. The corrected sign in the recursion matches the paper's root equation and discounted-revenue equation.

For stock x0x_0x0​, the model uses x0x_0x0​ independent rate-one exponential clocks. At inventory xxx, the next holding time is the next clock divided by λFˉ(π(x))\lambda\bar F(\pi(x))λFˉ(π(x)). Sales are numbered 1,…,x01,\ldots,x_01,…,x0​, and sale kkk earns π(x0−k+1)\pi(x_0-k+1)π(x0​−k+1). The law is the myopic sales process with intensity λFˉ(π(x))\lambda\bar F(\pi(x))λFˉ(π(x)). Expected revenues are lower integrals in the extended nonnegative reals; nonnegative prices make their real-to-extended conversion exact. An independent X∼Exp⁡(β)X\sim\operatorname{Exp}(\beta)X∼Exp(β) is represented by an outer integral over its law. The factor λ\lambdaλ missing from the paper's printed compensator display is restored, and the paper's intermediate 1+e−βT/T1+e^{-\beta T}/T1+e−βT/T is read with the missing β\betaβ restored.

The strategic-customer types, stopping rules, perfect Bayesian equilibrium, and the identification of the strategic policy's revenue with myopic-sales revenue are not formalized here. Neither are Theorem 1's OPT\mathrm{OPT}OPT and J∗J^*J∗ benchmarks: their specification is not pinned down sufficiently in the paper for this proposal. In particular, JJJ is constructed from sale times and payments, never from VVV, and VVV is specified by the Bellman equation, never as the policy's revenue. Contributions toward the goal include the individual milestones and reusable facts about the clock process, its revenue, and the exponential horizon.

Selected references

  • Y. Chen and V. F. Farias, Robust Dynamic Pricing with Strategic Customers, Mathematics of Operations Research 43(4):1119–1142, 2018. DOI: 10.1287/moor.2017.0897.
  • G. Gallego and G. van Ryzin, Optimal Dynamic Pricing of Inventories with Stochastic Demand over Finite Horizons, Management Science 40(8):999–1020, 1994. DOI: 10.1287/mnsc.40.8.999.
  • V. F. Farias and B. Van Roy, Dynamic Pricing with a Prior on Market Response, Operations Research 58(1):16–29, 2010. DOI: 10.1287/opre.1090.0729.
  • P. Brémaud, Point Processes and Queues: Martingale Dynamics, Springer, 1981. DOI: 10.1007/978-1-4684-9477-8.
11 thms1 active userReviewed
Algorithmic Game TheoryOperations Research·Captain: mikedeng1

Prophet Inequalities Made Easy: Stochastic Optimization by Pricing Nonstochastic Inputs I: Expected (α, β)-Balanced Prices Scaled by α/(1 + αβ) Earn 1/(1 + αβ) of the Expected Welfare of ALGResearch Paper

Motivation

A prophet inequality compares an online decision maker, who sees random values one at a time and must decide irrevocably, with a "prophet" who sees all values in advance. The classic result of Krengel, Sucheston and Garling (1977) and Samuel-Cahn (1984) says that for a single item and independent values a single threshold earns at least half of the expected maximum. In mechanism design the same inequalities describe posted-price mechanisms: a seller approaches buyers in sequence and offers each a menu of prices; each buyer takes the outcome it likes best. Such mechanisms are truthful, simple and online, so a prophet inequality for a feasibility constraint is a welfare guarantee for a practical selling procedure.

Before this paper, prophet inequalities were proved one constraint at a time (matroids, matroid intersections, polymatroids, unit-demand and combinatorial auctions), each with its own probabilistic argument. Dütting, Feldman, Kesselheim and Lucier (SIAM J. Comput. 49(3), 2020; conference version FOCS 2017) showed that a single extension theorem reduces all of them to a deterministic question about full-information instances.

Setting

There are nnn agents N={1,…,n}N=\{1,\dots,n\}N={1,…,n}. Agent iii has an outcome space XiX_iXi​ containing a null outcome ∅\varnothing∅. An outcome profile is x=(x1,…,xn)\mathbf x=(x_1,\dots,x_n)x=(x1​,…,xn​); x[i−1]\mathbf x_{[i-1]}x[i−1]​ is x\mathbf xx with the outcomes of agents i,…,ni,\dots,ni,…,n replaced by ∅\varnothing∅. The feasible profiles form a downward-closed set F\mathcal FF: removing outcomes keeps a profile feasible.

Agent iii's valuation vi:Xi→[0,1]v_i:X_i\to[0,1]vi​:Xi​→[0,1] is drawn independently from a known distribution Di\mathcal D_iDi​; D=∏iDi\mathcal D=\prod_i\mathcal D_iD=∏i​Di​. The welfare of x\mathbf xx is v(x)=∑ivi(xi)\mathbf v(\mathbf x)=\sum_iv_i(x_i)v(x)=∑i​vi​(xi​), and v(OPT(v,S))\mathbf v(\mathrm{OPT}(\mathbf v,S))v(OPT(v,S)) is the largest welfare over a set SSS of profiles. An outcome rule ALG\mathrm{ALG}ALG maps each valuation profile to a feasible profile.

A pricing rule assigns a price pi(xi∣y)∈[0,∞]p_i(x_i\mid\mathbf y)\in[0,\infty]pi​(xi​∣y)∈[0,∞] to outcome xix_ixi​ offered to agent iii when the partial allocation is y\mathbf yy, with price ∞\infty∞ when xix_ixi​ cannot be feasibly added to y\mathbf yy. In the posted-price mechanism the agents are approached in index order; agent iii sees the prices given the allocation so far and takes an outcome maximizing its utility vi(xi)−pi(xi∣y)v_i(x_i)-p_i(x_i\mid\mathbf y)vi​(xi​)−pi​(xi​∣y).

A set H\mathcal HH of profiles is exchange compatible with x\mathbf xx if (yi,x−i)∈F(y_i,\mathbf x_{-i})\in\mathcal F(yi​,x−i​)∈F for every y∈H\mathbf y\in\mathcal Hy∈H and every iii. For a fixed valuation profile v\mathbf vv, a pricing rule pvp^{\mathbf v}pv is (α,β)(\alpha,\beta)(α,β)-balanced (Definition 3.1) with respect to ALG\mathrm{ALG}ALG and an exchange-compatible family (Fx)x(\mathcal F_{\mathbf x})_{\mathbf x}(Fx​)x​ if for every x∈F\mathbf x\in\mathcal Fx∈F

(a)∑ipiv(xi∣x[i−1]) ≥ 1α(v(ALG(v))−v(OPT(v,Fx))),\text{(a)}\quad\sum_{i}p^{\mathbf v}_i(x_i\mid\mathbf x_{[i-1]})\ \ge\ \tfrac1\alpha\bigl(\mathbf v(\mathrm{ALG}(\mathbf v))-\mathbf v(\mathrm{OPT}(\mathbf v,\mathcal F_{\mathbf x}))\bigr),(a)i∑​piv​(xi​∣x[i−1]​) ≥ α1​(v(ALG(v))−v(OPT(v,Fx​))), (b)∑ipiv(xi′∣x[i−1]) ≤ β⋅v(OPT(v,Fx))for all x′∈Fx.\text{(b)}\quad\sum_{i}p^{\mathbf v}_i(x'_i\mid\mathbf x_{[i-1]})\ \le\ \beta\cdot\mathbf v(\mathrm{OPT}(\mathbf v,\mathcal F_{\mathbf x}))\quad\text{for all }\mathbf x'\in\mathcal F_{\mathbf x}.(b)i∑​piv​(xi′​∣x[i−1]​) ≤ β⋅v(OPT(v,Fx​))for all x′∈Fx​.

Property (a) says the prices paid by x\mathbf xx cover the welfare x\mathbf xx blocks; (b) says that what remains available is not overpriced.

Formalization targets

Goal: Theorem 3.2 (p. 549)

If every pv~p^{\tilde{\mathbf v}}pv~ is (α,β)(\alpha,\beta)(α,β)-balanced with respect to one exchange-compatible family, post δp\delta pδp with pi(xi∣y)=Ev~[piv~(xi∣y)]p_i(x_i\mid\mathbf y)=\mathbb E_{\tilde{\mathbf v}}[p^{\tilde{\mathbf v}}_i(x_i\mid\mathbf y)]pi​(xi​∣y)=Ev~​[piv~​(xi​∣y)] and δ=α/(1+αβ)\delta=\alpha/(1+\alpha\beta)δ=α/(1+αβ). Then the welfare x(v)\mathbf x(\mathbf v)x(v) of the mechanism satisfies

Ev[v(x(v))] ≥ 11+αβ Ev[v(ALG(v))].\mathbb E_{\mathbf v}\bigl[\mathbf v(\mathbf x(\mathbf v))\bigr]\ \ge\ \frac1{1+\alpha\beta}\,\mathbb E_{\mathbf v}\bigl[\mathbf v(\mathrm{ALG}(\mathbf v))\bigr].Ev​[v(x(v))] ≥ 1+αβ1​Ev​[v(ALG(v))].

Milestones

The steps of the paper's proof (pp. 549–550), in order: the prefix x[i−1](v)\mathbf x_{[i-1]}(\mathbf v)x[i−1]​(v) does not depend on viv_ivi​; the per-agent utility bound against a deviation drawn from an independent sample v′\mathbf v'v′; (3.1); the bound (3.2) from property (b); (3.3); the pointwise revenue bound from property (a); (3.4); and welfare = utilities + revenue. Two applications complete the list: Claim 3.3 (the single-item price max⁡ℓvℓ\max_\ell v_\ellmaxℓ​vℓ​ is (1,1)(1,1)(1,1)-balanced, which with the goal gives the classic factor 2) and the balancedness half of Theorem 4.2 for knapsack.

Significance

The theorem turns prophet inequalities into a deterministic design problem: to obtain a 1/(1+αβ)1/(1+\alpha\beta)1/(1+αβ) guarantee against any benchmark ALG\mathrm{ALG}ALG, it suffices to exhibit balanced prices for each fixed valuation profile. The paper uses it to rederive the single-item factor 2, to obtain O(d)O(d)O(d) guarantees for combinatorial auctions with bundles of size ddd, a constant for knapsack, and bounds for MPH-kkk valuations, sparse packing problems and multidimensional matroids. Because ALG\mathrm{ALG}ALG can be any (for example polynomial-time) algorithm, the guarantee is relative to what can be computed, not only to the optimum.

The result is proved in the paper. As far as is known, no part of it has a machine-checked proof; the platform holds no posted-price framework of this generality. A formal proof would provide a reusable model of sequential posted-price mechanisms over arbitrary outcome spaces and a checked reduction that later formalizations of the applications can invoke instead of repeating a probabilistic argument.

Difficulty

The obvious argument fixes the realized valuations and charges each agent's purchase against the welfare it blocks. That fails because the prices are averages over a hypothetical profile v~\tilde{\mathbf v}v~, not the prices of the realized instance, and the set of outcomes still available depends on everyone's random values. The proof needs an independent copy v′\mathbf v'v′ of the valuations and the fact that the allocation before agent iii does not depend on viv_ivi​, so that viv_ivi​ and vi′v'_ivi′​ can be exchanged without changing the distribution. Formally, the run of the mechanism is a recursion over agents, the exchange is a measure-preserving swap of coordinates of a product measure, the optimal value v(OPT(v,Fx))\mathbf v(\mathrm{OPT}(\mathbf v,\mathcal F_{\mathbf x}))v(OPT(v,Fx​)) is a supremum that need not be attained, and the family Fx\mathcal F_{\mathbf x}Fx​ may be empty.

Formalization scope

Agents are Fin n, 0-based; the index order is the order of Fin n. Outcomes are (X : Fin n → Type*) with a null profile nul. Valuation types V i are countable with the discrete σ-algebra (a disclosed restriction: the paper allows any distribution); this makes every function measurable and every bounded one integrable. A type w : V i values outcomes through val i w : X i → ℝ, with values in [0,1][0,1][0,1] as on p. 547. The distribution is Measure.pi μ with every μ i a probability measure.

Prices are ℝ≥0∞-valued, ⊤ standing for ∞\infty∞; the expected price is a Lebesgue integral in [0,∞][0,\infty][0,∞]. Property (a) passes its real right side through ENNReal.ofReal, and no price sum is ever converted to a real number, so an infinite price is never read as 000. v(OPT(v,S))\mathbf v(\mathrm{OPT}(\mathbf v,S))v(OPT(v,S)) is the real supremum sSup of the welfare over SSS, equal to 000 on S=∅S=\varnothingS=∅ as in the proof of Claim 3.3. Exchange compatibility of the family is required at every profile, not only feasible ones. The pricing-rule condition "∞\infty∞ when infeasible" is imposed at feasible partial allocations, which is where the paper defines prices.

Two standing assumptions are disclosed. First, the null outcome is free: piv~(∅∣y)=0p^{\tilde{\mathbf v}}_i(\varnothing\mid\mathbf y)=0piv~​(∅∣y)=0 for feasible y\mathbf yy. Every pricing rule of the paper satisfies it, and without it the proof's deviation does not exist when Fx(v)=∅\mathcal F_{\mathbf x(\mathbf v)}=\varnothingFx(v)​=∅. Second, the agents' behaviour is a choice rule that sees only the agent's own type and the allocation so far, assumed utility-maximizing (finite price, best utility among finitely priced outcomes) at every feasible partial allocation. The constant δ=α/(1+αβ)\delta=\alpha/(1+\alpha\beta)δ=α/(1+αβ) is written into the posted prices, not passed as a parameter.

A trivializing formalization is ruled out: the goal mentions only the run of the mechanism, the expected prices scaled by δ\deltaδ, balancedness and the two expectations. It does not assume the utility or revenue bounds, it does not let an agent's choice see other agents' types, and a supplementary check confirms that its hypotheses are jointly satisfiable. Milestones that involve the proof's deviation OPT(v′,Fx(v))\mathrm{OPT}(\mathbf v',\mathcal F_{\mathbf x(\mathbf v)})OPT(v′,Fx(v)​) are stated for an arbitrary selector from Fx\mathcal F_{\mathbf x}Fx​, a generalization that avoids assuming that the supremum is attained. Theorem 4.2 is stated as (2,1)(2,1)(2,1)-balanced, which is what its proof establishes; the printed "(1,2)(1,2)(1,2)" is false.

Useful infrastructure, reusable beyond this mission: measure-preserving coordinate swaps for Measure.pi and products of product measures, ε-maximizer selection for bounded suprema, and induction principles for the sequential run. Contributions of proofs of individual milestones are welcome.

Selected references

  • P. Dütting, M. Feldman, T. Kesselheim, B. Lucier, Prophet inequalities made easy: Stochastic optimization by pricing nonstochastic inputs, SIAM J. Comput. 49(3):540–582, 2020. https://doi.org/10.1137/20M1323850
  • U. Krengel, L. Sucheston, Semiamarts and finite values, Bull. Amer. Math. Soc. 83(4):745–747, 1977. https://doi.org/10.1090/S0002-9904-1977-14378-4
  • E. Samuel-Cahn, Comparison of threshold stop rules and maximum for independent nonnegative random variables, Ann. Probab. 12(4):1213–1216, 1984. https://doi.org/10.1214/aop/1176993150
  • R. Kleinberg, S. M. Weinberg, Matroid prophet inequalities and applications to multi-dimensional mechanism design, Games Econom. Behav. 113:97–115, 2019. https://doi.org/10.1016/j.geb.2014.11.002
  • S. Chawla, J. D. Hartline, D. Malec, B. Sivan, Multi-parameter mechanism design and sequential posted pricing, STOC 2010, 311–320. https://doi.org/10.1145/1806689.1806733
13 thms1 active userReviewed
Bandit AlgorithmsMachine LearningStatistics·Captain: mikedeng1

The Multi-Armed Bandit Problem with Covariates I: Successive Elimination Has Expected Regret at Most 392γ²(1 + n/T)(K/Δ)log(TΔ²/18γ²) + nΔ⁻ Under a Random HorizonResearch Paper

Motivation

In a multi-armed bandit, a decision maker repeatedly chooses an arm, observes a reward, and tries to learn which arm has the highest average reward while losing little reward during learning. A policy that requires the number of future pulls in advance cannot be used unchanged when stopping time is uncertain. Perchet and Rigollet study successive elimination (SE), which removes poorly performing arms from a shrinking active set and can run without knowing the realized horizon. Their Theorem 2.1 bounds its expected regret when the horizon is random and independent of rewards Perchet and Rigollet, 2013, Theorem 2.1.

The static-bandit theorem also supplies a concentration and elimination analysis used in the paper's broader study of bandits with covariates. This mission isolates that theorem and the intermediate estimates that make its quantitative constant explicit. The paper proves the result in §2, with a concentration lemma in the Appendix; the goal here is a machine-checked statement and, eventually, a machine-checked proof of that known result.

Setting

There are K+1≥2K+1\ge2K+1≥2 arms. For each arm iii, the reward stack Yi,1,Yi,2,…Y_{i,1},Y_{i,2},\ldotsYi,1​,Yi,2​,… lists the rewards that successive pulls of that arm will reveal. Rewards lie in [0,1][0,1][0,1], are identically distributed along each arm, and have mean fif_ifi​. The stacks are mutually independent. Arms are indexed so that f1≤⋯≤fK+1f_1\le\cdots\le f_{K+1}f1​≤⋯≤fK+1​, and the last arm, denoted ∗*∗, is uniquely best. Its gap from arm iii is Δi=f∗−fi\Delta_i=f_*-f_iΔi​=f∗​−fi​. A gap is zero only for the best arm.

The policy works in rounds. At the start of round τ\tauτ, a nonempty set IτI_\tauIτ​ is active. It pulls each arm in IτI_\tauIτ​ once, in increasing index order, and computes the mean Yˉi(τ)\bar Y_i(\tau)Yˉi​(τ) from that arm's first τ\tauτ rewards. It retains an arm when its mean is within γU(τ,T)\gamma U(\tau,T)γU(τ,T) of the largest active mean, where T≥1T\ge1T≥1 is an integer, γ≥1\gamma\ge1γ≥1, and

log⁡‾(x)=max⁡{log⁡x,1},U(τ,T)=22log⁡‾(T/τ)τ.\overline{\log}(x)=\max\{\log x,1\},\qquad U(\tau,T)=2\sqrt{\frac{2\overline{\log}(T/\tau)}{\tau}}.log​(x)=max{logx,1},U(τ,T)=2τ2log​(T/τ)​​.

Concatenating the rounds gives the pull sequence π^1,π^2,…\widehat\pi_1,\widehat\pi_2,\ldotsπ1​,π2​,…. For a horizon NNN, its regret is RN(π^)=∑t=1NΔπ^tR_N(\widehat\pi)=\sum_{t=1}^{N}\Delta_{\widehat\pi_t}RN​(π)=∑t=1N​Δπt​​. If NNN is random, independent of the reward stacks, and has mean nnn, the paper studies ERN(π^)\mathbb E R_N(\widehat\pi)ERN​(π). For a threshold Δ>0\Delta>0Δ>0, let Δ−\Delta^-Δ− be the largest arm gap strictly smaller than Δ\DeltaΔ, or zero if none exists.

Formalization targets

The goal is Theorem 2.1, for every positive threshold Δ\DeltaΔ:

ERN(π^)≤392γ2(1+nT)KΔ log⁡‾ ⁣(TΔ218γ2)+nΔ−.\mathbb E R_N(\widehat\pi) \le 392\gamma^2\left(1+\frac nT\right) \frac K\Delta\, \overline{\log}\!\left(\frac{T\Delta^2}{18\gamma^2}\right) +n\Delta^-.ERN​(π)≤392γ2(1+Tn​)ΔK​log​(18γ2TΔ2​)+nΔ−.

The KKK in this formula counts suboptimal arms; there are K+1K+1K+1 arms in total. The paper writes Δ≥0\Delta\ge0Δ≥0 under a convention that division by zero is infinite. The real-valued Lean goal states the informative positive-threshold case.

The milestones include the Appendix's maximal Hoeffding–Azuma inequality and Lemma A.1, which controls all empirical averages up to a fixed time. In §2 they continue with the stopping-round estimate (2.2), the empirical-gap tail (2.5), the bound Φj(τ)≤4τ/T\Phi_j(\tau)\le4\tau/TΦj​(τ)≤4τ/T for a false elimination of the best arm, and the final comparison for ϕ(x)=log⁡‾(ax2)/x\phi(x)=\overline{\log}(a x^2)/xϕ(x)=log​(ax2)/x. Each has its own statement so that a proof of the goal can import the exact estimate it needs.

Significance

The bound separates the cost of learning arms with gaps above Δ\DeltaΔ from the remaining regret represented by nΔ−n\Delta^-nΔ−. It applies when the policy uses a tuning value TTT but the number of pulls actually made is the random variable NNN. Thus it gives a quantitative guarantee without asking the algorithm to know its stopping time in advance. Its constants and dependence on TTT, nnn, KKK, and γ\gammaγ are explicit Perchet and Rigollet, 2013, §2.

The result is proved in the paper, while the Lean declarations in this proposal are open statements. Completing them would add a reusable formal account of elimination on a reward stack and of the Appendix's time-uniform concentration estimate. Mathlib and existing platform theorems contain related Hoeffding–Azuma tools, but this paper's peeling inequality, its particular stopping rule, and the random-horizon regret conclusion still require proofs in the stated model.

Difficulty

An ordinary concentration bound at one sample size does not control every round in which an arm could eliminate the best arm. Taking a naive union bound over all rounds loses the dependence on TTT needed for the displayed regret constant. The Appendix instead states a maximal inequality and a bound that holds across a range of sample sizes. The regret calculation must then account for both kinds of error: removing the best arm and retaining a suboptimal arm for too long.

The stopping round depends on the gap through U(τ,T)U(\tau,T)U(τ,T), which itself contains log⁡‾(T/τ)\overline{\log}(T/\tau)log​(T/τ). Its upper bound, the event decomposition across ordered arms, and the final comparison of gap-dependent terms all retain the paper's exact constants. A proof that only establishes the asymptotic rate would not close this goal.

Formalization scope

Lean uses Fin (K+1) with zero-based indices, so its last element is the paper's arm K+1K+1K+1. The published stochastic-bandit definition supplies the mutually independent reward stack and arm means. The [0,1] bound, sorted means, and unique best arm are explicit assumptions. A separate martingale-difference definition supplies the Appendix's adapted, integrable sequence, including centering of its first term; the page's one-line definition leaves that first-step condition unstated. The concentration milestones require a nondegenerate interval a<ba<ba<b: the displayed maximal exponent divides by (b−a)2(b-a)^2(b−a)2, and Lemma A.1 is false at a=b=0a=b=0a=b=0 for a zero sequence and δ<1\delta<1δ<1.

The SE policy is a function of observed rewards in the stack. It uses the end-of-round test described in the text and proof of §2, with equality retained, since Policy 1's pseudocode tests at a different point in the round. The active set starts with all arms and remains nonempty; the total Lean pull function has a fallback value that is unreachable under the theorem's parameters. The policy never reads the unknown means. Empirical means and UUU are used only at positive rounds.

The horizon is measurable and integrable, has mean nnn, and is independent of the entire reward stack. Since regret is bounded by the horizon, its expectation is a genuine integral rather than the default value of a non-integrable Lean integral. The positive-threshold condition avoids Lean's zero-valued division at Δ=0\Delta=0Δ=0. Finite maxima are taken over finite arm sets, and log⁡‾\overline{\log}log​ is kept distinct from the ordinary logarithm in the Appendix's bounds. The paper's printed nnn in the logarithm of (2.7) and the final ϕ\phiϕ comparison is read as TTT, consistently with (2.2) and Theorem 2.1. Contributions to the concentration lemmas, stopping-round analysis, elimination events, and final regret proof all fit this scope.

Selected references

  • V. Perchet and P. Rigollet, The multi-armed bandit problem with covariates, Annals of Statistics 41(2), 2013, 693–721. arXiv:1110.6084v3; DOI:10.1214/13-AOS1101.
  • S. Bubeck and N. Cesa-Bianchi, Regret Analysis of Stochastic and Nonstochastic Multi-armed Bandit Problems, Foundations and Trends in Machine Learning 5(1), 2012. arXiv:1204.5721.
10 thms1 active userReviewed
PreviousPage 21 of 38Next

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