Prove2Me
Navigate
DiscoverCollectionsFormalpediaBlogsUsersMomentumMy Missions+
Prove2Me
⌕
Log in

Get started

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

Operations Research

2,633 missions · 1,357 completed

The discipline of applying mathematical analysis to complex decision problems in operations: allocating scarce resources, scheduling, routing, inventory, and the design of service and production systems. Drawing on mathematical programming, stochastic modeling, queueing, simulation, and game-theoretic reasoning, it seeks policies that perform provably well in systems shaped by constraints, congestion, and uncertainty.

Missions

Open1276Completed1357All2633
Algorithmic Game TheoryCombinatoricsProbability·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
Algorithmic Game TheoryCombinatoricsProbability·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
Graph TheoryLinear OptimizationOptimization·Captain: mikedeng1

Minimum Cost Flows in Graphs with Unit Capacities: When the Second Phase of a Cost Scale Starts, at Most 3λ Excess Is Left for λ = √m and 36λ for λ = n^(2/3)Research Paper

Motivation

Minimum-cost flow is a standard model for routing a quantity through a network while paying a cost on each arc. Unit capacities arise when each input arc can carry at most one unit; they include flow formulations of several matching and assignment problems. Goldberg, Kaplan, Hed and Tarjan analyze a cost-scaling algorithm for this setting. Their second-phase argument depends on a bound on the excess left in a pseudoflow when every excess vertex has undergone enough potential increases. This mission formalizes that bound, Lemma 5 of their STACS 2015 paper. The published paper gives an O(min⁡{m1/2,n2/3} mlog⁡(nC))O(\min\{m^{1/2},n^{2/3}\}\,m\log(nC))O(min{m1/2,n2/3}mlog(nC)) running-time bound for its unit-capacity algorithm; the running-time statement itself is outside this mission.

Setting

A network has a finite vertex set VVV and an arc set EEE that contains both orientations of each input arc. Write n=∣V∣n=|V|n=∣V∣. Every input arc has capacity one; its reverse has capacity zero and opposite cost. Write mmm for the number of input arcs, so ∣E∣=2m|E|=2m∣E∣=2m. The paper assumes no parallel or antiparallel input arcs, allowing an arc to be identified by its ordered endpoints. Costs are antisymmetric between orientations.

A pseudoflow fff assigns an antisymmetric quantity to each arc and respects its capacity. Its excess at vvv is the sum ef(v)e_f(v)ef​(v) of incoming flow. Positive excess means that vvv has flow to discharge; negative excess is a deficit. A circulation is a pseudoflow with zero excess at every vertex. An arc is residual for fff when its residual capacity uf(v,w)=u(v,w)−f(v,w)u_f(v,w)=u(v,w)-f(v,w)uf​(v,w)=u(v,w)−f(v,w) is positive. For vertex prices ppp, its reduced cost is cp(v,w)=c(v,w)−p(v)+p(w)c_p(v,w)=c(v,w)-p(v)+p(w)cp​(v,w)=c(v,w)−p(v)+p(w). A pseudoflow is ε\varepsilonε-optimal at ppp when every residual arc has reduced cost at least −ε-\varepsilon−ε. These are the paper's conventions on p. 408.

At the start of a cost scale, let f′f'f′ be a circulation that is 2ε2\varepsilon2ε-optimal at prices p′p'p′. During the scale, fff is a pseudoflow that is ε\varepsilonε-optimal at prices ppp. The potential increase in scale units is d(v)=(p(v)−p′(v))/εd(v)=(p(v)-p'(v))/\varepsilond(v)=(p(v)−p′(v))/ε. Put λ=min⁡{m,n2/3}\lambda=\min\{\sqrt m,n^{2/3}\}λ=min{m​,n2/3}. The second phase begins when every excess vertex has d(v)≥λd(v)\ge\lambdad(v)≥λ; the paper also states that every deficit then has d(v)=0d(v)=0d(v)=0 and that the levels d(v)d(v)d(v) are integral. The remaining excess is X(f)=∑v:ef(v)>0ef(v)X(f)=\sum_{v:e_f(v)>0}e_f(v)X(f)=∑v:ef​(v)>0​ef​(v).

The comparison graph has arcs E+={(v,w)∈E:f′(v,w)>f(v,w)}E^+=\{(v,w)\in E:f'(v,w)>f(v,w)\}E+={(v,w)∈E:f′(v,w)>f(v,w)} and vertices VVV. An excess–deficit cut Y⊆E+Y\subseteq E^+Y⊆E+ meets every path in this graph from an excess vertex to a deficit vertex. Thus a cut is measured among comparison arcs, not among arbitrary arcs of the original network. Lemmas 3 and 4 of the paper establish the level and cut bounds on which the excess estimate rests.

Formalization targets

The goal is Lemma 5 on p. 410, stated with both constants obtained in its proof. Under the cost-scale conditions above, and the divisibility conventions in footnotes 5–6, its two cases are

λ=m⟹X(f)≤3λ,\lambda=\sqrt m\quad\Longrightarrow\quad X(f)\le 3\lambda,λ=m​⟹X(f)≤3λ, λ=n2/3⟹X(f)≤36λ.\lambda=n^{2/3}\quad\Longrightarrow\quad X(f)\le36\lambda.λ=n2/3⟹X(f)≤36λ.

The second case also assumes that six divides the integer λ\lambdaλ; the first uses divisibility by three. The goal keeps the two conclusions separate so that the stronger 3λ3\lambda3λ bound is retained when m\sqrt mm​ determines the scale. Its milestone list follows the paper: Lemma 3 bounds d(v)−d(w)d(v)-d(w)d(v)−d(w) by three on an arc of E+E^+E+; Lemma 4 bounds total excess by a cut; and the proof of Lemma 5 supplies the cut bands and the two size estimates. The bands use the exact three-level and six-level ranges printed on p. 410.

Significance

The result connects a condition on price increases to a numerical bound on the amount of flow still needing discharge. Combined with the paper's time per unit of remaining excess, this is the bound needed for its second-phase running-time estimate. The two regimes show how the same cost-scale state can be controlled by either the number of input arcs or the number of vertices. Without the excess bound, the paper's argument gives no comparable limit on the remaining work at the phase transition.

The paper proves Lemma 5, but this particular unit-capacity excess argument has no matching theorem in the platform catalog checked for this mission. The network, circulation and pseudoflow interfaces are already published as reusable definitions from earlier cost-scaling work. A completed machine-checked development would add the cut argument and band counts to those interfaces. Those facts can also support other analyses that compare an entry circulation with an intermediate pseudoflow; this mission does not claim to formalize the algorithm's entire execution or its running-time model.

Difficulty

The main obstacle is linking local price inequalities to a global bound on excess. A reduced-cost inequality controls a single residual arc, while excess may move through a path and is measured over many vertices. The comparison graph G+G^+G+ is essential: on a residual arc outside E+E^+E+, the level difference used by Lemma 3 is not available. Counting all residual arcs as if they were comparison arcs would therefore not justify the cut bound. There is a separate counting issue: the represented network has both directions of each input arc, so using ∣E∣|E|∣E∣ for mmm changes the constant in the square-root case.

Formalization scope

Lean represents the network by a finite symmetric arc set with real capacities and antisymmetric real costs. A local unit-capacity predicate marks each reverse pair as capacities (1,0)(1,0)(1,0) or (0,1)(0,1)(0,1); mmm counts the capacity-one orientations. Flows are real valued as on p. 408. Integer costs on that page are generalized to real costs because Lemmas 3–5 use only antisymmetry and reduced-cost inequalities. No result here depends on cost integrality.

The formalization takes the cost-scale state as hypotheses: ε>0\varepsilon>0ε>0, entry circulation and 2ε2\varepsilon2ε-optimality, current ε\varepsilonε-optimality, integral levels, level at least λ\lambdaλ at excess vertices, and zero level at deficits. It takes λ\lambdaλ as a natural number equal to min⁡{m,n2/3}\min\{\sqrt m,n^{2/3}\}min{m​,n2/3} and records the paper's divisibility simplifications. The indexed pigeonhole milestones also require λ>0\lambda>0λ>0 so their index ranges are nonempty; the goal includes the zero-arc boundary without requiring an index. Cuts use reachability after deleting a subset of E+E^+E+, equivalent to meeting every excess-to-deficit path. The proof of Lemma 5 prints “cut in GfG_fGf​”; the formalized claim uses G+G^+G+, the graph for which Lemma 3 gives the required level control and Lemma 4 states its bound.

The definition layer, path-cut lemma and finite counting arguments are the reusable infrastructure. Contributions that establish these exact statements, their necessary elementary flow identities, or a machine-checked goal proof are in scope. Counting mmm with reverse arcs, allowing cuts outside E+E^+E+, using arbitrary real levels in place of integral ones, or replacing both cases by the weaker 36λ36\lambda36λ bound does not express this target.

Selected references

  • A. V. Goldberg, H. Kaplan, S. Hed and R. E. Tarjan, Minimum Cost Flows in Graphs with Unit Capacities, 32nd International Symposium on Theoretical Aspects of Computer Science, LIPIcs 30, pp. 406–419, 2015. doi:10.4230/LIPIcs.STACS.2015.406.
9 thms1 active userReviewed
Bandit AlgorithmsMachine LearningStatistics·Captain: mikedeng1

Personalized Dynamic Pricing with Machine Learning: High-Dimensional Features and Heterogeneous Elasticity 2: Under Unknown Sparsity, Lasso-Based ILQX Pricing Has Regret at Most Cs√T(log d + log T)Research Paper

Motivation

Online sellers of loans, insurance, travel and retail goods observe data about each arriving customer before quoting a price. Using that data to price individually is attractive only if the seller can learn, from its own sales, how demand responds to price for customers with given characteristics. Two features make this hard in practice. The number of recorded characteristics ddd is often large, possibly larger than the selling horizon TTT, while only a few of them matter. And price sensitivity is itself heterogeneous: different customers respond to the same price change differently.

Ban and Keskin (Management Science 67(9), 2021) model both features and ask how fast the revenue lost to learning can grow. Their answer is that a seller who learns with lasso-regularized quasi-likelihood estimation and occasional price experiments loses revenue of order sTs\sqrt TsT​ up to logarithmic factors, where sss is the number of relevant characteristics, and that no policy can do better than order sTs\sqrt TsT​. This mission formalizes the upper bound for the general (nonlinear-link) model with unknown sparsity, Theorem 3 of the paper.

Setting

In each period t=1,2,…,Tt = 1, 2, \dots, Tt=1,2,…,T a customer arrives with a raw feature vector Zt∈RdZ_t \in \mathbb R^dZt​∈Rd. The ZtZ_tZt​ are i.i.d., take values in a compact set Z\mathcal ZZ inside a ball of radius zmax⁡z_{\max}zmax​, have mean zero and positive definite covariance ΣZ\Sigma_ZΣZ​. The seller sees the augmented vector Xt=[1;Zt]∈Rd+1X_t = [1; Z_t] \in \mathbb R^{d+1}Xt​=[1;Zt​]∈Rd+1, charges a price pt∈[ℓ,u]p_t \in [\ell, u]pt​∈[ℓ,u] with 0<ℓ<u0 < \ell < u0<ℓ<u, and observes the demand

Dt=g(α⋅Xt+(β⋅Xt) pt)+εt.D_t = g\big(\alpha \cdot X_t + (\beta \cdot X_t)\, p_t\big) + \varepsilon_t .Dt​=g(α⋅Xt​+(β⋅Xt​)pt​)+εt​.

Here θ=(α,β)\theta = (\alpha, \beta)θ=(α,β) is an unknown parameter in a compact rectangle Θ⊂R2(d+1)\Theta \subset \mathbb R^{2(d+1)}Θ⊂R2(d+1), the link ggg is known, differentiable and increasing with ℓ~≤g′≤u~\tilde\ell \le g' \le \tilde uℓ~≤g′≤u~ on the relevant domain, and the demand shocks εt\varepsilon_tεt​ form a sub-Gaussian martingale difference sequence. The term β⋅Xt\beta \cdot X_tβ⋅Xt​ is the customer's own price sensitivity. With u(p,x)=[1;p]⊗xu(p, x) = [1; p] \otimes xu(p,x)=[1;p]⊗x the argument of ggg is θ⋅u(p,x)\theta \cdot u(p, x)θ⋅u(p,x).

The sparsity of θ\thetaθ is s=∣S∣s = |\mathcal S|s=∣S∣, S={i:αi≠0}∪{i:βi≠0}\mathcal S = \{i : \alpha_i \neq 0\} \cup \{i : \beta_i \neq 0\}S={i:αi​=0}∪{i:βi​=0}. The seller does not know S\mathcal SS.

The expected revenue of price ppp is r(p,θ,x)=p g(θ⋅u(p,x))r(p, \theta, x) = p\, g(\theta \cdot u(p,x))r(p,θ,x)=pg(θ⋅u(p,x)), and the clairvoyant price φ(θ,x)\varphi(\theta, x)φ(θ,x) maximizes it; it is assumed to lie in the interior of [ℓ,u][\ell, u][ℓ,u] for θ∈Θ\theta \in \Thetaθ∈Θ, x∈X={1}×Zx \in \mathcal X = \{1\} \times \mathcal Zx∈X={1}×Z. The regret of a policy over TTT periods is

Δθπ(T)=E[∑t=1Tr∗(θ,Xt)−r(pt,θ,Xt)],r∗(θ,x)=r(φ(θ,x),θ,x).\Delta^\pi_\theta(T) = \mathbb E\Big[\sum_{t=1}^T r^*(\theta, X_t) - r(p_t, \theta, X_t)\Big], \qquad r^*(\theta,x) = r(\varphi(\theta,x),\theta,x).Δθπ​(T)=E[t=1∑T​r∗(θ,Xt​)−r(pt​,θ,Xt​)],r∗(θ,x)=r(φ(θ,x),θ,x).

The policy ILQX(m1,m2,λ)(m_1, m_2, \lambda)(m1​,m2​,λ) charges the experimental price m1m_1m1​ in periods M1={1,4,9,… }M_1 = \{1, 4, 9, \dots\}M1​={1,4,9,…} and m2m_2m2​ in periods M2={2,5,10,… }M_2 = \{2, 5, 10, \dots\}M2​={2,5,10,…}, and otherwise charges φ(PΘθ^t,Xt)\varphi(\mathcal P_\Theta \hat\theta_t, X_t)φ(PΘ​θ^t​,Xt​), where θ^t\hat\theta_tθ^t​ maximizes the lasso-penalized quasi-likelihood

Qt−1(θ~,λt)=∑k=1t−1χk∫Dkg(θ~⋅uk)Dk−yg′(g−1(y)) dy−λt∥θ~∥1Q_{t-1}(\tilde\theta, \lambda_t) = \sum_{k=1}^{t-1} \chi_k \int_{D_k}^{g(\tilde\theta \cdot u_k)} \frac{D_k - y}{g'(g^{-1}(y))}\, dy - \lambda_t \|\tilde\theta\|_1Qt−1​(θ~,λt​)=k=1∑t−1​χk​∫Dk​g(θ~⋅uk​)​g′(g−1(y))Dk​−y​dy−λt​∥θ~∥1​

over the experimental periods (χk=I{k∈M1∪M2}\chi_k = \mathbb I\{k \in M_1 \cup M_2\}χk​=I{k∈M1​∪M2​}), and PΘ\mathcal P_\ThetaPΘ​ is the projection onto Θ\ThetaΘ. The regularization level is λt+1=c~ t1/4log⁡d+log⁡t\lambda_{t+1} = \tilde c\, t^{1/4} \sqrt{\log d + \log t}λt+1​=c~t1/4logd+logt​ for a constant c~>0\tilde c > 0c~>0.

Formalization targets

Goal: Theorem 3

There is a finite constant C~>0\tilde C > 0C~>0 such that

Δθπ(T)≤C~ sT (log⁡d+log⁡T)for all θ∈Θ, T≥2.\Delta^\pi_\theta(T) \le \tilde C\, s \sqrt T\, (\log d + \log T) \qquad \text{for all } \theta \in \Theta,\ T \ge 2 .Δθπ​(T)≤C~sT​(logd+logT)for all θ∈Θ, T≥2.

C~\tilde CC~ is fixed with the model and is uniform over Θ\ThetaΘ; since s=s(θ)s = s(\theta)s=s(θ) varies over Θ\ThetaΘ, the factor sss carries content.

Milestones

  1. The counting property of the schedule (7): for t≥5t \ge 5t≥5 each experimental price is charged at least 14t\tfrac14 \sqrt t41​t​ times in periods 1,…,t1, \dots, t1,…,t.
  2. Lemma 3, the estimation error of the lasso quasi-likelihood estimate:
P{∥θ^t+1(lasso)(λt+1)−θ∥2≤ρ3s(log⁡d+log⁡t)t}≥1−κ3s(log⁡d+log⁡t)t,t≥t1.\mathbb P\Big\{\|\hat\theta^{(\mathrm{lasso})}_{t+1}(\lambda_{t+1}) - \theta\|^2 \le \rho_3 \frac{s(\log d + \log t)}{\sqrt t}\Big\} \ge 1 - \kappa_3 \frac{s(\log d + \log t)}{\sqrt t}, \qquad t \ge t_1 .P{∥θ^t+1(lasso)​(λt+1​)−θ∥2≤ρ3​t​s(logd+logt)​}≥1−κ3​t​s(logd+logt)​,t≥t1​.

Significance

Theorem 3 says that a seller facing many potentially relevant customer characteristics, without knowing which matter and under a nonlinear demand link such as logit, loses revenue at a rate governed by the number sss of relevant characteristics and only logarithmically by the total number ddd. Combined with the paper's sTs\sqrt TsT​ lower bound (Theorem 1, for the linear model), it makes ILQX first-order optimal up to logarithmic factors. An unregularized policy would pay order dTd\sqrt TdT​, which is linear in TTT when ddd exceeds T\sqrt TT​.

The paper's proofs are in its electronic companion, which this mission does not reproduce; no machine-checked version of the result exists. A complete formalization would also produce a reusable high-probability error bound for lasso-regularized quasi-likelihood estimation under adaptively collected data (Lemma 3), and a machine-checked model of feature-based dynamic pricing on which related policies can be stated.

Difficulty

The estimation error must be controlled on data collected by the policy itself. The design vectors uk=[1;mi(k)]⊗Xku_k = [1; m_{i(k)}] \otimes X_kuk​=[1;mi(k)​]⊗Xk​ use only two prices, so the information about the price sensitivity accrues only from the Θ(t)\Theta(\sqrt t)Θ(t​) experimental periods. Standard lasso error bounds assume a fixed or i.i.d. design and a restricted-eigenvalue condition on it; here a restricted-eigenvalue-type property of the random Kronecker design must be derived from ΣZ≻0\Sigma_Z \succ 0ΣZ​≻0 and m1≠m2m_1 \ne m_2m1​=m2​, and the noise enters through a martingale whose conditional law may depend on past prices. The quasi-likelihood is not a least-squares objective, so the curvature comes from the bounds ℓ~≤g′≤u~\tilde\ell \le g' \le \tilde uℓ~≤g′≤u~, which hold only on the bounded domain {θ⋅u(p,x):θ∈Θ}\{\theta \cdot u(p,x) : \theta \in \Theta\}{θ⋅u(p,x):θ∈Θ} and not at an unconstrained estimate. Finally, the regret bound has to turn a squared-error bound holding with probability 1−O(slog⁡(dt)/t)1 - O(s \log(dt)/\sqrt t)1−O(slog(dt)/t​) into an expected revenue loss, which requires the revenue loss to be second order in the parameter error at an interior optimum.

Formalization scope

  • Coordinates. x∈Rd+1x \in \mathbb R^{d+1}x∈Rd+1 is Fin (d+1) → ℝ with the intercept at index 0; θ\thetaθ is Fin 2 × Fin (d+1) → ℝ with θ(0,i)=αi\theta(0,i) = \alpha_iθ(0,i)=αi​, θ(1,i)=βi\theta(1,i) = \beta_iθ(1,i)=βi​. Periods are 1-based. d≥1d \ge 1d≥1.
  • Objective. The integrated form ∑kχk(Dk (θ~⋅uk)−∫0θ~⋅ukg)−λ∥θ~∥1\sum_k \chi_k \big(D_k\,(\tilde\theta \cdot u_k) - \int_0^{\tilde\theta \cdot u_k} g\big) - \lambda \|\tilde\theta\|_1∑k​χk​(Dk​(θ~⋅uk​)−∫0θ~⋅uk​​g)−λ∥θ~∥1​, which differs from the paper's (17) by a term not depending on θ~\tilde\thetaθ~ and therefore has the same maximizers.
  • Estimator. Any measurable selection that is a maximizer whenever one exists; the theorems hold for every such selection. The paper's claim that (17) is strictly concave with a unique maximizer fails when 2(d+1)2(d+1)2(d+1) exceeds the rank of the design.
  • Policy. The price is defined by strong recursion over the realized history; the estimator sees θ\thetaθ only through the realized demands.
  • Regret is a lower Lebesgue integral of nonnegative losses; conditional moments of the shocks use the conditional Lebesgue expectation; probabilities are compared in [0,∞][0, \infty][0,∞].
  • Standing assumptions formalized: i.i.d. measurable features in a compact Z\mathcal ZZ inside a Euclidean ball, mean zero, positive definite ΣZ\Sigma_ZΣZ​; fresh customers (Zt+1Z_{t+1}Zt+1​ independent of the past); sub-Gaussian martingale-difference shocks; m1≠m2m_1 \ne m_2m1​=m2​ in [ℓ,u][\ell, u][ℓ,u]; ggg differentiable and increasing with derivative bounds on the relevant domain (endnote 2 states them as a consequence, which is false; here they are a hypothesis); φ\varphiφ an unconstrained revenue maximizer in (ℓ,u)(\ell, u)(ℓ,u) whose pricing map after projection and feature augmentation is measurable. The processes are constrained only at the used periods t≥1t\ge1t≥1.
  • Added hypotheses, per θ\thetaθ: s(θ)≥1s(\theta) \ge 1s(θ)≥1 (the paper's s∈{1,…,d+1}s \in \{1, \dots, d+1\}s∈{1,…,d+1}), and almost surely every realized demand lies in the closure of the range of ggg, which guarantees that a maximizer of (17) exists (true for the linear link and for binary demand under a logit link).
  • Not formalized: the paper's support condition on continuous features; Remark 7's independence of C~\tilde CC~ from ddd; dependence of the shock law on θ\thetaθ (the shocks are one fixed process).
  • Ruled-out trivializations. The regret cannot vanish by non-integrability; the estimator cannot be an arbitrary function, since it must maximize the objective; C~\tilde CC~ and κ3,ρ3,t1\kappa_3, \rho_3, t_1κ3​,ρ3​,t1​ are chosen before θ\thetaθ, so they cannot absorb sss; the hypotheses are jointly satisfiable (a linear-link instance with d=1d = 1d=1 is checked).

Contributions welcome: existence of a measurable maximizer of the objective, the counting lemma, and supporting lemmas for Lemma 3 and the regret bound.

Selected references

  • G.-Y. Ban and N. B. Keskin, Personalized Dynamic Pricing with Machine Learning: High-Dimensional Features and Heterogeneous Elasticity, Management Science 67(9):5549–5568, 2021. https://doi.org/10.1287/mnsc.2020.3680
  • A. V. den Boer and B. Zwart, Simultaneously Learning and Optimizing Using Controlled Variance Pricing, Management Science 60(3):770–783, 2014. https://doi.org/10.1287/mnsc.2013.1788
  • N. B. Keskin and A. Zeevi, Dynamic Pricing with an Unknown Demand Model: Asymptotically Optimal Semi-Myopic Policies, Operations Research 62(5):1142–1167, 2014. https://doi.org/10.1287/opre.2014.1294
  • A. Javanmard and H. Nazerzadeh, Dynamic Pricing in High-Dimensions, Journal of Machine Learning Research 20(9):1–49, 2019. https://jmlr.org/papers/v20/17-357.html
6 thms1 active userReviewed
Dynamic ProgrammingMarkov ChainProbability·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
Algorithmic Game TheoryFunctional Analysis·Captain: mikedeng1

Stochastic Graphon Games: I. The Static Case 1: If the State Map b Is Bounded and √c_z‖W‖ < 1, the Graphon Game Has at Least One Nash EquilibriumResearch Paper

Motivation

Games with many interacting players appear in economics, epidemiology and network engineering: firms competing in a market, individuals deciding whether to vaccinate, agents in a social network choosing an effort level. When every player interacts with every other player symmetrically, mean field game theory replaces the finite game by a limit with a continuum of players. Real interactions are rarely symmetric: who influences whom is described by a network, and large networks are described in the limit by a graphon (Lovász, Large networks and graph limits, 2012).

Graphon games combine the two ideas. Parise and Ozdaglar (2018) studied deterministic graphon games with continuous strategy sets. Carmona, Cooney, Graves and Laurière (2019) introduced idiosyncratic random shocks, giving a stochastic graphon game in which each player's state is perturbed by independent noise, and proved existence, uniqueness and stability of Nash equilibria as well as convergence of finite network games to the graphon limit. This mission formalizes their existence theorem (Theorem 3.15).

Setting

The players form the interval I=[0,1]I=[0,1]I=[0,1] with Lebesgue measure λI\lambda_IλI​. Each player x∈Ix\in Ix∈I chooses an action αx∈A=R\alpha_x\in A=\mathbb Rαx​∈A=R; a strategy profile is a function α∈L2(I)\alpha\in L^2(I)α∈L2(I), defined up to λI\lambda_IλI​-null sets.

A graphon is a symmetric, Borel-measurable function w:I×I→Rw:I\times I\to\mathbb Rw:I×I→R with ∫I×Iw(x,y)2 dx dy<∞\int_{I\times I}w(x,y)^2\,dx\,dy<\infty∫I×I​w(x,y)2dxdy<∞. It acts on L2(I)L^2(I)L2(I) through the graphon operator

[Wg]x=∫Iw(x,y) g(y) dy,∥W∥=sup⁡∥φ∥L2(I)=1∥Wφ∥L2(I).[\mathbf Wg]_x=\int_I w(x,y)\,g(y)\,dy,\qquad \|\mathbf W\|=\sup_{\|\varphi\|_{L^2(I)}=1}\|\mathbf W\varphi\|_{L^2(I)}.[Wg]x​=∫I​w(x,y)g(y)dy,∥W∥=∥φ∥L2(I)​=1sup​∥Wφ∥L2(I)​.

A state map b:R×R→Rb:\mathbb R\times\mathbb R\to\mathbb Rb:R×R→R and a noise law μ0\mu_0μ0​ on R\mathbb RR define the state Xα,z,ξ=b(α,z)+ξX_{\alpha,z,\xi}=b(\alpha,z)+\xiXα,z,ξ​=b(α,z)+ξ of a player who takes action α\alphaα, feels aggregate zzz, and receives shock ξ∼μ0\xi\sim\mu_0ξ∼μ0​. Assumption 1 asks that ∣b(α,z)−b(α′,z′)∣2≤cα∣α−α′∣2+cz∣z−z′∣2|b(\alpha,z)-b(\alpha',z')|^2\le c_\alpha|\alpha-\alpha'|^2+c_z|z-z'|^2∣b(α,z)−b(α′,z′)∣2≤cα​∣α−α′∣2+cz​∣z−z′∣2 and that μ0\mu_0μ0​ has mean zero and a finite second moment.

The aggregate Zα\mathbf Z\alphaZα felt by the players under a profile α\alphaα is the solution z∈L2(I)z\in L^2(I)z∈L2(I) of

zx=∫Iw(x,y) b(αy,zy) dyfor λI-a.e. x.(6)z_x=\int_I w(x,y)\,b(\alpha_y,z_y)\,dy\qquad\text{for }\lambda_I\text{-a.e. }x.\tag{6}zx​=∫I​w(x,y)b(αy​,zy​)dyfor λI​-a.e. x.(6)

A player choosing α\alphaα against aggregate zzz pays J(α,z)=∫f(b(α,z)+ξ,α,z) μ0(dξ)J(\alpha,z)=\int f(b(\alpha,z)+\xi,\alpha,z)\,\mu_0(d\xi)J(α,z)=∫f(b(α,z)+ξ,α,z)μ0​(dξ). Assumption 4 asks that J(⋅,z)J(\cdot,z)J(⋅,z) be continuously differentiable and ℓc\ell_cℓc​-strongly convex uniformly in zzz, and that ∂αJ\partial_\alpha J∂α​J be ℓJ\ell_JℓJ​-Lipschitz in zzz. Assumption 5 asks that bbb be bounded: ∣b(α,z)∣≤c0|b(\alpha,z)|\le c_0∣b(α,z)∣≤c0​ for some finite c0>0c_0>0c0​>0. The best response [Bz]x[\mathbf Bz]_x[Bz]x​ is the unique minimizer of α↦J(α,zx)\alpha\mapsto J(\alpha,z_x)α↦J(α,zx​).

A Nash equilibrium is a profile α^∈L2(I)\hat\alpha\in L^2(I)α^∈L2(I) such that, for λI\lambda_IλI​-a.e. player xxx and every action β\betaβ,

J(α^x,(Zα^)x)≤J(β,(Zα^)x).J(\hat\alpha_x,(\mathbf Z\hat\alpha)_x)\le J(\beta,(\mathbf Z\hat\alpha)_x).J(α^x​,(Zα^)x​)≤J(β,(Zα^)x​).

Formalization targets

Goal: Theorem 3.15 (existence)

Assumptions 1, 4, 5 and cz ∥W∥<1 ⟹ there is at least one Nash equilibrium.\text{Assumptions 1, 4, 5 and }\sqrt{c_z}\,\|\mathbf W\|<1\ \Longrightarrow\ \text{there is at least one Nash equilibrium.}Assumptions 1, 4, 5 and cz​​∥W∥<1 ⟹ there is at least one Nash equilibrium.

The statement fixes no constants and asserts no uniqueness.

Milestones

  1. Proposition 3.1: under cz∥W∥<1\sqrt{c_z}\|\mathbf W\|<1cz​​∥W∥<1, (6) has a unique solution in L2(I)L^2(I)L2(I) for every α∈L2(I)\alpha\in L^2(I)α∈L2(I).
  2. Lemma 3.7: ∥Bz1−Bz2∥L2≤(ℓJ/ℓc)∥z1−z2∥L2\|\mathbf Bz^1-\mathbf Bz^2\|_{L^2}\le(\ell_J/\ell_c)\|z^1-z^2\|_{L^2}∥Bz1−Bz2∥L2​≤(ℓJ​/ℓc​)∥z1−z2∥L2​.
  3. Lemma 3.16, (15): ∥Zα1−Zα2∥L2≤cα∥W∥1−cz∥W∥∥α1−α2∥L2\|\mathbf Z\alpha^1-\mathbf Z\alpha^2\|_{L^2}\le\dfrac{\sqrt{c_\alpha}\|\mathbf W\|}{1-\sqrt{c_z}\|\mathbf W\|}\|\alpha^1-\alpha^2\|_{L^2}∥Zα1−Zα2∥L2​≤1−cz​​∥W∥cα​​∥W∥​∥α1−α2∥L2​.
  4. Proof of Theorem 3.15, first claim: with χ=c0\chi=c_0χ=c0​ and Bχ={g∈L2(I):∥g∥≤χ}B_\chi=\{g\in L^2(I):\|g\|\le\chi\}Bχ​={g∈L2(I):∥g∥≤χ}, every aggregate lies in W(Bχ)\mathbf W(B_\chi)W(Bχ​).
  5. Proof of Theorem 3.15, second claim: W(Bχ)\mathbf W(B_\chi)W(Bχ​) is relatively compact in L2(I)L^2(I)L2(I).

Significance

The existence theorem holds under mild hypotheses: unlike the uniqueness result of the same paper (Proposition 3.17), it needs no smallness condition linking ℓJ/ℓc\ell_J/\ell_cℓJ​/ℓc​, cαc_\alphacα​ and ∥W∥\|\mathbf W\|∥W∥. Boundedness of the state map takes its place. Existence of a graphon equilibrium is the starting point for the paper's §4, where graphon equilibria are shown to be limits of equilibria of finite network games and to yield ϵ\epsilonϵ-Nash equilibria for them.

The result is proved in the paper; no part of it is machine-checked. A complete formalization yields reusable pieces beyond this mission: square-integrable kernel operators on L2[0,1]L^2[0,1]L2[0,1] with their operator norm, the Hilbert–Schmidt compactness of such operators, Lipschitz dependence of minimizers of uniformly strongly convex functions on a parameter, and Schauder's fixed-point theorem on a Banach space, which is not in Mathlib.

Difficulty

A contraction argument does not apply: the best-response-of-aggregate map BZ\mathbf B\mathbf ZBZ is Lipschitz with constant ℓJℓc⋅cα∥W∥1−cz∥W∥\frac{\ell_J}{\ell_c}\cdot\frac{\sqrt{c_\alpha}\|\mathbf W\|}{1-\sqrt{c_z}\|\mathbf W\|}ℓc​ℓJ​​⋅1−cz​​∥W∥cα​​∥W∥​, which may exceed 111 under the hypotheses of the theorem. Existence therefore rests on a compactness fixed-point theorem in the infinite-dimensional space L2(I)L^2(I)L2(I). That requires compactness of the integral operator W\mathbf WW, which Mathlib does not provide for kernel operators, and Schauder's theorem itself, which Mathlib does not have. A bounded convex set in L2(I)L^2(I)L2(I) is not compact, so Brouwer-type finite-dimensional arguments do not transfer directly.

Formalization scope

All statements live in the namespace GraphonGames.Existence and share one definitions file.

  • III is Mathlib's unitInterval with volume. Profiles and aggregates are functions I→RI\to\mathbb RI→R with MemLp _ 2 volume; equalities between them are almost-everywhere equalities.
  • The action set is R\mathbb RR.
  • ∥W∥\|\mathbf W\|∥W∥ is the supremum, in [0,∞][0,\infty][0,∞], of ∥Wφ∥L2\|\mathbf W\varphi\|_{L^2}∥Wφ∥L2​ over ∥φ∥L2=1\|\varphi\|_{L^2}=1∥φ∥L2​=1; the graphon condition makes it finite, and the hypothesis cz∥W∥<1\sqrt{c_z}\|\mathbf W\|<1cz​​∥W∥<1 uses its real value.
  • The aggregate is a predicate: zzz is an aggregate of α\alphaα when z∈L2(I)z\in L^2(I)z∈L2(I) and (6) holds a.e. Remark 3.3's identity Zα=W[b(α,Zα)]\mathbf Z\alpha=\mathbf W[b(\alpha,\mathbf Z\alpha)]Zα=W[b(α,Zα)] is then the definition itself, the noise having mean zero.
  • The best response is a choice of minimizer, which exists and is unique under Assumption 4.
  • Pinned definition. The paper defines equilibrium (Definition 3.4) as a fixed point of a best response built on a rich Fubini extension, and states in Proposition 3.6 that under Assumption 4 this is equivalent to the a.e. best-response inequality above. The formalization takes that inequality as the definition. Fubini extensions and the exact law of large numbers (Proposition 3.2) are not formalized.
  • Norm inequalities are stated in [0,∞][0,\infty][0,∞], so an infinite left-hand side cannot make a bound vacuous.

Trivializing encodings are excluded: the operator norm is not a real supremum that could collapse to 000; an equilibrium must come with its aggregate solving (6), be square-integrable, and beat every action β∈R\beta\in\mathbb Rβ∈R for almost every player; and no norm is converted to a real number where it could be infinite.

The results of §4 (Theorems 4.3, 4.7, 4.9 on finite network games) and the uniqueness and stability results (Proposition 3.17, Theorem 3.20) are not posed in this mission. Contributions welcome include Schauder's fixed-point theorem, the compactness of Hilbert–Schmidt integral operators on L2L^2L2, and the contraction argument behind Proposition 3.1.

Selected references

  • R. Carmona, D. B. Cooney, C. V. Graves, M. Laurière, Stochastic Graphon Games: I. The Static Case, arXiv:1911.10664v1, 2019; Mathematics of Operations Research 47(1), 2022. https://arxiv.org/abs/1911.10664
  • F. Parise, A. Ozdaglar, Graphon games, arXiv:1802.00080, 2018. https://arxiv.org/abs/1802.00080
  • L. Lovász, Large Networks and Graph Limits, AMS Colloquium Publications 60, 2012. https://doi.org/10.1090/coll/060
  • R. Carmona, F. Delarue, Probabilistic Theory of Mean Field Games with Applications I, Springer, 2018. https://doi.org/10.1007/978-3-319-58920-6
7 thms1 active userReviewed
Bandit AlgorithmsProbabilityStatistics·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
Probability·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
Bandit AlgorithmsMachine LearningStatistics·Captain: mikedeng1

Personalized Dynamic Pricing with Machine Learning: High-Dimensional Features and Heterogeneous Elasticity 1: Under Known Sparsity, ILSX Pricing Has Regret at Most Cs√T log TResearch Paper

Motivation

A seller who changes prices while learning demand faces two costs: a price experiment may earn less revenue now, while insufficient experimentation leaves later prices poorly calibrated. Customer features make the balance harder because both willingness to buy and sensitivity to price may vary between customers. Ban and Keskin study this problem when the number of recorded features can be large relative to the selling horizon. Their known-support policy, ILSX, provides the upper bound formalized here Ban and Keskin, 2021, §§2 and 4.1.

The paper's result is a statement about learning a demand model while serving sequential customers, rather than about fitting that model after all observations are available. A price at period ttt must use information then available, including the arriving customer's features but excluding that customer's demand shock. The seller therefore needs a policy whose expected revenue loss can be bounded uniformly over the unknown parameters in a specified model Ban and Keskin, 2021, pp. 5552–5556.

Setting

In each period t=1,2,…t=1,2,\ldotst=1,2,…, a customer arrives with a feature vector Zt∈RdZ_t\in\mathbb R^dZt​∈Rd. The seller observes the augmented vector Xt=[1;Zt]X_t=[1;Z_t]Xt​=[1;Zt​], posts a price pt∈[ℓ,u]p_t\in[\ell,u]pt​∈[ℓ,u] with 0<ℓ<u0<\ell<u0<ℓ<u, and then observes demand. In the linear demand model, demand equals

Dt=α⋅Xt+(β⋅Xt)pt+εt.D_t=\alpha\cdot X_t+(\beta\cdot X_t)p_t+\varepsilon_t.Dt​=α⋅Xt​+(β⋅Xt​)pt​+εt​.

The intercept and slope vectors α,β\alpha,\betaα,β are unknown. The seller knows the sparsity structure S\mathcal SS: outside this set of coordinate indices, both parameter vectors vanish. The compressed parameter θS=(αS,βS)\theta_{\mathcal S}=(\alpha_{\mathcal S},\beta_{\mathcal S})θS​=(αS​,βS​) lies in a compact coordinate rectangle ΘS\Theta_{\mathcal S}ΘS​, and s=∣S∣≥1s=|\mathcal S|\ge1s=∣S∣≥1. The compressed customer vector is XS,t=(Xt,i)i∈SX_{\mathcal S,t}=(X_{t,i})_{i\in\mathcal S}XS,t​=(Xt,i​)i∈S​ Ban and Keskin, 2021, pp. 5552–5555.

At price ppp, expected single-period revenue is r(p,θS,x)=p[αS⋅x+(βS⋅x)p]r(p,\theta_{\mathcal S},x)=p[\alpha_{\mathcal S}\cdot x+(\beta_{\mathcal S}\cdot x)p]r(p,θS​,x)=p[αS​⋅x+(βS​⋅x)p]. The clairvoyant price is φ(θS,x)=−(αS⋅x)/(2βS⋅x)\varphi(\theta_{\mathcal S},x)=-(\alpha_{\mathcal S}\cdot x)/(2\beta_{\mathcal S}\cdot x)φ(θS​,x)=−(αS​⋅x)/(2βS​⋅x) when the slope is negative. The paper assumes this price is in (ℓ,u)(\ell,u)(ℓ,u) throughout its model domain. Expected regret ΔθSπ(T)\Delta^{\pi}_{\theta_{\mathcal S}}(T)ΔθS​π​(T) is the expected sum, through period TTT, of the clairvoyant revenue minus the revenue at the policy's price Ban and Keskin, 2021, (2)–(4), p. 5553–5554.

The ILSX policy charges one fixed experimental price m1m_1m1​ at periods 1,4,9,…1,4,9,\ldots1,4,9,… and a distinct price m2m_2m2​ at periods 2,5,10,…2,5,10,\ldots2,5,10,…. At every other period it estimates the compressed demand parameter by least squares using experimental observations, projects the estimate onto ΘS\Theta_{\mathcal S}ΘS​, and applies φ\varphiφ to the projected estimate and the current customer vector. Both prices lie in [ℓ,u][\ell,u][ℓ,u] Ban and Keskin, 2021, (7), (8), (12), pp. 5555–5556.

Formalization targets

The goal is Theorem 2: for the stated model and ILSX policy, a finite positive constant CCC exists such that

ΔθSILSX(m1,m2)(T)≤CsTlog⁡T(θS∈ΘS, T≥2).\Delta^{\mathrm{ILSX}(m_1,m_2)}_{\theta_{\mathcal S}}(T)\le Cs\sqrt T\log T\qquad(\theta_{\mathcal S}\in\Theta_{\mathcal S},\ T\ge2).ΔθS​ILSX(m1​,m2​)​(T)≤CsT​logT(θS​∈ΘS​, T≥2).

The constant is chosen after the model, shock process, experimental prices, and estimator selection are fixed, but before the parameter and horizon. The theorem and Remark 5 are on pp. 5556–5557 of Ban and Keskin, 2021.

The milestones reproduce the paper's route into this claim: the experimental-period count following (7); the growth bound (15) for the price-variation statistic JtJ_tJt​; Lemma 1's minimum-eigenvalue probability bound for the experimental information matrix J~S,t\widetilde{\mathcal J}_{\mathcal S,t}J​S,t​; and Lemma 2's high-probability bound for its noise-driven estimation error. The least-squares error identity (11) is included as a supporting item on invertible samples, not as a milestone. The constants and thresholds in these targets remain existential where the article leaves them existential Ban and Keskin, 2021, pp. 5555–5556.

Significance

The theorem gives a guarantee for a concrete pricing policy under feature-dependent price sensitivity. The bound grows with the horizon at order Tlog⁡T\sqrt T\log TT​logT for a fixed known support, with an explicit factor sss in the published statement. It lets a seller compare the cost of this policy with the revenue of a seller who knew the demand coefficients before pricing. The paper presents its result as matching its lower-bound benchmark up to logarithmic factors, under the paper's respective hypotheses Ban and Keskin, 2021, Theorems 1–2, pp. 5554–5556.

The mathematical result was published in 2021. This mission seeks a machine-checked statement and proof of the known-support guarantee and its stated intermediate results. The checked Prove2Me prior-art catalog contains featureless dynamic-pricing results and separate least-squares results, but no matching formalization of ILSX. A completed development would add reusable definitions for a feature-dependent price regressor, experimental information matrix, nonanticipating estimator selection, and expected regret. The article places its proofs in an electronic companion, which is not included in the held source PDF; this mission does not reproduce those proofs as source material.

Difficulty

The policy learns from a sparse set of deliberately chosen periods, while its regret accumulates at every period. A pointwise least-squares identity alone does not control the regret: the information matrix can be singular at early periods, and the customer features and shocks are random. In particular, the printed inverse identity (11) says “for all t≥2t\ge2t≥2”, although at t=2t=2t=2 the matrix has rank at most two and cannot be invertible when s≥2s\ge2s≥2. The formal statement therefore confines the inverse identity to invertible samples. Lemma 2's event also includes invertibility so a default inverse of a singular matrix cannot make its error bound automatic Ban and Keskin, 2021, (9)–(11), Lemma 2, p. 5556.

Formalization scope

Lean uses one-based periods and a support subtype of the augmented feature coordinates. Coordinate zero is the paper's intercept coordinate one. The parameter has coordinates (j,i)(j,i)(j,i) for demand intercept or price slope j∈{0,1}j\in\{0,1\}j∈{0,1} and support index i∈Si\in\mathcal Si∈S. The regressor [1;p]⊗XS,t[1;p]\otimes X_{\mathcal S,t}[1;p]⊗XS,t​ gives the information matrix in (10). Euclidean squared norms are explicit sums of squares, and the minimum-eigenvalue claim is the equivalent Rayleigh inequality. The parameter rectangle is nonempty, and its Euclidean projection is coordinatewise clamping.

The standing assumptions are d≥1d\ge1d≥1, s≥1s\ge1s≥1, 0<ℓ<u0<\ell<u0<ℓ<u, distinct experimental prices in [ℓ,u][\ell,u][ℓ,u], bounded i.i.d. features with mean zero and positive-definite covariance, and a fixed shock process with conditional mean zero, bounded conditional second moment, and finite conditional exponential moments near zero. The feature process is measurable and integrable. The model additionally states independence of the next customer's features from the past primitives, the reading of the paper's joint probability construction used here. The slope is negative and the clairvoyant price lies in (ℓ,u)(\ell,u)(ℓ,u) for every allowed parameter and realized feature. These conditions ensure that the price formula represents a genuine revenue maximum. The phrase about continuous features having positive measure in the interior of their domains is omitted because the page does not specify the domains or a reference measure Ban and Keskin, 2021, §2, pp. 5552–5554.

The estimator is any measurable minimizer of the experimental least-squares objective for histories with at least two observations. Its type depends only on earlier observations; at nonexperimental periods the recorded pseudo-observations are ignored by the objective. The shock law is fixed while the parameter varies. Expected regret uses a nonnegative integral, and conditional second and exponential moments use nonnegative conditional expectations. The theorem binds CCC before the parameter and horizon, while leaving support size fixed with the model. Consequently, the stronger uniformity in sss suggested by Remark 5 is outside this statement. These choices exclude vacuous integrals, unconstrained estimators, and accidental bounds from a singular matrix inverse. Work on the probability estimates, least-squares identity, and regret bound is welcome under these exact definitions.

Selected references

  • G.-Y. Ban and N. B. Keskin, Personalized Dynamic Pricing with Machine Learning: High-Dimensional Features and Heterogeneous Elasticity, Management Science 67(9):5549–5568, 2021. DOI: 10.1287/mnsc.2020.3680.
6 thms1 active userReviewed
Algorithmic Game TheoryProbability·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
OptimizationProbabilityStatistics·Captain: mikedeng1

Robust Sample Average Approximation 1: A Data-Driven DRO Converges in Objective, Optimal Value and Optimal Solutions for Every Admissible Cost If and Only If Its Test Is Uniformly ConsistentResearch Paper

Motivation

A decision maker who must choose xxx to minimize an expected cost EF[c(x;ξ)]\mathbb E_F[c(x;\xi)]EF​[c(x;ξ)] rarely knows the distribution FFF of the uncertain parameter ξ\xiξ; what is available is a sample ξ1,…,ξN\xi^1,\dots,\xi^Nξ1,…,ξN. Sample average approximation (SAA) replaces FFF by the empirical distribution. Data-driven distributionally robust optimization (DRO) instead minimizes the worst case of the expected cost over a set FN\mathcal F_NFN​ of distributions built from the data, which protects against the error of the empirical distribution at small NNN. Many such sets have been proposed: moment sets (Delage and Ye, 2010), ϕ\phiϕ-divergence balls (Ben-Tal et al., 2013), Wasserstein balls (Mohajerin Esfahani and Kuhn, 2018).

Bertsimas, Gupta and Kallus (arXiv:1408.4445v3; Mathematical Programming, 2018, doi:10.1007/s10107-017-1174-z) observe that each such set can be read as the confidence region of a statistical goodness-of-fit (GoF) test, and call the resulting method Robust SAA. A basic question for any such method is whether, as N→∞N\to\inftyN→∞, its objective, optimal value and optimal solutions converge to those of the full-information problem. The paper answers it for all data-driven DROs at once, with a single property of the test.

Setting

The support Ξ⊆Rd\Xi\subseteq\mathbb R^dΞ⊆Rd of ξ\xiξ is closed, and P(Ξ)\mathcal P(\Xi)P(Ξ) is the set of Borel probability distributions on Ξ\XiΞ with the topology of weak convergence. The data ξ1,ξ2,…\xi^1,\xi^2,\dotsξ1,ξ2,… are i.i.d. from an unknown F∈P(Ξ)F\in\mathcal P(\Xi)F∈P(Ξ). A data-driven uncertainty set (DUS) map assigns to each sample (ξ1,…,ξN)(\xi^1,\dots,\xi^N)(ξ1,…,ξN) a set FN⊆P(Ξ)\mathcal F_N\subseteq\mathcal P(\Xi)FN​⊆P(Ξ); it is the confidence region of the test that rejects a hypothesis F0F_0F0​ exactly when F0∉FNF_0\notin\mathcal F_NF0​∈/FN​. Decisions lie in a closed set X⊆RdxX\subseteq\mathbb R^{d_x}X⊆Rdx​, the cost is c(x;ξ)c(x;\xi)c(x;ξ), and the robust objective is the worst-case expectation

C(x;FN)=sup⁡F0∈FNEF0[c(x;ξ)].\mathcal C(x;\mathcal F_N)=\sup_{F_0\in\mathcal F_N}\mathbb E_{F_0}[c(x;\xi)] .C(x;FN​)=F0​∈FN​sup​EF0​​[c(x;ξ)].

A test is uniformly consistent (Definition 3) if, for every data-generating FFF, almost surely every sequence FNF_NFN​ that does not converge weakly to FFF lies outside FN\mathcal F_NFN​ for infinitely many NNN. The admissible instances are those satisfying three assumptions: (1) c(⋅;ξ)c(\cdot;\xi)c(⋅;ξ) is equicontinuous in xxx uniformly over ξ\xiξ; (2) XXX is bounded, or ccc is coercive in xxx uniformly over a set DDD of positive FFF-probability and bounded below off DDD; (3) Ξ\XiΞ is bounded, or there is a growth function ϕ≥0\phi\ge0ϕ≥0 with c(x;ξ)=O(ϕ(ξ))c(x;\xi)=O(\phi(\xi))c(x;ξ)=O(ϕ(ξ)) whose expectation the DUS tracks: sup⁡F0∈FN∣EF0ϕ−1N∑iϕ(ξi)∣→0\sup_{F_0\in\mathcal F_N}|\mathbb E_{F_0}\phi-\frac1N\sum_i\phi(\xi^i)|\to0supF0​∈FN​​∣EF0​​ϕ−N1​∑i​ϕ(ξi)∣→0 almost surely. The convergence conditions are (18) C(x;FN)→EF[c(x;ξ)]\mathcal C(x;\mathcal F_N)\to\mathbb E_F[c(x;\xi)]C(x;FN​)→EF​[c(x;ξ)] uniformly on compact subsets of XXX; (19) convergence of optimal values; and (20) every sequence of robust minimizers has a limit point, and all its limit points are true minimizers.

Formalization targets

Goal: Theorem 2

For a DUS map with nonempty sets,

FN is uniformly consistent  ⟺  [∀F,X,c satisfying Assumptions 1–3: (18),(19),(20) hold a.s.].\mathcal F_N\ \text{is uniformly consistent}\iff\Big[\forall F,X,c\ \text{satisfying Assumptions 1–3}:\ (18),(19),(20)\ \text{hold a.s.}\Big].FN​ is uniformly consistent⟺[∀F,X,c satisfying Assumptions 1–3: (18),(19),(20) hold a.s.].

The statement carries no rates and no constants; it is a characterization, so it does not depend on any particular test.

Milestones

  1. Proposition 8: under uniform consistency and Assumptions 1 and 3, almost surely EFN[c(x;ξ)]→EF[c(x;ξ)]\mathbb E_{F_N}[c(x;\xi)]\to\mathbb E_F[c(x;\xi)]EFN​​[c(x;ξ)]→EF​[c(x;ξ)] for every x∈Xx\in Xx∈X and every sequence FN∈FNF_N\in\mathcal F_NFN​∈FN​.
  2. Proposition 9: deterministically, convergence of C(xN;FN)\mathcal C(x_N;\mathcal F_N)C(xN​;FN​) along every convergent xN→xx_N\to xxN​→x implies (18).
  3. Theorem 2, "only if" side: uniform consistency implies (18)–(20) a.s. for every admissible instance.
  4. Theorem 2, "if" side: if the test is not uniformly consistent, some admissible instance fails (18)–(20) with positive probability.
  5. Theorem 3: if FN\mathcal F_NFN​ always contains the empirical distribution and the test is ccc-consistent (Definition 5, consistency with respect to the expected costs only), then (18)–(20) hold a.s.
  6. Proposition 4 (first sentence): a uniformly consistent test is consistent in the classical sense, P(F0∉FN)→1\mathbb P(F_0\notin\mathcal F_N)\to1P(F0​∈/FN​)→1 for each F0≠FF_0\neq FF0​=F.

Significance

The result. Theorem 2 reduces the asymptotic analysis of a data-driven DRO to a property of a statistical test, which can be checked once per test rather than once per optimization model. The paper uses it to classify DUSs: the χ2\chi^2χ2 and G-tests on a finite support and the Kolmogorov–Smirnov, Kuiper, Cramér–von Mises, Watson and Anderson–Darling tests on the line are uniformly consistent (Theorems 4–5, separate missions of this series), while the test of marginals is not even consistent and so cannot give convergence for every cost. The necessity direction makes the property sharp: a weaker test always admits a well-behaved cost on which the DRO fails to converge. Theorem 3 explains the remaining cases where a DRO converges for one given cost.

Formalizing it. The result is proved in the paper (§10.5), but its proof leans on several standard facts used without detail: the portmanteau theorem, uniform integrability arguments from Billingsley's Convergence of Probability Measures, the Lévy continuity theorem with the Cramér–Wold device, and coercivity arguments carried over from the paper's Theorem 1. A machine-checked version would fix exactly which regularity of the cost these steps need; three such conditions are stated explicitly here (see the scope section). No part of the result is formalized elsewhere to our knowledge; the platform has no notion of a data-driven uncertainty set or of uniform consistency.

Difficulty

Uniform consistency is a statement about weak convergence, while (18)–(20) are statements about expectations of an unbounded, xxx-dependent cost under the members of FN\mathcal F_NFN​. Weak convergence controls expectations only of bounded continuous functions, so the "only if" side must recover convergence of EFN[c(x;ξ)]\mathbb E_{F_N}[c(x;\xi)]EFN​​[c(x;ξ)] from Assumption 3, which requires uniform integrability, and then upgrade pointwise convergence to uniform convergence on compacts and to convergence of minimizers when XXX is unbounded. The tempting route through a single limiting minimization problem fails: the sets FN\mathcal F_NFN​ change with the sample, the worst case is a supremum over infinitely many distributions, and the minimizers of C(⋅;FN)\mathcal C(\cdot;\mathcal F_N)C(⋅;FN​) could escape to infinity unless coercivity on DDD transfers to every member of FN\mathcal F_NFN​. The "if" side needs, from the mere failure of uniform consistency, a single admissible cost whose worst-case objective detects every failure of weak convergence.

Formalization scope

Points are EuclideanSpace ℝ (Fin d), distributions are ProbabilityMeasure ↥Ξ (whose topology is weak convergence), and the data live on the canonical product space ΞN\Xi^{\mathbb N}ΞN with the product law, 0-indexed. "Almost surely" is ∀ᵐ, and in Definition 3 one null set serves every sequence FNF_NFN​. An expectation over a member of FN\mathcal F_NFN​ is +∞+\infty+∞ when the integrand is not integrable, never Lean's default 000. Minima in (19)–(20) are infima in the extended reals, so the empty and unattained cases need no special treatment; (20) applies to sequences that are minimizers for all large NNN.

Four hypotheses are added to the page, each used by the paper's own proof: every FN\mathcal F_NFN​ is nonempty (otherwise the worst case is −∞-\infty−∞ and uniform consistency is vacuous); c(x;⋅)c(x;\cdot)c(x;⋅) is continuous on Ξ\XiΞ; the set DDD of Assumption 2b is open; the growth function ϕ\phiϕ of Assumption 3b is continuous, with EFϕ<∞\mathbb E_F\phi<\inftyEF​ϕ<∞ as fixed on p. 10 of the paper. Theorem 3 additionally assumes Assumption 2, which its proof uses. Proposition 4 assumes the rejection events are measurable. A formalization in which the worst case uses the raw Bochner integral, the supremum ranges over a possibly empty set, or (20) requires minimizers at N=0N=0N=0 would make the statements false or vacuous, and is ruled out.

A full development needs weak convergence and the portmanteau theorem (in Mathlib), uniform integrability and its relation to convergence of means under weak convergence, the strong law of large numbers for the growth function and the empirical measure, the Lévy continuity theorem in several dimensions (in Mathlib, via characteristic functions), and a countable convergence-determining class of bounded continuous functions. The uniform-integrability step and the path characterization of local uniform convergence (Proposition 9) are reusable beyond this mission. Contributions to any milestone are welcome; Proposition 9 is self-contained and deterministic.

Selected references

  • D. Bertsimas, V. Gupta, N. Kallus, Robust Sample Average Approximation, arXiv:1408.4445v3, 2016; Mathematical Programming 171 (2018). https://arxiv.org/abs/1408.4445v3, https://doi.org/10.1007/s10107-017-1174-z
  • P. Billingsley, Convergence of Probability Measures, 2nd ed., Wiley, 1999. https://doi.org/10.1002/9780470316962
  • E. Delage, Y. Ye, Distributionally Robust Optimization Under Moment Uncertainty with Application to Data-Driven Problems, Operations Research 58(3), 2010. https://doi.org/10.1287/opre.1090.0741
  • A. Ben-Tal, D. den Hertog, A. De Waegenaere, B. Melenberg, G. Rennen, Robust Solutions of Optimization Problems Affected by Uncertain Probabilities, Management Science 59(2), 2013. https://doi.org/10.1287/mnsc.1120.1641
  • P. Mohajerin Esfahani, D. Kuhn, Data-Driven Distributionally Robust Optimization Using the Wasserstein Metric, Mathematical Programming 171, 2018. https://doi.org/10.1007/s10107-017-1172-1
9 thms1 active userReviewed
ProbabilityStatistics·Captain: mikedeng1

Robust Sample Average Approximation 3: The Kolmogorov–Smirnov, Kuiper, Cramér–von Mises, Watson and Anderson–Darling Tests with O(N^(−1/2)) Thresholds Are Uniformly ConsistentResearch Paper

Motivation

Robust Sample Average Approximation (Robust SAA) is a data-driven approach to stochastic optimization proposed by Bertsimas, Gupta and Kallus (arXiv:1408.4445; Mathematical Programming (2018), doi:10.1007/s10107-017-1174-z). Instead of minimizing the empirical average of a cost, a decision maker minimizes the worst-case expected cost over a distributional uncertainty set: the set of distributions that a statistical goodness-of-fit (GoF) test does not reject at a chosen significance level, given the observed data. Each classical GoF test therefore produces its own robust optimization model.

Whether such a model behaves well as data accumulate depends on a property of the test. Theorem 2 of the paper shows that Robust SAA converges, in objective, optimal value and optimal solutions, for every admissible cost if and only if the underlying test is uniformly consistent. This mission concerns the univariate case: it asks for a proof that the five standard tests on the real line, those of Kolmogorov–Smirnov, Kuiper, Cramér–von Mises, Watson and Anderson–Darling, are uniformly consistent (Theorem 5, p. 14). The companion missions of the series treat the characterization (Theorem 2) and the discrete χ² and G-tests (Theorem 4).

Setting

Let FFF be a probability distribution on R\mathbb RR, and let ξ1,ξ2,…\xi^1,\xi^2,\dotsξ1,ξ2,… be independent draws from FFF. As in §3.2 of the paper, ξ\xiξ is a continuous random variable: FFF has no atoms. For a candidate distribution F0F_0F0​ write F0(t)=F0((−∞,t])F_0(t)=F_0((-\infty,t])F0​(t)=F0​((−∞,t]) for its cdf, and write ξ(1)≤⋯≤ξ(N)\xi^{(1)}\le\dots\le\xi^{(N)}ξ(1)≤⋯≤ξ(N) for the sorted sample of size NNN. With ui=F0(ξ(i))u_i=F_0(\xi^{(i)})ui​=F0​(ξ(i)), the five statistics (8) are

DN=max⁡1≤i≤Nmax⁡{iN−ui, ui−i−1N},VN=max⁡i(ui−i−1N)+max⁡i(iN−ui),WN=(112N2+1N∑i=1N(2i−12N−ui)2)1/2,UN=(WN2−(1N∑i=1Nui−12)2)1/2,\begin{aligned} D_N&=\max_{1\le i\le N}\max\Big\{\tfrac iN-u_i,\ u_i-\tfrac{i-1}N\Big\},& V_N&=\max_i\Big(u_i-\tfrac{i-1}N\Big)+\max_i\Big(\tfrac iN-u_i\Big),\\ W_N&=\Big(\tfrac1{12N^2}+\tfrac1N\sum_{i=1}^N\Big(\tfrac{2i-1}{2N}-u_i\Big)^2\Big)^{1/2},& U_N&=\Big(W_N^2-\Big(\tfrac1N\sum_{i=1}^Nu_i-\tfrac12\Big)^2\Big)^{1/2}, \end{aligned}DN​WN​​=1≤i≤Nmax​max{Ni​−ui​, ui​−Ni−1​},=(12N21​+N1​i=1∑N​(2N2i−1​−ui​)2)1/2,​VN​UN​​=imax​(ui​−Ni−1​)+imax​(Ni​−ui​),=(WN2​−(N1​i=1∑N​ui​−21​)2)1/2,​ AN=(−1−∑i=1N2i−1N2(log⁡ui+log⁡(1−uN+1−i)))1/2.A_N=\Big(-1-\sum_{i=1}^N\tfrac{2i-1}{N^2}\big(\log u_i+\log(1-u_{N+1-i})\big)\Big)^{1/2}.AN​=(−1−i=1∑N​N22i−1​(logui​+log(1−uN+1−i​)))1/2.

For a statistic SNS_NSN​ and a threshold QSN(α)Q_{S_N}(\alpha)QSN​​(α), the confidence region (7) is FSNα={F0:SN(F0;ξ1,…,ξN)≤QSN(α)}\mathcal F^\alpha_{S_N}=\{F_0: S_N(F_0;\xi^1,\dots,\xi^N)\le Q_{S_N}(\alpha)\}FSN​α​={F0​:SN​(F0​;ξ1,…,ξN)≤QSN​​(α)}: the distributions the test does not reject. A test is uniformly consistent (Definition 3, p. 13) if, for every data-generating FFF, almost surely every sequence FNF_NFN​ that does not converge weakly to FFF is rejected infinitely often, i.e. FN∉FNF_N\notin\mathcal F_NFN​∈/FN​ for infinitely many NNN. The Lévy metric between cdfs is dLeˊvy(G,G′)=inf⁡{ϵ>0:G(ξ−ϵ)−ϵ≤G′(ξ)≤G(ξ+ϵ)+ϵ ∀ξ}d_{\text{Lévy}}(G,G')=\inf\{\epsilon>0: G(\xi-\epsilon)-\epsilon\le G'(\xi)\le G(\xi+\epsilon)+\epsilon\ \forall\xi\}dLeˊvy​(G,G′)=inf{ϵ>0:G(ξ−ϵ)−ϵ≤G′(ξ)≤G(ξ+ϵ)+ϵ ∀ξ}, and F^N\hat F_NF^N​ denotes the empirical cdf.

Formalization targets

Goal: Theorem 5

For thresholds τN\tau_NτN​ with τN≤q/N\tau_N\le q/\sqrt NτN​≤q/N​ for some constant qqq and all N≥1N\ge1N≥1,

the regions {F0:SN(F0)≤τN},SN∈{DN,VN,WN,UN,AN}, are uniformly consistent.\text{the regions }\{F_0: S_N(F_0)\le\tau_N\},\quad S_N\in\{D_N,V_N,W_N,U_N,A_N\},\ \text{are uniformly consistent.}the regions {F0​:SN​(F0​)≤τN​},SN​∈{DN​,VN​,WN​,UN​,AN​}, are uniformly consistent.

The statement fixes only the rate O(N−1/2)O(N^{-1/2})O(N−1/2) of the thresholds, which is all the paper's argument uses; it holds for any sequence with that rate, in particular for the tabulated quantiles QSN(α)Q_{S_N}(\alpha)QSN​​(α).

Milestones (§10.7, pp. 38–39)

  1. dLeˊvyd_{\text{Lévy}}dLeˊvy​ metrizes weak convergence on R\mathbb RR; almost surely dLeˊvy(F^N,F)→0d_{\text{Lévy}}(\hat F_N,F)\to0dLeˊvy​(F^N​,F)→0.
  2. dLeˊvy(F^N,F0)≤DN(F0)d_{\text{Lévy}}(\hat F_N,F_0)\le D_N(F_0)dLeˊvy​(F^N​,F0​)≤DN​(F0​), and DN≤VND_N\le V_NDN​≤VN​.
  3. DN2(F0)≤max⁡{1/N+3/(2N), N WN2(F0)+2/N}D_N^2(F_0)\le\max\{1/\sqrt N+3/(2N),\ \sqrt N\,W_N^2(F_0)+2/\sqrt N\}DN2​(F0​)≤max{1/N​+3/(2N), N​WN2​(F0​)+2/N​}, and the same with WN2W_N^2WN2​ replaced by AN2A_N^2AN2​.
  4. WN2−UN2=(1N∑iF0(ξ(i))−12)2≤max⁡{… }W_N^2-U_N^2=(\tfrac1N\sum_iF_0(\xi^{(i)})-\tfrac12)^2\le\max\{\dots\}WN2​−UN2​=(N1​∑i​F0​(ξ(i))−21​)2≤max{…}, a bound in terms of DN′(F0)=max⁡i∣F0(ξ(i))−2i−12N∣D'_N(F_0)=\max_i|F_0(\xi^{(i)})-\tfrac{2i-1}{2N}|DN′​(F0​)=maxi​∣F0​(ξ(i))−2N2i−1​∣.

Significance

The result. Together with Theorem 2 of the paper, Theorem 5 guarantees that Robust SAA built on any of these five univariate tests converges for every cost that is equicontinuous in the decision and satisfies the paper's growth conditions. It also feeds Proposition 5, which extends the guarantee to multivariate problems with separable costs by testing marginals. Without it, the univariate tests would carry only the finite-sample guarantee of Proposition 1 and no asymptotic one.

The formalization. The theorem is proved in the paper, but the published argument for the Watson test has a gap. Its last step asserts that (1N∑imin⁡{1,2i−12N+DN′}−12)2=O(1/N)(\tfrac1N\sum_i\min\{1,\tfrac{2i-1}{2N}+D'_N\}-\tfrac12)^2=O(1/N)(N1​∑i​min{1,2N2i−1​+DN′​}−21​)2=O(1/N) whenever DN′≥1/N+1/ND'_N\ge1/\sqrt N+1/NDN′​≥1/N​+1/N, which fails when DN′D'_NDN′​ stays bounded away from 000 (at DN′=1/2D'_N=1/2DN′​=1/2 the expression tends to 9/649/649/64). The Watson clause is believed true, but a complete proof needs a different argument. The Anderson–Darling step relies on a cited integral representation rather than an argument on the page. A machine-checked proof would close both points. No part of the theorem is formalized elsewhere, as far as a search of the platform shows.

Difficulty

The Kolmogorov–Smirnov case reduces to deterministic inequalities once the Lévy-metric picture is in place, but each of the other statistics controls DND_NDN​ only indirectly. The Cramér–von Mises statistic is an L2L^2L2 quantity, and an L2L^2L2 bound does not by itself control a supremum. What closes the gap is monotonicity of F0F_0F0​ together with the rate of the threshold, so a threshold that merely tends to 000 is not enough. The Watson statistic is invariant under cyclic rotation of the transformed sample u1,…,uNu_1,\dots,u_Nu1​,…,uN​ on the circle [0,1)[0,1)[0,1), so a small UNU_NUN​ leaves the centring term 1N∑iF0(ξ(i))−12\tfrac1N\sum_iF_0(\xi^{(i)})-\tfrac12N1​∑i​F0​(ξ(i))−21​ uncontrolled until the boundary constraints 0≤F0≤10\le F_0\le10≤F0​≤1 are used. The obvious attempt, bounding the centring term by DN′D'_NDN′​ alone, is the step that fails on p. 39.

Formalization scope

Distributions are MeasureTheory.ProbabilityMeasure ℝ, with weak convergence as their topology. The cdf is Mathlib's ProbabilityTheory.cdf. The data are the coordinates of the canonical product space ℕ → ℝ under Measure.infinitePi, so the sample of size NNN is (ω0,…,ωN−1)(\omega_0,\dots,\omega_{N-1})(ω0​,…,ωN−1​). Order statistics are the sorted sample s ∘ Tuple.sort s, and every formula of (8) is written with the paper's 1-based index (iii in Lean is the paper's i+1i+1i+1; ξ(N+1−i)\xi^{(N+1-i)}ξ(N+1−i) is Fin.rev i).

Committed conventions:

  • distributions range over all of R\mathbb RR; a support [ξ‾,ξ‾][\underline\xi,\overline\xi][ξ​,ξ​] with infinite endpoints is this case;
  • the data law has no atoms (§3.2's standing assumption, expressed by the predicate IsContinuousLaw), although the paper's proof does not use it;
  • ANA_NAN​ takes values in [0,∞][0,\infty][0,∞] and is +∞+\infty+∞ when some F0(ξ(i))∈{0,1}F_0(\xi^{(i)})\in\{0,1\}F0​(ξ(i))∈{0,1}, so Lean's log⁡0=0\log0=0log0=0 never makes a distribution look close;
  • the threshold hypothesis is an explicit constant, τN≤q/N\tau_N\le q/\sqrt NτN​≤q/N​ for N≥1N\ge1N≥1;
  • maxima over an empty index set are 000, and every fixed-NNN milestone assumes N≥1N\ge1N≥1.

Uniform consistency is quantified over all sequences FNF_NFN​ inside the almost-sure event, as in Definition 3. A pointwise consistency statement for a fixed F0F_0F0​ is a different and weaker theorem and does not count.

A full development needs the Lévy metric and its metrization of weak convergence on R\mathbb RR, the Glivenko–Cantelli theorem in Lévy form, and elementary order-statistic inequalities. The first two are reusable well beyond this mission. Contributions of independent value include a corrected proof of the Watson step, and a proof of the goal without the atomless hypothesis.

Selected references

  • D. Bertsimas, V. Gupta, N. Kallus, Robust Sample Average Approximation, arXiv:1408.4445v3 (2016); Mathematical Programming (2018). https://arxiv.org/abs/1408.4445, https://doi.org/10.1007/s10107-017-1174-z
  • R. M. Dudley, Real Analysis and Probability, Cambridge University Press, 2002, Theorem 11.4.1. https://doi.org/10.1017/CBO9780511755347
  • R. B. D'Agostino, M. A. Stephens (eds.), Goodness-of-Fit Techniques, Marcel Dekker, New York, 1986 (reference [14] of the paper; the source of the formulas (8)).
9 thms1 active userReviewed
ProbabilityStatistics·Captain: mikedeng1

Robust Sample Average Approximation 2: Pearson's χ² Test and the G-Test on a Known Finite Support, with Thresholds √(Q/N), Are Uniformly ConsistentResearch Paper

Motivation

A decision made from sampled data can be sensitive to the distribution chosen to represent future uncertainty. Robust sample average approximation addresses this by optimizing against every distribution that passes a goodness of fit test, rather than against one fitted distribution. The resulting set of plausible laws is a data driven uncertainty set. Its finite sample coverage and its behavior as the sample grows serve different purposes: coverage protects against sampling error at a chosen sample size, while convergence asks whether the uncertainty set eventually stops admitting laws far from the true one. Bertsimas, Gupta, and Kallus distinguish these properties in their Robust Sample Average Approximation paper.

This mission concerns the two tests the paper assigns to a known finite support: Pearson's χ2\chi^2χ2 test and the G test. Both compare the empirical frequencies of the observed support points with the probabilities under a proposed law. Their confidence regions give concrete uncertainty sets for discrete data. The paper's Theorem 4 states that both tests have the stronger convergence property it calls uniform consistency. That property is used in its general characterization of when robust sample average approximation converges for every admissible cost function (Theorem 2 and Table 1, pp. 13–14).

Setting

Fix a known set of nnn possible observations, labelled j=1,…,nj=1,\ldots,nj=1,…,n. A probability law FFF on this support is specified by probabilities pF(j)p_F(j)pF​(j), each nonnegative and summing to one. A sample of size NNN consists of independent observations ξ1,…,ξN\xi^1,\ldots,\xi^Nξ1,…,ξN from an unknown true law FFF. Its empirical frequency at label jjj is

p^N(j)=1N∑i=1N1{ξi=j}.\widehat p_N(j)=\frac1N\sum_{i=1}^N\mathbf 1\{\xi^i=j\}.p​N​(j)=N1​i=1∑N​1{ξi=j}.

To assess a hypothetical law F0F_0F0​ with probabilities p0(j)p_0(j)p0​(j), Pearson's statistic and the G statistic are, in the paper's normalization,

XN(F0)=(∑j=1n(p0(j)−p^N(j))2p0(j))1/2,GN(F0)=(2∑j=1np^N(j)log⁡p^N(j)p0(j))1/2.X_N(F_0)=\left(\sum_{j=1}^n\frac{(p_0(j)-\widehat p_N(j))^2}{p_0(j)}\right)^{1/2},\qquad G_N(F_0)=\left(2\sum_{j=1}^n\widehat p_N(j)\log\frac{\widehat p_N(j)}{p_0(j)}\right)^{1/2}.XN​(F0​)=(j=1∑n​p0​(j)(p0​(j)−p​N​(j))2​)1/2,GN​(F0​)=(2j=1∑n​p​N​(j)logp0​(j)p​N​(j)​)1/2.

The confidence region for each test is the set of hypothetical laws whose statistic is at most Q/N\sqrt{Q/N}Q/N​, where Q≥0Q\ge0Q≥0 is fixed across sample sizes. An observed point to which F0F_0F0​ assigns zero probability makes either statistic infinite. A term with zero empirical frequency contributes zero to the G sum, including when its hypothetical probability is zero. These conventions make the regions describe the same tests at the boundary of the probability simplex as they do in its interior (§3.1, p. 8; §10.6, p. 38).

A sequence of laws FNF_NFN​ is weakly convergent to FFF when expectations of every bounded continuous function converge. On this finite support, weak convergence can be measured by the total variation distance

dTV(q,q′)=12∑j=1n∣q(j)−q′(j)∣.d_{\mathrm{TV}}(q,q')=\frac12\sum_{j=1}^n|q(j)-q'(j)|.dTV​(q,q′)=21​j=1∑n​∣q(j)−q′(j)∣.

The paper's Definition 3 calls a test uniformly consistent when, for every true law FFF, almost surely every sequence FNF_NFN​ that fails to converge weakly to FFF lies outside the confidence region infinitely often (Definition 3, p. 13). The almost sure event applies to all such sequences at once; this is stronger than a statement about one fixed false hypothesis.

Formalization targets

Goal: both tests are uniformly consistent

For the Pearson and G confidence regions FNχ2\mathcal F_N^{\chi^2}FNχ2​ and FNG\mathcal F_N^GFNG​ with fixed Q≥0Q\ge0Q≥0, the target is

UC⁡(FNχ2)andUC⁡(FNG).\operatorname{UC}(\mathcal F_N^{\chi^2})\quad\text{and}\quad\operatorname{UC}(\mathcal F_N^G).UC(FNχ2​)andUC(FNG​).

Here UC⁡\operatorname{UC}UC has the quantifier order of Definition 3: every generating FFF, almost every infinite sample, and every nonconvergent sequence FNF_NFN​. This is Theorem 4, p. 14. Its four milestone claims are the finite support description of weak convergence, almost sure convergence of empirical frequencies, and the paper's total variation bounds for XNX_NXN​ and GNG_NGN​ (§10.6, p. 38).

Significance

The result says that distributions retained by either test cannot remain persistently separated from the data generating law as the sample grows. In the paper's framework, that is the statistical condition needed for its general theorem to turn test based confidence regions into convergent robust optimization values and decisions under its separate assumptions on the cost and feasible set (Theorem 2, p. 13). Theorem 4 itself concerns the tests; it does not assume a cost function or assert optimization convergence directly.

The paper proves Theorem 4, so the open work here is its Lean formalization. A complete development would also leave reusable statements for finite probability vectors, their empirical frequencies under an infinite product law, total variation, and the boundary aware Pearson and G divergences. These pieces can support other discrete goodness of fit and distributionally robust optimization results. The mission does not claim that Theorem 4 already has a machine checked proof.

Difficulty

Pointwise convergence of empirical frequencies does not by itself state uniform consistency. The latter quantifies over all moving hypothetical laws FNF_NFN​ on one almost sure sample event, including laws chosen differently at every NNN. A direct argument about one fixed F0≠FF_0\ne FF0​=F therefore does not address the target. The finite support makes weak convergence concrete, but the formal statement must still relate that topology to a finite total variation sum and keep the quantifiers in the order required by Definition 3.

The boundary of the probability simplex is another source of difficulty. Lean's real division and logarithm have default values at zero. Using those defaults in the test statistics would accept some laws that assign zero probability to an actually observed category. Both statistics must instead retain the paper's infinite penalty there. The G expression also contains individual logarithmic terms that may be negative before summation, so its extended value cannot be built by coercing each term directly to a nonnegative extended real.

Formalization scope

The known support is represented by Fin n. Its members are labels for the paper's distinct points ξ^1,…,ξ^n\hat\xi^1,\ldots,\hat\xi^nξ^​1,…,ξ^​n; coordinates in Euclidean space are irrelevant to these tests. Laws use Mathlib's ProbabilityMeasure (Fin n) and its weak topology. The data are coordinates of the canonical infinite product law, with Lean index zero corresponding to the paper's first observation. The empty support has no probability law; all substantive instances have n≥1n\ge1n≥1. The empty sample N=0N=0N=0 is defined by total operations but has no bearing on an eventual or infinitely often assertion.

Pearson's statistic lives in [0,+∞][0,+\infty][0,+∞]. The squared G statistic uses extended real arithmetic so finite negative summands can be combined with an infinite boundary term. Its region compares GN2G_N^2GN2​ with Q/NQ/NQ/N, equivalent to the paper's GN≤Q/NG_N\le\sqrt{Q/N}GN​≤Q/N​ at positive NNN. The free Q≥0Q\ge0Q≥0 stands for a fixed threshold constant, including the paper's chi square quantile; no particular quantile value is built into the theorem. This makes explicit the only property of that threshold used in §10.6: independence from NNN.

Contributions may prove the four milestone statements, connect Mathlib's finite probability and strong law interfaces to the empirical frequencies, or establish finite support versions of Cauchy–Schwarz and Pinsker with the correct infinite boundary cases. A proof of the goal must use the stated regions and preserve simultaneous quantification over every moving sequence of laws; replacing uniform consistency by convergence of p^N\widehat p_Np​N​ would remove its main claim.

Selected references

  • Dimitris Bertsimas, Vishal Gupta, Nathan Kallus, Robust Sample Average Approximation, arXiv:1408.4445v3, 2016; published in Mathematical Programming (2018). Preprint · Publication DOI.
6 thms1 active userReviewed
Algorithmic Game TheoryProbability·Captain: mikedeng1

Posted Price Mechanisms and Optimal Threshold Strategies for Random Arrivals II: For I.I.D. Variables, Adaptive Thresholds Achieve E[max Xᵢ] ≤ β*·E[X_t] with β* ≈ 1.341Research Paper

Motivation

A prophet inequality compares a gambler who sees nonnegative random rewards X1,…,XnX_1,\ldots,X_nX1​,…,Xn​ one at a time, and must accept or reject each on arrival, with a prophet who knows all values in advance and simply takes max⁡iXi\max_i X_imaxi​Xi​. For independent and identically distributed (i.i.d.) rewards the question is the best constant β\betaβ such that some online rule always earns at least E(max⁡iXi)/βE(\max_i X_i)/\betaE(maxi​Xi​)/β. The question matters in operations research and mechanism design because a seller who posts take-it-or-leave-it prices to arriving customers faces exactly this stopping problem on the customers' virtual values; prophet inequalities translate into revenue guarantees for posted price mechanisms relative to Myerson's optimal auction (Chawla, Hartline, Malec and Sivan 2010).

Timeline (as recounted in §1 of the paper):

  • Hill and Kertz (1982) characterized, through a recursion, the smallest constants ana_nan​ with E(max⁡iXi)≤ansup⁡tE(Xt)E(\max_iX_i)\le a_n\sup_tE(X_t)E(maxi​Xi​)≤an​supt​E(Xt​) for nnn i.i.d. nonnegative variables, showed 1.1<an<1.61.1<a_n<1.61.1<an​<1.6, gave instances where ana_nan​ cannot be beaten, and conjectured that (an)(a_n)(an​) is monotone; the computation of its limit was left open.
  • Kertz (1986) proved an→β∗≈1.341a_n\to\beta^*\approx1.341an​→β∗≈1.341, the solution of the integral equation (2) below, and conjectured that β∗\beta^*β∗ bounds every ana_nan​; the best uniform bound then known was an≤e/(e−1)≈1.582a_n\le e/(e-1)\approx1.582an​≤e/(e−1)≈1.582.
  • Abolhassani et al. (2017) improved the uniform bound to 1/0.738≈1.3551/0.738\approx1.3551/0.738≈1.355.
  • Correa, Foncea, Hoeksma, Oosterwijk and Vredeveld (EC 2017; Math. Oper. Res. 2021) proved E(max⁡iXi)≤β∗E(Xt)E(\max_iX_i)\le\beta^*E(X_t)E(maxi​Xi​)≤β∗E(Xt​) for a threshold rule ttt, so an≤β∗a_n\le\beta^*an​≤β∗ for all nnn and 1/β∗≈0.7451/\beta^*\approx0.7451/β∗≈0.745 is the tight i.i.d. guarantee.

Setting

Let X1,…,XnX_1,\ldots,X_nX1​,…,Xn​ be i.i.d. nonnegative random variables with distribution function FFF. A threshold rule fixes real numbers τ1,…,τn\tau_1,\ldots,\tau_nτ1​,…,τn​ and stops at

t=min⁡{i∈{1,…,n}:Xi≥τi},t=\min\{i\in\{1,\ldots,n\}: X_i\ge\tau_i\},t=min{i∈{1,…,n}:Xi​≥τi​},

receiving XtX_tXt​. Because the thresholds may change with the step iii, the rule is adaptive; for i.i.d. variables the number of variables already seen is all the history the rule needs.

The analysis works with quantiles. The generalized inverse is F−1(q)=inf⁡{x≥0:F(x)≥q}F^{-1}(q)=\inf\{x\ge0: F(x)\ge q\}F−1(q)=inf{x≥0:F(x)≥q}, and τ(q)=F−1(1−q)\tau(q)=F^{-1}(1-q)τ(q)=F−1(1−q) is the threshold that accepts with probability qqq (with a coin toss to break ties at an atom). The function

R(q)=∫0qF−1(1−θ) dθR(q)=\int_0^qF^{-1}(1-\theta)\,d\thetaR(q)=∫0q​F−1(1−θ)dθ

is the expected reward from a variable accepted with probability qqq. The quantile stopping rule of Algorithm 2 partitions [0,1][0,1][0,1] into intervals Ai=[εi−1,εi]A_i=[\varepsilon_{i-1},\varepsilon_i]Ai​=[εi−1​,εi​], draws an acceptance probability qi∈Aiq_i\in A_iqi​∈Ai​ with density proportional to ψ(q)=(n−1)(1−q)n−2\psi(q)=(n-1)(1-q)^{n-2}ψ(q)=(n−1)(1−q)n−2, and stops at step iii with probability qiq_iqi​; its value is E(Xr)E(X_r)E(Xr​). The analysis then passes through the normalizations γi=∫Aiψ\gamma_i=\int_{A_i}\psiγi​=∫Ai​​ψ and ρi\rho_iρi​, the recursion (8)–(9) on xi=1−εix_i=1-\varepsilon_ixi​=1−εi​, and the differential equation

y′(t)=y(t)(ln⁡y(t)−1)−(β−1),y(0)=1.(ODE)y'(t)=y(t)\big(\ln y(t)-1\big)-(\beta-1),\qquad y(0)=1.\qquad\text{(ODE)}y′(t)=y(t)(lny(t)−1)−(β−1),y(0)=1.(ODE)

Formalization targets

Goal: Theorem 2

Let β∗>1\beta^*>1β∗>1 be the unique solution of

∫01dyy(1−ln⁡y)+(β−1)=1.(2)\int_0^1\frac{dy}{y(1-\ln y)+(\beta-1)}=1.\qquad(2)∫01​y(1−lny)+(β−1)dy​=1.(2)

For every n≥1n\ge1n≥1 and every law of nonnegative i.i.d. X1,…,XnX_1,\ldots,X_nX1​,…,Xn​ there are thresholds τ1,…,τn\tau_1,\ldots,\tau_nτ1​,…,τn​ with

E(max⁡{X1,…,Xn})≤β∗E(Xt).E(\max\{X_1,\ldots,X_n\})\le\beta^*E(X_t).E(max{X1​,…,Xn​})≤β∗E(Xt​).

Milestones

In order: the quantile sandwich P(X≥τ(q))≥q≥P(X>τ(q))P(X\ge\tau(q))\ge q\ge P(X>\tau(q))P(X≥τ(q))≥q≥P(X>τ(q)); the representation (7) E(max⁡iXi)=n∫01ψ(q)R(q) dqE(\max_iX_i)=n\int_0^1\psi(q)R(q)\,dqE(maxi​Xi​)=n∫01​ψ(q)R(q)dq; the identity "expected reward of the quantile rule at qqq equals R(q)R(q)R(q)"; Lemma 4 (E(Xr)=∑iρi∫AiψRE(X_r)=\sum_i\rho_i\int_{A_i}\psi RE(Xr​)=∑i​ρi​∫Ai​​ψR); Lemma 5 (E(max⁡iXi)=nγ1E(Xr)E(\max_iX_i)=n\gamma_1E(X_r)E(maxi​Xi​)=nγ1​E(Xr​) when all ρi\rho_iρi​ agree); the equivalence of equal ρi\rho_iρi​ with the recursion (8) and γ1=1−x1n−1\gamma_1=1-x_1^{n-1}γ1​=1−x1n−1​; the closed form (9)–(10); Lemma 6 on (ODE); Proposition D.1; Lemma 7 (xin−1<y(i/n)x_i^{n-1}<y(i/n)xin−1​<y(i/n)); uniqueness of β∗\beta^*β∗ in (2); optimality of the dynamic-programming thresholds; domination of the randomized rule by them; and n(1−x1n−1)≤β∗n(1-x_1^{n-1})\le\beta^*n(1−x1n−1​)≤β∗.

A companion statement records the numerical guarantee 1/β∗>0.7451/\beta^*>0.7451/β∗>0.745.

Significance

The result. Theorem 2 gives the optimal constant of the i.i.d. prophet inequality and, through an=n(1−x1n−1)a_n=n(1-x_1^{n-1})an​=n(1−x1n−1​), shows an≤β∗a_n\le\beta^*an​≤β∗ for every nnn, which settles the problem posed by Hill and Kertz in 1982. Since an→β∗a_n\to\beta^*an​→β∗ (Kertz 1986), no rule of any kind achieves a better constant uniformly in nnn. Through the virtual-value connection it implies that, for customers with i.i.d. valuations, an adaptive posted price mechanism earns at least a 1/β∗>0.7451/\beta^*>0.7451/β∗>0.745 fraction of the optimal auction's revenue (Corollary 2 of the paper).

Formalizing it. The theorem is proved on paper; no machine-checked proof of it or of any of its lemmas is known. The proof combines a measure-theoretic part (quantile functions, randomized tie-breaking at atoms, Tonelli and integration by parts), a combinatorial recursion, and real analysis of an autonomous ODE with a third-derivative Taylor argument. Formal versions of these parts are reusable: the quantile representation of E(max⁡iXi)E(\max_iX_i)E(maxi​Xi​), the expected reward of a quantile rule, and the optimality of dynamic-programming thresholds are standard tools across prophet inequalities and optimal stopping.

Difficulty

The obvious attempt, computing the optimal dynamic-programming thresholds and bounding their value against E(max⁡iXi)E(\max_iX_i)E(maxi​Xi​), does not lead to a closed-form constant: the optimal values satisfy a nonlinear recursion that depends on the whole distribution. The paper avoids it by analysing a randomized rule whose value is exactly a ψ\psiψ-weighted average of RRR, which reduces the problem to a distribution-free recursion. Two steps then carry the weight. First, randomizing correctly at atoms of FFF so that every step accepts with probability exactly qqq; without this, R(q)R(q)R(q) is not the step reward. Second, comparing the discrete recursion (9) with the continuous (ODE): the comparison needs strict convexity and y′′′>0y'''>0y′′′>0, and the latter holds only because β>1.25\beta>1.25β>1.25, a standing assumption validated at β∗\beta^*β∗ at the end. The existence of the partition with all ρi\rho_iρi​ equal (a solution of (8) with x0=1x_0=1x0​=1, xn=0x_n=0xn​=0, xi∈(0,1)x_i\in(0,1)xi​∈(0,1)) is taken from Hill and Kertz and is not proved in the paper.

Formalization scope

Lean conventions:

  • The common law is a probability measure μ\muμ on R\mathbb RR with μ((−∞,0))=0\mu((-\infty,0))=0μ((−∞,0))=0; the sample is the coordinate vector of SamuelCahnProphet.IID.iidLaw μ n (the published i.i.d. setting), and E(max⁡iXi)E(\max_iX_i)E(maxi​Xi​) is its Emax. Coordinates are 0-based; the partition εi\varepsilon_iεi​, γi\gamma_iγi​, ρi\rho_iρi​ and the recursion keep the paper's 1-based index.
  • All expectations of the random variables, R(q)R(q)R(q) and E(Xr)E(X_r)E(Xr​) are lower Lebesgue integrals in [0,∞][0,\infty][0,∞]; no integrability is assumed in the goal, so infinite means are covered. Only bounded deterministic integrals (γi\gamma_iγi​, ρi\rho_iρi​, (2)) are interval integrals.
  • XtX_tXt​ is 000 when no XiX_iXi​ reaches its threshold (the page leaves it undefined; 000 is the least favourable choice).
  • The quantile rule's coin at an atom is an independent uniform variable on [0,1][0,1][0,1]; acceptance probabilities are drawn independently with density ψ/γi\psi/\gamma_iψ/γi​ on AiA_iAi​.
  • β∗\beta^*β∗ enters through (2) and β∗>1\beta^*>1β∗>1, not through a choice function; uniqueness is a milestone.
  • Statements involving ψ\psiψ, γi\gamma_iγi​, ρi\rho_iρi​ or the recursion assume n≥2n\ge2n≥2 (for n=1n=1n=1, ψ≡0\psi\equiv0ψ≡0). Lemmas 6 and 7 carry §4's standing assumption β>1.25\beta>1.25β>1.25; Proposition D.1 carries β>1.25\beta>1.25β>1.25. Lemma 6 does not assert existence of a solution of (ODE) on [0,1][0,1][0,1], which fails for β>β∗\beta>\beta^*β>β∗.
  • The dynamic-programming milestones assume EX<∞E X<\inftyEX<∞ and use the optimal-stopping recurrence Vi=E[max⁡(X,Vi+1)]V_i=E[\max(X,V_{i+1})]Vi​=E[max(X,Vi+1​)]; the page prints Vi=E(X∣X≥τi)V_i=E(X\mid X\ge\tau_i)Vi​=E(X∣X≥τi​), which is a slip.
  • The page's "β∗≈1.341>1/0.745\beta^*\approx1.341>1/0.745β∗≈1.341>1/0.745" has its inequality reversed and is not formalized.

A degenerate β\betaβ cannot trivialize the goal: β\betaβ is pinned by (2) and β>1\beta>1β>1, and the goal is stated for every n≥1n\ge1n≥1 and every nonnegative law, including laws with infinite mean.

A complete development needs the quantile function of a law on R\mathbb RR, Tonelli's theorem on product measures, uniqueness for ODEs with locally Lipschitz right-hand side (ODE_solution_unique), and Taylor's theorem; the existence of the Hill–Kertz partition is a further proof obligation. Contributions to any milestone, and reusable lemmas on quantile functions and finite-horizon optimal stopping, are welcome.

Selected references

  • J. Correa, P. Foncea, R. Hoeksma, T. Oosterwijk, T. Vredeveld, Posted price mechanisms and optimal threshold strategies for random arrivals, Mathematics of Operations Research 46(4):1452–1478, 2021. https://doi.org/10.1287/moor.2020.1105
  • T. P. Hill, R. P. Kertz, Comparisons of stop rule and supremum expectations of i.i.d. random variables, Annals of Probability 10(2):336–345, 1982. https://doi.org/10.1214/aop/1176993861
  • R. P. Kertz, Stop rule and supremum expectations of i.i.d. random variables: a complete comparison by conjugate duality, Journal of Multivariate Analysis 19(1):88–112, 1986. https://doi.org/10.1016/0047-259X(86)90095-3
  • M. Abolhassani, S. Ehsani, H. Esfandiari, M. T. HajiAghayi, R. Kleinberg, B. Lucier, Beating 1−1/e for ordered prophets, Proc. 49th ACM STOC, 61–71, 2017. https://doi.org/10.1145/3055399.3055479
  • S. Chawla, J. D. Hartline, D. L. Malec, B. Sivan, Multi-parameter mechanism design and sequential posted pricing, Proc. 42nd ACM STOC, 311–320, 2010. https://doi.org/10.1145/1806689.1806733
17 thms1 active userReviewed
Mechanism DesignProbability·Captain: mikedeng1

Pricing and Matching with Forward-Looking Buyers and Sellers 3: Every Incentive-Compatible Two-Sided Dynamic Mechanism Earns at Most the Expected Virtual-Surplus Assignment ValueResearch Paper

Motivation

Ride-hailing, freelancing and delivery platforms are intermediaries that post a price to buyers and a price to sellers and then match the two sides. Both sides are forward-looking: an agent who arrives can wait before requesting, and an agent who requests waits until she is matched, at a cost per unit of time. Chen and Hu, Pricing and Matching with Forward-Looking Buyers and Sellers (SSRN 2859864; Manufacturing & Service Operations Management 22(4), 2020, doi:10.1287/msom.2018.0769), measure simple pricing and matching policies against an upper bound that holds for every policy. This mission formalizes that upper bound in its mechanism-design form: no incentive-compatible two-sided dynamic mechanism earns more than the expected value of a clairvoyant assignment problem in virtual values.

The bound extends the virtual-value approach of Myerson (1981) and its bilateral-trade version by Myerson and Satterthwaite (1983) from a single static trade to a market with random arrivals on both sides, waiting costs and a timing constraint that matched units change hands at the same instant.

Setting

Over a horizon [0,T][0,T][0,T], buyers arrive as a Poisson process of rate λd>0\lambda^d>0λd>0. A buyer has type ϕ=(tϕ,vϕ)\phi=(t_\phi,v_\phi)ϕ=(tϕ​,vϕ​): arrival time tϕt_\phitϕ​ and valuation vϕ∈[v‾,vˉ]v_\phi\in[\underline v,\bar v]vϕ​∈[v​,vˉ], drawn with continuous positive density fdf^dfd and c.d.f. FdF^dFd, independently of the arrival time. Sellers arrive as an independent Poisson process of rate λs>0\lambda^s>0λs>0, with types ψ=(tψ,cψ)\psi=(t_\psi,c_\psi)ψ=(tψ​,cψ​) and costs cψ∈[c‾,cˉ]c_\psi\in[\underline c,\bar c]cψ​∈[c​,cˉ] of continuous positive density fsf^sfs and c.d.f. FsF^sFs. The paper assumes c‾≤vˉ\underline c\le\bar vc​≤vˉ and that the virtual value Vd(v)=v−(1−Fd(v))/fd(v)V^d(v)=v-(1-F^d(v))/f^d(v)Vd(v)=v−(1−Fd(v))/fd(v) and the virtual cost Vs(c)=c+Fs(c)/fs(c)V^s(c)=c+F^s(c)/f^s(c)Vs(c)=c+Fs(c)/fs(c) are increasing (Assumptions 1–2). Write HTH^THT for the realized arrivals.

A direct mechanism yTy^TyT assigns to every buyer an outcome (sϕ,mϕ,pϕ)(s_\phi,m_\phi,p_\phi)(sϕ​,mϕ​,pϕ​): the time sϕ∈[tϕ,T]s_\phi\in[t_\phi,T]sϕ​∈[tϕ​,T] at which she leaves, the indicator mϕ∈{0,1}m_\phi\in\{0,1\}mϕ​∈{0,1} that she receives a unit, and her payment pϕp_\phipϕ​. Every seller receives (sψ,mψ,pψ)(s_\psi,m_\psi,p_\psi)(sψ​,mψ​,pψ​) likewise, pψp_\psipψ​ now paid to her. The balancing condition requires that at every t∈[0,T]t\in[0,T]t∈[0,T] the number of buyers who receive a unit at ttt equals the number of sellers who deliver one at ttt. The buyer's utility is Ud=vϕmϕ−pϕ−b(sϕ−tϕ)U^d=v_\phi m_\phi-p_\phi-b(s_\phi-t_\phi)Ud=vϕ​mϕ​−pϕ​−b(sϕ​−tϕ​) and the seller's is Us=pψ−cψmψ−h(sψ−tψ)U^s=p_\psi-c_\psi m_\psi-h(s_\psi-t_\psi)Us=pψ​−cψ​mψ​−h(sψ​−tψ​), with waiting costs b,h≥0b,h\ge0b,h≥0. The intermediary's profit is

Π(yT)=∑ϕ∈HTpϕ−∑ψ∈HTpψ.\Pi(y^T)=\sum_{\phi\in H^T}p_\phi-\sum_{\psi\in H^T}p_\psi .Π(yT)=ϕ∈HT∑​pϕ​−ψ∈HT∑​pψ​.

Incentive compatibility (ICd) asks that no buyer gain, in the interim expectation E−ϕE_{-\phi}E−ϕ​ over all other agents, by reporting another valuation or a later arrival time tϕ^∈[tϕ,T]t_{\hat\phi}\in[t_\phi,T]tϕ^​​∈[tϕ​,T]; individual rationality (IRd) asks that her interim utility be nonnegative. (ICs) and (IRs) are the seller counterparts. Problem (B′) maximizes E[Π(yT)]E[\Pi(y^T)]E[Π(yT)] over feasible mechanisms subject to these constraints; its value is J∗J^*J∗.

The comparison object is the clairvoyant assignment problem (B): knowing HTH^THT in advance, choose a partial matching xxx of buyers to sellers maximizing

∑ϕ,ψ(Vd(vϕ)−Vs(cψ)−b (tψ−tϕ)+−h (tϕ−tψ)+)xϕψ;\sum_{\phi,\psi}\Big(V^d(v_\phi)-V^s(c_\psi)-b\,(t_\psi-t_\phi)^+-h\,(t_\phi-t_\psi)^+\Big)x_{\phi\psi};ϕ,ψ∑​(Vd(vϕ​)−Vs(cψ​)−b(tψ​−tϕ​)+−h(tϕ​−tψ​)+)xϕψ​;

its optimal value is Jˉ(HT)≥0\bar J(H^T)\ge0Jˉ(HT)≥0.

Formalization targets

Goal: the mechanism-design bound

For every feasible, well-defined mechanism satisfying (ICd), (IRd), (ICs) and (IRs),

E[Π(yT)]≤E[Jˉ(HT)],E\big[\Pi(y^T)\big]\le E\big[\bar J(H^T)\big],E[Π(yT)]≤E[Jˉ(HT)],

with Jˉ(HT)\bar J(H^T)Jˉ(HT) integrable. This is the last display of the proof of Lemma 1 (Online Appendix p. 4), J∗≤Jˉ1≤Jˉ2≤E[Jˉ(HT)]J^*\le\bar J^1\le\bar J^2\le E[\bar J(H^T)]J∗≤Jˉ1≤Jˉ2≤E[Jˉ(HT)].

Milestones

  1. Lemma S.1 (Supplemental Note p. 1): E−ϕ[pϕ]≤vϕE−ϕ[mϕ]−∫v‾vϕE−ϕ[mϕv′] dv′−bE−ϕ[sϕ−tϕ]E_{-\phi}[p_\phi]\le v_\phi E_{-\phi}[m_\phi]-\int_{\underline v}^{v_\phi}E_{-\phi}[m_{\phi_{v'}}]\,dv'-bE_{-\phi}[s_\phi-t_\phi]E−ϕ​[pϕ​]≤vϕ​E−ϕ​[mϕ​]−∫v​vϕ​​E−ϕ​[mϕv′​​]dv′−bE−ϕ​[sϕ​−tϕ​] under (ICd′) and (IRd).
  2. Lemma S.2: E[∑ϕpϕ]≤E[∑ϕVd(vϕ)mϕ−b(sϕ−tϕ)]E[\sum_\phi p_\phi]\le E[\sum_\phi V^d(v_\phi)m_\phi-b(s_\phi-t_\phi)]E[∑ϕ​pϕ​]≤E[∑ϕ​Vd(vϕ​)mϕ​−b(sϕ​−tϕ​)].
  3. Lemma S.4: E[∑ψpψ]≥E[∑ψVs(cψ)mψ+h(sψ−tψ)]E[\sum_\psi p_\psi]\ge E[\sum_\psi V^s(c_\psi)m_\psi+h(s_\psi-t_\psi)]E[∑ψ​pψ​]≥E[∑ψ​Vs(cψ​)mψ​+h(sψ​−tψ​)].
  4. The pathwise step Jˉ2≤E[Jˉ(HT)]\bar J^2\le E[\bar J(H^T)]Jˉ2≤E[Jˉ(HT)]: for every profile and every balanced assignment of outcomes, buyer virtual surplus net of waiting, less seller virtual cost plus waiting, is at most the value of (B).

Lemma S.3 (Supplemental Note p. 3), the seller envelope bound with information rent ∫cψcˉE−ψ[mψc′] dc′\int_{c_\psi}^{\bar c}E_{-\psi}[m_{\psi_{c'}}]\,dc'∫cψ​cˉ​E−ψ​[mψc′​​]dc′, is posed as a further theorem of the mission, the mirror image of Lemma S.1; it is not listed as a milestone.

Here (ICd′) and (ICs′) are the one-dimensional constraints in which the reported arrival time is the true one; Lemmas S.1–S.4 are stated under these weaker hypotheses, as in the paper.

Significance

Lemma 1 of the paper, Jπ,M≤E[Jˉ(HT)]J^{\pi,M}\le E[\bar J(H^T)]Jπ,M≤E[Jˉ(HT)] for every pricing policy π\piπ and matching policy MMM, is the benchmark against which the paper's main theorem measures its waiting-adjusted fixed-price policy, and E[Jˉ(HT)]E[\bar J(H^T)]E[Jˉ(HT)] is in turn bounded by a deterministic fluid value (Lemma 2). Lemma 1 is the composition of two facts: the revelation principle Jπ,M≤J∗J^{\pi,M}\le J^*Jπ,M≤J∗ (Lemma 3) and the mechanism-design bound J∗≤E[Jˉ(HT)]J^*\le E[\bar J(H^T)]J∗≤E[Jˉ(HT)] posed here. The second carries the analytic content: it is where the envelope theorem, the exchange of integrals that turns payments into virtual values, and the timing argument appear.

The result is proved in the paper; its proofs are in the Online Appendix (Section A.2) and the Supplemental Note (Section A). No machine-checked proof of it, or of a dynamic two-sided Myerson bound of this kind, is known. A formal proof would also yield a reusable two-sided, Poisson-arrival version of the interim envelope and revenue-equivalence identities.

Difficulty

Each agent's type is two-dimensional, an arrival time and a valuation or cost, and arrival can only be misreported later. The relaxation to (ICd′)/(ICs′) makes the screening problem one-dimensional, but the interim expectation E−ϕE_{-\phi}E−ϕ​ ranges over a random population whose size is itself random. Turning E[∑ϕ∈HTpϕ]E[\sum_{\phi\in H^T}p_\phi]E[∑ϕ∈HT​pϕ​] into an integral of interim payments requires the Mecke formula for the Poisson process, and the envelope argument must be carried out for interim utilities that are only known to be suprema of affine functions, using no differentiability of the mechanism. The pathwise step is combinatorial: the balancing condition at each instant pairs buyers with sellers discharged together, and those pairs must be compared with an optimal matching of (B), with duplicate types and unmatched agents handled.

A frequent first idea, bounding payments by valuations directly, gives only the first-best surplus E[∑(vϕ−cψ)]E[\sum (v_\phi-c_\psi)]E[∑(vϕ​−cψ​)], which is not the target; the information rents must be extracted.

Formalization scope

The primitives are the published bilateral-trade environment MechanismDesign.BilateralTrade.Environment (densities on intervals, cdfB, cdfS, psiB =Vd=V^d=Vd, psiS =Vs=V^s=Vs, Regular === Assumptions 1–2), and (B) is AssignmentGame.CoreLP.worth of the published assignment game, with sellers as rows and buyers as columns. The conventions committed to are:

  1. Arrivals. Each side is a Poisson(λT\lambda TλT) number of agents with i.i.d. types, arrival time uniform on [0,T][0,T][0,T] and independent of the mark; this is equal in law to the Poisson process restricted to [0,T][0,T][0,T]. Profiles are multisets of types: agents are named by their types, as in the paper, and coincide only with probability zero.
  2. Mechanisms map the reported buyer and seller profiles to an outcome for each type, (s,m,p)(s,m,p)(s,m,p) with t≤s≤Tt\le s\le Tt≤s≤T and m∈{0,1}m\in\{0,1\}m∈{0,1} for listed agents and the balancing condition at every t∈[0,T]t\in[0,T]t∈[0,T]. The causality conditions of the paper (stopping times with respect to the filtration Ht\mathcal H_tHt​) and the timing variables τ,a\tau,aτ,a are dropped. This enlarges the class of mechanisms; it is the clairvoyant class through which the paper's own proof passes (Jˉ1≤Jˉ2\bar J^1\le\bar J^2Jˉ1≤Jˉ2), and the bound for the larger class implies the paper's.
  3. Interim expectations. E−ϕ[g(yϕ^)]E_{-\phi}[g(y_{\hat\phi})]E−ϕ​[g(yϕ^​​)] is the expectation, over an independent copy of HTH^THT, of ggg at the outcome of the report ϕ^\hat\phiϕ^​ inserted among the buyers (the Slivnyak–Mecke reading of "expectation taken with respect to HT∖{ϕ}H^T\setminus\{\phi\}HT∖{ϕ}").
  4. Incentive constraints are quantified over every type and every report in the type spaces, not almost everywhere, with reported arrivals tϕ^∈[tϕ,T]t_{\hat\phi}\in[t_\phi,T]tϕ^​​∈[tϕ​,T].
  5. Added hypotheses. Continuity of fdf^dfd, fsf^sfs on their supports, without which Vd,VsV^d,V^sVd,Vs are defined only almost everywhere; and WellDefined: outcomes are Borel in the listed types, interim payments are integrable, and total absolute payments are integrable. Integrability of Jˉ(HT)\bar J(H^T)Jˉ(HT) and of the virtual-surplus sums is part of each conclusion.

A trivializing formalization is ruled out: the mechanism class is not restricted beyond the paper's feasibility, the constraints are not vacuous (the null mechanism satisfies them, so J∗≥0J^*\ge0J∗≥0), and expectations are of integrable functions. The goal quantifies over mechanisms rather than taking a supremum, which is equivalent because the class is nonempty.

A complete development needs the Mecke formula for the Poisson-count representation, an envelope theorem for suprema of affine functions on an interval, Fubini on the type space, and an exchange argument for matchings. All of these are reusable beyond this mission. Contributions of any of them, of the milestones in any order, or of alternative proofs are welcome.

Selected references

  • Y. Chen and M. Hu, Pricing and Matching with Forward-Looking Buyers and Sellers, Manufacturing & Service Operations Management 22(4):717–734, 2020; submitted manuscript SSRN 2859864. https://ssrn.com/abstract=2859864, https://doi.org/10.1287/msom.2018.0769
  • R. B. Myerson, Optimal Auction Design, Mathematics of Operations Research 6(1):58–73, 1981. https://doi.org/10.1287/moor.6.1.58
  • R. B. Myerson and M. A. Satterthwaite, Efficient Mechanisms for Bilateral Trading, Journal of Economic Theory 29(2):265–281, 1983. https://doi.org/10.1016/0022-0531(83)90048-0
  • P. Milgrom and I. Segal, Envelope Theorems for Arbitrary Choice Sets, Econometrica 70(2):583–601, 2002. https://doi.org/10.1111/1468-0262.00296
  • G. Last and M. Penrose, Lectures on the Poisson Process, Cambridge University Press, 2017 (Mecke equation, Theorem 4.1). https://doi.org/10.1017/9781316104477
9 thms1 active userReviewed
Optimization·Captain: mikedeng1

Pricing and Matching with Forward-Looking Buyers and Sellers 2: Both Fixed Prices Rise with the Buyer Arrival Rate and Fall with the Seller Arrival Rate; Matches and Fluid Profit Rise with BothResearch Paper

Why market size changes prices

Two-sided intermediaries charge buyers and pay sellers while trying to match the two groups. A larger arrival stream on one side changes both the number of possible matches and the prices that clear the market. Chen and Hu study this in a model of buyers and sellers who arrive over time and may wait before requesting a match. Their market-size theorem concerns the deterministic fluid benchmark used to choose fixed prices for the dynamic system. This mission formalizes that theorem from the submitted manuscript, main text §6, rather than the performance of the waiting-adjusted policy itself.

The fluid benchmark matters because it supplies an interpretable price and quantity for each pair of arrival rates. When buyers become more numerous, the model predicts that both the buyer price and the seller payment rise. When sellers become more numerous, both prices fall. In either direction, more arrivals raise the maximum matching quantity and the fluid profit. These claims concern the same valuation and cost distributions as the rates vary, so they isolate the effect of market size.

The two-sided fluid market

Buyer valuations lie in [v‾,vˉ][\underline v,\bar v][v​,vˉ], with density fdf^dfd and distribution function FdF^dFd. Seller costs lie in [c‾,cˉ][\underline c,\bar c][c​,cˉ], with density fsf^sfs and distribution function FsF^sFs. Write Fˉd=1−Fd\bar F^d=1-F^dFˉd=1−Fd. Buyer and seller arrival rates are λd>0\lambda^d>0λd>0 and λs>0\lambda^s>0λs>0, and the horizon has length T>0T>0T>0. A buyer accepts a fixed ask price ppp when the valuation is at least ppp; a seller accepts a fixed bid price www when the cost is at most www. Thus λdFˉd(p)\lambda^d\bar F^d(p)λdFˉd(p) and λsFs(w)\lambda^s F^s(w)λsFs(w) are the corresponding acceptance rates.

The buyer virtual value is Vd(v)=v−Fˉd(v)/fd(v)V^d(v)=v-\bar F^d(v)/f^d(v)Vd(v)=v−Fˉd(v)/fd(v), and the seller virtual cost is Vs(c)=c+Fs(c)/fs(c)V^s(c)=c+F^s(c)/f^s(c)Vs(c)=c+Fs(c)/fs(c). Assumptions 1 and 2 of the manuscript say these functions are weakly increasing on their respective supports. The supports overlap in the sense c‾≤vˉ\underline c\le\bar vc​≤vˉ. The density assumptions ensure that the survival and cost distribution functions have inverse quantiles on [0,1][0,1][0,1]. For a candidate volume μ\muμ, the virtual surplus is

V(μ)=Vd ⁣(Fˉd,−1 ⁣(μλdT))−Vs ⁣(Fs,−1 ⁣(μλsT)).V(\mu)=V^d\!\left(\bar F^{d,-1}\!\left(\frac{\mu}{\lambda^dT}\right)\right)-V^s\!\left(F^{s,-1}\!\left(\frac{\mu}{\lambda^sT}\right)\right).V(μ)=Vd(Fˉd,−1(λdTμ​))−Vs(Fs,−1(λsTμ​)).

The optimal fluid matching quantity μ∗\mu^*μ∗ is the largest μ\muμ in [0,min⁡{λdT,λsT}][0,\min\{\lambda^dT,\lambda^sT\}][0,min{λdT,λsT}] for which V(μ)≥0V(\mu)\ge0V(μ)≥0. The corresponding fixed prices are p∗=Fˉd,−1(μ∗/(λdT))p^*=\bar F^{d,-1}(\mu^*/(\lambda^dT))p∗=Fˉd,−1(μ∗/(λdT)) and w∗=Fs,−1(μ∗/(λsT))w^*=F^{s,-1}(\mu^*/(\lambda^sT))w∗=Fs,−1(μ∗/(λsT)).

The fluid program (D) permits any measurable nonnegative ask and bid price path (π^td,π^ts)(\hat\pi^d_t,\hat\pi^s_t)(π^td​,π^ts​) on [0,T][0,T][0,T] that clears the market at every time: λdFˉd(π^td)=λsFs(π^ts)\lambda^d\bar F^d(\hat\pi^d_t)=\lambda^sF^s(\hat\pi^s_t)λdFˉd(π^td​)=λsFs(π^ts​). Its value Jˉ∗\bar J^*Jˉ∗ is the supremum of buyer revenue minus seller payments, integrated over the horizon. Proposition 1 states that the constant path (p∗,w∗)(p^*,w^*)(p∗,w∗) attains this value and that Jˉ∗=(p∗−w∗)μ∗\bar J^*=(p^*-w^*)\mu^*Jˉ∗=(p∗−w∗)μ∗.

Formalization targets

The goal is Theorem 3, the eight weak comparisons of prices, volume, and fluid profit. For 0<λ1d≤λ2d0<\lambda^d_1\le\lambda^d_20<λ1d​≤λ2d​ with λs\lambda^sλs fixed, it asserts

p1∗≤p2∗,w1∗≤w2∗,μ1∗≤μ2∗,Jˉ1∗≤Jˉ2∗.p^*_1\le p^*_2,\quad w^*_1\le w^*_2,\quad \mu^*_1\le\mu^*_2,\quad \bar J^*_1\le\bar J^*_2.p1∗​≤p2∗​,w1∗​≤w2∗​,μ1∗​≤μ2∗​,Jˉ1∗​≤Jˉ2∗​.

For 0<λ1s≤λ2s0<\lambda^s_1\le\lambda^s_20<λ1s​≤λ2s​ with λd\lambda^dλd fixed, the two price inequalities reverse while the volume and profit inequalities retain their directions. The goal leaves the rates and the two distributions as variables, so it states the market-size result across all systems meeting the paper's assumptions.

The milestones follow the source's dependency structure: Proposition 1 identifies the fluid optimizer; the proof of Lemma S.5 notes that V(μ)V(\mu)V(μ) decreases in μ\muμ; Lemma S.5 establishes the volume and profit comparisons; and the two parts of the proof of Theorem 3 compare the fractions μ∗/(λdT)\mu^*/(\lambda^dT)μ∗/(λdT) and μ∗/(λsT)\mu^*/(\lambda^sT)μ∗/(λsT) served on the growing side. The paper's proofs are in the Supplemental Note, with related proofs in the Online Appendix.

What the result provides

The theorem gives directional predictions without differentiating optimal prices with respect to arrival rates. A platform comparing two otherwise identical markets can infer which way its fluid ask and bid prices move, whether it matches more units, and whether the benchmark profit increases. Proposition 1 also connects the optimization over time-varying paths to a fixed price pair, so the static claims apply to the entire feasible fluid program.

The result is proved in the manuscript; this mission asks for a machine-checked proof of the same statement. The formal development would make the inverse-quantile construction, attainment of the optimal volume, and the nonempty bounded fluid optimization problem explicit. Those objects can be reused in other two-sided pricing or matching models. This mission does not claim that the corresponding stochastic system has identical comparative statics; the manuscript itself focuses on the fluid regime in §6.

Where the proof is difficult

Increasing one arrival rate expands the feasible number of matches, but the two optimal prices do not follow from that fact alone. The buyer price depends on the fraction of buyers served, and the seller price depends on the fraction of sellers served. A growing side can have a larger matched volume yet a smaller matched fraction. The proof must control both quantities under the appropriate inverse distribution function. Profit also needs a comparison of two separate fluid optimization problems, rather than only a comparison of their candidate fixed-price formulas.

Formalization scope

The Lean model reuses the published bilateral-trade environment for the two positive densities, interval supports, cumulative distributions, and virtual values. The new definitions clamp both cumulative distributions to their supports so a price outside the support still has the correct acceptance probability. Inverse quantiles are supported level-set infima, used only at arguments in [0,1][0,1][0,1]. All rates and the horizon are positive; both densities are continuous on their supports; the virtual functions satisfy the paper's regularity assumptions; and c‾≤vˉ\underline c\le\bar vc​≤vˉ. Density continuity is a disclosed addition needed for the asserted inverse-quantile and maximum properties.

The claims involving the value of (D) assume a nonnegative seller cost support. This disclosed domain correction puts the fixed seller price in (D), whose price paths take values in R+2\mathbb R_+^2R+2​; the imported bilateral-trade environment itself permits negative supports. The virtual-surplus and served-fraction inequalities hold without this condition. Every feasible path is measurable, nonnegative on [0,T][0,T][0,T], clears the market at every point of that interval, and has integrable buyer and seller revenue terms. The fluid value is defined as the supremum over those paths. It is never defined by Proposition 1's fixed-price identity, which would erase the optimization claim. The volume set is restricted to [0,min⁡{λdT,λsT}][0,\min\{\lambda^dT,\lambda^sT\}][0,min{λdT,λsT}], so its supremum has the source's intended domain. The source uses “increasing” and “decreasing” in the weak sense.

Useful contributions include lemmas showing the clamped distributions are monotone, their supported inverses have the stated identities, the candidate volume set has a greatest member, the fluid feasible set is nonempty and bounded above, and the two fraction comparisons. These are substantive mathematical obligations: no arbitrary inverse, empty feasible set, non-integrable path, or preassigned fluid value should satisfy the goal by a default value.

Selected references

  • Y. Chen and M. Hu, Pricing and Matching with Forward-Looking Buyers and Sellers, submitted manuscript, SSRN 2859864, 2020; published in Manufacturing & Service Operations Management 22(4), 717–734. SSRN manuscript.
  • T. Börgers, An Introduction to the Theory of Mechanism Design, Oxford University Press, 2015, §3.4.1. Publisher.
8 thms1 active userReviewed
Algorithmic Game TheoryOptimization·Captain: mikedeng1

A game-theoretic approach to computation offloading in mobile cloud computing III: Variational Solutions Are the Nash Equilibria of the Extended Cloudlet-Pricing GameResearch Paper

Why a cloudlet price matters

Mobile devices can divide work between local execution, a nearby shared cloudlet, and a remote cloud. A cloudlet has limited capacity, so one user's choice changes the delay faced by the others. Cardellini and coauthors model this as a generalized Nash equilibrium problem, in which each user's feasible choices depend on the aggregate cloudlet load. The shared utilization cap makes a direct distributed solution awkward: a user must know how much capacity all other users consume before choosing a feasible deviation. The paper therefore introduces an additional player who sets a cloudlet price. This mission formalizes the relationship between the original coupled game and that extended, ordinary Nash game Cardellini et al., accepted manuscript, pp. 10–16.

The offloading game

There are NNN users and nnn cloudlet servers. User uuu chooses fractions xu,mx_{u,m}xu,m​, xu,cletx_{u,\mathrm{clet}}xu,clet​, and xu,cloudx_{u,\mathrm{cloud}}xu,cloud​ for local, cloudlet, and remote-cloud execution. These fractions are nonnegative and sum to one. At most a fraction χ\chiχ can be offloaded, and a user-specific linear power constraint applies. The set K~u\widetilde K_uKu​ contains these private feasible choices. The cloudlet load and the shared feasible set are

L(x)=1n∑u=1Nδuxu,clet,Ω={x:L(x)≤Umax⁡},K=(∏u=1NK~u)∩Ω.L(x)=\frac1n\sum_{u=1}^{N}\delta_u x_{u,\mathrm{clet}},\qquad \Omega=\{x:L(x)\le U_{\max}\},\qquad K=\left(\prod_{u=1}^{N}\widetilde K_u\right)\cap\Omega.L(x)=n1​u=1∑N​δu​xu,clet​,Ω={x:L(x)≤Umax​},K=(u=1∏N​Ku​)∩Ω.

Each user's cost λuRu\lambda_u R_uλu​Ru​ is the response-time expression on p. 10, with coefficients αu,βu,γu,δu\alpha_u,\beta_u,\gamma_u,\delta_uαu​,βu​,γu​,δu​. The VI map FFF is the stack of the user's partial gradients, written in closed form on p. 11. A variational solution xˉ\bar xxˉ belongs to KKK and satisfies ⟨F(xˉ),x−xˉ⟩≥0\langle F(\bar x),x-\bar x\rangle\ge0⟨F(xˉ),x−xˉ⟩≥0 for every x∈Kx\in Kx∈K. The pairing sums all three coordinates for each user. Section 5 assumes 0<Umax⁡<10<U_{\max}<10<Umax​<1 and 0<αu,δu<10<\alpha_u,\delta_u<10<αu​,δu​<1 for every user; the model also has n>0n>0n>0 and 0<χ≤10<\chi\le10<χ≤1 Cardellini et al., pp. 7, 9–11.

The extended game replaces the joint cap in each user's feasible set by a price. User uuu minimizes λuRu(x)+ρ(δu/n)xu,clet\lambda_uR_u(x)+\rho(\delta_u/n)x_{u,\mathrm{clet}}λu​Ru​(x)+ρ(δu​/n)xu,clet​ over K~u\widetilde K_uKu​. A cloudlet manager chooses ρ≥0\rho\ge0ρ≥0 to maximize ρ(L(x)−Umax⁡)\rho(L(x)-U_{\max})ρ(L(x)−Umax​). The manager's choice is part of the Nash equilibrium condition, not a preimposed complementary-slackness equation. The extended game's VI set is Ke=(∏uK~u)×R+K_e=(\prod_u\widetilde K_u)\times\mathbb R_+Ke​=(∏u​Ku​)×R+​ and its map is Fe(x,ρ)=(F(x)+price⁡(ρ),Umax⁡−L(x))F_e(x,\rho)=(F(x)+\operatorname{price}(\rho),U_{\max}-L(x))Fe​(x,ρ)=(F(x)+price(ρ),Umax​−L(x)) Cardellini et al., p. 15.

Formalization targets

The goal is the two-part assertion of Theorem 2. Under the standing assumptions and a strict load bound on the full product of private feasible sets, its equilibrium characterization is

xˉ∈SOL⁡(K,F)⟺∃ρˉ,  (xˉ,ρˉ)∈NE⁡ext.\bar x\in\operatorname{SOL}(K,F) \quad\Longleftrightarrow\quad \exists\bar\rho,\; (\bar x,\bar\rho)\in\operatorname{NE}_{\mathrm{ext}}.xˉ∈SOL(K,F)⟺∃ρˉ​,(xˉ,ρˉ​)∈NEext​.

The same goal states that monotonicity of FFF on ∏uK~u\prod_u\widetilde K_u∏u​Ku​ implies monotonicity of FeF_eFe​ on KeK_eKe​. The three milestones follow the paper's path: the extended Nash game is equivalent to its VI; the shared affine constraint has a scalar multiplier at a VI solution; and the price-load coupling cancels in the monotonicity pairing. The source states the first as a classical result and gives the last two in its proof Cardellini et al., pp. 15–16.

What the result provides

The characterization lets a variational solution of a game with coupled feasible sets be found as an equilibrium of a game whose NNN users have separate feasible sets. The additional scalar has a concrete interpretation: it is the multiplier of the shared cloudlet-utilization cap. The monotonicity assertion connects the paper's analysis of the original VI to methods for the extended one. The paper uses this relationship as the entry point for its distributed algorithms in the next section Cardellini et al., pp. 16–17.

The result is proved in the paper but is not yet machine checked in this mission. A complete development would give reusable Lean interfaces for a finite-dimensional VI on an intersection of a product polytope and one affine half-space, a multiplier characterization for that intersection, and an equilibrium-to-VI conversion for separable player constraints. Those components could support other capacity-sharing games. The theorem drafts remain proof obligations; the definitions compile locally.

Where the proof is delicate

A Nash equilibrium requires each user to be optimal against every privately feasible deviation, including deviations that violate the original shared cap. Consequently, the user's cloudlet cost must remain well defined throughout K~u\widetilde K_uKu​ for the equivalence to hold. A load bound only on KKK does not cover those deviations. The multiplier comparison must also account for the manager's maximization condition and the sign of the load term in FeF_eFe​. For monotonicity, FeF_eFe​ evaluates FFF on the whole product ∏uK~u\prod_u\widetilde K_u∏u​Ku​, a larger domain than KKK. These domain distinctions explain why the statement needs two explicit repairs to the printed theorem Cardellini et al., pp. 15–16.

Formalization scope

Lean represents the users as Fin N and the three processing locations as a finite tier type. Indices therefore start at zero. The costs, the map FFF, the private sets, Ω\OmegaΩ, and KKK are separate definitions. The Euclidean pairing is written as a finite sum; no Pi-type sup norm is used. The extended Nash predicate quantifies over user deviations in K~u\widetilde K_uKu​ and over every nonnegative manager price. The VI predicates are in their own module so the game is not built into their meaning. The theorem's statement remains meaningful at N=0N=0N=0; no positive-user hypothesis is added to Section 5's assumptions.

The first repair requires L(x)<1L(x)<1L(x)<1 at every x∈∏uK~ux\in\prod_u\widetilde K_ux∈∏u​Ku​. The paper's Theorem 1 condition (18) implies this bound, but is stronger than needed. Without the bound, a private deviation can cross the pole of δuxu,clet/(1−L(x))\delta_u x_{u,\mathrm{clet}}/(1-L(x))δu​xu,clet​/(1−L(x)); Lean's total division then assigns a value at a zero denominator even though the paper's cost has none. In a two-user instance with n=1n=1n=1, δu=0.9\delta_u=0.9δu​=0.9 and Umax⁡=0.5U_{\max}=0.5Umax​=0.5, the pole can lie within one user's private feasible interval. The second repair reads “the original game is monotone” on ∏uK~u\prod_u\widetilde K_u∏u​Ku​, because that is precisely the xxx-domain of KeK_eKe​. Monotonicity on KKK alone is insufficient. Both repairs are disclosed in the item notes, and the source wording remains visible in the milestone quotations.

An empty KKK would make a solution statement vacuous; the formal development should retain concrete feasible parameter instances. The map FFF is the paper's displayed formula, rather than a derivative chosen to make a milestone tautological. The VI uses all user coordinates, and the extended game does not silently reinstate the shared constraint in user deviations. Contributions toward the stated multiplier theorem, Nash-to-VI equivalence, and the cancellation identity are welcome, along with reusable finite-dimensional convex-analysis lemmas.

Selected references

  • V. Cardellini, V. De Nitto Personé, V. Di Valerio, F. Facchinei, V. Grassi, F. Lo Presti, and V. Piccialli, A game-theoretic approach to computation offloading in mobile cloud computing, accepted manuscript for Mathematical Programming, 2015, DOI 10.1007/s10107-015-0881-6.
8 thms1 active userReviewed
Mechanism DesignOptimization·Captain: mikedeng1

Double Counting in Supply Chain Carbon Footprinting II: A Carbon Leader Paying Linear Emission-Contingent Rules That Double-Count Internally Earns Its Effort-Contracting ProfitResearch Paper

Motivation

Firms increasingly account for, and offset, the greenhouse-gas emissions of their entire supply chain, not only their own operations. Emissions in a supply chain are often jointly produced: the footprint of one process (transport, packaging, a molten-bulk delivery) depends on the abatement efforts of several firms at once. When each firm is charged for "its" emissions, the question is how to allocate the footprint of a jointly produced process so that every firm has the right incentive to abate.

Caro, Corbett, Tan and Zuidwijk (working paper, 2012; published in Manufacturing & Service Operations Management 15(4), 2013) model this as a moral-hazard problem in teams. Their §4 shows that a social planner who wants the first-best efforts must, under joint production, double-count emissions: the marginal payments across firms for one process exceed the carbon price. This mission formalizes §5, which looks at a different institutional arrangement observed in practice. One firm, the carbon leader (the paper's example is the cosmetics company Natura), voluntarily pays for all supply-chain emissions at a carbon price ppp and compensates the other firms to induce abatement.

Timeline. Holmström (1982) showed that a budget-balanced sharing rule cannot implement efficient efforts in a team. Strausz (1999) showed that sequential moves can restore efficiency. Battaglini (2006) studied joint production with multidimensional output. Caro et al. (2012) placed the carbon leader in this line: a Stackelberg leader whose commitment plays the role of Strausz's sequential structure.

Setting

There are finite sets of firms N\mathcal NN and processes I\mathcal II. Firm nnn chooses efforts en∈[0,A]mn\mathbf e_n\in[0,A]^{m_n}en​∈[0,A]mn​, one for each of its abatement actions, and the effort profile e\mathbf ee lies in the box [0,A]M[0,A]^M[0,A]M. Firm nnn earns profit Vn(en)V_n(\mathbf e_n)Vn​(en​), which is differentiable, concave and componentwise decreasing: effort is costly. Process iii emits fi(e)≥0f_i(\mathbf e)\ge0fi​(e)≥0, which is differentiable, convex and componentwise decreasing on the box. The influence matrix BBB has bn,i=1b_{n,i}=1bn,i​=1 if firm nnn influences process iii, that is, if ∑j∂fi/∂en,j<0\sum_j\partial f_i/\partial e_{n,j}<0∑j​∂fi​/∂en,j​<0, and bn,i=0b_{n,i}=0bn,i​=0 otherwise.

The first best at carbon price ccc maximizes the social value ∑nVn(en)−c∑ifi(e)\sum_n V_n(\mathbf e_n)-c\sum_i f_i(\mathbf e)∑n​Vn​(en​)−c∑i​fi​(e) over the box. At the societal cost c=pSc=p_Sc=pS​ its maximizer is e∗\mathbf e^*e∗, assumed unique.

The carbon leader is a firm NNN, written L in Lean. Every other firm nnn has a reservation profit πˉn\bar\pi_nπˉn​. The signed transfer gng_ngn​ is recorded as a payment by firm nnn; a negative value means compensation from the leader. The leader faces two problems.

  • PEP_EPE​ (6)–(8): it sets gn(en)g_n(\mathbf e_n)gn​(en​), contingent on firm nnn's own effort. It maximizes VN(eN)−p∑ifi(e)+∑n≠Ngn(en)V_N(\mathbf e_N)-p\sum_i f_i(\mathbf e)+\sum_{n\ne N}g_n(\mathbf e_n)VN​(eN​)−p∑i​fi​(e)+∑n=N​gn​(en​) subject to participation, Vn−gn≥πˉnV_n-g_n\ge\bar\pi_nVn​−gn​≥πˉn​, and incentive compatibility, en\mathbf e_nen​ being a best response. Its optimal value is zEz^EzE, attained at efforts eE\mathbf e^EeE.
  • PFP_FPF​ (9)–(11): the same problem, except that the payments gn(f)g_n(\mathbf f)gn​(f) may depend only on the emissions vector f=f(e)\mathbf f=\mathbf f(\mathbf e)f=f(e), because efforts cannot be verified.

Formalization targets

Goal: Proposition 5 (p. 16)

Let (gE,eE)(g^E,\mathbf e^E)(gE,eE) be optimal for PEP_EPE​ and define, for n≠Nn\neq Nn=N,

gn(f)=p∑i∈Ibn,ifi+kn,kn=Vn(enE)−p∑i∈Ibn,ifi(eE)−πˉn.(12)g_n(\mathbf f)=p\sum_{i\in\mathcal I}b_{n,i}f_i+k_n,\qquad k_n=V_n(\mathbf e^E_n)-p\sum_{i\in\mathcal I}b_{n,i}f_i(\mathbf e^E)-\bar\pi_n.\qquad(12)gn​(f)=pi∈I∑​bn,i​fi​+kn​,kn​=Vn​(enE​)−pi∈I∑​bn,i​fi​(eE)−πˉn​.(12)

Then (g,eE)(g,\mathbf e^E)(g,eE) is feasible for PFP_FPF​, its value equals zEz^EzE, and every feasible solution of PFP_FPF​ has value at most zEz^EzE:

max⁡PF  =  zE  =  max⁡PE.\max P_F \;=\; z^E \;=\; \max P_E .maxPF​=zE=maxPE​.

Milestones (proof in App. A.4, p. 26)

  1. Reduction of PEP_EPE​ (p. 15). zE=max⁡e[∑nVn(en)−p∑ifi(e)]−∑n≠Nπˉnz^E=\max_{\mathbf e}\big[\sum_n V_n(\mathbf e_n)-p\sum_i f_i(\mathbf e)\big]-\sum_{n\ne N}\bar\pi_nzE=maxe​[∑n​Vn​(en​)−p∑i​fi​(e)]−∑n=N​πˉn​, attained at eE\mathbf e^EeE.
  2. Participation under (12). Vn(enE)−gn(f(eE))=πˉnV_n(\mathbf e^E_n)-g_n(\mathbf f(\mathbf e^E))=\bar\pi_nVn​(enE​)−gn​(f(eE))=πˉn​ for n≠Nn\ne Nn=N.
  3. First-order identity. ∂Vn(eE)/∂en,j=p∑ibn,i ∂fi(eE)/∂en,j\partial V_n(\mathbf e^E)/\partial e_{n,j}=p\sum_i b_{n,i}\,\partial f_i(\mathbf e^E)/\partial e_{n,j}∂Vn​(eE)/∂en,j​=p∑i​bn,i​∂fi​(eE)/∂en,j​ for n≠Nn\neq Nn=N.
  4. Incentive compatibility under (12). enE\mathbf e^E_nenE​ maximizes Vn(x)−gn(f(e−nE,x))V_n(\mathbf x)-g_n(\mathbf f(\mathbf e^E_{-n},\mathbf x))Vn​(x)−gn​(f(e−nE​,x)) over firm nnn's box.
  5. Upper bound. Every feasible solution of PFP_FPF​ has value at most zEz^EzE.

Companion: Lemma 6 (p. 17)

If 0≤p<pS0\le p<p_S0≤p<pS​, each VnV_nVn​ is supermodular and each fif_ifi​ is submodular, then eE≤e∗\mathbf e^E\le\mathbf e^*eE≤e∗ componentwise.

Significance

The result. Proposition 5 says that, when efforts are unobservable, a carbon leader that can commit loses nothing compared with contracting on efforts. A linear rule, the form footprint accounting already takes in practice, achieves the benchmark. The rule charges each firm the full price ppp on every footprint it influences, so emissions are double-counted internally, yet the supply chain as a whole pays for carbon once. This is how the paper reconciles its §4 impossibility (a planner must double-count) with a footprint-balanced scheme: the leader moves first and is excluded from the incentive constraints. Lemma 6 then says that a leader paying a price below the societal cost induces weakly less effort than the first best.

Formalizing it. Both results are proved in the paper. As far as the platform's catalog shows, neither is machine-checked anywhere. The proof of Proposition 5 combines a reduction of a bilevel problem with a first-order argument for concave games on a box. The appendix compresses these into three sentences, one of which contains a misprint ("since f is concave" for convex). Formalizing it makes the leader–follower conventions explicit: the leader's tie-breaking choice of follower efforts, and the unilateral deviation in (11).

Difficulty

The paper calls PFP_FPF​ a constrained version of PEP_EPE​, but their feasible sets are not literally nested: PFP_FPF​'s payments are functions of the emission vector and PEP_EPE​'s of a firm's own effort. Establishing the common upper bound therefore requires care with the two payment domains and the participation constraint (10).

The incentive compatibility of (12) is the substantive step. It is a global best-response claim on a closed effort box, while the appendix states only a first-order identity at an interior point. The influence indicator bn,i=0b_{n,i}=0bn,i​=0 also has to agree with every action-level derivative, rather than only their sum. For Lemma 6, the optimum at price ppp need not be unique.

Formalization scope

  • Indices. Firms, actions and processes are finite types. Efforts are e : (n : ι) → act n → ℝ, and a unilateral deviation is Function.update.
  • Partial derivatives. A partial derivative is the one-variable deriv along one coordinate. Every theorem that uses one assumes differentiability.
  • Optimality. "Optimal" means feasible with value at least that of every feasible point, never sSup. zEz^EzE is the PEP_EPE​ objective at an optimal pair.
  • Influence matrix. BBB is a real 0/1 matrix, and its defining sign condition is assumed at every profile in the box.
  • Interiority. eE\mathbf e^EeE is assumed interior, the page's "A is a sufficiently large bound to guarantee interior solutions"; the box itself stays closed.
  • Followers' efforts. In PEP_EPE​ and PFP_FPF​ the leader chooses the followers' efforts subject to the best-response constraints, as the page assumes.
  • Hypotheses not imposed. The row and column sums of BBB, the strictness of cnc_ncn​ and gn≤0g_n\le0gn​≤0 are unused and not imposed.

Rule (12) is a function of the emission vector only, and the deviation in (11) moves only firm nnn's efforts. A formalization in which payments may depend on efforts, or in which a deviation also moves the leader, would make the goal trivial and is ruled out. A sorry-free instance with two firms checks that every hypothesis of the goal can be satisfied. Supermodularity reuses the published Supermodularity.Monotonicity.SupermodularOn. All other definitions live in CarbonDoubleCount.Leader.Setting.

Proofs of the milestones, or alternative routes such as a Lagrangian argument for incentive compatibility, are welcome. A general lemma that a concave differentiable function on a box is maximized at an interior stationary point would be reusable beyond this mission.

Selected references

  • F. Caro, C. J. Corbett, T. Tan, R. Zuidwijk, Double-Counting in Supply Chain Carbon Footprinting, working paper dated December 21, 2012; Manufacturing & Service Operations Management 15(4), 2013. https://doi.org/10.1287/msom.2013.0443
  • B. Holmström, Moral Hazard in Teams, Bell Journal of Economics 13(2), 1982. https://doi.org/10.2307/3003457
  • R. Strausz, Efficiency in Sequential Partnerships, Journal of Economic Theory 85, 1999. https://doi.org/10.1006/jeth.1998.2488
  • M. Battaglini, Joint Production in Teams, Journal of Economic Theory 130(1), 2006. https://doi.org/10.1016/j.jet.2005.03.003
  • D. M. Topkis, Minimizing a Submodular Function on a Lattice, Operations Research 26(2), 1978. https://doi.org/10.1287/opre.26.2.305
  • R. K. Sundaram, A First Course in Optimization Theory, Cambridge University Press, 1996. https://doi.org/10.1017/CBO9780511804526
9 thms1 active userReviewed
Algorithmic Game TheoryCombinatorics·Captain: mikedeng1

A Little Charity Guarantees Almost Envy-Freeness I: Under Monotone Valuations an EFX Allocation Exists Whose Unenvied Pool Has Fewer Goods Than There Are Unenvied AgentsResearch Paper

Motivation

Allocating indivisible goods fairly is difficult because a single good can change how an agent compares two bundles. Exact envy-freeness can be impossible (one good, two agents), so the field works with relaxations. Envy-freeness up to any good (EFX) requires each agent to prefer her own bundle to every other agent's bundle after any one good of that other bundle is removed. Whether EFX allocations of all goods always exist is open; the question motivates weakening the requirement that every good be assigned, and asking how little must be left over. Chaudhury, Kavitha, Mehlhorn and Sgouritsa show that leaving a small, unenvied set of goods unallocated (donated to "charity") suffices for every normalized monotone valuation profile (arXiv:1907.04596v3, Theorem 8, pp. 11–12).

Timeline (as recounted on pp. 2–3 of the source):

  • 2004: Lipton, Markakis, Mossel and Saberi introduce the envy-graph and cycle-elimination technique and compute allocations now called EF1 (DOI 10.1145/988772.988792).
  • 2016: Caragiannis, Kurokawa, Moulin, Procaccia, Shah and Wang define EFX (DOI 10.1145/2940716.2940726).
  • 2018: Plaut and Roughgarden prove EFX allocations of all goods exist for two agents, or for any number of agents with identical valuations (arXiv:1707.04769).
  • 2019: Caragiannis, Gravin and Huang introduce EFX with charity: a partial EFX allocation in which each agent gets at least half her value in a Nash-optimal allocation, with no bound on the number or value of donated goods (arXiv:1902.04319).
  • 2020: Chaudhury, Garg and Mehlhorn prove existence of complete EFX allocations for three agents with additive valuations (arXiv:2002.05119).
  • 2020/2021: the present paper bounds the charity: fewer goods than unenvied agents, and no agent envies the donated set. Its SODA 2020 version assumed gross-substitute valuations; the journal version (arXiv v3) needs only normalized monotone valuations (p. 8, §1.4).

The two bounds matter separately. The empty allocation is EFX and donates everything; a cardinality bound alone still lets the donated goods carry most of an agent's value; the value bound alone allows almost all goods to be donated. The paper's later maximin-share and groupwise-maximin-share guarantees (Theorems 14 and 16, treated in the next missions of this series) are derived from exactly the allocation produced here.

Setting

There are n≥1n\ge1n≥1 agents, indexed by N={1,…,n}N=\{1,\ldots,n\}N={1,…,n} in the paper, and a finite set MMM of mmm indivisible goods. Agent iii has a valuation vi:2M→R≥0v_i:2^M\to\mathbb R_{\ge0}vi​:2M→R≥0​. It is normalized if vi(∅)=0v_i(\varnothing)=0vi​(∅)=0 and monotone if S⊆TS\subseteq TS⊆T implies vi(S)≤vi(T)v_i(S)\le v_i(T)vi​(S)≤vi​(T). No additivity, submodularity or substitutes condition is assumed. These are the general valuations of §1.1.1, p. 3, and §2 of the source.

A partial allocation X=(Xi)i∈NX=(X_i)_{i\in N}X=(Xi​)i∈N​ consists of pairwise disjoint bundles Xi⊆MX_i\subseteq MXi​⊆M. Its pool is the goods not assigned to anyone,

P(X)=M∖⋃i∈NXi.P(X)=M\setminus\bigcup_{i\in N}X_i.P(X)=M∖i∈N⋃​Xi​.

The allocation is EFX when

vi(Xj∖{g})≤vi(Xi)for all i,j∈N and every g∈Xj.v_i(X_j\setminus\{g\})\le v_i(X_i) \qquad\text{for all }i,j\in N\text{ and every }g\in X_j.vi​(Xj​∖{g})≤vi​(Xi​)for all i,j∈N and every g∈Xj​.

The paper calls allocations partial by default (p. 3, footnote 5). The envy graph GXG_XGX​ has the agents as vertices and an edge i→ji\to ji→j when vi(Xi)<vi(Xj)v_i(X_i)<v_i(X_j)vi​(Xi​)<vi​(Xj​). A source is a vertex with no incoming edge. Its reachability component C(s)C(s)C(s) includes sss itself. The social welfare of XXX is ϕ(X)=∑ivi(Xi)\phi(X)=\sum_i v_i(X_i)ϕ(X)=∑i​vi​(Xi​). These notions are introduced on pp. 8–9 of the source.

For a set SSS of goods, an envied set is one some agent values more than her own bundle. An inclusion-wise minimal envied subset ZZZ is envied, but no strict subset of ZZZ is envied by anyone. An agent who envies such a ZZZ is a most envious agent of SSS; the paper allows ties. The three update rules U0, U1 and U2 use these sets to change an allocation while maintaining EFX (Definition 3 and Algorithm 2, p. 9).

Formalization targets

The goal is the allocation guarantee of Theorem 8:

∃X partial:X is EFX,vi(P(X))≤vi(Xi) (i∈N),∣P(X)∣<∣sources⁡(GX)∣≤n.\exists X\text{ partial}:\quad X\text{ is EFX},\qquad v_i(P(X))\le v_i(X_i)\ (i\in N),\qquad |P(X)|<|\operatorname{sources}(G_X)|\le n.∃X partial:X is EFX,vi​(P(X))≤vi​(Xi​) (i∈N),∣P(X)∣<∣sources(GX​)∣≤n.

Both the stronger source-count bound and its consequence ∣P(X)∣<n|P(X)|<n∣P(X)∣<n are stated. The theorem asserts existence for every normalized, monotone valuation profile, including nonadditive profiles. It does not require the returned envy graph to be acyclic; acyclicity is an invariant used during construction, rather than an additional advertised property of the final allocation.

The attack path consists of the paper's numbered Lemmas 2, 4, 5, 6 and 7. Lemma 2 concerns removal of envy cycles; its Lean form keeps the page's conditional "if XXX is EFX then X′X'X′ is EFX" and the unconditional strict welfare increase. Lemma 4 records what U0 does and what its failure implies. Lemma 5 records U1's welfare and EFX guarantees. Lemma 6 supplies the finite source, good and agent data required for U2. Lemma 7 states U2's welfare and EFX guarantees. Their source wording and page locations appear in the milestone list; the Lean statements include the exact allocation updates where those lemmas rely on them.

Significance

The theorem guarantees a partial EFX allocation even for general monotone valuations, while limiting the unallocated pool more tightly than the number of agents: its size is below the number of unenvied agents. Every agent also values that entire pool no more than her own bundle. These two bounds distinguish the allocation from a merely EFX partial allocation. They provide the starting point for the paper's later statements on approximate maximin shares under additive valuations (Theorems 14 and 16).

The mathematical result is proved in the paper. The mission asks for a machine-checked proof of its allocation guarantee and of the supporting lemmas, with a reusable Lean account of finite partial allocations, the unallocated pool, envy graphs, minimal envied subsets and the update rules. Such definitions are useful for other fair-division questions involving unallocated indivisible goods.

Difficulty

EFX must hold after the allocation is changed, for every agent pair and for removal of every good from the other agent's new bundle. Simply giving a pool good to an agent may violate that condition, even though it makes the pool smaller. Conversely, replacing a bundle with an envied subset can improve welfare while returning goods to the pool. The source-count bound depends on controlling both kinds of progress. Rule U2 adds a further complication: its bundle moves follow envy paths, so an overlap between selected paths would make the new allocation ambiguous. These are the points at which the straightforward approach of distributing the remaining goods one at a time stops working.

Formalization scope

Lean represents agents by Fin n, goods by Fin m, bundles by finite sets of goods, and valuations by real-valued functions on bundles. The paper's one-based subscripts become zero-based Lean indices. Nonnegative values need no separate hypothesis: normalization and monotonicity imply them. The hypothesis n>0n>0n>0 is explicit because the paper presumes agents and the strict bound ∣P∣<n|P|<n∣P∣<n is false when n=0n=0n=0. The pool is computed as the complement of the union of pairwise disjoint assigned bundles. It is never a free set. EFX ranges over every pair i,ji,ji,j, including i=ji=ji=j, and every good of XjX_jXj​.

The minimal-envied-subset definition checks every strict subset against every agent. The most envious agent is a relation, allowing the paper's arbitrary tie choice. Reachability includes a path of length zero. U2's new allocation is calculated from its paths and chosen subsets. Lemma 6 includes Algorithm 1's acyclic-graph invariant and the ordering property established in its proof. It replaces the printed threshold ∣P∣≥n|P|\ge n∣P∣≥n with ∣P∣≥∣sources⁡(GX)∣|P|\ge|\operatorname{sources}(G_X)|∣P∣≥∣sources(GX​)∣, the exact number of fresh goods its construction needs and the threshold used by the goal. Lemma 7 takes that ordering property to make path assignments unambiguous.

The printed U2 pool expression in Algorithm 2, line 13, is malformed; the formal statement uses the pool of the constructed allocation and the corrected union of leftover goods. Theorem 8's final sentence also asserts an nmV/ΔnmV/\DeltanmV/Δ iteration bound. This mission poses the allocation conclusions, while a future run-level formalization would be needed to state and check the iteration count. The paper's preceding count of U1/U2 steps and successive U0 steps does not by itself imply the printed number. Definitions and the numbered lemma statements are welcome contributions; a purported proof that uses a free empty pool, overlapping bundles, zero agents or an arbitrary replacement allocation would not establish this target.

Selected references

  • B. R. Chaudhury, T. Kavitha, K. Mehlhorn and A. Sgouritsa, A Little Charity Guarantees Almost Envy-Freeness, arXiv:1907.04596v3, 2020; published in SIAM Journal on Computing, 2021. Pinned preprint.
  • R. J. Lipton, E. Markakis, E. Mossel and A. Saberi, On approximately fair allocations of indivisible goods, ACM EC 2004. https://doi.org/10.1145/988772.988792
  • I. Caragiannis, D. Kurokawa, H. Moulin, A. D. Procaccia, N. Shah and J. Wang, The unreasonable fairness of maximum Nash welfare, ACM EC 2016. https://doi.org/10.1145/2940716.2940726
  • B. Plaut and T. Roughgarden, Almost envy-freeness with general valuations, SODA 2018. https://arxiv.org/abs/1707.04769
  • I. Caragiannis, N. Gravin and X. Huang, Envy-freeness up to any item with high Nash welfare: the virtue of donating items, ACM EC 2019. https://arxiv.org/abs/1902.04319
  • B. R. Chaudhury, J. Garg and K. Mehlhorn, EFX exists for three agents, ACM EC 2020. https://arxiv.org/abs/2002.05119
8 thms1 active userReviewed
Partial Differential EquationsProbabilityStochastic Systems·Captain: mikedeng1

Steady-State Analysis of the Join-the-Shortest-Queue Model in the Halfin-Whitt Regime 2: The JSQ Diffusion Limit Satisfies a Foster-Lyapunov Drift ConditionResearch Paper

Motivation

Join-the-shortest-queue (JSQ) is the basic load-balancing rule for a system of nnn parallel single-server queues: every arriving customer joins a queue of minimal length. It is optimal in several senses for homogeneous servers and is the reference policy against which cheaper rules (power-of-ddd choices, join-the-idle-queue) are measured. In the Halfin–Whitt regime, where the arrival rate is nλn\lambdanλ with λ=1−β/n\lambda = 1-\beta/\sqrt nλ=1−β/n​ for a fixed β>0\beta>0β>0, Eschenfeldt and Gamarnik (arXiv:1502.00999) proved that the suitably centred and scaled numbers of queues of length at least one and at least two converge, on finite time intervals, to a two-dimensional reflected diffusion. Braverman (arXiv:1801.05121) justified the interchange of limits (the stationary distributions of the nnn-server chains converge to the stationary distribution of the diffusion). One ingredient of that result is that the diffusion is exponentially ergodic, so it has a unique stationary distribution that is approached at a geometric rate. This mission formalizes the analytic core of that ingredient: a Foster–Lyapunov drift condition for the diffusion's generator.

Timeline. Eschenfeldt and Gamarnik (2015) proved the process-level diffusion limit, the first analysis of JSQ in this regime. Mukherjee, Borst, van Leeuwaarden and Whiting (J. Appl. Probab. 53, 2016) showed that a class of load-balancing schemes shares the same diffusion limit. Braverman (arXiv 2018, v2 2019; Math. Oper. Res. 45(3), 2020) proved tightness of the scaled stationary distributions and exponential ergodicity of the limit, which together give the interchange of limits.

Setting

The state space is the closed quadrant Ω=(−∞,0]×[0,∞)\Omega = (-\infty,0]\times[0,\infty)Ω=(−∞,0]×[0,∞). For f:Ω→Rf:\Omega\to\mathbb Rf:Ω→R write f1=∂f/∂x1f_1 = \partial f/\partial x_1f1​=∂f/∂x1​, f2=∂f/∂x2f_2 = \partial f/\partial x_2f2​=∂f/∂x2​, f11=∂f1/∂x1f_{11} = \partial f_1/\partial x_1f11​=∂f1​/∂x1​. On the boundary ∂Ω\partial\Omega∂Ω these are one-sided derivatives, and C2(Ω)C^2(\Omega)C2(Ω) denotes the twice continuously differentiable functions in this sense.

The diffusion limit Y=(Y1,Y2)Y=(Y_1,Y_2)Y=(Y1​,Y2​) solves

Y1(t)=Y1(0)+2W(t)−βt+∫0t(−Y1(s)+Y2(s)) ds−U(t),Y2(t)=Y2(0)+U(t)−∫0tY2(s) ds,Y_1(t)=Y_1(0)+\sqrt2W(t)-\beta t+\int_0^t(-Y_1(s)+Y_2(s))\,ds-U(t),\qquad Y_2(t)=Y_2(0)+U(t)-\int_0^tY_2(s)\,ds,Y1​(t)=Y1​(0)+2​W(t)−βt+∫0t​(−Y1​(s)+Y2​(s))ds−U(t),Y2​(t)=Y2​(0)+U(t)−∫0t​Y2​(s)ds,

with WWW a standard Brownian motion and UUU a nondecreasing regulator that increases only when Y1=0Y_1=0Y1​=0. For f∈C2(Ω)f\in C^2(\Omega)f∈C2(Ω) with the reflection condition f1(0,x2)=f2(0,x2)f_1(0,x_2)=f_2(0,x_2)f1​(0,x2​)=f2​(0,x2​), its generator acts by

GYf(x)=(−x1+x2−β)f1(x)−x2f2(x)+f11(x),x∈Ω.G_Yf(x)=(-x_1+x_2-\beta)f_1(x)-x_2f_2(x)+f_{11}(x),\qquad x\in\Omega .GY​f(x)=(−x1​+x2​−β)f1​(x)−x2​f2​(x)+f11​(x),x∈Ω.

The construction of the Lyapunov function uses an auxiliary integer n≥1n\ge1n≥1 and the fluid operator

Lf(x)=(−x1+x2−βn)f1(x)−x2f2(x),Lf(x)=\Big(-x_1+x_2-\frac{\beta}{\sqrt n}\Big)f_1(x)-x_2f_2(x),Lf(x)=(−x1​+x2​−n​β​)f1​(x)−x2​f2​(x),

the generator of the deterministic fluid path v1′=−v1+v2−β/nv_1'=-v_1+v_2-\beta/\sqrt nv1′​=−v1​+v2​−β/n​, v2′=−v2v_2'=-v_2v2′​=−v2​. A smoothed indicator ϕ(ℓ,u)\phi^{(\ell,u)}ϕ(ℓ,u) (a piecewise cubic that rises from 000 at ℓ\ellℓ to 111 at uuu) and two levels β<κ1<κ2\beta<\kappa_1<\kappa_2β<κ1​<κ2​ define the PDEs

Lf(1)(x)=−ϕ(κ1/n,κ2/n)(−x1),Lf(2)(x)=−ϕ(κ1/n,κ2/n)(x2),x∈Ω,Lf^{(1)}(x)=-\phi^{(\kappa_1/\sqrt n,\kappa_2/\sqrt n)}(-x_1),\qquad Lf^{(2)}(x)=-\phi^{(\kappa_1/\sqrt n,\kappa_2/\sqrt n)}(x_2),\qquad x\in\Omega,Lf(1)(x)=−ϕ(κ1​/n​,κ2​/n​)(−x1​),Lf(2)(x)=−ϕ(κ1​/n​,κ2​/n​)(x2​),x∈Ω,

each with the reflection condition. Their solutions are written explicitly in terms of the fluid hitting times τ~(κ)\tilde\tau^{(\kappa)}τ~(κ) (through the Lambert W function) and τ\tauτ, and the curves Γ(κ)\Gamma^{(\kappa)}Γ(κ) that split Ω\OmegaΩ into regions S0,…,S3S_0,\dots,S_3S0​,…,S3​.

Formalization targets

Goal: Theorem 4

For every β>0\beta>0β>0 there exist c>0c>0c>0, d>0d>0d>0, a compact set KKK and V∈C2(Ω)V\in C^2(\Omega)V∈C2(Ω) with V≥1V\ge1V≥1 on Ω\OmegaΩ, V(x)→∞V(x)\to\inftyV(x)→∞ as ∣x∣→∞|x|\to\infty∣x∣→∞ in Ω\OmegaΩ, and V1(0,x2)=V2(0,x2)V_1(0,x_2)=V_2(0,x_2)V1​(0,x2​)=V2​(0,x2​), such that

GYV(x)≤−cV(x)+d 1(x∈K),x∈Ω.G_YV(x)\le -cV(x)+d\,1(x\in K),\qquad x\in\Omega .GY​V(x)≤−cV(x)+d1(x∈K),x∈Ω.

The constants are unspecified and depend only on β\betaβ; the statement involves no nnn.

Milestones

  1. (5.6)–(5.7): ϕ(ℓ,u)\phi^{(\ell,u)}ϕ(ℓ,u) has an absolutely continuous derivative with ϕ′(ℓ)=ϕ′(u)=0\phi'(\ell)=\phi'(u)=0ϕ′(ℓ)=ϕ′(u)=0, ∣ϕ′∣≤4/(u−ℓ)|\phi'|\le 4/(u-\ell)∣ϕ′∣≤4/(u−ℓ) and ∣ϕ′′∣≤12/(u−ℓ)2|\phi''|\le12/(u-\ell)^2∣ϕ′′∣≤12/(u−ℓ)2.
  2. Lemma 10: τ~(κ)\tilde\tau^{(\kappa)}τ~(κ) extends to {x2=0}\{x_2=0\}{x2​=0}, vanishes on {x1=−κ/n}\{x_1=-\kappa/\sqrt n\}{x1​=−κ/n​}, and has explicit one-sided partial derivatives.
  3. Lemma 11: the explicit f(1)f^{(1)}f(1) is nonnegative, lies in C2(Ω)C^2(\Omega)C2(Ω), and has explicit first partials.
  4. Lemma 12: the explicit six-case f(2)f^{(2)}f(2) is well defined, lies in C2(Ω)C^2(\Omega)C2(Ω), and has the partials (C.8)–(C.9).
  5. Lemma 8: for ϵ>0\epsilon>0ϵ>0, κ1=β+ϵ\kappa_1=\beta+\epsilonκ1​=β+ϵ, κ2=β+2ϵ\kappa_2=\beta+2\epsilonκ2​=β+2ϵ, both PDEs have C2(Ω)C^2(\Omega)C2(Ω) solutions with
f(1)≤log⁡2,f(2)≤log⁡2+ϵβ on a box,∣f1(1)∣≤4nϵlog⁡2, ∣f11(1)∣≤12nϵ2log⁡2, ∣f1(2)∣≤nβ, ∣f11(2)∣≤nβϵ(1+4β+ϵϵ).f^{(1)}\le\log2,\quad f^{(2)}\le\log2+\tfrac{\epsilon}{\beta}\ \text{on a box},\qquad |f^{(1)}_1|\le\tfrac{4\sqrt n}{\epsilon}\log2,\ |f^{(1)}_{11}|\le\tfrac{12n}{\epsilon^2}\log2,\ |f^{(2)}_1|\le\tfrac{\sqrt n}{\beta},\ |f^{(2)}_{11}|\le\tfrac{n}{\beta\epsilon}\big(1+4\tfrac{\beta+\epsilon}{\epsilon}\big).f(1)≤log2,f(2)≤log2+βϵ​ on a box,∣f1(1)​∣≤ϵ4n​​log2, ∣f11(1)​∣≤ϵ212n​log2, ∣f1(2)​∣≤βn​​, ∣f11(2)​∣≤βϵn​(1+4ϵβ+ϵ​).

Significance

By the criterion of Down, Meyn and Tweedie (Ann. Probab. 23, 1995, Theorem 5.2), a drift inequality of the form (5.2) with compact KKK and V→∞V\to\inftyV→∞ implies that the diffusion is positive Harris recurrent and exponentially ergodic in the VVV-norm (Theorem 3 and Corollary 1 of the paper). Together with tightness of the prelimit stationary distributions (mission 1 of this series), this gives convergence of the stationary distributions of the JSQ chains to that of the diffusion. That justifies approximating steady-state JSQ performance by the diffusion.

The construction itself is reusable. Solving the fluid-model PDE Lf=−(smoothed indicator)Lf=-(\text{smoothed indicator})Lf=−(smoothed indicator) and exponentiating the solution is a general route from fluid stability to exponential ergodicity of a diffusion. This mission writes that route out for one reflected diffusion with explicit hitting-time formulas.

Status: the result is proved in the paper; no part of it has a machine-checked proof. Theorem 3 and Corollary 1 are not posed, because they rest on the cited Down–Meyn–Tweedie theory and on the existence of the reflected process, neither of which Mathlib has.

Difficulty

The reflecting boundary rules out the first candidates. The condition V1(0,x2)=V2(0,x2)V_1(0,x_2)=V_2(0,x_2)V1​(0,x2​)=V2​(0,x2​) excludes exponentials of linear functions such as ea(x2−x1)e^{a(x_2-x_1)}ea(x2​−x1​) unless a=0a=0a=0. The natural candidate from the fluid model, the exponential of the fluid time needed to reach a box, has a discontinuous integrand and is not in C2(Ω)C^2(\Omega)C2(Ω), so GYG_YGY​ cannot be applied to it. The milestone functions f(1)f^{(1)}f(1), f(2)f^{(2)}f(2) are given piecewise across the curves Γ(κ1)\Gamma^{(\kappa_1)}Γ(κ1​), Γ(κ2)\Gamma^{(\kappa_2)}Γ(κ2​) and the lines x1=−κi/nx_1=-\kappa_i/\sqrt nx1​=−κi​/n​, through implicitly defined hitting times. Three things are hard: C2C^2C2 regularity across these interfaces, the Lambert W asymptotics as x2↓0x_2\downarrow0x2​↓0, and second-derivative bounds that hold uniformly on the unbounded domain Ω\OmegaΩ with the stated dependence on nnn and ϵ\epsilonϵ.

Formalization scope

  • Points of Ω\OmegaΩ are pairs ℝ × ℝ, and Ω\OmegaΩ is Set.Iic 0 ×ˢ Set.Ici 0. Partials are fderivWithin ℝ f Ω in the coordinate directions, which are one-sided at boundary points. C2(Ω)C^2(\Omega)C2(Ω) is ContDiffOn ℝ 2 f Ω.
  • GYG_YGY​ is the generator formula of p. 16, and every statement applies it only to C2(Ω)C^2(\Omega)C2(Ω) functions with the reflection condition. The process (2.1) and the extended generator are not constructed. V∈C2(Ω)V\in C^2(\Omega)V∈C2(Ω) is part of the conclusion of Theorem 4, so a VVV whose derivatives within Ω\OmegaΩ are undefined (and default to 000) cannot satisfy (5.2) vacuously. Dropping V≥1V\ge1V≥1, V→∞V\to\inftyV→∞, the compactness of KKK, or c,d>0c,d>0c,d>0 would make (5.2) trivial, and all four are kept.
  • The auxiliary nnn is an integer n≥1n\ge1n≥1; no relation between β\betaβ and n\sqrt nn​ is assumed.
  • The Lambert W function (principal branch) is defined locally, since Mathlib has none. τ\tauτ takes values in WithTop ℝ, with ∞\infty∞ when no solution exists, and its real value is used only where it is finite. ν∗\nu^*ν∗ (the curve Γ(κ)\Gamma^{(\kappa)}Γ(κ)) is the solution of a nonlinear system whose uniqueness is Lemma 5.
  • Two printed misprints are corrected. Lemma 8's "for every ϵ>0\epsilon>0ϵ>0" is formalized with κ1=β+ϵ\kappa_1=\beta+\epsilonκ1​=β+ϵ, κ2=β+2ϵ\kappa_2=\beta+2\epsilonκ2​=β+2ϵ, the choice its proof makes; read literally, the printed statement is false. In (C.5) of Lemma 10 the identity τ~2=−τ~1τ~\tilde\tau_2=-\tilde\tau_1\tilde\tauτ~2​=−τ~1​τ~ is formalized as τ~2=τ~1τ~\tilde\tau_2=\tilde\tau_1\tilde\tauτ~2​=τ~1​τ~, which is what implicit differentiation gives.
  • Lemmas 5, 6 and 9 (properties of Γ(κ)\Gamma^{(\kappa)}Γ(κ), τ\tauτ and WWW) are milestones of mission 1; here their objects are only definitions. Contributions proving those facts in this namespace, or general lemmas on the Lambert W function, are welcome.

Selected references

  • A. Braverman, Steady-State Analysis of the Join-the-Shortest-Queue Model in the Halfin-Whitt Regime, Math. Oper. Res. 45(3), 2020; preprint arXiv:1801.05121v2, 2019. https://arxiv.org/abs/1801.05121, https://doi.org/10.1287/moor.2019.1023
  • P. Eschenfeldt, D. Gamarnik, Join the Shortest Queue with Many Servers. The Heavy-Traffic Asymptotics, 2015 (Math. Oper. Res. 43(3), 2018). https://arxiv.org/abs/1502.00999
  • D. Down, S. P. Meyn, R. L. Tweedie, Exponential and Uniform Ergodicity of Markov Processes, Ann. Probab. 23, 1671–1691, 1995. https://doi.org/10.1214/aop/1176987798
  • D. Mukherjee, S. C. Borst, J. S. H. van Leeuwaarden, P. A. Whiting, Universality of Load Balancing Schemes on the Diffusion Scale, J. Appl. Probab. 53, 1111–1124, 2016. https://projecteuclid.org/euclid.jap/1481132840
9 thms1 active userReviewed
Control TheoryMachine LearningProbability+1·Captain: mikedeng1

Convergence Analysis of Machine Learning for Mean Field Control, Finite Horizon: The N-Agent Discrete-Time Neural-Net Value Is Within O(N^{-1/max(d,4)} + n_in^{-1/(3(d+1))} + √Δt) of the MKV OptimumResearch Paper

Motivation

Mean field control studies the optimal control of a very large population of interacting agents through the limit in which the population is replaced by its distribution. The limit problem is a McKean–Vlasov (MKV) control problem: the dynamics and the cost of a representative agent depend on the law of its own state. Such problems arise in systemic risk, crowd motion, energy management and the planning of large distributed systems; the reference account is Carmona and Delarue, Probabilistic Theory of Mean Field Games with Applications I–II (Springer, 2018).

MKV control problems are rarely solvable in closed form, and grid-based numerical methods suffer from the dimension of the state. A practical alternative is to simulate a finite population on a time grid, parametrize the feedback control by a neural network, and minimize the simulated cost by stochastic gradient descent. Carmona and Laurière, Convergence analysis of machine learning algorithms for the numerical solution of mean field control and games: II — the finite horizon case (Ann. Appl. Probab. 32(6), 2022, doi:10.1214/21-AAP1715), quantify how close this computable proxy is to the true optimum. This mission formalizes the statement of their main result and of the propositions and lemmas its proof is built from.

Setting

Fix a horizon T>0T>0T>0, a state dimension ddd, a control dimension kkk (controls take values in A=RkA=\mathbb R^kA=Rk), an initial law μ0\mu_0μ0​ on Rd\mathbb R^dRd, a drift b(t,x,μ,α)b(t,x,\mu,\alpha)b(t,x,μ,α), a volatility σ(t,x,μ)\sigma(t,x,\mu)σ(t,x,μ) (not controlled), a running cost f(t,x,μ,α)f(t,x,\mu,\alpha)f(t,x,μ,α) and a terminal cost g(x,μ)g(x,\mu)g(x,μ); μ\muμ ranges over P2(Rd)\mathcal P_2(\mathbb R^d)P2​(Rd), the probability measures with finite second moment, compared with the Wasserstein distance W2W_2W2​.

Problem 1 (the MKV problem). On a probability space carrying a ddd-dimensional Wiener process WWW and an independent X0∼μ0X_0\sim\mu_0X0​∼μ0​, minimize over progressively measurable square-integrable controls α∈A\alpha\in\mathbb Aα∈A

J(α)=E[∫0Tf(t,Xt,L(Xt),αt) dt+g(XT,L(XT))],dXt=b(t,Xt,L(Xt),αt) dt+σ(t,Xt,L(Xt)) dWt,J(\alpha)=\mathbb E\Big[\int_0^Tf(t,X_t,\mathcal L(X_t),\alpha_t)\,dt+g(X_T,\mathcal L(X_T))\Big],\qquad dX_t=b(t,X_t,\mathcal L(X_t),\alpha_t)\,dt+\sigma(t,X_t,\mathcal L(X_t))\,dW_t,J(α)=E[∫0T​f(t,Xt​,L(Xt​),αt​)dt+g(XT​,L(XT​))],dXt​=b(t,Xt​,L(Xt​),αt​)dt+σ(t,Xt​,L(Xt​))dWt​,

where L(Xt)\mathcal L(X_t)L(Xt​) is the law of XtX_tXt​.

Problem 3 (the NNN-agent problem). NNN agents with independent Wiener processes WiW^iWi and i.i.d. initial states X0i∼μ0X^i_0\sim\mu_0X0i​∼μ0​ use a common feedback function v(t,x)v(t,x)v(t,x); the law is replaced by the empirical measure μtN=1N∑jδXtj\mu^N_t=\frac1N\sum_j\delta_{X^j_t}μtN​=N1​∑j​δXtj​​, and the cost is JN(v)=1N∑iE[∫0Tf(t,Xti,μtN,v(t,Xti)) dt+g(XTi,μTN)]J^N(v)=\frac1N\sum_i\mathbb E[\int_0^Tf(t,X^i_t,\mu^N_t,v(t,X^i_t))\,dt+g(X^i_T,\mu^N_T)]JN(v)=N1​∑i​E[∫0T​f(t,Xti​,μtN​,v(t,Xti​))dt+g(XTi​,μTN​)].

Problem 2 (the computable proxy). With Δt=T/NT\Delta t=T/N_TΔt=T/NT​ and tn=nΔtt_n=n\Delta ttn​=nΔt, the NNN agents follow the Euler scheme

Xˇtn+1i=Xˇtni+b(tn,Xˇtni,μˇtn,φ(tn,Xˇtni))Δt+σ(tn,Xˇtni,μˇtn)ΔWˇni\check X^i_{t_{n+1}}=\check X^i_{t_n}+b(t_n,\check X^i_{t_n},\check\mu_{t_n},\varphi(t_n,\check X^i_{t_n}))\Delta t+\sigma(t_n,\check X^i_{t_n},\check\mu_{t_n})\Delta\check W^i_nXˇtn+1​i​=Xˇtn​i​+b(tn​,Xˇtn​i​,μˇ​tn​​,φ(tn​,Xˇtn​i​))Δt+σ(tn​,Xˇtn​i​,μˇ​tn​​)ΔWˇni​

with i.i.d. Gaussian increments ΔWˇni∼N(0,Δt Id)\Delta\check W^i_n\sim\mathcal N(0,\Delta t\,I_d)ΔWˇni​∼N(0,ΔtId​), and JˇN(φ)\check J^N(\varphi)JˇN(φ) is the expected time-discretized cost. The feedback φ\varphiφ ranges over Nd+1,nin,kψ\mathbf N^\psi_{d+1,n_{\rm in},k}Nd+1,nin​,kψ​, the one-hidden-layer neural networks from (t,x)∈Rd+1(t,x)\in\mathbb R^{d+1}(t,x)∈Rd+1 to Rk\mathbb R^kRk with ninn_{\rm in}nin​ hidden neurons and activation ψ\psiψ; ψ\psiψ is 2π2\pi2π-periodic, of class C3\mathcal C^3C3, with ψ^1=∫−ππψ(x)e−ix dx≠0\hat\psi_1=\int_{-\pi}^{\pi}\psi(x)e^{-ix}\,dx\ne0ψ^​1​=∫−ππ​ψ(x)e−ixdx=0.

The standing assumptions (§2.2 and Appendix A of the paper) make b,σb,\sigmab,σ affine in (x,μˉ,α)(x,\bar\mu,\alpha)(x,μˉ​,α) (A1), make f,gf,gf,g differentiable in (x,α,μ)(x,\alpha,\mu)(x,α,μ) (in μ\muμ in the sense of the L-derivative) with Lipschitz derivatives (A2)–(A3), make fff strongly convex in α\alphaα (A4), and add the regularity (B1)–(B3), (C1)–(C3). Under them the reduced Hamiltonian b⋅y+fb\cdot y+fb⋅y+f has a minimizer α^(t,x,μ,y)\hat\alpha(t,x,\mu,y)α^(t,x,μ,y), the Pontryagin forward–backward system has a solution (X,Y,Z)(X,Y,Z)(X,Y,Z) with marginal flow μt\mu_tμt​, and Yt=V(t,Xt)Y_t=V(t,X_t)Yt​=V(t,Xt​) for a regular decoupling field VVV. The optimal feedback is v^(t,x)=α^(t,x,μt,V(t,x))\hat v(t,x)=\hat\alpha(t,x,\mu_t,V(t,x))v^(t,x)=α^(t,x,μt​,V(t,x)).

Formalization targets

Goal: Theorem 3 (p. 4070)

inf⁡α∈AJ(α) ≥ inf⁡φ∈Nd+1,nin,kψJˇN(φ)−C(N−1/max⁡(d,4)1+ln⁡(N)1{d=4}+nin−1/(3(d+1))+Δt)\inf_{\alpha\in\mathbb A}J(\alpha)\ \ge\ \inf_{\varphi\in\mathbf N^\psi_{d+1,n_{\rm in},k}}\check J^N(\varphi)-C\Big(N^{-1/\max(d,4)}\sqrt{1+\ln(N)\mathbf 1_{\{d=4\}}}+n_{\rm in}^{-1/(3(d+1))}+\sqrt{\Delta t}\Big)α∈Ainf​J(α) ≥ φ∈Nd+1,nin​,kψ​inf​JˇN(φ)−C(N−1/max(d,4)1+ln(N)1{d=4}​​+nin−1/(3(d+1))​+Δt​)

for a constant CCC depending only on the data and on ψ\psiψ, for all N,nin,NT≥1N,n_{\rm in},N_T\ge1N,nin​,NT​≥1.

The three steps

  1. Proposition 7 (finite population): inf⁡αJ(α)≥JN(v^)−c1N−2/max⁡(d,4)(1+ln⁡(N)1{d=4})\inf_\alpha J(\alpha)\ge J^N(\hat v)-c_1\sqrt{N^{-2/\max(d,4)}(1+\ln(N)\mathbf 1_{\{d=4\}})}infα​J(α)≥JN(v^)−c1​N−2/max(d,4)(1+ln(N)1{d=4}​)​.
  2. Proposition 8 (network class): some φ^∈Nd+1,nin,kψ\hat\varphi\in\mathbf N^\psi_{d+1,n_{\rm in},k}φ^​∈Nd+1,nin​,kψ​ with bounded Lipschitz constants has JN(v^)≥JN(φ^)−K2nin−1/(3(d+1))J^N(\hat v)\ge J^N(\hat\varphi)-K_2n_{\rm in}^{-1/(3(d+1))}JN(v^)≥JN(φ^​)−K2​nin−1/(3(d+1))​.
  3. Proposition 15 (time step): ∣JN(φ)−JˇN(φ)∣≤CΔt|J^N(\varphi)-\check J^N(\varphi)|\le C\sqrt{\Delta t}∣JN(φ)−JˇN(φ)∣≤CΔt​ for regular feedbacks.

Supporting results

Lemma 18 (Lipschitz continuity of α^\hat\alphaα^), Proposition 10 (network approximation of a function and its first two xxx-derivatives at rate nin−1/(2(d+1))n_{\rm in}^{-1/(2(d+1))}nin−1/(2(d+1))​ on [0,T]×Bˉd(0,R)[0,T]\times\bar B_d(0,R)[0,T]×Bˉd​(0,R)), (C.2) (moment bounds uniform in NNN), Proposition 13 (stability of JNJ^NJN under feedbacks close on a ball), Lemma 19 (time regularity of the particle system) and Lemma 14 (strong error CNΔtCN\Delta tCNΔt of the Euler scheme).

Significance

The theorem certifies the standard deep-learning approach to mean field control: the value computed by minimizing the simulated finite-population, discrete-time cost over networks is, up to an explicit error with separate rates for population size, network width and time step, no larger than the true MKV value. The rates also expose where the method loses: the nin−1/(3(d+1))n_{\rm in}^{-1/(3(d+1))}nin−1/(3(d+1))​ term is a curse of dimensionality in the network step (Remark 11).

The result is proved in the paper (§3 and Appendices B–D); the proof of Proposition 7 is a modification of Carmona–Delarue, Vol. II, Theorem 6.17. None of it is formalized. A formal development would provide the first machine-checked bridge between a stochastic control problem of McKean–Vlasov type and a concrete learning architecture, and its intermediate results (moment and stability estimates for interacting particle systems, the strong error of an Euler scheme with constants independent of the number of particles, simultaneous approximation of a function and its derivatives by periodic networks) are reusable on their own.

Difficulty

Each step rests on research-level analysis. Proposition 7 needs the Pontryagin principle for McKean–Vlasov control, the propagation of chaos for the optimally controlled system, and the Fournier–Guillin rate for empirical measures, whose max⁡(d,4)\max(d,4)max(d,4) exponent and logarithmic correction at d=4d=4d=4 appear in the bound. Proposition 8 requires approximating the optimal feedback together with its first two derivatives in xxx, because the time-discretization step needs regularity of the approximating network; standard universal approximation results give no control of derivatives, and the available simultaneous approximation results hold only for periodic functions, which forces a localization to a ball whose radius must be balanced against the approximation rate. Proposition 15 cannot be obtained from classical Euler estimates for a system of NdNdNd equations, whose constants grow with NNN: the estimates must be done agent by agent through the empirical measure.

Formalization scope

Rd\mathbb R^dRd is Fin d → ℝ with the sup norm; every constant is existential, so all statements are norm-independent. Time is ℝ≥0 for processes. Solutions of (2.2) and (3.2) use the published Itô-process predicate Peng1990.SMP.IsItoProcess (for X−X0X-X_0X−X0​, with continuous paths and the natural filtration of the initial data and the noise); the backward equation of the Pontryagin system uses Peng1990.SMP.SolvesBSDE; W2W_2W2​ is the published WassersteinDRO.Duality.wassersteinDistance. Problem 2 is law-level: the Euler scheme is a deterministic map of initial positions and increments, integrated against μ0⊗N⊗N(0,ΔtId)⊗NNT\mu_0^{\otimes N}\otimes\mathcal N(0,\Delta t I_d)^{\otimes NN_T}μ0⊗N​⊗N(0,ΔtId​)⊗NNT​; Lemma 14 couples it with the continuous system through the Brownian increments, as the paper does.

The hypotheses are bundled as Coeff ((A1)–(A4), (B1), (B3), (C1), with witnesses for every derivative; the L-derivative is the Fréchet derivative of the lift on L2([0,1])L^2([0,1])L2([0,1])) and DecouplingField ((2.6), existence of a solution of (2.7) on some Problem-1 space, Yt=V(t,Xt)Y_t=V(t,X_t)Yt​=V(t,Xt​), (B2), (C2), (C3)). Disclosed readings: ggg's L-convexity is (A4) without α\alphaα and with right-hand side 000; unique solvability of (2.7), uniqueness and Lipschitz continuity of α^\hat\alphaα^ are not assumed (the last is Lemma 18); Proposition 13 lets its constant depend on bounds for ∣v(0,0)∣,∣w(0,0)∣|v(0,0)|,|w(0,0)|∣v(0,0)∣,∣w(0,0)∣, as its proof does; Proposition 15 reads the first bound of (3.16) as C1(1+∣(t,x)∣)C_1(1+|(t,x)|)C1​(1+∣(t,x)∣). Infima are stated pointwise (∀α ∀η>0 ∃φ\forall\alpha\,\forall\eta>0\,\exists\varphi∀α∀η>0∃φ), never as a real ⨅, and N,nin,NT≥1N,n_{\rm in},N_T\ge1N,nin​,NT​≥1 throughout. The hypothesis bundle is not vacuous: a sorry-free check exhibits the model b=σ=0b=\sigma=0b=σ=0, f=∣α∣2f=|\alpha|^2f=∣α∣2, g=0g=0g=0, μ0=δ0\mu_0=\delta_0μ0​=δ0​ with α^=0\hat\alpha=0α^=0, V=0V=0V=0 satisfying every assumption on any Problem-1 space. Junk values cannot trivialize the goal: the costs are Bochner integrals whose integrability follows from the assumptions, so a non-integrable cost would only make the claimed inequality harder, never vacuous.

Contributions are welcome at every level: proofs of the milestones, in particular the self-contained Lemmas 14 and 19, (C.2) and Proposition 13; infrastructure for the L-derivative, Itô's formula for the published Itô processes, empirical-measure convergence rates, and trigonometric approximation of Cr\mathcal C^rCr periodic functions.

Selected references

  • R. Carmona, M. Laurière, Convergence analysis of machine learning algorithms for the numerical solution of mean field control and games: II — the finite horizon case, Ann. Appl. Probab. 32(6) (2022), 4065–4105. https://doi.org/10.1214/21-AAP1715
  • R. Carmona, F. Delarue, Probabilistic Theory of Mean Field Games with Applications I–II, Springer (2018). https://doi.org/10.1007/978-3-319-58920-6
  • N. Fournier, A. Guillin, On the rate of convergence in Wasserstein distance of the empirical measure, Probab. Theory Relat. Fields 162 (2015), 707–738. https://doi.org/10.1007/s00440-014-0583-7
  • S. Peng, A general stochastic maximum principle for optimal control problems, SIAM J. Control Optim. 28(4) (1990), 966–979. https://doi.org/10.1137/0328054
15 thms1 active userReviewed
OptimizationProbability·Captain: mikedeng1

Retail Assortment Planning in the Presence of Consumer Search IV: Under Overlapping-Assortment Search, the Full Assortment Is Optimal Once Search Is Cheap EnoughResearch Paper

Assortment planning when consumers can shop elsewhere

A retailer that decides which variants of a product to stock (its assortment) trades off the revenue each variant brings against the operational cost of carrying it. The classical model of this trade-off uses the multinomial logit (MNL) model of consumer choice: a consumer who does not like what the store offers simply buys nothing (van Ryzin & Mahajan 1999). In practice a consumer who does not like what the store offers can also go elsewhere. Cachon, Terwiesch and Xu (working paper, December 2002; published in MSOM 7(4), 2005) add this consumer search to the MNL model and ask how it changes the retailer's optimal assortment.

They study two search models. In the overlapping assortment model (§3.1), a consumer who searches pays a search cost bbb and then finds a competitor carrying every variant, so search can only reveal variants the store does not carry. This mission formalizes the paper's Theorem 6 about that model: when search is cheap, the store's best response is to carry everything.

Setting

There are nnn variants, N={1,…,n}N = \{1,\dots,n\}N={1,…,n}, and an assortment is a subset S⊆NS \subseteq NS⊆N; write Sˉ=N−S\bar S = N - SSˉ=N−S for the variants the store does not carry. A consumer's utility for variant iii is Ui=(ui−pi)+ζiU_i = (u_i - p_i) + \zeta_iUi​=(ui​−pi​)+ζi​, where ui−piu_i - p_iui​−pi​ is the variant's expected net utility and the shocks ζi\zeta_iζi​ are independent with the zero-mean Gumbel distribution function F(x)=exp⁡[−exp⁡(−(x/μ+γ))]F(x) = \exp[-\exp(-(x/\mu + \gamma))]F(x)=exp[−exp(−(x/μ+γ))], with scale μ>0\mu > 0μ>0 and γ\gammaγ Euler's constant. The no-purchase option is a "faux variant" 000 with utility U0U_0U0​. Its preference is vi=exp⁡((ui−pi)/μ)>0v_i = \exp((u_i - p_i)/\mu) > 0vi​=exp((ui​−pi​)/μ)>0, and v0>0v_0 > 0v0​>0 is the preference of not buying.

Without search, variant i∈Si \in Si∈S has the MNL demand (2), qim(S)=vi/(∑k∈Svk+v0)q_i^m(S) = v_i / (\sum_{k\in S} v_k + v_0)qim​(S)=vi​/(∑k∈S​vk​+v0​). Following Theorem 1, put λ(u)=exp⁡[−(u/μ+γ)]\lambda(u) = \exp[-(u/\mu + \gamma)]λ(u)=exp[−(u/μ+γ)] and

H(u,S)=exp⁡(−λ(u)(v0+∑j∈Svj)),H(u, S) = \exp\Big(-\lambda(u)\big(v_0 + \textstyle\sum_{j\in S} v_j\big)\Big),H(u,S)=exp(−λ(u)(v0​+∑j∈S​vj​)),

the probability that every option the store offers, including not buying, has utility at most uuu. The best utility outside the store, max⁡i∈SˉUi\max_{i\in\bar S} U_imaxi∈Sˉ​Ui​, has a density w(yˉ,S)w(\bar y, S)w(yˉ​,S). A consumer whose best in-store utility is yyy searches when the expected gain exceeds the cost,

ΦS(y)=∫y∞(yˉ−y) w(yˉ,S) dyˉ ≥ b,(5)\Phi_S(y) = \int_y^\infty (\bar y - y)\, w(\bar y, S)\, d\bar y \ \ge\ b, \qquad (5)ΦS​(y)=∫y∞​(yˉ​−y)w(yˉ​,S)dyˉ​ ≥ b,(5)

and the search threshold Uˉ(S)\bar U(S)Uˉ(S) is the solution of ΦS(Uˉ(S))=b\Phi_S(\bar U(S)) = bΦS​(Uˉ(S))=b, equation (4). Theorem 3 gives the resulting demand,

qiso(S)=qim(S)(1−H(Uˉ(S),S)),i∈S, S≠N.q_i^{so}(S) = q_i^m(S)\big(1 - H(\bar U(S), S)\big), \qquad i \in S,\ S \ne N.qiso​(S)=qim​(S)(1−H(Uˉ(S),S)),i∈S, S=N.

When the store carries everything there is nothing to search for (p. 28: "If the retailer offers the full assortment, there will be no consumer search"), so qiso(N)=qim(N)q_i^{so}(N) = q_i^m(N)qiso​(N)=qim​(N).

Variant iii earns margin mim_imi​ per unit and costs c(qi)c(q_i)c(qi​) to carry, where ccc is concave and increasing (economies of scale). The retailer's profit is

π(S)=∑i∈S[mi qiso(S)−c(qiso(S))],(1), (13).\pi(S) = \sum_{i\in S}\big[m_i\, q_i^{so}(S) - c(q_i^{so}(S))\big], \qquad (1),\ (13).π(S)=i∈S∑​[mi​qiso​(S)−c(qiso​(S))],(1), (13).

As printed, the second line of (13) writes qiso(S)(1−H(Uˉ(S),S))q_i^{so}(S)(1 - H(\bar U(S), S))qiso​(S)(1−H(Uˉ(S),S)) where qim(S)(1−H(Uˉ(S),S))q_i^m(S)(1 - H(\bar U(S), S))qim​(S)(1−H(Uˉ(S),S)) is meant; the definition uses the latter.

Formalization targets

Goal: Theorem 6 (p. 18)

If the full assortment is profitable, π(N)>0\pi(N) > 0π(N)>0, then there is a bˉ>0\bar b > 0bˉ>0 such that

πb(S)≤πb(N)for every S⊆N and every search cost 0<b≤bˉ.\pi_b(S) \le \pi_b(N) \qquad \text{for every } S \subseteq N \text{ and every search cost } 0 < b \le \bar b .πb​(S)≤πb​(N)for every S⊆N and every search cost 0<b≤bˉ.

The threshold bˉ\bar bbˉ depends on the data v,v0,μ,m,cv, v_0, \mu, m, cv,v0​,μ,m,c but is one number for all assortments SSS.

Milestones

  1. w(⋅,S)w(\cdot, S)w(⋅,S) is the density of max⁡i∈SˉUi\max_{i\in\bar S} U_imaxi∈Sˉ​Ui​ (proof of Theorem 3, p. 11).
  2. For S≠NS \ne NS=N, the gain ΦS(y)\Phi_S(y)ΦS​(y) is finite and strictly decreasing in yyy (proof of Theorem 3, p. 11).
  3. For S≠NS \ne NS=N and b>0b > 0b>0, Uˉ(S)\bar U(S)Uˉ(S) is the unique solution of (4) (Theorem 3, p. 11).
  4. Uˉ(S)\bar U(S)Uˉ(S) is strictly decreasing in bbb (proof of Theorem 6, p. 18).
  5. Uˉ(S)→−∞\bar U(S) \to -\inftyUˉ(S)→−∞ as b→∞b \to \inftyb→∞ and Uˉ(S)→∞\bar U(S) \to \inftyUˉ(S)→∞ as b→0+b \to 0^+b→0+ (p. 18).
  6. H(Uˉ(S),S)→1H(\bar U(S), S) \to 1H(Uˉ(S),S)→1 and qiso(S)→0q_i^{so}(S) \to 0qiso​(S)→0 as b→0+b \to 0^+b→0+ (p. 18).

Significance

Theorem 6 separates the overlapping model sharply from the no-search MNL model and from the paper's independent-assortment model. In those two models every variant in an optimal assortment must earn a strictly positive profit, because dropping a loss-making variant only shifts demand to the remaining ones. In the overlapping model the assortment also controls whether consumers leave: carrying an unprofitable variant can pay because it removes a reason to search. Theorem 6 is the extreme form of this effect, stating that the motive to prevent search can override the profitability of individual variants entirely (§4.3, p. 17).

The paper gives only a short sketch of the proof. The steps it asserts without proof are the existence and uniqueness of Uˉ(S)\bar U(S)Uˉ(S), its monotonicity and limits in bbb, and the explicit density www. This mission states each of them as a milestone. No machine-checked version of any of these results is known.

Difficulty

The search threshold Uˉ(S)\bar U(S)Uˉ(S) is defined only implicitly, through an improper integral of the density of a maximum of Gumbel variables. Each property the theorem needs (that the equation has exactly one root, that the root moves monotonically with bbb, and that it escapes to +∞+\infty+∞ as b→0b \to 0b→0) has to be read off the integral. That means controlling ΦS\Phi_SΦS​ on the whole real line: its convergence, its strict decrease, and its limits at ±∞\pm\infty±∞.

The final comparison is not a direct consequence of qiso(S)→0q_i^{so}(S) \to 0qiso​(S)→0 either. A variant's profit miq−c(q)m_i q - c(q)mi​q−c(q) tends to −c(0)-c(0)−c(0) rather than to 000, so the conclusion depends on the sign of c(0)c(0)c(0), and bˉ\bar bbˉ must be chosen uniformly over all 2n2^n2n assortments.

Formalization scope

The shared model is the definition AssortSearch.FullAssort.Model, which builds on the published MNL share RetailVariety.Structure.share. Its conventions:

  • Variants are Fin n, numbered from 0; the paper's variant iii is index i−1i-1i−1. The no-purchase option is not a variant: v0v_0v0​ is a separate positive real.
  • Preferences vi>0v_i > 0vi​>0, v0>0v_0 > 0v0​>0 and the scale μ>0\mu > 0μ>0 are free data ("for any given set of preference"). γ\gammaγ is Real.eulerMascheroniConstant.
  • Assortments are arbitrary Finset (Fin n), including ∅\emptyset∅, whose profit is 000.
  • w(yˉ,S)w(\bar y, S)w(yˉ​,S) is given by its explicit formula (VˉS/μ) λ(yˉ)exp⁡(−λ(yˉ)VˉS)(\bar V_S/\mu)\,\lambda(\bar y)\exp(-\lambda(\bar y)\bar V_S)(VˉS​/μ)λ(yˉ​)exp(−λ(yˉ​)VˉS​) with VˉS=∑i∈Sˉvi\bar V_S = \sum_{i\in\bar S} v_iVˉS​=∑i∈Sˉ​vi​. Milestone 1 proves that it is the density of the maximum.
  • ΦS(u)\Phi_S(u)ΦS​(u) is the Lebesgue integral over (u,∞)(u, \infty)(u,∞). Its integrability is part of milestone 2, so no result relies on Lean's value 000 for a non-integrable function.
  • Uˉ(S)\bar U(S)Uˉ(S) is defined as inf⁡{u:ΦS(u)≤b}\inf\{u : \Phi_S(u) \le b\}inf{u:ΦS​(u)≤b}; milestone 3 proves that it solves (4). The demand for S=NS = NS=N is the no-search MNL share, and Uˉ(N)\bar U(N)Uˉ(N) is never used.
  • The cost c:R→Rc : \mathbb R \to \mathbb Rc:R→R is concave and monotone on [0,1][0,1][0,1], where all demands lie. The hypothesis c(0)≥0c(0) \ge 0c(0)≥0 is added to Theorem 6; without it the theorem fails (a counterexample is in the goal's statement).
  • Search costs are positive: 0<b≤bˉ0 < b \le \bar b0<b≤bˉ.
  • The margins mim_imi​ are reals that satisfy §3's standing assumption of monotone margins (p. 6): mj≥mkm_j \ge m_kmj​≥mk​ whenever vj≥vkv_j \ge v_kvj​≥vk​ (equivalently uj−pj≥uk−pku_j - p_j \ge u_k - p_kuj​−pj​≥uk​−pk​). §3's labelling of variants by decreasing net utility is a labelling without loss of generality and is not encoded, since assortments are arbitrary subsets.

Several formalizations would make the goal trivial or empty, and each is excluded:

  • taking Uˉ(S)\bar U(S)Uˉ(S) as a free parameter instead of the solution of (4);
  • concluding ∃bˉ\exists \bar b∃bˉ without bˉ>0\bar b > 0bˉ>0;
  • defining qso(N)q^{so}(N)qso(N) through a junk value of Uˉ(N)\bar U(N)Uˉ(N);
  • dropping the hypothesis π(N)>0\pi(N) > 0π(N)>0;
  • fixing a specific cost function.

A complete development needs the integral calculus of the Gumbel tail (∫u∞(yˉ−u)w\int_u^\infty (\bar y - u) w∫u∞​(yˉ​−u)w and its derivative in uuu) and the product-measure computation of the law of a maximum of independent variables. Both are reusable beyond this paper, for instance in other search and extreme-value models. Contributions are welcome on any milestone, and on helper lemmas such as the continuity of ΦS\Phi_SΦS​ and its limits at ±∞\pm\infty±∞.

Selected references

  • G. P. Cachon, C. Terwiesch, Y. Xu, Retail Assortment Planning in the Presence of Consumer Search, Wharton working paper, December 20, 2002; published in Manufacturing & Service Operations Management 7(4):330–346, 2005. https://doi.org/10.1287/msom.1050.0088
  • G. van Ryzin, S. Mahajan, On the Relationship Between Inventory Costs and Variety Benefits in Retail Assortments, Management Science 45(11):1496–1509, 1999. https://doi.org/10.1287/mnsc.45.11.1496
  • D. McFadden, Conditional Logit Analysis of Qualitative Choice Behavior, in P. Zarembka (ed.), Frontiers in Econometrics, Academic Press, 1974.
10 thms1 active userReviewed
Graph TheoryProbabilityStochastic Systems·Captain: mikedeng1

Graphon Mean Field Systems V: On Percolated Graphs Sampled From a Lipschitz Graphon the Mean-Square Error Is at Most κ(q)/(nβₙ)^{1/q}Research Paper

Motivation

Large systems of interacting diffusions — neurons, oscillators, agents in a network game, particles in a queueing or load-balancing model — are usually analysed through their mean-field limit: as the population grows, each particle interacts with the empirical law of the others, and the system is replaced by a single nonlinear (McKean–Vlasov) equation. The classical theory assumes that every particle interacts with every other one with the same strength. Real networks are heterogeneous and sparse: particle iii interacts with particle jjj only through an edge of a random graph, and the number of neighbours of a particle is much smaller than the population size.

Bayraktar, Chakraborty and Wu (Ann. Appl. Probab. 33(5), 2023) describe the heterogeneity by a graphon and prove laws of large numbers and rates of convergence for such systems. This mission is the fifth of a series formalizing that paper. Its subject is the paper's quantitative result for not-so-dense graphs: when the edges are sampled from a Lipschitz graphon and then thinned by a sparsity factor βn\beta_nβn​, the average mean-square distance between each particle and its graphon limit is of order (nβn)−1/q(n\beta_n)^{-1/q}(nβn​)−1/q for every q>1q>1q>1.

Setting

Let I=[0,1]I=[0,1]I=[0,1] with Lebesgue measure. A graphon is a measurable symmetric function G:I×I→[0,1]G:I\times I\to[0,1]G:I×I→[0,1]. Fix a horizon T>0T>0T>0 and write Cd=C([0,T]:Rd)\mathcal C_d=C([0,T]:\mathbb R^d)Cd​=C([0,T]:Rd) with the norm ∥x∥∗,T=sup⁡0≤s≤T∣xs∣\|x\|_{*,T}=\sup_{0\le s\le T}|x_s|∥x∥∗,T​=sup0≤s≤T​∣xs​∣.

One probability space carries, for every label u∈Iu\in Iu∈I, an initial state Xu(0)X_u(0)Xu​(0) with law μu(0)\mu_u(0)μu​(0) and a standard ddd-dimensional Brownian motion BuB_uBu​, the whole family {Xu(0),Bu:u∈I}\{X_u(0),B_u:u\in I\}{Xu​(0),Bu​:u∈I} being independent. The limit system (4.2) is a continuum of diffusions, one per label:

Xu(t)=Xu(0)+∫0t ⁣∫I ⁣∫Rdb(Xu(s),x) G(u,v) μv,s(dx) dv ds+∫0tσ(Xu(s)) dBu(s),μu,t=L(Xu(t)).X_u(t)=X_u(0)+\int_0^t\!\int_I\!\int_{\mathbb R^d}b(X_u(s),x)\,G(u,v)\,\mu_{v,s}(dx)\,dv\,ds+\int_0^t\sigma(X_u(s))\,dB_u(s),\qquad \mu_{u,t}=\mathcal L(X_u(t)).Xu​(t)=Xu​(0)+∫0t​∫I​∫Rd​b(Xu​(s),x)G(u,v)μv,s​(dx)dvds+∫0t​σ(Xu​(s))dBu​(s),μu,t​=L(Xu​(t)).

The particles are independent but not identically distributed; their laws are coupled through the dvdvdv-integral.

The not-so-dense nnn-particle system (4.1) uses the same initial states and Brownian motions at the labels i/ni/ni/n:

Xin(t)=Xi/n(0)+∫0t1nβn∑j=1nξijn b(Xin(s),Xjn(s)) ds+∫0tσ(Xin(s)) dBi/n(s),i=1,…,n.X^n_i(t)=X_{i/n}(0)+\int_0^t\frac{1}{n\beta_n}\sum_{j=1}^n\xi^n_{ij}\,b(X^n_i(s),X^n_j(s))\,ds+\int_0^t\sigma(X^n_i(s))\,dB_{i/n}(s),\qquad i=1,\dots,n.Xin​(t)=Xi/n​(0)+∫0t​nβn​1​j=1∑n​ξijn​b(Xin​(s),Xjn​(s))ds+∫0t​σ(Xin​(s))dBi/n​(s),i=1,…,n.

The edges are percolated samples of GGG (Condition 4.3): ξijn=ξjin∼Bernoulli(βnG(i/n,j/n))\xi^n_{ij}=\xi^n_{ji}\sim\mathrm{Bernoulli}(\beta_nG(i/n,j/n))ξijn​=ξjin​∼Bernoulli(βn​G(i/n,j/n)), independent for 1≤i≤j≤n1\le i\le j\le n1≤i≤j≤n and independent of the noise. A particle has about nβnn\beta_nnβn​ neighbours, which is why the interaction is scaled by 1/(nβn)1/(n\beta_n)1/(nβn​).

The standing hypotheses are Condition 4.1 — u↦μu(0)u\mapsto\mu_u(0)u↦μu​(0) measurable with bounded second moments; bbb bounded and Lipschitz; σ\sigmaσ bounded, Lipschitz and invertible with bounded inverse; βn∈(0,1]\beta_n\in(0,1]βn​∈(0,1] and nβn→∞n\beta_n\to\inftynβn​→∞ — and Condition 2.3: on finitely many intervals I1,…,INI_1,\dots,I_NI1​,…,IN​ covering III, the map u↦μu(0)u\mapsto\mu_u(0)u↦μu​(0) is Lipschitz in the Wasserstein distance W2W_2W2​, and GGG is Lipschitz on every block Ii×IjI_i\times I_jIi​×Ij​.

Formalization targets

Goal: Theorem 4.2

Under Conditions 2.3, 4.1 and 4.3, for each q∈(1,∞)q\in(1,\infty)q∈(1,∞) there is κ(q)∈(0,∞)\kappa(q)\in(0,\infty)κ(q)∈(0,∞) with

1n∑i=1nE∥Xin−Xi/n∥∗,T2≤κ(q)(nβn)1/qfor all n∈N.\frac1n\sum_{i=1}^n\mathbb E\big\|X^n_i-X_{i/n}\big\|_{*,T}^2\le\frac{\kappa(q)}{(n\beta_n)^{1/q}}\qquad\text{for all }n\in\mathbb N .n1​i=1∑n​E​Xin​−Xi/n​​∗,T2​≤(nβn​)1/qκ(q)​for all n∈N.

The constant depends on qqq and on the data, never on nnn. The bound is on the average over particles, not on the maximum.

Milestones

  1. Lemma 7.1: all moments of the coupling error are bounded uniformly in nnn and iii.
  2. (7.10): E[ξijn∣Xjn(s)−Xj/n(s)∣2]≤(2E∣Xjn(s)−Xj/n(s)∣2+κ(q)(nβn)−1/q)βnG(i/n,j/n)\mathbb E[\xi^n_{ij}|X^n_j(s)-X_{j/n}(s)|^2]\le\big(2\mathbb E|X^n_j(s)-X_{j/n}(s)|^2+\kappa(q)(n\beta_n)^{-1/q}\big)\beta_nG(i/n,j/n)E[ξijn​∣Xjn​(s)−Xj/n​(s)∣2]≤(2E∣Xjn​(s)−Xj/n​(s)∣2+κ(q)(nβn​)−1/q)βn​G(i/n,j/n).
  3. (7.15): the interaction error Rsn,2R^{n,2}_sRsn,2​ is at most (κ/n)∑jE∣Xjn(s)−Xj/n(s)∣2+κ(q)(nβn)−1/q(\kappa/n)\sum_j\mathbb E|X^n_j(s)-X_{j/n}(s)|^2+\kappa(q)(n\beta_n)^{-1/q}(κ/n)∑j​E∣Xjn​(s)−Xj/n​(s)∣2+κ(q)(nβn​)−1/q.
  4. Theorem 2.1(b) for (4.2): W2,T(μu,μv)≤κ∣u−v∣W_{2,T}(\mu_u,\mu_v)\le\kappa|u-v|W2,T​(μu​,μv​)≤κ∣u−v∣ for u,vu,vu,v in one interval IiI_iIi​.
  5. §7.3 display: the discretization error Rsn,4R^{n,4}_sRsn,4​ is at most κ/n2\kappa/n^2κ/n2.

Significance

Theorem 4.2 gives an explicit rate for the convergence of a particle system on a sparse random graph to its graphon limit. Its companion, Theorem 4.1 (mission IV), proves the law of large numbers under weaker hypotheses but without a rate. The rate says that the effective sample size is nβnn\beta_nnβn​, the typical degree, not nnn: graphs with βn→0\beta_n\to0βn​→0 still approximate the graphon system as long as the degree grows. Results of this type justify replacing a large heterogeneous network by its graphon limit in control and game problems, such as graphon mean-field games, where the error of the replacement must be quantified.

The paper proves the theorem by hand; none of its results has a machine-checked proof. The formalization produces precise statements of the not-so-dense model, of the estimates of §7.1 that hold for any sparsity, and of the discretization bound, each usable on its own. The milestones (7.10) and (7.15) are the paper's key estimates for sparse interactions and are shared with mission IV.

Difficulty

The obvious approach is to compare each XinX^n_iXin​ with Xi/nX_{i/n}Xi/n​ by Gronwall's inequality, after bounding the difference of drifts. That step fails at the interaction term: the edge ξijn\xi^n_{ij}ξijn​ and the error Xjn−Xj/nX^n_j-X_{j/n}Xjn​−Xj/n​ are dependent, because the edge enters the dynamics of particle jjj. A plain Cauchy–Schwarz bound loses a factor 1/βn1/\beta_n1/βn​ and gives nothing when βn→0\beta_n\to0βn​→0. Estimates (7.10) and (7.15) quantify how weak this dependence is, and they are where the hypotheses that σ\sigmaσ is invertible with bounded inverse and that nβn→∞n\beta_n\to\inftynβn​→∞ are needed. The discretization term needs the Lipschitz regularity of the limit laws in the label (Theorem 2.1(b)), itself a fixed-point estimate on a space of measure-valued maps.

Formalization scope

The Lean development uses Mathlib's unitInterval for III, Fin d → ℝ with its sup norm for Rd\mathbb R^dRd (all constants are existential, so the choice of norm does not matter), time ℝ≥0, and the path space Cd\mathcal C_dCd​ as continuous maps on [0,T][0,T][0,T] with the sup norm. Brownian motions, Itô integrals and Itô processes are those of the published Peng1990.SMP.Stochastic; Wasserstein distances are the published WassersteinDRO.Duality.wassersteinDistance, valued in [0,∞][0,\infty][0,∞]. Particle iii is Lean's i : Fin n with label (i+1)/n(i+1)/n(i+1)/n.

Committed conventions:

  • A solution of (4.2) has continuous paths, AE-measurable path maps, and path laws in the class M\mathcal MM (measurable in uuu, bounded second moments); each XuX_uXu​ is a strong solution for the filtration of (Xu(0),Bu)(X_u(0),B_u)(Xu​(0),Bu​). A solution of (4.1) is a strong solution for the filtration of the edges, the initial states and the Brownian motions of the nnn labels.
  • The edges are {0,1}\{0,1\}{0,1}-valued, symmetric, with the diagonal included, mutually independent over i≤ji\le ji≤j, and independent of the σ-algebra of the noise.
  • Condition 4.1(d) is required for n≥1n\ge1n≥1; the invertibility of σ\sigmaσ is that of a d×dd\times dd×d matrix; the intervals of Condition 2.3 are order-connected subsets of III.
  • Expectations of nonnegative quantities are lower Lebesgue integrals, so no expectation is silently 000 for a non-integrable integrand.

The solution predicates are not vacuous: the zero-coefficient system Xu(t)=Xu(0)X_u(t)=X_u(0)Xu​(t)=Xu​(0) satisfies both, Condition 4.1 and Condition 4.3 are satisfiable, and every statement quantifies over solutions rather than assuming the conclusion. The constant κ(q)\kappa(q)κ(q) is chosen after qqq and before nnn, so it cannot depend on nnn, and the bound is the paper's average over iii, not a weaker or a stronger form. At n=0n=0n=0 both sides of every rate are 000.

The paper's proofs are in its §5 (Theorem 2.1), §7.1 (Lemma 7.1, (7.10), (7.15)) and §7.3 (Theorem 4.2). Stochastic analysis at that level — Burkholder–Davis–Gundy and Rosenthal inequalities, Girsanov's theorem for these SDEs, Gronwall arguments for measure-valued fixed points — is mostly absent from Mathlib, and contributions of such infrastructure are welcome, as are proofs of any milestone.

Selected references

  • E. Bayraktar, S. Chakraborty, R. Wu, Graphon mean field systems, Ann. Appl. Probab. 33(5):3587–3619, 2023. https://doi.org/10.1214/22-AAP1901
  • L. Lovász, Large Networks and Graph Limits, AMS Colloquium Publications 60, 2012. https://doi.org/10.1090/coll/060
  • B. Bollobás, C. Borgs, J. Chayes, O. Riordan, Percolation on dense graph sequences, Ann. Probab. 38(1):150–183, 2010. https://doi.org/10.1214/09-AOP478
  • S. Peng, A general stochastic maximum principle for optimal control problems, SIAM J. Control Optim. 28(4):966–979, 1990. https://doi.org/10.1137/0328054
11 thms1 active userReviewed
ProbabilityStochastic Systems·Captain: mikedeng1

Steady-State Analysis of the Join-the-Shortest-Queue Model in the Halfin-Whitt Regime 1: Diffusion-Scaled Stationary Queue Lengths Have Expectations Bounded Uniformly in nResearch Paper

Motivation

Join-the-shortest-queue (JSQ) is the load-balancing rule in which each arriving job is sent to a server with the fewest jobs. It is a baseline in the design of server farms and data centers, and in the many-server Halfin–Whitt regime, where the load per server approaches one at rate 1/n1/\sqrt n1/n​, it combines near-full utilization with vanishing waiting. Eschenfeldt and Gamarnik (arXiv:1502.00999, 2015) proved that the suitably scaled JSQ process converges, over bounded time intervals, to a two-dimensional reflected diffusion; Mukherjee, Borst, van Leeuwaarden and Whiting (J. Appl. Probab. 53, 2016) gave the form of the result that Braverman quotes as Theorem 1. Process-level convergence over bounded time intervals says nothing about the stationary distribution. To conclude that the stationary distributions converge as well, one needs tightness of the scaled stationary distributions, uniformly in the number of servers nnn.

Braverman (arXiv:1801.05121, Math. Oper. Res. 45(3), 2020) supplied this, together with positive recurrence of the diffusion limit. This mission formalizes the first of the two main results, Theorem 2, and the chain of lemmas the paper uses to prove it.

Setting

There are nnn identical servers, each with its own infinite buffer. Jobs arrive as a Poisson process of rate nλn\lambdanλ and service times are i.i.d. exponential with mean one. An arriving job joins a server with the fewest jobs, ties broken arbitrarily. The Halfin–Whitt regime is

λ=1−β/n,β>0 fixed.\lambda = 1 - \beta/\sqrt n,\qquad \beta > 0 \text{ fixed}.λ=1−β/n​,β>0 fixed.

For i≥1i \ge 1i≥1, QiQ_iQi​ is the number of servers with at least iii jobs. The process Q=(Q1,Q2,… )Q = (Q_1, Q_2, \dots)Q=(Q1​,Q2​,…) is a continuous-time Markov chain on

S={q∈{0,1,…,n}∞ ∣ qi≥qi+1, ∑iqi<∞}.S = \Big\{ q \in \{0,1,\dots,n\}^\infty \ \Big|\ q_i \ge q_{i+1},\ \textstyle\sum_i q_i < \infty \Big\}.S={q∈{0,1,…,n}∞ ​ qi​≥qi+1​, ∑i​qi​<∞}.

Its generator GQG_QGQ​ sends qqq to q+e(i)q + e^{(i)}q+e(i) at rate nλn\lambdanλ when q1=⋯=qi−1=n>qiq_1 = \dots = q_{i-1} = n > q_iq1​=⋯=qi−1​=n>qi​, and to q−e(i)q - e^{(i)}q−e(i) at rate qi−qi+1q_i - q_{i+1}qi​−qi+1​. The fluid-scaled coordinates are X1=(Q1−n)/n≤0X_1 = (Q_1 - n)/n \le 0X1​=(Q1​−n)/n≤0, the negative of the fraction of idle servers, and Xi=Qi/nX_i = Q_i/nXi​=Qi​/n for i≥2i \ge 2i≥2. The diffusion-scaled coordinates are nXi\sqrt n X_in​Xi​. Expectations E\mathbb EE are taken under a stationary distribution π\piπ of the chain.

The proof works with functions on Ω=(−∞,0]×[0,∞)\Omega = (-\infty, 0] \times [0, \infty)Ω=(−∞,0]×[0,∞) and with the first-order operator

Lf(x)=(−x1+x2−β/n)f1(x)−x2f2(x),Lf(x) = (-x_1 + x_2 - \beta/\sqrt n) f_1(x) - x_2 f_2(x),Lf(x)=(−x1​+x2​−β/n​)f1​(x)−x2​f2​(x),

where fif_ifi​ is the partial derivative in xix_ixi​, taken one-sided on ∂Ω\partial\Omega∂Ω.

Formalization targets

Goal: Theorem 2

For each β>0\beta > 0β>0 there is a constant C(β)C(\beta)C(β) such that for all n≥1n \ge 1n≥1 with β<n\beta < \sqrt nβ<n​ and every stationary distribution,

E∣nXi∣≤C(β),  i=1,2,E∣nXi∣≤C(β),  i≥3.\mathbb E\big|\sqrt n X_i\big| \le C(\beta),\ \ i = 1, 2,\qquad \mathbb E\big|n X_i\big| \le C(\beta),\ \ i \ge 3.E​n​Xi​​≤C(β),  i=1,2,E​nXi​​≤C(β),  i≥3.

The constant depends on β\betaβ only, not on nnn, the stationary distribution, or iii.

Milestones, in the order the proof uses them

  1. Lemma 1: E GQf(Q)=0\mathbb E\,G_Q f(Q) = 0EGQ​f(Q)=0 whenever E∣f(Q)∣<∞\mathbb E|f(Q)| < \inftyE∣f(Q)∣<∞.
  2. Lemma 2: EQ1=nλ\mathbb E Q_1 = n\lambdaEQ1​=nλ and EQi=nλ P(Q1=⋯=Qi−1=n)\mathbb E Q_i = n\lambda\,\mathbb P(Q_1 = \dots = Q_{i-1} = n)EQi​=nλP(Q1​=⋯=Qi−1​=n).
  3. Lemma 3: the generator, applied to f(x1,x2)f(x_1, x_2)f(x1​,x2​), equals Lf(x)+(f2−f1)(x)λ1(x1=0)+ε(x)Lf(x) + (f_2 - f_1)(x)\lambda 1(x_1 = 0) + \varepsilon(x)Lf(x)+(f2​−f1​)(x)λ1(x1​=0)+ε(x) with an explicit remainder ε\varepsilonε built from second weak derivatives.
  4. Lemmas 5 and 6: the curve Γ(κ)\Gamma^{(\kappa)}Γ(κ) and the hitting time τ(x)\tau(x)τ(x) of the fluid model; existence and uniqueness, finiteness, derivative formulas and monotonicity in κ\kappaκ.
  5. Lemma 7 and Lemma 4: for κ>β\kappa > \betaκ>β, the explicit function f∗f^*f∗ of (4.14) solves
Lf(x)=−((x2−κ/n)∨0) on Ω,f1(0,x2)=f2(0,x2),Lf(x) = -\big((x_2 - \kappa/\sqrt n)\vee 0\big)\ \text{on } \Omega,\qquad f_1(0, x_2) = f_2(0, x_2),Lf(x)=−((x2​−κ/n​)∨0) on Ω,f1​(0,x2​)=f2​(0,x2​),

with f11∗,f12∗,f22∗≥0f^*_{11}, f^*_{12}, f^*_{22} \ge 0f11∗​,f12∗​,f22∗​≥0, vanishing second derivatives for x2≤κ/nx_2 \le \kappa/\sqrt nx2​≤κ/n​, and explicit upper bounds of order n/β\sqrt n/\betan​/β above that level. 6. (3.18): E((X2−κ/n)∨0)≤1βn(12+6κκ−β)P(X2≥κ/n−1/n)\mathbb E\big((X_2 - \kappa/\sqrt n)\vee 0\big) \le \frac{1}{\beta\sqrt n}\big(12 + \frac{6\kappa}{\kappa - \beta}\big)\mathbb P(X_2 \ge \kappa/\sqrt n - 1/n)E((X2​−κ/n​)∨0)≤βn​1​(12+κ−β6κ​)P(X2​≥κ/n​−1/n). 7. Theorem 2, (2.2): the i=1,2i = 1, 2i=1,2 half of the goal.

Significance

Theorem 2 shows that in steady state the number of idle servers and the number of servers holding two or more jobs are O(n)O(\sqrt n)O(n​) on average, and the number of servers holding i≥3i \ge 3i≥3 jobs is O(1)O(1)O(1) on average, uniformly in nnn. Combined with the process-level limit and the positive recurrence of the limiting diffusion (Theorem 3, the companion mission), it yields Proposition 1: the diffusion-scaled stationary distribution converges to the stationary distribution of the diffusion. This justifies the diffusion as a steady-state approximation for JSQ.

The method is the generator-comparison or Stein's-method approach: a Lyapunov function solving a first-order PDE for the fluid model, with explicit bounds on its derivatives. It applies to other many-server systems, and the paper is an example of that approach on an infinite-dimensional chain whose diffusion limit is two-dimensional.

The result is proved in the paper; no machine-checked proof is known to exist. The work remaining is to formalize the known proof: the adjoint relation for a countable-state chain with bounded rates, the Taylor expansion of the generator, and the explicit analysis of f∗f^*f∗ through the fluid curves.

Difficulty

The obvious route is to bound P(Q1=⋯=Qi−1=n)\mathbb P(Q_1 = \dots = Q_{i-1} = n)P(Q1​=⋯=Qi−1​=n) in Lemma 2 directly. The paper notes that no good handle on these probabilities exists. The PDE route needs a solution whose second derivatives are bounded by O(n)O(\sqrt n)O(n​). General results that give such bounds require a continuous fluid vector field, and the JSQ fluid model is discontinuous because of the reflection at x1=0x_1 = 0x1​=0. So f∗f^*f∗ is constructed explicitly from the fluid paths, and its derivatives must be bounded by hand through the implicitly defined τ(x)\tau(x)τ(x) and ν∗(x1)\nu^*(x_1)ν∗(x1​).

A second difficulty is dimensional: the chain is infinite-dimensional while the PDE is two-dimensional, so the remainder of Lemma 3 contains a term with q3q_3q3​ that is controlled only by combining Lemma 2 with the monotonicity of f2∗f^*_2f2∗​.

Formalization scope

  • Model. States are functions ℕ → ℕ with 0-based indices: Lean q.1 i is the paper's qi+1q_{i+1}qi+1​, so the goal's "i≥3i \ge 3i≥3" is Lean index ≥2\ge 2≥2. Each state is bounded by nnn, nonincreasing and eventually zero. The generator genQ is the displayed formula on p. 6.
  • Stationarity. No process is constructed. A stationary distribution is a probability π\piπ on SSS satisfying global balance πGQ=0\pi G_Q = 0πGQ​=0, which characterizes stationarity for this chain with bounded rates. Such a π\piπ exists, since the chain is positive recurrent for λ<1\lambda < 1λ<1 (Bramson 2011), so statements quantifying over it are not vacuous. Expectations are sums against π\piπ. Lemma 1 keeps its integrability hypothesis, and all other expectations are of bounded functions.
  • Parameters. λ=1−β/n\lambda = 1 - \beta/\sqrt nλ=1−β/n​ with 0<β<n0 < \beta < \sqrt n0<β<n​, so that λ>0\lambda > 0λ>0; λ\lambdaλ is free in Lemmas 1–2. Constants are quantified as ∃C\exists C∃C after β\betaβ and before nnn.
  • Misprint. The printed Theorem 2 omits E\mathbb EE. Read as an almost-sure bound it is false, so the expectation version, which the paper proves, is the target.
  • Derivatives. Points of Ω\OmegaΩ are pairs in ℝ × ℝ. Partial derivatives are named functions asserted to be one-sided derivatives along coordinate lines (HasDerivWithinAt on Set.Iic 0 or Set.Ici 0). Weak second derivatives are a.e. derivatives of absolutely continuous functions, in integral form on compact intervals. A "solution" of the PDE whose derivatives are the junk values of a derivative operator is thereby excluded. The PDE needs genuine derivatives, and Lemma 7 is about the explicit f∗f^*f∗ of (4.14), not about an arbitrary function.
  • Fluid objects. ν∗,η∗\nu^*, \eta^*ν∗,η∗ are the components of the unique solution of (4.8). τ\tauτ is valued in WithTop ℝ, with ∞\infty∞ as ⊤, so a missing solution is never read as 000. f∗f^*f∗ is an if cascade, and Lemma 7 states the agreement of its branches. f∗f^*f∗ is defined by the closed form (4.14), not through the fluid model (4.1).
  • Not included. Lemma 9 (via the Lambert W function, which Mathlib lacks), the bootstrapping argument of App. A.4 (not a numbered result), Theorems 1 and 3, and Proposition 1.

Contributions of reusable infrastructure are welcome: the adjoint relation for countable-state chains with bounded rates, and FTC/Taylor lemmas for absolutely continuous first derivatives.

Selected references

  • A. Braverman, Steady-State Analysis of the Join-the-Shortest-Queue Model in the Halfin–Whitt Regime, Math. Oper. Res. 45(3), 2020; arXiv:1801.05121v2. https://arxiv.org/abs/1801.05121
  • P. Eschenfeldt, D. Gamarnik, Join the Shortest Queue with Many Servers. The Heavy Traffic Asymptotics, arXiv preprint, 2015. https://arxiv.org/abs/1502.00999
  • D. Mukherjee, S. C. Borst, J. S. H. van Leeuwaarden, P. A. Whiting, Universality of Load Balancing Schemes on the Diffusion Scale, J. Appl. Probab. 53, 1111–1124, 2016. https://projecteuclid.org/euclid.jap/1481132840
  • M. Bramson, Stability of Join the Shortest Queue Networks, Ann. Appl. Probab. 21, 1568–1625, 2011. https://doi.org/10.1214/10-AAP726
  • S. Halfin, W. Whitt, Heavy-Traffic Limits for Queues with Many Exponential Servers, Oper. Res. 29(3), 1981. https://doi.org/10.1287/opre.29.3.567
17 thms1 active userReviewed
Graph TheoryProbabilityStochastic Systems·Captain: mikedeng1

Graphon Mean Field Systems IV: Law of Large Numbers for Particle Systems on Percolated Not-So-Dense Graphs With nβₙ → ∞Research Paper

Motivation

Large networks are often modeled through a graphon, a measurable function that records how likely two labeled vertices are to interact. For interacting diffusions, this lets one ask whether a finite random network has a predictable population limit. Bayraktar, Chakraborty and Wu study that question for graphon weighted systems, including graphs obtained by independently retaining potential edges with a sparsity parameter βn\beta_nβn​. Their not-so-dense regime allows βn\beta_nβn​ to decrease, provided the typical interaction scale nβnn\beta_nnβn​ diverges. Theorem 4.1 is the paper's law of large numbers for this regime; its two conclusions distinguish convergence of the population's empirical law from mean-square closeness of each particle to a coupled graphon limit. Bayraktar, Chakraborty and Wu, 2023.

The authors present three main convergence results: continuity of the continuum graphon system in Section 2, a dense-system law of large numbers in Section 3, and the percolated-system law of large numbers in Section 4. Section 7 establishes the latter. The graphon and cut-metric framework used here follows the graph-limit treatment cited by the paper. Bayraktar, Chakraborty and Wu, 2023; Lovász, 2012.

Setting

Let I=[0,1]I=[0,1]I=[0,1]. A graphon is a symmetric measurable G:I2→[0,1]G:I^2\to[0,1]G:I2→[0,1]. The cut norm ∥W∥□\|W\|_\square∥W∥□​ measures the largest absolute integral of WWW over a measurable rectangle S×US\times US×U. A sequence of step graphons GnG_nGn​ approaches GGG when ∥Gn−G∥□→0\|G_n-G\|_\square\to0∥Gn​−G∥□​→0. The related operator norm ∥W∥∞→1\|W\|_{\infty\to1}∥W∥∞→1​ measures the effect of WWW on bounded measurable functions. Remark 2.1 says cut convergence implies convergence in this operator norm. Bayraktar, Chakraborty and Wu, 2023, Remark 2.1.

Fix a dimension d≥1d\ge1d≥1 and horizon T>0T>0T>0. A label u∈Iu\in Iu∈I has an initial state Xu(0)∈RdX_u(0)\in\mathbb R^dXu​(0)∈Rd and an independent ddd-dimensional Brownian motion BuB_uBu​. The initial states are independent of each other and of the Brownian motions. The limiting graphon particle XuX_uXu​ is a continuous path whose drift averages b(Xu(s),x)G(u,v)b(X_u(s),x)G(u,v)b(Xu​(s),x)G(u,v) against the time-sss laws of all labels vvv, and whose diffusion matrix is σ(Xu(s))\sigma(X_u(s))σ(Xu​(s)). The time-sss law at label vvv is μv,s=L(Xv(s))\mu_{v,s}=\mathcal L(X_v(s))μv,s​=L(Xv​(s)). The family of path laws μu=L(Xu)\mu_u=\mathcal L(X_u)μu​=L(Xu​) is measurable in uuu and has uniformly bounded second moments. Bayraktar, Chakraborty and Wu, 2023, (4.2).

For nnn particles, label iii means i/ni/ni/n, with 1≤i≤n1\le i\le n1≤i≤n. Each undirected edge weight ξijn=ξjin\xi^n_{ij}=\xi^n_{ji}ξijn​=ξjin​ is Bernoulli with mean βnGn(i/n,j/n)\beta_nG_n(i/n,j/n)βn​Gn​(i/n,j/n); diagonal weights are included. These edge variables are independent of the Brownian motions and initial states. The finite particle XinX_i^nXin​ uses the same initial state and Brownian motion as Xi/nX_{i/n}Xi/n​. Its drift is the sum of ξijnb(Xin,Xjn)\xi^n_{ij}b(X_i^n,X_j^n)ξijn​b(Xin​,Xjn​) divided by nβnn\beta_nnβn​; its diffusion is σ(Xin)\sigma(X_i^n)σ(Xin​). The coefficients bbb and σ\sigmaσ are bounded and Lipschitz, and σ(x)\sigma(x)σ(x) is invertible with uniformly bounded inverse. Initial laws have uniformly bounded second moments, 0<βn≤10<\beta_n\le10<βn​≤1, and nβn→∞n\beta_n\to\inftynβn​→∞. Bayraktar, Chakraborty and Wu, 2023, Conditions 4.1–4.2.

Formalization targets

Let Cd=C([0,T];Rd)\mathcal C_d=C([0,T];\mathbb R^d)Cd​=C([0,T];Rd), μn=n−1∑i=1nδXin\mu^n=n^{-1}\sum_{i=1}^{n}\delta_{X_i^n}μn=n−1∑i=1n​δXin​​, and μˉ=∫Iμu du\bar\mu=\int_I\mu_u\,duμˉ​=∫I​μu​du. The main target, Theorem 4.1, first asserts convergence in probability in the weak topology on P(Cd)\mathcal P(\mathcal C_d)P(Cd​):

μn→Pμˉ.\mu^n\xrightarrow{\mathbb P}\bar\mu.μnP​μˉ​.

This uses Condition 2.2(a): on each of finitely many intervals covering III, the initial law varies continuously in Wasserstein distance W2W_2W2​. If the graphon also satisfies Condition 2.2(b), its stated sectional continuity outside null sets, the theorem adds

1n∑i=1nE∥Xin−Xi/n∥∗,T2⟶0,∥x∥∗,T=sup⁡0≤s≤T∣xs∣.\frac1n\sum_{i=1}^{n}\mathbb E\|X_i^n-X_{i/n}\|_{*,T}^2\longrightarrow0, \qquad \|x\|_{*,T}=\sup_{0\le s\le T}|x_s|.n1​i=1∑n​E∥Xin​−Xi/n​∥∗,T2​⟶0,∥x∥∗,T​=0≤s≤Tsup​∣xs​∣.

The milestones record the paper's cut-to-operator observation, uniform coupling moments, the one-edge and two-edge estimates, the bound on the second interaction-error term, the quantitative coupling lemma, and the empirical law of large numbers for the independent limiting particles. The second displayed limit is conditional on 2.2(b); the first is not. Bayraktar, Chakraborty and Wu, 2023, Theorem 4.1.

Significance

The first conclusion identifies the population-level limit even when the sampled network has a decreasing edge-retention probability. It gives a deterministic law on path space, rather than only a limit at one time. The second conclusion is stronger: it compares each finite particle with the limiting particle built from the same initial state and Brownian motion, averaged over all labels. This comparison is available when the extra graphon continuity assumption holds. Bayraktar, Chakraborty and Wu, 2023, Theorem 4.1 and Remark 4.1.

The mathematical results were proved in the published paper. The present mission asks for machine-checked Lean proofs of its specified statements; the draft theorems themselves currently contain sorry. Reusable outputs would include the percolated-edge model, a path-law formulation of graphon diffusions, and estimates for random interactions that depend on the particle states they influence. The paper's proofs appear in Sections 5–7, with the sparse-system argument in Section 7. Bayraktar, Chakraborty and Wu, 2023.

Difficulty

The edge variables are independent of the driving noise, but a particle's state depends on its incident edges. Therefore an expression such as E[ξijn∣Xjn−Xj/n∣2]\mathbb E[\xi^n_{ij}|X_j^n-X_{j/n}|^2]E[ξijn​∣Xjn​−Xj/n​∣2] cannot be factored into the edge mean times the error moment. The same dependence affects pairs of edges in the squared interaction sum. Equations (7.10), (7.14), and (7.15) quantify the discrepancy. A second difficulty is that cut convergence controls integrals against fixed bounded tests, while the dynamics and their path laws supply tests that vary with the system. The paper's regularity conditions and quantitative estimates are needed to connect those forms of control. Bayraktar, Chakraborty and Wu, 2023, Sections 7.1–7.2.

Formalization scope

Lean represents III by Mathlib's unit interval, Rd\mathbb R^dRd by Fin d → ℝ, and continuous paths by maps from [0,T][0,T][0,T] into that vector space. The chosen vector and matrix norms are fixed sup-type norms; finite-dimensional norm equivalence preserves the limit statements and the estimates with existential constants. The horizon and dimension are explicitly positive. Particle index i : Fin n denotes the paper's (i+1)(i+1)(i+1)st particle, and empirical laws use positive sizes n+1n+1n+1 to avoid division by zero. Wasserstein distances use extended nonnegative reals, while nonnegative expectations use lower integrals.

The noise lives on one probability space. The continuum equation uses the natural filtration of its own initial state and Brownian motion. The finite equation uses the natural filtration of its edge array, sampled initial states, and Brownian motions. A solution carries a measurable continuous-path law family in the paper's class M\mathcal MM; this prevents a nonmeasurable pushforward or a default-zero integral from silently satisfying an equation. Edge symmetry, Bernoulli marginals, independence, and diagonal self-loops are explicit. The zero-coefficient system has a separate sorry-free Lean witness for both solution predicates, ruling out a vacuous solution concept.

The formalization uses the published Peng1990.SMP.Stochastic Itô interface and WassersteinDRO.Duality.wassersteinDistance definition. Contributions can prove the listed estimates, strengthen their reusable stochastic and measure-theoretic infrastructure, or connect this mission's local setting to the other graphon missions after their definitions are published. The auxiliary changed-measure particles of Section 7.1 are outside the initial milestone list. Bayraktar, Chakraborty and Wu, 2023, Sections 4 and 7.

Selected references

  • E. Bayraktar, S. Chakraborty and R. Wu, Graphon mean field systems, Annals of Applied Probability 33(5):3587–3619, 2023. DOI.
  • L. Lovász, Large Networks and Graph Limits, American Mathematical Society Colloquium Publications 60, 2012. DOI.
13 thms1 active userReviewed
PreviousPage 95 of 106Next

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