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,413 missions · 1,137 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

Open1276Completed1137All2413
CombinatoricsGraph Theory·Captain: mikedeng1

Bounding Residence Times for Atomic Dynamic Routings: If No Later Entrant Leaves the Subnetwork First, an Agent Exits by Time r0 plus the Longest Path Length plus the Earlier Agents Still InsideResearch Paper

Motivation

In atomic dynamic routing, indivisible agents (vehicles, packets, travellers) enter a network over time and travel along fixed routes. Each road admits a limited number of agents per unit of time, so agents queue at its entrance. A basic question for traffic models, packet-switching networks and dynamic routing games is how long an agent can stay in the network, its residence time. Upper bounds on residence times are what make it possible to bound queue lengths and the price of anarchy in dynamic routing games.

Greedy packet scheduling, where an edge never idles while packets wait for it, was studied by Cidon, Kutten, Mansour and Peleg (SIAM J. Comput. 1995), who showed that cycles can make delays large, and by Mansour and Patt-Shamir (J. Algorithms 1993), who bounded delivery times when every route is a shortest path. Cao, Chen, Chen and Wang (MOR 2022) prove a bound for arbitrary routes on acyclic networks. It holds under any rule that decides which waiting agent may leave next, and agents may keep entering forever. Their Theorem 1 is the goal of this mission. The paper then uses it to show that equilibria of three dynamic routing games on series-parallel networks have bounded queues (§6).

Setting

A network G=(V,E)G=(V,E)G=(V,E) is a finite directed multigraph (parallel edges allowed) that is acyclic. Every edge has unit length and unit capacity. A path is a sequence of edges, each starting where the previous one ends. LvwL_{vw}Lvw​ is the number of edges of a longest vvv-www path.

Time is discrete, 0,1,2,…0,1,2,\dots0,1,2,…. Agents form a set Δ\DeltaΔ, possibly infinite. At each time finitely many agents enter. Agent iii enters at time rir_iri​ at his origin oio_ioi​ and follows a fixed route, a nonempty oio_ioi​-did_idi​ path. The tail of each edge holds a buffer. Two rules govern the dynamics:

  • (A1) an agent reaching a vertex that is not his destination enters the buffer of his next edge at once;
  • (A2) at every integer time, whenever a buffer is nonempty, exactly one agent leaves it and reaches the head of the edge one time unit later; the others wait.

No rule says which agent leaves. Every choice consistent with (A2) is an admissible run. An agent exits GGG on reaching his destination. His residence time is his exit time minus his entry time (Definition 1).

Fix a focal agent ζ\zetaζ with entry time rrr. GζG_\zetaGζ​ is the subnetwork formed by all oζo_\zetaoζ​-dζd_\zetadζ​ paths. For an integer time r0r_0r0​:

  • Aζ,pA_{\zeta,p}Aζ,p​ is the set of agents who enter after time rrr and use an edge of GζG_\zetaGζ​;
  • Bζ,p(r0)B_{\zeta,p}(r_0)Bζ,p​(r0​) is the set of agents i≠ζi\neq\zetai=ζ who entered by time rrr, are still inside GGG at time r0r_0r0​, and use an edge of GζG_\zetaGζ​.

An agent exits GζG_\zetaGζ​ when he reaches the head of the last edge of his route lying in GζG_\zetaGζ​.

Formalization targets

Goal: Theorem 1 (p. 7)

For every admissible run and every integer r0≥rr_0\ge rr0​≥r: if no agent of Aζ,pA_{\zeta,p}Aζ,p​ exits GζG_\zetaGζ​ before ζ\zetaζ exits GGG, then

tζ ≤ r0+Loζdζ+∣Bζ,p(r0)∣,t_\zeta\ \le\ r_0+L_{o_\zeta d_\zeta}+\bigl|B_{\zeta,p}(r_0)\bigr|,tζ​ ≤ r0​+Loζ​dζ​​+​Bζ,p​(r0​)​,

where tζt_\zetatζ​ is the exit time of ζ\zetaζ.

Milestones (§3.2, p. 9)

When ζ\zetaζ is still inside GGG at r0r_0r0​, the proof follows tokens: imaginary markers that ζ\zetaζ produces from time r0r_0r0​ on, and that travel with the agents leaving buffers ((s1), (s2), p. 8). The milestones are the paper's statements about them, with L=LoζdζL=L_{o_\zeta d_\zeta}L=Loζ​dζ​​:

  1. (p1) each agent holds at most one token, and a buffer never contains more tokens than agents;
  2. (p2) the agent leaving ζ\zetaζ's buffer always holds exactly one token;
  3. (p3) after r0r_0r0​, ζ\zetaζ produces a new token only while he waits;
  4. Claim 1 the induction step of Lemma 1 at the vertices where token kkk departs right after a smaller token;
  5. Lemma 1 token kkk departs every vertex vvv of its trail inside GζG_\zetaGζ​ by time r0+(L−Lvdζ+k−1)r_0+(L-L_{vd_\zeta}+k-1)r0​+(L−Lvdζ​​+k−1);
  6. token exit token kkk exits GζG_\zetaGζ​ by time r0+(L+k−1)r_0+(L+k-1)r0​+(L+k−1).

Significance

The result. The bound does not depend on the rule that selects the leaving agent, so it covers every queueing discipline at once: FIFO, global agent priorities, and the edge-priority rules of earlier dynamic routing games. With r0=rr_0=rr0​=r and no later entrants on GζG_\zetaGζ​, it gives Corollary 1: the residence time of ζ\zetaζ is at most Loζdζ+∣Bζ,p(r)∣L_{o_\zeta d_\zeta}+|B_{\zeta,p}(r)|Loζ​dζ​​+∣Bζ,p​(r)∣. Example 1 (p. 8) shows that the bound can be attained. The freedom in r0r_0r0​ is what §6 uses to show that queue lengths stay bounded at equilibria of dynamic routing games on series-parallel networks.

Formalizing it. The theorem is proved on paper. No machine-checked version is known to exist, and the platform has no rule-free model of atomic dynamic routing. The mission adds that model: a run is a set of leave times constrained by (A1)–(A2), not a fixed simulation. It also adds the multi-token process and the paper's chain of statements leading to Theorem 1.

Difficulty

The obvious argument is wrong. One might reason that ζ\zetaζ travels at most LoζdζL_{o_\zeta d_\zeta}Loζ​dζ​​ edges and each earlier agent delays him at most once, giving the bound at once. The paper notes (p. 7) that this intuition fails even with two agents: a single agent can delay ζ\zetaζ more than ∣Δ∣|\Delta|∣Δ∣ times, because delays propagate through other buffers. Bounding the total therefore means tracking chains of delays across the network rather than charging each delay to one agent.

Formally, the hard parts are:

  • the run is only constrained, not determined;
  • the token process is defined by recursion over time, using the unique leaver of each buffer;
  • (p1)–(p3) are invariants proved by induction on time;
  • Lemma 1 combines induction on the token index with longest-path inequalities in an acyclic graph.

Formalization scope

  • Model. Network V E (tail and head maps) with finite V, E and a hypothesis Acyclic: no nonempty closed walk. Paths are edge lists. longest v w is the supremum of path lengths in N\mathbb NN; it is a maximum wherever it is used.

  • Runs. RunData holds entry times, routes, origins and destinations, and leave times leave i j in N∪{∞}\mathbb N\cup\{\infty\}N∪{∞} for the jjj-th edge of iii's route, numbered from 000. IsADR states that:

    • routes are nonempty origin–destination paths;
    • finitely many agents enter per time;
    • nobody leaves a buffer before entering it;
    • (A2) holds, as unique existence of the leaver of each nonempty buffer at each time.

    Agents form an arbitrary type and may be infinitely many. An agent may never be served (leave time ∞\infty∞), as in the paper, so the theorem also asserts that ζ\zetaζ exits.

  • Derived objects. Exit times, GζG_\zetaGζ​, Aζ,pA_{\zeta,p}Aζ,p​, Bζ,p(r0)B_{\zeta,p}(r_0)Bζ,p​(r0​) and the exit time from GζG_\zetaGζ​ are definitions computed from the run, never free parameters. ∣Bζ,p(r0)∣|B_{\zeta,p}(r_0)|∣Bζ,p​(r0​)∣ is Set.ncard of a set that is finite under IsADR.

  • Tokens. The token milestones assume that ζ\zetaζ is inside GGG at r0r_0r0​, as §3.1 requires when starting the construction. Tokens are indexed from 111. A new token gets index one plus the number produced so far. Bounds with subtraction are written additively over N\mathbb NN, e.g. s+Lvdζ+1≤r0+L+ks+L_{vd_\zeta}+1\le r_0+L+ks+Lvdζ​​+1≤r0​+L+k, which is the same inequality over the integers.

  • Ruled out. None of these is used, since each proves a different statement: a fixed tie-breaking rule (FIFO or priorities), exit times as free parameters, finitely many agents, or dropping acyclicity, which Remark 2 shows is necessary.

Contributions welcome: proofs of the token invariants, of Lemma 1, and of the goal. The network and run definitions can be reused for the paper's other results (Theorem 2 on uniform capacities, the routing games of §5–6).

Selected references

  • Z. Cao, B. Chen, X. Chen, C. Wang, Bounding residence times for atomic dynamic routings, Mathematics of Operations Research, 2022 (accepted manuscript: WRAP 160244). https://doi.org/10.1287/moor.2021.1242 , https://wrap.warwick.ac.uk/160244
  • Y. Mansour, B. Patt-Shamir, Greedy packet scheduling on shortest paths, Journal of Algorithms 14(3):449–465, 1993.
  • I. Cidon, S. Kutten, Y. Mansour, D. Peleg, Greedy packet scheduling, SIAM Journal on Computing 24(1):148–157, 1995.
  • T. Harks, B. Peis, D. Schmand, B. Tauer, L. Vargas Koch, Competitive packet routing with priority lists, ACM Transactions on Economics and Computation 6(1):4:1–4:26, 2018.
10 thms1 active userReviewed
Economics·Captain: mikedeng1

A Theory of ICOs: Diversification, Agency, and Information Asymmetry: Tokens Have a Signaling Advantage over Equity Exactly for Moderate Venture Heterogeneity and for Moderate Payoff VolatilityResearch Paper

Why financing can convey venture quality

An entrepreneur who knows more about a venture than its investors can convey information through the share of the venture that she keeps. A larger retained share exposes her to more of its uncertain payoff. The cost of this signal depends on what investors buy: an equity claim receives revenue after production costs, while a token claim receives a share of gross revenue. Chod and Lyandres examine how that difference affects the ability of high- and low-quality ventures to separate in an initial coin offering. Their analysis identifies a range of venture heterogeneity and payoff volatility in which tokens support separation and equity does not. The source for this mission is their SSRN working paper, version of December 14, 2019, §4 and Proposition 6.

This result concerns a specific information problem. Investors cannot observe venture quality before buying; an entrepreneur knows her own type. If the same financing contract is used regardless of type, investors cannot identify quality from the offering. A separating outcome instead assigns different sale fractions to the two types and makes the resulting investor beliefs consistent with each type's incentives. The paper compares the feasibility of such outcomes for the two securities, under a model in which the agency effects studied in its earlier sections are set aside.

The two-type financing game

There are two venture types, high HHH and low LLL. The entrepreneur observes the type; investors observe the fraction α\alphaα offered for sale. Type iii has uncertain shock Xi=ξi2X_i=\xi_i^2Xi​=ξi2​ with a normal approximation of mean μi\mu_iμi​ and standard deviation σ\sigmaσ, where μH>μL\mu_H>\mu_LμH​>μL​. Revenue is Ri=Xi/4R_i=X_i/4Ri​=Xi​/4. The entrepreneur pays fixed production cost CCC and initial investment III. Her CARA utility at terminal wealth ω\omegaω is

uγ(ω)=γ−1−γ−1exp⁡(−γω),γ>0.u_\gamma(\omega)=\gamma^{-1}-\gamma^{-1}\exp(-\gamma\omega),\qquad\gamma>0.uγ​(ω)=γ−1−γ−1exp(−γω),γ>0.

Under a belief that the venture has type jjj, investors value the revenue using that type's distribution. An actual type iii who sells fraction α\alphaα therefore receives the cash value αEj[Rj]\alpha\mathbb E_j[R_j]αEj​[Rj​], retains fraction 1−α1-\alpha1−α of her own realized revenue, and has expected utility

Uij(α)=Ei[uγ ⁣(αEj[Rj]+(1−α)Ri−C−I)].U_{ij}(\alpha)=\mathbb E_i\left[u_\gamma\!\left(\alpha\mathbb E_j[R_j]+(1-\alpha)R_i-C-I\right)\right].Uij​(α)=Ei​[uγ​(αEj​[Rj​]+(1−α)Ri​−C−I)].

The financing method changes the smallest fraction she must sell. The equity financing threshold is αˉjE=I/(Ej[Rj]−C)\bar\alpha_j^E=I/(\mathbb E_j[R_j]-C)αˉjE​=I/(Ej​[Rj​]−C), and the token financing threshold is αˉjT=I/Ej[Rj]\bar\alpha_j^T=I/\mathbb E_j[R_j]αˉjT​=I/Ej​[Rj​]. The cost CCC appears in the equity denominator because equity buyers receive a claim after that cost. With a positive cost, the token threshold is strictly smaller.

Investors use a threshold belief: they infer HHH when α≤t\alpha\le tα≤t and LLL when α>t\alpha>tα>t. A separating equilibrium requires both types to prefer the side of the threshold assigned to them, subject to their financing constraints. The paper's conditions (23) and (24) compare the best feasible payoffs on the two sides. A least-cost separating equilibrium is a separating equilibrium that is not Pareto-dominated by another such equilibrium. The paper defines an interior share α∗\alpha^*α∗ in equation (25); Proposition 5 identifies its role as the least-cost threshold whenever separation is feasible.

Formalization targets

The central question is when tokens have a signaling advantage: a separating equilibrium exists with tokens, and none exists with equity. Proposition 6 gives two parameter windows. For fixed values of the other parameters, there are thresholds μl<μu\mu_l<\mu_uμl​<μu​ and σl<σu\sigma_l<\sigma_uσl​<σu​ such that

token signaling advantage⟺μl≤μL<μu,\text{token signaling advantage}\quad\Longleftrightarrow\quad\mu_l\le\mu_L<\mu_u,token signaling advantage⟺μl​≤μL​<μu​,

and, in a separate comparison,

token signaling advantage⟺σl≤σ<σu.\text{token signaling advantage}\quad\Longleftrightarrow\quad\sigma_l\le\sigma<\sigma_u.token signaling advantage⟺σl​≤σ<σu​.

The first window varies the low-type mean while keeping the high-type mean fixed. Increasing μL\mu_LμL​ reduces the difference between the types. The second varies payoff volatility. Both claims concern the game's equilibrium conditions, rather than only a numerical inequality between closed-form expressions.

The milestone sequence begins with the CARA-normal expected-utility formula (38) and the interior range of α∗\alpha^*α∗. It then records the low type's indifference at α∗\alpha^*α∗ and the threshold form of condition (24). Proposition 5 characterizes the equilibrium under each security. The remaining milestones compare the financing thresholds and record the monotonicity of α∗\alpha^*α∗ with respect to μL\mu_LμL​ and σ\sigmaσ. These are the paper's stated results and intermediate claims in §4 and its appendix.

What the result establishes

The theorem isolates a case in which the form of the security changes what investors can infer. A token offering can permit the high type to sell enough to finance the investment while retaining a stake that the low type does not want to imitate. Equity may require a larger sale fraction, closing that opportunity. The window statement also says the advantage is limited: when types are too different or too similar, or when volatility is outside the relevant range, tokens do not uniquely enable separation under this model.

A machine-checked development would connect an explicit probabilistic payoff model, incentive comparisons over feasible actions, the equilibrium characterization, and comparative statics in one chain. The cited paper proves these results informally; the theorem files here are open proof targets. The CARA-normal calculation and the representation of threshold equilibria are reusable in other signaling models with normally distributed payoffs.

Where the mathematical work lies

The direct comparison of financing thresholds is insufficient for the goal. Separation also depends on whether a low-quality entrepreneur would imitate a high-quality one, and whether a high-quality entrepreneur prefers to reveal her type. These conditions optimize utility over different feasible sale fractions, including a strict boundary on the low-belief side of the threshold. The paper's closed form for α∗\alpha^*α∗ expresses the point at which the low type becomes indifferent, but the full result must also establish the high type's incentive and the uniqueness of the refined equilibrium.

The comparative statics require control of two moving parts: the threshold α∗\alpha^*α∗ changes with the low-type mean and volatility, while the high-type financing thresholds remain fixed in those comparisons. The interval endpoints in Proposition 6 are existential, so their placement must reflect the feasible parameter domain and the strict or weak inequalities at each boundary.

Formalization scope

Lean represents the normal approximation for Xi=ξi2X_i=\xi_i^2Xi​=ξi2​ directly with gaussianReal of mean μi\mu_iμi​ and variance σ2\sigma^2σ2. This follows the paper's approximation even though a literal square is nonnegative and a nondegenerate normal variable is not. Revenue is Xi/4X_i/4Xi​/4. Expected utility remains a Gaussian integral; equation (38) is a theorem target, not a definition of that utility. All payoffs and thresholds are real numbers.

The parameter domain has γ>0\gamma>0γ>0, σ>0\sigma>0σ>0, μL<μH\mu_L<\mu_HμL​<μH​, C>0C>0C>0, I>0I>0I>0, and I≤μL/4−CI\le\mu_L/4-CI≤μL​/4−C. The last two substantive additions to the page's stated model are C>0C>0C>0, needed for the strict token-equity threshold comparison, and the low-type financing condition, needed for the paper's action αL=1\alpha_L=1αL​=1 to be feasible. An equilibrium threshold satisfies αˉHi≤t<1\bar\alpha_H^i\le t<1αˉHi​≤t<1. The least-cost refinement is represented by the paper's operational exclusion of Pareto-dominated separating equilibria, without claiming a general formalization of the Cho–Kreps intuitive criterion.

The high and low incentive conditions quantify only over feasible fractions. Equilibrium payoffs use real suprema on the corresponding action sets; within an equilibrium these sets are nonempty and the utility is bounded above. The goal defines signaling advantage through existence and nonexistence of separating equilibria, so it cannot be discharged by redefining the advantage as the closed-form inequality. Useful contributions include Gaussian expected-utility lemmas, optimization over the feasible intervals, equilibrium existence and uniqueness, and the real-analysis results needed for the parameter windows.

Selected references

  • Jacob Chod and Evgeny Lyandres, A Theory of ICOs: Diversification, Agency, and Information Asymmetry, SSRN working paper 3159528, version of December 14, 2019, §4 and Appendix; published in Management Science. Working paper. Published article DOI.
10 thms1 active userReviewed
Algorithmic Game TheoryControl TheoryProbability+1·Captain: mikedeng1

Dynkin Games with Incomplete and Asymmetric Information: With Linear Payoffs, a Threshold Stop and a Belief-Reflecting Randomised Stop Form a Saddle PointResearch Paper

Why this game matters

An optimal stopping game asks two players to choose when to stop a shared process. The first player receives one payoff if they stop first and another if the second player stops first. Such a game changes when one player knows a parameter that the other can only infer from observations. The informed player can reveal information through the timing of a stop, and can conceal it by randomising. De Angelis, Ekström and Glover study this interaction for a diffusion with one of two possible drifts. Their preprint gives a verification theorem for general payoffs and solves a linear geometric Brownian example explicitly. This mission targets the latter result, Theorem 6.10. The result is proved in the paper; the work here is to express its model, strategy classes, supporting results and conclusion in Lean.

The game and its information

The asset value starts at x>0x>0x>0 and evolves as a geometric Brownian motion Xt=xexp⁡((μθ−σ2/2)t+σWt)X_t=x\exp((\mu_\theta-\sigma^2/2)t+\sigma W_t)Xt​=xexp((μθ​−σ2/2)t+σWt​), where σ>0\sigma>0σ>0, μ0<0<μ1\mu_0<0<\mu_1μ0​<0<μ1​, and the hidden type θ∈{0,1}\theta\in\{0,1\}θ∈{0,1} selects the drift. Initially P(θ=1)=π∈(0,1)\mathbb P(\theta=1)=\pi\in(0,1)P(θ=1)=π∈(0,1); write ϕ=π/(1−π)\phi=\pi/(1-\pi)ϕ=π/(1−π) for the initial odds. Player 1 sees the history of XXX and chooses an almost surely finite stopping time τ\tauτ for its augmented observation filtration FX\mathcal F^XFX. Player 2 also knows θ\thetaθ and chooses an adapted, nondecreasing control Γ\GammaΓ for the augmented filtration FX,θ\mathcal F^{X,\theta}FX,θ. An independent uniform variable UUU turns this control into the randomised stopping time γ=inf⁡{t≥0:Γt>U}\gamma=\inf\{t\ge0:\Gamma_t>U\}γ=inf{t≥0:Γt​>U}. The initial value Γ0\Gamma_0Γ0​ may be positive, representing a jump at time zero.

If τ<γ\tau<\gammaτ<γ, Player 1 receives f(Xτ)f(X_\tau)f(Xτ​); if γ≤τ\gamma\le\tauγ≤τ, Player 1 receives g(Xγ)g(X_\gamma)g(Xγ​). Thus a finite tie pays ggg. The expected payoff is J(τ,Γ)J(\tau,\Gamma)J(τ,Γ). Player 1 maximises it and Player 2 minimises it. A saddle point (τ∗,Γ∗)(\tau^*,\Gamma^*)(τ∗,Γ∗) satisfies J(τ,Γ∗)≤J(τ∗,Γ∗)≤J(τ∗,Γ)J(\tau,\Gamma^*)\le J(\tau^*,\Gamma^*)\le J(\tau^*,\Gamma)J(τ,Γ∗)≤J(τ∗,Γ∗)≤J(τ∗,Γ) for every admissible pair (τ,Γ)(\tau,\Gamma)(τ,Γ). In the linear example f(y)=yf(y)=yf(y)=y and g(y)=(1+ε)yg(y)=(1+\varepsilon)yg(y)=(1+ε)y, where ε>0\varepsilon>0ε>0.

The uninformed player's posterior probability Πt=P(θ=1∣FtX)\Pi_t=\mathbb P(\theta=1\mid\mathcal F^X_t)Πt​=P(θ=1∣FtX​) is expressed through the posterior odds Φt=Πt/(1−Πt)\Phi_t=\Pi_t/(1-\Pi_t)Φt​=Πt​/(1−Πt​). The informed strategy is a pair (Γ0,Γ1)(\Gamma^0,\Gamma^1)(Γ0,Γ1), one control for each type. Conditional laws Pi=P(⋅∣θ=i)P^i=P(\cdot\mid\theta=i)Pi=P(⋅∣θ=i) define the interim payoffs J0,J1J^0,J^1J0,J1. The scaled ex-ante payoff is J^=(1+ϕ)J\widehat J=(1+\phi)JJ=(1+ϕ)J. These quantities permit both an ex-ante saddle-point statement and the agent-form equilibrium statement of Proposition 3.6.

Formalization targets

The supporting targets begin with the filtering and payoff identities of Section 3, the quasi-variational verification result of Theorem 5.1, and the reflected-likelihood and boundary-value facts of Section 6. They culminate in the explicit linear-game equilibrium. Two real numbers 0<A<B0<A<B0<A<B are determined by equations (90) and (93); piecewise power functions V1,VV_1,VV1​,V are defined by (89) and (91), extended by (94)–(97), and V0(ϕ)=V(ϕ)−ϕV1(ϕ)V_0(\phi)=V(\phi)-\phi V_1(\phi)V0​(ϕ)=V(ϕ)−ϕV1​(ϕ). The type-one control ΓB\Gamma^BΓB keeps the adjusted odds ΦB=Φ(1−ΓB)\Phi^B=\Phi(1-\Gamma^B)ΦB=Φ(1−ΓB) below BBB, while Player 1 stops at τA=inf⁡{t:ΦtB≤A}\tau_A=\inf\{t:\Phi^B_t\le A\}τA​=inf{t:ΦtB​≤A}.

The goal is Theorem 6.10: (τA,(0,ΓB))(\tau_A,(0,\Gamma^B))(τA​,(0,ΓB)) is a saddle point, and for every x,ϕ>0x,\phi>0x,ϕ>0,

J^x,ϕ(τA,(0,ΓB))=xV(ϕ),Jx,ϕ0(τA,0)=xV0(ϕ),Jx,ϕ1(τA,ΓB)=xV1(ϕ).\widehat J_{x,\phi}(\tau_A,(0,\Gamma^B))=xV(\phi),\qquad J^0_{x,\phi}(\tau_A,0)=xV_0(\phi),\qquad J^1_{x,\phi}(\tau_A,\Gamma^B)=xV_1(\phi).Jx,ϕ​(τA​,(0,ΓB))=xV(ϕ),Jx,ϕ0​(τA​,0)=xV0​(ϕ),Jx,ϕ1​(τA​,ΓB)=xV1​(ϕ).

The mission's milestone list follows the paper's stated results. It includes the payoff reductions, Theorem 5.1, reflection at BBB, the ODE bounds, the unique ratio A/BA/BA/B, and the verification of the explicit value functions.

What the result establishes

The theorem identifies actual strategies for both players, rather than only showing that the game's upper and lower values coincide. It specifies how the informed player's randomisation depends on the unobserved type and how the uninformed player responds to the adjusted odds. It also gives the ex-ante and both type-conditional payoffs. Without the saddle inequalities over the full strategy classes, the formulas would describe candidate strategies rather than an equilibrium.

A Lean development must capture the paper's information restrictions and time-zero behaviour before the analytical proof can be meaningful. Reusable components include a model-independent two-filtration randomised Dynkin frame, the augmented observation filtrations, randomised stopping times, nonnegative extended payoffs, conditional-law decomposition, and Stieltjes integration with an initial jump. The paper has a written proof; the statements here are open Lean goals, not machine-checked proofs.

Main difficulty

The verification theorem couples a partial differential inequality to a singular control and a stopping time. The adjusted odds can jump when the control jumps, so ordinary continuous-path differential conditions do not by themselves describe the admissible action. The support conditions and transversality limits in Theorem 5.1 are needed to connect the analytic candidate values to payoffs for strategies that can stop at any time. In the linear example, obtaining the power-function candidates is only one part of the target: their free boundaries must be valid, the reflected control must be admissible, and the saddle inequalities must hold against all admissible alternatives.

Formalization scope and conventions

The Lean model uses the explicit strong solution of the §6 GBM on a complete probability space carrying a Brownian motion and independent uniform devices U,ΘU,\ThetaU,Θ. It sets θπ=1{Θ<π}\theta^\pi=\mathbf1_{\{\Theta<\pi\}}θπ=1{Θ<π}​. Paths of the chosen Brownian modification are continuous. Player 1's filtration is generated by XXX plus null sets; Player 2's also contains θ\thetaθ; neither contains UUU. Stopping times live in [0,∞][0,\infty][0,∞], with Player 1 required to stop almost surely. Payoffs are nonnegative extended integrals and equal zero if both stopping times are infinite. The Stieltjes convention counts Γ0\Gamma_0Γ0​.

The posterior odds Φ\PhiΦ use the GBM closed form. The running-maximum formula defines ΓB\Gamma^BΓB and permits its initial jump. The tilted law used in Lemma 6.5 is represented by an equivalent likelihood-path law. Theorem 5.1 is stated for the GBM instance with general continuous 0≤f≤g0\le f\le g0≤f≤g and relative closures inside (0,∞)2(0,\infty)^2(0,∞)2. Its localisation integral retains the factor printed in the preprint. The explicit boundaries, controls and values in the goal are definitions, never free variables or hypotheses; the saddle inequalities quantify over the full player strategy classes.

There is a normalization slip in the paper's (23): it prints Φt=dP1/dP0\Phi_t=dP^1/dP^0Φt​=dP1/dP0 on FtX\mathcal F^X_tFtX​, although Φ0=ϕ\Phi_0=\phiΦ0​=ϕ and a Radon–Nikodym derivative at time zero equals one. Equations (20)–(21) give dP1/dP0=Φt/ϕdP^1/dP^0=\Phi_t/\phidP1/dP0=Φt​/ϕ. The formal filtering item uses that corrected identity. Equation (50) similarly identifies the ex-ante value with the scaled payoff; the formal target records J^=(1+ϕ)J\widehat J=(1+\phi)JJ=(1+ϕ)J and states the payoff identities without confusing the two. The source's domain R+2\mathbb R_+^2R+2​ in Theorem 6.10 is read as x,ϕ>0x,\phi>0x,ϕ>0, as required by the model. These readings are documented in the local moderation notes.

Selected references

  • T. De Angelis, E. Ekström and K. Glover, Dynkin Games with Incomplete and Asymmetric Information, arXiv preprint 1810.07674v5, 2020. Preprint.
19 thms1 active userReviewed
EconomicsOptimization·Captain: mikedeng1

Non-Ownership Business Models for Solar Energy: When a Solar Power Company Optimally Offers Power Purchase Agreements, Leasing and SalesResearch Paper

Why the business-model choice matters

Residential solar panels can reach customers through a purchase, a lease, or a power purchase agreement (PPA), under which the customer pays for the electricity the panels generate. These arrangements allocate panel ownership, subsidy receipts, and generation risk differently. A company choosing prices for all three therefore also chooses which customers adopt each arrangement. Agrawal, Toktay and Yücel study that choice in a model with net-metered electricity, investment and generation subsidies, and customers who differ in expected solar generation. Their main structural result identifies the company's optimal mix as panel cost and the investment subsidy change. The subsequent results about subsidy effects, renewable adoption, and restrictions on leasing use this mix as their starting point. Agrawal, Toktay and Yücel (2022).

Customers, prices and profit

One solar power company sets a normalized monthly sales price sss, a monthly lease fee lll, and a PPA price ppp per unit of generation. A customer can choose sale SSS, lease LLL, PPA PPP, or continued electricity supply from the utility UUU. The customer has expected monthly solar generation η\etaη and expected electricity demand ddd. The retail electricity price is rrr; α\alphaα is the fraction of panel investment cost covered by an investment subsidy, β\betaβ is the generation subsidy per unit, ttt is the company's PPA monitoring cost per unit, and ccc is its normalized panel cost. Under net metering, the customer's utility bill is r(q−n)r(q-n)r(q−n) when usage is qqq and generation is nnn. The paper's expected monthly costs of the four choices are

CS=(1−α)s+r(d−η)−βη,CL=l+r(d−η),CP=pη+r(d−η),CU=rd.C_S=(1-\alpha)s+r(d-\eta)-\beta\eta,\quad C_L=l+r(d-\eta),\quad C_P=p\eta+r(d-\eta),\quad C_U=rd.CS​=(1−α)s+r(d−η)−βη,CL​=l+r(d−η),CP​=pη+r(d−η),CU​=rd.

The customer chooses a least-cost option. Expected generation is uniform on [0,g][0,g][0,g] and independent of expected demand, whose distribution has positive support. Let Ni(s,l,p)N^i(s,l,p)Ni(s,l,p) be the population share choosing option iii and Gi(s,l,p)G^i(s,l,p)Gi(s,l,p) the total expected generation among those customers. An option is offered when its share is positive. The company's expected monthly profit, equation (1) of the source, is

Π(s,l,p)=(s−c)NS+(l−(1−α)c)NL+βGL+(p+β−t)GP−(1−α)cNP.\Pi(s,l,p)=(s-c)N^S+\bigl(l-(1-\alpha)c\bigr)N^L+\beta G^L +(p+\beta-t)G^P-(1-\alpha)cN^P.Π(s,l,p)=(s−c)NS+(l−(1−α)c)NL+βGL+(p+β−t)GP−(1−α)cNP.

These amounts come from actual customer choices, so their adoption shares are nonnegative by definition. The sales-only benchmark, in which customers choose between SSS and UUU, yields the panel-cost bound C0=g(r+β)/(1−α)C_0=g(r+\beta)/(1-\alpha)C0​=g(r+β)/(1−α). The source restricts its main analysis to c<C0c<C_0c<C0​. Agrawal, Toktay and Yücel (2022), §§3–4 and Online Appendix A.

Formalization targets

The goal is Proposition 1, the global choice of business models. Write A2=β/(β+r)A_2=\beta/(\beta+r)A2​=β/(β+r). When α<A2\alpha<A_2α<A2​, there are panel-cost thresholds C1≤C2C_1\le C_2C1​≤C2​ such that every profit-maximizing price triple has the offer set

{P,S}(c<C1),{P,L,S}(C1<c<C2),{L,S}(C2<c<C0).\{P,S\}\quad(c<C_1),\qquad \{P,L,S\}\quad(C_1<c<C_2),\qquad \{L,S\}\quad(C_2<c<C_0).{P,S}(c<C1​),{P,L,S}(C1​<c<C2​),{L,S}(C2​<c<C0​).

When A2≤αA_2\le\alphaA2​≤α, a threshold C3C_3C3​ separates {P,S}\{P,S\}{P,S} for c<C3c<C_3c<C3​ from {S}\{S\}{S} for C3<c<C0C_3<c<C_0C3​<c<C0​. The same thresholds apply as ccc and the admissible demand distribution vary, and a maximizing price triple exists for each admissible ccc. A zero-width middle regime is allowed when C1=C2C_1=C_2C1​=C2​. The equality points at C1,C2,C3C_1,C_2,C_3C1​,C2​,C3​ are excluded because the appendix's candidate cases become degenerate or have tied global optima there, so the printed assignment at those points does not hold for every maximizer. Agrawal, Toktay and Yücel (2022), Proposition 1 and proof.

The milestones establish the common optimal electricity usage of Lemma 1, customer segmentation by η\etaη in Lemma 2, and the sales-only profitability limit of Lemma A.1. They then identify the best profits in the regions p≤rp\le rp≤r and p>rp>rp>r and compare those profits across the investment-subsidy regimes. The mission follows the author's source indices and case names, including the appendix's closed-form candidate prices and values. Agrawal, Toktay and Yücel (2022), pp. 11, 14 and OA1–OA8.

What the result enables

The proposition states when a company can profitably maintain PPA, leasing, and sales side by side, and when it optimally drops one or both non-ownership models. This matters for the paper's later claims: a change in a subsidy or a legal restriction may alter not only prices within one model but which models serve customers at all. The paper's comparisons of total adoption, generation, and utility-firm profit therefore depend on the equilibrium classification. Agrawal, Toktay and Yücel (2022), Propositions 2–6.

The mathematical result is already proved in the source. This mission supplies compiled Lean statements for its model and central claims; the theorem proofs remain open. Completing them would provide a machine-checked route from individual cost comparisons and population integrals to the global pricing decision. The customer-choice and adoption definitions can also support formal work on the paper's later subsidy and restriction results. No machine-checked proof of these draft theorems is claimed here.

Main difficulty

Comparing customer costs at a fixed price triple gives generation intervals, but it does not by itself identify the company's best prices. The company optimizes three prices while the adopting population changes at interval boundaries. In the appendix, the resulting problem has distinct candidates when p≤rp\le rp≤r and p>rp>rp>r, and its preferred candidate changes with both ccc and α\alphaα. A direct comparison of two displayed candidate profits is therefore insufficient unless the candidate is optimal throughout its stated region and all boundary cases are accounted for. At a cost where two regional values agree, different maximizing price triples can have different offer sets. Agrawal, Toktay and Yücel (2022), proof of Proposition 1, OA4–OA8.

Formalization scope

Lean represents prices by three real coordinates and parameters by one structure. The standing conditions are 0<α0<\alpha0<α, 0<β0<\beta0<β, 0<t<r0<t<r0<t<r, β<r\beta<rβ<r, 0<g0<g0<g, and α<A\alpha<Aα<A, where A=β(r+t)/(r(r+β)+βt)A=\sqrt{\beta(r+t)/(r(r+\beta)+\beta t)}A=β(r+t)/(r(r+β)+βt)​ is the paper's concavity threshold. Lemma 2 and Lemma A.1, which the paper states before introducing β<r\beta<rβ<r and α<A\alpha<Aα<A on p. 15, assume only the §3 conditions (and α<1\alpha<1α<1 for Lemma A.1). The goal also assumes 0<c<C00<c<C_00<c<C0​. Its demand distribution is an arbitrary probability measure concentrated on (0,∞)(0,\infty)(0,∞), which generalizes the paper's density h(d)h(d)h(d) without changing the claims, because the choice comparisons do not depend on ddd. Generation uses the uniform probability law on [0,g][0,g][0,g] and integrates η\etaη over the actual choice sets.

The four realized and expected costs are defined separately. Ties among S,L,PS,L,PS,L,P use a fixed priority order; UUU is selected only if strictly cheaper than all three solar options. This implements the appendix's allocation to PPA at p=rp=rp=r, where PPA and utility costs tie on a set of positive mass. The theorem ranges over all real price triples. Its offer set is calculated from positive adoption at a maximizer, and it asserts existence of that maximizer; neither the offer regime nor optimality is assigned by a case label. The cost thresholds are chosen before quantifying over the demand distribution and ccc. A milestone using η2=l/p\eta_2=l/pη2​=l/p requires p>0p>0p>0, and the closed-form denominators are used only under the stated parameter conditions. Population measures are finite, and the choice regions are measurable; these facts matter when converting adoption mass to a real number and integrating generation. The sales-only benchmark and the appendix formulas are separate definitions, not substitutes for equation (1).

Contributions are welcome on customer-choice measurability, integral identities for adoption and generation, the regional maximization claims, and the algebraic value comparisons. Propositions 2–6 and the Technical Appendix's alternative policy models are outside this mission.

Selected references

  • V. V. Agrawal, L. B. Toktay and Ş. Yücel, Non-Ownership Business Models for Solar Energy, Manufacturing & Service Operations Management, 2022; authors' manuscript dated 18 September 2021 used for the page and theorem citations here. DOI: 10.1287/msom.2021.1049.
12 thms1 active userReviewed
Algorithmic Game TheoryMechanism Design·Captain: mikedeng1

Procurement Mechanisms with Post-Auction Pre-Award Cost-Reduction Investigations: Symmetric Suppliers Earn Weakly Higher Ex Ante Profits When Investigations Precede the AwardResearch Paper

Motivation

Large manufacturers procure components by competitive bidding, and many of them also run cost-reduction investigations: the buyer's engineers visit a supplier's plant and look for process changes that lower the cost of the contracted part. The savings are passed to the buyer as a lower contract price. In practice such investigations usually happen after the contract is awarded, on the winner only. Chen, Beil and Duenyas (Operations Research 70(6), 2022) compare that convention with a mechanism in which the buyer investigates before the award, possibly several suppliers, and then awards the contract on the investigated prices.

Investigating before the award obviously helps the buyer for fixed bids: she sees more cost reductions before choosing. The question is whether suppliers lose. The buyer keeps every identified saving, and investigations of competitors make their prices more aggressive, so suppliers seem to have reason to resist. The paper shows that, when suppliers are ex ante symmetric, they do not: every supplier's ex ante expected profit is at least as high when investigations come before the award. This mission formalizes that result.

Setting

A buyer procures one contract from suppliers N={1,…,N}\mathcal N=\{1,\dots,N\}N={1,…,N} and must transact with one of them. Supplier iii has a private cost ci=Δic+ϵicc_i=\Delta^c_i+\epsilon^c_ici​=Δic​+ϵic​, where Δic≥0\Delta^c_i\ge0Δic​≥0 is public and ϵic\epsilon^c_iϵic​ has density fff on [c‾,cˉ][\underline c,\bar c][c​,cˉ] with mean 000; so cic_ici​ has density fi(c)=f(c−Δic)f_i(c)=f(c-\Delta^c_i)fi​(c)=f(c−Δic​), c.d.f. FiF_iFi​, and support Ci=[c‾i,cˉi]C_i=[\underline c_i,\bar c_i]Ci​=[c​i​,cˉi​]. Costs are independent.

An investigation of supplier iii costs the buyer d≥0d\ge0d≥0 and reveals a cost reduction ti=Δit+ϵit≥0t_i=\Delta^t_i+\epsilon^t_i\ge0ti​=Δit​+ϵit​≥0, where ϵit\epsilon^t_iϵit​ has density ggg on [t‾,tˉ][\underline t,\bar t][t​,tˉ] with mean 000, independently of everything else. At most N^∈{1,…,N}\hat N\in\{1,\dots,N\}N^∈{1,…,N} suppliers can be investigated, and costs always exceed reductions: Δit+tˉ<Δic+c‾\Delta^t_i+\bar t<\Delta^c_i+\underline cΔit​+tˉ<Δic​+c​. Each supplier submits a sealed bid bi∈Bi=[0,cˉi]b_i\in B_i=[0,\bar c_i]bi​∈Bi​=[0,cˉi​]. If an investigated supplier wins, he is paid his updated bid bi−tib_i-t_ibi​−ti​ and his cost becomes ci−tic_i-t_ici​−ti​, so his margin is bi−cib_i-c_ibi​−ci​ in every case.

  • BAI (bid, award, investigate). The buyer awards the contract to the lowest λi=bi−ri\lambda_i=b_i-r_iλi​=bi​−ri​, where ri=max⁡{0,Δit−d}r_i=\max\{0,\Delta^t_i-d\}ri​=max{0,Δit​−d} is the expected saving from investigating the winner afterwards. Ties are broken at random.
  • BIA (bid, investigate, award). The buyer investigates sequentially. Each round she uses the bids and the reductions revealed so far to stop or to investigate one more supplier, and stopping is final. She then awards the contract to the lowest updated bid. Her policy minimizes the expected payment plus investigation cost, problem (1). Theorem 1 solves (1) by an index policy: with τ\tauτ the root of d=∫τ∞(ϵt−τ) dG(ϵt)d=\int_\tau^\infty(\epsilon^t-\tau)\,dG(\epsilon^t)d=∫τ∞​(ϵt−τ)dG(ϵt), each supplier gets the index μi(bi)=bi−max⁡{0,τ+Δit}\mu_i(b_i)=b_i-\max\{0,\tau+\Delta^t_i\}μi​(bi​)=bi​−max{0,τ+Δit​}. The buyer investigates in increasing order of index while the current lowest updated bid is at least the next index, up to a cap Nˉ≤N^\bar N\le\hat NNˉ≤N^.

Supplier iii's winning probability at bids bbb is Wi(b)W_i(b)Wi​(b), over ttt and the tie-break. With the others bidding β−i(c−i)\beta_{-i}(c_{-i})β−i​(c−i​), his expected profit from bidding bbb at cost ccc is ui(b;c)=(b−c) Ec−i[Wi(b,β−i(c−i))]u_i(b;c)=(b-c)\,\mathbf E_{c_{-i}}[W_i(b,\beta_{-i}(c_{-i}))]ui​(b;c)=(b−c)Ec−i​​[Wi​(b,β−i​(c−i​))], equation (5). A Bayesian Nash equilibrium is a profile of strategies βi:Ci→Bi\beta_i:C_i\to B_iβi​:Ci​→Bi​ that are best replies at every type. Supplier iii's ex ante expected profit is Eci[ui(βi(ci);ci)]\mathbf E_{c_i}[u_i(\beta_i(c_i);c_i)]Eci​​[ui​(βi​(ci​);ci​)].

Formalization targets

Goal: Theorem 2 (weak form)

Suppose that suppliers are ex ante symmetric (Δic=Δc\Delta^c_i=\Delta^cΔic​=Δc, Δit=Δt\Delta^t_i=\Delta^tΔit​=Δt), that Fi/fiF_i/f_iFi​/fi​ is nondecreasing (Assumption 1), and that d<max⁡iΔitd<\max_i\Delta^t_id<maxi​Δit​ (Assumption 2). Let βBIA\beta^{\mathrm{BIA}}βBIA be a symmetric equilibrium of BIA, and βBAI\beta^{\mathrm{BAI}}βBAI a symmetric nondecreasing equilibrium of BAI. Then for every supplier iii,

Eci[uiBAI(βiBAI(ci);ci)]  ≤  Eci[uiBIA(βiBIA(ci);ci)].\mathbf E_{c_i}\big[u^{\mathrm{BAI}}_i(\beta^{\mathrm{BAI}}_i(c_i);c_i)\big]\;\le\;\mathbf E_{c_i}\big[u^{\mathrm{BIA}}_i(\beta^{\mathrm{BIA}}_i(c_i);c_i)\big].Eci​​[uiBAI​(βiBAI​(ci​);ci​)]≤Eci​​[uiBIA​(βiBIA​(ci​);ci​)].

Milestones

  1. Proposition 1. There exists an equilibrium βˇ\check\betaβˇ​ of the first-price auction with costs ci−ric_i-r_ici​−ri​ whose shifted strategy β~i(ci)=βˇi(ci−ri)+ri\tilde\beta_i(c_i)=\check\beta_i(c_i-r_i)+r_iβ~​i​(ci​)=βˇ​i​(ci​−ri​)+ri​ is a BAI equilibrium.
  2. Theorem 1(a). The index policy (3) minimizes the buyer's expected cost (1) over all admissible sequential policies, for every bid vector.
  3. Zero profit at the top. UiM(cˉi)=0U^{\mathcal M}_i(\bar c_i)=0UiM​(cˉi​)=0 for M∈{BIA,BAI}\mathcal M\in\{\mathrm{BIA},\mathrm{BAI}\}M∈{BIA,BAI}.
  4. Envelope formula (A.2). UiM(ci)=UiM(cˉi)+∫cicˉiϕiM(y) dyU^{\mathcal M}_i(c_i)=U^{\mathcal M}_i(\bar c_i)+\int_{c_i}^{\bar c_i}\phi^{\mathcal M}_i(y)\,dyUiM​(ci​)=UiM​(cˉi​)+∫ci​cˉi​​ϕiM​(y)dy.
  5. Ex ante identity. Eci[UiM(ci)]=UiM(cˉi)+Ec[AiM(c) Fi(ci)/fi(ci)]\mathbf E_{c_i}[U^{\mathcal M}_i(c_i)]=U^{\mathcal M}_i(\bar c_i)+\mathbf E_{c}\big[A^{\mathcal M}_i(c)\,F_i(c_i)/f_i(c_i)\big]Eci​​[UiM​(ci​)]=UiM​(cˉi​)+Ec​[AiM​(c)Fi​(ci​)/fi​(ci​)].
  6. Pointwise comparison. At cost vector ccc and realized reduction vector ttt, ∑iAiBAI(c)F1(ci)/f1(ci)=min⁡iF1(ci)/f1(ci)≤∑iAiBIA(c,t)F1(ci)/f1(ci)\sum_i A^{\mathrm{BAI}}_i(c) F_1(c_i)/f_1(c_i)=\min_i F_1(c_i)/f_1(c_i)\le\sum_i A^{\mathrm{BIA}}_i(c,t) F_1(c_i)/f_1(c_i)∑i​AiBAI​(c)F1​(ci​)/f1​(ci​)=mini​F1​(ci​)/f1​(ci​)≤∑i​AiBIA​(c,t)F1​(ci​)/f1​(ci​) almost everywhere, with ∑iAiM=1\sum_iA^{\mathcal M}_i=1∑i​AiM​=1.

Significance

The theorem identifies a win-win procurement design. With symmetric suppliers, a buyer who chooses the cheaper of BAI and BIA weakly lowers her own cost, and this choice never lowers any supplier's ex ante profit. The gain does not come from sharing savings, since the buyer keeps all identified savings. It comes from strategic bidding: under BIA a supplier need not submit the lowest initial bid to win, so price competition at the bid stage softens. The paper's Figure 1 shows the higher BIA bid curves. Theorem 1 is of independent interest: it is an index solution of a capped sequential search problem in which an uninvestigated bid may still win.

The paper proves the result on paper. To our knowledge none of it is machine-checked. A formalization adds three things. It makes the mechanisms precise, including the tie-breaking and the stopping rule that the text leaves partly implicit. It pins down the right inequality: the printed "higher" is not strict, and with a single supplier both mechanisms give the same profit. It also produces reusable components: an envelope-theorem layer for bidding games with an indirect mechanism, and a capped-search index theorem.

Difficulty

The proof of Theorem 2 needs the BIA winning probability only through its total, ∑iAiBIA=1\sum_i A^{\mathrm{BIA}}_i=1∑i​AiBIA​=1. Everything else rests on three technical steps. The first is the envelope formula (A.2) for an indirect mechanism whose allocation WiW_iWi​ has no closed form and need not be continuous in the bid. The usual argument assumes a differentiable direct mechanism and does not apply as stated. The second is the claim that in BAI the lowest-cost supplier wins almost surely. That needs a symmetric equilibrium that is not merely nondecreasing but free of pooling, and this must be derived from the equilibrium conditions. The third is the "by symmetry" step that equates each supplier's ex ante profit with the average over suppliers. It needs both the equilibria and the mechanisms to be invariant under relabelling the suppliers.

Theorem 1(a) is a dynamic program over adaptive policies with a cap and an outside option. A Weitzman-style argument must handle a cap Nˉ\bar NNˉ and suppliers whose index makes investigation worthless (τi=0\tau_i=0τi​=0). It cannot simply be imported from Pandora's-box results.

Formalization scope

Suppliers are Fin N with N≥1N\ge1N≥1; all values are real numbers. A single definitions module PreAwardInvest.WinWin.Model contains the following.

  • Parameters and assumptions. The parameters, and the standing assumptions of §3.1 (densities with mean zero, nonnegative costs and reductions, footnote 6, 1≤N^≤N1\le\hat N\le N1≤N^≤N, d≥0d\ge0d≥0), with f>0f>0f>0 on (c‾,cˉ)(\underline c,\bar c)(c​,cˉ) added. Assumptions 1 and 2 are separate predicates.
  • Laws. The cost and reduction laws are densities with respect to Lebesgue measure, and their vectors have product laws.
  • Policies. A policy reads the bids, the investigated set and the revealed reductions. An admissible policy is nonanticipative and measurable in the reductions, and stopping is final.
  • Index policy. It takes τ\tauτ as any root of its equation. Ties in the ranking and in both awards are broken by one uniformly random priority order, so that BIA and BAI are symmetric under relabelling. Theorem 1(a) compares policies for bids in Bi=[0,cˉi]B_i=[0,\bar c_i]Bi​=[0,cˉi​].
  • Profits and equilibria. Profits follow (5), with the expectation over c−ic_{-i}c−i​ taken over the whole cost vector with coordinate iii overwritten. Equilibria require measurable strategies into Bi=[0,cˉi]B_i=[0,\bar c_i]Bi​=[0,cˉi​]. Ex ante profit is the expectation of the equilibrium profit, not the Myerson virtual-cost expression.

Several conventions are disclosed hypotheses or corrections. The goal uses ≤\le≤ in place of the printed "higher". The equilibria in the goal are assumed symmetric, and the BAI one nondecreasing (Lemma 1's class). Strict monotonicity is not assumed, because with one supplier the equilibrium bid is the constant cˉ1\bar c_1cˉ1​.

The goal is not trivialized. BAI awards on bids through λi=bi−ri\lambda_i=b_i-r_iλi​=bi​−ri​ and BIA runs the index policy on bids and revealed reductions. Neither mechanism has "lowest cost wins" built in, which is milestone 6. Ex ante profit is defined from (5), so milestones 4–5 are not definitional. A sorry-free check exhibits a one-supplier instance that satisfies every hypothesis of the goal. In it both ex ante profits equal 1/21/21/2.

Not formalized: Theorem 1(b)–(c), Lemma 1 (existence via Athey 2001), the closed form (4) of WiW_iWi​, the asymmetric results (Propositions 2–4), the buyer-cost comparison (Theorem 3), the optimal mechanism (Theorem 4), and the extensions of §9. Contributions welcome: the envelope theorem for monotone-selection bidding games, measurability of the winning probabilities, and the index policy's optimality.

Selected references

  • Q. Chen, D. R. Beil, I. Duenyas, Procurement mechanisms with post-auction pre-award cost-reduction investigations, Operations Research 70(6):3054–3075, 2022. https://doi.org/10.1287/opre.2022.2349 (accepted manuscript: LBS Research Online eprint 2544).
  • S. Athey, Single crossing properties and the existence of pure strategy equilibria in games of incomplete information, Econometrica 69(4):861–889, 2001. https://doi.org/10.1111/1468-0262.00223
  • R. B. Myerson, Optimal auction design, Mathematics of Operations Research 6(1):58–73, 1981. https://doi.org/10.1287/moor.6.1.58
  • P. Milgrom, I. Segal, Envelope theorems for arbitrary choice sets, Econometrica 70(2):583–601, 2002. https://doi.org/10.1111/1468-0262.00296
  • M. L. Weitzman, Optimal search for the best alternative, Econometrica 47(3):641–654, 1979. https://doi.org/10.2307/1910412
  • V. Krishna, Auction Theory, 2nd ed., Academic Press, 2009. https://doi.org/10.1016/C2009-0-22474-3
8 thms1 active userReviewed
AnalysisOptimization·Captain: mikedeng1

Personalized Robo-Advising: Enhancing Investment Through Client Interaction: The Approximate Personalization Measure Has a Unique Optimal Time Between Interactions, Increasing in Behavioral BiasResearch Paper

How often should a robo-advisor talk to its client?

A robo-advisor is an automated portfolio manager. It invests a client's wealth according to a model of the client's risk preferences, and it learns those preferences by interacting with the client: questionnaires, check-ins, requests to confirm or revise a risk profile. The design question of how often such a system should query its client is practical: each interaction updates the advisor's view of a preference that drifts over time, but each interaction is also exposed to the client's momentary reaction to recent market moves.

Capponi, Ólafsson and Zariphopoulou (Management Science 68(4), 2022; read here as arXiv:1911.01391v2) model this trade-off. Their Proposition 5.1 gives a closed-form approximation of how well the portfolio tracks the client's true risk aversion as a function of the time between interactions, and shows that this approximation has a unique optimal interaction time whose position responds monotonically to the model's parameters. This mission formalizes part (ii) of that proposition and the comparative statics (5.6).

Setting

Time is discrete. The client's true risk aversion γnC\gamma^C_nγnC​ evolves over the investment horizon; in each period it receives, with probability pϵ∈[0,1]p_\epsilon \in [0,1]pϵ​∈[0,1], an idiosyncratic shock σϵWn−σϵ2/2\sigma_\epsilon W_n - \sigma_\epsilon^2/2σϵ​Wn​−σϵ2​/2 with WnW_nWn​ standard Gaussian and σϵ>0\sigma_\epsilon > 0σϵ​>0, and no shock otherwise (p. 17). The robo-advisor and client interact every ϕ≥1\phi \ge 1ϕ≥1 periods. At an interaction the client communicates a risk-aversion parameter that is distorted by a behavioral bias: it is the true value multiplied by e−β⋅(average excess market return since the last interaction)e^{-\beta \cdot(\text{average excess market return since the last interaction})}e−β⋅(average excess market return since the last interaction), where the coefficient β≥0\beta \ge 0β≥0 measures the strength of the bias (eq. (5.2), pp. 17–18). Between interactions the advisor keeps using the last communicated value. Market returns have volatility σ(y)>0\sigma(y) > 0σ(y)>0 in economic state yyy, and σ0:=σ(Y0)\sigma_0 := \sigma(Y_0)σ0​:=σ(Y0​) is the volatility in the initial state.

The personalization measure (5.4) is the expected time-average relative error between the advisor's risk tolerance 1/γn1/\gamma_n1/γn​ and the client's 1/γnC1/\gamma^C_n1/γnC​,

R(ϕ,β)=E[1T∑n=0T−1∣1/γn−1/γnC1/γnC∣].\mathcal R(\phi,\beta) = \mathbb E\Big[\frac1T\sum_{n=0}^{T-1}\Big|\frac{1/\gamma_n - 1/\gamma^C_n}{1/\gamma^C_n}\Big|\Big].R(ϕ,β)=E[T1​n=0∑T−1​​1/γnC​1/γn​−1/γnC​​​].

Lower is better. Proposition 5.1(i) shows that R\mathcal RR is, up to error terms, a multiple of the approximate personalization measure

R~(ϕ,β)=2π(βσ0ϕ(1−ϕ−12pϵ)+β2σ02ϕ+σϵ2  ϕ−12pϵ).\widetilde{\mathcal R}(\phi,\beta) = \sqrt{\frac{2}{\pi}}\left(\frac{\beta\sigma_0}{\sqrt{\phi}}\Big(1-\frac{\phi-1}{2}p_\epsilon\Big) + \sqrt{\frac{\beta^2\sigma_0^2}{\phi}+\sigma_\epsilon^2}\;\frac{\phi-1}{2}p_\epsilon\right).R(ϕ,β)=π2​​​ϕ​βσ0​​(1−2ϕ−1​pϵ​)+ϕβ2σ02​​+σϵ2​​2ϕ−1​pϵ​​.

The first term is the bias cost, which shrinks as interactions become rarer (the bias averages returns over a longer window); the second is the staleness cost, which grows with the number of shocks missed between interactions. The mission is about R~\widetilde{\mathcal R}R, as Proposition 5.1(ii) is. In Lean the parameters are a structure Params with fields sig0 (σ0\sigma_0σ0​), pe (pϵp_\epsilonpϵ​), se (σϵ\sigma_\epsilonσϵ​), β, and R~\widetilde{\mathcal R}R is Rtilde P φ.

Formalization targets

Goal: Proposition 5.1(ii) and (5.6)

Assume σ0>0\sigma_0 > 0σ0​>0, σϵ>0\sigma_\epsilon > 0σϵ​>0, pϵ∈[0,1]p_\epsilon \in [0,1]pϵ​∈[0,1], β≥0\beta \ge 0β≥0 and β+pϵ>0\beta + p_\epsilon > 0β+pϵ​>0. Then

arg⁡min⁡ϕ≥1R~(ϕ,β)={1,β=0,ϕ0∈[1,∞),β>0, pϵ>0,∞,pϵ=0,\arg\min_{\phi\ge1}\widetilde{\mathcal R}(\phi,\beta) = \begin{cases} 1, & \beta = 0,\\ \phi_0\in[1,\infty), & \beta>0,\ p_\epsilon>0,\\ \infty, & p_\epsilon = 0,\end{cases}argϕ≥1min​R(ϕ,β)=⎩⎨⎧​1,ϕ0​∈[1,∞),∞,​β=0,β>0, pϵ​>0,pϵ​=0,​

and ϕ0\phi_0ϕ0​ is increasing in β\betaβ and σ0\sigma_0σ0​ and decreasing in pϵp_\epsilonpϵ​ and σϵ\sigma_\epsilonσϵ​ (5.6). "arg⁡min⁡=∞\arg\min = \inftyargmin=∞" means that no ϕ≥1\phi \ge 1ϕ≥1 is a minimizer, R~\widetilde{\mathcal R}R decreases strictly on [1,∞)[1,\infty)[1,∞) and tends to its infimum 000.

Milestones

The proof (Appendix C, p. 46) works with the bracket f=π/2 R~f = \sqrt{\pi/2}\,\widetilde{\mathcal R}f=π/2​R and an auxiliary function ggg with f′(ϕ)=12βσ0pϵϕ−2g(ϕ)f'(\phi) = \tfrac12\beta\sigma_0p_\epsilon\phi^{-2}g(\phi)f′(ϕ)=21​βσ0​pϵ​ϕ−2g(ϕ). The milestones are its displayed claims, in order: f≥0f \ge 0f≥0; the two degenerate cases pϵ=0p_\epsilon = 0pϵ​=0 and β=0\beta = 0β=0; the factorization of f′f'f′; the value of f′(1)f'(1)f′(1) and the criterion f′(1)≥0  ⟺  pϵσϵ2/(β2σ02)≥2+1/pϵf'(1) \ge 0 \iff p_\epsilon\sigma_\epsilon^2/(\beta^2\sigma_0^2) \ge 2 + 1/p_\epsilonf′(1)≥0⟺pϵ​σϵ2​/(β2σ02​)≥2+1/pϵ​; the sufficient condition pϵσϵ/(βσ0)<1⇒∂ϕR~∣ϕ=1<0p_\epsilon\sigma_\epsilon/(\beta\sigma_0) < 1 \Rightarrow \partial_\phi\widetilde{\mathcal R}|_{\phi=1} < 0pϵ​σϵ​/(βσ0​)<1⇒∂ϕ​R∣ϕ=1​<0 of p. 22; the lower bound g′(ϕ)≥12ϕ−1/2g(1)g'(\phi) \ge \tfrac12\phi^{-1/2}g(1)g′(ϕ)≥21​ϕ−1/2g(1), which settles the corner case g(1)≥0g(1)\ge 0g(1)≥0; and the interior case g(1)<0g(1) < 0g(1)<0, where ggg is strictly convex, tends to ∞\infty∞, and has a unique zero at which fff is minimized.

Significance

The result gives the model's two qualitative predictions in exact form. Without behavioral bias the advisor should interact at every period; without idiosyncratic shocks it should never interact; when both are present there is a single best interaction frequency, and it moves in the economically expected direction: more bias or more volatile markets call for rarer interaction, more frequent or larger preference shocks call for more frequent interaction. The p. 22 criterion gives a checkable condition, in terms of pϵσϵp_\epsilon\sigma_\epsilonpϵ​σϵ​ against βσ0\beta\sigma_0βσ0​, under which interacting at every period is not optimal.

The proof in the paper is a sketch. The convexity of ggg in the interior case is asserted ("it can be shown"), and the comparative statics (5.6) are stated as following from "the form of the function fff". As printed, (5.6) asserts strict signs of partial derivatives of ϕ0\phi_0ϕ0​, which fail on the open parameter region where the minimizer sits at the corner ϕ0=1\phi_0 = 1ϕ0​=1 (there ϕ0\phi_0ϕ0​ is locally constant). A machine-checked proof therefore supplies the missing convexity argument and the corrected comparative statics. No part of the proposition is known to be formalized.

Difficulty

The function R~\widetilde{\mathcal R}R mixes ϕ−1/2\phi^{-1/2}ϕ−1/2, ϕ1/2\phi^{1/2}ϕ1/2 and a square root of β2σ02/ϕ+σϵ2\beta^2\sigma_0^2/\phi + \sigma_\epsilon^2β2σ02​/ϕ+σϵ2​. The sign of its derivative is not visible from the formula, and the paper's sign analysis passes through ggg; the paper's claim that ggg is strictly convex has no written proof. The comparative statics are not a differentiable implicit-function computation either: the minimizer switches between the corner ϕ=1\phi = 1ϕ=1 and an interior root of ggg, and the monotonicity must hold across that switch. The obvious move of differentiating the first-order condition g(ϕ0)=0g(\phi_0) = 0g(ϕ0​)=0 in the parameters covers only the interior region and says nothing at the corner.

Formalization scope

All quantities are real. The time between interactions ϕ\phiϕ is a real variable on [1,∞)[1,\infty)[1,∞): the model uses integer ϕ\phiϕ, but Proposition 5.1 differentiates in ϕ\phiϕ and places ϕ0\phi_0ϕ0​ in [1,∞)[1,\infty)[1,∞). The standing ranges σ0>0\sigma_0 > 0σ0​>0, σϵ>0\sigma_\epsilon > 0σϵ​>0, pϵ∈[0,1]p_\epsilon \in [0,1]pϵ​∈[0,1], β≥0\beta \ge 0β≥0 are bundled in Params.Standing. "ϕ\phiϕ minimizes R~\widetilde{\mathcal R}R" is the predicate IsOptInteraction P φ: ϕ≥1\phi \ge 1ϕ≥1 and R~(ϕ)≤R~(ψ)\widetilde{\mathcal R}(\phi) \le \widetilde{\mathcal R}(\psi)R(ϕ)≤R(ψ) for all ψ≥1\psi \ge 1ψ≥1; there is no argmin function, so no junk value can satisfy the goal. Square roots are Real.sqrt, ϕ3/2\phi^{3/2}ϕ3/2 is ϕϕ\phi\sqrt\phiϕϕ​, and every statement restricts ϕ\phiϕ so that no radicand is negative and no denominator vanishes. Derivatives are stated with HasDerivAt and an explicit value. Rtilde P φ takes ϕ\phiϕ last, while the page writes both R~(ϕ,β)\widetilde{\mathcal R}(\phi,\beta)R(ϕ,β) and R~(β,ϕ)\widetilde{\mathcal R}(\beta,\phi)R(β,ϕ).

Two readings are fixed. The case "arg⁡min⁡=∞\arg\min = \inftyargmin=∞" is: no minimizer exists, R~\widetilde{\mathcal R}R is strictly decreasing on [1,∞)[1,\infty)[1,∞), and R~(ϕ)→0\widetilde{\mathcal R}(\phi) \to 0R(ϕ)→0. The comparative statics (5.6) are stated as monotone comparative statics: if a parameter increases (keeping pϵ≤1p_\epsilon \le 1pϵ​≤1), the minimizer moves weakly in the stated direction, and strictly whenever the original minimizer exceeds 111. The monotonicity clauses quantify over minimizers for both parameter sets, and the existence clause of the goal ensures both exist, so they are not vacuous. Stating the goal about fff alone, or defining ϕ0\phi_0ϕ0​ by sInf or by a formula, would trivialize the claim; the goal is about R~\widetilde{\mathcal R}R and its minimizers.

Part (i) of Proposition 5.1, which relates R\mathcal RR to R~\widetilde{\mathcal R}R with unquantified O(⋅)O(\cdot)O(⋅) error terms, is not formalized; neither is the optimal investment strategy of Proposition 4.1. The development needs one-variable real calculus from Mathlib, including square-root and quotient derivatives, monotonicity and convexity on intervals, and limits at infinity. Proofs of the milestones, of the convexity of ggg on all of [1,∞)[1,\infty)[1,∞) (which appears to hold without the hypothesis g(1)<0g(1) < 0g(1)<0), and of the derivative form of (5.6) on the interior region are all welcome.

Selected references

  • A. Capponi, S. Ólafsson, T. Zariphopoulou, Personalized Robo-Advising: Enhancing Investment through Client Interaction, Management Science 68(4), 2022; preprint arXiv:1911.01391v2 (the version cited for all page and equation numbers). https://arxiv.org/abs/1911.01391
10 thms1 active userReviewed
EconomicsOptimization·Captain: mikedeng1

Labor Welfare in On-Demand Service Platforms: With Delay-Sensitive Customers, Average Labor Welfare First Increases and Then Decreases in the Labor Pool SizeResearch Paper

Motivation

On-demand service platforms (ride-hailing, delivery, home services) do not employ their workers. They set a price for customers and a wage per unit of work, and independent agents decide how much time to make themselves available. A recurring policy question, raised by regulators and labor advocates alike, is whether platforms should be allowed to recruit labor pools of unlimited size: a larger pool is often argued to depress what each worker earns. Benjaafar, Ding, Kong and Taylor (SSRN 3103447; published in Manufacturing & Service Operations Management, doi:10.1287/msom.2020.0964) study this question in a model where customers care about delay. Their central result, Theorem 2, is that the answer depends on delay sensitivity: in a standard market a larger pool lowers average labor welfare, but when customers are delay-sensitive, average labor welfare first rises and then falls in the pool size.

Setting

The model has two populations, both treated as continua.

Agents. A labor pool of mass N>0N>0N>0. An agent with opportunity-cost parameter aaa who is available for time xxx incurs cost ax2a x^2ax2; aaa has distribution FFF on [a‾,aˉ][\underline a,\bar a][a​,aˉ] with 0<a‾0<\underline a0<a​. Given the wage www, demand DDD and labor supply SSS, work is shared in proportion to availability, so the agent's effective wage is we=wD/Sw_e=wD/Swe​=wD/S and his best response is x∗(a)=wD/(2aS)x^*(a)=wD/(2aS)x∗(a)=wD/(2aS). Aggregation gives the labor supply equation

S=wDN E[1/a]/2.S=\sqrt{wDN\,E[1/a]/2}.S=wDNE[1/a]/2​.

Customers. Needs arise at rate Λ\LambdaΛ (the market size); a customer with valuation v∼Unif[0,vˉ]v\sim\mathrm{Unif}[0,\bar v]v∼Unif[0,vˉ] seeks service if v−p−hW≥0v-p-hW\ge 0v−p−hW≥0, where h≥0h\ge0h≥0 is the delay sensitivity and

W=(D/S)βS−D,β≥1,W=\frac{(D/S)^\beta}{S-D},\qquad \beta\ge1,W=S−D(D/S)β​,β≥1,

is the delay. Demand is D=Λ[1−G(p+hW)]D=\Lambda[1-G(p+hW)]D=Λ[1−G(p+hW)] with GGG the uniform cdf, and an equilibrium requires D<SD<SD<S. In the benchmark h=0h=0h=0 demand is Λ[1−G(p)]\Lambda[1-G(p)]Λ[1−G(p)] and the capacity constraint is D≤SD\le SD≤S.

Platform. The platform chooses p≥0p\ge0p≥0, w≥0w\ge0w≥0 to maximize (p−w)D(p-w)D(p−w)D subject to these equilibrium conditions; the optimum is (p∗,w∗,D∗,S∗)(p^*,w^*,D^*,S^*)(p∗,w∗,D∗,S∗). The agent utilization is ρ∗=D∗/S∗\rho^*=D^*/S^*ρ∗=D∗/S∗, and average labor welfare is

LW‾=∫[w∗D∗S∗x∗(a)−a x∗(a)2]dF(a).\overline{LW}=\int\Big[\frac{w^*D^*}{S^*}x^*(a)-a\,x^*(a)^2\Big]dF(a).LW=∫[S∗w∗D∗​x∗(a)−ax∗(a)2]dF(a).

The pool sizes for which a positive optimum exists are N>N‾N>\underline NN>N​, with N‾=0\underline N=0N​=0 for β>1\beta>1β>1 and N‾=8h/(vˉ2E[1/a])\underline N=8h/(\bar v^2E[1/a])N​=8h/(vˉ2E[1/a]) for β=1\beta=1β=1.

Formalization targets

Goal: Theorem 2 (p. 11)

  1. If h=0h=0h=0: w∗w^*w∗, we∗w^*_ewe∗​ and LW‾\overline{LW}LW are strictly decreasing in NNN.
  2. If h>0h>0h>0: w∗w^*w∗ is strictly decreasing in NNN on (N‾,∞)(\underline N,\infty)(N​,∞), and for each Q∈{D∗/S∗, we∗, LW‾}Q\in\{D^*/S^*,\ w^*_e,\ \overline{LW}\}Q∈{D∗/S∗, we∗​, LW} there is N0>max⁡(0,N‾)N_0>\max(0,\underline N)N0​>max(0,N​) with
Q strictly increasing on (N‾,N0),Q strictly decreasing on (N0,∞).Q\ \text{strictly increasing on}\ (\underline N,N_0),\qquad Q\ \text{strictly decreasing on}\ (N_0,\infty).Q strictly increasing on (N​,N0​),Q strictly decreasing on (N0​,∞).

The thresholds are not given in closed form; the goal asserts only their existence.

Milestones (in attack order)

  • Lemma 11 (p. OA1): unique positive market equilibrium for each admissible (p,w)(p,w)(p,w).
  • (8)–(9) (p. OA1): price, wage and profit as functions of (D,S)(D,S)(D,S); (10): LW‾=E[1/a] we2/4\overline{LW}=E[1/a]\,w_e^2/4LW=E[1/a]we2​/4.
  • Lemma 12 (p. OA2): the revenue-maximizing demand D~(S)\tilde D(S)D~(S) increases in SSS and Π(D~(S),S)\Pi(\tilde D(S),S)Π(D~(S),S) is single-peaked at S∗S^*S∗.
  • Lemma 1 (pp. 9–10): a unique optimum exists iff β>1\beta>1β>1 or β=1\beta=1β=1 and N>8h/(vˉ2E[1/a])N>8h/(\bar v^2E[1/a])N>8h/(vˉ2E[1/a]).
  • (17)–(18) (p. OA7): first-order conditions at the optimum.
  • Lemma 13 (p. OA4): S∗S^*S∗, D∗D^*D∗ increase in NNN; S∗→∞S^*\to\inftyS∗→∞, D∗→Λ/2D^*\to\Lambda/2D∗→Λ/2; S∗,D∗,ρ∗→0S^*,D^*,\rho^*\to0S∗,D∗,ρ∗→0 as N↓N‾N\downarrow\underline NN↓N​.
  • Case h=0h=0h=0 of the proof (pp. OA4–OA5): D∗=S∗=vˉ/(2(vˉ/Λ+2/(NE[1/a])))D^*=S^*=\bar v/(2(\bar v/\Lambda+2/(NE[1/a])))D∗=S∗=vˉ/(2(vˉ/Λ+2/(NE[1/a]))), w∗=vˉΛ/(vˉNE[1/a]+2Λ)w^*=\bar v\Lambda/(\bar vNE[1/a]+2\Lambda)w∗=vˉΛ/(vˉNE[1/a]+2Λ).
  • Parts (i)–(iv) of the proof (pp. OA5–OA7): the threshold statements for we∗w^*_ewe∗​, ρ∗\rho^*ρ∗ and LW‾\overline{LW}LW, the identity w∗=12vˉ−vˉD∗/Λw^*=\tfrac12\bar v-\bar vD^*/\Lambdaw∗=21​vˉ−vˉD∗/Λ, and the monotonicity of w∗w^*w∗.

Significance

The result. Theorem 2 separates two regimes. In a congestion regime (small pools, large delay cost hW∗hW^*hW∗) extra agents reduce delay, stimulate demand faster than supply grows, and raise utilization, the effective wage and average labor welfare: agents exert a positive externality on one another. In a cost regime (large pools) the demand effect is exhausted and extra agents only dilute work. The paper uses this to argue against blunt caps on platform labor pools and as the baseline for its analysis of nominal and effective wage floors (§5). Without delay sensitivity the classical conclusion, that a larger pool lowers welfare, is recovered exactly.

Formalizing it. The theorem is proved in the paper, with an online appendix of implicit-function and envelope arguments. None of it is formalized. This mission produces a machine-checked model of a platform with endogenous supply and congestion-dependent demand, the reduced (D,S)(D,S)(D,S) problem, and the comparative statics in NNN. The paper's Lemma 11 needs a correction at β=1\beta=1β=1 (see below), which is the kind of gap a formal development exposes.

Difficulty

The optimum is defined implicitly by two nonlinear first-order conditions, and there is no closed form for h>0h>0h>0. The obvious route, differentiating we∗=2S∗/(NE[1/a])w^*_e=2S^*/(NE[1/a])we∗​=2S∗/(NE[1/a]) in NNN, gives a quantity whose sign is the competition between the growth of S∗S^*S∗ and the factor 1/N1/N1/N, and that sign changes with NNN. Monotonicity arguments of the usual kind (supermodularity in (S,N)(S,N)(S,N), which does give that S∗S^*S∗ increases) say nothing about such ratios. The obstacle is single crossing of a quantity defined only implicitly through the optimum, not differentiability. The boundary behaviour at N↓N‾N\downarrow\underline NN↓N​ and N→∞N\to\inftyN→∞ (Lemma 13) is also needed and is itself nontrivial at β=1\beta=1β=1, where N‾>0\underline N>0N​>0.

Formalization scope

Everything lives in the namespace LaborWelfare.PoolSize, with one definitions file Setting.

  • Parameters. vˉ>0\bar v>0vˉ>0 (vbar), Λ>0\Lambda>0Λ>0 (Lam), h≥0h\ge0h≥0, real β≥1\beta\ge1β≥1 (beta), and FFF a probability measure on R\mathbb RR carried by [a‾,aˉ][\underline a,\bar a][a​,aˉ] with 0<a‾≤aˉ0<\underline a\le\bar a0<a​≤aˉ. The density of FFF is not required: only E[1/a]=∫a−1dFE[1/a]=\int a^{-1}dFE[1/a]=∫a−1dF enters any result. (D/S)β(D/S)^\beta(D/S)β is a real power.
  • NNN is a real number; all results are about functions of NNN on (N‾,∞)(\underline N,\infty)(N​,∞) (on (0,∞)(0,\infty)(0,∞) for h=0h=0h=0, where N‾=0\underline N=0N​=0).
  • Equilibrium. D>0D>0D>0, S>0S>0S>0, (2), (3)–(4), and D<SD<SD<S for h>0h>0h>0, D≤SD\le SD≤S for h=0h=0h=0; the uniform cdf is clamped to [0,1][0,1][0,1].
  • The optimum is an outcome maximizing (p−w)D(p-w)D(p−w)D over every price p≥0p\ge0p≥0 and wage w≥0w\ge0w≥0 subject to the market equilibrium. Statements quantify over a selection N↦N\mapstoN↦ optimal outcome, which Lemma 1 (and the h=0h=0h=0 milestone) show exists and is unique on the regime, so the monotonicity claims concern the platform's actual optimum and cannot be met by a degenerate outcome or a vacuous hypothesis.
  • Labor welfare is the integral of the page, not its closed form; the closed form is milestone (10).
  • Strictness. "Increasing" and "decreasing" are strict (StrictMonoOn, StrictAntiOn), as the paper fixes on p. 11; the thresholds are three separate numbers.
  • Correction. For w>0w>0w>0 and p<vˉp<\bar vp<vˉ, Lemma 11's positive equilibrium exists exactly when β>1\beta>1β>1 or p+2h/(wNE[1/a])<vˉp+2h/(wNE[1/a])<\bar vp+2h/(wNE[1/a])<vˉ. At β=1\beta=1β=1, the delay stays bounded away from zero as utilization vanishes; the paper omits that condition.

A complete development needs real analysis on implicitly defined functions (implicit-function theorem in one variable, envelope arguments, monotone limits) and an argmax framework for the two-variable profit; these are reusable for other platform-pricing papers. Proofs of any milestone are welcome, as are alternative routes to the threshold statements.

Selected references

  • S. Benjaafar, J.-Y. Ding, G. Kong, T. Taylor, Labor Welfare in On-Demand Service Platforms, Manufacturing & Service Operations Management, 2022 (online 2021). SSRN manuscript with online appendix: https://ssrn.com/abstract=3103447; https://doi.org/10.1287/msom.2020.0964
  • J. Bai, K. C. So, C. S. Tang, X. Chen, H. Wang, Coordinating supply and demand on an on-demand service platform with impatient customers, Manufacturing & Service Operations Management 21(3), 2019. https://doi.org/10.1287/msom.2018.0707
  • T. A. Taylor, On-demand service platforms, Manufacturing & Service Operations Management 20(4), 2018. https://doi.org/10.1287/msom.2017.0678
15 thms1 active userReviewed
Linear algebra·Captain: mikedeng1

Controlling Epidemic Spread: Reducing Economic Losses with Targeted Closures: The Effective Reproduction Rate Is Below One if and only if a Linear Matrix Inequality in the Activity Levels HoldsResearch Paper

Why targeted closures matter

A city can reduce contacts by limiting activity in selected neighborhoods. Such a policy affects more than the residents of the closed neighborhood: people travel across neighborhood boundaries, meet there, and carry infection back to their home neighborhoods. A citywide closure ignores these differences in movement and economic cost. Birge, Candogan, and Feng model these links and ask how to curb infection while preserving activity where it is most valuable. Their 2020 working paper gives separate optimization formulations for a large infection regime and a small infection regime; this mission concerns the latter. The paper later uses movement and economic data for New York City, while the theorem here is independent of those data. Birge, Candogan, and Feng, BFI Working Paper 2020-57.

The spatial epidemic model

There are nnn neighborhoods. Neighborhood iii has population Ni>0N_i>0Ni​>0, and its residents spend a fraction τij≥0\tau_{ij}\geq0τij​≥0 of their time in neighborhood jjj. The row sum ∑jτij\sum_j\tau_{ij}∑j​τij​ is at most one; remaining time is spent in isolation. An activity level xj∈[0,1]x_j\in[0,1]xj​∈[0,1] controls contact at destination jjj. This makes the number of infectious contacts depend on where residents of every neighborhood meet, not only where they live.

Each neighborhood has susceptible, exposed, clinically infected, subclinically infected, and recovered populations. The effective contact rate is β\betaβ, the progression rate from exposed to infected is κ\kappaκ, the clinical fraction is ρ\rhoρ, the relative infectiousness of subclinical cases is α\alphaα, the removal rate is γ\gammaγ, and the natural birth and death rate is μ\muμ. The model uses β,κ,γ,μ>0\beta,\kappa,\gamma,\mu>0β,κ,γ,μ>0 and 0≤ρ,α≤10\leq\rho,\alpha\leq10≤ρ,α≤1. At the disease-free equilibrium, the exposed and infected populations vanish and the susceptible population in neighborhood iii equals NiN_iNi​.

The paper separates the infected-coordinate dynamics into a vector F\mathcal FF of new infections and a vector V\mathcal VV of other transfers. In the coordinate order exposed, clinical, subclinical, their Jacobians at the disease-free equilibrium are F(x)F(x)F(x) and VVV, each a 3n×3n3n\times3n3n×3n real matrix. Define λ1(A)\lambda_1(A)λ1​(A) as the greatest modulus of a complex eigenvalue of AAA. The effective reproduction rate in condition (6) is Re=λ1(F(x)V−1)\mathcal R_e=\lambda_1(F(x)V^{-1})Re​=λ1​(F(x)V−1). This is a property of the full infected-coordinate system. Model and condition (6), pp. 7–13.

Formalization targets

The intermediate targets identify the displayed inverse of VVV, the smaller matrix sharing the nonzero eigenvalues of F(x)V−1F(x)V^{-1}F(x)V−1, and its diagonal similarity to a symmetric matrix. They connect condition (6) with a bound on the spectral radius of the reduced matrix. Each target is taken from the proof of Lemma 2.1 in the working paper.

The goal is the paper's Lemma 2.1. Let T=(τij)T=(\tau_{ij})T=(τij​), X=diag⁡(xi)X=\operatorname{diag}(x_i)X=diag(xi​), N^=diag⁡(Ni)\widehat N=\operatorname{diag}(N_i)N=diag(Ni​), and Λ=diag⁡((∑lNlτlk)−1)\Lambda=\operatorname{diag}((\sum_l N_l\tau_{lk})^{-1})Λ=diag((∑l​Nl​τlk​)−1). With A≺BA\prec BA≺B meaning that B−AB-AB−A is positive definite, the target is

Re<1⟺TXΛTT≺N^−1(μ+κ)(μ+γ)(ρ+α(1−ρ))βκ.\mathcal R_e<1\quad\Longleftrightarrow\quad TX\Lambda T^{\mathsf T}\prec \widehat N^{-1}\frac{(\mu+\kappa)(\mu+\gamma)}{(\rho+\alpha(1-\rho))\beta\kappa}.Re​<1⟺TXΛTT≺N−1(ρ+α(1−ρ))βκ(μ+κ)(μ+γ)​.

The strict inequality is the stability condition stated in the paper. The semidefinite program labelled (P2) uses a non-strict inequality to close its feasible region; that optimization problem is outside this mission. Lemma 2.1 and (P2), pp. 13–14.

What the result provides

The left side of Lemma 2.1 is a spectral test on a 3n×3n3n\times3n3n×3n next generation matrix. The right side is a linear matrix inequality in the activity variables because XXX is diagonal and depends linearly on xxx. It therefore supplies the mathematical constraint used by the paper's small infection optimization problem. The theorem states equivalence for the paper's spatial epidemic model; it does not assert that the optimization problem has a solution or that a chosen activity vector minimizes economic loss. Birge, Candogan, and Feng, pp. 13–14.

The working paper proves the lemma in ordinary mathematics. The formalization task is to give its matrices and hypotheses precise Lean definitions, then prove the spectral and positive definiteness statements. The definition of the reproduction rate and the intermediate matrix identities can also support later formal statements about related compartmental models. No machine-checked proof of these draft items is claimed here.

Where the difficulty lies

The matrix F(x)V−1F(x)V^{-1}F(x)V−1 is generally nonsymmetric and has three blocks of infected coordinates. A direct comparison between its entries and a symmetric positive definite matrix does not express the theorem. The paper first relates its nonzero eigenvalues to those of an n×nn\times nn×n matrix, then uses a diagonal similarity to obtain a symmetric matrix. The final equivalence relies on the activity and mobility bounds that make this symmetric matrix positive semidefinite. A matrix statement that starts directly from the reduced symmetric form would omit the spectral reduction at the heart of the lemma.

Formalization scope

Lean uses Fin n for neighborhoods and Fin 3 × Fin n for infected coordinates, in exposed, clinical, subclinical order. All model matrices are real; λ1\lambda_1λ1​ takes the spectral radius after mapping a matrix to complex entries, so complex eigenvalues are included. The inverse in Re\mathcal R_eRe​ is Mathlib's matrix inverse, and A≺BA\prec BA≺B is represented by positive definiteness of B−AB-AB−A. The only numerical instance in the sanity file has one neighborhood, but the theorem quantifies over arbitrary finite nnn.

The standing assumptions include positive populations and rates, nonnegative travel fractions with row sums at most one, clinical and subclinical fractions in [0,1][0,1][0,1], and activity levels in [0,1][0,1][0,1]. The formalization also requires ρ+α(1−ρ)>0\rho+\alpha(1-\rho)>0ρ+α(1−ρ)>0: the paper divides by this quantity without stating the exclusion of ρ=α=0\rho=\alpha=0ρ=α=0. If no one visits destination kkk, its mixing denominator is zero. Under the standing assumptions that entire column of TTT vanishes, so Lean's value 1/0=01/0=01/0=0 has no effect on the theorem. These conventions are part of the model, not optional proof shortcuts.

The proposal includes the full F(x)F(x)F(x) and VVV, their displayed inverse, the paper's reduced matrix, the diagonal similarity, and the Loewner form. Contributions may develop spectral identities for block matrices, diagonal similarity, and positive definite matrix bounds. The large infection linear program, the stability theorem cited from earlier epidemic literature, and the New York City numerical application are beyond this target.

Selected references

  • John R. Birge, Ozan Candogan, and Yiding Feng, Controlling Epidemic Spread: Reducing Economic Losses with Targeted Closures, Becker Friedman Institute Working Paper No. 2020-57, May 18, 2020 (preliminary), especially §2.1, §2.3, and Lemma 2.1. Working paper PDF.
  • John R. Birge, Ozan Candogan, and Yiding Feng, Controlling Epidemic Spread: Reducing Economic Losses with Targeted Closures, Management Science 68(5), 2022, pp. 3175–3195. DOI. This mission's statements and page numbers follow the 2020 working paper.
7 thms1 active userReviewed
Algorithmic Game TheoryEconomics·Captain: mikedeng1

Revenue-Sharing Allocation Strategies for Two-Sided Media Platforms: Pro-Rata vs. User-Centric: The Ratios π2j/π1j for Which Each Rule Sustains All Artists Form Two Explicit IntervalsResearch Paper

Motivation

Music and podcast streaming services such as Spotify, Apple Music and Deezer charge every subscriber the same monthly fee and pass a share of the collected revenue to the artists. How that share is split is contested. The industry default, pro-rata, pools all fees and pays each artist her share of total streams. The alternative, user-centric, splits each subscriber's fee according to that subscriber's own listening. The two rules pay identical amounts when all users stream the same volume, and differ whenever heavy users and light users listen to different artists; which artists gain and which lose under each rule is the subject of an active industry debate.

Alaei, Makhdoumi, Malekian and Pekeč (SSRN 3645521; published in Management Science 68(12), 2022, doi:10.1287/mnsc.2022.4307) model the question as a two-sided market in which every artist has an outside option: she can leave the platform and sell her content directly. A rule is only viable if it pays each artist at least what she would earn alone. The paper's first main result, Theorem 1, says exactly when each rule can keep every artist on the platform. This mission formalizes that theorem and the chain of results its proof rests on, citing the authors' submitted manuscript (SSRN 3645521) for all page and result numbers.

Setting

There is a finite set M\mathcal MM of artists and two user types i∈{1,2}i\in\{1,2\}i∈{1,2}. Type iii users consume content at usage rate λi\lambda_iλi​ and form a population of mass qiq_iqi​; a type iii user wants artist jjj with probability πij\pi_{ij}πij​ per use. Throughout,

λ2>λ1>0,q1,q2∈(0,1), q1+q2=1,πij≥0, ∑jπij=1.\lambda_2>\lambda_1>0,\qquad q_1,q_2\in(0,1),\ q_1+q_2=1,\qquad \pi_{ij}\ge0,\ \sum_{j}\pi_{ij}=1 .λ2​>λ1​>0,q1​,q2​∈(0,1), q1​+q2​=1,πij​≥0, j∑​πij​=1.

Type 2 users are the heavy users.

Going solo. If artist jjj sells directly at price rrr, a type iii user buys when her value λiπij\lambda_i\pi_{ij}λi​πij​ is at least rrr, so the demand d(j,r)d(j,r)d(j,r) is 000, q1q_1q1​, q2q_2q2​ or 111 according to which of λ1π1j≥r\lambda_1\pi_{1j}\ge rλ1​π1j​≥r and λ2π2j≥r\lambda_2\pi_{2j}\ge rλ2​π2j​≥r hold. The solo revenue is

Rj=max⁡rr d(j,r).R_j=\max_{r} r\,d(j,r).Rj​=rmax​rd(j,r).

The platform. With all artists on the platform and both types subscribing, the fee is Fee=min⁡iλi∑jπij\mathrm{Fee}=\min_i\lambda_i\sum_j\pi_{ij}Fee=mini​λi​∑j​πij​ and the revenue is Rev=Fee (q1+q2)\mathrm{Rev}=\mathrm{Fee}\,(q_1+q_2)Rev=Fee(q1​+q2​). A payout rate β∈[0,1]\beta\in[0,1]β∈[0,1] is the fraction of revenue paid to artists. The pro-rata payment (Definition 1) and the user-centric payment (Definition 2) to artist jjj are

pjpr(β)=β Rev ∑iqiλiπij∑iqiλi∑kπik,pjuc(β)=β∑iqi Fee λiπij∑kλiπik.p^{\mathrm{pr}}_j(\beta)=\beta\,\mathrm{Rev}\,\frac{\sum_i q_i\lambda_i\pi_{ij}}{\sum_i q_i\lambda_i\sum_k\pi_{ik}},\qquad p^{\mathrm{uc}}_j(\beta)=\beta\sum_i q_i\,\mathrm{Fee}\,\frac{\lambda_i\pi_{ij}}{\sum_k\lambda_i\pi_{ik}} .pjpr​(β)=βRev∑i​qi​λi​∑k​πik​∑i​qi​λi​πij​​,pjuc​(β)=βi∑​qi​Fee∑k​λi​πik​λi​πij​​.

The platform is sustainable with a rule if some β∈[0,1]\beta\in[0,1]β∈[0,1] makes every artist's payment at least RjR_jRj​.

Assumption 1 is λ1≥λ2q2\lambda_1\ge\lambda_2q_2λ1​≥λ2​q2​: it rules out the case where the platform would rather serve heavy users alone, in which the two rules coincide. Artists are sorted by the ratio π2j/π1j\pi_{2j}/\pi_{1j}π2j​/π1j​ into four regions

W=[0,q1λ1λ2),Z=[q1λ1λ2,λ1λ2),Y=[λ1λ2,λ1λ2q2),X=[λ1λ2q2,∞).W=[0,\tfrac{q_1\lambda_1}{\lambda_2}),\quad Z=[\tfrac{q_1\lambda_1}{\lambda_2},\tfrac{\lambda_1}{\lambda_2}),\quad Y=[\tfrac{\lambda_1}{\lambda_2},\tfrac{\lambda_1}{\lambda_2q_2}),\quad X=[\tfrac{\lambda_1}{\lambda_2q_2},\infty).W=[0,λ2​q1​λ1​​),Z=[λ2​q1​λ1​​,λ2​λ1​​),Y=[λ2​λ1​​,λ2​q2​λ1​​),X=[λ2​q2​λ1​​,∞).

Formalization targets

Goal: Theorem 1 (p. 17)

Under Assumption 1 (and π1j>0\pi_{1j}>0π1j​>0 for all jjj), the platform is sustainable with pro-rata if and only if, for all jjj,

π2jπ1j∈[q1(λ2−λ1)λ2,q1λ12λ22q2+λ1λ2(q1−q2)]∪[1,q1λ12λ2q22(λ2−λ1)],\frac{\pi_{2j}}{\pi_{1j}}\in\Bigl[\tfrac{q_1(\lambda_2-\lambda_1)}{\lambda_2},\tfrac{q_1\lambda_1^2}{\lambda_2^2q_2+\lambda_1\lambda_2(q_1-q_2)}\Bigr]\cup\Bigl[1,\tfrac{q_1\lambda_1^2}{\lambda_2q_2^2(\lambda_2-\lambda_1)}\Bigr],π1j​π2j​​∈[λ2​q1​(λ2​−λ1​)​,λ22​q2​+λ1​λ2​(q1​−q2​)q1​λ12​​]∪[1,λ2​q22​(λ2​−λ1​)q1​λ12​​],

and sustainable with user-centric if and only if, for all jjj,

π2jπ1j∈[0,q1λ1λ2−q2λ1]∪[1,q1λ1q2(λ2−λ1)].\frac{\pi_{2j}}{\pi_{1j}}\in\Bigl[0,\tfrac{q_1\lambda_1}{\lambda_2-q_2\lambda_1}\Bigr]\cup\Bigl[1,\tfrac{q_1\lambda_1}{q_2(\lambda_2-\lambda_1)}\Bigr].π1j​π2j​​∈[0,λ2​−q2​λ1​q1​λ1​​]∪[1,q2​(λ2​−λ1​)q1​λ1​​].

Milestones

  1. Lemma 1 (p. 14): RjR_jRj​ equals q1λ1π1jq_1\lambda_1\pi_{1j}q1​λ1​π1j​, λ2π2j\lambda_2\pi_{2j}λ2​π2j​, λ1π1j\lambda_1\pi_{1j}λ1​π1j​, q2λ2π2jq_2\lambda_2\pi_{2j}q2​λ2​π2j​ on W,Z,Y,XW,Z,Y,XW,Z,Y,X, attained at the stated fee with the stated subscribers.
  2. Eq. (10) (p. 16): every artist is paid at least RjR_jRj​ under pro-rata with payout β\betaβ iff β≥βpr=q1λ1+q2λ2λ1max⁡jRjq1λ1π1j+q2λ2π2j\beta\ge\beta^{\mathrm{pr}}=\frac{q_1\lambda_1+q_2\lambda_2}{\lambda_1}\max_j\frac{R_j}{q_1\lambda_1\pi_{1j}+q_2\lambda_2\pi_{2j}}β≥βpr=λ1​q1​λ1​+q2​λ2​​maxj​q1​λ1​π1j​+q2​λ2​π2j​Rj​​.
  3. Lemma 6 (p. 45): pro-rata sustainability is equivalent to four region-wise bounds on π2j/π1j\pi_{2j}/\pi_{1j}π2j​/π1j​.
  4. Eq. (11) (p. 16): the same for user-centric with βuc=1λ1max⁡jRjq1π1j+q2π2j\beta^{\mathrm{uc}}=\frac1{\lambda_1}\max_j\frac{R_j}{q_1\pi_{1j}+q_2\pi_{2j}}βuc=λ1​1​maxj​q1​π1j​+q2​π2j​Rj​​.
  5. Lemma 7 (p. 45): user-centric sustainability is equivalent to three region-wise bounds.

Significance

Theorem 1 separates the two rules by the kind of artist they can retain. Both rules lose artists whose audience is concentrated among heavy users, since such an artist captures those users' full fee by going solo; pro-rata tolerates a larger ratio π2j/π1j\pi_{2j}/\pi_{1j}π2j​/π1j​ before losing her. User-centric, in turn, can retain artists listened to mainly by light users, which pro-rata pays too little. The later results of the paper (which rule artists and the platform prefer, Propositions 2–3, and the optimal artist sets of §3.4) all start from this characterization.

The result is proved in the paper. No part of it, nor the underlying model of subscription fees, outside options and payout rates, has a machine-checked proof. The platform already has an axiomatic treatment of pro-rata and user-centric indices on stream-count tables; this mission adds the economic model with prices, outside options and payout rates, and its sustainability theorem. A complete development gives a checked version of the paper's case analysis, including its boundary cases.

Difficulty

Each step is elementary, and the difficulty lies in the bookkeeping. The solo revenue is a supremum of a step function times the price, and identifying its value requires comparing three candidates on four regions with half-open boundaries. The payments must be simplified from their general form (with sums and a minimum fee) to the closed forms of (10)–(11), which uses the standing assumptions and Assumption 1. Finally, combining four region-wise conditions into two intervals does not go region by region: when λ2>2λ1\lambda_2>2\lambda_1λ2​>2λ1​ the first pro-rata interval is empty and no artist in WWW or ZZZ can be sustained, and the upper intervals agree with the region boundary λ1/(λ2q2)\lambda_1/(\lambda_2q_2)λ1​/(λ2​q2​) only because of Assumption 1. A transcription that treats the printed intervals as non-empty, or drops Assumption 1, gives a false statement.

Formalization scope

Everything sits in the namespace ProRataUC.Sustain, in one definitions file and six theorem files. The setting is a structure Setting M over a type M of artists with [Fintype M], with real fields lam1 lam2 q1 q2 and listening probabilities pi1 pi2 : M → ℝ; the standing assumptions are one predicate Setting.Standing, and Assumption 1 is a separate predicate. Conventions:

  • The solo revenue is the real supremum of r d(j,r)r\,d(j,r)rd(j,r) over all real rrr, with ddd the printed four-case demand; the family is bounded above, so the supremum is genuine.
  • The fee and revenue are those with all artists on the platform and both types subscribing, with the sums ∑jπij\sum_j\pi_{ij}∑j​πij​ as printed, not hard-coded to λ1\lambda_1λ1​.
  • Sustainability is the existence of a payout rate β∈[0,1]\beta\in[0,1]β∈[0,1] under which every artist's payment is at least RjR_jRj​, not the inequality βpr≤1\beta^{\mathrm{pr}}\le1βpr≤1 and not the interval condition; equations (10)–(11) are the bridge.
  • βpr\beta^{\mathrm{pr}}βpr and βuc\beta^{\mathrm{uc}}βuc use a maximum over a nonempty finite set, so (10)–(11) assume at least one artist.
  • The intervals are closed in Theorem 1 and half-open in Lemma 1 and the regions, as printed.

Added hypotheses, each disclosed in the statement: π1j>0\pi_{1j}>0π1j​>0 (the paper's ratio presupposes it; Lean's division by zero would put an artist heard only by heavy users in the region [0,⋅)[0,\cdot)[0,⋅) and make part (b) false); Assumption 1 in equations (10)–(11) and Lemmas 6–7, whose proofs use it although the lemmas say "for a given λ\lambdaλ and qqq".

Formalizations that would trivialize the goal are ruled out: sustainability is not defined as the printed intervals or as βpr≤1\beta^{\mathrm{pr}}\le1βpr≤1, RjR_jRj​ is not defined by Lemma 1's case formula, π1j>0\pi_{1j}>0π1j​>0 is explicit, and Assumption 1 is kept.

Contributions welcome: proofs of the milestones in order, reusable lemmas on suprema of step functions times price (the two-type monopoly pricing problem) and on Finset.sup' of ratios, and the interval bookkeeping of Theorem 1.

Selected references

  • S. Alaei, A. Makhdoumi, A. Malekian, S. Pekeč, Revenue-sharing allocation strategies for two-sided media platforms: Pro-rata vs. user-centric, submitted manuscript, SSRN 3645521 (2020); published in Management Science 68(12):8699–8721, 2022. https://ssrn.com/abstract=3645521 , https://doi.org/10.1287/mnsc.2022.4307
  • G. Bergantiños, J. D. Moreno-Ternero, Revenue sharing at music streaming platforms, Management Science, 2025 (axiomatic comparison of pro-rata and user-centric); arXiv:2310.11861. https://arxiv.org/abs/2310.11861
7 thms1 active userReviewed
Algorithmic Game TheoryMechanism Design·Captain: mikedeng1

Preference Elicitation for Participatory Budgeting: Randomized Threshold Approval Votes Achieve Distortion O(log² m)Research Paper

Motivation

In participatory budgeting, residents of a city vote on how to spend part of a public budget on a list of proposed projects, each with a known cost. The practice has been adopted by cities in North America and Europe, including New York City, Chicago, Madrid and Paris, the last two each spending more than 100 million euros through participatory budgets in 2017 (Cabannes, 2004; paper §1). Every deployment must choose a ballot: voters may be asked to approve some projects, to rank them, or to fill a knapsack within the budget. The ballot determines how much information about the residents' preferences reaches the aggregation step, and so how much welfare the chosen set of projects can guarantee.

The implicit utilitarian framework of Procaccia and Rosenschein (2006) makes this question quantitative. Voters are assumed to hold cardinal utilities, which the ballot reveals only partially; the quality of a ballot-plus-aggregation pair is measured by its distortion, the worst-case ratio between the best achievable welfare and the welfare it actually obtains. Benadè, Nath, Procaccia and Shah (Management Science, 2021; AAAI 2017 conference version) apply this lens to participatory budgeting and compare the input formats used in practice. Their main positive result is that threshold approval votes with a randomly drawn threshold reach distortion O(log⁡2m)O(\log^2 m)O(log2m), against mmm for a rule that uses no information at all and Ω(m)\Omega(\sqrt m)Ω(m​) for any deterministic threshold.

Setting

There are nnn voters N=[n]N=[n]N=[n] and a set AAA of mmm alternatives (projects). Alternative aaa costs ca∈[0,1]c_a\in[0,1]ca​∈[0,1], and the budget is 111. Write c(S)=∑a∈Scac(S)=\sum_{a\in S}c_ac(S)=∑a∈S​ca​, and let

Fc={S⊆A: c(S)≤1, c(T)>1 for all S⊊T⊆A}\mathcal F_c=\{S\subseteq A:\ c(S)\le 1,\ c(T)>1\text{ for all }S\subsetneq T\subseteq A\}Fc​={S⊆A: c(S)≤1, c(T)>1 for all S⊊T⊆A}

be the family of inclusion-maximal budget-feasible sets.

Voter iii has a utility vi(a)≥0v_i(a)\ge 0vi​(a)≥0 for each alternative, utilities are additive over sets, and they are normalized: ∑a∈Avi(a)=1\sum_{a\in A}v_i(a)=1∑a∈A​vi​(a)=1. The social welfare of a set is sw(S,v⃗)=∑a∈S∑i∈Nvi(a)\mathrm{sw}(S,\vec v)=\sum_{a\in S}\sum_{i\in N}v_i(a)sw(S,v)=∑a∈S​∑i∈N​vi​(a).

A threshold approval vote at threshold ttt is the set τi={a∈A:vi(a)≥t}\tau_i=\{a\in A: v_i(a)\ge t\}τi​={a∈A:vi​(a)≥t}. A randomized threshold rule chooses a threshold ttt at random and maps the resulting profile ρ⃗(v⃗,t)=(τ1,…,τn)\vec\rho(\vec v,t)=(\tau_1,\dots,\tau_n)ρ​(v,t)=(τ1​,…,τn​) to a (possibly random) set in Fc\mathcal F_cFc​. In the formalization, a rule is a finite mixture of arms; each arm has a probability, a threshold, and a deterministic aggregation map from approval profiles to Fc\mathcal F_cFc​ (ThresholdRule), and its expected welfare is the weighted sum over arms (ThresholdRule.expWelfare). The rule never sees the utilities, only the approval sets. The distortion of a rule is bounded by DDD when

max⁡T∈Fcsw(T,v⃗) ≤ D⋅E[sw(rule,v⃗)]for every normalized profile v⃗,\max_{T\in\mathcal F_c}\mathrm{sw}(T,\vec v)\ \le\ D\cdot\mathbb E[\mathrm{sw}(\text{rule},\vec v)]\quad\text{for every normalized profile }\vec v,T∈Fc​max​sw(T,v) ≤ D⋅E[sw(rule,v)]for every normalized profile v,

and the distortion associated with an input format is that of the best rule for it.

Formalization targets

Goal: Theorem 3.5

∃ C>0  ∀n, ∀A with m≥2, ∀c∈[0,1]A: ∃ rule R  ∀v⃗, ∀T∈Fc:sw(T,v⃗)≤Clog⁡22m⋅E[sw(R,v⃗)].\exists\,C>0\ \ \forall n,\ \forall A\text{ with }m\ge 2,\ \forall c\in[0,1]^A:\ \exists\text{ rule }R\ \ \forall\vec v,\ \forall T\in\mathcal F_c:\quad \mathrm{sw}(T,\vec v)\le C\log_2^2 m\cdot\mathbb E[\mathrm{sw}(R,\vec v)].∃C>0  ∀n, ∀A with m≥2, ∀c∈[0,1]A: ∃ rule R  ∀v, ∀T∈Fc​:sw(T,v)≤Clog22​m⋅E[sw(R,v)].

The goal fixes only the shape O(log⁡2m)O(\log^2 m)O(log2m), with an absolute constant chosen before nnn, AAA and the costs, and it covers every m≥2m\ge 2m≥2.

Milestones

For m=2Km=2^Km=2K, with ℓj=2j−1/m2\ell_j=2^{j-1}/m^2ℓj​=2j−1/m2, uj=2j/m2u_j=2^j/m^2uj​=2j/m2, Ij=(ℓj,uj]I_j=(\ell_j,u_j]Ij​=(ℓj​,uj​], cost buckets Sj={a:ca∈Ij}S_j=\{a:c_a\in I_j\}Sj​={a:ca​∈Ij​} (S0S_0S0​ uses I0=[0,1/m2]I_0=[0,1/m^2]I0​=[0,1/m2]) and counts nja=∣{i:vi(a)∈Ij}∣n^a_j=|\{i:v_i(a)\in I_j\}|nja​=∣{i:vi​(a)∈Ij​}∣. In the statements involving S∗S^*S∗, it is an optimal member of Fc\mathcal F_cFc​ for the given utility profile (IsWelfareOptimal).

  1. Observation 3.1: some rule has distortion at most mmm.
  2. Eq. (6): ∑j≥1njaℓj≤sw(a,v⃗)≤n/m2+∑j≥1njauj\sum_{j\ge1}n^a_j\ell_j\le\mathrm{sw}(a,\vec v)\le n/m^2+\sum_{j\ge1}n^a_ju_j∑j≥1​nja​ℓj​≤sw(a,v)≤n/m2+∑j≥1​nja​uj​.
  3. At most 1/uj1/u_j1/uj​ alternatives from SjS_jSj​ are budget-feasible.
  4. Eq. (7): the same sandwich for sw(S∗∩Sk,v⃗)\mathrm{sw}(S^*\cap S_k,\vec v)sw(S∗∩Sk​,v).
  5. Eq. (8): a greedy top set Bj,k⊆SkB_{j,k}\subseteq S_kBj,k​⊆Sk​ at threshold ℓj\ell_jℓj​ has sw(Bj,k,v⃗)≥14nj,k∗uj\mathrm{sw}(B_{j,k},\vec v)\ge\frac14 n^*_{j,k}u_jsw(Bj,k​,v)≥41​nj,k∗​uj​.
  6. Mechanism A's expected welfare is at least (sw(S∗∖S0,v⃗)−n/m)/(16log⁡2m)(\mathrm{sw}(S^*\setminus S_0,\vec v)-n/m)/(16\log^2 m)(sw(S∗∖S0​,v)−n/m)/(16log2m).
  7. Mechanism B: S0S_0S0​ is feasible and covers S∗∩S0S^*\cap S_0S∗∩S0​.
  8. The explicit form for m=2Km=2^Km=2K: distortion at most 2+16log⁡2m2+16\log^2 m2+16log2m.

Significance

The theorem separates randomized from deterministic elicitation: with a fixed threshold no aggregation rule beats Ω(m)\Omega(\sqrt m)Ω(m​) (Theorem 3.4 of the paper), while drawing the threshold at random brings the distortion to O(log⁡2m)O(\log^2 m)O(log2m), close to the Ω(log⁡m/log⁡log⁡m)\Omega(\log m/\log\log m)Ω(logm/loglogm) lower bound of Theorem 3.6. It is the paper's argument for threshold approval as a practical, low-burden ballot, and it supports the empirical comparison in the paper's §5 between input formats on data from real participatory budgeting elections.

The result is proved in the paper; no machine-checked version is known. A formalization yields a reusable layer for implicit utilitarian voting with budgets (feasible-set families, normalized additive profiles, randomized rules over partial ballots, distortion), on which the paper's other bounds (Theorem 3.3 for rankings, Theorems 3.4 and 3.6 for lower bounds) and the related literature on distortion can be stated. It also makes explicit the reading of distortion under a random threshold that the proof actually uses.

Difficulty

The obvious rule picks the threshold and then funds the most-approved projects greedily. It fails because one threshold cannot serve both high-utility and low-utility voters, and because popular projects may be expensive: approval counts say nothing about costs. Mixing over thresholds alone does not help either; the ratio for a single badly chosen threshold can be of order mmm. What has to be shown is that one fixed mixture, chosen before any utilities are seen, recovers a 1/O(log⁡2m)1/O(\log^2 m)1/O(log2m) share of the optimum simultaneously for every profile, despite rounding of utilities into intervals, budget constraints within each cost class, and the loss of everything below the smallest threshold.

Formalization scope

All objects live in PBElicit.RandThresh (definition file PBElicit.RandThresh.Setting). Voters are Fin n; alternatives are a Fintype with mmm = Fintype.card A. Utilities are real, nonnegative and normalized (IsUtilityProfile); costs satisfy 0≤ca≤10\le c_a\le 10≤ca​≤1, which the paper uses when it calls a single alternative budget-feasible. Rules must return sets of Fc\mathcal F_cFc​ (IsMaxFeasible); the paper's mechanisms return feasible sets, which are extended to maximal ones. The maximum over Fc\mathcal F_cFc​ is stated as "for every T∈FcT\in\mathcal F_cT∈Fc​"; distortion bounds are multiplicative. Logarithms are base 222. The milestones take m=2Km=2^Km=2K, as the paper's proof does; the goal does not.

Distortion is read as the optimum over the total expected welfare (threshold and aggregation together), which is what the proof of Theorem 3.5 bounds; the paper's displayed definition puts the expectation over the threshold outside the ratio. The interval bound (8) is formalized as its outer inequality, because the printed middle step fails when a utility equals a threshold exactly.

The rule is fixed before the utilities, sees only threshold approval votes, and must return an inclusion-maximal feasible set; choosing the rule after the profile (which makes distortion 111 trivially) is excluded by the quantifier order, and the uniform rule of Observation 3.1 only gives distortion mmm, so a Clog⁡2mC\log^2 mClog2m bound is a genuine improvement.

Welcome contributions: proofs of the milestones, the reduction from general mmm to the next power of 222, and further results of the paper stated on the same definitions.

Selected references

  • G. Benadè, S. Nath, A. D. Procaccia, N. Shah, Preference Elicitation for Participatory Budgeting, Management Science, 2021 (AAAI 2017 conference version). https://doi.org/10.1287/mnsc.2020.3666
  • A. D. Procaccia, J. S. Rosenschein, The Distortion of Cardinal Preferences in Voting, Cooperative Information Agents (CIA), LNCS 4149, pp. 317–331, 2006. https://doi.org/10.1007/11839354_23
  • I. Caragiannis, A. D. Procaccia, Voting Almost Maximizes Social Welfare Despite Limited Communication, Artificial Intelligence 175(9–10), pp. 1655–1671, 2011. https://doi.org/10.1016/j.artint.2011.03.005
  • C. Boutilier, I. Caragiannis, S. Haber, T. Lu, A. D. Procaccia, O. Sheffet, Optimal Social Choice Functions: A Utilitarian View, Artificial Intelligence 227, pp. 190–213, 2015. https://doi.org/10.1016/j.artint.2015.06.003
  • Y. Cabannes, Participatory Budgeting: A Significant Contribution to Participatory Democracy, Environment and Urbanization 16(1), pp. 27–46, 2004. https://doi.org/10.1177/095624780401600104
10 thms1 active userReviewed
Algorithmic Game TheoryCombinatorics·Captain: mikedeng1

Best of Both Worlds: Ex-Ante and Ex-Post Fairness in Resource Allocation: Recursive Probabilistic Serial Is Ex-Ante Envy-Free and Ex-Post Envy-Free up to One GoodResearch Paper

Motivation

Many allocation problems hand out indivisible goods without money: course seats, shifts, donated food, inherited items. An integral allocation of such goods can rarely be envy-free. With one good and two agents, someone always envies. The standard remedies are of two kinds. Randomization restores fairness in expectation: the probabilistic serial rule of Bogomolnaia and Moulin (J. Econ. Theory 2001) yields a fractional allocation that is envy-free and can be implemented as a lottery over integral allocations. Approximation restores fairness after the draw: envy-freeness up to one good (EF1), introduced by Lipton, Markakis, Mossel and Saberi (EC 2004) and named by Budish (J. Political Economy 2011), asks only that any envy disappear once a single good is removed from the envied bundle.

The two remedies pull apart. Giving all goods to one agent chosen uniformly at random is envy-free ex-ante, but after the draw one agent holds everything. Round-robin picking with a uniformly random order is EF1 after every draw, and proportional but not envy-free in expectation (Proposition 2 and Example 1 of the paper). Freeman, Shah and Vaish (arXiv:2005.14122, 2020; with Aziz, Operations Research 2024) ask whether one lottery can be fair in both senses at once. They show that it can, and that the answer depends on how the lottery is built, not only on which fractional allocation it implements.

Setting

There is a finite set NNN of n≥1n \ge 1n≥1 agents and a finite set MMM of mmm goods. Agent iii has an additive valuation with values vi,g≥0v_{i,g} \ge 0vi,g​≥0. A fractional allocation is a matrix XXX with Xi,g∈[0,1]X_{i,g} \in [0,1]Xi,g​∈[0,1] the share of good ggg given to agent iii, and vi(Xi)=∑gXi,gvi,gv_i(X_i) = \sum_g X_{i,g} v_{i,g}vi​(Xi​)=∑g​Xi,g​vi,g​. An integral allocation gives each agent a bundle Ai⊆MA_i \subseteq MAi​⊆M, with bundles pairwise disjoint, and vi(Ai)=∑g∈Aivi,gv_i(A_i) = \sum_{g \in A_i} v_{i,g}vi​(Ai​)=∑g∈Ai​​vi,g​. A randomized allocation is a finite lottery {(pk,Ak)}k\{(p^k, A^k)\}_k{(pk,Ak)}k​ over integral allocations. It implements the fractional allocation ∑kpkAk\sum_k p^k A^k∑k​pkAk.

XXX is envy-free (EF) if vi(Xi)≥vi(Xh)v_i(X_i) \ge v_i(X_h)vi​(Xi​)≥vi​(Xh​) for all i,hi, hi,h. An integral allocation AAA is EF1 if for all i,hi, hi,h with Ah≠∅A_h \neq \emptysetAh​=∅ some j∈Ahj \in A_hj∈Ah​ satisfies vi(Ai)≥vi(Ah∖{j})v_i(A_i) \ge v_i(A_h \setminus \{j\})vi​(Ai​)≥vi​(Ah​∖{j}). A lottery is ex-ante EF if the allocation it implements is EF. It is ex-post EF1 if every allocation it draws is EF1.

The eating protocol Eating(M′,1)\mathrm{Eating}(M', 1)Eating(M′,1) runs for one unit of time on the goods M′M'M′, each of mass one. At every moment each agent eats her favorite good among those not yet fully consumed, at rate one. Recursive Probabilistic Serial (Algorithm 1) proceeds in ⌈m/n⌉\lceil m/n \rceil⌈m/n⌉ stages, starting from M1=MM^1 = MM1=M. Stage ttt computes Xt=Eating(Mt,1)X^t = \mathrm{Eating}(M^t, 1)Xt=Eating(Mt,1) and decomposes it as Xt=∑kwt,kAt,kX^t = \sum_k w^{t,k} A^{t,k}Xt=∑k​wt,kAt,k by the generalized Birkhoff–von Neumann theorem, Theorem 1 of the paper. It then draws Bt=At,kB^t = A^{t,k}Bt=At,k with probability wt,kw^{t,k}wt,k and continues on Mt+1=Mt∖⋃iBitM^{t+1} = M^t \setminus \bigcup_i B^t_iMt+1=Mt∖⋃i​Bit​. The final bundles are Bi=⋃tBitB_i = \bigcup_t B^t_iBi​=⋃t​Bit​.

Formalization targets

Goal: Theorem 2 (p. 11)

For every instance, every way of breaking ties between equally valued goods, and every decomposition rule satisfying the conclusions of Theorem 1,

the lottery of Algorithm 1 is ex-ante envy-free and ex-post EF1.\text{the lottery of Algorithm 1 is ex-ante envy-free and ex-post EF1.}the lottery of Algorithm 1 is ex-ante envy-free and ex-post EF1.

Milestones

  1. Theorem 1 (p. 9). A matrix with entries in [0,1][0,1][0,1] and row and column sums at most 111 is a convex combination of 0/10/10/1 allocations. Each of their row and column sums is the floor or the ceiling of the corresponding sum of XXX.
  2. One eating stage (§3.1, p. 9). Eating(M′,1)\mathrm{Eating}(M', 1)Eating(M′,1) satisfies the hypotheses of Theorem 1 and is supported on M′M'M′. If ∣M′∣≥n|M'| \ge n∣M′∣≥n, every agent eats a total mass of exactly one.
  3. One item per stage (§3.1, p. 9). Every allocation of the decomposition uses only goods of M′M'M′. If ∣M′∣≥n|M'| \ge n∣M′∣≥n, it gives each agent exactly one good.
  4. Each eating stage is envy-free (proof of Theorem 2, p. 11).
  5. Lemma 1 (p. 11). In iteration t<⌈m/n⌉t < \lceil m/n \rceilt<⌈m/n⌉, each agent weakly prefers the good she receives to every good still unallocated.
  6. Lemma 2 (p. 11). Every allocation returned by Algorithm 1 is EF1.
  7. Algorithm 1 returns a randomized allocation (§3.1, p. 9).
  8. The expectation recursion (proof of Theorem 2, p. 11; Algorithm 3, p. 33). The implemented allocation is the first eating stage plus the weighted sum of the allocations implemented on the remaining goods.

Significance

Theorem 2 shows that ex-ante envy-freeness and ex-post EF1 are compatible for additive goods. Algorithm 1 is a constructive witness. The paper builds on this in two directions. A polynomial-time variant (§3.2) computes such a lottery. An impossibility result (Theorem 3) shows that adding ex-post fractional Pareto optimality is too much. The paper also extends the approach to bads and mixed manna (§7), where the natural extension guarantees only ex-post EF2.

The result is proved on paper. No formalization of probabilistic serial, the eating protocol, EF1 of randomized allocations or the rectangular Birkhoff–von Neumann theorem with floor and ceiling bounds is known to us in any proof assistant. Mathlib has the square doubly stochastic case (doublyStochastic_eq_convexHull_permMatrix), and Prove2Me has it as FamousTheorems.birkhoff_von_neumann. This mission asks for a machine-checked proof of Theorem 2 about an explicit definition of the eating protocol and of Algorithm 1.

Difficulty

The proof is short on paper but hides two pieces of continuous-time bookkeeping. First, the eating protocol is a process in time. The stage statements need a description of it in which "agent iii ate a positive share of ggg while g∗g^*g∗ was still available" can be stated and checked. Second, the ex-ante argument conditions on the history of earlier draws. Because the decomposition may have exponentially many branches, it has to be turned into an induction over the recursive lottery. The obvious shortcut, proving the result for "the" BvN decomposition of each stage, fails because the decomposition is not unique and the theorem must hold for all of them. Weight-zero terms of a decomposition can violate Lemma 1. They are never drawn, but any argument over the whole support must exclude them explicitly.

Formalization scope

Agents and goods are finite types with decidable equality. The good set MMM is the whole type of goods, and n≥1n \ge 1n≥1 is [Nonempty N]. Values are real and non-negative. All quantities are in R\mathbb RR, and ⌈m/n⌉\lceil m/n \rceil⌈m/n⌉ is Nat.ceil of a real quotient. Integral allocations are bundle maps N → Finset G, and lotteries are finite lists of (probability, allocation) pairs.

The eating protocol is a definition computed event by event. Between events every agent eats a fixed good chosen by a choice rule. The next event comes when the time budget ends or a good being eaten is exhausted. The minimum runs only over goods actually being eaten, so no division by zero occurs. ∣M∣+1|M| + 1∣M∣+1 events suffice, and later iterations change nothing. Algorithm 1 is a definition: a recursion over stages producing a lottery over stage histories.

The only free parameters of the goal are the tie-breaking rule, constrained to choose a most valuable available good, and the decomposition rule, constrained to satisfy Theorem 1's conclusions. The eating outcome, the lottery and its support are computed by the definitions, never assumed. A statement taking XtX^tXt or the lottery as free variables with envy-freeness as a hypothesis would trivialize the goal and is not this mission.

Conventions that differ from the page:

  • Ex-post EF1 is over allocations drawn with positive probability.
  • The without-loss-of-generality assumption that every good is valued positively by some agent is omitted. The proof does not use it, so the goal is stronger.
  • Theorem 1 is stated as existence, without its running-time clause.
  • Algorithm 1's printed line 2 initializes M0M^0M0, while the loop reads M1M^1M1. Its line 8 is a union over all stages inside the loop. Both are read as the intended algorithm, consistent with Algorithm 3 and with the proof of Theorem 2.
  • Histories are indexed from 000.

Not included are Proposition 1, the §3.2 polynomial-time variant, the SD-EF remark after Theorem 2, and §§4–7 (Theorem 3, the group-fairness results, the MNW characterization, bads). A proof of Theorem 1 in this rectangular form, with the floor and ceiling bounds, would be reusable well beyond this mission, for example in other fair-division and assignment results. So would a clean theory of the eating protocol. Contributions to either are welcome.

Selected references

  • R. Freeman, N. Shah, R. Vaish, Best of Both Worlds: Ex-Ante and Ex-Post Fairness in Resource Allocation, arXiv:2005.14122v1, 2020. https://arxiv.org/abs/2005.14122 (journal version with H. Aziz, Operations Research, 2024).
  • A. Bogomolnaia, H. Moulin, A New Solution to the Random Assignment Problem, Journal of Economic Theory 100(2), 2001. https://doi.org/10.1006/jeth.2000.2710
  • R. J. Lipton, E. Markakis, E. Mossel, A. Saberi, On Approximately Fair Allocations of Indivisible Goods, ACM EC 2004. https://doi.org/10.1145/988772.988792
  • E. Budish, The Combinatorial Assignment Problem: Approximate Competitive Equilibrium from Equal Incomes, Journal of Political Economy 119(6), 2011. https://doi.org/10.1086/664613
  • E. Budish, Y.-K. Che, F. Kojima, P. Milgrom, Designing Random Allocation Mechanisms: Theory and Applications, American Economic Review 103(2), 2013. https://doi.org/10.1257/aer.103.2.585
10 thms1 active userReviewed
ProbabilityStochastic Systems·Captain: mikedeng1

When Service Times Depend on Customers’ Delays: A Relationship Between Two Models of Dependence: Every Exogenous Service–Patience Dependence Has an Equivalent Endogenous DependenceResearch Paper

Motivation

In call centers, hospitals and other service systems, customers who wait too long leave: each has a patience time, and abandons the queue once the wait in queue exceeds it. Empirical work has found that service times and waiting are dependent. Customers who waited longer often need longer (or shorter) service, and those who were served are a self-selected group: they had enough patience to wait. Two models of this dependence are in use.

  • In the exogenous model, each customer arrives carrying a pair (service time, patience time) with a joint density. The dependence is a property of the customer population.
  • In the endogenous model, the service time a customer requires depends on the delay that customer actually experiences in queue, and the patience is independent of everything else.

The two models suggest different statistical procedures. Estimating an exogenous dependence means fitting a two-sided censored bivariate distribution. Estimating an endogenous one means a univariate patience estimate (Kaplan–Meier) plus a conditional service-time distribution per delay. Wu, Bassamboo & Perry (2022) prove that, as far as the law of the queueing dynamics is concerned, every exogenous dependence is an endogenous one. The simpler estimation procedure therefore never loses generality.

Setting

A service system has n≥1n \ge 1n≥1 identical servers. Customers arrive one at a time, at times κ1<κ2<⋯\kappa_1 < \kappa_2 < \cdotsκ1​<κ2​<⋯, and are served first come, first served. A customer who finds a free server starts service at once. Otherwise the customer waits, and abandons if the wait exceeds the customer's patience.

The offered wait ZiZ_iZi​ of customer iii is the time the customer would wait with infinite patience. With the sorting operator R\mathcal RR (which sorts a vector of Rn\mathbb R^nRn in increasing order), e1=(1,0,…,0)e_1 = (1,0,\dots,0)e1​=(1,0,…,0), 1=(1,…,1)\mathbf 1 = (1,\dots,1)1=(1,…,1) and interarrival times αi+1=κi+1−κi\alpha_{i+1} = \kappa_{i+1} - \kappa_iαi+1​=κi+1​−κi​, the offered-wait vectors follow the recursion

Zi+1=R([Zi+I{Zi(1)≤Ti}Si(Zi(1)) e1−αi+11]+),Zi=Zi(1).(2)\mathbf Z_{i+1} = \mathcal R\big([\mathbf Z_i + \mathbb I_{\{\mathbf Z_i(1) \le T_i\}} S_i(\mathbf Z_i(1))\,e_1 - \alpha_{i+1}\mathbf 1]^+\big), \qquad Z_i = \mathbf Z_i(1). \tag{2}Zi+1​=R([Zi​+I{Zi​(1)≤Ti​}​Si​(Zi​(1))e1​−αi+1​1]+),Zi​=Zi​(1).(2)

Here TiT_iTi​ is the patience of customer iii and Si(z)S_i(z)Si​(z) the service time it requires at offered wait zzz.

The queue is described by the process X(t)=(Y(t),Q(t),U(t),V(t),B(t))\mathbf X(t) = (Y(t), Q(t), U(t), V(t), B(t))X(t)=(Y(t),Q(t),U(t),V(t),B(t)), t≥0t \ge 0t≥0. Its components are:

  • Y(t)Y(t)Y(t) and Q(t)Q(t)Q(t), the numbers of customers in service and in queue;

  • U(t)∈R+nU(t) \in \mathbb R^n_+U(t)∈R+n​, which lists the idle servers as zeros, followed by the remaining service times in increasing order;

  • V(t)V(t)V(t), the elapsed waiting times of the customers in queue, in decreasing order and padded with zeros;

  • B(t)B(t)B(t), the time to the next arrival.

  • Exogenous dependence fexf^{ex}fex. The pairs (Siex,Tiex)(S^{ex}_i, T^{ex}_i)(Siex​,Tiex​) are i.i.d. with joint density fexf^{ex}fex on R+2\mathbb R_+^2R+2​, and independent of the arrivals and the initial state. The service time does not depend on the wait. The patience marginal is fTex(y)=∫0∞fex(x,y) dxf^{ex}_T(y) = \int_0^\infty f^{ex}(x,y)\,dxfTex​(y)=∫0∞​fex(x,y)dx.

  • Endogenous dependence (fTen,ψen)(f^{en}_T, \psi^{en})(fTen​,ψen). The patience times TienT^{en}_iTien​ are i.i.d. with density fTenf^{en}_TfTen​. Customer iii has a virtual service-time process z↦Sien(z)z \mapsto S^{en}_i(z)z↦Sien​(z) whose value at zzz has density ψen(⋅,z)\psi^{en}(\cdot, z)ψen(⋅,z). These are independent across customers, of the patience and of the arrivals.

Lex(A,n;fex)\mathcal L^{ex}(A, n; f^{ex})Lex(A,n;fex) and Len(A,n;fTen,ψen)\mathcal L^{en}(A, n; f^{en}_T, \psi^{en})Len(A,n;fTen​,ψen) denote the laws of the process X\mathbf XX. The two dependencies are equivalent (Definition 1) if these laws agree for every arrival process AAA and every nnn, whenever the initial states agree in law.

Formalization targets

Goal: Theorem 1 (p. 3348)

For every joint density fexf^{ex}fex on R+2\mathbb R_+^2R+2​, the pair

fTen(z)=fTex(z),ψen(x,z)=∫z∞fex(x,y) dy∫z∞fTex(y) dy(1)f^{en}_T(z) = f^{ex}_T(z), \qquad \psi^{en}(x,z) = \frac{\int_z^\infty f^{ex}(x,y)\,dy}{\int_z^\infty f^{ex}_T(y)\,dy} \tag{1}fTen​(z)=fTex​(z),ψen(x,z)=∫z∞​fTex​(y)dy∫z∞​fex(x,y)dy​(1)

is an endogenous dependence (a patience density and a conditional service-time density), and it is equivalent to fexf^{ex}fex:

Lex(A,n;fex)=Len(A,n;fTen,ψen)for every n≥1 and every arrival process A.\mathcal L^{ex}(A, n; f^{ex}) = \mathcal L^{en}(A, n; f^{en}_T, \psi^{en}) \quad \text{for every } n \ge 1 \text{ and every arrival process } A.Lex(A,n;fex)=Len(A,n;fTen​,ψen)for every n≥1 and every arrival process A.

Milestones

  1. Formula (1) defines densities.
  2. Lemma 1, sufficiency: the law of Xen\mathbf X^{en}Xen is determined by (fT,ψ)(f_T, \psi)(fT​,ψ), the arrivals and nnn.
  3. The pathwise coupling: two systems in which every customer has the same outcome at its offered wait have the same offered waits and the same path X\mathbf XX.
  4. The display after the proof: Ψˉen(x,z)FˉTen(z)=∫z∞∫x∞fex(u,v) du dv\bar\Psi^{en}(x,z)\bar F^{en}_T(z) = \int_z^\infty\int_x^\infty f^{ex}(u,v)\,du\,dvΨˉen(x,z)FˉTen​(z)=∫z∞​∫x∞​fex(u,v)dudv, decreasing in xxx and in zzz.
  5. The coupling kernel: the pair generated by steps (ii)–(iii) of the proof has joint density fexf^{ex}fex, whatever the offered wait.
  6. The law of Xex\mathbf X^{ex}Xex is determined by fexf^{ex}fex, the arrivals and nnn.

Significance

The theorem shows that the systems with exogenous dependence form a subclass of those with endogenous dependence. Two consequences follow:

  • Any statement about the law of a queue with exogenous dependence (delay distribution, abandonment probability, occupancy) is a statement about some endogenous system.
  • Statistical analysis can be carried out with the endogenous model, a Kaplan–Meier patience estimate and binned conditional service-time distributions, without assuming the dependence is truly endogenous.

The formula for ψen\psi^{en}ψen also explains the self-selection effect. It is the service-time density of a customer, conditioned on having patience exceeding the offered wait.

The paper's proof omits the verification step "for brevity". It also states the law-determination facts (milestones 2 and 6) without proof, and treats an initially empty system. A formal proof would supply the per-customer density computation, the measure-theoretic construction of the coupled system, and the measurability of the path map.

Difficulty

The obvious argument would compare the two systems customer by customer with the same random draws. It fails: in the exogenous system the service time and the patience are drawn together before the customer meets the queue, while in the endogenous system the service time is drawn after the offered wait is known. Equality holds only in law, through a coupling built on a new probability space, and the coupled pairs must be shown to be i.i.d. with density fexf^{ex}fex even though they are constructed from the state of the queue.

Two technical points are not visible in the paper. First, the law of the whole path t↦X(t)t \mapsto \mathbf X(t)t↦X(t) must be shown to be determined by the primitives. This requires measurability of a map built from sorting, counting and an infimum over arrival indices. Second, the conditional densities in (1) and in steps (ii)–(iii) are ratios whose denominators can vanish, so null sets need care.

Formalization scope

Everything lives in the namespace DelayDep.Equiv. The conventions are:

  • Paths and time. Customers are numbered from 000 (customer iii is the paper's customer i+1i+1i+1). Time is NNReal. The law of X\mathbf XX is the image measure on NNReal → State n with the product σ-algebra, which is the σ-algebra of the finite-dimensional distributions.
  • Densities. Densities are real-valued, nonnegative and measurable on [0,∞)[0,\infty)[0,∞) or [0,∞)2[0,\infty)^2[0,∞)2, and enter measures through withDensity.
  • Ties. Recursion (2)'s convention is used: a customer is served iff offered wait ≤\le≤ patience. Ties are null events.
  • Arrivals and initial state. The arrivals are any measurable positive interarrival times whose sums tend to infinity. The initial state is an empty queue with arbitrary nonnegative residual service times. "Any arrival process AAA with Xex(0)=dXen(0)\mathbf X^{ex}(0) \overset{d}{=} \mathbf X^{en}(0)Xex(0)=dXen(0)" becomes "any joint law of interarrival times and initial residual service times, the same in both systems".
  • Assumption 1. Assumption 1 (differentiability of Ψen\Psi^{en}Ψen) is not imposed. The ψen\psi^{en}ψen of (1) need not satisfy it, and Theorem 1 does not use it.
  • Tail guard. The endogenous service-time law is prescribed only at offered waits with positive patience tail. Elsewhere a customer abandons almost surely, and (1) is 0/00/00/0.

Every law equality includes the almost-everywhere measurability of both path maps as part of its conclusion. Without it, the equality could hold trivially through Lean's convention that the image of a measure under a non-measurable map is 000. For the same reason, the goal also asserts that (1) defines densities, so that the endogenous systems it quantifies over are not vacuous.

Out of scope are Theorem 2, the necessity half of Lemma 1, Example 1, Theorem A.1 and the estimation section §4. Useful contributions are:

  • a general measurability lemma for sample-path functionals of this kind;
  • the kernel computation of milestone 5;
  • the construction of the coupled system;
  • the law-determination lemmas, which apply to any queue driven by i.i.d. customer primitives.

Selected references

  • C. A. Wu, A. Bassamboo, O. Perry, When service times depend on customers' delays: A relationship between two models of dependence, Operations Research 70(6):3345–3354, 2022. https://doi.org/10.1287/opre.2021.2179
9 thms1 active userReviewed
ProbabilityStochastic Systems·Captain: mikedeng1

Technical Note–A Simple Heuristic Policy for Stochastic Distribution Inventory Systems with Fixed Shipment Costs: An Upper Bound on the Long-Run Average Cost of Modified Echelon (r,Q) PoliciesResearch Paper

Motivation

A one-warehouse, multi-retailer (OWMR) distribution system is the basic model of a supply chain in which a central warehouse buys from an outside supplier and ships to several retail locations. When every shipment carries a fixed cost and a positive lead time, the optimal policy of the stochastic system is unknown, and the classical approach is to search over a structured class of policies and bound its cost. Exact evaluation of even a simple policy is hard, because the warehouse can run out of stock and then the retailers' replenishments interact.

Zhu, Chen, Hu and Yang (SSRN 2921743; Operations Research, 2021, doi:10.1287/opre.2021.2101) study modified echelon (r,Q)(r,Q)(r,Q) policies for continuous-review OWMR systems with Poisson demand. Their central tool is an explicit upper bound on the long-run average cost of any such policy, written as the retailers' single-stage (r,Q)(r,Q)(r,Q) costs, a warehouse term that depends on the retailer policies, and a fixed-cost allowance. Their heuristic policy, its worst-case performance guarantee and its asymptotic optimality (Theorems 2–4 of the paper) are all built on this bound. The single-stage cost functions follow Zheng (1992) and the echelon cost accounting follows Chen and Zheng (1994).

Setting

There are N≥1N\ge1N≥1 retailers, indexed i=1,…,Ni=1,\dots,Ni=1,…,N, and one warehouse, indexed 000. Retailer iii faces Poisson demand of rate λi>0\lambda_i>0λi​>0, independent across retailers, so the warehouse sees Poisson demand of rate λ0=∑iλi\lambda_0=\sum_i\lambda_iλ0​=∑i​λi​. Installation iii has a constant lead time Li>0L_i>0Li​>0, a fixed cost Ki>0K_i>0Ki​>0 per shipment, and an echelon holding cost rate hi>0h_i>0hi​>0; unmet retailer demand is backlogged at rate pi>0p_i>0pi​>0. Write Kˉ=max⁡iKi\bar K=\max_iK_iKˉ=maxi​Ki​.

A modified echelon (r,Q)(r,Q)(r,Q) policy has reorder levels rir_iri​ and order quantities Qi>0Q_i>0Qi​>0 for every installation. When the echelon inventory position IPiIP_iIPi​ of installation iii is at or below rir_iri​ and the upstream installation has stock on hand, a shipment raises IPiIP_iIPi​ as close as possible to ri+Qir_i+Q_iri​+Qi​.

With Di∼Poisson(λiLi)D_i\sim\mathrm{Poisson}(\lambda_iL_i)Di​∼Poisson(λi​Li​) the lead-time demand of retailer iii, the single-stage holding-and-backlog rate and the single-stage (r,Q)(r,Q)(r,Q) cost are

Gi(y)=E[hi(y−Di)++(h0+pi)(y−Di)−],Ci(ri,Qi)=λiKi+∫riri+QiGi(y) dyQi.G_i(y)=\mathbb E\big[h_i(y-D_i)^++(h_0+p_i)(y-D_i)^-\big],\qquad C_i(r_i,Q_i)=\frac{\lambda_iK_i+\int_{r_i}^{r_i+Q_i}G_i(y)\,dy}{Q_i}.Gi​(y)=E[hi​(y−Di​)++(h0​+pi​)(y−Di​)−],Ci​(ri​,Qi​)=Qi​λi​Ki​+∫ri​ri​+Qi​​Gi​(y)dy​.

With Gi(wi)=sup⁡ri<z≤ri+QiGi(z)G_i(w_i)=\sup_{r_i<z\le r_i+Q_i}G_i(z)Gi​(wi​)=supri​<z≤ri​+Qi​​Gi​(z), Mj=max⁡{Gj(wj),Cj(rj,Qj)}M_j=\max\{G_j(w_j),C_j(r_j,Q_j)\}Mj​=max{Gj​(wj​),Cj​(rj​,Qj​)} and Si=∑j≠i(rj+Qj)S_i=\sum_{j\ne i}(r_j+Q_j)Si​=∑j=i​(rj​+Qj​), the envelope of Lemma 3 at warehouse echelon level xxx is

Γˉi(x)=∑j≠iMj+{Mix−Si>ri,max⁡{Gi(x−Si),Mi}otherwise,Γˉ(x)=max⁡iΓˉi(x).\bar\Gamma_i(x)=\sum_{j\ne i}M_j+\begin{cases}M_i & x-S_i>r_i,\\ \max\{G_i(x-S_i),M_i\} & \text{otherwise},\end{cases}\qquad \bar\Gamma(x)=\max_i\bar\Gamma_i(x).Γˉi​(x)=j=i∑​Mj​+{Mi​max{Gi​(x−Si​),Mi​}​x−Si​>ri​,otherwise,​Γˉ(x)=imax​Γˉi​(x).

With D0∼Poisson(λ0L0)D_0\sim\mathrm{Poisson}(\lambda_0L_0)D0​∼Poisson(λ0​L0​), set G^=Γˉ−∑iCi\hat G=\bar\Gamma-\sum_iC_iG^=Γˉ−∑i​Ci​, Λ0(y)=E[h0(y−D0)+G^(y−D0)]\Lambda_0(y)=\mathbb E[h_0(y-D_0)+\hat G(y-D_0)]Λ0​(y)=E[h0​(y−D0​)+G^(y−D0​)] and

C^0(r0,Q0)=1Q0[λ0K0+∫r0r0+Q0Λ0(y) dy].\hat C_0(r_0,Q_0)=\frac1{Q_0}\Big[\lambda_0K_0+\int_{r_0}^{r_0+Q_0}\Lambda_0(y)\,dy\Big].C^0​(r0​,Q0​)=Q0​1​[λ0​K0​+∫r0​r0​+Q0​​Λ0​(y)dy].

The long-run average system cost is (1) of the paper,

C(r,Q)=1Q0[λ0K0+∫r0r0+Q0E[h0(y−D0)+Ω(y)] dy],C(r,Q)=\frac1{Q_0}\Big[\lambda_0K_0+\int_{r_0}^{r_0+Q_0}\mathbb E\big[h_0(y-D_0)+\Omega(y)\big]\,dy\Big],C(r,Q)=Q0​1​[λ0​K0​+∫r0​r0​+Q0​​E[h0​(y−D0​)+Ω(y)]dy],

where Ω(y)\Omega(y)Ω(y) is the expected cost rate of all retailers when the warehouse inventory position is yyy.

Formalization targets

Goal: Theorem 1 (p. 13)

C(r,Q)≤CU(r,Q)≡∑i=1NCi(ri,Qi)+C^0(r0,Q0)+λ0KˉQ0.C(r,Q)\le C^U(r,Q)\equiv\sum_{i=1}^NC_i(r_i,Q_i)+\hat C_0(r_0,Q_0)+\frac{\lambda_0\bar K}{Q_0}.C(r,Q)≤CU(r,Q)≡i=1∑N​Ci​(ri​,Qi​)+C^0​(r0​,Q0​)+Q0​λ0​Kˉ​.

It is stated for every policy, with no assumption on the parameters beyond Qi,Q0>0Q_i,Q_0>0Qi​,Q0​>0.

Milestones

  1. App. A (3): the per-retailer bound: the charge of retailer iii is at most the own-retailer term of Γˉi(x)\bar\Gamma_i(x)Γˉi​(x).
  2. Lemma 3: the total retailer charge at a feasible state is at most Γˉ(x)\bar\Gamma(x)Γˉ(x).
  3. App. A (8): Γ^(x)≤Γˉ(x)=G^(x)+∑iCi\hat\Gamma(x)\le\bar\Gamma(x)=\hat G(x)+\sum_iC_iΓ^(x)≤Γˉ(x)=G^(x)+∑i​Ci​.
  4. App. A (9): Γ0(y)≤Λ0(y)+∑iCi\Gamma_0(y)\le\Lambda_0(y)+\sum_iC_iΓ0​(y)≤Λ0​(y)+∑i​Ci​.
  5. App. A (10): 1Q0[λ0K0+∫r0r0+Q0Γ0]≤C^0(r0,Q0)+∑iCi\frac1{Q_0}[\lambda_0K_0+\int_{r_0}^{r_0+Q_0}\Gamma_0]\le\hat C_0(r_0,Q_0)+\sum_iC_iQ0​1​[λ0​K0​+∫r0​r0​+Q0​​Γ0​]≤C^0​(r0​,Q0​)+∑i​Ci​.

Significance

The bound replaces the cost of an (N+1)(N+1)(N+1)-installation stochastic system, which has no closed form, by computable retailer costs and a warehouse term. This is what makes the paper's heuristic policy computable and what its performance guarantees compare against: each retailer uses its single-stage optimal (ri∗,Qi∗)(r_i^*,Q_i^*)(ri∗​,Qi∗​), and the warehouse minimizes the remaining term. Bounds of this decomposition type are a standard device in multi-echelon inventory theory, and Lemma 3's envelope, which charges each retailer at the lowest inventory position the others can leave it, is the step specific to the distribution setting.

The result is proved in the paper's online appendix. It has not been formalized. A formal development adds to the platform a single-stage (r,Q)(r,Q)(r,Q) cost with Poisson lead-time demand, the echelon envelope of Lemma 3 for general NNN, and an explicit OWMR cost decomposition, all of which are reusable by later work on distribution systems, including the paper's own performance bounds.

Difficulty

The core is Lemma 3 for general NNN. Bounding each retailer separately by its worst case max⁡{Gi(wi),Ci}\max\{G_i(w_i),C_i\}max{Gi​(wi​),Ci​} fails, because a retailer below its reorder level can have a much larger cost rate Gi(IPi)G_i(IP_i)Gi​(IPi​); the bound must use that the inventory positions are coupled through the warehouse echelon level. The printed proof is an induction on NNN whose step needs convexity of a sum in which one term is only bounded by a convex function, so the order of bounding and convexity arguments matters. Convexity and continuity of GiG_iGi​ (a Poisson mixture of piecewise-linear functions) must be established, and every expectation against the Poisson law needs an integrability argument, since the integrands grow linearly.

Formalization scope

All declarations live in the namespace DistribRQ.UpperBound. Retailers are indexed by Fin N with [NeZero N]; the warehouse data L0,K0,h0,r0,Q0L_0,K_0,h_0,r_0,Q_0L0​,K0​,h0​,r0​,Q0​ are separate. The primitives are bundled in Primitives N (all rates and costs positive) and the policy in Policy N (real ri,r0r_i,r_0ri​,r0​; Qi,Q0>0Q_i,Q_0>0Qi​,Q0​>0), following the paper's continuous approximation (p. 7). Poisson expectations are Bochner integrals over N\mathbb NN against poissonMeasure, evaluated at real arguments y−ny-ny−n. Gi(wi)G_i(w_i)Gi​(wi​) is the supremum of GiG_iGi​ over (ri,ri+Qi](r_i,r_i+Q_i](ri​,ri​+Qi​], since the paper's arg⁡max⁡\arg\maxargmax need not be attained on a half-open interval.

The sample-path model (cycles, shipment intervals, the paper's Lemmas 1–2) has no substrate, so the system is encoded as the paper's proof uses it:

  • C(r,Q)C(r,Q)C(r,Q) is (1), with Ω\OmegaΩ any interval-integrable function satisfying the bound of p. 12, Ω(y)≤λ0Kˉ/Q0+E[Γ^(y−D0)]\Omega(y)\le\lambda_0\bar K/Q_0+\mathbb E[\hat\Gamma(y-D_0)]Ω(y)≤λ0​Kˉ/Q0​+E[Γ^(y−D0​)]; this bound is where Lemmas 1–2 enter and is a hypothesis.
  • Γ^(x)=∫∑ici(zi) κ(x)(dz)\hat\Gamma(x)=\int\sum_ic_i(z_i)\,\kappa(x)(dz)Γ^(x)=∫∑i​ci​(zi​)κ(x)(dz) for any Markov kernel κ\kappaκ concentrated on the feasible states of p. A1 (zj≤rj+Qjz_j\le r_j+Q_jzj​≤rj​+Qj​; ∑jzj≤x\sum_jz_j\le x∑j​zj​≤x; and ∑jzj<x\sum_jz_j<x∑j​zj​<x forces zi>riz_i>r_izi​>ri​ for all iii), and any measurable, nonnegative admissible charges (ci(z)≤max⁡{Gi(z),Ci}c_i(z)\le\max\{G_i(z),C_i\}ci​(z)≤max{Gi​(z),Ci​} for z>riz>r_iz>ri​, ci(z)≤Gi(z)c_i(z)\le G_i(z)ci​(z)≤Gi​(z) for z≤riz\le r_iz≤ri​).
  • The theorem holds for every such κ\kappaκ, ccc and Ω\OmegaΩ, so the real system is an instance.

Γ^\hat\GammaΓ^ is never a free function bounded by Γˉ\bar\GammaΓˉ, the feasible set keeps the clause for positive warehouse stock, and GiG_iGi​, CiC_iCi​, Λ0\Lambda_0Λ0​ are computed from the Poisson law rather than taken as parameters; formalizations that relax any of these, or that rely on a junk value of sSup, an integral or a division, trivialize the statement and are not accepted.

Contributions welcome: convexity and continuity of GiG_iGi​, Poisson integrability lemmas, Lemma 3 for general NNN, and the averaging steps of App. A (8)–(10).

Selected references

  • H. Zhu, Y. (F.) Chen, M. Hu, Y. Yang, Technical Note—A Simple Heuristic Policy for Stochastic Distribution Inventory Systems with Fixed Shipment Costs, Operations Research, 2021. SSRN copy of 2020-09-08: https://ssrn.com/abstract=2921743; doi:10.1287/opre.2021.2101.
  • Y.-S. Zheng, On properties of stochastic inventory systems, Management Science 38(1), 1992. https://doi.org/10.1287/mnsc.38.1.87
  • F. Chen, Y.-S. Zheng, Lower bounds for multi-echelon stochastic inventory systems, Management Science 40(11), 1994. https://doi.org/10.1287/mnsc.40.11.1426
7 thms1 active userReviewed
Control TheoryStochastic Systems·Captain: mikedeng1

New Venture Creation: A Drift-Variance Diffusion Control Model: The Optimal Policy Is at Most a Two-Interval PolicyResearch Paper

Motivation

An entrepreneur developing a venture repeatedly chooses between strategies that differ in speed, risk and cost: a safe, slow approach or a fast, risky one. Wang and Zenios model this as a drift–variance diffusion control problem (SSRN 3667132; published in Operations Research, 2021, doi:10.1287/opre.2021.2171): the controller chooses not only the drift of a Brownian state but also its variance, and pays a running cost for each choice. Earlier models in which a decision maker switches controls as a diffusion evolves, Radner and Shepp (1996) on corporate risk choice and Harrison and Sunar (2015) on multiple modes of learning, describe their optimal policies through numerical procedures; here the controller chooses drift and variance together, the trade-off between moving fast and moving reliably is the object of study, and the optimal policy is in closed form.

The paper's motivating example (p. 10) suggests that the optimal policy might be complicated: use the safe control near 000 to avoid failure, the cheaper risky control in the middle, and the safe control again near UUU to secure the reward. The paper's main result says this never happens.

Setting

There are two controls i∈{1,2}i \in \{1,2\}i∈{1,2}. Control iii gives drift μi>0\mu_i > 0μi​>0, volatility σi>0\sigma_i > 0σi​>0 (variance σi2\sigma_i^2σi2​) and costs hi>0h_i > 0hi​>0 per unit time. The state follows

dX(t)=μ(A(t)) dt+σ(A(t)) dB(t)dX(t) = \mu(A(t))\,dt + \sigma(A(t))\,dB(t)dX(t)=μ(A(t))dt+σ(A(t))dB(t)

on [0,U][0,U][0,U], U>0U > 0U>0, and stops on hitting 000 (no payment) or UUU (terminal cost M<0M < 0M<0, i.e. a reward ∣M∣|M|∣M∣). The entrepreneur minimizes the expected total cost.

Write vi=2μi/σi2v_i = 2\mu_i/\sigma_i^2vi​=2μi​/σi2​ for the effective drift, hi∗=−Mμi/Uh_i^* = -M\mu_i/Uhi∗​=−Mμi​/U for the critical cost (the cost at which control iii alone breaks even with no noise), and Δhi=hi−hi∗\Delta h_i = h_i - h_i^*Δhi​=hi​−hi∗​. As in the paper (p. 12), v1<v2v_1 < v_2v1​<v2​ throughout.

A kkk-interval policy (Definition 1) cuts [0,U)[0,U)[0,U) at 0=S0<S1<⋯<Sk=U0 = S_0 < S_1 < \dots < S_k = U0=S0​<S1​<⋯<Sk​=U and uses control iji_jij​ on [Sj−1,Sj)[S_{j-1},S_j)[Sj−1​,Sj​), with ij≠ij+1i_j \ne i_{j+1}ij​=ij+1​. Its value is the function WWW that is continuously differentiable on [0,U][0,U][0,U], satisfies W(0)=0W(0) = 0W(0)=0, W(U)=MW(U) = MW(U)=M, and on each piece solves

hij+μijW′(x)+σij22W′′(x)=0.h_{i_j} + \mu_{i_j}W'(x) + \frac{\sigma_{i_j}^2}{2}W''(x) = 0 .hij​​+μij​​W′(x)+2σij​2​​W′′(x)=0.

A policy is optimal if its value lies below the value of every interval policy on [0,U][0,U][0,U].

Formalization targets

Goal: Theorem 1 (p. 10)

With v1<v2v_1 < v_2v1​<v2​:

∃ π optimal with k(π)≤2,and(h1,h2)≠(h1∗,h2∗)  ⟹  ∀ π optimal, k(π)≤2.\exists\,\pi \text{ optimal with } k(\pi) \le 2, \qquad\text{and}\qquad (h_1,h_2) \ne (h_1^*,h_2^*) \implies \forall\,\pi \text{ optimal},\ k(\pi) \le 2 .∃π optimal with k(π)≤2,and(h1​,h2​)=(h1∗​,h2∗​)⟹∀π optimal, k(π)≤2.

In words: an optimal policy either uses one control throughout or switches exactly once, and away from the critical point no optimal policy switches more than once.

Milestones

  • Theorem 2 (pp. 12–13): four explicit slopes d1L,d1U,d2L,d2Ud_1^L, d_1^U, d_2^L, d_2^Ud1L​,d1U​,d2L​,d2U​ cut the (Δh1,Δh2)(\Delta h_1, \Delta h_2)(Δh1​,Δh2​)-plane into regions (7)–(10) on which "always 1", "always 2", "1 then 2" and "2 then 1" are optimal.
  • Exhaustiveness (p. 13): one of (7)–(10) always holds.
  • Lemma 1 (p. 14): d1L<d2L<σ22/σ12<d1U<d2U<μ2/μ1d_1^L < d_2^L < \sigma_2^2/\sigma_1^2 < d_1^U < d_2^U < \mu_2/\mu_1d1L​<d2L​<σ22​/σ12​<d1U​<d2U​<μ2​/μ1​ and d2L>0d_2^L > 0d2L​>0, with the elementary Lemmas 20–21.
  • Lemma 12, analytic form: a C1C^1C1, piecewise C2C^2C2 function with the boundary values that satisfies the HJB inequality lies below every policy value.
  • The ingredients of Theorem 2 (a) and (c): the static value (54) and inequality (55); the two-interval value (74)/(87); Lemmas 4, 5, 6, 7, 9, 10 and 22 on the test functions f1,f2,F1,F2f_1, f_2, F_1, F_2f1​,f2​,F1​,F2​.

Significance

The result reduces an infinite-dimensional control problem to a two-parameter family of candidates (which control first, and where to switch), and Theorem 2 identifies the optimum explicitly. Downstream, the paper reads managerial conclusions off the region picture: when both controls are inefficient, when to begin with the safe control and when with the risky one (§4, §8).

The theorem is proved in the paper; nothing here is open. No machine-checked version of any part of it is known. The mission produces a checked version of the analytic core: the closed-form value functions, the sign analysis of the switching functions, a one-dimensional comparison principle (the analytic content of the verification lemma), and the correction of the uniqueness half of Theorem 1 at the critical cost point.

Difficulty

The existence half is a long chain of explicit but delicate calculus. The slopes involve U/(1−e−viU)U/(1 - e^{-v_iU})U/(1−e−vi​U) terms whose ordering (Lemma 1) is not visible from the formulas, and the switching point S∗S^*S∗ of the two-interval policy is the root of a function f1f_1f1​ built from the derivatives of a piecewise exponential solution whose coefficients are ratios of exponential sums. Showing that F1,F2F_1, F_2F1​,F2​ keep their signs away from S∗S^*S∗ requires the sign facts of Lemma 22 in exactly the right region.

The uniqueness half looks immediate (two optimal policies have the same value, which cannot solve two different ODEs on an interval) but is false as printed: at (h1,h2)=(h1∗,h2∗)(h_1,h_2) = (h_1^*,h_2^*)(h1​,h2​)=(h1∗​,h2∗​) the affine function Mx/UMx/UMx/U solves both ODEs, and every interval policy is optimal. The correct statement must exclude that point, and its proof must show that, elsewhere, the optimal value is affine on no interval.

Formalization scope

Mathlib has Brownian motion but no stochastic integral or Itô formula, so the problem is stated analytically, in the objects the paper's proofs manipulate. A policy's value is the unique solution of its ODE boundary-value system with C1C^1C1 pasting; the identity with the expected cost (Dynkin's formula, the second half of Lemma 12) is not formalized, and optimality is over interval policies, to which the paper restricts (p. 8). Conventions: L=0L = 0L=0, m=0m = 0m=0, M<0M < 0M<0, minimization; σi\sigma_iσi​ is the volatility and σi2\sigma_i^2σi2​ enters the ODE; v1<v2v_1 < v_2v1​<v2​ is an explicit hypothesis; controls are a two-constructor type, policies are structures carrying kkk, the controls, strictly increasing breakpoints from 000 to UUU and distinct adjacent controls.

Corrections relative to the page: Theorem 1's uniqueness half excludes Δh1=Δh2=0\Delta h_1 = \Delta h_2 = 0Δh1​=Δh2​=0; Lemma 22 (204) is stated strictly; Lemmas 9–10 carry the hypotheses of condition (9) that their proofs use; the exclusivity claim of p. 13 is dropped (only exhaustiveness is stated); (73)'s σ12\sigma_1^2σ12​ in the control-2 equation is read as σ22\sigma_2^2σ22​.

The encoding cannot be trivialized: the value fixes both boundary values and C1C^1C1 gluing, so no policy has arbitrarily low "values"; optimality quantifies over policies of every length and all their values at every point of [0,U][0,U][0,U]; and the policy structure forbids empty pieces and repeated adjacent controls, so kkk is the true number of intervals.

Needed infrastructure: elementary real analysis (exponential inequalities, derivatives of explicit functions, the intermediate value theorem), a one-dimensional maximum principle for μW′+σ22W′′≥0\mu W' + \frac{\sigma^2}{2}W'' \ge 0μW′+2σ2​W′′≥0 with piecewise-constant coefficients, and linear independence of 1,x,e−v1x,e−v2x1, x, e^{-v_1x}, e^{-v_2x}1,x,e−v1​x,e−v2​x. The comparison principle and the interval-policy layer are reusable for other one-dimensional switching problems. Proofs of any milestone are welcome, as is a formal Dynkin link once Mathlib has stochastic integration.

Selected references

  • Z. Wang, S. Zenios, New Venture Creation: A Drift-Variance Diffusion Control Model, SSRN working paper 3667132; Operations Research, 2021. https://ssrn.com/abstract=3667132, https://doi.org/10.1287/opre.2021.2171
  • R. Radner, L. Shepp, Risk vs. profit potential: A model for corporate strategy, Journal of Economic Dynamics and Control 20(8), 1373–1393, 1996 (as cited on p. 33 of the source).
  • J. M. Harrison, N. Sunar, Investment timing with incomplete information and multiple means of learning, Operations Research 63(2), 442–457, 2015 (as cited on p. 33 of the source).
19 thms1 active userReviewed
ProbabilityStatistics·Captain: mikedeng1

Optimal Sequential Multi-class Diagnosis: Under Exponential Tilting, Every Reachable Posterior Belief Is a Closed-Form Function of an r-Dimensional Diagnostic Statistic, r = rank(B)Research Paper

Motivation

In sequential multi-class diagnosis, a decision maker observes a system repeatedly, then chooses when to stop and which hidden class to report. Each new observation can improve the diagnosis, but waiting has a cost. A direct dynamic program uses the vector of posterior class probabilities as its state. With N+1N+1N+1 possible classes, that state lies in an NNN-dimensional probability simplex. Wang's model asks when the observed sequence can instead be summarized by fewer coordinates, so that the same posterior can be recovered without retaining a full belief vector as an independent state Wang, §2–§3.

The paper considers class distributions related by exponential tilting. This relationship appears in common one-dimensional distribution families, including the normal, binomial, Poisson, and gamma examples collected in its Table 1 Wang, p. 8. The result targeted here is the belief reconstruction identity of Theorem 1. Later parts of the paper use the reduced statistic in optimality equations and study stopping regions; those results are outside this mission Wang, pp. 12–18.

Setting

The hidden class is j∈{0,…,N}j\in\{0,\ldots,N\}j∈{0,…,N}. Before observing data, the decision maker has a prior belief θ=(θ0,…,θN)\theta=(\theta_0,\ldots,\theta_N)θ=(θ0​,…,θN​): every component is nonnegative and their sum is one. Conditional on class jjj, observations y1,…,yty_1,\ldots,y_ty1​,…,yt​ are independent and have the same density fjf_jfj​ on the observation space Y\mathcal YY. Class 000 supplies a baseline density f0f_0f0​. The one-step Bayesian operator B\mathcal BB multiplies each current class weight by that class's likelihood at the new observation and divides by the sum of all weighted likelihoods. Starting at Π0=θ\Pi_0=\thetaΠ0​=θ, the belief Πt\Pi_tΠt​ is obtained by applying B\mathcal BB once for each observation in arrival order Wang, pp. 6–7.

In the exponential tilting model, a transformed observation h(y)∈Rph(y)\in\mathbb R^ph(y)∈Rp and a parameter βj∈Rp\beta_j\in\mathbb R^pβj​∈Rp describe each nonbaseline density:

fj(y)=f0(y)exp⁡{βj0+βjTh(y)},βj0=−log⁡∫Yexp⁡{βjTh(y)}f0(y) dy.f_j(y)=f_0(y)\exp\{\beta_{j0}+\beta_j^{\mathsf T}h(y)\},\qquad \beta_{j0}=-\log\int_{\mathcal Y}\exp\{\beta_j^{\mathsf T}h(y)\}f_0(y)\,dy.fj​(y)=f0​(y)exp{βj0​+βjT​h(y)},βj0​=−log∫Y​exp{βjT​h(y)}f0​(y)dy.

Here j=1,…,Nj=1,\ldots,Nj=1,…,N. Place the βjT\beta_j^{\mathsf T}βjT​ in the rows of the diagnostic parameter matrix B\mathbb BB. Its rank decomposition is B=RP\mathbb B=\mathbb R\mathbb PB=RP, where R\mathbb RR has NNN rows and rrr linearly independent columns, while P\mathbb PP has rrr linearly independent rows and ppp columns. Thus r=rank⁡(B)r=\operatorname{rank}(\mathbb B)r=rank(B). The paper calls R\mathbb RR the reconstruction matrix and P\mathbb PP the projection matrix Wang, Definition 1, p. 8.

Formalization targets

Theorem 1: posterior reconstruction

For t≥1t\ge1t≥1, define the diagnostic statistic xt=P(∑m=1th(ym))/t∈Rrx_t=\mathbb P\big(\sum_{m=1}^t h(y_m)\big)/t\in\mathbb R^rxt​=P(∑m=1t​h(ym​))/t∈Rr. For classes j=1,…,Nj=1,\ldots,Nj=1,…,N, the target identity is

Πt(y1,…,yt,θ)j=Tjt(xt,θ)=θjexp⁡{t(βj0+ejRxt)}θ0+∑i=1Nθiexp⁡{t(βi0+eiRxt)}.\Pi_t(y_1,\ldots,y_t,\theta)_j =\mathcal T_j^t(x_t,\theta) =\frac{\theta_j\exp\{t(\beta_{j0}+e_j\mathbb R x_t)\}} {\theta_0+\sum_{i=1}^{N}\theta_i\exp\{t(\beta_{i0}+e_i\mathbb R x_t)\}}.Πt​(y1​,…,yt​,θ)j​=Tjt​(xt​,θ)=θ0​+∑i=1N​θi​exp{t(βi0​+ei​Rxt​)}θj​exp{t(βj0​+ej​Rxt​)}​.

The baseline coordinate is T0t=1−∑j=1NTjt\mathcal T_0^t=1-\sum_{j=1}^{N}\mathcal T_j^tT0t​=1−∑j=1N​Tjt​. The normalizer appearing in Tt\mathcal T^tTt is evaluated from the row eiRPe_i\mathbb R\mathbb Pei​RP, as in equation (3), and the theorem identifies it with the model normalizer through B=RP\mathbb B=\mathbb R\mathbb PB=RP. The target also identifies the set of posteriors reachable after exactly ttt observations with Tt(Ωt,θ)\mathcal T^t(\Omega_t,\theta)Tt(Ωt​,θ), where Ωt\Omega_tΩt​ is the set of attainable diagnostic statistics. These statements give a precise reading of the paper's claim that reachable beliefs admit an rrr-coordinate representation Wang, Theorem 1, p. 9; Definition 2, p. 10.

The goal additionally states that the map x↦Tt(x,θ)x\mapsto\mathcal T^t(x,\theta)x↦Tt(x,θ) is injective on Rr\mathbb R^rRr when every prior component is positive. This is an explicit, qualified reading of the paper's phrase that the dimensionality is equal to rrr. For discrete observations, the reachable set Ωt\Omega_tΩt​ can be finite, so the phrase cannot always mean its geometric dimension is literally rrr. The five milestones follow the source's Bayes operator and the displayed identities in its proof of Theorem 1 Wang, p. ec1.

Significance

The formula determines every posterior component from rrr real coordinates, the prior, and the observation count. When rrr is smaller than NNN, this supplies the representation used by the paper to reformulate a diagnosis problem in a smaller state space. The paper's examples show cases with one or two diagnostic coordinates; Proposition 1 then rewrites its optimality equations using those coordinates Wang, pp. 9–12. This mission formalizes the reconstruction claim itself, including the exact reachable-set statement, while leaving the decision costs and Bellman equations for later work.

A completed development would provide reusable definitions for a finite-class Bayes update and its sequential iterate, an exponential-tilting model over densities, and a rank-factorized reconstruction map. The result is proved in the cited paper. The Lean declarations in this proposal compile as statements, but their sorry bodies are open proof obligations; this mission seeks machine-checked proofs of them.

Difficulty

The main verification burden is connecting several expressions for the same posterior without quietly replacing one definition by another. A sequential posterior is an iterated one-step update; the paper's closed form is a normalized product of likelihoods; and equation (3) is an exponential expression in a projected sum of transformed observations. Each form has its own denominator, so the positivity conditions that make division meaningful must be stated. The rank claim also requires full-rank assumptions on both factors of B\mathbb BB. For a prior with zero mass on a class, the posterior formula remains meaningful, but the injectivity assertion needs a separate full-support condition.

Formalization scope

Lean represents classes by Fin (N + 1), with class zero as the baseline. Row i : Fin N of the diagnostic matrix is class i.succ. Observations form a List in arrival order; ttt is its length, and the posterior is the list fold of B\mathcal BB. The observation space is a measurable type with a reference measure μ\muμ, covering continuous and discrete densities with one notation. This generalizes the paper's Y⊆R\mathcal Y\subseteq\mathbb RY⊆R; no order or arithmetic operation on observations is used by Theorem 1. Integrals are real-valued, and N=0N=0N=0 and p=0p=0p=0 are permitted. The reconstruction identity is stated for t≥1t\ge1t≥1, matching the paper and avoiding the division by zero in xtx_txt​.

The standing assumptions say that every fjf_jfj​ is measurable, f0f_0f0​ is positive at every admitted observation, its integral is one, and each tilted baseline density is integrable. These make the normalizers and the Bayesian denominators meaningful. Positivity, normalization, and integrability are explicit refinements of the paper's use of densities and normalizing constants. The prior lies in the standard simplex; θ0\theta_0θ0​ need not be positive for reconstruction. The additional full-support condition appears only in the injectivity clause. The paper's conditional independence assumption is represented by multiplying class likelihoods along the list.

The generic Bayesian update and its reachable set are separate from the paper-specific model. In particular, the posterior is an actual fold of the one-step update, and the class densities are given data constrained by exponential tilting. Defining the posterior directly as equation (3), or defining the densities to make the tilting identity true by construction, would erase the result being formalized. A complete proof needs finite-sum identities, exponential and logarithm facts, integrability facts, matrix rank and row identities, and the algebra of probability normalization. The generic Bayes definitions and list-level likelihood identity are reusable beyond this paper; contributions that establish those interfaces or close any source-derived milestone are welcome.

Selected references

  • J. Wang, Optimal Sequential Multi-class Diagnosis, Operations Research 70(5), 2022; arXiv:1506.08915v4. The mission cites the v4 preprint's pages and indices.
8 thms1 active userReviewed
ProbabilityStochastic Systems·Captain: mikedeng1

Fluid Limits for Multiclass Many-Server Queues with General Reneging Distributions and Head-of-the-Line Scheduling: Fluid-Scaled States Are Tight and Their Limit Points Are Fluid Model SolutionsResearch Paper

Motivation

Call centers, hospital wards and other service systems with many agents serve several customer classes, and customers who wait too long hang up or leave. The standard model is the multiclass many-server queue with reneging: NNN identical servers, JJJ classes, class-dependent service-time and patience-time distributions, and a control policy that decides when an idle server takes a customer and from which class. Exact analysis is out of reach once service and patience times are not exponential, so performance is studied through fluid limits: the limit of the system scaled by 1/N1/N1/N as N→∞N\to\inftyN→∞.

A. L. Puha and A. R. Ward (MOR 2022, doi:10.1287/moor.2021.1166) prove such a limit theorem for general service and patience distributions, batch arrivals and an arbitrary admissible head-of-the-line policy, including policies that idle or randomize. Earlier results covered one class without reneging (Kaspi and Ramanan, Ann. Appl. Probab. 2011), one class with reneging and unit arrivals (Kang and Ramanan, Ann. Appl. Probab. 2010), and static priority policies (Atar, Kaspi and Shimkin, Math. Oper. Res. 2014).

Setting

Fix J≥1J\ge1J≥1 classes. Class jjj has a service-time distribution GjsG^s_jGjs​ on R+\mathbb R_+R+​ with density gjsg^s_jgjs​, and a patience-time distribution GjrG^r_jGjr​ on [0,∞][0,\infty][0,∞] with density gjrg^r_jgjr​ on R+\mathbb R_+R+​, so a customer may never abandon. For a distribution GGG with density ggg, H=sup⁡{x:G(x)<1}H=\sup\{x:G(x)<1\}H=sup{x:G(x)<1} is the right edge of its support and h=g/(1−G)h=g/(1-G)h=g/(1−G) on [0,H)[0,H)[0,H) is its hazard function.

In the NNN-th queue, class-jjj customers arrive by a counting process EjNE^N_jEjN​ that may jump by more than one. The state is (αN,XN,νN,ηN)(\alpha^N,X^N,\nu^N,\eta^N)(αN,XN,νN,ηN): XjNX^N_jXjN​ counts class-jjj customers in system; νjN\nu^N_jνjN​ is the age-in-service measure, a unit atom at the service age of each class-jjj customer in service; and ηjN\eta^N_jηjN​ is the potential-queue measure, a unit atom at the time in system of each class-jjj customer whose patience has not run out. The control enters only through the entry-into-service times. A head-of-the-line (HL) policy serves each class in order of arrival. It is admissible when the entry-into-service process KNK^NKN is adapted to the filtration generated by the state history and the policy's own independent uniform random variables.

The fluid model works on the state space X=R+J×∏jM[0,Hjs)×∏jM[0,Hjr)\mathbb X=\mathbb R^J_+\times\prod_j\mathbf M[0,H^s_j)\times\prod_j\mathbf M[0,H^r_j)X=R+J​×∏j​M[0,Hjs​)×∏j​M[0,Hjr​) of finite measures with the weak topology. Given an arrival function EEE (nondecreasing, rcll, E(0)=0E(0)=0E(0)=0), a path (X,ν,η)(X,\nu,\eta)(X,ν,η) defines the auxiliary functions Bj=⟨1,νj⟩B_j=\langle1,\nu_j\rangleBj​=⟨1,νj​⟩, Qj=Xj−BjQ_j=X_j-B_jQj​=Xj​−Bj​, Dj(t)=∫0t⟨hjs,νj(u)⟩duD_j(t)=\int_0^t\langle h^s_j,\nu_j(u)\rangle duDj​(t)=∫0t​⟨hjs​,νj​(u)⟩du, Kj=Bj+Dj−Bj(0)K_j=B_j+D_j-B_j(0)Kj​=Bj​+Dj​−Bj​(0), I=1−∑jBjI=1-\sum_jB_jI=1−∑j​Bj​, and the reneging function

Rj(t)=∫0t∫0Qj(u)hjr((Fj,u)−1(y)) dy du,Fj,u(x)=⟨1[0,x],ηj(u)⟩.R_j(t)=\int_0^t\int_0^{Q_j(u)}h^r_j\big((F_{j,u})^{-1}(y)\big)\,dy\,du,\qquad F_{j,u}(x)=\langle1_{[0,x]},\eta_j(u)\rangle .Rj​(t)=∫0t​∫0Qj​(u)​hjr​((Fj,u​)−1(y))dydu,Fj,u​(x)=⟨1[0,x]​,ηj​(u)⟩.

A fluid model solution for EEE satisfies the constraints (35)–(37), the monotonicity of KKK, the balance law Xj=Xj(0)+Ej−Rj−DjX_j=X_j(0)+E_j-R_j-D_jXj​=Xj​(0)+Ej​−Rj​−Dj​, and the transport equations (46)–(47) for νj\nu_jνj​ and ηj\eta_jηj​ driven by dKjdK_jdKj​ and dEjdE_jdEj​.

Formalization targets

Goal: Theorem 1 (p. 17)

Under Assumptions 4 and 5, for every sequence of admissible HL policies:

{(XˉN,νˉN,ηˉN)}N is tight in D(X),\{(\bar X^N,\bar\nu^N,\bar\eta^N)\}_{N}\ \text{is tight in }\mathbf D(\mathbb X),{(XˉN,νˉN,ηˉ​N)}N​ is tight in D(X),

and every distributional limit point (X,ν,η)(X,\nu,\eta)(X,ν,η) is, almost surely, a continuous fluid model solution for EEE started from (X0,ν0,η0)(X^0,\nu^0,\eta^0)(X0,ν0,η0), with continuous auxiliary functions and each ηj(t)\eta_j(t)ηj​(t) charging no points. Assumptions 1–3 ask that the scaled arrivals and initial states converge. Assumption 4 asks that the hazards be bounded or lower semicontinuous near the right edges Hjs,HjrH^s_j,H^r_jHjs​,Hjr​. Assumption 5 adds continuity of EEE, an atomless η0\eta^0η0, and control of KˉN\bar K^NKˉN by C-tightness or by inequality (55).

Milestones

  1. Lemma 1, (24) and (26): the prelimit dynamic equations for ⟨f,νjN⟩\langle f,\nu^N_j\rangle⟨f,νjN​⟩ and ⟨ς,ηjN⟩\langle\varsigma,\eta^N_j\rangle⟨ς,ηjN​⟩.
  2. Lemma 6, (61): ⟨1[0,x]hjr,ηjN(t)⟩=∫0Fj,tN(x)hjr((Fj,tN)−1(y)) dy\langle1_{[0,x]}h^r_j,\eta^N_j(t)\rangle=\int_0^{F^N_{j,t}(x)}h^r_j((F^N_{j,t})^{-1}(y))\,dy⟨1[0,x]​hjr​,ηjN​(t)⟩=∫0Fj,tN​(x)​hjr​((Fj,tN​)−1(y))dy.
  3. Lemma 7: the compensator bounds ∣As,jN(φ,t)∣≤∥φ∥∞(XjN(0)+EjN(t))∫0mhjs|\mathcal A^N_{s,j}(\varphi,t)|\le\|\varphi\|_\infty(X^N_j(0)+E^N_j(t))\int_0^mh^s_j∣As,jN​(φ,t)∣≤∥φ∥∞​(XjN​(0)+EjN​(t))∫0m​hjs​ and the reneging analogue.
  4. Lemma 8: relative compactness in D(R+)\mathbf D(\mathbb R_+)D(R+​) of the fluid-scaled real processes, and their C-tightness under Assumption 5.
  5. Lemma 2 (i) and (iii): continuity and the absence of atoms propagate along fluid model solutions.

Significance

Theorem 1 characterizes fluid limits for every admissible HL policy at once. To get the fluid limit of a specific policy, one adds the equations that policy imposes and proves uniqueness for the resulting system; the paper does this for weighted random buffer selection (its Theorems 3–5). The reneging equation (40) is the main new feature: abandonment depends on the whole distribution of waiting times in queue, not only on the queue length. Allowing batch arrivals means the reneging count can only be sandwiched between two compensators, which is why Lemma 6 is needed.

The results are proved in the paper, with parts of the tightness argument delegated to Kang–Ramanan and Kaspi–Ramanan. None of them is formalized. The mission defines from scratch a measure-valued queue with abandonment, the Skorokhod space D(X)\mathbf D(\mathbb X)D(X) with J1J_1J1​ convergence, and the fluid model of §3. These definitions are reusable for the single-class models of Kaspi–Ramanan and Kang–Ramanan, for the WRBS and static-priority policies, and for other measure-valued fluid limits.

Difficulty

Each milestone is a known argument, but the obvious route to the goal breaks at the reneging term. A naive limit passage in the prelimit reneging count RjNR^N_jRjN​ fails: with batch arrivals, RjNR^N_jRjN​ is not a functional of ηjN\eta^N_jηjN​ and QjNQ^N_jQjN​ alone, because customers with equal waiting times may be split between queue and service. Identifying the limit as the nonlinear function (40) requires bounds from above and below by marked point processes. It also requires showing that the compensators of these processes converge to a common limit, which uses that η\etaη charges no points along the limit. The tightness part needs Kurtz–Jakubowski criteria for measure-valued processes on spaces [0,H)[0,H)[0,H) whose right end may carry escaping mass. Lemma 6 is a pathwise identity for finite counting measures, and Lemma 7 is a pathwise bound.

Formalization scope

Classes are Fin J with J≥1J\ge1J≥1. Time is ℝ≥0. Customers are indexed by integers as in the paper. Event times are extended reals: a patience of ∞\infty∞ keeps a customer in the potential queue forever, and an entry time of ∞\infty∞ means the customer is never served. M[0,H)\mathbf M[0,H)M[0,H) is FiniteMeasure on the subtype [0,H)[0,H)[0,H) of ℝ≥0. J1J_1J1​ convergence uses one time change for all coordinates, with the maximum of the sup-norm and Lévy–Prokhorov distances on X\mathbb XX. Convergence in distribution, distributional limit points, relative compactness and C-tightness are stated in coupling form; on the Polish space D(X)\mathbf D(\mathbb X)D(X) these are equivalent to the usual notions. Hazards take values in [0,∞][0,\infty][0,∞] and are 000 off [0,H)[0,H)[0,H). The integrals in (37), (40), (41), Lemma 6 and Lemma 7 are integrals of nonnegative functions. The derivative indicators that define the departure, reneging and marked processes are written as explicit counts.

Four readings are disclosed in the items. The goal's initial condition is stated through a coupling carrying (E,X0,ν0,η0)(E,X^0,\nu^0,\eta^0)(E,X0,ν0,η0) jointly, because the limit point lives on its own probability space. Gs,GrG^s,G^rGs,Gr in (46)–(47) are read as Gjs,GjrG^s_j,G^r_jGjs​,Gjr​. In Lemma 8, the (K.2) branch of Assumption 2 is read as relative compactness of KˉjN\bar K^N_jKˉjN​, because the printed (K.2) does not imply it. The goal encodes Assumption 2 as printed. The head-of-the-line order is the HL condition stated in words on pp. 8–9 rather than inequality (10) alone, which would let a policy skip a waiting customer for good and serve younger ones.

The statement admits trivializing formalizations, and they are excluded here. The HL predicate contains only the defining constraints of Definitions 1–3: entry no earlier than arrival and before abandonment (9), the head-of-the-line order (the HL condition of pp. 8–9, which contains (10): a customer can enter service only after every earlier customer of its class has entered or abandoned), the state constraints (S.1), (S.2), (15) and (16), and rcll paths. The identities (18)–(26) are not in it. All standing independence and distributional hypotheses and the adaptedness of KNK^NKN are explicit. The auxiliary functions R,D,B,Q,K,IR,D,B,Q,K,IR,D,B,Q,K,I are computed from the path, never chosen. The policy is not specialized: idling and randomization are allowed.

Not posed: Lemmas 3–5, 9 and 10, Theorem 2, and Lemma 2 (ii) and (iv). Lemma 4 needs a server-assignment filtration, Lemma 5 needs predictable quadratic variations, and Theorem 2, Lemma 9 and Lemma 10 need measure-valued marked point processes on [0,H)×R+[0,H)\times\mathbb R_+[0,H)×R+​. Formalizing these, or a sorry-free construction of an admissible system satisfying Assumptions 1–5, is welcome.

Selected references

  • A. L. Puha, A. R. Ward, Fluid limits for multiclass many-server queues with general reneging distributions and head-of-the-line scheduling, Mathematics of Operations Research, Articles in Advance, 2021 (2022). https://doi.org/10.1287/moor.2021.1166
  • H. Kaspi, K. Ramanan, Law of large numbers limits for many-server queues, Annals of Applied Probability 21 (2011) 33–114. https://doi.org/10.1214/10-AAP681
  • W. Kang, K. Ramanan, Fluid limits of many-server queues with reneging, Annals of Applied Probability 20 (2010) 2204–2260. https://doi.org/10.1214/10-AAP683
  • R. Atar, H. Kaspi, N. Shimkin, Fluid limits for many-server systems with reneging under a priority policy, Mathematics of Operations Research 39 (2014) 672–696. https://doi.org/10.1287/moor.2013.0632
  • S. N. Ethier, T. G. Kurtz, Markov Processes: Characterization and Convergence, Wiley, 1986. https://doi.org/10.1002/9780470316658
10 thms1 active userReviewed
ProbabilityStochastic Systems·Captain: mikedeng1

Exploiting Random Lead Times for Significant Inventory Cost Savings: With Exponential Lead Times, the Generalized Base-Stock Policy's Scaled Stationary Inventory Level Converges to N(0, 1/(βγ))Research Paper

Motivation

In continuous-review inventory systems with backlogging, the classical constant base-stock (CBS) policy orders one unit whenever one unit of demand arrives, so that the inventory position (net inventory plus in-transit inventory) stays constant. Under CBS with Poisson demand of rate rrr and i.i.d. lead times with mean 1/β1/\beta1/β, the in-transit inventory has a Poisson law with mean r/βr/\betar/β, and the steady-state net inventory level Yr(∞)Y^r(\infty)Yr(∞) fluctuates on the scale r\sqrt rr​: E∣Yr(∞)∣/r→2(2πβ)−1/2\mathbb E|Y^r(\infty)|/\sqrt r \to 2(2\pi\beta)^{-1/2}E∣Yr(∞)∣/r​→2(2πβ)−1/2. Average holding and backlog costs therefore grow like r\sqrt rr​.

Stolyar and Wang (arXiv:1801.02646, Operations Research 2022) observe that when lead times are random and orders may cross, an inventory manager can exploit the randomness: by ordering more aggressively when the inventory level is low and less when it is high, the variance of the inventory level can be reduced. Their generalized base-stock (GBS) policy implements this with a single gain parameter γ\gammaγ. The main theoretical result, Theorem 1, identifies the diffusion-scale steady-state law of the inventory level under GBS when lead times are exponential: it is normal with variance 1/(βγ)1/(\beta\gamma)1/(βγ), a factor γ\gammaγ below the CBS variance. With γ=γ(r)→∞\gamma = \gamma(r) \to \inftyγ=γ(r)→∞ (Corollary 2), the scaled average cost E∣Yr(∞)∣/r\mathbb E|Y^r(\infty)|/\sqrt rE∣Yr(∞)∣/r​ tends to 000, so GBS becomes infinitely better than CBS in the limit.

Setting

Demand units (customers) arrive as a Poisson process of rate rrr. Supply units (items) are ordered by the policy and each arrives after an independent exponential lead time with mean 1/β1/\beta1/β. Let YYY be the net inventory level (positive or negative, backlog allowed) and Z≥0Z \ge 0Z≥0 the number of items in transit. Write a∨b=max⁡(a,b)a\vee b=\max(a,b)a∨b=max(a,b), a∧b=min⁡(a,b)a\wedge b=\min(a,b)a∧b=min(a,b), a+=a∨0a^+=a\vee 0a+=a∨0, a−=(−a)+a^-=(-a)^+a−=(−a)+.

The GBS policy has an integer gain γ≥1\gamma\ge 1γ≥1, a truncation level f>0f>0f>0 and a centering x∗x_*x∗​. With X∗=r/βX_* = r/\betaX∗​=r/β it uses the target and the truncated target

X=X∗−γ(Y−x∗),T=[X∧(X∗+f)]∨0.X = X_* - \gamma(Y - x_*), \qquad T = [X\wedge(X_*+f)]\vee 0 .X=X∗​−γ(Y−x∗​),T=[X∧(X∗​+f)]∨0.

At every arrival: a customer sets Y:=Y−1Y := Y-1Y:=Y−1; an item sets Y:=Y+1Y := Y+1Y:=Y+1, Z:=Z−1Z := Z-1Z:=Z−1. Then TTT is recomputed and A=⌈T−Z⌉∨0A = \lceil T - Z\rceil\vee 0A=⌈T−Z⌉∨0 items are ordered, Z:=Z+AZ := Z + AZ:=Z+A. Under exponential lead times (Y,Z)(Y,Z)(Y,Z) is a continuous-time Markov chain on Z×Z+\mathbb Z\times\mathbb Z_+Z×Z+​ with customer rate rrr and item rate βZ\beta ZβZ.

The systems are indexed by r→∞r\to\inftyr→∞, with f=f(r)f = f(r)f=f(r) satisfying

f(r)/r→0,f(r)/r→∞,f(r)/r\to 0, \qquad f(r)/\sqrt r\to\infty,f(r)/r→0,f(r)/r​→∞,

and x∗=x∗(r)x_* = x_*(r)x∗​=x∗​(r) arbitrary. A stationary distribution ppp of a countable chain with rates q(s,s′)q(s,s')q(s,s′) and exit rates q(s)q(s)q(s) is a probability mass function with finite mean exit rate satisfying global balance ∑sp(s)q(s,s′)=p(s′)q(s′)\sum_s p(s)q(s,s') = p(s')q(s')∑s​p(s)q(s,s′)=p(s′)q(s′).

Formalization targets

Goal: Theorem 1

For every r>0r>0r>0 some state is reachable from every state, for all large rrr the chain has a unique stationary distribution, and for every choice of stationary distributions

Yr(∞)−x∗r⇒N(0,(βγ)−1),E[Yr(∞)−x∗r]±→(2πβγ)−1/2.\frac{Y^r(\infty)-x_*}{\sqrt r}\Rightarrow\mathcal N\big(0,(\beta\gamma)^{-1}\big),\qquad \mathbb E\Big[\frac{Y^r(\infty)-x_*}{\sqrt r}\Big]^{\pm}\to(2\pi\beta\gamma)^{-1/2}.r​Yr(∞)−x∗​​⇒N(0,(βγ)−1),E[r​Yr(∞)−x∗​​]±→(2πβγ)−1/2.

Milestones

In proof order: the stability part of Theorem 1; the conservation law βEZr(∞)=r\beta\mathbb E Z^r(\infty)=rβEZr(∞)=r (14); the mean bound for the reflected random walk QnQ_nQn​ that dominates the order gap (p. 16, EC.3); Lemma 4, EDr(∞)≤ED~r(∞)≤C\mathbb E D^r(\infty)\le\mathbb E\tilde D^r(\infty)\le CEDr(∞)≤ED~r(∞)≤C; the conservation law (19) and Lemma 5 for the artificial process, a birth–death chain Y^r\hat Y^rY^r that keeps Z^=⌈T^⌉\hat Z = \lceil\hat T\rceilZ^=⌈T^⌉ exactly; Corollary 6; the stochastic comparison Lemma 8, Y^r(∞)≤stYr(∞)\hat Y^r(\infty)\le_{st}Y^r(\infty)Y^r(∞)≤st​Yr(∞) and T^r(∞)≥stTr(∞)\hat T^r(\infty)\ge_{st}T^r(\infty)T^r(∞)≥st​Tr(∞); the target-gap bound (31); the transfer Lemma 9; and the reduction (36), EYr(∞)/r−EY^r(∞)/r→0\mathbb E Y^r(\infty)/\sqrt r - \mathbb E\hat Y^r(\infty)/\sqrt r\to 0EYr(∞)/r​−EY^r(∞)/r​→0. For §4 milestones the centering is x∗≡0x_*\equiv 0x∗​≡0, as in the paper (p. 12).

Significance

The result. Theorem 1 quantifies the value of random lead times: the GBS policy divides the asymptotic variance of the steady-state inventory level by γ\gammaγ, and the half-normal limits of E[⋅]±\mathbb E[\cdot]^\pmE[⋅]± give the asymptotic holding and backlog costs, which is what makes the cost comparison with CBS (bounds (12)–(13) and Corollary 2) possible. The artificial-process comparison and Lemma 9 form a reusable template for transferring a diffusion limit from a tractable one-dimensional chain to a higher-dimensional controlled chain.

Formalizing it. The result is proved in the paper; no machine-checked version exists. The mission formalizes the GBS chain itself as explicit transition rates, so the theorem cannot be discharged by a model that only resembles the policy. The stability part is only sketched in the paper ("We omit further details of the stability proof", p. 13), so a complete formal proof also fills that gap. Lemma 5 is a stationary Halfin–Whitt-type limit for a birth–death chain with state-dependent up-rates, including uniform integrability of the scaled tails.

Difficulty

The two-dimensional chain has no product form, and its in-transit inventory does not track the target: after an item arrival YYY rises by one while the target XXX falls by γ\gammaγ, but only one item leaves transit, so ZZZ overshoots its target and the gap D=Z−⌈T⌉D = Z - \lceil T\rceilD=Z−⌈T⌉ is random. The natural argument, "the chain behaves like the birth–death chain with Z=⌈T⌉Z = \lceil T\rceilZ=⌈T⌉", fails on two counts: the steady-state gap must be bounded uniformly in rrr (Lemma 4), although the rates of the chain grow with rrr and ZZZ is unbounded on the state space, and the truncations of TTT at 000 and at X∗+fX_*+fX∗​+f produce error terms that must vanish after 1/r1/\sqrt r1/r​ scaling (Corollary 6). Stochastic order alone gives only an inequality between the two laws; closing it needs the mean identity (36), which rests on the exact conservation laws (14) and (19).

Formalization scope

All objects live in the namespace RandomLeadGBS.Diffusion. The state is (Y,Z)∈Z×N(Y,Z)\in\mathbb Z\times\mathbb N(Y,Z)∈Z×N; the demand rate rrr is real and r→∞r\to\inftyr→∞ along atTop. The parameters are β>0\beta>0β>0, γ∈N\gamma\in\mathbb Nγ∈N with γ≥1\gamma\ge 1γ≥1, f:R→Rf:\mathbb R\to\mathbb Rf:R→R with f>0f>0f>0 and (10), and arbitrary x∗:R→Rx_*:\mathbb R\to\mathbb Rx∗​:R→R. The policy is encoded literally: target (4), truncated target (5), order step Z+toNat⌈T−Z⌉Z + \mathrm{toNat}\lceil T - Z\rceilZ+toNat⌈T−Z⌉ applied after every arrival; the item rate is βZ\beta ZβZ and the cap is X∗+fX_*+fX∗​+f (not X∗+γx∗+fX_*+\gamma x_*+fX∗​+γx∗​+f).

Conventions committed to:

  • a stationary distribution is a probability solution of global balance with finite mean exit rate (the state space is unbounded in ZZZ);
  • expectations and probabilities are series over states, and every expectation bound also asserts the summability it needs;
  • weak convergence ⇒\Rightarrow⇒ is stated by its definition, through every bounded continuous test function, against Mathlib's gaussianReal 0 (1/(βγ));
  • the paper's "irreducible" is replaced by "some state is reachable from every state", because states with Z<⌈T(Y)⌉Z<\lceil T(Y)\rceilZ<⌈T(Y)⌉ are never entered;
  • "uniformly in rrr" and "sup⁡r\sup_rsupr​" (Lemma 4, Corollary 6, (31)) are stated for all r≥r0r\ge r_0r≥r0​.

The goal concerns the two-dimensional GBS chain; stating it for the artificial chain, or only for an empty family of stationary laws, would trivialize it. The uniqueness clause and the existence clause of Lemma 5 rule out the vacuous reading.

A complete development needs a theory of stationary laws of countable continuous-time chains (existence via drift/Foster criteria, embedded-chain time change), birth–death stationary laws and their local normal asymptotics, and stochastic comparison of chains by coupling. The countable-chain stationarity layer and Lemma 9 are reusable beyond this mission. Proofs of any milestone, and alternative routes to the goal, are welcome. Propositions 3 and 7 (generic comparison results used for Lemma 4 and Lemma 8) are not posed as milestones.

Selected references

  • A. L. Stolyar and Q. Wang, Exploiting random lead times for significant inventory cost savings, Operations Research, 2022; preprint arXiv:1801.02646v6, 2021. https://arxiv.org/abs/1801.02646
  • S. M. Disney, A. Maltz, X. Wang and R. D. H. Warburton, Inventory management for stochastic lead times with order crossovers, European Journal of Operational Research 248(2), 473–486, 2016. https://doi.org/10.1016/j.ejor.2015.07.047
  • S. Halfin and W. Whitt, Heavy-traffic limits for queues with many exponential servers, Operations Research 29(3), 567–588, 1981. https://doi.org/10.1287/opre.29.3.567
  • T. Lindvall, Lectures on the Coupling Method, Dover, 1992.
  • S. P. Meyn and R. L. Tweedie, Markov Chains and Stochastic Stability, Springer, 1993. https://doi.org/10.1007/978-1-4471-3267-7
14 thms1 active userReviewed
Algorithmic Game TheoryEconomics·Captain: mikedeng1

Implications of Price Matching in Supply Chain Negotiation: Under Symmetric Bargaining Power the Seller Prefers Price Matching to Simultaneous Negotiation and Both Buyers Prefer the ReverseResearch Paper

Motivation

A seller that supplies several competing buyers may negotiate a wholesale price with each of them independently and simultaneously (simultaneous negotiation, SN), or may commit to charge all of them the same price (price-matching negotiation, PM). PM is the bargaining counterpart of the most-favoured-customer clause, a contractual guarantee that a customer receives the best price the vendor gives anyone. Retailers such as Walmart and Costco demand such clauses from common suppliers, and buyers command them expecting to obtain the same price as their rivals (Shang and Cai, §1). Shang and Cai (MSOM, 2021; SSRN 3021563) ask who benefits from PM relative to SN when wholesale prices are the outcome of Nash bargaining between the seller and each buyer, and show that with symmetric buyers PM benefits the seller and hurts both buyers.

The model builds on the bilateral-bargaining framework of Horn and Wolinsky (1988), in which a common supplier negotiates simultaneously with competing downstream firms, and on the analysis of negotiated wholesale contracts by O'Brien and Shaffer (2005, RAND Journal of Economics 36, 573–595) and Dukes, Gal-Or and Srinivasan (2006).

Setting

A seller supplies an intermediate good to Buyer 1 and Buyer 2, who each turn one unit of it into one unit of a substitutable end product. With demand base a>0a > 0a>0 and substitutability γ∈[0,1)\gamma \in [0,1)γ∈[0,1), retail prices and quantities satisfy the reverse demand system pi=a−xi−γxjp_i = a - x_i - \gamma x_jpi​=a−xi​−γxj​, j=3−ij = 3 - ij=3−i.

The game has two stages.

  1. Retail stage. Given wholesale prices (wi,wj)(w_i, w_j)(wi​,wj​), the buyers set retail prices in a differentiated Bertrand game. Buyer iii earns πi=(pi−wi)xi\pi_i = (p_i - w_i)x_iπi​=(pi​−wi​)xi​ and the seller earns πs=wixi+wjxj\pi_s = w_i x_i + w_j x_jπs​=wi​xi​+wj​xj​. Writing p^i(wi,wj)\hat p_i(w_i,w_j)p^​i​(wi​,wj​) and x^i(wi,wj)\hat x_i(w_i,w_j)x^i​(wi​,wj​) for the equilibrium prices and quantities, the second-stage profits are πi(wi,wj)\pi_i(w_i,w_j)πi​(wi​,wj​) and πs(wi,wj)\pi_s(w_i,w_j)πs​(wi​,wj​).
  2. Negotiation stage. Wholesale prices are set by generalized Nash bargaining: Buyer iii has bargaining power θi∈(0,1)\theta_i \in (0,1)θi​∈(0,1) against the seller, who has 1−θi1-\theta_i1−θi​. A bargain over a price w≥0w \ge 0w≥0 with buyer payoff U(w)U(w)U(w), seller payoff V(w)V(w)V(w) and disagreement point (0,d)(0, d)(0,d) selects the individually rational www maximizing U(w)θi(V(w)−d)1−θiU(w)^{\theta_i}(V(w) - d)^{1-\theta_i}U(w)θi​(V(w)−d)1−θi​.

If a negotiation fails, the seller can sell to the other buyer alone. Then Buyer jjj faces demand a−pja - p_ja−pj​, and the seller earns the monopoly profit πsM(wj)\pi^M_s(w_j)πsM​(wj​).

The two mechanisms are:

  • Simultaneous negotiation (SN). The seller bargains with both buyers at once. A pair (w1SN,w2SN)(w^{SN}_1, w^{SN}_2)(w1SN​,w2SN​) is an equilibrium if each wiSNw^{SN}_iwiSN​ solves dyad iii's bargain with payoffs πi(⋅,wjSN)\pi_i(\cdot, w^{SN}_j)πi​(⋅,wjSN​), πs(⋅,wjSN)\pi_s(\cdot, w^{SN}_j)πs​(⋅,wjSN​) and disagreement point (0,πsM(wjSN))(0, \pi^M_s(w^{SN}_j))(0,πsM​(wjSN​)).
  • Price matching (PM). The seller bargains with one buyer first; the other buyer accepts the same price. If the first bargain fails, the seller bargains with the other buyer alone, reaching w^j\hat w_jw^j​. The PM price wPMw^{PM}wPM solves the bargain with payoffs πi(w,w)\pi_i(w,w)πi​(w,w), πs(w,w)\pi_s(w,w)πs​(w,w) and disagreement point (0,πsM(w^j))(0, \pi^M_s(\hat w_j))(0,πsM​(w^j​)).

The Lean development uses this notation: pHat, xHat, piB (πi\pi_iπi​), piS (πs\pi_sπs​), piMb, piMs, IsNashBargain, IsSNEq, IsPMPrice.

Formalization targets

Goal: Theorem 1 (p. 10)

With symmetric power θ1=θ2=θ∈(0,1)\theta_1 = \theta_2 = \theta \in (0,1)θ1​=θ2​=θ∈(0,1), an SN equilibrium and a PM equilibrium price exist, and for every SN equilibrium (w1,w2)(w_1, w_2)(w1​,w2​) and every PM price www

πs(w1,w2)<πs(w,w),π1(w,w)<π1(w1,w2),π2(w,w)<π2(w2,w1).\pi_s(w_1,w_2) < \pi_s(w,w), \qquad \pi_1(w,w) < \pi_1(w_1,w_2), \qquad \pi_2(w,w) < \pi_2(w_2,w_1).πs​(w1​,w2​)<πs​(w,w),π1​(w,w)<π1​(w1​,w2​),π2​(w,w)<π2​(w2​,w1​).

The goal is stated about the equilibrium predicates, not about the closed-form prices, so that it asserts the comparison of the two mechanisms as the paper defines them.

Milestones

  1. The retail equilibrium p^i\hat p_ip^​i​, x^i\hat x_ix^i​ and the profit formulas (1)–(2), §3.1.1.
  2. The single-channel monopoly price and profits, §3.1.2.
  3. The single-channel bargain w^j=(1−θj)a/2\hat w_j = (1-\theta_j)a/2w^j​=(1−θj​)a/2, with πsM(w^j)=(1−θj2)a2/8\pi^M_s(\hat w_j) = (1-\theta_j^2)a^2/8πsM​(w^j​)=(1−θj2​)a2/8, §3.2.2.
  4. The SN equilibrium is unique and equals wSN=(1−θ)(1−γ)(2+γ)a2(1−γ)(2+γ)+γ2(1+γ)(2−γ)θw^{SN} = \dfrac{(1-\theta)(1-\gamma)(2+\gamma)a}{2(1-\gamma)(2+\gamma) + \gamma^2(1+\gamma)(2-\gamma)\theta}wSN=2(1−γ)(2+γ)+γ2(1+γ)(2−γ)θ(1−θ)(1−γ)(2+γ)a​, §3.2.1.
  5. The PM equilibrium price is unique and equals wiPM=a4(3−θi−1−γθi+γ2θi+θi2+(1+γ)(2−γ)θiθj2)w^{PM}_i = \tfrac a4\bigl(3 - \theta_i - \sqrt{1 - \gamma\theta_i + \gamma^2\theta_i + \theta_i^2 + (1+\gamma)(2-\gamma)\theta_i\theta_j^2}\bigr)wiPM​=4a​(3−θi​−1−γθi​+γ2θi​+θi2​+(1+γ)(2−γ)θi​θj2​​), §3.2.2, for general θi,θj\theta_i, \theta_jθi​,θj​.
  6. 0<wSN<wPM≤a/20 < w^{SN} < w^{PM} \le a/20<wSN<wPM≤a/2 under symmetric power, §3.2.2.
  7. At a common price, πs(w,w)\pi_s(w,w)πs​(w,w) increases on [0,a/2][0, a/2][0,a/2] and πi(w,w)\pi_i(w,w)πi​(w,w) decreases on [0,a][0,a][0,a] (proof of Theorem 1, Online Supplement p. 1).

Significance

Theorem 1 separates the effect of the negotiation mechanism from any asymmetry between buyers: under symmetric power the buyers are interchangeable, SN gives them the same price, and still the mechanism alone moves surplus from the buyers to the seller. The paper's explanation is that price matching lowers a buyer's incentive to negotiate the price down, so the matched wholesale price is higher (p. 10). The paper builds its asymmetric-power, side-payment and spot-market results (§§4–6) on this comparison.

The result is proved in the paper; nothing here is open. Neither the result nor its ingredients have been formalized. The mission produces a machine-checked account of a two-stage bargaining-over-competition model: the Bertrand stage, the Nash-bargaining stage with endogenous disagreement points, and the closed-form equilibria. Two steps the paper only asserts become formal obligations: uniqueness of the SN equilibrium, and the claim that the closed-form PM price is the global maximizer of the Nash product.

Difficulty

The comparison step is short once the closed forms are known: by milestones 6 and 7 the PM price is higher, the seller's diagonal profit increases, and the buyers' decreases. The work is in the equilibrium characterizations.

  • In SN each dyad's disagreement payoff πsM(wj)\pi^M_s(w_j)πsM​(wj​) depends on the other dyad's outcome, so an equilibrium is a fixed point of two coupled Nash-product maximizations. The paper asserts uniqueness. The best-response map is not a contraction for γ\gammaγ near 111, so a contraction argument does not give uniqueness directly.
  • The Nash products involve real powers with exponents θ\thetaθ and 1−θ1-\theta1−θ. Showing that the first-order condition picks out the global maximizer over the individually rational set requires a concavity or single-peakedness argument, not just a derivative computation.
  • The PM closed form contains a square root. Comparing it with the rational function wSNw^{SN}wSN is an inequality in three parameters (γ,θ,a)(\gamma, \theta, a)(γ,θ,a).

Formalization scope

All quantities are real numbers. The standing assumptions are a>0a > 0a>0, γ∈[0,1)\gamma \in [0,1)γ∈[0,1) and every bargaining power in (0,1)(0,1)(0,1). The sign a>0a > 0a>0 is added; the paper never states it.

  • Bargaining. A Nash bargaining solution maximizes Uθ(V−d)1−θU^{\theta}(V-d)^{1-\theta}Uθ(V−d)1−θ over the individually rational prices {w≥0:U(w)≥0, V(w)≥d}\{w \ge 0 : U(w) \ge 0,\ V(w) \ge d\}{w≥0:U(w)≥0, V(w)≥d}. Lean's Real.rpow is the real power only on non-negative bases, and the individually rational set guarantees non-negative bases.
  • Feasible domain. Equilibrium wholesale prices are required to lie in [0,a][0,a][0,a], the paper's "entire feasible domain" (p. 9; Online Supplement p. 1); bargaining deviations range over all individually rational w≥0w \ge 0w≥0.
  • Second-stage profits. πi\pi_iπi​ and πs\pi_sπs​ are defined from the demand system at (p^i,p^j)(\hat p_i,\hat p_j)(p^​i​,p^​j​), not by the closed forms (1)–(2), so milestone 1 is a genuine identity.
  • PM outcome. In PM the matched price applies to both buyers, and the seller negotiates first with Buyer 1. With symmetric power the choice of first buyer is immaterial.
  • Preference. "Prefers" is strict; θ∈(0,1)\theta \in (0,1)θ∈(0,1) excludes the paper's equality cases.

Substituting the closed-form prices wSNw^{SN}wSN, wPMw^{PM}wPM into the profit comparison would turn Theorem 1 into the algebra of milestones 6–7. The goal is therefore stated for every equilibrium of each mechanism, together with existence of both equilibria, so neither the closed forms nor a vacuous equilibrium set can discharge it.

A complete development needs:

  • concavity of the Nash product along the individually rational interval (reusable for other one-dimensional Nash-bargaining models);
  • the Bertrand equilibrium of a linear differentiated duopoly;
  • real-analysis inequalities involving square roots.

Contributions to any milestone, or to reusable lemmas on maximizing uθv1−θu^{\theta} v^{1-\theta}uθv1−θ over an interval, are welcome.

Selected references

  • W. Shang and G. Cai, Implications of Price Matching in Supply Chain Negotiation, Manufacturing & Service Operations Management, 2021 (accepted manuscript, SSRN 3021563). https://ssrn.com/abstract=3021563 — DOI https://doi.org/10.1287/msom.2021.0992
  • H. Horn and A. Wolinsky, Bilateral Monopolies and Incentives for Merger, RAND Journal of Economics 19(3), 408–419, 1988. https://doi.org/10.2307/2555721
  • D. P. O'Brien and G. Shaffer, Bargaining, Bundling, and Clout: The Portfolio Effects of Horizontal Mergers, RAND Journal of Economics 36, 573–595, 2005.
  • A. Dukes, E. Gal-Or and K. Srinivasan, Channel Bargaining with Retailer Asymmetry, Journal of Marketing Research 43, 84–97, 2006. https://doi.org/10.1509/jmkr.43.1.84
  • E. Kalai, Nonsymmetric Nash Solutions and Replications of 2-Person Bargaining, International Journal of Game Theory 6(3), 129–133, 1977. https://doi.org/10.1007/BF01770215
10 thms1 active userReviewed
Algorithmic Game TheoryTheoretical Computer Science·Captain: mikedeng1

Dynamic Fair Resource Division: The Optimal Fairness Ratio of d-Uniform Control Vectors Is 1/((d+1)·ln((d+1)/d)), Attained by the Frugally Fair AlgorithmResearch Paper

Motivation

Shared computing infrastructure — a data-centre cluster, a cloud platform's memory — is divided among users who arrive over time. Two requirements pull against each other. Proportional fairness asks that, when ttt users are present, each holds at least a 1/t1/t1/t fraction of the resource. Limited disruption asks that the system not take resource away from running jobs at every arrival, because preemption and reallocation are costly; production cluster managers such as Google's Borg impose an explicit cap on the number of task disruptions an update may cause (Verma et al., 2015, as quoted in Vardi–Psomas–Friedman 2022, pp. 945–946). If reallocation were free, the resource could be re-split evenly at every arrival; if fairness did not matter, the first user could keep everything. The question is how much fairness survives a hard cap on disruptions.

Earlier work on dynamic fair division (Walsh 2011, online cake cutting; Kash, Procaccia and Shah 2014) did not allow reallocation at all. Friedman, Psomas and Vardi (EC 2017) introduced control vectors and gave non-optimal algorithms. Vardi, Psomas and Friedman (Mathematics of Operations Research, 2022) give an optimal algorithm for every control vector once the target ratio is known, and compute the optimal ratio exactly for uniform disruption budgets. This mission formalizes that exact computation.

Setting

One unit of a homogeneous, divisible resource is shared. Users arrive one at a time and never leave (the paper treats arrivals and departures separately; §§3–5, which this mission covers, are arrivals only). Users are labelled by arrival order, and XktX^t_kXkt​ is the share of the (k+1)(k+1)(k+1)-th user when ttt users are present.

A control vector ψ\psiψ of length n∈{1,2,… }∪{∞}n \in \{1,2,\dots\}\cup\{\infty\}n∈{1,2,…}∪{∞} fixes, for each arrival 1≤t≤n1 \le t \le n1≤t≤n, a budget ψ[t]\psi[t]ψ[t]: when the ttt-th user arrives, at most ψ[t]\psi[t]ψ[t] of the t−1t-1t−1 users already present may have their share strictly reduced (disrupted). Increasing a share is free. The ddd-uniform control vector has ψ[t]=d\psi[t] = dψ[t]=d for every ttt, of any length, finite or infinite.

An allocation sequence X=(X1,X2,… )X = (X^1, X^2, \dots)X=(X1,X2,…) for ψ\psiψ is feasible if at every time 1≤t≤n1 \le t \le n1≤t≤n the shares are nonnegative, sum to at most 111, and respect the disruption budget. Its fairness ratio is

inf⁡1≤t≤n  t⋅min⁡k<tXkt,\inf_{1 \le t \le n}\; t \cdot \min_{k<t} X^t_k ,1≤t≤ninf​t⋅k<tmin​Xkt​,

the smallest share relative to the proportional share 1/t1/t1/t, in the worst time step. The fairness ratio of the control vector, FAIRNESS(ψ)\mathrm{FAIRNESS}(\psi)FAIRNESS(ψ), is the supremum of this quantity over all feasible allocation sequences, and

FAIRNESS(d-uniform)=inf⁡ψ d-uniformFAIRNESS(ψ).\mathrm{FAIRNESS}(d\text{-uniform}) = \inf_{\psi\ d\text{-uniform}} \mathrm{FAIRNESS}(\psi).FAIRNESS(d-uniform)=ψ d-uniforminf​FAIRNESS(ψ).

The frugally fair algorithm FFA(ψ,σ)\mathrm{FFA}(\psi,\sigma)FFA(ψ,σ) takes a target ratio σ\sigmaσ: when user ttt arrives, it gives the newcomer σ/t\sigma/tσ/t and cuts the ψ[t]\psi[t]ψ[t] largest existing shares to σ/t\sigma/tσ/t. It may allocate more than one unit; StS^tSt denotes its total allocation at time ttt.

Formalization targets

Goal: Theorem 4

For every integer d≥1d \ge 1d≥1,

FAIRNESS(d-uniform)  =  σd∗  :=  1(d+1)ln⁡ ⁣(d+1d).\mathrm{FAIRNESS}(d\text{-uniform}) \;=\; \sigma^*_d \;:=\; \frac{1}{(d+1)\ln\!\big(\frac{d+1}{d}\big)} .FAIRNESS(d-uniform)=σd∗​:=(d+1)ln(dd+1​)1​.

The goal states this in three parts: the infinite ddd-uniform vector admits one feasible allocation sequence with ratio at least σd∗\sigma^*_dσd∗​; for every σ>σd∗\sigma > \sigma^*_dσ>σd∗​ some finite ddd-uniform vector admits no feasible allocation sequence with ratio σ\sigmaσ; and the infimum over all lengths of FAIRNESS\mathrm{FAIRNESS}FAIRNESS equals σd∗\sigma^*_dσd∗​. For d=1d = 1d=1, σ1∗=1/(2ln⁡2)≈0.721\sigma^*_1 = 1/(2\ln 2) \approx 0.721σ1∗​=1/(2ln2)≈0.721.

Milestones

  1. Feasible runs of FFA yield feasible allocation sequences with the same shares and ratio σ\sigmaσ (§3, p. 949).
  2. Theorem 1 (p. 950): if any algorithm achieves a nonnegative ratio σ\sigmaσ on ψ\psiψ, then FFA(ψ,σ)\mathrm{FFA}(\psi,\sigma)FFA(ψ,σ) is feasible at every time.
  3. Lemma 1 (p. 951): FFA(ψ,σ)\mathrm{FFA}(\psi,\sigma)FFA(ψ,σ) is FFA(ψ,σ′)\mathrm{FFA}(\psi,\sigma')FFA(ψ,σ′) scaled by σ/σ′\sigma/\sigma'σ/σ′.
  4. Proposition 2 (p. 952): FAIRNESS(ψ)≤FAIRNESS(ψ′)\mathrm{FAIRNESS}(\psi) \le \mathrm{FAIRNESS}(\psi')FAIRNESS(ψ)≤FAIRNESS(ψ′) when ψ′\psi'ψ′ is a prefix of ψ\psiψ.
  5. Lemma E.1 (p. 959): the explicit allocation of FFA on the infinite ddd-uniform vector at time t=ℓ(d+1)+it = \ell(d+1)+it=ℓ(d+1)+i.
  6. StS^tSt is nondecreasing in ttt on the ddd-uniform vector (p. 959).
  7. Display (E.1) (p. 959): lim⁡ℓ→∞Sℓ(d+1)=σ(d+1)ln⁡d+1d\lim_{\ell\to\infty} S^{\ell(d+1)} = \sigma(d+1)\ln\frac{d+1}{d}limℓ→∞​Sℓ(d+1)=σ(d+1)lndd+1​.

Significance

The theorem gives a closed-form price of limited disruption: with ddd reallocations per arrival, a system can guarantee every user σd∗/t\sigma^*_d/tσd∗​/t of the resource forever, and no algorithm can guarantee more uniformly over all horizons. Since σd∗=1−12d+O(d−2)\sigma^*_d = 1 - \frac{1}{2d} + O(d^{-2})σd∗​=1−2d1​+O(d−2), a modest budget already recovers most of proportional fairness, which is the quantitative content a system designer needs to choose ddd. Combined with Theorem 1, the value also turns FFA into a concrete optimal algorithm for the infinite ddd-uniform vector, running online with no knowledge of the horizon.

The result is proved in the paper; to the best of a search of the Prove2Me catalogue, nothing about control vectors, disruption-limited fair division or this constant is formalized there. The mission produces a machine-checked model of disruption-limited dynamic allocation, a formal optimality proof for a greedy online algorithm (Theorem 1) that is reusable for other control vectors (the paper's finite-vector algorithm and ccc-gap bounds rest on the same theorem), and an explicit harmonic-sum asymptotic.

Difficulty

The impossibility half must rule out every algorithm, including ones that over-allocate early, raise shares opportunistically and choose whom to disrupt adaptively. The first idea, comparing total allocations time by time, fails: an arbitrary algorithm's allocation is not entrywise comparable with any fixed reference allocation, and a single well-chosen disruption can change which users are largest. The attainment half needs one allocation sequence for an unbounded horizon, so it cannot be obtained from finitely many feasibility checks; and the value σd∗\sigma^*_dσd∗​ is not the outcome of any finite computation but a limit of harmonic-number differences. The paper's own computation of the one-step change in total allocation contains a slip (a missing factor ℓ\ellℓ; the sign is unaffected), so that step cannot be transcribed as printed.

Formalization scope

Everything is real-valued and lives in DynFairDiv.Uniform. A control vector is ψ : ℕ → ℕ (index t≥1t \ge 1t≥1) with a horizon n : ℕ∞, where ⊤ is the infinite vector; uniformCV d is the constant vector. Allocation sequences are X : ℕ → ℕ → ℝ with users labelled by arrival order, which is equivalent to the paper's sorted vectors because feasibility uses only sums, fairness only minima, and a disruption is a strict decrease of a given user's share; the newcomer is never counted as disrupted. Fairness at least σ\sigmaσ is σ / t ≤ X t k for all 1≤t≤n1 \le t \le n1≤t≤n and k<tk < tk<t. FAIRNESS(ψ)\mathrm{FAIRNESS}(\psi)FAIRNESS(ψ) is a real sSup, genuine for n≥1n \ge 1n≥1. FFA's state is a list of shares kept in nonincreasing order, so "the ψ[t]\psi[t]ψ[t] largest shares" is the front of the list. Departures (the paper's §6) are not modelled.

Added hypotheses, each disclosed in its item: d≥1d \ge 1d≥1 throughout §5.1 (the paper's own assumption); σ≥0\sigma \ge 0σ≥0 for milestones 1, 2 and 6 (the paper's fairness ratios lie in [0,1][0,1][0,1]); σ′≠0\sigma' \ne 0σ′=0 in Lemma 1; a prefix of length at least 111 in Proposition 2. Milestone 1 identifies the resulting sequence's shares with FFA's shares at each time.

A trivializing reading is ruled out: the impossibility half quantifies over every feasible allocation sequence, not only FFA-like or sorted ones, and the attainment half asks for a single sequence on the infinite horizon rather than one per finite length.

A complete development needs basic facts about FFA's list (length, order, scaling) and the harmonic asymptotics already in Mathlib (Real.tendsto_harmonic_sub_log). The FFA model and Theorem 1 apply to arbitrary control vectors and are reusable for the paper's finite-vector and ccc-gap results; contributions of auxiliary lemmas about FFA (length, sortedness, scaling) are welcome.

Selected references

  • S. Vardi, A. Psomas, E. Friedman, Dynamic Fair Resource Division, Mathematics of Operations Research 47(2):945–968, 2022. https://doi.org/10.1287/moor.2021.1155
  • E. Friedman, C.-A. Psomas, S. Vardi, Controlled Dynamic Fair Division, Proc. 18th ACM Conference on Economics and Computation (EC), 461–478, 2017.
  • I. A. Kash, A. D. Procaccia, N. Shah, No Agent Left Behind: Dynamic Fair Division of Multiple Resources, Journal of Artificial Intelligence Research 51:579–603, 2014.
  • T. Walsh, Online Cake Cutting, Algorithmic Decision Theory (ADT), Springer, 292–305, 2011.
  • A. Verma, L. Pedrosa, M. Korupolu, D. Oppenheimer, E. Tune, J. Wilkes, Large-Scale Cluster Management at Google with Borg, Proc. 10th European Conference on Computer Systems (EuroSys), 2015.
9 thms1 active userReviewed
Markov ChainStochastic Systems·Captain: mikedeng1

Smart Policies for Multisource Inventory Systems and General Tandem Queues with Order Tracking and Expediting: State Space and Product-Form Stationary Distribution under Policy-MResearch Paper

Why state-dependent sourcing matters

An inventory system can replenish demand from several stages of a supply process. When an order is placed, the manager chooses where it enters; when work finishes at one stage, the manager may send it past intermediate stages to speed delivery. These choices change the future congestion of every stage. A stationary occupancy distribution is needed to calculate how often stock is available, how often demand waits, and how often expediting is used. Song, Xiao, Zhang, and Zipkin study such policies through a tandem network whose arrival and routing decisions depend on the current occupancies. They identify Policy-M as one family with a product-form stationary distribution, despite this state dependence (Song et al., SSRN 3702068, §§3–4).

The paper's Theorem 1 is the target here. Its formula permits the system's long-run inventory cost to be evaluated by summing over states with known product weights. The same paper develops a distinct Policy-E family, but its theorem has a different routing rule and is outside this mission (Song et al., Theorems 1–2).

The tandem supply system

There are K≥1K\ge1K≥1 service nodes in series, numbered 1,…,K1,\ldots,K1,…,K. A state n=(n1,…,nK)∈NKn=(n_1,\ldots,n_K)\in\mathbb N^Kn=(n1​,…,nK​)∈NK records the number of units at each node. The segment occupancy n[a,b]=∑l=abnln[a,b]=\sum_{l=a}^{b}n_ln[a,b]=∑l=ab​nl​ counts units between nodes aaa and bbb, inclusive. Demand arrives at positive rate λ\lambdaλ. Node kkk completes service at positive rate μk\mu_kμk​ whenever nk>0n_k>0nk​>0. A completion at node KKK leaves the supply system.

Policy parameters are thresholds uiju_{ij}uij​ for 0≤i<j≤K0\le i<j\le K0≤i<j≤K, with values in N∪{∞}\mathbb N\cup\{\infty\}N∪{∞}. For fixed origin iii, thresholds are nondecreasing as jjj grows; for fixed destination jjj, they are nonincreasing as iii grows. Every terminal threshold uiKu_{iK}uiK​ equals ∞\infty∞. An arriving demand places a new unit at node 111 when no prefix n[1,k]n[1,k]n[1,k] reaches u0ku_{0k}u0k​; otherwise it places the unit at one node beyond the largest such kkk. A service completion at node m<Km<Km<K follows the analogous Policy-M rule using downstream segments n[m+1,k]n[m+1,k]n[m+1,k] and thresholds umku_{mk}umk​. An empty set of reached thresholds sends the unit to m+1m+1m+1 (Song et al., (5), (7), §4).

Assumption A imposes two inequalities on these thresholds: a triangle inequality ukj≤uki+uiju_{kj}\le u_{ki}+u_{ij}ukj​≤uki​+uij​ for k<i<jk<i<jk<i<j, and supermodularity uij+ukl≤ukj+uilu_{ij}+u_{kl}\le u_{kj}+u_{il}uij​+ukl​≤ukj​+uil​ for k<i<j<lk<i<j<lk<i<j<l. The indices may include the demand source 000. The admissible state set SM(u)S_M(u)SM​(u) consists of all occupancies whose segments obey n[k+1,l]≤ukln[k+1,l]\le u_{kl}n[k+1,l]≤ukl​ for every 0≤k<l≤K0\le k<l\le K0≤k<l≤K (Song et al., Assumption A and (11)).

Formalization targets

The first part of Theorem 1 identifies the states reachable from the empty system through positive-rate Policy-M transitions:

Reach⁡M(0)=SM(u)={n∈NK:n[k+1,l]≤ukl for 0≤k<l≤K}.\operatorname{Reach}_{M}(0)=S_M(u) =\{n\in\mathbb N^K:n[k+1,l]\le u_{kl}\text{ for }0\le k<l\le K\}.ReachM​(0)=SM​(u)={n∈NK:n[k+1,l]≤ukl​ for 0≤k<l≤K}.

For λ<μK\lambda<\mu_Kλ<μK​ and finite u0iu_{0i}u0i​ whenever 1≤i<K1\le i<K1≤i<K, define w(n)=∏k=1K(λ/μk)nkw(n)=\prod_{k=1}^{K}(\lambda/\mu_k)^{n_k}w(n)=∏k=1K​(λ/μk​)nk​ and G(u)=∑n∈SM(u)w(n)G(u)=\sum_{n\in S_M(u)}w(n)G(u)=∑n∈SM​(u)​w(n). The second part gives a normalized product law:

0<G(u)<∞,πM(n)=1SM(u)(n)G(u)∏k=1K(λμk)nk.0<G(u)<\infty,\qquad \pi_M(n)=\frac{\mathbf1_{S_M(u)}(n)}{G(u)} \prod_{k=1}^{K}\left(\frac{\lambda}{\mu_k}\right)^{n_k}.0<G(u)<∞,πM​(n)=G(u)1SM​(u)​(n)​k=1∏K​(μk​λ​)nk​.

The mission also requires πM\pi_MπM​ to have total mass one and to satisfy global balance for the rates generated by the stated routing rules. Its milestones are the closure claims (A-2) and (A-5), equality of the Appendix proof sets in (A-30), the indicator identity (A-25), and the local balance equations (A-12)–(A-15) (Song et al., Theorem 1 and Appendix A).

What the result gives

The formula describes the joint occupancy of all service nodes from KKK rates and the policy thresholds. Once its normalizer is known, probabilities of stock levels and congestion events can be summed over SM(u)S_M(u)SM​(u), supplying the inputs to the paper's long-run cost calculation (9). The state-space result determines which occupancies contribute to those sums (Song et al., (9), Theorem 1).

Theorem 1 is proved in the paper. This mission asks for machine-checked proofs of its formal statements, not a new queueing result. A complete development would provide reusable finite-segment reasoning, routing-set identities, and countable-state balance calculations. The declarations staged here are open Lean statements, so they do not yet constitute machine-checked proofs.

The obstacle in the balance equations

With fixed routing, a proposed product weight can be compared across neighboring states using only the node rates. Here the route of a newly arriving or completed unit changes with several segment occupancies. An incoming transition to a state may have skipped nodes, and the condition that selected its destination must hold in the predecessor state. Checking only the product ratio therefore leaves out the routing indicators. The source's local equations (A-12)–(A-15) and its proof-set identity (A-25) make those indicators explicit. The paper also gives an example showing that state-dependent routing does not yield product form for every policy (Song et al., Example 1 and Appendix A).

Formalization scope

Lean represents nnn as a function from Fin K to natural numbers, while the functions occ and seg retain the paper's one-based node numbers. The threshold array takes values in extended natural numbers so uiK=∞u_{iK}=\inftyuiK​=∞ is representable; a strict comparison with a threshold expresses the paper's “at most u−1u-1u−1” condition. Both routing rules specify the destination when no threshold is reached. A service transition is enabled only when its origin node is occupied. The theorem states K≥1K\ge1K≥1, which is implicit in the paper's numbered-node model, along with positive arrival and service rates.

The source defines π\piπ as the limiting law of a continuous-time Markov chain and invokes known ergodicity. This Lean development encodes its stationary claim by a summable, normalized product weight and the global balance equations for explicit Policy-M rates. It does not formalize the CTMC limit or uniqueness of that limit. The normalizer is a countable sum; the goal asserts summability and G(u)>0G(u)>0G(u)>0, so a default value for a divergent sum or division by zero cannot certify the claim. The state space is constructed from uuu and proved reachable rather than supplied as a hypothesis. Empty-node departures are guarded to avoid natural-number subtraction at zero. The special K=1K=1K=1 case remains in the goal; the milestone combining four separate local equations assumes K≥2K\ge2K≥2 because its node-1 and node-KKK formulas must refer to distinct nodes.

The staged definition supplies segments, thresholds, routing, transition rates, product weights, and the source's Appendix sets. Contributions may prove the listed milestones, build general lemmas about segment sums and reachability, or establish summability and the passage from local to global balance. These supporting tools can be reused for other countable-state queueing models.

Selected references

  • J.-S. Song, L. Xiao, H. Zhang, and P. Zipkin, Smart Policies for Multisource Inventory Systems and General Tandem Queues with Order Tracking and Expediting, SSRN 3702068, version of September 28, 2020; subsequently published in Operations Research, DOI 10.1287/opre.2021.2124. Source preprint.
7 thms1 active userReviewed
OptimizationTheoretical Computer Science·Captain: mikedeng1

Dynamic Resource Allocation in the Cloud with Near-Optimal Efficiency: The Multiplicative Weight Algorithm Nearly Matches the Offline Dynamic Optimum's Work While Approximately Satisfying Every SLAResearch Paper

Motivation

Cloud platforms divide a shared physical resource, such as a pool of processors, among many users at once. Each user is promised a service level agreement (SLA): a fraction β(i)\beta(i)β(i) of the resource that it should effectively receive whenever it has work to do. The operator wants two things that pull against each other: to keep the resource busy (total work done should be large) and to honour every SLA. The difficulty is informational. The operator does not see future demand, and in the setting of Perez-Salazar, Menache, Singh and Toriello (arXiv:1809.02688) it does not even see current demand: at each step it only learns which users have a non-empty queue.

The paper answers with a multiplicative weight rule. Multiplicative weight methods go back to online learning (Freund and Schapire, 1997; survey: Arora, Hazan and Kale, 2012) and can be read as entropic mirror descent (Beck and Teboulle, 2003). Their classical guarantees compare an online learner with the best fixed decision in hindsight. This paper compares the algorithm with the best time-varying offline allocation instead, and proves the comparison directly, without the standard regret analysis.

Setting

There are N≥2N \ge 2N≥2 users and discrete time steps t=1,2,…t = 1, 2, \dotst=1,2,…. Each user iii has an SLA β(i)≥0\beta(i) \ge 0β(i)≥0, with ∑iβ(i)≤1\sum_i \beta(i) \le 1∑i​β(i)≤1. Loads Lt(i)≥0L_t(i) \ge 0Lt​(i)≥0 arrive online. User iii keeps a queue Qt(i)Q_t(i)Qt​(i) of unfinished work. At each step:

  1. the decision maker observes the active users At={i:Qt(i)>0}A_t = \{i : Q_t(i) > 0\}At​={i:Qt​(i)>0};
  2. it chooses an allocation hth_tht​ with ht(i)≥0h_t(i) \ge 0ht​(i)≥0 and ∑iht(i)≤1\sum_i h_t(i) \le 1∑i​ht​(i)≤1;
  3. the load LtL_tLt​ arrives, user iii completes wt(i)=min⁡{ht(i), Lt(i)+Qt(i)}w_t(i) = \min\{h_t(i),\, L_t(i) + Q_t(i)\}wt​(i)=min{ht​(i),Lt​(i)+Qt​(i)}, and Qt+1(i)=max⁡{0, Lt(i)+Qt(i)−ht(i)}Q_{t+1}(i) = \max\{0,\, L_t(i) + Q_t(i) - h_t(i)\}Qt+1​(i)=max{0,Lt​(i)+Qt​(i)−ht​(i)}.

The total work of an allocation sequence hhh up to time ttt is workh(t)=∑i∑τ≤twτ(i)\mathrm{work}_h(t) = \sum_i \sum_{\tau \le t} w_\tau(i)workh​(t)=∑i​∑τ≤t​wτ​(i).

Algorithm 1 has parameters 0<ε≤1/100 < \varepsilon \le 1/100<ε≤1/10 and 0<η<1/30 < \eta < 1/30<η<1/3, and sets λ=ε2/(8N)\lambda = \varepsilon^2/(8N)λ=ε2/(8N). It keeps hth_tht​ in the truncated simplex Δε={x:∑ix(i)=1, x(i)≥ε/N}\Delta_\varepsilon = \{x : \sum_i x(i) = 1,\ x(i) \ge \varepsilon/N\}Δε​={x:∑i​x(i)=1, x(i)≥ε/N}, starting anywhere in Δε\Delta_\varepsilonΔε​. At step ttt it assigns the gain gt(i)=1+λg_t(i) = 1 + \lambdagt​(i)=1+λ to active users below their SLA (At1A^1_tAt1​: ht(i)<β(i)h_t(i) < \beta(i)ht​(i)<β(i)), gt(i)=1g_t(i) = 1gt​(i)=1 to the other active users (At2A^2_tAt2​) and gt(i)=0g_t(i) = 0gt​(i)=0 to inactive users. It then sets h^t+1(i)=ht(i)eηgt(i)\widehat h_{t+1}(i) = h_t(i)e^{\eta g_t(i)}ht+1​(i)=ht​(i)eηgt​(i) and takes ht+1=πΔε(h^t+1)h_{t+1} = \pi_{\Delta_\varepsilon}(\widehat h_{t+1})ht+1​=πΔε​​(ht+1​), the KL projection: the minimizer over Δε\Delta_\varepsilonΔε​ of ∑ix(i)ln⁡(x(i)/h^t+1(i))\sum_i x(i)\ln(x(i)/\widehat h_{t+1}(i))∑i​x(i)ln(x(i)/ht+1​(i)).

Formalization targets

Goal: Theorem 1 with the explicit constants of p. 10

Let β(i)≥2ε/N\beta(i) \ge 2\varepsilon/Nβ(i)≥2ε/N for all iii and let hhh be any run of Algorithm 1. For every offline allocation sequence h∗h^*h∗ and every time ttt,

workalg(t)≥(1−4ε) workh∗(t)−2000 Nε2ηln⁡Nε,\mathrm{work}_{\mathrm{alg}}(t) \ge (1 - 4\varepsilon)\,\mathrm{work}_{h^*}(t) - 2000\,\frac{N}{\varepsilon^2\eta}\ln\frac{N}{\varepsilon},workalg​(t)≥(1−4ε)workh∗​(t)−2000ε2ηN​lnεN​,

and for every user iii, every comparator h′h'h′ with 0≤ht′(i)≤β(i)0 \le h'_t(i) \le \beta(i)0≤ht′​(i)≤β(i) and every ttt,

∑τ≤twτ(i)≥(1−2ε)∑τ≤twτ′(i)−β(i) p~,p~=32N2ε3ηln⁡Nε.\sum_{\tau \le t} w_\tau(i) \ge (1 - 2\varepsilon)\sum_{\tau \le t} w'_\tau(i) - \beta(i)\,\tilde p, \qquad \tilde p = \frac{32N^2}{\varepsilon^3\eta}\ln\frac{N}{\varepsilon}.τ≤t∑​wτ​(i)≥(1−2ε)τ≤t∑​wτ′​(i)−β(i)p~​,p~​=ε3η32N2​lnεN​.

Intermediate targets (milestones)

  • Theorem 3: against every sequence of (1−4ε)(1-4\varepsilon)(1−4ε)-allocations, workh∗(t)−workalg(t)≤8Nε2ηln⁡Nε+2\mathrm{work}_{h^*}(t) - \mathrm{work}_{\mathrm{alg}}(t) \le \frac{8N}{\varepsilon^2\eta}\ln\frac{N}{\varepsilon} + 2workh∗​(t)−workalg​(t)≤ε2η8N​lnεN​+2 at every ttt.
  • Theorem 4: if β(j)>eη(1+λ)ε/N\beta(j) > e^{\eta(1+\lambda)}\varepsilon/Nβ(j)>eη(1+λ)ε/N for all jjj, each user's work by time t+⌈p~⌉t + \lceil\tilde p\rceilt+⌈p~​⌉ is at least what the constant share (1−2ε)β(i)(1-2\varepsilon)\beta(i)(1−2ε)β(i) completes by time ttt.
  • Lemma 1 and Lemma 2: per-step growth factors 1+c1 + c1+c, 1−εc1 - \varepsilon c1−εc and 1+c′1 + c'1+c′ of a busy user's allocation, with c=εη/(4N)c = \varepsilon\eta/(4N)c=εη/(4N) and c′=εηλ/(2N)c' = \varepsilon\eta\lambda/(2N)c′=εηλ/(2N).
  • Proposition 1(b): the closed form of the KL projection onto Δε\Delta_\varepsilonΔε​.
  • Proposition 2 (weak-duality half), Claims 3.1 and 3.2 (first conclusion): the amortization behind Theorem 3.
  • The scaling step of p. 10: scaling an allocation by κ≤1\kappa \le 1κ≤1 keeps at least a κ\kappaκ fraction of its work, and work is monotone in the allocation.
  • Existence and uniqueness of the KL projection, hence of the run from each start.

Significance

The result. Theorem 1 is a bicriteria guarantee against a dynamic benchmark. The work benchmark may change its allocation arbitrarily over time with full knowledge of all loads, and the algorithm, which sees only which queues are non-empty, loses only an additive amount that does not grow with time, at the price of a 4ε4\varepsilon4ε fraction of capacity. At the same time every user receives, up to a fixed delay, at least the work of a dedicated (1−2ε)(1-2\varepsilon)(1−2ε) fraction of its SLA. The paper also shows (arXiv:1809.02688, Theorem 2) that against the unscaled offline optimum every deterministic online algorithm loses Ω(T)\Omega(\sqrt T)Ω(T​) work on some load sequence, so some form of slack is unavoidable.

Formalizing it. The result is proved on paper; no machine-checked proof exists. A formalization would check the constant bookkeeping, where the printed text has slips (an ε\varepsilonε versus ε/4\varepsilon/4ε/4 in Theorem 1, a real block length s~\tilde ss~ used as an integer, an off-by-one in Claim 3.1). It would also produce a reusable deterministic queueing model and a verified closed form of the KL projection onto a truncated simplex.

Difficulty

The obvious first idea, the regret analysis of mirror descent, does not apply. Regret bounds compare with a fixed decision, but here the benchmark changes over time. The algorithm's gains also depend on its own past decisions through the queues, so they are not an exogenous sequence of losses. The paper notes (p. 10) that its guarantees do not follow from the standard analysis of mirror descent. A second obstacle is the feedback: a busy user is not the same as a user that uses its whole allocation. A user whose queue empties during a step completes less than its allocation, so utilization and work done can differ. Finally, the constants matter, because the goal states them explicitly. The printed derivation mixes ε\varepsilonε and ε/4\varepsilon/4ε/4, uses the real number s~\tilde ss~ as a block length, and is off by one step when it identifies the work done with the arrived load.

Formalization scope

  • Model. Everything is real-valued. Time is 0-based: Lean's step ttt is the paper's step t+1t + 1t+1, the queue at Lean time ttt is the paper's Qt+1Q_{t+1}Qt+1​, and ∑τ=1t\sum_{\tau=1}^{t}∑τ=1t​ is a sum over the first ttt steps. Queues start empty, a convention the paper uses but does not state. Loads form an infinite sequence; the horizon TTT is not a parameter and every bound holds at every ttt.
  • Runs. The KL projection is a predicate for positive inputs, the domain where the logarithm is defined; Algorithm 1 always projects a positive vector. A run is any sequence satisfying the recursion. The gains are computed from the queues the run itself produces before step ttt, so ht+1h_{t+1}ht+1​ never uses the load of step ttt or later. A milestone proves that each start in Δε\Delta_\varepsilonΔε​ has exactly one run, and a sorry-free sanity check exhibits one. Statements about runs are therefore not vacuous.
  • Benchmarks. "The optimal offline sequence" is replaced by every allocation sequence with ∑iht∗(i)≤1\sum_i h^*_t(i) \le 1∑i​ht∗​(i)≤1 (or ≤1−4ε\le 1 - 4\varepsilon≤1−4ε), which is equivalent and avoids a supremum. The SLA comparator ranges over every sequence with ht′(i)≤β(i)h'_t(i) \le \beta(i)ht′​(i)≤β(i), not only the constant one.
  • Disclosed readings and corrections. Part 1 of the goal uses the (1−4ε)(1 - 4\varepsilon)(1−4ε) form of p. 10 for the same ε\varepsilonε, not the headline's (1−ε)(1-\varepsilon)(1−ε). Theorem 3 carries a +2+2+2 for integer blocks. Claim 3.1 is stated in its corrected 0-based form. Theorem 4 assumes the SLA bound strictly. 0<η<1/30 < \eta < 1/30<η<1/3 (Algorithm 1's input) is assumed throughout. ∑iβ(i)≤1\sum_i \beta(i) \le 1∑i​β(i)≤1 replaces the =1= 1=1 used in the proof of Lemma 2.
  • Ruled out. The goal is stated for every run of the algorithm, with the active sets read from the run's own queues, against every offline sequence. A formalization that fixes a closed-form allocation, decouples the queues from the allocations, or restricts the benchmark to static allocations does not count.
  • Needed and reusable. The deterministic single-user queue with its work and monotonicity facts, KKT conditions for an entropic objective over a box-constrained simplex, and elementary exponential and logarithm inequalities. The queue model and the projection closed form are reusable beyond this mission. Contributions to any milestone are welcome.

Selected references

  • S. Perez-Salazar, I. Menache, M. Singh, A. Toriello, Dynamic Resource Allocation in the Cloud with Near-Optimal Efficiency, arXiv:1809.02688v3, 2021 (Operations Research, 2022). https://arxiv.org/abs/1809.02688
  • Y. Freund, R. E. Schapire, A decision-theoretic generalization of on-line learning and an application to boosting, Journal of Computer and System Sciences 55(1), 1997. https://doi.org/10.1006/jcss.1997.1504
  • S. Arora, E. Hazan, S. Kale, The multiplicative weights update method: a meta-algorithm and applications, Theory of Computing 8, 2012. https://doi.org/10.4086/toc.2012.v008a006
  • A. Beck, M. Teboulle, Mirror descent and nonlinear projected subgradient methods for convex optimization, Operations Research Letters 31(3), 2003. https://doi.org/10.1016/S0167-6377(02)00231-6
13 thms1 active userReviewed
ProbabilityStochastic Systems·Captain: mikedeng1

Large Fork-Join Queues with Nearly Deterministic Arrival and Service Times: Fluid Limit of the Maximum Queue Length with a Non-Zero Initial ConditionResearch Paper

Motivation

A fork-join queue splits every incoming job into tasks, sends one task to each of NNN parallel servers, and declares the job complete only when all tasks are done. The model describes assembly systems in which a manufacturer orders parts from NNN suppliers and can only start production when every part has arrived, as well as parallel computing and storage systems in which a request is answered only after all sub-requests return. Because a job waits for its slowest task, the performance of the system is governed by the maximum queue length over the NNN servers, and the question of interest is how this maximum behaves when NNN is large and the system is heavily loaded.

Schol, Vlasiou and Zwart (arXiv:1912.11661, Mathematics of Operations Research 47(2), 2022) study this question in a nearly deterministic regime: arrivals and services happen in almost every time slot, and the traffic intensity of each server tends to 111 at rate N−2N^{-2}N−2. Even in the empty-start case the queue lengths of different servers are dependent, because every server sees the same arrival process; the maximum of NNN dependent random walks has to be analysed directly.

Timeline. Early work treats two service stations and obtains exact asymptotics (Flatto and Hahn, 1984; Baccelli, 1985; Wright, 1992). For more stations the literature turned to approximations and bounds: Nelson and Tantawi (1988) and Baccelli and Makowski (1989) bound the mean response time, and Ko and Serfozo (2004) approximate the response-time distribution in M/M/s fork-join queues. Heavy-traffic process limits for a fixed number of stations were obtained by Varma (1990) and Nguyen (1993, 1994), with reflected Brownian motion limits, and by Lu and Pang (2015–2017) for many-server stations. Sigman and Whitt (2011) analysed nearly deterministic G/G/1 and G/D/N queues. The paper of this mission gives a process-level (fluid) limit for the largest queue when the number of servers NNN itself grows, in the nearly deterministic regime and with a non-zero initial condition.

Setting

Fix α,β>0\alpha,\beta>0α,β>0. For each NNN the system runs in discrete time slots j=1,2,…j=1,2,\dotsj=1,2,…. In slot jjj a task arrives at all NNN servers with probability p(N)=1−α/N−β/N2p^{(N)}=1-\alpha/N-\beta/N^2p(N)=1−α/N−β/N2 (indicator X(N)(j)X^{(N)}(j)X(N)(j)), and server iii completes a task with probability q(N)=1−α/Nq^{(N)}=1-\alpha/Nq(N)=1−α/N (indicator Yi(N)(j)Y_i^{(N)}(j)Yi(N)​(j)), all indicators being independent. The arrival process is A(N)(n)=∑j=1nX(N)(j)A^{(N)}(n)=\sum_{j=1}^{n}X^{(N)}(j)A(N)(n)=∑j=1n​X(N)(j) and the service process of server iii is Si(N)(n)=∑j=1nYi(N)(j)S_i^{(N)}(n)=\sum_{j=1}^{n}Y_i^{(N)}(j)Si(N)​(n)=∑j=1n​Yi(N)​(j).

Server iii starts with Qi(N)(0)=⌊rNUi⌋Q_i^{(N)}(0)=\lfloor r_NU_i\rfloorQi(N)​(0)=⌊rN​Ui​⌋ tasks, where U1,U2,…U_1,U_2,\dotsU1​,U2​,… are i.i.d. random variables, the initial queue lengths are non-negative, and rN≥0r_N\ge0rN​≥0 is a scaling sequence. By Lindley's recursion the maximum queue length at time nnn is

Q(α,β)(N)(n)=max⁡i≤Nmax⁡(sup⁡0≤k≤n[(A(N)(n)−A(N)(k))−(Si(N)(n)−Si(N)(k))],  Qi(N)(0)+A(N)(n)−Si(N)(n)).Q^{(N)}_{(\alpha,\beta)}(n)=\max_{i\le N}\max\Big(\sup_{0\le k\le n}\big[(A^{(N)}(n)-A^{(N)}(k))-(S_i^{(N)}(n)-S_i^{(N)}(k))\big],\;Q_i^{(N)}(0)+A^{(N)}(n)-S_i^{(N)}(n)\Big).Q(α,β)(N)​(n)=i≤Nmax​max(0≤k≤nsup​[(A(N)(n)−A(N)(k))−(Si(N)​(n)−Si(N)​(k))],Qi(N)​(0)+A(N)(n)−Si(N)​(n)).

The first term alone is the maximum queue length of the system started empty. The initial condition is assumed to satisfy Q(α,β)(N)(0)/(Nlog⁡N)→q(0)≥0Q^{(N)}_{(\alpha,\beta)}(0)/(N\log N)\to q(0)\ge0Q(α,β)(N)​(0)/(NlogN)→q(0)≥0 in probability (Assumption 2.3), and either UiU_iUi​ has a finite right endpoint (Assumption 2.4) or UiU_iUi​ is continuous with −log⁡P(Ui>vt)/−log⁡P(Ui>t)→h(v)-\log\mathbb P(U_i>vt)/-\log\mathbb P(U_i>t)\to h(v)−logP(Ui​>vt)/−logP(Ui​>t)→h(v) as t→∞t\to\inftyt→∞ for all v∈[0,1]v\in[0,1]v∈[0,1] (Assumption 2.5, the Gumbel-domain case).

The fluid limit is

q(t)=max⁡((2αt−βt)1(t<α2β2)+α2β1(t≥α2β2),  g(t,q(0))−βt),q(t)=\max\Big(\big(\sqrt{2\alpha t}-\beta t\big)\mathbb 1\big(t<\tfrac{\alpha}{2\beta^2}\big)+\tfrac{\alpha}{2\beta}\mathbb 1\big(t\ge\tfrac{\alpha}{2\beta^2}\big),\;g(t,q(0))-\beta t\Big),q(t)=max((2αt​−βt)1(t<2β2α​)+2βα​1(t≥2β2α​),g(t,q(0))−βt),

with g(t,q(0))=q(0)+2αtg(t,q(0))=q(0)+\sqrt{2\alpha t}g(t,q(0))=q(0)+2αt​ under Assumption 2.4 and g(t,q(0))=sup⁡{2αt u+q(0)v:u2+h(v)≤1, 0≤u,v≤1}g(t,q(0))=\sup\{\sqrt{2\alpha t}\,u+q(0)v : u^2+h(v)\le1,\ 0\le u,v\le1\}g(t,q(0))=sup{2αt​u+q(0)v:u2+h(v)≤1, 0≤u,v≤1} under Assumption 2.5. In Lean these are fluid, fluidEmpty, gFinite and gGumbel; the queue is maxQueue and the empty-start queue maxQueueEmpty.

Formalization targets

Goal: Theorem 2.1 (fluid limit)

For every T>0T>0T>0 and ϵ>0\epsilon>0ϵ>0,

P(sup⁡0≤t≤T∣Q(α,β)(N)(tN3log⁡N)Nlog⁡N−q(t)∣>ϵ)→N→∞0,\mathbb P\left(\sup_{0\le t\le T}\left|\frac{Q^{(N)}_{(\alpha,\beta)}(tN^3\log N)}{N\log N}-q(t)\right|>\epsilon\right)\xrightarrow{N\to\infty}0,P​0≤t≤Tsup​​NlogNQ(α,β)(N)​(tN3logN)​−q(t)​>ϵ​N→∞​0,

under Assumptions 2.1–2.3 and either Assumption 2.4 or Assumption 2.5 (one conjunct per case).

Milestones

  1. Section 4.2 tools. The pathwise oscillation bound of Lemma 4.1; the non-uniform Berry–Esseen rate of Lemma 4.2 for the standardized service process; the moment bound of Lemma 4.3; the extreme-value Lemma B.1 for sums of independent relatively stable variables and its consequence Lemma 4.4 for the Gaussian surrogate of the initial condition.
  2. Pointwise convergence (Section 4.3). Lemmas 4.6 and 4.8 (bounds for the empty-start system), Lemma 4.5 (its pointwise limit), Lemma 4.9 (the initial-condition part) and Theorem 4.1: Q(α,β)(N)(tN3log⁡N)/(Nlog⁡N)→q(t)Q^{(N)}_{(\alpha,\beta)}(tN^3\log N)/(N\log N)\to q(t)Q(α,β)(N)​(tN3logN)/(NlogN)→q(t) in probability for each t>0t>0t>0.
  3. Tightness (Section 4.4). Lemma 4.10 (maxima preserve uniform convergence), Lemmas 4.11–4.12 (moduli of continuity of the two parts) and Corollary 4.1 (tightness of the scaled process on [0,T][0,T][0,T]).

Significance

The result. Theorem 2.1 identifies the scaling under which the maximum queue length of a large, heavily loaded fork-join system follows a deterministic curve: time of order N3log⁡NN^3\log NN3logN and space of order Nlog⁡NN\log NNlogN. The limit separates the contribution of the empty-start dynamics, which saturate at the steady-state level α/(2β)\alpha/(2\beta)α/(2β) after time α/(2β2)\alpha/(2\beta^2)α/(2β2), from the contribution of the initial backlog, which is governed by the extreme-value class of the initial condition. For instance, with h(v)=vh(v)=vh(v)=v and q(0)=α/(2β)q(0)=\alpha/(2\beta)q(0)=α/(2β) the system "starts and stays in steady state" on an initial interval (p. 12). The theorem gives the transient behaviour needed to judge how long a backlog at the slowest supplier persists.

Formalizing it. The result is proved in the paper; no part of it has a machine-checked proof. The mission produces a formal statement of the model and of the fluid limit, and the intermediate results of Section 4. The proof combines a non-uniform Berry–Esseen inequality, extreme-value theory for maxima of sums, Doob's maximal inequality and a tightness criterion; formalizing it would give reusable Lean versions of several of these tools.

Difficulty

The obvious argument treats the NNN queues as independent reflected random walks and applies extreme-value theory to their maximum. It fails because all servers share the arrival process, so the queues are dependent. The paper handles this by splitting each queue into a common arrival term and an individual service term, but the individual terms are maxima of NNN centred random walks evaluated at a time of order N3log⁡NN^3\log NN3logN, and their Gaussian approximation must be accurate at the scale of the maximum: a uniform Berry–Esseen bound is too weak, and a non-uniform bound with the right rate in NNN is needed (Lemma 4.2). Upgrading pointwise convergence to uniform convergence on [0,T][0,T][0,T] requires a modulus-of-continuity estimate for a maximum of NNN processes (Lemmas 4.11–4.12), where the number of processes grows with NNN. The initial-condition term requires an extreme-value result for maxima of sums of two independent coordinates with different tails (Lemma B.1).

Formalization scope

All objects are defined in one module, LargeForkJoin.FluidLimit.Model. Times are natural numbers, and Q(tN3log⁡N)Q(tN^3\log N)Q(tN3logN) is evaluated at ⌊tN3log⁡N⌋\lfloor tN^3\log N\rfloor⌊tN3logN⌋, which is the paper's own convention Q(n)=Q(⌊n⌋)Q(n)=Q(\lfloor n\rfloor)Q(n)=Q(⌊n⌋). Servers are indexed 0,…,N−10,\dots,N-10,…,N−1. All variables live on one probability space; the Bernoulli laws are imposed for N≥N0N\ge N_0N≥N0​ (for small NNN, p(N)p^{(N)}p(N) can be negative), and within each N≥N0N\ge N_0N≥N0​ the arrival indicators, service indicators and initial-condition variables UiU_iUi​ used by that system are mutually independent. The independence of (Ui)(U_i)(Ui​) from arrivals and services is an added assumption: the paper uses it in its proofs but does not state it. Initial queue lengths are non-negative almost surely. Assumption 2.5 adds P(Ui>t)>0\mathbb P(U_i>t)>0P(Ui​>t)>0 for all ttt, without which its ratio is undefined; "continuous" means "no atoms". Convergence in probability to a constant is written with ϵ\epsilonϵ-events, and sup⁡0≤t≤T∣⋅∣>ϵ\sup_{0\le t\le T}|\cdot|>\epsilonsup0≤t≤T​∣⋅∣>ϵ as "some t∈[0,T]t\in[0,T]t∈[0,T] has ∣⋅∣>ϵ|\cdot|>\epsilon∣⋅∣>ϵ", the same event. Lemmas 4.4 and 4.9 speak of Gaussian variables outside the model; they are stated for the law of NNN independent pairs (standard normal, law of U1U_1U1​).

The goal is stated about the maximum queue length built from the Bernoulli indicators through (2.3), uniformly on [0,T][0,T][0,T]. A statement about a Gaussian or reflected-Brownian surrogate (the heuristics of Sections 2.2–2.3), or a pointwise statement (Theorem 4.1), would not be Theorem 2.1.

Not included: Propositions 2.1–2.3 (other scalings and the steady state), the variant of Theorem 2.1 under Assumption 2.6, Theorem 4.2 (finite-dimensional distributions, a generic fact), and Lemma 4.7, whose statement includes the stationary case t=∞t=\inftyt=∞ and needs a stationary-queue object. Contributions of general tools are welcome: a non-uniform Berry–Esseen inequality for Bernoulli sums, extreme-value results for maxima of i.i.d. sequences (relative stability, Gumbel-domain tails), and the tightness criterion of Billingsley's Theorem 7.3.

Selected references

  • D. Schol, M. Vlasiou, B. Zwart, Large fork-join queues with nearly deterministic arrival and service times, Mathematics of Operations Research 47(2), 2022; arXiv:1912.11661v3. https://arxiv.org/abs/1912.11661v3
  • L. Flatto, S. Hahn, Two parallel queues created by arrivals with two demands I, SIAM Journal on Applied Mathematics 44(5), 1984. https://doi.org/10.1137/0144072
  • F. Baccelli, A. M. Makowski, Queueing models for systems with synchronization constraints, Proceedings of the IEEE 77(1), 1989. https://doi.org/10.1109/5.21076
  • R. Nelson, A. N. Tantawi, Approximate analysis of fork/join synchronization in parallel queues, IEEE Transactions on Computers 37(6), 1988. https://doi.org/10.1109/12.2212
  • K. Sigman, W. Whitt, Heavy-traffic limits for nearly deterministic queues, Journal of Applied Probability 48(3), 2011. https://doi.org/10.1239/jap/1316796905
  • S.-S. Ko, R. F. Serfozo, Response times in M/M/s fork-join networks, Advances in Applied Probability 36(3), 2004. https://doi.org/10.1239/aap/1093962240
  • P. Billingsley, Convergence of Probability Measures, 2nd ed., Wiley, 1999. https://doi.org/10.1002/9780470316962
16 thms1 active userReviewed
EconomicsProbability·Captain: mikedeng1

Costly Interpretation of Asset Prices: With Costly Price Interpretation There Is a Unique Linear Financial Market Equilibrium, and Its Price Carries Endogenous NoiseResearch Paper

Motivation

A central question of financial economics is how much of the information dispersed among traders ends up in market prices. The workhorse model is the noisy rational expectations equilibrium (REE): investors with private signals trade a risky asset, the price aggregates their information, and each investor learns from the price as well as from her own signal (Grossman and Stiglitz 1980; Hellwig 1980). In that tradition the price is a free public signal: every investor reads off from it exactly what it contains.

Vives and Yang (SSRN 2763813; published with Mondria as Mondria, Vives and Yang 2021) drop that assumption. Extracting information from a price takes effort, so each investor observes the price signal only through her own interpretation error, and can buy sophistication to reduce that error. The interpretation errors share a common component, and the paper's first result is that this common error re-enters the price. Even with no exogenous noise traders, costly interpretation makes the price noisy. This mission formalizes that first result: the existence, uniqueness and closed form of the linear financial-market equilibrium (Proposition 1, with the remark au>0a_u > 0au​>0 that follows it).

Setting

A risky asset pays V~=v~+ξ~\tilde V = \tilde v + \tilde\xiV~=v~+ξ~​, with independent v~∼N(0,τv−1)\tilde v \sim N(0,\tau_v^{-1})v~∼N(0,τv−1​) (learnable) and ξ~∼N(0,τξ−1)\tilde\xi \sim N(0,\tau_\xi^{-1})ξ~​∼N(0,τξ−1​) (unlearnable). A continuum i∈[0,1]i\in[0,1]i∈[0,1] of investors have CARA utility with risk aversion γ\gammaγ. Investor iii observes:

  1. a private signal s~i=v~+ε~i\tilde s_i = \tilde v + \tilde\varepsilon_is~i​=v~+ε~i​, with ε~i∼N(0,τε−1)\tilde\varepsilon_i \sim N(0,\tau_\varepsilon^{-1})ε~i​∼N(0,τε−1​);
  2. a personalized price signal s~p,i=s~p+u~+e~i\tilde s_{p,i} = \tilde s_p + \tilde u + \tilde e_is~p,i​=s~p​+u~+e~i​, where s~p\tilde s_ps~p​ is the information in the price, u~∼N(0,τu−1)\tilde u \sim N(0,\tau_u^{-1})u~∼N(0,τu−1​) is a common interpretation error and e~i∼N(0,τe−1)\tilde e_i \sim N(0,\tau_e^{-1})e~i​∼N(0,τe−1​) an idiosyncratic one;
  3. an interpretation signal z~i=u~+e~i+η~i\tilde z_i = \tilde u + \tilde e_i + \tilde\eta_iz~i​=u~+e~i​+η~​i​, with η~i∼N(0,τη−1)\tilde\eta_i \sim N(0,\tau_\eta^{-1})η~​i​∼N(0,τη−1​), where τη\tau_\etaτη​ is the common sophistication level.

All shocks are independent and centred, and all precisions and γ\gammaγ are positive. A linear price function is p~=avv~+auu~\tilde p = a_v\tilde v + a_u\tilde up~​=av​v~+au​u~; when av≠0a_v\ne 0av​=0 it is equivalent to the signal s~p=v~+αu~\tilde s_p = \tilde v + \alpha\tilde us~p​=v~+αu~ with α=au/av\alpha = a_u/a_vα=au​/av​. Each investor demands Di=(E(V~∣s~i,s~p,i,z~i)−p~)/(γ Var(V~∣s~i,s~p,i,z~i))D_i = (E(\tilde V\mid \tilde s_i,\tilde s_{p,i},\tilde z_i) - \tilde p)/(\gamma\,\mathrm{Var}(\tilde V\mid \tilde s_i,\tilde s_{p,i},\tilde z_i))Di​=(E(V~∣s~i​,s~p,i​,z~i​)−p~​)/(γVar(V~∣s~i​,s~p,i​,z~i​)), computed by Bayes' rule from the jointly Gaussian signals. A linear equilibrium is a pair (av,au)(a_v,a_u)(av​,au​), av≠0a_v\neq0av​=0, such that aggregate demand ∫01Di di\int_0^1 D_i\,di∫01​Di​di is zero for every realization of (v~,u~)(\tilde v,\tilde u)(v~,u~).

The paper's closed forms use the precision of the combined price information,

τp=τeτu+(τe+τu)τητu+τe(α+1)2+α2τη,\tau_p = \frac{\tau_e\tau_u + (\tau_e+\tau_u)\tau_\eta}{\tau_u + \tau_e(\alpha+1)^2 + \alpha^2\tau_\eta},τp​=τu​+τe​(α+1)2+α2τη​τe​τu​+(τe​+τu​)τη​​,

and the cubic

f(α)=(τeτε+τετη)α3+2τeτεα2+(τeτε+τeτη+τuτε)α−τeτu.f(\alpha) = (\tau_e\tau_\varepsilon + \tau_\varepsilon\tau_\eta)\alpha^3 + 2\tau_e\tau_\varepsilon\alpha^2 + (\tau_e\tau_\varepsilon + \tau_e\tau_\eta + \tau_u\tau_\varepsilon)\alpha - \tau_e\tau_u.f(α)=(τe​τε​+τε​τη​)α3+2τe​τε​α2+(τe​τε​+τe​τη​+τu​τε​)α−τe​τu​.

Formalization targets

Goal: Proposition 1 and the noise remark

There is exactly one linear equilibrium (av,au)(a_v,a_u)(av​,au​). In it

av=τε+τpτv+τε+τp,au=τpτv+τε+τp⋅τu(τe+ατe+ατη)τeτu+τeτη+τuτη,a_v = \frac{\tau_\varepsilon+\tau_p}{\tau_v+\tau_\varepsilon+\tau_p},\qquad a_u = \frac{\tau_p}{\tau_v+\tau_\varepsilon+\tau_p}\cdot\frac{\tau_u(\tau_e+\alpha\tau_e+\alpha\tau_\eta)}{\tau_e\tau_u+\tau_e\tau_\eta+\tau_u\tau_\eta},av​=τv​+τε​+τp​τε​+τp​​,au​=τv​+τε​+τp​τp​​⋅τe​τu​+τe​τη​+τu​τη​τu​(τe​+ατe​+ατη​)​,

with τp\tau_pτp​ evaluated at α=au/av\alpha = a_u/a_vα=au​/av​; α\alphaα lies in (0,αˉ)(0,\bar\alpha)(0,αˉ), αˉ=τeτu/(τeτε+τeτη+τuτε)\bar\alpha = \tau_e\tau_u/(\tau_e\tau_\varepsilon+\tau_e\tau_\eta+\tau_u\tau_\varepsilon)αˉ=τe​τu​/(τe​τε​+τe​τη​+τu​τε​), and is the unique positive root of fff; and au>0a_u > 0au​>0.

Milestones, in the order of the paper's proof

  1. The Bayes coefficients (βs,βp,βz)(\beta_s,\beta_p,\beta_z)(βs​,βp​,βz​) of E(v~∣⋅)E(\tilde v\mid\cdot)E(v~∣⋅) in closed form (proof of Proposition 1, p. 55).
  2. Equation (17): Var(V~∣⋅)=(τv+τε+τp)−1+τξ−1\mathrm{Var}(\tilde V\mid\cdot) = (\tau_v+\tau_\varepsilon+\tau_p)^{-1} + \tau_\xi^{-1}Var(V~∣⋅)=(τv​+τε​+τp​)−1+τξ−1​ (p. 18).
  3. Market clearing holds iff av=βs+βpa_v = \beta_s+\beta_pav​=βs​+βp​ and au=βp(α+1)+βza_u = \beta_p(\alpha+1)+\beta_zau​=βp​(α+1)+βz​ (p. 56).
  4. These sums equal the closed forms of ava_vav​, aua_uau​ (p. 56).
  5. α av(α)=au(α)\alpha\,a_v(\alpha) = a_u(\alpha)αav​(α)=au​(α) iff f(α)=0f(\alpha)=0f(α)=0 (p. 56).
  6. f(0)<0<f(αˉ)f(0) < 0 < f(\bar\alpha)f(0)<0<f(αˉ), hence a root in (0,αˉ)(0,\bar\alpha)(0,αˉ) (p. 56).
  7. The discriminant of fff is the printed 25-term expression Δ<0\Delta<0Δ<0, hence fff has one real root (p. 56).

Significance

Proposition 1 is the foundation of the rest of the paper. Every comparative static on price informativeness, volatility, return predictability, disagreement and trading volume (Propositions 2–3), and the analysis of the sophistication-acquisition game (Propositions 4–6), is computed from the closed forms of ava_vav​, aua_uau​ and the root α\alphaα. The remark au>0a_u>0au​>0 is the paper's mechanism: the common interpretation error is a source of price noise generated by the traders themselves, which vanishes only in the limit of full sophistication (Corollary 1).

The result is proved in the paper; nothing here is open. To our knowledge no part of it, nor any noisy-REE equilibrium with Bayesian price learning, has been machine-checked. The mission produces a formal account of the equilibrium construction of the noisy REE literature, from Gaussian conditioning to market clearing to a polynomial fixed point. The same pattern (conjecture a linear price, condition, match coefficients, solve for the fixed point) recurs throughout that literature.

Difficulty

Each step is elementary, but each is a large symbolic computation. Inverting the 3×33\times33×3 covariance matrix and simplifying Σ−1c\Sigma^{-1}cΣ−1c to the printed fractions involves rational functions in seven variables. The discriminant of fff is a polynomial of degree ten with 25 terms, and its sign has to be established for all positive parameters.

The conceptual pitfall is uniqueness. Showing that fff has one positive root is not enough: a linear equilibrium may a priori have α=au/av≤0\alpha = a_u/a_v \le 0α=au​/av​≤0, so uniqueness among all linear equilibria needs fff to have one real root, which is what the discriminant gives. The fixed-point condition α=au(α)/av(α)\alpha = a_u(\alpha)/a_v(\alpha)α=au​(α)/av​(α) also needs av(α)≠0a_v(\alpha)\neq0av​(α)=0 at every candidate, which has to be proved rather than assumed.

Formalization scope

All quantities are real. A parameter structure carries τv,τξ,τε,τe,τu,τη,γ\tau_v,\tau_\xi,\tau_\varepsilon,\tau_e,\tau_u,\tau_\eta,\gammaτv​,τξ​,τε​,τe​,τu​,τη​,γ; the standing assumption is that all seven are positive, and nothing else is assumed. The model is encoded in reduced form:

  • Conditioning. The conditional mean coefficients are the linear projection Σ(α)−1c\Sigma(\alpha)^{-1}cΣ(α)−1c, with Σ(α)\Sigma(\alpha)Σ(α) the covariance matrix of (s~i,s~p,i,z~i)(\tilde s_i,\tilde s_{p,i},\tilde z_i)(s~i​,s~p,i​,z~i​) and ccc their covariance with v~\tilde vv~. The conditional variance is τv−1−c⊤Σ−1c+τξ−1\tau_v^{-1} - c^\top\Sigma^{-1}c + \tau_\xi^{-1}τv−1​−c⊤Σ−1c+τξ−1​. No probability space is built.
  • Market clearing. It is stated for every (v~,u~)∈R2(\tilde v,\tilde u)\in\mathbb R^2(v~,u~)∈R2 rather than almost surely, which is equivalent for linear functions of a nondegenerate Gaussian pair. Aggregate demand is the demand at the cross-sectional mean (v~, v~+(α+1)u~, u~)(\tilde v,\ \tilde v+(\alpha+1)\tilde u,\ \tilde u)(v~, v~+(α+1)u~, u~) of the signals, the usual continuum convention.
  • The equilibrium predicate. It requires av≠0a_v\neq0av​=0, as on p. 17 of the paper.
  • Equation (4). As printed, it attaches the law N(0,τη−1)N(0,\tau_\eta^{-1})N(0,τη−1​) to z~i\tilde z_iz~i​; (5) and (15) show it belongs to η~i\tilde\eta_iη~​i​, which is the reading used.

The equilibrium predicate is defined from Bayesian demands and market clearing only. Defining it by the closed forms of ava_vav​, aua_uau​, or as "α\alphaα is a root of (18)", would make Proposition 1 a statement about a cubic and is ruled out. Likewise the Bayes coefficients are Σ−1c\Sigma^{-1}cΣ−1c, never the printed fractions.

A complete development needs 3×33\times33×3 matrix inversion over R\mathbb RR (Mathlib's Matrix.inv), rational-function simplification, Mathlib's Cubic and its discriminant, and the intermediate value theorem. A reusable lemma "a real cubic with negative discriminant has exactly one real root" would serve beyond this mission. Proofs of any milestone, and alternative proofs of uniqueness, are welcome.

Selected references

  • X. Vives and L. Yang, Costly Interpretation of Asset Prices, SSRN working paper 2763813, September 2018 version (the version cited throughout). https://ssrn.com/abstract=2763813
  • J. Mondria, X. Vives and L. Yang, Costly Interpretation of Asset Prices, Management Science, 2021. https://doi.org/10.1287/mnsc.2020.3871
  • S. J. Grossman and J. E. Stiglitz, On the Impossibility of Informationally Efficient Markets, American Economic Review 70(3), 1980. https://www.jstor.org/stable/1805228
  • M. F. Hellwig, On the Aggregation of Information in Competitive Markets, Journal of Economic Theory 22(3), 1980. https://doi.org/10.1016/0022-0531(80)90065-5
9 thms1 active userReviewed
OptimizationProbability·Captain: mikedeng1

Stochastic Knapsack Revisited: The Service Level Perspective: Anticipative Policies Need at Most 2Δ More Capacity Than Adaptive Policies to Meet Type-I Service LevelsResearch Paper

Motivation

A central planner holds a single pool of QQQ units of a product and must supply NNN customers whose demands d~1,…,d~N\tilde d_1, \dots, \tilde d_Nd~1​,…,d~N​ are random. Each customer iii asks to be fully served, i.e. to receive her entire demand, with probability at least a target αi\alpha_iαi​. This Type-I service level (the persistency value of the stochastic knapsack) is the standard requirement in supply-chain allocation and in humanitarian logistics. The planner's question is how much capacity the targets require.

The answer depends on when the planner learns the demands. Lyu, Chou, Teo, Zheng and Zhong (SSRN 3674463, 2021) compare three information regimes:

  • a responsive policy sees every demand before allocating;
  • an adaptive policy learns a customer's demand only once it has decided to serve her, and may choose the next customer in light of what it has seen (Dean, Goemans and Vondrák, 2008);
  • an anticipative policy fixes a priority order before any demand is realized (Alptekinoğlu, Banerjee, Paul and Jain, M&SOM 15(1), 2013).

Anticipative rules are the ones used in practice, because they are simple to announce and to audit. The paper's main result quantifies what they cost in capacity compared with adaptive rules.

Timeline. Alptekinoğlu et al. (2013) introduced capacity pooling with differentiated service probabilities under anticipative priority rules. Dean, Goemans and Vondrák (2008) bounded the adaptivity gap of the stochastic knapsack for revenue maximization by a constant factor. Zhong, Zheng, Chou and Teo (Management Science 64(4), 2018) treated the Type-II (fill-rate) version, in which partial service counts in proportion. They showed that its achievable region is a polymatroid Γ(Q)\Gamma(Q)Γ(Q) and that anticipative policies are optimal there. Lyu et al. (2021) use this Type-II system to compare the Type-I classes.

Setting

Customers are i∈{1,…,N}i \in \{1, \dots, N\}i∈{1,…,N}. The demands are random variables on a probability space (Ω,F,P)(\Omega, \mathcal F, \mathbf P)(Ω,F,P). An allocation run with capacity Q≥0Q \ge 0Q≥0 serves customers one at a time. A chosen customer is fully served if her demand does not exceed the remaining capacity, and her demand is then deducted. Otherwise she is not fully served, and the run stops.

  • A deterministic adaptive policy chooses the next customer as a function of the customers chosen so far and their realized demands.
  • A deterministic anticipative policy is a priority order σ\sigmaσ, fixed in advance.
  • A randomized policy in either class is a finite mixture ∑kwkAk\sum_k w_k A_k∑k​wk​Ak​ of deterministic ones, drawn independently of the demands.

The optimal Type-I capacities are

QIAdap=inf⁡{Q≥0:some randomized adaptive policy has P(i fully served)≥αi ∀i},Q^{\mathrm{Adap}}_I = \inf\{Q \ge 0 : \text{some randomized adaptive policy has } \mathbf P(i \text{ fully served}) \ge \alpha_i\ \forall i\},QIAdap​=inf{Q≥0:some randomized adaptive policy has P(i fully served)≥αi​ ∀i},

and QIAntiQ^{\mathrm{Anti}}_IQIAnti​ is defined in the same way with anticipative policies. These are the optimal values of problem (CPP-A) of the paper.

On the Type-II side, Γ(Q)\Gamma(Q)Γ(Q) is the set of α\alphaα with ∑i∈UαiE[d~i]≤E[min⁡(Q,∑i∈Ud~i)]\sum_{i \in U} \alpha_i \mathbf E[\tilde d_i] \le \mathbf E[\min(Q, \sum_{i\in U} \tilde d_i)]∑i∈U​αi​E[d~i​]≤E[min(Q,∑i∈U​d~i​)] for every U⊆{1,…,N}U \subseteq \{1,\dots,N\}U⊆{1,…,N}. QIIQ_{II}QII​ is the least Q≥0Q \ge 0Q≥0 with α∈Γ(Q)\alpha \in \Gamma(Q)α∈Γ(Q). For a priority order σ\sigmaσ, the greedy fill xi∗(σ,d~,Q)=min⁡{d~i,max⁡(0,Q−∑j ahead of id~j)}x^*_i(\sigma, \tilde d, Q) = \min\{\tilde d_i, \max(0, Q - \sum_{j \text{ ahead of } i} \tilde d_j)\}xi∗​(σ,d~,Q)=min{d~i​,max(0,Q−∑j ahead of i​d~j​)} solves the fractional knapsack (CPP-F).

Standing assumptions (§3.3): the demands are independent and positive, and almost surely max⁡id~i≤Δ\max_i \tilde d_i \le \Deltamaxi​d~i​≤Δ for a constant Δ>0\Delta > 0Δ>0. The targets satisfy 0≤αi≤10 \le \alpha_i \le 10≤αi​≤1.

Formalization targets

Goal: Theorem 5

QIAnti−QIAdap≤2Δ.Q^{\mathrm{Anti}}_I - Q^{\mathrm{Adap}}_I \le 2\Delta.QIAnti​−QIAdap​≤2Δ.

The bound is uniform in NNN and in the joint law of the demands, subject only to the standing assumptions.

The two halves: Theorems 3 and 4

QIAnti≤QII+Δ,QIAdap≥QII−Δ.Q^{\mathrm{Anti}}_I \le Q_{II} + \Delta, \qquad Q^{\mathrm{Adap}}_I \ge Q_{II} - \Delta.QIAnti​≤QII​+Δ,QIAdap​≥QII​−Δ.

Milestones inside the proofs

  1. Achievability of Γ(Q)\Gamma(Q)Γ(Q) in fill rate by randomized priority orders (p. 21, from Zhong et al. 2018).
  2. Structure of the greedy fill: it is fractional for at most one customer, and positive exactly when the demand ahead is below QQQ (p. 17).
  3. E[xi∗]/E[d~i]≤P(xi∗>0)\mathbf E[x^*_i]/\mathbf E[\tilde d_i] \le \mathbf P(x^*_i > 0)E[xi∗​]/E[d~i​]≤P(xi∗​>0) (p. 17).
  4. Theorem 3.
  5. The rounded adaptive policy A′A'A′, which also serves the one customer the run turns away, satisfies E[∑k∈Ud~kykA′]≥∑k∈UE[d~k]αk\mathbf E[\sum_{k\in U}\tilde d_k y^{A'}_k] \ge \sum_{k \in U} \mathbf E[\tilde d_k]\alpha_kE[∑k∈U​d~k​ykA′​]≥∑k∈U​E[d~k​]αk​ (p. 37).
  6. Adaptive feasibility of QQQ implies α∈Γ(Q+Δ)\alpha \in \Gamma(Q + \Delta)α∈Γ(Q+Δ) (p. 37). The page writes QIAdptQ^{\mathrm{Adpt}}_IQIAdpt​, a misprint for QIAdapQ^{\mathrm{Adap}}_IQIAdap​.
  7. Theorem 4.

Significance

The result. Theorem 5 says that the information an adaptive planner gains costs at most two maximal demands of capacity, independently of the number of customers. Anticipative priority rules are therefore near-optimal for Type-I service whenever individual demands are small relative to the pool. This contrasts with the responsive regime: Example 1 of the paper shows that QIAdapQ^{\mathrm{Adap}}_IQIAdap​ can exceed QIRespQ^{\mathrm{Resp}}_IQIResp​ by an unbounded factor. In Theorems 3–4 the Type-II capacity QIIQ_{II}QII​ serves as a common yardstick, so each policy class is compared with an explicitly computable polymatroid quantity.

Formalizing it. The theorem is proved on paper. To our knowledge none of it is machine-checked. A formal development has to make precise what the paper leaves informal: what a randomized adaptive policy is, why the event "customer iii is reached" is independent of d~i\tilde d_id~i​, and why the infima defining the optimal capacities are over non-empty sets. The milestones on the greedy fill and the polymatroid Γ(Q)\Gamma(Q)Γ(Q) can be reused for other fill-rate and persistency problems.

Difficulty

The obvious comparison fails because the two classes face different constraints. A Type-I requirement is an all-or-nothing event, so a deterministic priority order cannot hedge. With three identical unit-demand customers and targets 1/31/31/3, any fixed order needs Q=3Q = 3Q=3, while QII=1Q_{II} = 1QII​=1. Randomization over orders is therefore essential on the anticipative side, and the link between a Type-I requirement and the Type-II region Γ(Q)\Gamma(Q)Γ(Q) is not a statement about any single order. On the adaptive side, the decisions depend on revealed demands. Whether customer kkk is fully served is correlated with d~k\tilde d_kd~k​, so a Type-I service probability does not translate directly into a bound on the expected demand served, which is what Γ(Q)\Gamma(Q)Γ(Q) constrains. Conditioning arguments that ignore this correlation give false inequalities.

Formalization scope

  • Demands are a vector d : Fin N → Ω → ℝ on a probability space. Probabilities are μ.real, expectations are Bochner integrals. The standing assumptions are bundled in Standing μ d Δ α: measurable, iIndepFun, almost surely positive, almost surely ≤ Δ, 0 < Δ, and α i ∈ [0, 1].
  • A run has at most NNN steps. A demand equal to the remaining capacity is served. Every run, anticipative ones included, stops at the first customer not fully served.
  • Adaptive policies are functions of the chosen list and the demand vector that depend only on the chosen customers' demands, with measurable decision events. Anticipative policies are permutations, embedded as adaptive policies that ignore the demands.
  • Randomized policies are finite mixtures. Service levels are probabilities, not the long-run scenario frequencies of p. 11.
  • QIAdapQ^{\mathrm{Adap}}_IQIAdap​, QIAntiQ^{\mathrm{Anti}}_IQIAnti​ and QIIQ_{II}QII​ are real infima. Under the standing assumptions every feasible set is non-empty (Q=NΔQ = N\DeltaQ=NΔ works) and bounded below by 000, so no value is a junk default. QIIQ_{II}QII​ uses the subset form Γ(Q)\Gamma(Q)Γ(Q), which p. 17 shows equals the Type-II problem (CPP-B).
  • The greedy fill is stated for a priority order σ\sigmaσ rather than for weights θ′\theta'θ′. It equals x∗(θ′,d~,Q)x^*(\theta', \tilde d, Q)x∗(θ′,d~,Q) for every θ′\theta'θ′ whose decreasing order, ties broken by index, is σ\sigmaσ.
  • Two encodings would trivialize the goal, and both are excluded. Restricting to deterministic priority orders makes Theorem 3 false. Letting adaptive runs skip an oversized customer and continue enlarges the adaptive class beyond the paper's.

The infrastructure the proofs need: measurability of runs, independence of functions of disjoint coordinates, continuity of Q↦E[min⁡(Q,S)]Q \mapsto \mathbf E[\min(Q, S)]Q↦E[min(Q,S)] (so that QIIQ_{II}QII​ is attained), and the greedy description of a polymatroid's base polytope. Proofs of any milestone are welcome. So are alternative proofs of milestone 1 that avoid the general polymatroid theorem.

Selected references

  • G. Lyu, M. C. Chou, C.-P. Teo, Z. Zheng, Y. Zhong, Stochastic Knapsack Revisited: The Service Level Perspective, SSRN working paper 3674463, 2021. https://ssrn.com/abstract=3674463
  • Y. Zhong, Z. Zheng, M. C. Chou, C.-P. Teo, Resource Pooling and Allocation Policies to Deliver Differentiated Service, Management Science 64(4), 1555–1573, 2018.
  • A. Alptekinoğlu, A. Banerjee, A. Paul, N. Jain, Inventory Pooling to Deliver Differentiated Service, Manufacturing & Service Operations Management 15(1), 33–44, 2013.
  • B. C. Dean, M. X. Goemans, J. Vondrák, Approximating the Stochastic Knapsack Problem: The Benefit of Adaptivity, Mathematics of Operations Research 33(4), 945–964, 2008. https://doi.org/10.1287/moor.1080.0330
9 thms1 active userReviewed
Dynamic ProgrammingProbabilityStochastic Systems·Captain: mikedeng1

Discrete Dividend Payments in Continuous Time: The Value Function Is the Unique Fixed Point of the One-Period Dividend-and-Issuance OperatorResearch Paper

Why periodic dividends matter

A firm can return cash to shareholders by paying dividends, but a payment reduces the reserves available to survive future losses. If cash can be distributed only at scheduled dates while equity can be raised between those dates, the firm must coordinate a discrete dividend decision with a continuous-time financing decision. Keppo, Reppen and Soner study this combination for a firm whose cash flow is a Brownian motion with drift. Their preprint identifies the firm's optimal value through an operator that first solves the financing problem between dividend dates and then chooses the next dividend. The mission concerns that characterization, Theorem 3.4, together with the regularity results on which it depends.

The cash reserve model

Let x≥0x\ge0x≥0 be the initial reserve and let Ct=μt+σWtC_t=\mu t+\sigma W_tCt​=μt+σWt​, where WWW is a one-dimensional Brownian motion, μ>0\mu>0μ>0 is profitability, and σ>0\sigma>0σ>0 is volatility. Dividends of amounts ℓn≥0\ell_n\ge0ℓn​≥0 may be paid at dates nTnTnT, with T>0T>0T>0. Equity issues have sizes ζk≥0\zeta_k\ge0ζk​≥0 and occur at stopping times τk\tau_kτk​; the issue times are locally finite. A positive issue costs

c(ζ)=λf+(1+λp)ζ,λf,λp>0,c(\zeta)=\lambda_f+(1+\lambda_p)\zeta,\qquad \lambda_f,\lambda_p>0,c(ζ)=λf​+(1+λp​)ζ,λf​,λp​>0,

while no issue costs zero. The cash reserve after actions at time ttt is Xt=x+Ct−Lt+ItX_t=x+C_t-L_t+I_tXt​=x+Ct​−Lt​+It​. The ruin time θ\thetaθ is the first strictly positive time that Xt<0X_t<0Xt​<0. A dividend cannot exceed the reserve immediately before its date, and no action is taken after ruin. These conditions define the admissible controls Ax\mathcal A_xAx​.

The firm's value function is the supremum, over admissible controls, of expected discounted dividends less expected discounted issue costs:

V(x)=sup⁡(L,I)∈AxEx[∑n≥0e−ρnTℓn−∑ke−ρτkc(ζk)],ρ>0.V(x)=\sup_{(L,I)\in\mathcal A_x}\mathbb E_x\left[\sum_{n\ge0}e^{-\rho nT}\ell_n-\sum_k e^{-\rho\tau_k}c(\zeta_k)\right],\qquad \rho>0.V(x)=(L,I)∈Ax​sup​Ex​[n≥0∑​e−ρnTℓn​−k∑​e−ρτk​c(ζk​)],ρ>0.

The first sum includes a possible dividend at time zero, following the displayed objective in §3.1 of the preprint. An issue at time TTT belongs to the next period. Thus the one-period financing operator Lϕ\mathcal L\phiLϕ charges issues on [0,T)[0,T)[0,T) and evaluates its continuation reward ϕ\phiϕ at the reserve immediately before TTT. The dividend operator is

Dϕ(x)=sup⁡0≤ℓ≤x{ϕ(x−ℓ)+ℓ},T=D∘L.\mathcal D\phi(x)=\sup_{0\le\ell\le x}\{\phi(x-\ell)+\ell\},\qquad \mathcal T=\mathcal D\circ\mathcal L.Dϕ(x)=0≤ℓ≤xsup​{ϕ(x−ℓ)+ℓ},T=D∘L.

The continuation class CA\mathcal C_ACA​ consists of the continuous, nondecreasing functions ϕ:[0,∞)→[0,∞)\phi:[0,\infty)\to[0,\infty)ϕ:[0,∞)→[0,∞) that satisfy the issuance boundary condition of Definition 3.2, are bounded above by x+Ax+Ax+A, and have ϕ(x)−x\phi(x)-xϕ(x)−x nondecreasing and bounded in absolute value by AAA. Its larger class C^A\widehat{\mathcal C}_ACA​ omits the last condition. The intervention value in the boundary condition is Iϕ(x)=sup⁡ζ>0{ϕ(x+ζ)−c(ζ)}\mathcal I\phi(x)=\sup_{\zeta>0}\{\phi(x+\zeta)-c(\zeta)\}Iϕ(x)=supζ>0​{ϕ(x+ζ)−c(ζ)}.

Formalization targets

The first target is the mapping statement in Theorem 3.3: D:C^A→CA\mathcal D:\widehat{\mathcal C}_A\to\mathcal C_AD:CA​→CA​ for A≥0A\ge0A≥0. With

A∗=μTeρT−1,A^*=\frac{\mu T}{e^{\rho T}-1},A∗=eρT−1μT​,

the next targets are the one-period bound 0≤Lϕ(x)≤e−ρTx+A0\le\mathcal L\phi(x)\le e^{-\rho T}x+A0≤Lϕ(x)≤e−ρTx+A for A≥A∗A\ge A^*A≥A∗, the classical one-period characterization of Theorem 3.6, and the remaining mapping and contraction assertions of Theorem 3.3. In particular, T\mathcal TT contracts the uniform distance between functions in CA\mathcal C_ACA​ by the factor e−ρTe^{-\rho T}e−ρT.

The goal is Theorem 3.4:

V∈CA∗,TV=V,ϕ∈CA∗ and Tϕ=ϕ ⟹ ϕ=V.V\in\mathcal C_{A^*},\qquad \mathcal T V=V,\qquad \phi\in\mathcal C_{A^*}\ \text{and}\ \mathcal T\phi=\phi\ \Longrightarrow\ \phi=V.V∈CA∗​,TV=V,ϕ∈CA∗​ and Tϕ=ϕ ⟹ ϕ=V.

This identifies the fixed point with the value defined from admissible controls. It also asserts that the value is finite at every nonnegative initial reserve.

What the characterization gives

Theorem 3.4 turns an infinite-horizon control problem into repeated one-period calculations. Since Theorem 3.3 makes T\mathcal TT a strict contraction on CA∗\mathcal C_{A^*}CA∗​, iterating it determines the value in that class. The source connects this characterization to a numerical method when the one-period financing problem can be computed efficiently. A proof in Lean would also provide reusable statements about locally finite impulse controls, periodic actions, ruin, and discounted extended-real objectives. The preprint proves the mathematical results; this mission records their statements as open Lean theorems for machine-checked proofs.

Where the difficulty lies

The dividend operator is an elementary supremum over one payment, but the financing operator ranges over adapted impulses before ruin. Its regularity is central: the paper invokes Theorem 3.6 to show that Lϕ\mathcal L\phiLϕ is continuous and satisfies the issuance boundary condition. That theorem characterizes the arbitrary-maturity issuance value as the unique classical solution of a parabolic equation with a nonlocal boundary condition. The fixed-point conclusion alone would leave an additional gap: Theorem 3.4 must identify that fixed point with the original control supremum, including the limit over dividend periods and its terminal reward.

Formalization scope

Lean represents nonnegative reserves by R≥0\mathbb R_{\ge0}R≥0​ and the discount, costs and function values by real numbers. Control values and one-period rewards are suprema in extended real numbers, with dividends and issue costs integrated separately as nonnegative quantities. The L\mathcal LL bound proves the finiteness needed before converting to real values. An expression with both expected dividends and expected costs infinite takes the extended-real value −∞-\infty−∞ and cannot improve a supremum with a finite admissible alternative.

The paper's cumulative processes are represented by nonnegative dividend amounts and a locally finite impulse sequence of nonnegative issues. This makes explicit that equity cannot be supplied by an uncharged continuous component. Simultaneous impulses are allowed; merging them cannot increase the cost. The Brownian cash flow uses the published AtaBarjesteh23.Workload.IsDriftBM definition. The model requires a probability measure, adaptation to a filtration, and positive ρ,μ,σ,λf,λp,T\rho,\mu,\sigma,\lambda_f,\lambda_p,Tρ,μ,σ,λf​,λp​,T. The filtration is right-continuous and its initial σ-algebra contains all subsets of ambient null measurable sets, as required by the usual conditions. A dividend at time zero follows the objective's sum. In L\mathcal LL, the continuation value uses XT−X_{T-}XT−​, resolving the free-issue-at-TTT ambiguity in the printed formula. The arbitrary-maturity problem retains its separately printed non-strict survival indicator.

The value VVV is defined solely from admissible controls, never from T\mathcal TT or its iterates. Theorem 3.6 uses time to maturity and states smoothness only on the interior rectangle; uniqueness is restricted to the closed nonnegative time-state rectangle. Reusable work includes stopping-time measurability, impulse summation and parabolic regularity. Contributions proving the goal or the four listed milestones are within scope; §4's random-profitability model and §3.4's numerical experiments are outside this mission.

Selected references

  • J. Keppo, A. M. Reppen and H. M. Soner, Discrete Dividend Payments in Continuous Time, Mathematics of Operations Research, 2021; preprint arXiv:1805.05077v2, 2019. https://arxiv.org/abs/1805.05077v2
7 thms1 active userReviewed
PreviousPage 62 of 97Next

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