Prove2Me
Navigate
DiscoverCollectionsFormalpediaBlogsUsersMomentumMy Missions+
Prove2Me
⌕
Log in
← Collections

Revenue Management and Choice Models

Dynamic pricing, assortment optimization, discrete choice models, and airline seat control.

23 completed missions

Missions

1–20 of 23
OpenCompletedAll
🏆Completed
Operations ResearchOptimization·Captain: naimengye

The Theory and Practice of Revenue Management I: Single-Resource Capacity ControlTextbook

Which fares to open, and when to close them

An airline sells one flight, a hotel one night, a car-rental firm one day of one car: a fixed capacity, perishable at a deadline, sold to customers who arrive over time and are willing to pay different amounts. Chapter 2 of Talluri and van Ryzin's The Theory and Practice of Revenue Management (2004) is the theory of that single resource. Its three models answer the same question with increasing generality: Littlewood's two-class rule, the nnn-class static model and its dynamic-arrival version give the seller a protection level per class, a booking limit or a bid price, all three equivalent; the discrete-choice model, in which customers buy down when a cheaper fare is open, replaces classes by offer sets and shows that only the efficient sets, ordered by their purchase probability, are ever offered, with a higher set the more capacity or the less time remains. This mission formalizes that last result, Theorem 2.3, together with the structural results of the two earlier models that it generalizes.

Setting

Static model. Classes 1,…,n1, \dots, n1,…,n with prices p1≥⋯≥pn≥0p_1 \ge \dots \ge p_n \ge 0p1​≥⋯≥pn​≥0 arrive in stages, lowest class first, with demands DjD_jDj​ distributed on N\mathbb NN. With xxx units left at stage jjj the seller observes DjD_jDj​ and accepts u≤min⁡{Dj,x}u \le \min\{D_j, x\}u≤min{Dj​,x} units; the value function is the Bellman equation (2.3), Vj(x)=E[max⁡u{pju+Vj−1(x−u)}]V_j(x) = \mathbb E[\max_u \{p_j u + V_{j-1}(x-u)\}]Vj​(x)=E[maxu​{pj​u+Vj−1​(x−u)}], V0=0V_0 = 0V0​=0 (staticValue), and ΔVj(x)=Vj(x)−Vj(x−1)\Delta V_j(x) = V_j(x) - V_j(x-1)ΔVj​(x)=Vj​(x)−Vj​(x−1) is the marginal value of capacity. The protection level yj∗=max⁡{x:pj+1<ΔVj(x)}y_j^* = \max\{x : p_{j+1} < \Delta V_j(x)\}yj∗​=max{x:pj+1​<ΔVj​(x)} (protLevel), the booking limit bj∗=C−yj−1∗b_j^* = C - y_{j-1}^*bj∗​=C−yj−1∗​ (bookLimit) and the bid price πj+1(x)=ΔVj(x)\pi_{j+1}(x) = \Delta V_j(x)πj+1​(x)=ΔVj​(x) (bidPrice) define the three controls of Theorem 2.1.

Dynamic model. Over TTT periods at most one request arrives per period, of class jjj with probability λj(t)\lambda_j(t)λj​(t); the value function (2.17) is Vt(x)=Vt+1(x)+E[max⁡u∈{0,1}(R(t)−ΔVt+1(x))u]V_t(x) = V_{t+1}(x) + \mathbb E[\max_{u \in \{0,1\}} (R(t) - \Delta V_{t+1}(x))u]Vt​(x)=Vt+1​(x)+E[maxu∈{0,1}​(R(t)−ΔVt+1​(x))u] (dynValue), with time-dependent protection levels (2.19), booking limits (2.20) and bid prices (2.18).

Choice model. When the set SSS of classes is open an arriving customer buys class j∈Sj \in Sj∈S with probability Pj(S)P_j(S)Pj​(S); Q(S)=∑j∈SPj(S)Q(S) = \sum_{j \in S} P_j(S)Q(S)=∑j∈S​Pj​(S) is the purchase probability and R(S)=∑j∈SPj(S)pjR(S) = \sum_{j \in S} P_j(S) p_jR(S)=∑j∈S​Pj​(S)pj​ the expected revenue (purchaseProb, expRevenue). The value function (2.26) is Vt(x)=max⁡Sλt(R(S)−Q(S)ΔVt+1(x))+Vt+1(x)V_t(x) = \max_{S} \lambda_t (R(S) - Q(S)\Delta V_{t+1}(x)) + V_{t+1}(x)Vt​(x)=maxS​λt​(R(S)−Q(S)ΔVt+1​(x))+Vt+1​(x) (choiceValue). A set TTT is inefficient (Definition 2.1, IsInefficient) if a randomization α\alphaα over the subsets has ∑Sα(S)Q(S)≤Q(T)\sum_S \alpha(S) Q(S) \le Q(T)∑S​α(S)Q(S)≤Q(T) and ∑Sα(S)R(S)>R(T)\sum_S \alpha(S) R(S) > R(T)∑S​α(S)R(S)>R(T), and efficient otherwise.

Formalization targets

Goal: Theorem 2.3

In every period with capacity left, some efficient set maximizes (2.26); and, the efficient sets being ordered by QQQ, the largest optimal set is nondecreasing in the remaining capacity xxx and nondecreasing in the period ttt: choice_optimal_policy. Monotonicity is stated as "every efficient optimal set at (t,x)(t, x)(t,x) is matched by one at (t,x′)(t, x')(t,x′), x′≥xx' \ge xx′≥x, with at least as large a purchase probability", and likewise in ttt.

Supporting targets

Littlewood's rule (2.1), ΔV1(x)=p1P(D1≥x)\Delta V_1(x) = p_1 \mathbb P(D_1 \ge x)ΔV1​(x)=p1​P(D1​≥x) and the acceptance criterion; Proposition 2.1, the marginal values of the static model are decreasing in xxx and increasing in the stages remaining; Theorem 2.1, nested protection levels, nested booking limits and bid-price tables each attain the Bellman maximum at every stage; Proposition 2.2 and Theorem 2.2, the same two results for the dynamic model, with marginal values now decreasing in time; Proposition 2-2.A.4 of the appendix, the marginal values of the choice model are decreasing in xxx and in ttt; Proposition 2.3, an inefficient set is never optimal; and the ordering of efficient sets, Q(S)≤Q(S′)Q(S) \le Q(S')Q(S)≤Q(S′) implies R(S)≤R(S′)R(S) \le R(S')R(S)≤R(S′) when S′S'S′ is efficient.

The continuous-demand optimality conditions (2.9) of Sect. 2.2.2.3, stated without proof, the computational and heuristic methods of Sects. 2.2.3-2.2.4, the overbooking models of Sect. 2.7 and the nested-policy characterization of Sect. 2.6.2.5 are not targets.

Significance

Theorem 2.3 is the structural result behind choice-based revenue management: it reduces the 2n2^n2n offer sets to the efficient frontier of (Q(S),R(S))(Q(S), R(S))(Q(S),R(S)), orders that frontier, and shows the optimal policy walks up it as capacity grows or the deadline nears. It was the analytical core of Talluri and van Ryzin's (2004) choice-model paper and is the reason the efficient sets, not the fare classes, are the unit of control when customers substitute between fares. The static and dynamic results, from Littlewood (1972) and Brumelle and McGill (1993) to Lee and Hersh (1993), are the foundation of every airline seat inventory control system; the equivalence of protection levels, booking limits and bid prices is what lets the same optimal policy be implemented on any of the three kinds of reservation system. None of these results has a machine-checked proof.

Difficulty

The two marginal-value propositions are inductions in which the inductive step is the discrete concavity of a max-plus convolution, Lemma 2-2.A.1 of the appendix: x↦max⁡0≤a≤m{ap+g(x−a)}x \mapsto \max_{0 \le a \le m}\{ap + g(x-a)\}x↦max0≤a≤m​{ap+g(x−a)} is concave when ggg is, which in Lean requires reasoning about the argmax on N\mathbb NN and the truncated subtraction. The static model's expectation is a tsum against a pmf, so every step also needs summability of a bounded family. The protection-level theorems then need the down-set structure of {x:pj+1<ΔVj(x)}\{x : p_{j+1} < \Delta V_j(x)\}{x:pj+1​<ΔVj​(x)} under monotonicity of ΔVj\Delta V_jΔVj​, and the three controls have to be shown to coincide unit by unit. For the choice model, Proposition 2.3 is a one-line convexity argument once ΔV≥0\Delta V \ge 0ΔV≥0 is known, and the monotonicity in Theorem 2.3 is a monotone comparative-statics argument on the objective R(S)−Q(S)ΔR(S) - Q(S)\DeltaR(S)−Q(S)Δ, which is easy in Δ\DeltaΔ but must be combined with Proposition 2-2.A.4 in both xxx and ttt; the existence of an efficient maximizer uses Proposition 2.3 and the finiteness of the subsets.

Formalization scope

Capacities, stages and periods are natural numbers, the value functions recurse on the stage or on the number of periods to go, and the book's ranges (x≤Cx \le Cx≤C, t≤Tt \le Tt≤T, j≤nj \le nj≤n) are hypotheses of the theorems. Demand in the static model is a pmf on N\mathbb NN rather than a random variable, so the expectation in (2.3) is a tsum; the dynamic model's expectation over R(t)R(t)R(t) is written out, including the no-arrival term, which vanishes under nonnegative prices. The choice model is defined by its compact form (2.26), and the maximization includes the empty offer set. Optimality of a control means attaining the inner maximum of the Bellman equation at every state, which is what the book's proofs establish. The bid-price control is formalized with the bid price πj+1(x+1−z)\pi_{j+1}(x + 1 - z)πj+1​(x+1−z) of the zzz-th unit allocated; the book prints x−zx - zx−z, which is one unit off from (2.5). The appendix's Proposition 2-2.A.4 prints its time monotonicity in the reverse direction; the formal statement is the direction consistent with Proposition 2.2 and Theorem 2.3. The ordering of efficient sets is stated with non-strict inequalities, since Definition 2.1 admits ties in revenue.

Selected references

  • K. T. Talluri and G. J. van Ryzin, The Theory and Practice of Revenue Management, Kluwer/Springer, 2004, Chapter 2. https://doi.org/10.1007/b139000
  • K. Littlewood, Forecasting and control of passenger bookings, AGIFORS Symposium Proceedings 12, 1972; reprinted in Journal of Revenue and Pricing Management 4(2), 2005. https://doi.org/10.1057/palgrave.rpm.5170134
  • S. L. Brumelle and J. I. McGill, Airline seat allocation with multiple nested fare classes, Operations Research 41(1), 1993. https://doi.org/10.1287/opre.41.1.127
  • T. C. Lee and M. Hersh, A model for dynamic airline seat inventory control with multiple seat bookings, Transportation Science 27(3), 1993. https://doi.org/10.1287/trsc.27.3.252
  • K. T. Talluri and G. J. van Ryzin, Revenue management under a general discrete choice model of consumer behavior, Management Science 50(1), 2004. https://doi.org/10.1287/mnsc.1030.0147
  • C. J. Lautenbacher and S. Stidham, The underlying Markov decision process in the single-leg airline yield-management problem, Transportation Science 33(2), 1999. https://doi.org/10.1287/trsc.33.2.136
11 thms4 active usersReviewed
🏆Completed
Operations ResearchOptimization·Captain: naimengye

The Theory and Practice of Revenue Management II: OverbookingTextbook

How far to oversell

Every airline, hotel and car-rental firm sells more reservations than it has capacity, because some customers cancel or do not show. Chapter 4 of Talluri and van Ryzin's The Theory and Practice of Revenue Management (2004) is the theory of that decision. Its static models pick one overbooking limit from the show distribution; its dynamic model, a simplification of Chatwin's (1998), follows reservations, cancellations and refunds period by period and proves that the optimal control is still a limit, one that declines toward the deadline and falls when more demand is expected; and its substitutable-capacity model, from Karaesmen and van Ryzin (2004), sets joint limits for several classes whose oversold customers can be moved between resources, showing the expected net revenue is concave in each limit and submodular across them. This mission formalizes the chapter's four propositions and the corollary the book draws from the last.

Setting

Dynamic overbooking (Sect. 4.3.1). With yyy reservations on hand in period ttt, DtD_tDt​ new requests arrive; the firm books up to x∈[y,y+Dt]x \in [y, y + D_t]x∈[y,y+Dt​] at revenue p(t)p(t)p(t) each, and every reservation survives the period with probability qtq_tqt​, a cancellation refunding r(t)r(t)r(t). At the deadline T+1T + 1T+1 the firm pays the convex denied-service cost c(y−C)c(y - C)c(y−C) on reservations beyond capacity CCC, Eq. (4.11). The recursion is vt+1(x)=E[Vt+1(Zt(x))−(x−Zt(x)) r(t)]v_{t+1}(x) = \mathbb E[V_{t+1}(Z_t(x)) - (x - Z_t(x))\,r(t)]vt+1​(x)=E[Vt+1​(Zt​(x))−(x−Zt​(x))r(t)] with Zt(x)∼Bin(x,qt)Z_t(x) \sim \mathrm{Bin}(x, q_t)Zt​(x)∼Bin(x,qt​) and Vt(y)=E[max⁡y≤x≤y+Dt{vt+1(x)+(x−y) p(t)}]V_t(y) = \mathbb E[\max_{y \le x \le y + D_t}\{v_{t+1}(x) + (x - y)\,p(t)\}]Vt​(y)=E[maxy≤x≤y+Dt​​{vt+1​(x)+(x−y)p(t)}] (value, postValue). The greatest optimal overbooking limit x∗(t)x^*(t)x∗(t) (overbookingLimit) is the largest level at which vt+1(x)+x p(t)v_{t+1}(x) + x\,p(t)vt+1​(x)+xp(t) is at least its value at every smaller level, an element of N∪{∞}\mathbb N \cup \{\infty\}N∪{∞}; the limit policy books min⁡{y+Dt,max⁡{y,x∗}}\min\{y + D_t, \max\{y, x^*\}\}min{y+Dt​,max{y,x∗}} (limitPolicy).

Substitutable capacity (Sect. 4.5). Classes j=1,…,nj = 1, \dots, nj=1,…,n hold yjy_jyj​ reservations and are overbooked to levels xjx_jxj​; in the service period Zj∼Poisson(qjxj)Z_j \sim \mathrm{Poisson}(q_j x_j)Zj​∼Poisson(qj​xj​) customers show and are assigned to resources i=1,…,mi = 1, \dots, mi=1,…,m of capacities CiC_iCi​, or to the virtual resource 000 (denied service), at net benefit hjih_{ji}hji​, by the transportation problem (TP) with value V(z,C)V(z, C)V(z,C) (serviceValue). The expected net revenue (4.21) is G(x)=p⊤(x−y)−E[s⊤(x−Z(x))]+E[V(Z(x),C)]G(x) = p^\top(x - y) - \mathbb E[s^\top(x - Z(x))] + \mathbb E[V(Z(x), C)]G(x)=p⊤(x−y)−E[s⊤(x−Z(x))]+E[V(Z(x),C)] (expNetRevenue), and jointLimit is the greatest optimal limit of one class with the others fixed.

Formalization targets

Goal: Proposition 4.4

With Poisson show demands, GGG has decreasing first differences in every direction: G(x+ei+ej)−G(x+ei)≤G(x+ej)−G(x)G(x + e_i + e_j) - G(x + e_i) \le G(x + e_j) - G(x)G(x+ei​+ej​)−G(x+ei​)≤G(x+ej​)−G(x) for all xxx and all classes i,ji, ji,j, which is component-wise concavity (i=ji = ji=j) and submodularity (i≠ji \ne ji=j): joint_overbooking_concave_submodular.

Supporting targets

Proposition 4.1, with a convex denied-service cost the limit policy with the greatest optimal limit attains the maximum of the recursion at every state; Proposition 4.2, under qt(p(t)−p(t+1))+(1−qt)(p(t)−r(t))≥0q_t(p(t) - p(t+1)) + (1 - q_t)(p(t) - r(t)) \ge 0qt​(p(t)−p(t+1))+(1−qt​)(p(t)−r(t))≥0 the greatest optimal limits decline with time; Proposition 4.3, stochastically larger demand to come gives limits that are no larger; and the corollary of Sect. 4.5.2, the greatest optimal limit of class iii is nonincreasing in the level of any other class.

The static overbooking models of Sect. 4.2 (binomial, normal and Gram-Charlier approximations, Type 1 and Type 2 service levels), the net-bookings heuristics of Sect. 4.3.2, the combined capacity-control models of Sect. 4.4 and the stochastic-gradient algorithm of Appendix 4.A carry no numbered results and are not targets.

Significance

Proposition 4.4 is the structural fact that makes joint overbooking of related resources tractable: concavity gives each class a critical booking level and submodularity makes those levels move in opposite directions, so a stochastic-gradient or coordinate search on the limits is well behaved, and the pattern of Example 4.5, overbooking an early flight aggressively because its oversold passengers can be moved to later ones, is a consequence rather than a heuristic. The dynamic propositions are the theoretical support for the overbooking curves that reservation systems post, limits that fall as departure approaches, and they quantify the sense in which a static model, which ignores future demand, overbooks too much. The proof of Proposition 4.4 passes through the discrete concavity of the transportation problem's value in its supply vector, an M-natural-concavity fact in the sense of Murota, and the Poisson-expectation identity for second differences; none of this has a machine-checked proof.

Difficulty

The dynamic model needs the concavity of VtV_tVt​ on N\mathbb NN to be propagated through two operations, the binomial thinning x↦E[V(Bin(x,q))]x \mapsto \mathbb E[V(\mathrm{Bin}(x, q))]x↦E[V(Bin(x,q))] and the windowed maximum y↦max⁡y≤x≤y+Dg(x)y \mapsto \max_{y \le x \le y + D} g(x)y↦maxy≤x≤y+D​g(x), both of which preserve discrete concavity but require explicit manipulation of binomial sums and of the argmax; Propositions 4.2 and 4.3 then compare greatest maximizers of concave sequences through lower bounds on marginal values, with the value ∞\infty∞ handled in ℕ∞. The substitutable-capacity goal is harder: the value of (TP) as a function of the integer supply vector must be shown to have decreasing differences, which is the submodularity of a max-weight transportation value in its supplies, a linear programming duality argument (or Murota's M-natural-concavity of min-cost flow), and the Poisson expectation of it, a tsum over Nn\mathbb N^nNn, must be differenced in two coordinates using the identity E[f(Nμ+δ)]−E[f(Nμ)]\mathbb E[f(N_{\mu + \delta})] - \mathbb E[f(N_\mu)]E[f(Nμ+δ​)]−E[f(Nμ​)] for Poisson pmfs. The linear terms of GGG cancel in second differences and the refund term is linear in xxx.

Formalization scope

Periods are natural numbers with value t the value with T+1−tT + 1 - tT+1−t periods to go, and the book's ranges 1≤t≤T1 \le t \le T1≤t≤T are hypotheses. The denied-service cost is normalized, c(0)=0c(0) = 0c(0)=0 and c≥0c \ge 0c≥0, as a cost "penalizing denied service" is. Convexity of the sequence alone is not enough, because (4.11) never reads c(0)c(0)c(0). Demands are pmfs on N\mathbb NN and cancellations exact binomial sums. The greatest optimal limit lives in N∪{∞}\mathbb N \cup \{\infty\}N∪{∞} because a mild denied-service cost can make accepting every request optimal, in which case the book's critical value is +∞+\infty+∞; the limit policy then accepts everything. Proposition 4.3 is stated for two demand families ordered by first-order stochastic dominance rather than a parametrized family. In the substitutable-capacity model the virtual resource is uncapacitated, the book's "finite but very high" C0C_0C0​ taken as infinite so that (TP) is feasible for every Poisson realization, and (TP) is over real assignments, whose optimum at integer supplies is integral. Eq. (4.21) is printed with −E[V(Z(x),C)]-\mathbb E[V(Z(x), C)]−E[V(Z(x),C)]; VVV being the maximum net benefit, the expected net revenue adds it, and the definition uses +++, without which Proposition 4.4 fails numerically on every sampled instance.

Selected references

  • K. T. Talluri and G. J. van Ryzin, The Theory and Practice of Revenue Management, Kluwer/Springer, 2004, Chapter 4. https://doi.org/10.1007/b139000
  • R. E. Chatwin, Multiperiod airline overbooking with a single fare class, Operations Research 46(6), 1998. https://doi.org/10.1287/opre.46.6.805
  • I. Karaesmen and G. J. van Ryzin, Overbooking with substitutable inventory classes, Operations Research 52(1), 2004. https://doi.org/10.1287/opre.1030.0079
  • M. Rothstein, OR and the airline overbooking problem, Operations Research 33(2), 1985. https://doi.org/10.1287/opre.33.2.237
  • K. Murota, Discrete Convex Analysis, SIAM, 2003. https://doi.org/10.1137/1.9780898718508
7 thms3 active usersReviewed
🏆Completed
Operations ResearchOptimization·Captain: naimengye

The Theory and Practice of Revenue Management III: Dynamic PricingTextbook

Prices that respond to inventory

A retailer marking down a seasonal line, an airline raising fares as seats sell, a manufacturer pricing while restocking: each sets prices over time against a finite and changing inventory. Chapter 5 of Talluri and van Ryzin's The Theory and Practice of Revenue Management (2004) collects the structural theory of that problem. Without replenishment, the Bernoulli-arrival model of Gallego and van Ryzin (1994) gives a marginal value of capacity that falls with inventory and with time, hence prices that jump up at each sale and drift down between sales, and the deterministic fluid model bounds it from above. With replenishment, the model of Federgruen and Heching (1999) has a jointly concave and supermodular continuation value, from which the base-stock, posted-price policy follows: below a base-stock level order up to it and post a fixed price, above it order nothing and discount, more the higher the inventory. This mission formalizes those results with Proposition 5.3 as the goal.

Setting

Bernoulli demand (Sect. 5.2.2.2). One customer arrives per period with a random willingness to pay; the firm's decision is the demand rate d∈[0,1]d \in [0, 1]d∈[0,1], the probability of a sale at the inverse-demand price p(t,d)p(t, d)p(t,d), with revenue rate r(t,d)=d p(t,d)r(t, d) = d\,p(t, d)r(t,d)=dp(t,d) (revenueRate). The value function (5.12) is Vt(x)=max⁡d{r(t,d)−d ΔVt+1(x)}+Vt+1(x)V_t(x) = \max_{d}\{r(t, d) - d\,\Delta V_{t+1}(x)\} + V_{t+1}(x)Vt​(x)=maxd​{r(t,d)−dΔVt+1​(x)}+Vt+1​(x), VT+1=0V_{T+1} = 0VT+1​=0, Vt(0)=0V_t(0) = 0Vt​(0)=0 (bernoulliValue), and ΔVt(x)=Vt(x)−Vt(x−1)\Delta V_t(x) = V_t(x) - V_t(x-1)ΔVt​(x)=Vt​(x)−Vt​(x−1) (bernoulliDelta). The deterministic model (5.1) maximizes ∑tr(t,d(t))\sum_t r(t, d(t))∑t​r(t,d(t)) over rates with ∑td(t)≤C\sum_t d(t) \le C∑t​d(t)≤C (deterministicValue).

Pricing with replenishment (Sect. 5.3.2). Inventory may be negative (backorders). In period ttt with inventory xxx the firm orders up to y≥xy \ge xy≥x at unit cost ctc_tct​, chooses a rate d∈[0,dˉ]d \in [0, \bar d]d∈[0,dˉ], sells the random demand D(t,d,ξt)=at(ξt) d+bt(ξt)D(t, d, \xi_t) = a_t(\xi_t)\,d + b_t(\xi_t)D(t,d,ξt​)=at​(ξt​)d+bt​(ξt​) (additive, multiplicative or mixed noise), and pays the convex cost hth_tht​ on ending inventory. The value function (5.20) is Vt(x)=sup⁡y≥x, d{r(t,d)−ct(y−x)+Gt+1(y,d)}V_t(x) = \sup_{y \ge x,\, d}\{r(t, d) - c_t(y - x) + G_{t+1}(y, d)\}Vt​(x)=supy≥x,d​{r(t,d)−ct​(y−x)+Gt+1​(y,d)} with the continuation value Gt+1(y,d)=E[Vt+1(y−D(t,d,ξt))−ht(y−D(t,d,ξt))]G_{t+1}(y, d) = \mathbb E[V_{t+1}(y - D(t, d, \xi_t)) - h_t(y - D(t, d, \xi_t))]Gt+1​(y,d)=E[Vt+1​(y−D(t,d,ξt​))−ht​(y−D(t,d,ξt​))] (ReplPricing.value, contValue). Assumption 7.2, marginal revenue decreasing, is the concavity of r(t,⋅)r(t, \cdot)r(t,⋅).

Formalization targets

Goal: Proposition 5.3

For every period, Gt+1G_{t+1}Gt+1​ is jointly concave on R×[0,dˉ]\mathbb R \times [0, \bar d]R×[0,dˉ], VtV_tVt​ is concave on R\mathbb RR, and Gt+1G_{t+1}Gt+1​ has increasing differences in (y,d)(y, d)(y,d), the supermodularity that the book states through its partial derivatives: replenishment_concave_supermodular.

Supporting targets

Proposition 5.2, the marginal value of capacity in the Bernoulli model decreases in ttt and in xxx; the upper bound of Sect. 5.2.2.3, the optimal deterministic revenue dominates the optimal expected stochastic revenue; and the base-stock, posted-price structure of Sect. 5.3.2.1, derived from Proposition 5.3: below the unconstrained optimum y0y^0y0 order up to it and use d0d^0d0, above it order nothing and use a rate at least d0d^0d0 that is nondecreasing in the inventory.

Proposition 5.1 and Lemma 5-5.A.1 (the continuous-demand model without replenishment) are not targets; see the formalization scope. The deterministic sections (efficient prices, discrete price sets), the asymptotic optimality of the deterministic heuristic, the infinite-horizon stationary problem and the multiproduct and finite-population models carry no numbered results.

Significance

Proposition 5.3 is the structural core of joint pricing and inventory control: joint concavity makes the period problem a concave program, and supermodularity is what turns its solution into a policy, the base-stock, posted-price rule that Federgruen and Heching showed optimal and that later work on pricing with inventory builds on. Proposition 5.2 is the reason optimal dynamic prices in the stochastic single-item model rise at every sale and fall while inventory sits, the behaviour of Figure 5.5, and the deterministic upper bound is what justifies the fluid model as a benchmark and a heuristic, the pattern quantified in Table 5.6. None of these has a machine-checked proof; the replenishment result in particular needs the interplay of concavity, expectation and partial maximization on all of R\mathbb RR.

Difficulty

Proposition 5.2 is an induction whose step compares suprema over the rate interval, with the boundary condition Vt(0)=0V_t(0) = 0Vt​(0)=0 breaking the recursion at x=1x = 1x=1 and requiring r(t,0)=0r(t, 0) = 0r(t,0)=0. The deterministic bound is an induction on periods that uses the concavity of the deterministic value in the inventory (a concave program's value) to absorb the two branches of the Bernoulli recursion. The goal needs: integrability and continuity of Vt+1(y−D)−ht(y−D)V_{t+1}(y - D) - h_t(y - D)Vt+1​(y−D)−ht​(y−D) under bounded noise; that the expectation of a concave function of an affine map is jointly concave, and its increasing differences from those of the concave integrand; that a partial supremum of a jointly concave function over the convex feasible set {y≥x}\{y \ge x\}{y≥x} is concave in xxx; and the boundedness of the objective so that every supremum is a real number. The base-stock item is the segment argument that moves an unconstrained maximizer onto the boundary y=xy = xy=x and a monotone comparative-statics argument on the supermodular objective, with maxima attained by continuity on the compact rate interval.

Formalization scope

Periods are natural numbers with value t the value with T+1−tT + 1 - tT+1−t periods to go, the maxima are suprema, and the book's ranges are hypotheses. The demand is affine in the rate, which is the additive and multiplicative models the book names; with a merely convex demand (Assumption 5.1) the joint concavity of Proposition 5.3 fails when Vt+1−htV_{t+1} - h_tVt+1​−ht​ is not monotone, and the noise has bounded support, strengthening Assumption 7.6. The partial-derivative statements (iii)-(iv) are in difference form. Proposition 5.1 is not formalized: its model (5.11) evaluates Vt+1(x−D)V_{t+1}(x - D)Vt+1​(x−D) at negative inventories the model does not define while truncating revenue at xxx, and its Lemma 5-5.A.1 (joint concavity of r+r^+r+) is false as stated, its Hessian argument mistaking an indefinite matrix for a negative definite one; a counterexample is in the mission's check. The deterministic model restricts rates to [0,1][0, 1][0,1], the rates the Bernoulli model can realize.

Selected references

  • K. T. Talluri and G. J. van Ryzin, The Theory and Practice of Revenue Management, Kluwer/Springer, 2004, Chapter 5. https://doi.org/10.1007/b139000
  • G. Gallego and G. J. van Ryzin, Optimal dynamic pricing of inventories with stochastic demand over finite horizons, Management Science 40(8), 1994. https://doi.org/10.1287/mnsc.40.8.999
  • A. Federgruen and A. Heching, Combined pricing and inventory control under uncertainty, Operations Research 47(3), 1999. https://doi.org/10.1287/opre.47.3.454
  • W. Elmaghraby and P. Keskinocak, Dynamic pricing in the presence of inventory considerations, Management Science 49(10), 2003. https://doi.org/10.1287/mnsc.49.10.1287.17315
  • D. M. Topkis, Supermodularity and Complementarity, Princeton University Press, 1998. https://doi.org/10.1515/9781400822539
5 thms2 active usersReviewed
🏆Completed
Mechanism DesignOperations Research·Captain: naimengye

The Theory and Practice of Revenue Management IV: AuctionsTextbook

Why a reserve price, and why it does not matter which auction

Airlines selling last seats, Priceline's name-your-own-price, procurement of supply contracts: Chapter 6 of Talluri and van Ryzin's The Theory and Practice of Revenue Management (2004) treats auctions as pricing mechanisms and asks what revenue they earn and how to design them. Its centre is Myerson's (1981) theory for independent private values: whatever the mechanism, so long as bidders with higher valuations are more likely to win and the lowest type gains nothing, the firm's expected revenue is the expected virtual value ∑iJ(vi)yi(v)\sum_i J(v_i) y_i(v)∑i​J(vi​)yi​(v) of the winners, with J(v)=v−(1−F(v))/f(v)J(v) = v - (1 - F(v))/f(v)J(v)=v−(1−F(v))/f(v) (Theorem 6.1, the revenue equivalence theorem). Maximizing that expression pointwise gives the optimal auction: the standard first- or second-price auction with a reserve price v∗v^*v∗ at the zero of JJJ (Theorem 6.2). This mission formalizes the second-price form of Theorem 6.2 as its goal, with the dominant-strategy and first-price equilibria of the informal analysis, Theorem 6.1, the optimal allocation and Proposition 6.1 on list prices as supporting results.

Setting

NNN customers have i.i.d. valuations on [0,vˉ][0, \bar v][0,vˉ] with a continuously differentiable, strictly increasing distribution FFF and positive density fff (PrivateValues, IsRegular); the joint law is the product measure (joint). A direct-revelation mechanism (Mechanism) maps reported valuations to allocations yi(v)∈{0,1}y_i(v) \in \{0, 1\}yi​(v)∈{0,1}, at most CCC units in total, and payments pi(v)p_i(v)pi​(v). For a report www by customer iii, Pi(w)P_i(w)Pi​(w) is the win probability, Ri(w)R_i(w)Ri​(w) the expected payment and Si(w)=wPi(w)−Ri(w)S_i(w) = w P_i(w) - R_i(w)Si​(w)=wPi​(w)−Ri​(w) the surplus (winProb, expPayment, expSurplus); incentive compatibility, Si(w)≥wPi(w′)−Ri(w′)S_i(w) \ge w P_i(w') - R_i(w')Si​(w)≥wPi​(w′)−Ri​(w′), is the equilibrium condition of the direct mechanism (IsIncentiveCompatible). The chapter's mechanisms are the CCC-unit second-price auction with reserve price rrr (secondPriceReserve: the CCC highest valuations above rrr win and pay the larger of rrr and the highest losing valuation), the list-price mechanism for N≤CN \le CN≤C (listPrice), and the single-unit first-price auction with its equilibrium bid b∗(v)=v−∫0vP(s) ds/P(v)b^*(v) = v - \int_0^v P(s)\,ds / P(v)b∗(v)=v−∫0v​P(s)ds/P(v), P=FN−1P = F^{N-1}P=FN−1 (firstPriceBid).

Formalization targets

Goal: Theorem 6.2

With JJJ strictly increasing (Assumption 7.2) and v∗v^*v∗ its zero, the CCC-unit second-price auction with reserve price v∗v^*v∗ is a feasible, incentive-compatible mechanism with monotone allocations and zero surplus at zero, and its expected revenue is at least that of every such mechanism: reserve_price_auction_optimal.

Supporting targets

Bidding one's valuation is dominant in the second-price auction (Sect. 6.2.2.1); the bid (6.4) solves the first-order condition (6.3), is a symmetric equilibrium of the first-price auction and shades below the valuation (Sect. 6.2.2.2); Theorem 6.1, revenue equals expected virtual surplus and each expected payment is wPi(w)−∫0wPiw P_i(w) - \int_0^w P_iwPi​(w)−∫0w​Pi​; the pointwise optimal allocation of Sect. 6.2.5; and Proposition 6.1, a list price at v∗v^*v∗ is optimal when N≤CN \le CN≤C.

Proposition 6.2 (asymptotic optimality of list prices, a law-of-large-numbers statement about scaled auctions), the first-price form of Theorem 6.2 with its equilibrium (6.9) stated without proof, and the dynamic, replenishment and network auctions of Sects. 6.3-6.5 (Propositions 6.3-6.11, from Vulcano, van Ryzin and Maglaras and from Cooper and Menich) are not targets of this mission.

Significance

Theorem 6.1 is the tool that lets revenue be computed from allocations alone, which is why the first- and second-price auctions of Examples 6.1-6.3 earn the same (N−1)/(N+1)(N-1)/(N+1)(N−1)/(N+1) and why any dynamic pricing scheme that ends with the same winners earns the same as the optimal auction (Sect. 6.2.6.3). Theorem 6.2 says a firm with private-value customers cannot do better than a standard auction with the right reserve price, and Proposition 6.1 that with enough capacity a list price already does it: auctions are a small-numbers phenomenon. These are the foundations on which the chapter's dynamic auctions and the list-price comparisons of Sects. 6.3-6.4 rest, and Myerson's optimal auction has no machine-checked proof in its multi-unit form.

Difficulty

Theorem 6.1 is an envelope argument in measure-theoretic clothing: incentive compatibility gives the two-sided inequalities of Appendix 6.A, monotonicity of PiP_iPi​ makes SiS_iSi​ convex with derivative PiP_iPi​ almost everywhere, so Si(w)=∫0wPiS_i(w) = \int_0^w P_iSi​(w)=∫0w​Pi​, and then an integration by parts against the density converts ∫(wPi(w)−Si(w))f(w) dw\int (w P_i(w) - S_i(w)) f(w)\,dw∫(wPi​(w)−Si​(w))f(w)dw into ∫J(w)Pi(w)f(w) dw\int J(w) P_i(w) f(w)\,dw∫J(w)Pi​(w)f(w)dw; the win probabilities are integrals over a product measure with one coordinate replaced, and Fubini is needed to return to E[J(vi)yi(v)]\mathbb E[J(v_i) y_i(v)]E[J(vi​)yi​(v)]. The goal then needs the reserve-price auction shown incentive compatible (a dominant-strategy argument on the threshold payment), measurable, monotone and with zero surplus at zero, and the pointwise optimal allocation integrated. The first-price item is calculus on an interval integral with a vanishing denominator at 000 and a monotone comparative-statics argument for the equilibrium inequality.

Formalization scope

Mechanisms are direct-revelation mechanisms on [0,vˉ]N[0, \bar v]^N[0,vˉ]N, as the book reduces to in Sect. 6.2.3.1; expectations over the other customers are integrals over the joint law with customer iii's coordinate overwritten by the report. Payments are assumed bounded on reports in [0,vˉ]N[0, \bar v]^N[0,vˉ]N (not on all of RN\mathbb R^NRN, where the second-price payment is unbounded) and the rules measurable. Ties in the second-price auction are broken by index, a null event, and when every customer wins the losing supremum is 000 so the winner pays the reserve. Theorem 6.2 is stated for the second-price auction; the first-price version with reserve price, whose equilibrium (6.9) the book asserts without proof, is left out and noted. Optimality is over mechanisms satisfying conditions (i) and (ii) of Theorem 6.1 and incentive compatibility, which is the class the book compares against. The virtual value's zero v∗v^*v∗ is a parameter with J(v∗)=0J(v^*) = 0J(v∗)=0 rather than the maximum of (6.8), which under strict monotonicity is the same point.

Selected references

  • K. T. Talluri and G. J. van Ryzin, The Theory and Practice of Revenue Management, Kluwer/Springer, 2004, Chapter 6. https://doi.org/10.1007/b139000
  • R. B. Myerson, Optimal auction design, Mathematics of Operations Research 6(1), 1981. https://doi.org/10.1287/moor.6.1.58
  • J. G. Riley and W. F. Samuelson, Optimal auctions, American Economic Review 71(3), 1981. https://www.jstor.org/stable/1802786
  • P. Klemperer, Auction theory: a guide to the literature, Journal of Economic Surveys 13(3), 1999. https://doi.org/10.1111/1467-6419.00083
  • W. Vickrey, Counterspeculation, auctions, and competitive sealed tenders, Journal of Finance 16(1), 1961. https://doi.org/10.1111/j.1540-6261.1961.tb02789.x
  • E. Maskin and J. Riley, Optimal multi-unit auctions, in The Economics of Missing Markets, Information, and Games, Oxford University Press, 1989.
7 thms3 active usersReviewed
🏆Completed
Algorithmic Game TheoryOperations Research·Captain: naimengye

The Theory and Practice of Revenue Management V: CompetitionTextbook

When competing firms settle on prices and allocations

Revenue management is practised by firms that compete: airlines matching fares while allocating seats, retailers ordering stock and then pricing to clear it, hotels protecting rooms for late high-paying guests. Chapter 8 of Talluri and van Ryzin's The Theory and Practice of Revenue Management (2004) surveys the economics behind these situations: monopoly pricing and mechanism design, the Coase problem, advance-purchase discounts, and oligopoly in quantities (Cournot), in prices (Bertrand, Bertrand–Edgeworth with capacities), in newsvendor capacities and in RM allocations. Most of its numbered results are quoted from the literature; what the book itself establishes is a family of equilibrium arguments built on Talluri's equilibrium graph: when each firm's best response moves monotonically with the rival's action, following best-response arcs must end in an equilibrium pair. This mission formalizes those arguments, with Proposition 8.3, existence of an equilibrium in offer sets under the multinomial-logit choice model, as the goal.

Setting

A two-firm game on finite chains of strategies has payoffs u1,u2u_1, u_2u1​,u2​, best responses and pure Nash equilibria (IsBestResponse1, IsNashEquilibrium); its equilibrium graph has no crossing arcs (NoCrossing1) when best-response arcs (k2,k1)(k_2, k_1)(k2​,k1​), (l2,l1)(l_2, l_1)(l2​,l1​) with k2<l2k_2 < l_2k2​<l2​ always have k1≤l1k_1 \le l_1k1​≤l1​. The chapter's games are: the RM duopoly of Sect. 8.4.1.3, two firms with capacity CCC and fares pL<pHp_L < p_HpL​<pH​ whose high-fare demand spills over to the rival beyond its protection level (spilloverDemand), each firm answering with Littlewood's protection level (littlewoodResponse); the duopoly newsvendor of Example 8.16 with effective demands R1=D1+(D2−x2)+R_1 = D_1 + (D_2 - x_2)^+R1​=D1​+(D2​−x2​)+ (effectiveDemand); linear Cournot (cournotPayoff); Bertrand–Edgeworth price competition with capacity CCC per firm, linear demand and efficient rationing (beSales, bePayoff, IsBEEquilibrium); the advance-purchase model of Sect. 8.3.6 (peakLoadRevenue, advancePurchaseRevenue); and the one-period offer-set game of Sect. 8.4.3.2, where a firm offering its first kkk products earns g(Ck)/(W1(k)+W2(l)+w0)g(C_k)/(W_1(k) + W_2(l) + w_0)g(Ck​)/(W1​(k)+W2​(l)+w0​) with g(Ck)=∑j≤kwj(pj−Δ)−w0δg(C_k) = \sum_{j \le k} w_j(p_j - \Delta) - w_0\deltag(Ck​)=∑j≤k​wj​(pj​−Δ)−w0​δ (OfferFirm, offerPayoff1, CaseI, CaseII).

Formalization targets

Goal: Proposition 8.3

If both firms are in Case I (g(G∗)≥0g(G^*) \ge 0g(G∗)≥0) or both in Case II (g<0g < 0g<0 on every complete set), the offer-set game has a pure-strategy equilibrium in complete sets: mnl_offer_set_equilibrium.

Supporting targets

The equilibrium-graph lemma, monotone best responses for both firms give an equilibrium (Sect. 8.4.1.3); Proposition 8.1, Littlewood responses in the RM duopoly are monotone in the rival's protection level, and the resulting equilibrium of the RM duopoly game; Example 8.16, duopoly newsvendor capacities total at least the monopoly capacity; Example 8.15, the symmetric Cournot equilibrium and its price; Theorem 8.5 (i)-(ii), the pure-strategy Bertrand–Edgeworth equilibria; and the advance-purchase comparison of Sect. 8.3.6, w2<w^w_2 < \hat ww2​<w^ and VAPD(w^)≥Vpeak(w2)V_{APD}(\hat w) \ge V_{peak}(w_2)VAPD​(w^)≥Vpeak​(w2​).

Not targets: Theorems 8.1-8.2 (the Coase problem, subgame-perfect equilibria of an infinite-horizon game, proofs in von der Fehr and Kühn), Theorem 8.3 (Harris–Raviv priority pricing, stated with a garbled price formula), Theorem 8.4 (existence via quasiconcavity, cited), Theorem 8.5 (iii), 8.6 and 8.7-8.8 (mixed-strategy and supergame equilibria of Kreps and Scheinkman and Benoit and Krishna), Proposition 8.2 and Proposition 8.4 (the dynamic offer-set game under condition (8.32), proved in Talluri's paper), and the Kreps–Scheinkman derivation of Sect. 8.4.1.6, which rests on Theorem 8.6.

Significance

The equilibrium graph is the chapter's own contribution: a bipartite picture of best responses on chains in which monotonicity, the absence of crossing arcs, forces an equilibrium. It is the finite, combinatorial form of the monotone comparative-statics route to Nash equilibrium (Tarski's fixed point on a chain), and it is what makes RM allocation games tractable: the two-class duopoly with Littlewood responses always has an equilibrium, and the offer-set duopoly does whenever the two firms face the same sign of ggg, while Example 8.18 shows a best-response cycle when they do not. The Bertrand–Edgeworth, Cournot and newsvendor results are the benchmarks the book uses to interpret RM competition: capacity constraints soften Bertrand's zero-profit outcome, competing in allocations tends to raise total capacity above the monopoly level, and advance-purchase discounts dominate peak-load pricing as a self-selection mechanism. None of these has a machine-checked proof.

Difficulty

The lemma and the goal are fixed-point arguments on finite chains: the largest best response is a monotone map of the rival's index, the composition of two monotone (or two antitone) maps on a finite chain has a fixed point, and for the offer-set game the monotonicity itself must be extracted from the ratio structure g(Ck)/(W(k)+a)g(C_k)/(W(k) + a)g(Ck​)/(W(k)+a) as in the appendix's inequalities (8.A.3)-(8.A.5), separately in the two cases. Proposition 8.1 is a monotonicity of tail probabilities under the pointwise order of effective demands. Example 8.16 is a short probabilistic argument that needs the identity {D>x1+x2}={R1>x1}∩{R2>x2}\{D > x_1 + x_2\} = \{R_1 > x_1\} \cap \{R_2 > x_2\}{D>x1​+x2​}={R1​>x1​}∩{R2​>x2​} and the strict monotonicity of the tail. Theorem 8.5 requires computing efficient-rationing sales for every unilateral deviation from a symmetric profile, through the recursive definition of residual demand, and a quadratic inequality for upward deviations in part (ii). The advance-purchase item reduces to the identity VAPD(w)−Vpeak(w)=(1−α)wV_{APD}(w) - V_{peak}(w) = (1 - \alpha)wVAPD​(w)−Vpeak​(w)=(1−α)w and the strict decrease of the first-order-condition function.

Formalization scope

Products, protection levels and complete sets are natural numbers; the offer-set game's strategies are the complete sets C1,…,CnC_1, \dots, C_nC1​,…,Cn​ of the book, and its payoff is (8.31) up to the positive factor λ\lambdaλ and the terms independent of both offer sets. Prices decreasing in the product index are a hypothesis, as the nested-by-revenue order of Sect. 8.4.3.2. The RM duopoly is modeled through Littlewood's response, as the book's appendix argues, rather than through expected revenues. Efficient rationing is defined recursively over the firms priced strictly below a given firm, with equal sharing among firms at the same price. Example 8.16 is stated with the equilibrium conditions (8.23) and the strictly increasing distribution of aggregate demand as hypotheses. The advance-purchase item takes the first-order conditions (8.13) and (8.15) and the book's uniqueness assumption as hypotheses, the latter as (v−w)−F(w)/f(w)(v - w) - F(w)/f(w)(v−w)−F(w)/f(w) strictly decreasing in www (the page prints "increasing", but its footnote identifies it with the monotone marginal-revenue assumption, which is decreasing in the waiting cost). Theorem 8.5 is stated for its pure-strategy parts (i) and (ii) only.

Selected references

  • K. T. Talluri and G. J. van Ryzin, The Theory and Practice of Revenue Management, Kluwer/Springer, 2004, Chapter 8. https://doi.org/10.1007/b139000
  • D. M. Kreps and J. A. Scheinkman, Quantity precommitment and Bertrand competition yield Cournot outcomes, Bell Journal of Economics 14(2), 1983. https://doi.org/10.2307/3003636
  • S. A. Lippman and K. F. McCardle, The competitive newsboy, Operations Research 45(1), 1997. https://doi.org/10.1287/opre.45.1.54
  • S. Netessine and R. A. Shumsky, Revenue management games: horizontal and vertical competition, Management Science 51(5), 2005. https://doi.org/10.1287/mnsc.1040.0356
  • I. L. Gale and T. J. Holmes, Advance-purchase discounts and monopoly allocation of capacity, American Economic Review 83(1), 1993. https://www.jstor.org/stable/2117500
  • D. M. Topkis, Supermodularity and Complementarity, Princeton University Press, 1998. https://doi.org/10.1515/9781400822539
8 thms2 active usersReviewed
🏆Completed
Algorithmic Game TheoryOperations ResearchOptimization+1·Captain: mikedeng1

Optimal Pricing of Seasonal Products in the Presence of Forward-Looking Consumers 1: Threshold Purchasing Policies under Contingent PricingResearch Paper

Motivation

Retailers of fashion and seasonal goods sell at a premium price early in the season and mark the remaining stock down later. When customers anticipate the markdown, some of them who would buy at the premium price instead wait, trading a lower price against the risk that the item sells out and against the decline of their own valuation over the season. How forward-looking ("strategic") customers respond to a markdown policy is the first question any model of such pricing has to answer, because the seller's optimal prices depend on it.

Aviv and Pazgal (MSOM 2008) model a seller with a fixed inventory, Poisson arrivals of customers with heterogeneous, exponentially declining valuations, and two pricing regimes: contingent pricing, where the discount depends on the inventory left at the markdown time, and announced fixed discounts. The first step of their analysis of contingent pricing is Theorem 1: whatever the other customers do, a customer's best response is a threshold rule on his current valuation, with a threshold that rises as the markdown approaches. Their numerical study of equilibria and of the value of price commitment (§§4.2–7) is built on this reduction.

Setting

A seller holds QQQ units over a season [0,H][0, H][0,H] split at a fixed time TTT with 0<T≤H0 < T \le H0<T≤H. On [0,T)[0, T)[0,T) the premium price p1p_1p1​ applies. At time TTT the seller observes the remaining inventory QT∈{0,1,…,Q}Q_T \in \{0, 1, \dots, Q\}QT​∈{0,1,…,Q} and charges the discount menu price p2(QT)p_2(Q_T)p2​(QT​), where p2(q)≤p1p_2(q) \le p_1p2​(q)≤p1​ for q=1,…,Qq = 1, \dots, Qq=1,…,Q. Customer jjj has a base valuation VjV_jVj​ and valuation Vj(t)=Vje−αtV_j(t) = V_j e^{-\alpha t}Vj​(t)=Vj​e−αt at time ttt, with a common decline factor α≥0\alpha \ge 0α≥0.

A customer arriving at t<Tt < Tt<T either buys immediately at p1p_1p1​ or waits until TTT, when he requests a unit if the discounted price leaves him a nonnegative surplus. Waiting is uncertain in two ways: the remaining inventory QTQ_TQT​ is random, and when fewer units remain than customers request them, units are rationed at random. A belief is a probability mass function π\piπ of QTQ_TQT​ on {0,…,Q}\{0, \dots, Q\}{0,…,Q} together with allocation probabilities a(q)=Pr⁡{A∣QT=q}∈[0,1]a(q) = \Pr\{\mathcal A \mid Q_T = q\} \in [0,1]a(q)=Pr{A∣QT​=q}∈[0,1], a(0)=0a(0) = 0a(0)=0, where A\mathcal AA is the event that the customer is allocated a unit. It is determined by the other customers' strategies, which are arbitrary.

With δ=e−α(T−t)\delta = e^{-\alpha(T-t)}δ=e−α(T−t), the expected surplus of waiting of a customer with current valuation ψ\psiψ is

Wt(ψ)=EQT ⁣[max⁡{ψδ−p2(QT),0}⋅1{A∣QT}]=∑q=0Qπ(q) a(q) max⁡{ψδ−p2(q),0}.W_t(\psi) = \mathrm E_{Q_T}\!\left[\max\{\psi\delta - p_2(Q_T), 0\}\cdot \mathbf 1\{\mathcal A \mid Q_T\}\right] = \sum_{q=0}^{Q}\pi(q)\,a(q)\,\max\{\psi\delta - p_2(q), 0\}.Wt​(ψ)=EQT​​[max{ψδ−p2​(QT​),0}⋅1{A∣QT​}]=q=0∑Q​π(q)a(q)max{ψδ−p2​(q),0}.

The paper's purchase rule (p. 344): buy immediately iff the current surplus V(t)−p1V(t) - p_1V(t)−p1​ is nonnegative and at least Wt(V(t))W_t(V(t))Wt​(V(t)).

Formalization targets

Goal: Theorem 1 and Corollary 1

Assume p1≥0p_1 \ge 0p1​≥0, and α>0\alpha > 0α>0 or ∑qπ(q)a(q)<1\sum_q \pi(q)a(q) < 1∑q​π(q)a(q)<1. For every t∈[0,T)t \in [0,T)t∈[0,T) the equation

ψ−p1=Wt(ψ)(2)\psi - p_1 = W_t(\psi) \tag{2}ψ−p1​=Wt​(ψ)(2)

has a unique solution ψ(t)≥p1\psi(t) \ge p_1ψ(t)≥p1​; a customer arriving at ttt buys immediately under the purchase rule if and only if V(t)≥ψ(t)V(t) \ge \psi(t)V(t)≥ψ(t); and the threshold function ψ:[0,T)→[p1,∞)\psi : [0, T) \to [p_1, \infty)ψ:[0,T)→[p1​,∞) is nondecreasing in ttt.

Milestones

  1. The right-hand side of (2) is nonnegative and nondecreasing in ψ\psiψ, with increments bracketed by δ Pr⁡{ψδ≥p2(QT),A}\delta\,\Pr\{\psi\delta \ge p_2(Q_T), \mathcal A\}δPr{ψδ≥p2​(QT​),A} at the two endpoints, and this slope is below one.
  2. Equation (2) has a unique solution ψ≥p1\psi \ge p_1ψ≥p1​.

Significance

Theorem 1 reduces a customer's strategy, a function of arrival time and valuation, to one threshold function ψ\psiψ on [0,T)[0, T)[0,T). The segment sizes ΛI,ΛS,ΛW,ΛL\Lambda_I, \Lambda_S, \Lambda_W, \Lambda_LΛI​,ΛS​,ΛW​,ΛL​ of §4.2, the seller's menu problem (3), the equilibrium iteration (4) and the closed form of Proposition 2 are all written in terms of ψ\psiψ; without Theorem 1 none of them is defined. Corollary 1, that the threshold rises toward the markdown, is what the paper calls "useful in our analyses below"; the customer segments of Figure 1 are drawn with it.

The result is proved in the paper, with a short appendix argument. No machine-checked version exists. The mission produces a formal statement and proof of the reduction for an arbitrary belief, which fixes the exact hypotheses under which it holds: the paper's slope bound needs either valuation decline (α>0\alpha > 0α>0) or imperfect availability, and the monotonicity of the threshold needs a nonnegative premium price. A formal WtW_tWt​ and threshold are the starting point for formalizing the equilibrium and pricing results of the paper.

Difficulty

The mathematics is one-dimensional. The difficulty is in stating it exactly. WtW_tWt​ is piecewise linear with a kink wherever ψδ\psi\deltaψδ crosses a menu price, so the paper's derivative is only a one-sided derivative, and the uniqueness argument has to use increments. The paper's bound "slope <1< 1<1" is false when α=0\alpha = 0α=0 and a unit is allocated with certainty; then (2) has either no finite solution or a half-line of them. The threshold's monotonicity in ttt rests on Wt(ψ)W_t(\psi)Wt​(ψ) increasing in ttt for fixed ψ\psiψ, which needs ψ≥0\psi \ge 0ψ≥0; with a negative premium price the threshold can decrease. The naive reading of "optimal to use a threshold" as an abstract fixed-point fact about any monotone function with slope below one discards the model and is not the goal.

Formalization scope

Lean namespace SeasonalPricing.Contingent. Time, prices and valuations are real numbers. The belief is a pair pmf alloc : ℕ → ℝ restricted to {0, …, Q} (IsInventoryBelief), not a random variable on a probability space; only the law of (QT,1{A})(Q_T, \mathbf 1\{\mathcal A\})(QT​,1{A}) enters (2). The menu is p2 : ℕ → ℝ with p2(q)≤p1p_2(q) \le p_1p2​(q)≤p1​ required on {1,…,Q}\{1, \dots, Q\}{1,…,Q} only; p2(0)p_2(0)p2​(0) never matters because a(0)=0a(0) = 0a(0)=0. The belief does not depend on the arrival time, as in Eq. (4) of the paper. waitingSurplus is WtW_tWt​ with e−α(T−t)e^{-\alpha(T-t)}e−α(T−t) written Real.exp (-(α * (T - t))); buysNow is the purchase rule, stated on the current valuation V(t)V(t)V(t).

Readings of the paper's words:

  • "the unique solution" of (2): existence and uniqueness of a real ψ≥p1\psi \ge p_1ψ≥p1​ (∃!). The paper's "ψ∈[p1,∞]\psi \in [p_1, \infty]ψ∈[p1​,∞]" includes ∞\infty∞ only in the case excluded by the added hypothesis.
  • "it is optimal to base purchasing decisions on a threshold function": the purchase rule of p. 344 holds exactly when V(t)≥ψ(t)V(t) \ge \psi(t)V(t)≥ψ(t).
  • "derivative … <1< 1<1": a two-sided bracket on increments of WtW_tWt​, with right slope δPr⁡{ψδ≥p2(QT),A}\delta\Pr\{\psi\delta \ge p_2(Q_T), \mathcal A\}δPr{ψδ≥p2​(QT​),A}, below one.
  • "increasing" (Corollary 1): nondecreasing (MonotoneOn), since ψ\psiψ is constant on an initial interval whenever no menu price is reachable (p. 347).

Added hypotheses, both named in the statements: α>0\alpha > 0α>0 or ∑qπ(q)a(q)<1\sum_q \pi(q)a(q) < 1∑q​π(q)a(q)<1, the one hypothesis the paper's proof uses without stating it; and p1≥0p_1 \ge 0p1​≥0, the model's convention that prices are nonnegative. Only the branch 0≤t<T0 \le t < T0≤t<T of the threshold θ\thetaθ is stated: for t≥Tt \ge Tt≥T the paper's θ(t)=p2\theta(t) = p_2θ(t)=p2​ is the model's rule for late customers. The belief enters through the explicit sum; a formalization with an unspecified monotone WWW, or with ψ(t)\psi(t)ψ(t) defined by choice inside a definition, is not the target.

No new library is needed beyond finite sums, max and Real.exp. A lemma on unique roots of ψ↦ψ−c−f(ψ)\psi \mapsto \psi - c - f(\psi)ψ↦ψ−c−f(ψ) for fff with increments bounded by k(ψ′−ψ)k(\psi' - \psi)k(ψ′−ψ), k<1k < 1k<1, is reusable. Proofs of the milestones and the goal, in any order, are welcome.

Selected references

  • Y. Aviv and A. Pazgal, Optimal Pricing of Seasonal Products in the Presence of Forward-Looking Consumers, Manufacturing & Service Operations Management 10(3):339–359, 2008. https://doi.org/10.1287/msom.1070.0183
  • X. Su, Intertemporal Pricing with Strategic Customer Behavior, Management Science 53(5):726–741, 2007. https://doi.org/10.1287/mnsc.1060.0667
  • G. Gallego and G. van Ryzin, Optimal Dynamic Pricing of Inventories with Stochastic Demand over Finite Horizons, Management Science 40(8):999–1020, 1994. https://doi.org/10.1287/mnsc.40.8.999
5 thms2 active usersReviewed
🏆Completed
Algorithmic Game TheoryOperations ResearchOptimization+1·Captain: mikedeng1

Optimal Pricing of Seasonal Products in the Presence of Forward-Looking Consumers 2: A Threshold Nash Equilibrium under Announced Fixed-Discount PricingResearch Paper

Motivation

Retailers of fashion and seasonal goods sell at a premium price early in the season and mark down later. When customers anticipate the markdown, some of them wait, and the seller's pricing problem becomes a game between the seller and a population of forward-looking (strategic) customers. Aviv and Pazgal (MSOM 10(3), 2008) study this game in a model with limited inventory, stochastic arrivals and valuations that decline over the season, under two classes of seller policies: contingent pricing, where the discount depends on the inventory left, and announced fixed-discount pricing, where the seller commits to both prices upfront. Their numerical study (§7.3) compares the two classes and finds that precommitment can raise expected revenue by up to about 8%.

That comparison needs, for every announced price path, the customers' equilibrium response. Theorem 2 of the paper (p. 348) supplies it: a threshold purchasing policy, pinned down by a scalar fixed-point equation for the probability that a waiting customer is served. This mission formalizes Theorem 2. A companion mission of the same series formalizes Theorem 1, the contingent-pricing counterpart.

Setting

A seller has Q≥1Q \ge 1Q≥1 units to sell over a season [0,H][0, H][0,H], split at a fixed time TTT with 0<T≤H0 < T \le H0<T≤H. Customers arrive by a Poisson process with rate λ>0\lambda > 0λ>0. Customer jjj has a base valuation VjV_jVj​ drawn from a continuous distribution FFF (tail Fˉ=1−F\bar F = 1 - FFˉ=1−F), and at time ttt values the product at Vj(t)=Vje−αtV_j(t) = V_j e^{-\alpha t}Vj​(t)=Vj​e−αt, where the decline factor α≥0\alpha \ge 0α≥0 is common to all customers.

Under an announced price path the seller commits to a premium price p1p_1p1​ on [0,T)[0, T)[0,T) and a discount price p2≤p1p_2 \le p_1p2​≤p1​ from TTT on; p2p_2p2​ does not depend on the remaining inventory. Customers know the initial inventory but not the current one.

A customer arriving at t<Tt < Tt<T buys immediately if and only if (i) the current surplus V(t)−p1V(t) - p_1V(t)−p1​ is nonnegative and (ii) it is at least the expected surplus of waiting,

ω⋅max⁡{V(T)−p2,0},\omega\cdot\max\{V(T) - p_2, 0\},ω⋅max{V(T)−p2​,0},

where ω\omegaω is the probability that a unit will be allocated to the customer at time TTT. Units left at TTT are rationed at random among the customers who request one.

For a threshold function ψ\psiψ on [0,T)[0, T)[0,T) the paper defines three segment rates: ΛI(ψ)\Lambda_I(\psi)ΛI​(ψ), the expected number of customers who buy at p1p_1p1​; ΛS(ψ,p1,p2)\Lambda_S(\psi, p_1, p_2)ΛS​(ψ,p1​,p2​), those who could buy at p1p_1p1​ but wait and want to buy at p2p_2p2​; and ΛW(p1,p2)\Lambda_W(p_1, p_2)ΛW​(p1​,p2​), those whose valuation was below p1p_1p1​ and who want to buy at p2p_2p2​. Each is λ\lambdaλ times an integral over [0,T][0, T][0,T] of Fˉ\bar FFˉ at scaled prices (p. 345). With P(x∣Λ)P(x \mid \Lambda)P(x∣Λ) the Poisson probabilities, the allocation probability of qqq units is

A(q∣Λ)=∑y=0∞qmax⁡{1+y,q} P(y∣Λ).A(q \mid \Lambda) = \sum_{y=0}^{\infty} \frac{q}{\max\{1+y, q\}}\,P(y \mid \Lambda).A(q∣Λ)=y=0∑∞​max{1+y,q}q​P(y∣Λ).

Formalization targets

Goal: Theorem 2 (p. 348)

For w∈[0,1]w \in [0,1]w∈[0,1] let

ψA(t)=max⁡{p1,p1−wp21−we−α(T−t)},0≤t<T,(7)\psi_A(t) = \max\left\{p_1, \frac{p_1 - wp_2}{1 - we^{-\alpha(T-t)}}\right\},\qquad 0 \le t < T, \tag{7}ψA​(t)=max{p1​,1−we−α(T−t)p1​−wp2​​},0≤t<T,(7)

and suppose www solves

w=∑x=0Q−1P(x∣ΛI(ψA))⋅A(Q−x∣ΛS(ψA,p1,p2)+ΛW(p1,p2)).(8)w = \sum_{x=0}^{Q-1} P\big(x \mid \Lambda_I(\psi_A)\big)\cdot A\big(Q-x \mid \Lambda_S(\psi_A, p_1, p_2) + \Lambda_W(p_1, p_2)\big). \tag{8}w=x=0∑Q−1​P(x∣ΛI​(ψA​))⋅A(Q−x∣ΛS​(ψA​,p1​,p2​)+ΛW​(p1​,p2​)).(8)

Then, when all other customers use ψA\psi_AψA​ (so that a waiting customer is served with the probability on the right of (8)), every customer arriving at t∈[0,T)t \in [0, T)t∈[0,T) buys immediately if and only if V(t)≥ψA(t)V(t) \ge \psi_A(t)V(t)≥ψA​(t): the symmetric threshold profile is a Nash equilibrium.

Milestones: the two cases of the proof (p. 358)

  1. If e−α(T−t)≤p2/p1e^{-\alpha(T-t)} \le p_2/p_1e−α(T−t)≤p2​/p1​, the threshold is p1p_1p1​.
  2. If e−α(T−t)>p2/p1e^{-\alpha(T-t)} > p_2/p_1e−α(T−t)>p2​/p1​, the threshold is (p1−wp2)/(1−we−α(T−t))≥p1(p_1 - wp_2)/(1 - we^{-\alpha(T-t)}) \ge p_1(p1​−wp2​)/(1−we−α(T−t))≥p1​.

Significance

Theorem 2 reduces the customers' equilibrium under an announced path to a single scalar www. Everything downstream in §5 and §7 rests on it: the seller's expected revenue πA/S(p1,p2)\pi_{A/S}(p_1, p_2)πA/S​(p1​,p2​) (p. 348) is written in terms of ψA\psi_AψA​, the seller's optimal announced path maximizes it, and the comparison between announced and contingent pricing uses the resulting value πA/S∗\pi^*_{A/S}πA/S∗​. The theorem also explains the qualitative prediction of the model: the threshold exceeds p1p_1p1​ exactly when the announced discount is deep relative to the decline of valuations, and it rises with the perceived availability www.

The result is proved in the paper; to the best of our search it has no machine-checked proof. A formal development contributes the model objects (segment rates for threshold policies, the allocation probability for random rationing among Poisson requesters) in a form reusable by the rest of the series and by other strategic-customer pricing models, and a checked proof of the equilibrium property. The existence of a solution to (8) is not proved in the paper and is a natural further target.

Difficulty

The best-response part of the argument is elementary once the availability is known. The substance of the statement lies in the availability itself: the probability that a waiting customer is served is not a free parameter but the one generated, through (8), by the other customers' use of the same threshold. A formalization must connect the segment rates, the Poisson counts and random rationing into one expression and keep the fixed-point coupling between www and ψA\psi_AψA​ intact; dropping it turns the theorem into a one-line inequality about an arbitrary www. The division by 1−we−α(T−t)1 - we^{-\alpha(T-t)}1−we−α(T−t) also degenerates when w=1w = 1w=1 and α=0\alpha = 0α=0, and has to be excluded explicitly.

Formalization scope

The Lean development lives in namespace SeasonalPricing.Announced. Conventions:

  • Time is real; base valuations have law μ : Measure ℝ with IsProbabilityMeasure μ, FFF = ProbabilityTheory.cdf μ, and continuity of FFF (the paper's "continuous distribution") is a hypothesis of the goal. No support condition on [0,∞)[0,\infty)[0,∞) is imposed; the statement quantifies over every real base valuation VVV.
  • ΛI,ΛS,ΛW\Lambda_I, \Lambda_S, \Lambda_WΛI​,ΛS​,ΛW​ are interval integrals over [0,T][0, T][0,T] exactly as printed. P(x∣Λ)=e−ΛΛx/x!P(x \mid \Lambda) = e^{-\Lambda}\Lambda^x/x!P(x∣Λ)=e−ΛΛx/x! is written out; A(q∣Λ)A(q\mid\Lambda)A(q∣Λ) is the infinite series (tsum) as printed, not its closed form.
  • availability is the right-hand side of (8), with ψA\psi_AψA​ built from www by (7).

Readings of the paper's informal words:

  • "Nash equilibrium" is read as the best-response property the paper's proof checks: against the availability generated by (8), the immediate-purchase rule of p. 344 coincides with the threshold ψA\psi_AψA​ at every t∈[0,T)t \in [0, T)t∈[0,T) and every valuation. The paper defines no strategy space beyond threshold rules.
  • "www is a solution to (8)": the theorem is conditional on a solution; its existence is neither assumed elsewhere nor claimed. The conditional statement has content only when (8) has a solution, which the paper does not prove.
  • www as a likelihood: 0≤w≤10 \le w \le 10≤w≤1 is a hypothesis (it also follows from (8)).
  • Added hypothesis: α>0\alpha > 0α>0 or w<1w < 1w<1, which keeps 1−we−α(T−t)>01 - we^{-\alpha(T-t)} > 01−we−α(T−t)>0 for t<Tt < Tt<T; the paper's formula is undefined when it fails. In the milestones the same condition appears as we−α(T−t)<1we^{-\alpha(T-t)} < 1we−α(T−t)<1, and 0<p10 < p_10<p1​ is added so that p2/p1p_2/p_1p2​/p1​ is meaningful.
  • The rule on [T,H][T, H][T,H] (buy at TTT iff V(T)>p2V(T) > p_2V(T)>p2​) is part of the model and is not restated; HHH does not enter the statements.

A formalization in which www is an arbitrary number in [0,1][0,1][0,1], not tied to (8), is ruled out: it is the best-response lemma alone, not Theorem 2. Contributions welcome: proofs of the two milestones and the goal; lemmas such as 0≤A(q∣Λ)≤10 \le A(q\mid\Lambda) \le 10≤A(q∣Λ)≤1 and summability of its series; the closed form of A(q∣Λ)A(q \mid \Lambda)A(q∣Λ) printed on p. 346; and an existence result for (8).

Selected references

  • Y. Aviv and A. Pazgal, Optimal Pricing of Seasonal Products in the Presence of Forward-Looking Consumers, Manufacturing & Service Operations Management 10(3):339–359, 2008. https://doi.org/10.1287/msom.1070.0183
  • G. Gallego and G. van Ryzin, Optimal Dynamic Pricing of Inventories with Stochastic Demand over Finite Horizons, Management Science 40(8):999–1020, 1994. https://doi.org/10.1287/mnsc.40.8.999
  • X. Su, Intertemporal Pricing with Strategic Customer Behavior, Management Science 53(5):726–741, 2007. https://doi.org/10.1287/mnsc.1060.0667
5 thms3 active usersReviewed
🏆Completed
Operations ResearchOptimizationProbability·Captain: mikedeng1

Optimal Pricing of Seasonal Products in the Presence of Forward-Looking Consumers 3: Optimal Contingent-Pricing Revenue with Myopic Customers and Exponential ValuationsResearch Paper

Motivation

Retailers of seasonal goods (fashion, electronics, holiday items) sell a fixed stock over a short season and routinely cut prices toward its end. A markdown of this kind segments the market over time: customers with high valuations buy early at a premium price, and customers with lower valuations are served later at a discount price. Aviv and Pazgal (MSOM 2008) study how much such two-price schemes are worth when customers arrive over time, differ in their valuations, and may or may not anticipate the discount.

To measure the value of price segmentation, the paper compares every two-price scheme with the best fixed-price policy, a single price held for the whole season. Its benchmark is the case of myopic customers, who never delay a purchase strategically. Proposition 3 of the paper computes this benchmark in closed form in the simplest nontrivial setting: exponentially distributed valuations that do not decline over the season, and unlimited inventory. The resulting formula explains the pattern of the paper's Table 1, where the benefit of segmentation grows with the heterogeneity of valuations and with a late discount time.

Setting

A seller offers a product during the season [0,H][0, H][0,H]; throughout this mission H=1H = 1H=1, so time is measured as a fraction of the season. Customers arrive as a Poisson process with rate λ>0\lambda > 0λ>0. Customer jjj has a base valuation VjV_jVj​ drawn independently from a distribution FFF with tail Fˉ(x)=1−F(x)\bar F(x) = 1 - F(x)Fˉ(x)=1−F(x), and values the product at Vje−αtV_j e^{-\alpha t}Vj​e−αt at time ttt, where α≥0\alpha \ge 0α≥0 is the decline factor. The paper reparametrizes it as ρ=e−αH\rho = e^{-\alpha H}ρ=e−αH, the fraction of the base valuation left at the end of the season.

In the numerical study, FFF is a Gamma law with mean μ\muμ and coefficient of variation ccc (standard deviation over mean): shape 1/c21/c^21/c2 and rate 1/(μc2)1/(\mu c^2)1/(μc2). The paper sets μ=1\mu = 1μ=1. For c=1c = 1c=1 this is the exponential law with mean one, Fˉ(x)=e−x\bar F(x) = e^{-x}Fˉ(x)=e−x for x≥0x \ge 0x≥0.

A contingent two-price policy posts the premium price p1p_1p1​ on [0,T)[0, T)[0,T), where 0<T≤10 < T \le 10<T≤1 is fixed, and a discount price p2≤p1p_2 \le p_1p2​≤p1​ from time TTT on. A myopic customer arriving at t<Tt < Tt<T buys at p1p_1p1​ if his valuation is at least p1p_1p1​; otherwise he waits and buys at TTT if his valuation is then at least p2p_2p2​. Customers arriving at or after TTT buy if their valuation is at least p2p_2p2​. The numbers of customers in these groups are Poisson with means

ΛI(p1)=λ∫0TFˉ(p1eαt) dt,ΛW(p1,p2)=λ∫0T[Fˉ(min⁡{p1eαt,p2eαT})−Fˉ(p1eαt)]dt,ΛL(p2)=λ∫THFˉ(p2eαt) dt.\Lambda_I(p_1) = \lambda\int_0^T \bar F(p_1 e^{\alpha t})\,dt, \quad \Lambda_W(p_1,p_2) = \lambda\int_0^T \big[\bar F(\min\{p_1e^{\alpha t}, p_2e^{\alpha T}\}) - \bar F(p_1e^{\alpha t})\big]dt, \quad \Lambda_L(p_2) = \lambda\int_T^H \bar F(p_2e^{\alpha t})\,dt .ΛI​(p1​)=λ∫0T​Fˉ(p1​eαt)dt,ΛW​(p1​,p2​)=λ∫0T​[Fˉ(min{p1​eαt,p2​eαT})−Fˉ(p1​eαt)]dt,ΛL​(p2​)=λ∫TH​Fˉ(p2​eαt)dt.

With unlimited inventory, the expected revenue of the policy is

RC/N(p1,p2)=p1ΛI(p1)+p2(ΛW(p1,p2)+ΛL(p2)),R_{C/N}(p_1, p_2) = p_1\Lambda_I(p_1) + p_2\big(\Lambda_W(p_1,p_2) + \Lambda_L(p_2)\big),RC/N​(p1​,p2​)=p1​ΛI​(p1​)+p2​(ΛW​(p1​,p2​)+ΛL​(p2​)),

and the expected revenue of a single price ppp is RF(p)=p λ∫0HFˉ(peαt) dtR_F(p) = p\,\lambda\int_0^H \bar F(p e^{\alpha t})\,dtRF​(p)=pλ∫0H​Fˉ(peαt)dt (Eq. (9) of the paper). The optimal values are πC/N∗=max⁡p2≤p1RC/N(p1,p2)\pi^*_{C/N} = \max_{p_2 \le p_1} R_{C/N}(p_1,p_2)πC/N∗​=maxp2​≤p1​​RC/N​(p1​,p2​) and πF∗=max⁡pRF(p)\pi^*_F = \max_p R_F(p)πF∗​=maxp​RF​(p).

Formalization targets

Goal: Proposition 3

Suppose c=1c = 1c=1, ρ=1\rho = 1ρ=1 and Q/λ→∞Q/\lambda \to \inftyQ/λ→∞ (unlimited inventory), with μ=1\mu = 1μ=1 and H=1H = 1H=1. Then

πC/N∗=(λe−1)⋅eT/e=πF∗⋅eT/e.\pi^*_{C/N} = (\lambda e^{-1})\cdot e^{T/e} = \pi^*_F \cdot e^{T/e}.πC/N∗​=(λe−1)⋅eT/e=πF∗​⋅eT/e.

Both maxima are attained. The goal states the two optimal values; it does not fix the optimal prices.

Milestones from the paper's proof

  1. The reduced problem: for 0≤p2≤p10 \le p_2 \le p_10≤p2​≤p1​, RC/N(p1,p2)=p2⋅λe−p2+(p1−p2)⋅λTe−p1R_{C/N}(p_1,p_2) = p_2\cdot\lambda e^{-p_2} + (p_1-p_2)\cdot\lambda T e^{-p_1}RC/N​(p1​,p2​)=p2​⋅λe−p2​+(p1​−p2​)⋅λTe−p1​.
  2. Its solution: over p2≤p1p_2 \le p_1p2​≤p1​ the maximum is λe−1+T/e\lambda e^{-1+T/e}λe−1+T/e, attained exactly at p1∗=2−T/e≥1p_1^* = 2 - T/e \ge 1p1∗​=2−T/e≥1, p2∗=p1∗−1≤1p_2^* = p_1^* - 1 \le 1p2∗​=p1∗​−1≤1.
  3. The fixed-price optimum (a supporting item of the goal, stated in the proof on pp. 358–359): p∗=μ=1p^* = \mu = 1p∗=μ=1 is the unique optimal single price and πF∗=λe−1\pi^*_F = \lambda e^{-1}πF∗​=λe−1.

Significance

Proposition 3 gives the relative benefit of contingent pricing over a single price, eT/e−1e^{T/e} - 1eT/e−1, as a function of the discount time alone. It increases in TTT and is largest at T=1T = 1T=1, where it equals e1/e−1≈44.46%e^{1/e} - 1 \approx 44.46\%e1/e−1≈44.46%. This is the paper's analytic anchor for its numerical findings: segmentation is most valuable when valuations are heterogeneous and customers are carried to the discount at little cost, and a late discount exposes more customers to the premium price. Under strategic customers the same quantity serves as an upper bound on the benefit of segmentation (§6.1 of the paper).

The result is proved in the paper, in a short appendix argument that states the reduced problem and its solution without the calculus. No machine-checked version exists. Formalizing it produces a reusable Lean encoding of the paper's segment rates ΛI,ΛW,ΛL\Lambda_I, \Lambda_W, \Lambda_LΛI​,ΛW​,ΛL​ as integrals of a valuation tail, a Gamma valuation law through Mathlib's gammaMeasure, and a complete verification that the integral model reduces to the two-variable problem and that the stated prices are its unique maximizer.

Difficulty

The obvious route is to write the revenue in closed form and set the gradient to zero. Two steps of that route are not automatic. First, the reduction requires evaluating the three integrals with the piecewise tail of the exponential law, including the min⁡\minmin inside ΛW\Lambda_WΛW​, and the reduced formula is valid only for nonnegative prices; negative prices must be handled separately in the model itself, where the tail equals one. Second, the reduced objective p2λe−p2+(p1−p2)λTe−p1p_2\lambda e^{-p_2} + (p_1-p_2)\lambda T e^{-p_1}p2​λe−p2​+(p1​−p2​)λTe−p1​ is not concave on the region p2≤p1p_2 \le p_1p2​≤p1​, so a stationary point is not automatically a global maximizer, and the boundary p2=p1p_2 = p_1p2​=p1​ and unbounded directions have to be ruled out. Uniqueness of the maximizer, which the paper asserts, fails at T=0T = 0T=0 and needs T>0T > 0T>0.

Formalization scope

All declarations sit in the namespace SeasonalPricing.MyopicExp. Time, prices and rates are real numbers. The season is [0,1][0, 1][0,1] with 0<T≤10 < T \le 10<T≤1 and λ>0\lambda > 0λ>0. Integrals are interval integrals. The valuation tail is gammaValuationTail μ c x = 1 - cdf (gammaMeasure (1/c^2) (1/(μ c^2))) x, used at μ=c=1\mu = c = 1μ=c=1. The hypothesis ρ=1\rho = 1ρ=1 is decayRatio α 1 = 1 with α≥0\alpha \ge 0α≥0.

Readings of the paper's informal words:

  • "Q/λ→∞Q/\lambda \to \inftyQ/λ→∞" is read as unlimited inventory: the truncated Poisson mean N(q,Λ)N(q,\Lambda)N(q,Λ) of §4.2 is replaced by Λ\LambdaΛ and stock-outs never occur. This is what the proof computes, what p. 348 writes as Q=∞Q = \inftyQ=∞, and what §7.1 calls inventory that is "practically unlimited". A limit of finite-inventory optimal revenues is not stated.
  • "max" is an attained maximum (IsGreatest), not a supremum.
  • The optimum is taken over all real prices with p2≤p1p_2 \le p_1p2​≤p1​, as printed; the paper never restricts signs, and negative prices are never optimal in the model.
  • The seller's discount at TTT is a best response to p1p_1p1​ in the paper (R(q∣p1)R(q \mid p_1)R(q∣p1​), p. 349). With unlimited inventory it does not depend on the realized sales, and the nested maximum equals the joint maximum over (p1,p2)(p_1, p_2)(p1​,p2​), which is what the goal states.
  • "The solution … is" (milestone 2) and "the optimal single price is given by p∗=μ=1p^* = \mu = 1p∗=μ=1" (the fixed-price item) are read as unique maximizers.

The Gamma density printed on p. 349 has the exponent 1/(sc2−1)1/(sc^2-1)1/(sc2−1), a misprint for 1/c2−11/c^2 - 11/c2−1; at c=1c = 1c=1 the exponent is 000 either way.

A trivializing formalization would state the goal on the reduced two-variable function, dropping the model: the goal here is about RC/NR_{C/N}RC/N​ built from ΛI,ΛW,ΛL\Lambda_I, \Lambda_W, \Lambda_LΛI​,ΛW​,ΛL​ and the Gamma tail, and about RFR_FRF​ built from Eq. (9). The platform's BuyingToBundle.monopolyRevenue (definition monopoly_pricing) is a related object, sup⁡pp ν([p,∞))\sup_p p\,\nu([p,\infty))supp​pν([p,∞)); with ρ=1\rho = 1ρ=1 and H=1H = 1H=1, πF∗\pi^*_FπF∗​ equals λ\lambdaλ times it for the exponential law, but it is a supremum without arrivals or time and is not reused.

Contributions welcome: closed forms of the segment rates for the exponential tail, a general lemma that negative prices are dominated, and the two-variable maximization.

Selected references

  • Y. Aviv and A. Pazgal, Optimal Pricing of Seasonal Products in the Presence of Forward-Looking Consumers, Manufacturing & Service Operations Management 10(3):339–359, 2008. https://doi.org/10.1287/msom.1070.0183
  • D. Besanko and W. L. Winston, Optimal Price Skimming by a Monopolist Facing Rational Consumers, Management Science 36(5):555–567, 1990. https://doi.org/10.1287/mnsc.36.5.555
  • G. Gallego and G. van Ryzin, Optimal Dynamic Pricing of Inventories with Stochastic Demand over Finite Horizons, Management Science 40(8):999–1020, 1994. https://doi.org/10.1287/mnsc.40.8.999
6 thms2 active usersReviewed
🏆Completed
Operations ResearchPartial Differential EquationsProbability+1·Captain: mikedeng1

Revenue Management of a Make-to-Stock Queue: Exponential Stationary Density under Normal Reflection (Proposition 2)Research Paper

Motivation

A make-to-stock manufacturer who also sells on a spot market must decide, at every moment, whether to keep producing and whether to accept or reject incoming orders at the prevailing price. Caldentey and Wein (Revenue Management of a Make-to-Stock Queue, Operations Research 54(5), 2006) study this problem in heavy traffic. The limit is a two-dimensional singular control problem for a diffusion: the inventory level and the logarithm of the price move jointly as a correlated Brownian motion, and the controls push the inventory only when it reaches one of two free boundaries. The optimal boundaries are characterized by an elliptic free-boundary problem that the authors could not solve in closed form.

The paper's way forward is an approximation: change the direction of reflection on the boundary so that the stationary distribution of the controlled process becomes an explicit exponential. Proposition 2 states that exponential form, and it turns the free-boundary problem into a calculus-of-variations problem for the two boundary curves. Explicit stationary densities of reflected diffusions in two dimensions are rare; the classical condition for an exponential stationary density of a reflected Brownian motion, and the characterization of the stationary law by a basic adjoint relation, are due to Harrison and Williams, Multidimensional reflected Brownian motions having exponential stationary distributions, Annals of Probability 15, 1987, the reference the paper cites. This mission formalizes the analytic core of Proposition 2: the exponential density satisfies that relation for the reflection field the proposition singles out.

Setting

Points of the plane are (x,y)(x,y)(x,y), with xxx the inventory level and yyy the logarithm of the price. The limiting process (X,Y)(\mathcal X,\mathcal Y)(X,Y) has drift (θ,0)(\theta,0)(θ,0) and covariance matrix

Σ=(σ2σδϱσδϱδ2),σ>0, δ>0, −1<ϱ<1,\Sigma=\begin{pmatrix}\sigma^2&\sigma\delta\varrho\\ \sigma\delta\varrho&\delta^2\end{pmatrix},\qquad \sigma>0,\ \delta>0,\ -1<\varrho<1,Σ=(σ2σδϱ​σδϱδ2​),σ>0, δ>0, −1<ϱ<1,

so its generator is

Γ=θ∂∂x+σ22∂2∂x2+σδϱ∂2∂x ∂y+δ22∂2∂y2.\Gamma=\theta\frac{\partial}{\partial x}+\frac{\sigma^2}{2}\frac{\partial^2}{\partial x^2}+\sigma\delta\varrho\frac{\partial^2}{\partial x\,\partial y}+\frac{\delta^2}{2}\frac{\partial^2}{\partial y^2}.Γ=θ∂x∂​+2σ2​∂x2∂2​+σδϱ∂x∂y∂2​+2δ2​∂y2∂2​.

Two curves bound the region where the process lives: the rejection boundary x=η(y)x=\eta(y)x=η(y) (below it, orders are rejected) and the idleness boundary x=ξ(y)x=\xi(y)x=ξ(y) (above it, production stops). For ymin⁡<ymax⁡y_{\min}<y_{\max}ymin​<ymax​ the region is

Ω={(x,y): ymin⁡<y<ymax⁡, η(y)<x<ξ(y)},\Omega=\{(x,y):\ y_{\min}<y<y_{\max},\ \eta(y)<x<\xi(y)\},Ω={(x,y): ymin​<y<ymax​, η(y)<x<ξ(y)},

and its boundary splits into four pieces: x=η(y)x=\eta(y)x=η(y), x=ξ(y)x=\xi(y)x=ξ(y), y=ymin⁡y=y_{\min}y=ymin​, y=ymax⁡y=y_{\max}y=ymax​. Write n⃗\vec nn for the inward unit normal on ∂Ω\partial\Omega∂Ω and dldldl for arc length. A reflection field v⃗\vec vv on ∂Ω\partial\Omega∂Ω gives the direction in which the process is pushed back into Ω\OmegaΩ. The basic adjoint relation (BAR) of the paper, equation (43), is

∫ΩΓf πΩ ds+12∫∂Ωv⃗⋅∇f πΩ dl=0for all test functions f,\int_\Omega \Gamma f\,\pi_\Omega\,ds+\frac12\int_{\partial\Omega}\vec v\cdot\nabla f\,\pi_\Omega\,dl=0\quad\text{for all test functions } f,∫Ω​ΓfπΩ​ds+21​∫∂Ω​v⋅∇fπΩ​dl=0for all test functions f,

and the paper cites Harrison and Williams for the fact that the stationary distribution πΩ\pi_\OmegaπΩ​ of the reflected process satisfies it. Proposition 2 introduces the eigen-decomposition Σ=V′EV\Sigma=V'EVΣ=V′EV (VVV a rotation whose rows are eigenvectors, EEE diagonal), the whitening map T=E−1/2VT=E^{-1/2}VT=E−1/2V and Ω∗=T(Ω)\Omega^*=T(\Omega)Ω∗=T(Ω), and assumes that Tv⃗T\vec vTv is normal to ∂Ω∗\partial\Omega^*∂Ω∗. The exponents are

mx=2θσ2(1−ϱ2),my=−2ϱθσδ(1−ϱ2).(47)m_x=\frac{2\theta}{\sigma^2(1-\varrho^2)},\qquad m_y=\frac{-2\varrho\theta}{\sigma\delta(1-\varrho^2)}.\tag{47}mx​=σ2(1−ϱ2)2θ​,my​=σδ(1−ϱ2)−2ϱθ​.(47)

Formalization targets

Goal: the exponential density satisfies the BAR under conormal reflection

For η,ξ\eta,\xiη,ξ continuously differentiable with η<ξ\eta<\xiη<ξ on [ymin⁡,ymax⁡][y_{\min},y_{\max}][ymin​,ymax​], π(x,y)=emxx+myy\pi(x,y)=e^{m_xx+m_yy}π(x,y)=emx​x+my​y, and every C2C^2C2 function fff on R2\mathbb R^2R2,

∫ΩΓf  π ds+12∫∂Ω(Σn⃗)⋅∇f  π dl=0.\int_\Omega \Gamma f\;\pi\,ds+\frac12\int_{\partial\Omega}(\Sigma\vec n)\cdot\nabla f\;\pi\,dl=0 .∫Ω​Γfπds+21​∫∂Ω​(Σn)⋅∇fπdl=0.

The boundary integral is written out on the four pieces, with n⃗ dl\vec n\,dlndl equal to (1,−η′(y)) dy(1,-\eta'(y))\,dy(1,−η′(y))dy, (−1,ξ′(y)) dy(-1,\xi'(y))\,dy(−1,ξ′(y))dy, (0,1) dx(0,1)\,dx(0,1)dx and (0,−1) dx(0,-1)\,dx(0,−1)dx respectively. The normalizing constant is left out because the relation is linear in π\piπ.

Milestones

  1. The interior equation: Γ∗π=−θπx+σ22πxx+σδϱ πxy+δ22πyy=0\Gamma^*\pi=-\theta\pi_x+\frac{\sigma^2}{2}\pi_{xx}+\sigma\delta\varrho\,\pi_{xy}+\frac{\delta^2}{2}\pi_{yy}=0Γ∗π=−θπx​+2σ2​πxx​+σδϱπxy​+2δ2​πyy​=0 everywhere.
  2. The zero-flux identity: 12Σ∇π=(θ,0) π\frac12\Sigma\nabla\pi=(\theta,0)\,\pi21​Σ∇π=(θ,0)π everywhere.
  3. The meaning of the hypothesis: with T=E−1/2VT=E^{-1/2}VT=E−1/2V, (Tv)⋅(Tw)=v⋅Σ−1w(Tv)\cdot(Tw)=v\cdot\Sigma^{-1}w(Tv)⋅(Tw)=v⋅Σ−1w, and for n≠0n\neq0n=0, TvTvTv is orthogonal to TTT of every vector orthogonal to nnn exactly when vvv is a multiple of Σn\Sigma nΣn.
  4. The normalizing constant: π\piπ is integrable on Ω\OmegaΩ and a unique KΩ>0K_\Omega>0KΩ​>0 makes KΩπK_\Omega\piKΩ​π integrate to one.

Significance

For the operations model, Proposition 2 is what makes the problem computable. Once the stationary density is explicit, the long-run average cost of any pair of boundary curves is an explicit integral, and optimizing over (η,ξ)(\eta,\xi)(η,ξ) becomes a variational problem with Euler–Lagrange equations; the paper's proposed policy and its numerical comparisons all rest on it.

For formalization, the mission produces a machine-checked version of a statement whose proof the paper does not contain (it is in an online companion) and whose hypothesis is stated only in words. The formal statements fix exactly which reflection field makes the claim true, which the prose leaves ambiguous. None of the statements has, to our knowledge, a machine-checked proof anywhere; the result itself is classical in spirit (an integration by parts on a planar region), but no divergence theorem on a region between two graphs with an anisotropic operator is currently available as a ready-made statement.

Difficulty

The interior equation and the zero-flux identity are finite computations with the exponential. The difficulty is the goal: it is an integration-by-parts identity on a curved planar region with an anisotropic second-order operator. The obvious first step, "apply Green's identity", presupposes a divergence theorem on a region bounded by two graphs x=η(y)x=\eta(y)x=η(y), x=ξ(y)x=\xi(y)x=ξ(y) and two horizontal segments, with the boundary integral written in the parametrization of each piece and the orientation of every normal tracked. Mathlib has the divergence theorem on rectangular boxes, not on such regions, and the moving limits η(y)\eta(y)η(y), ξ(y)\xi(y)ξ(y) are exactly where the terms in η′\eta'η′ and ξ′\xi'ξ′ of the boundary integral come from.

The second trap is the reflection field. The page describes the modification as substituting the inward unit normal n⃗\vec nn for v⃗\vec vv; with v⃗=n⃗\vec v=\vec nv=n the identity is false as soon as Σ\SigmaΣ is not a multiple of the identity (on a random instance the residual is of order one). Only the conormal field Σn⃗\Sigma\vec nΣn, which is what the hypothesis of Proposition 2 selects, gives a true statement.

Formalization scope

Everything lives in the namespace MakeToStockRM.ExpDensity. The plane is ℝ × ℝ with the inventory first; partial derivatives are Fréchet derivatives applied to (1, 0) and (0, 1), and the mixed partial is ∂x(∂yf)\partial_x(\partial_y f)∂x​(∂y​f). Parameters satisfy σ>0\sigma>0σ>0, δ>0\delta>0δ>0, ∣ϱ∣<1|\varrho|<1∣ϱ∣<1; θ\thetaθ is any real number, and θ=0\theta=0θ=0 (then π≡1\pi\equiv1π≡1) is allowed.

This is the analytic, pinned-down content of Proposition 2. The identification "the BAR characterizes the stationary law of the reflected diffusion" (Harrison–Williams 1987) is out of scope: Mathlib has no reflected Brownian motion. Relative to the page, the formalization commits to the following:

  • The reflection field is v⃗=Σn⃗\vec v=\Sigma\vec nv=Σn with n⃗\vec nn the inward unit normal and dldldl arc length. The hypothesis "Tv⃗T\vec vTv is normal to ∂Ω∗\partial\Omega^*∂Ω∗" fixes only the direction of v⃗\vec vv (milestone 3); the length Σn⃗\Sigma\vec nΣn is the one for which the BAR holds. The page's phrase "substituting the inward unit normal n⃗\vec nn for v⃗\vec vv" is inconsistent with the proposition's own hypothesis and is not followed.
  • The boundary curves are C1C^1C1 on R\mathbb RR with η<ξ\eta<\xiη<ξ on [ymin⁡,ymax⁡][y_{\min},y_{\max}][ymin​,ymax​], and ymin⁡<ymax⁡y_{\min}<y_{\max}ymin​<ymax​, so Ω\OmegaΩ is a nonempty bounded region; the paper assumes this implicitly.
  • Test functions are all C2C^2C2 functions on R2\mathbb R^2R2, which are bounded with bounded derivatives on the closure of Ω\OmegaΩ (the paper's "twice continuous and bounded").
  • The constant KΩK_\OmegaKΩ​ is dropped from the goal and treated in milestone 4.

The goal quantifies over every C2C^2C2 test function; restricting to functions supported inside Ω\OmegaΩ would delete the boundary term and reduce the goal to milestone 1, and that trivialization is ruled out. The second half of Proposition 2 ("(45)–(46) is equivalent to (48)–(49)"), Proposition 1, the heavy-traffic limit, the HJB equation and the proposed policy are not formalized: their normalizations or proofs are only in the online companion.

A complete development needs a divergence theorem on regions between two C1C^1C1 graphs, which is reusable for any planar PDE statement on such regions. Contributions of that lemma, and of the four milestones, are welcome.

Selected references

  • R. Caldentey, L. M. Wein, Revenue Management of a Make-to-Stock Queue, Operations Research 54(5):859–875, 2006. https://doi.org/10.1287/opre.1060.0289
  • J. M. Harrison, R. J. Williams, Multidimensional reflected Brownian motions having exponential stationary distributions, Annals of Probability 15(1):115–137, 1987. https://doi.org/10.1214/aop/1176992259
  • F. John, Partial Differential Equations, 4th ed., Springer, 1982. https://doi.org/10.1007/978-1-4684-9333-7
10 thms2 active usersReviewed
🏆Completed
Dynamic ProgrammingOperations ResearchOptimization+1·Captain: mikedeng1

An Overview of Pricing Models for Revenue Management: Performance Guarantee of the Deterministic Price Heuristic in Periodic-Review PricingResearch Paper

Motivation

Dynamic pricing under limited inventory is a core problem of revenue management: a seller holds C0C_0C0​ units of a perishable product (airline seats, hotel rooms, seasonal goods) and must choose prices over a finite selling horizon while demand is random and responds to price. The exactly optimal policy solves a stochastic dynamic program whose state is the remaining inventory, and it changes the price after every sale. In practice sellers often prefer a simpler rule, fixed in advance: solve the deterministic version of the problem, in which random demand is replaced by its mean, and charge the resulting prices whatever happens.

Gallego and van Ryzin (Management Science 1994) showed that in continuous time with Poisson demand this fixed-price heuristic is asymptotically optimal, and bounded its relative loss by the coefficient of variation of demand. Bitran and Caldentey's survey (MSOM 2003, §3.2.1) extends this bound to a discrete-time, periodic-review model with general demand distributions, as Proposition 8, proved in the paper's Appendix. The proof combines three ingredients: a Lagrangian duality argument showing that the deterministic problem is an upper bound on the optimal expected revenue, a sample-path comparison of lost sales, and a distribution-free moment bound of Gallego (1992) on E[(X−C)+]E[(X-C)^+]E[(X−C)+].

Setting

A single product is sold over N≥1N \ge 1N≥1 periods n=1,…,Nn = 1, \dots, Nn=1,…,N starting from inventory C0C_0C0​. In period nnn the seller charges a price p≥0p \ge 0p≥0, and the demand Dn(p)D_n(p)Dn​(p) is a nonnegative random variable with finite mean E[Dn(p)]E[D_n(p)]E[Dn​(p)], whose law may depend on nnn and ppp arbitrarily. With inventory CCC the seller sells min⁡{D,C}\min\{D, C\}min{D,C} units; unmet demand is lost.

The optimal expected revenue V1(C0)V_1(C_0)V1​(C0​) is defined by the Bellman recursion VN+1≡0V_{N+1} \equiv 0VN+1​≡0,

Vn(C)=sup⁡p≥0E[pmin⁡{Dn(p),C}+Vn+1(C−min⁡{Dn(p),C})].V_n(C) = \sup_{p \ge 0} E\big[p\min\{D_n(p), C\} + V_{n+1}\big(C - \min\{D_n(p), C\}\big)\big].Vn​(C)=p≥0sup​E[pmin{Dn​(p),C}+Vn+1​(C−min{Dn​(p),C})].

The deterministic problem (32)–(33) is

V1det⁡(C0)=sup⁡p∈[0,∞)N∑n=1NpnE[Dn(pn)]subject to∑n=1NE[Dn(pn)]≤C0,V_1^{\det}(C_0) = \sup_{p \in [0,\infty)^N} \sum_{n=1}^N p_n E[D_n(p_n)] \quad \text{subject to} \quad \sum_{n=1}^N E[D_n(p_n)] \le C_0,V1det​(C0​)=p∈[0,∞)Nsup​n=1∑N​pn​E[Dn​(pn​)]subject ton=1∑N​E[Dn​(pn​)]≤C0​,

and pdet⁡p^{\det}pdet denotes an optimal solution. The deterministic-price heuristic charges pndet⁡p^{\det}_npndet​ in period nnn regardless of sales. Its period demands Dn(pndet⁡)D_n(p_n^{\det})Dn​(pndet​) are independent, Dndet⁡=∑i=1nDi(pidet⁡)\mathscr{D}_n^{\det} = \sum_{i=1}^n D_i(p_i^{\det})Dndet​=∑i=1n​Di​(pidet​) is the cumulative demand, and its expected revenue is

V1(pdet⁡,C0)=∑n=1Npndet⁡E[Dn(pndet⁡)−(Dn(pndet⁡)−(C0−Dn−1det⁡)+)+].V_1(p^{\det}, C_0) = \sum_{n=1}^N p_n^{\det} E\Big[D_n(p_n^{\det}) - \big(D_n(p_n^{\det}) - (C_0 - \mathscr{D}_{n-1}^{\det})^+\big)^+\Big].V1​(pdet,C0​)=n=1∑N​pndet​E[Dn​(pndet​)−(Dn​(pndet​)−(C0​−Dn−1det​)+)+].

With σn2=Var⁡(Dndet⁡)\sigma_n^2 = \operatorname{Var}(\mathscr{D}_n^{\det})σn2​=Var(Dndet​), eq. (34) defines

ηndet⁡(C0)=σn2+(C0−E[Dndet⁡])2−(C0−E[Dndet⁡])2,\eta_n^{\det}(C_0) = \frac{\sqrt{\sigma_n^2 + (C_0 - E[\mathscr{D}_n^{\det}])^2} - (C_0 - E[\mathscr{D}_n^{\det}])}{2},ηndet​(C0​)=2σn2​+(C0​−E[Dndet​])2​−(C0​−E[Dndet​])​,

and ν(C0)=σN/E[DNdet⁡]\nu(C_0) = \sigma_N / E[\mathscr{D}_N^{\det}]ν(C0​)=σN​/E[DNdet​] is the coefficient of variation of the total demand.

Formalization targets

Goal: Proposition 8, eq. (35)

Assume that p↦p E[Dn(p)]p \mapsto p\,E[D_n(p)]p↦pE[Dn​(p)] is concave and p↦E[Dn(p)]p \mapsto E[D_n(p)]p↦E[Dn​(p)] is convex on [0,∞)[0,\infty)[0,∞) for each nnn, that some price p∞≥0p^\infty \ge 0p∞≥0 has ∑nE[Dn(p∞)]<C0\sum_n E[D_n(p^\infty)] < C_0∑n​E[Dn​(p∞)]<C0​, and that pdet⁡p^{\det}pdet is optimal for (32)–(33). Then

1≥V1(pdet⁡,C0)V1(C0)≥1V1det⁡(C0)∑n=1Npndet⁡E[Dn(pndet⁡)](1−ηndet⁡(C0)E[Dn(pndet⁡)])≥1−max⁡nηndet⁡(C0)E[Dn(pndet⁡)].1 \ge \frac{V_1(p^{\det}, C_0)}{V_1(C_0)} \ge \frac{1}{V_1^{\det}(C_0)} \sum_{n=1}^N p_n^{\det} E[D_n(p_n^{\det})]\left(1 - \frac{\eta_n^{\det}(C_0)}{E[D_n(p_n^{\det})]}\right) \ge 1 - \max_n \frac{\eta_n^{\det}(C_0)}{E[D_n(p_n^{\det})]}.1≥V1​(C0​)V1​(pdet,C0​)​≥V1det​(C0​)1​n=1∑N​pndet​E[Dn​(pndet​)](1−E[Dn​(pndet​)]ηndet​(C0​)​)≥1−nmax​E[Dn​(pndet​)]ηndet​(C0​)​.

The constants are the paper's, and the goal is the full chain.

Milestones

  1. Gallego's bound (Appendix, (*)): for square-integrable XXX, E[(X−C)+]≤12(Var⁡X+(C−EX)2−(C−EX))E[(X - C)^+] \le \frac12\big(\sqrt{\operatorname{Var}X + (C - EX)^2} - (C - EX)\big)E[(X−C)+]≤21​(VarX+(C−EX)2​−(C−EX)), which is at most 12Var⁡X\frac12\sqrt{\operatorname{Var} X}21​VarX​ when EX≤CEX \le CEX≤C.
  2. Eq. (24): in one period, E[pmin⁡{D(p),C}]≤pmin⁡{E[D(p)],C}E[p\min\{D(p), C\}] \le p\min\{E[D(p)], C\}E[pmin{D(p),C}]≤pmin{E[D(p)],C}, so V(C)≤Vdet⁡(C)V(C) \le V^{\det}(C)V(C)≤Vdet(C).
  3. Proposition 6, eq. (26): in one period, V(C,pdet⁡)/V(C)≥1−νdet⁡/2V(C, p^{\det})/V(C) \ge 1 - \nu^{\det}/2V(C,pdet)/V(C)≥1−νdet/2.
  4. V1det⁡V_1^{\det}V1det​ is concave in the capacity (first sentence of the Appendix proof).
  5. Proposition 8, first assertion: V1(C0)≤V1det⁡(C0)V_1(C_0) \le V_1^{\det}(C_0)V1​(C0​)≤V1det​(C0​).
  6. The Appendix display bounding V1(pdet⁡,C0)V_1(p^{\det}, C_0)V1​(pdet,C0​) below by ∑npndet⁡E[Dn](1−E[(Dndet⁡−C0)+]/E[Dn])\sum_n p_n^{\det}E[D_n](1 - E[(\mathscr{D}_n^{\det} - C_0)^+]/E[D_n])∑n​pndet​E[Dn​](1−E[(Dndet​−C0​)+]/E[Dn​]).
  7. Eq. (36), after the goal: if E[Dn(p)]=Tnλ(p)E[D_n(p)] = T_n\lambda(p)E[Dn​(p)]=Tn​λ(p), one constant price solves (32)–(33) and 1≥V1(pdet⁡,C0)/V1(C0)≥1−ν(C0)/21 \ge V_1(p^{\det},C_0)/V_1(C_0) \ge 1 - \nu(C_0)/21≥V1​(pdet,C0​)/V1​(C0​)≥1−ν(C0​)/2.

Significance

The result gives a guarantee for a pricing policy that needs no inventory tracking: its relative loss is controlled by the first two moments of cumulative demand, with no distributional assumption beyond finite variance. When demand grows while its coefficient of variation shrinks, as for sums of independent period demands, the guarantee tends to one, which is the discrete-time form of asymptotic optimality of fixed prices. The upper bound V1≤V1det⁡V_1 \le V_1^{\det}V1​≤V1det​ is used throughout revenue management as the benchmark for heuristics (fluid or deterministic LP bounds).

The paper's proof is complete in the Appendix, and the mathematics is not in question beyond minor typos. None of it is machine-checked. A formal development adds a checked Bellman model of periodic-review pricing with general demand laws, a checked fluid upper bound for it, and a checked form of Gallego's moment bound, all reusable for other revenue-management statements. A related platform theorem, RevenueManagement.deterministic_upper_bound (mission The Theory and Practice of Revenue Management III), states the deterministic bound for a Bernoulli-arrival model with at most one sale per period; the model here is different (arbitrary demand laws, continuous inventory), and neither statement implies the other.

Difficulty

The upper bound V1(C0)≤V1det⁡(C0)V_1(C_0) \le V_1^{\det}(C_0)V1​(C0​)≤V1det​(C0​) is the hard part. The natural induction replaces V2V_2V2​ by V2det⁡V_2^{\det}V2det​ and applies Jensen's inequality, but V2det⁡(C)V_2^{\det}(C)V2det​(C) is −∞-\infty−∞ at capacities where the remaining deterministic problem is infeasible, while V2(C)V_2(C)V2​(C) stays nonnegative. The induction therefore fails as stated when mean demand never vanishes. The paper's argument also passes through the Lagrangian dual of (32)–(33), and strong duality for that program on [0,∞)N[0,\infty)^N[0,∞)N under the Slater point p∞p^\inftyp∞ is not available in Mathlib in this form.

The lower bound is less deep but technical: it needs the product law of the period demands, linearity of expectation for truncated sums, and variances of partial sums. Gallego's inequality is elementary once the right quadratic bound on (x−C)+(x - C)^+(x−C)+ is found, but it is not in Mathlib.

Formalization scope

  • Periods are Fin N, 0-based. Demand laws are μ n p : Measure ℝ, probability measures carried by [0,∞)[0,\infty)[0,∞) with finite mean, required at every price (laws at negative prices never enter a statement).
  • The Bellman value is computed in [0,∞][0,\infty][0,∞] with the lower Lebesgue integral and ⨆ over p≥0p \ge 0p≥0, so no supremum takes a junk value; the paper's VnV_nVn​ is valueToGo M (N - n + 1). Ratios use its real value.
  • V1det⁡V_1^{\det}V1det​ is an EReal supremum over the feasible set (−∞-\infty−∞ if infeasible). In the goal pdet⁡p^{\det}pdet is an optimal solution, so V1det⁡(C0)V_1^{\det}(C_0)V1det​(C0​) is its objective value.
  • The heuristic's demands are independent: the joint law is the product measure.
  • Implicit hypotheses made explicit: "concave objective and convex feasible region" is read as concavity of p E[Dn(p)]p\,E[D_n(p)]pE[Dn​(p)] and convexity of E[Dn(p)]E[D_n(p)]E[Dn​(p)] on [0,∞)[0,\infty)[0,∞), the latter making (33) convex for every capacity; existence of the optimal deterministic solution; finite variance and positive mean of each Dn(pndet⁡)D_n(p_n^{\det})Dn​(pndet​); and V1det⁡(C0)>0V_1^{\det}(C_0) > 0V1det​(C0​)>0, without which the ratios are 0/00/00/0.
  • The typo Dndet⁡:=∑i=1nDn(pdet⁡)\mathscr{D}_n^{\det} := \sum_{i=1}^n D_n(p^{\det})Dndet​:=∑i=1n​Dn​(pdet) on p. 221 is read as ∑i=1nDi(pidet⁡)\sum_{i=1}^n D_i(p_i^{\det})∑i=1n​Di​(pidet​).

Trivializing formalizations ruled out: V1V_1V1​ is defined by the Bellman recursion, not as a supremum over an unspecified policy class or as a variable constrained by hypotheses; ηndet⁡\eta_n^{\det}ηndet​ is the expression (34), not a hypothesis-supplied bound on expected overflow; the goal does not assume strong duality or a Lagrange multiplier, since that is the proof's key step; independence is built into the joint law, not assumed as an inequality; variances are taken only under square integrability, since Mathlib's variance is 000 for infinite variance.

Needed infrastructure: measurability of the Bellman integrand (monotonicity of VnV_nVn​ in inventory), Jensen's inequality for min⁡{⋅,C}\min\{\cdot, C\}min{⋅,C} (ConcaveOn.le_map_integral), Lagrangian strong duality for a separable concave program with one convex constraint, marginals and variances of sums under Measure.pi, and Gallego's bound. The duality and moment results are reusable beyond this mission. Proofs of any milestone, and alternative proofs of the upper bound, are welcome.

Selected references

  • G. R. Bitran, R. Caldentey, An Overview of Pricing Models for Revenue Management, Manufacturing & Service Operations Management 5(3):203–230, 2003. https://doi.org/10.1287/msom.5.3.203.16061
  • G. Gallego, G. van Ryzin, Optimal Dynamic Pricing of Inventories with Stochastic Demand over Finite Horizons, Management Science 40(8):999–1020, 1994. https://doi.org/10.1287/mnsc.40.8.999
  • G. Gallego, A Minmax Distribution Free Procedure for the (Q, R) Inventory Model, Operations Research Letters 11(1):55–60, 1992 (as cited in Bitran and Caldentey 2003).
  • M. S. Bazaraa, H. D. Sherali, C. M. Shetty, Nonlinear Programming: Theory and Algorithms, 2nd ed., Wiley, 1993.
9 thms3 active usersReviewed
🏆Completed
Operations ResearchOptimizationProbability·Captain: mikedeng1

Optimal Pricing of Seasonal Products in the Presence of Forward-Looking Consumers 4: Optimal Prices and Discount Time with Myopic Customers and Identical Declining ValuationsResearch Paper

Motivation

Retailers of fashion and seasonal goods sell a fixed stock over a short season and routinely cut the price part-way through it. The markdown trades off two effects: a late discount keeps early, high-valuation customers paying the full price, while an early discount reaches customers whose interest in the product fades as the season goes on. Aviv and Pazgal (MSOM 2008) build a two-price model of this trade-off with Poisson arrivals and valuations that decline exponentially over the season, and compare sellers facing myopic customers, who buy as soon as the current price is acceptable, with sellers facing strategic customers, who may wait for the discount.

This mission formalizes the benchmark of that comparison in which the problem can be solved in closed form: myopic customers who all share the same base valuation, so that the only source of price discrimination is the decline of valuations over time. Proposition 4 of the paper identifies the optimal premium price, discount price and discount time, and the paper's Proposition 5 and Example 1 then measure how much strategic behaviour costs the seller against it.

Setting

The season is [0,1][0, 1][0,1]. Customers arrive as a Poisson process with rate λ>0\lambda > 0λ>0, so λ\lambdaλ is the expected number of arrivals in the season. Every customer has base valuation 111, and a customer's valuation at time ttt is ρt\rho^tρt for a fixed decline parameter 0<ρ<10 < \rho < 10<ρ<1 (equivalently e−αte^{-\alpha t}e−αt with α=−ln⁡ρ\alpha = -\ln\rhoα=−lnρ; ρ\rhoρ is the fraction of the valuation left at the end of the season). In the paper's notation this is the case c=0c = 0c=0, μ=1\mu = 1μ=1, H=1H = 1H=1 of a family of Gamma-distributed base valuations with mean μ\muμ and coefficient of variation ccc; the tail of the base valuation is Fˉ(x)=1\bar F(x) = 1Fˉ(x)=1 for x≤1x \le 1x≤1 and 000 otherwise.

The seller posts a premium price p1p_1p1​ on [0,T)[0, T)[0,T) and a discount price p2≤p1p_2 \le p_1p2​≤p1​ from the discount time T∈[0,1]T \in [0, 1]T∈[0,1] on, and has unlimited inventory. A myopic customer arriving at t<Tt < Tt<T buys at once iff ρt≥p1\rho^t \ge p_1ρt≥p1​; otherwise the customer waits and buys at TTT iff ρT≥p2\rho^T \ge p_2ρT≥p2​. A customer arriving after TTT buys iff the current valuation is at least p2p_2p2​. The expected numbers of buyers in the three groups are the segment rates ΛI(p1)=λ∫0TFˉ(p1eαt) dt\Lambda_I(p_1) = \lambda\int_0^T \bar F(p_1e^{\alpha t})\,dtΛI​(p1​)=λ∫0T​Fˉ(p1​eαt)dt, ΛW(p1,p2)=λ∫0T[Fˉ(min⁡{p1eαt,p2eαT})−Fˉ(p1eαt)] dt\Lambda_W(p_1, p_2) = \lambda\int_0^T[\bar F(\min\{p_1e^{\alpha t}, p_2e^{\alpha T}\}) - \bar F(p_1 e^{\alpha t})]\,dtΛW​(p1​,p2​)=λ∫0T​[Fˉ(min{p1​eαt,p2​eαT})−Fˉ(p1​eαt)]dt and ΛL(p2)=λ∫T1Fˉ(p2eαt) dt\Lambda_L(p_2) = \lambda\int_T^1 \bar F(p_2 e^{\alpha t})\,dtΛL​(p2​)=λ∫T1​Fˉ(p2​eαt)dt, and the expected revenue is

Rρ(p1,p2;T)=p1 ΛI(p1)+p2 (ΛW(p1,p2)+ΛL(p2)).R_\rho(p_1, p_2; T) = p_1\,\Lambda_I(p_1) + p_2\,\big(\Lambda_W(p_1, p_2) + \Lambda_L(p_2)\big).Rρ​(p1​,p2​;T)=p1​ΛI​(p1​)+p2​(ΛW​(p1​,p2​)+ΛL​(p2​)).

For a price p∈[ρ,1]p \in [\rho, 1]p∈[ρ,1] let τ(p)=ln⁡p/ln⁡ρ\tau(p) = \ln p/\ln\rhoτ(p)=lnp/lnρ, the time at which the valuation has fallen to ppp, and write τ1=τ(p1)\tau_1 = \tau(p_1)τ1​=τ(p1​), τ2=τ(p2)\tau_2 = \tau(p_2)τ2​=τ(p2​). The reduced objective is

G(p1,p2)=(p1−p2) ln⁡p1ln⁡ρ+p2 ln⁡p2ln⁡ρ,ρ≤p2≤p1≤1.G(p_1, p_2) = (p_1 - p_2)\,\frac{\ln p_1}{\ln\rho} + p_2\,\frac{\ln p_2}{\ln\rho}, \qquad \rho \le p_2 \le p_1 \le 1 .G(p1​,p2​)=(p1​−p2​)lnρlnp1​​+p2​lnρlnp2​​,ρ≤p2​≤p1​≤1.

Formalization targets

Goal: Proposition 4 (p. 351)

πC/N∗=λ⋅max⁡ρ≤p2≤p1≤1G(p1,p2)=max⁡0<p2≤p1, 0≤T≤1Rρ(p1,p2;T),\pi^*_{C/N} = \lambda\cdot\max_{\rho \le p_2 \le p_1 \le 1} G(p_1, p_2) = \max_{0 < p_2 \le p_1,\ 0 \le T \le 1} R_\rho(p_1, p_2; T),πC/N∗​=λ⋅ρ≤p2​≤p1​≤1max​G(p1​,p2​)=0<p2​≤p1​, 0≤T≤1max​Rρ​(p1​,p2​;T),

every maximizer (p1∗,p2∗)(p_1^*, p_2^*)(p1∗​,p2∗​) of GGG together with every TTT with p2∗≤ρT≤p1∗p_2^* \le \rho^T \le p_1^*p2∗​≤ρT≤p1∗​ attains πC/N∗\pi^*_{C/N}πC/N∗​, and, if ρ≤e−2+e−1\rho \le e^{-2+e^{-1}}ρ≤e−2+e−1, the maximizer is unique,

p1∗=e−1+e−1,p2∗=p1∗/e,πC/N∗=−λ e−1+e−1ln⁡ρ,p_1^* = e^{-1+e^{-1}}, \qquad p_2^* = p_1^*/e, \qquad \pi^*_{C/N} = -\frac{\lambda\, e^{-1+e^{-1}}}{\ln\rho},p1∗​=e−1+e−1,p2∗​=p1∗​/e,πC/N∗​=−lnρλe−1+e−1​,

and every TTT with ρT∈[e−2+e−1,e−1+e−1]\rho^T \in [e^{-2+e^{-1}}, e^{-1+e^{-1}}]ρT∈[e−2+e−1,e−1+e−1] is optimal.

Milestones (Proof of Proposition 4, p. 359)

For ρ≤p2≤p1≤1\rho \le p_2 \le p_1 \le 1ρ≤p2​≤p1​≤1:

  1. Rρ(p1,p2;T)≤Rρ(p1,p2;τ1)R_\rho(p_1, p_2; T) \le R_\rho(p_1, p_2; \tau_1)Rρ​(p1​,p2​;T)≤Rρ​(p1​,p2​;τ1​) for T∈[0,τ1]T \in [0, \tau_1]T∈[0,τ1​];
  2. Rρ(p1,p2;T)≤Rρ(p1,p2;τ2)R_\rho(p_1, p_2; T) \le R_\rho(p_1, p_2; \tau_2)Rρ​(p1​,p2​;T)≤Rρ​(p1​,p2​;τ2​) for T∈[τ2,1]T \in [\tau_2, 1]T∈[τ2​,1];
  3. Rρ(p1,p2;T)=λ G(p1,p2)R_\rho(p_1, p_2; T) = \lambda\, G(p_1, p_2)Rρ​(p1​,p2​;T)=λG(p1​,p2​) for T∈[τ1,τ2]T \in [\tau_1, \tau_2]T∈[τ1​,τ2​];
  4. for ρ≤e−2+e−1\rho \le e^{-2+e^{-1}}ρ≤e−2+e−1, max⁡G=−e−1+e−1/ln⁡ρ\max G = -e^{-1+e^{-1}}/\ln\rhomaxG=−e−1+e−1/lnρ, attained only at (e−1+e−1,e−2+e−1)(e^{-1+e^{-1}}, e^{-2+e^{-1}})(e−1+e−1,e−2+e−1).

Significance

Proposition 4 gives an explicit optimal markdown policy in a model where segmentation happens purely by arrival time: it shows that the discount time is not pinned down but can be placed anywhere in the interval in which the valuation lies between the two prices, and that for strongly declining valuations the optimal prices do not depend on ρ\rhoρ at all. The paper uses it as the benchmark πC/N∗\pi^*_{C/N}πC/N∗​ against which the strategic-customer equilibrium of Proposition 5 and the losses of Example 1 are measured.

The result is proved in the paper by a short argument; nothing in it has been machine-checked. A formal development makes the three observations of the proof precise (in particular, that prices outside [ρ,1][\rho, 1][ρ,1] are dominated, which the paper leaves implicit) and supplies the omitted calculus for the special case.

Difficulty

The revenue is defined through integrals of a step function of time, and the reduction to GGG needs these integrals evaluated in every configuration of p1p_1p1​, p2p_2p2​ and TTT, including prices above 111 (nobody buys) and below ρ\rhoρ (everyone buys, at a needlessly low price). The paper's proof covers only ρ≤p2≤p1≤1\rho \le p_2 \le p_1 \le 1ρ≤p2​≤p1​≤1 and asserts the domination of the remaining prices without argument. The special case is a constrained two-variable maximization of a function that is not jointly concave; the unconstrained critical point must be shown to be feasible exactly when ρ≤e−2+e−1\rho \le e^{-2+e^{-1}}ρ≤e−2+e−1, and boundary points of the region must be excluded.

Formalization scope

Everything is over R\mathbb RR. Logarithms are Real.log, powers ρT\rho^TρT are real powers, the segment rates are interval integrals ∫ t in a..b of the tail Fˉ(x)=1{x≤1}\bar F(x) = \mathbf 1\{x \le 1\}Fˉ(x)=1{x≤1}, and α=−ln⁡ρ\alpha = -\ln\rhoα=−lnρ with H=1H = 1H=1. The model definitions (ΛI\Lambda_IΛI​, ΛW\Lambda_WΛW​, ΛL\Lambda_LΛL​ and the revenue) are stated for a general tail Fˉ\bar FFˉ, decline factor, season length and discount time and then specialized.

The following readings of the paper's words are fixed:

  • "c=0c = 0c=0": every base valuation equals μ=1\mu = 1μ=1 (the degenerate end of the paper's Gamma family, outside §3's "continuous distribution").
  • "Q/λ→∞Q/\lambda \to \inftyQ/λ→∞": unlimited inventory; the truncated Poisson mean N(q,Λ)N(q, \Lambda)N(q,Λ) is replaced by Λ\LambdaΛ. With unlimited inventory, choosing the contingent discount at time TTT and choosing both prices in advance give the same optimum.
  • Myopic waiting customers buy at TTT iff their valuation at TTT is at least p2p_2p2​, as in ΛW\Lambda_WΛW​.
  • "TTT could be optimally selected": T∈[0,1]T \in [0, 1]T∈[0,1] is a decision variable together with the prices, which range over all 0<p2≤p10 < p_2 \le p_10<p2​≤p1​, not only over [ρ,1][\rho, 1][ρ,1].
  • "Maximize his expected revenues": IsGreatest of the set of attainable revenues.
  • "Setting TTT to any value within the range p2∗≤ρT≤p1∗p_2^* \le \rho^T \le p_1^*p2∗​≤ρT≤p1∗​", and "it would be optimal to select TTT so that ρT∈[e−2+e−1,e−1+e−1]\rho^T \in [e^{-2+e^{-1}}, e^{-1+e^{-1}}]ρT∈[e−2+e−1,e−1+e−1]": every such TTT is optimal; it is not claimed that no other TTT is.
  • "The prices p1∗p_1^*p1∗​ and p2∗p_2^*p2∗​ that solve the problem" in the special case: the maximizer of GGG is unique.
  • "Never optimal" in the first two observations: a weak inequality between revenues.

The decimals 0.1960.1960.196 and 0.5320.5320.532 are not stated. A formalization that restricted prices to [ρ,1][\rho, 1][ρ,1] in the revenue maximization, or that fixed TTT in advance, would assume half of what the proposition proves and is ruled out. Welcome contributions include general lemmas evaluating interval integrals of indicator functions of intervals, and the domination argument for prices outside [ρ,1][\rho, 1][ρ,1].

Selected references

  • Y. Aviv and A. Pazgal, Optimal Pricing of Seasonal Products in the Presence of Forward-Looking Consumers, Manufacturing & Service Operations Management 10(3):339–359, 2008. https://doi.org/10.1287/msom.1070.0183
  • N. Stokey, Intertemporal Price Discrimination, Quarterly Journal of Economics 93(3):355–371, 1979. https://doi.org/10.2307/1883163
  • D. Besanko and W. L. Winston, Optimal Price Skimming by a Monopolist Facing Rational Consumers, Management Science 36(5):555–567, 1990. https://doi.org/10.1287/mnsc.36.5.555
9 thms2 active usersReviewed
🏆Completed
Linear OptimizationOperations ResearchOptimization·Captain: mikedeng1

Revenue Management Under the Markov Chain Choice Model I: The Dual Linear Program Yields an Optimal AssortmentResearch Paper

Motivation

A retailer or an airline decides which products to make available, and customers choose among what is offered. When a preferred product is missing, many customers substitute to another product instead of leaving. Assortment optimization asks which subset of products to offer so that the expected revenue from a customer is as large as possible, and its answer depends entirely on the choice model used to describe substitution.

The Markov chain choice model was introduced by Blanchet, Gallego and Goyal (EC 2013; Oper. Res. 64(4), 2016), who showed that it contains the multinomial logit model as a special case and proposed it as an approximation of general random-utility choice models. Feldman and Topaloglu (Oper. Res. 65(5), 2017) study revenue management under this model. Their first result, the subject of this mission, is that the assortment problem, a search over all 2n2^n2n offer sets, is solved by one linear program with nnn variables. The same paper uses this result for its dynamic single-resource and network results, which are the subjects of the companion missions II–IV of this series.

Timeline:

  • 2013/2016: Blanchet, Gallego and Goyal introduce the model and give a polynomial-time assortment algorithm.
  • 2017: Feldman and Topaloglu show that the optimal assortment is read off an optimal solution of a dual linear program (their Theorem 2), and derive structural and capacity-control consequences.

Setting

There are nnn products N={1,…,n}N=\{1,\dots,n\}N={1,…,n}. A customer arrives to purchase product jjj with probability λj\lambda_jλj​. If the product she visits is offered, she buys it. Otherwise she transitions to product iii with probability ρj,i\rho_{j,i}ρj,i​ and checks whether iii is offered, or leaves without buying with probability 1−∑i∈Nρj,i1-\sum_{i\in N}\rho_{j,i}1−∑i∈N​ρj,i​. Throughout, λj>0\lambda_j>0λj​>0, ρj,i≥0\rho_{j,i}\ge0ρj,i​≥0 and ∑i∈Nρj,i<1\sum_{i\in N}\rho_{j,i}<1∑i∈N​ρj,i​<1 for all jjj.

For an offer set S⊆NS\subseteq NS⊆N, let Pj,SP_{j,S}Pj,S​ be the expected number of visits to product jjj while it is offered (the probability that jjj is purchased) and Rj,SR_{j,S}Rj,S​ the expected number of visits to jjj while it is not offered. The pair (PS,RS)(P_S,R_S)(PS​,RS​) is the unique solution of the (Balance) equations

Pj,S+Rj,S=λj+∑i∈Nρi,jRi,S  ∀j∈N,Pj,S=0  ∀j∉S,Rj,S=0  ∀j∈S.P_{j,S}+R_{j,S}=\lambda_j+\sum_{i\in N}\rho_{i,j}R_{i,S}\ \ \forall j\in N,\qquad P_{j,S}=0\ \ \forall j\notin S,\qquad R_{j,S}=0\ \ \forall j\in S .Pj,S​+Rj,S​=λj​+i∈N∑​ρi,j​Ri,S​  ∀j∈N,Pj,S​=0  ∀j∈/S,Rj,S​=0  ∀j∈S.

With revenue rj∈Rr_j\in\mathbb Rrj​∈R for product jjj, the (Assortment) problem is max⁡S⊆N∑j∈NPj,Srj\max_{S\subseteq N}\sum_{j\in N}P_{j,S}r_jmaxS⊆N​∑j∈N​Pj,S​rj​. Two linear programs enter: the maximization of ∑jrjxj\sum_j r_jx_j∑j​rj​xj​ over the polyhedron

H={(x,z)∈R+2n: xj+zj=λj+∑i∈Nρi,jzi  ∀j∈N},\mathcal H=\Big\{(x,z)\in\mathbb R^{2n}_+:\ x_j+z_j=\lambda_j+\sum_{i\in N}\rho_{i,j}z_i\ \ \forall j\in N\Big\},H={(x,z)∈R+2n​: xj​+zj​=λj​+i∈N∑​ρi,j​zi​  ∀j∈N},

and its dual

min⁡v∈Rn{∑j∈Nλjvj: vj≥rj  ∀j∈N,  vj≥∑i∈Nρj,ivi  ∀j∈N}.(Dual)\min_{v\in\mathbb R^n}\Big\{\sum_{j\in N}\lambda_jv_j:\ v_j\ge r_j\ \ \forall j\in N,\ \ v_j\ge\sum_{i\in N}\rho_{j,i}v_i\ \ \forall j\in N\Big\}.\qquad\text{(Dual)}v∈Rnmin​{j∈N∑​λj​vj​: vj​≥rj​  ∀j∈N,  vj​≥i∈N∑​ρj,i​vi​  ∀j∈N}.(Dual)

Formalization targets

Goal: Theorem 2 (p. 1326)

For every optimal solution v^\hat vv^ of (Dual), the set S^={j∈N:v^j=rj}\hat S=\{j\in N:\hat v_j=r_j\}S^={j∈N:v^j​=rj​} is an optimal assortment:

∑j∈NPj,S rj ≤ ∑j∈NPj,S^ rjfor all S⊆N.\sum_{j\in N}P_{j,S}\,r_j\ \le\ \sum_{j\in N}P_{j,\hat S}\,r_j\qquad\text{for all } S\subseteq N .j∈N∑​Pj,S​rj​ ≤ j∈N∑​Pj,S^​rj​for all S⊆N.

Milestones, in the order the paper uses them

  1. (p. 1325) The (Balance) equations have a unique nonnegative solution for every SSS.
  2. Lemma 1 (p. 1326): for an extreme point (x^,z^)(\hat x,\hat z)(x^,z^) of H\mathcal HH and Sx^={j:x^j>0}S_{\hat x}=\{j:\hat x_j>0\}Sx^​={j:x^j​>0}, Pj,Sx^=x^jP_{j,S_{\hat x}}=\hat x_jPj,Sx^​​=x^j​ and Rj,Sx^=z^jR_{j,S_{\hat x}}=\hat z_jRj,Sx^​​=z^j​ for all jjj.
  3. (p. 1326) The linear program over H\mathcal HH has an optimal solution, and its optimal value equals the optimal value of (Assortment).
  4. (p. 1326) (Dual) has an optimal solution, and its optimal value equals that of the linear program over H\mathcal HH.
  5. (p. 1326, proof of Theorem 2) An optimal v^\hat vv^ satisfies v^j=rj\hat v_j=r_jv^j​=rj​ or v^j=∑iρj,iv^i\hat v_j=\sum_i\rho_{j,i}\hat v_iv^j​=∑i​ρj,i​v^i​ for each jjj.

Significance

Theorem 2 reduces a combinatorial problem over 2n2^n2n offer sets to a linear program, so the assortment problem under the Markov chain choice model is solvable in polynomial time. The same structure drives the rest of the paper: the dual variables v^j\hat v_jv^j​ are used to show that optimal offer sets shrink when all revenues fall by a common amount, to show that the optimal offer sets of the single-resource dynamic program are nested in the remaining capacity and time, and to reduce the choice-based network linear program to a compact one.

The result is proved in the paper. This mission produces a machine-checked version of it and of its supporting lemmas; to the knowledge of the mission author no formalization of the Markov chain choice model exists. The definition layer (the model, the (Balance) solution, H\mathcal HH and (Dual)) is the first formal encoding of this choice model and is shared, with the same encoding, by missions II–IV.

Difficulty

The obvious route, comparing ∑jPj,Srj\sum_jP_{j,S}r_j∑j​Pj,S​rj​ across offer sets directly, fails because Pj,SP_{j,S}Pj,S​ is defined only implicitly through a linear system whose coefficient matrix changes with SSS; there is no closed-form expression that can be compared across offer sets, and enumerating the 2n2^n2n sets is exponential. The connection to linear programming needs an exact correspondence between the vertices of H\mathcal HH and the (Balance) solutions, which is a statement about polyhedra, not about Markov chains. Even well-posedness is not free: existence, uniqueness and nonnegativity of (PS,RS)(P_S,R_S)(PS​,RS​) depend on the row sums ∑iρj,i\sum_i\rho_{j,i}∑i​ρj,i​ being strictly below one. Mathlib has extreme points of convex sets but no ready-made theory of vertices of polyhedra or of linear programming duality in this form.

Formalization scope

Products are Fin n (0-based), offer sets are Finset (Fin n) (the paper's S⊂NS\subset NS⊂N is non-strict inclusion, so every subset including ∅\emptyset∅ and NNN is an offer set), and rho j i is ρj,i\rho_{j,i}ρj,i​, the transition from jjj to iii. The model structure carries the paper's standing assumptions λj>0\lambda_j>0λj​>0 and ∑iρj,i<1\sum_i\rho_{j,i}<1∑i​ρj,i​<1 (p. 1325) and the implicit nonnegativity ρj,i≥0\rho_{j,i}\ge0ρj,i​≥0. Revenues are arbitrary reals with no sign assumption. R2n\mathbb R^{2n}R2n is (Fin n → ℝ) × (Fin n → ℝ), and extreme points are Mathlib's Set.extremePoints ℝ.

The pair (PS,RS)(P_S,R_S)(PS​,RS​) is a solution of (Balance) chosen by Classical.epsilon; milestone 1 states that it satisfies (Balance), is nonnegative and equals every solution, so all statements are about the paper's (PS,RS)(P_S,R_S)(PS​,RS​) and not about a junk value. Every optimum (of (Assortment), of the linear program over H\mathcal HH, of (Dual)) is stated as "feasible and at least as good as every feasible point", never through sSup or sInf. The goal is stated for every optimal v^\hat vv^ of (Dual), which is what the proof uses; uniqueness of v^\hat vv^ is neither assumed nor claimed. Replacing the hypothesis "v^\hat vv^ optimal for (Dual)" by "v^\hat vv^ feasible for (Dual)" would make the statement false, and an existential over v^\hat vv^ would weaken it; neither is the target.

No hypothesis beyond the paper's is added. Two printed index slips on pp. 1325–1326 (a garbled sum in the discussion after (Balance), and ∑iρi,jz^j\sum_i\rho_{i,j}\hat z_j∑i​ρi,j​z^j​ for ∑iρi,jz^i\sum_i\rho_{i,j}\hat z_i∑i​ρi,j​z^i​ in the proof of Lemma 1) are not part of any statement here.

Welcome contributions: the existence and uniqueness of (Balance) solutions via Neumann series for substochastic matrices, a characterization of vertices of polyhedra given by equality constraints and nonnegativity, and a strong duality statement usable for this primal–dual pair. These are reusable well beyond this mission.

Selected references

  • J. B. Feldman, H. Topaloglu, Revenue Management Under the Markov Chain Choice Model, Operations Research 65(5):1322–1342, 2017. https://doi.org/10.1287/opre.2017.1628
  • J. Blanchet, G. Gallego, V. Goyal, A Markov Chain Approximation to Choice Modeling, Operations Research 64(4):886–905, 2016. https://doi.org/10.1287/opre.2016.1505
  • M. L. Puterman, Markov Decision Processes: Discrete Stochastic Dynamic Programming, Wiley, 1994. https://doi.org/10.1002/9780470316887
8 thms2 active usersReviewed
🏆Completed
Dynamic ProgrammingOperations ResearchOptimization·Captain: mikedeng1

Revenue Management Under the Markov Chain Choice Model II: Optimal Offer Sets Grow with Remaining Capacity and Remaining TimeResearch Paper

Motivation

Airlines, hotels and rental firms sell a fixed stock of capacity (seats on a flight leg, rooms on a night) over a finite selling horizon, and control revenue mainly by deciding which products (fare classes, rate plans) to make available at each moment. When customers substitute between products, the decision in each period is an assortment: a subset of products to offer. Talluri and van Ryzin (Management Science 2004) formulated this single resource revenue management problem as a dynamic program for a general choice model, and showed that structural properties of the optimal policy depend strongly on how customers choose.

The Markov chain choice model of Blanchet, Gallego and Goyal (technical report 2013; Oper. Res. 2016) describes substitution by a Markov chain on the products and approximates a broad class of random-utility models. Feldman and Topaloglu (Oper. Res. 65(5), 2017) study assortment and revenue management problems under this model. This mission formalizes their structural result for the single resource problem (Theorem 5): there is an optimal policy whose offer sets are nested in the remaining capacity and in the remaining time, so that it can be implemented by protection levels, one capacity threshold per product and period.

Setting

There are nnn products N={1,…,n}N=\{1,\dots,n\}N={1,…,n}. A customer arrives to purchase product jjj with probability λj\lambda_jλj​. If jjj is offered she buys it; otherwise she moves to product iii with probability ρj,i\rho_{j,i}ρj,i​, and leaves with probability 1−∑i∈Nρj,i1-\sum_{i\in N}\rho_{j,i}1−∑i∈N​ρj,i​. The paper assumes λj>0\lambda_j>0λj​>0 and ∑i∈Nρj,i<1\sum_{i\in N}\rho_{j,i}<1∑i∈N​ρj,i​<1 for all jjj. For an offer set S⊆NS\subseteq NS⊆N, the purchase probabilities Pj,SP_{j,S}Pj,S​ and the visit counts Rj,SR_{j,S}Rj,S​ of unavailable products solve the (Balance) equations

Pj,S+Rj,S=λj+∑i∈Nρi,jRi,S  ∀j,Pj,S=0  (j∉S),Rj,S=0  (j∈S).P_{j,S}+R_{j,S}=\lambda_j+\sum_{i\in N}\rho_{i,j}R_{i,S}\ \ \forall j,\qquad P_{j,S}=0\ \ (j\notin S),\qquad R_{j,S}=0\ \ (j\in S).Pj,S​+Rj,S​=λj​+i∈N∑​ρi,j​Ri,S​  ∀j,Pj,S​=0  (j∈/S),Rj,S​=0  (j∈S).

With product revenues rjr_jrj​, the (Assortment) problem is max⁡S⊆N∑jPj,Srj\max_{S\subseteq N}\sum_{j}P_{j,S}r_jmaxS⊆N​∑j​Pj,S​rj​, and its (Dual) linear program is

min⁡v∈Rn{∑jλjvj:vj≥rj, vj≥∑iρj,ivi  ∀j}.\min_{v\in\mathbb R^n}\Big\{\sum_{j}\lambda_jv_j : v_j\ge r_j,\ v_j\ge\sum_{i}\rho_{j,i}v_i\ \ \forall j\Big\}.v∈Rnmin​{j∑​λj​vj​:vj​≥rj​, vj​≥i∑​ρj,i​vi​  ∀j}.

In the single resource problem there are TTT periods and ccc units of capacity. In each period at most one customer arrives; a sale of product jjj earns rjr_jrj​ and uses one unit. The optimal expected revenue Vt(x)V_t(x)Vt​(x) from period ttt on with xxx units left satisfies the (Single Resource) dynamic program

Vt(x)=max⁡S⊆N{∑j∈NPj,S{rj+Vt+1(x−1)−Vt+1(x)}}+Vt+1(x),V_t(x)=\max_{S\subseteq N}\Big\{\sum_{j\in N}P_{j,S}\{r_j+V_{t+1}(x-1)-V_{t+1}(x)\}\Big\}+V_{t+1}(x),Vt​(x)=S⊆Nmax​{j∈N∑​Pj,S​{rj​+Vt+1​(x−1)−Vt+1​(x)}}+Vt+1​(x),

with VT+1(x)=0V_{T+1}(x)=0VT+1​(x)=0 and Vt(0)=0V_t(0)=0Vt​(0)=0. An optimal subset S^t(x)\hat S_t(x)S^t​(x) is a maximizer of the problem on the right side. The marginal value of capacity is ΔVt(x)=Vt(x)−Vt(x−1)\Delta V_t(x)=V_t(x)-V_t(x-1)ΔVt​(x)=Vt​(x)−Vt​(x−1).

Formalization targets

Goal: Theorem 5 (p. 1329)

There exists an optimal policy, i.e. a choice of maximizers S^t(x)\hat S_t(x)S^t​(x) for all 1≤t≤T1\le t\le T1≤t≤T, 1≤x≤c1\le x\le c1≤x≤c, such that

S^t(x−1)⊆S^t(x)andS^t−1(x)⊆S^t(x).\hat S_t(x-1)\subseteq \hat S_t(x)\qquad\text{and}\qquad \hat S_{t-1}(x)\subseteq\hat S_t(x).S^t​(x−1)⊆S^t​(x)andS^t−1​(x)⊆S^t​(x).

The offer set shrinks as capacity runs down, and it is smaller when more periods remain.

Milestones

  1. (§2, p. 1325) The (Balance) equations have a unique nonnegative solution for every SSS.
  2. (Theorem 2, p. 1326) If v^\hat vv^ is optimal for (Dual), then {j:v^j=rj}\{j:\hat v_j=r_j\}{j:v^j​=rj​} is optimal for (Assortment).
  3. (Lemma 3, p. 1327) For η≥0\eta\ge0η≥0, with v^η\hat v^\etav^η optimal for (Dual) with revenues rj−ηr_j-\etarj​−η,
{j:v^jη=rj−η}⊆{j:v^j0=rj}.\{j:\hat v^\eta_j=r_j-\eta\}\subseteq\{j:\hat v^0_j=r_j\}.{j:v^jη​=rj​−η}⊆{j:v^j0​=rj​}.
  1. ((Single Resource), pp. 1328–1329) The value functions satisfy the printed recursion Vt(x)=max⁡S{∑jPj,S{rj+Vt+1(x−1)}+{1−∑jPj,S}Vt+1(x)}V_t(x)=\max_S\{\sum_jP_{j,S}\{r_j+V_{t+1}(x-1)\}+\{1-\sum_jP_{j,S}\}V_{t+1}(x)\}Vt​(x)=maxS​{∑j​Pj,S​{rj​+Vt+1​(x−1)}+{1−∑j​Pj,S​}Vt+1​(x)} and the boundary conditions.
  2. (Proof of Theorem 5, p. 1329) ΔVt+1(x)≤ΔVt+1(x−1)\Delta V_{t+1}(x)\le\Delta V_{t+1}(x-1)ΔVt+1​(x)≤ΔVt+1​(x−1) and ΔVt+1(x)≤ΔVt(x)\Delta V_{t+1}(x)\le\Delta V_t(x)ΔVt+1​(x)≤ΔVt​(x).

Significance

Theorem 5 makes the optimal policy a protection level policy: for each product jjj and period ttt there is a threshold xˉjt\bar x_{jt}xˉjt​ such that jjj is offered exactly when at least xˉjt\bar x_{jt}xˉjt​ units remain. The policy can then be stored as n Tn\,TnT numbers instead of a table of subsets, and each product can be controlled separately, which is how airline inventory systems are organized. The paper also shows (Table 2, p. 1330) that the products need not be closed in revenue order: a higher-fare product can be closed before a lower-fare one, so the result is genuinely about nested sets, not about nested fare classes. Talluri and van Ryzin (2004) show that under the multinomial logit model an optimal assortment consists of a number of products with the largest revenues; the Markov chain choice model does not have this property (Table 1, p. 1328), so their structure of the optimal policy does not carry over directly.

The result is proved in the paper. No part of it is formalized on Prove2Me. The monotonicity of marginal values for an abstract choice model is published and proved as RevenueManagement.choice_marginal_values, stated over the definitions RevenueManagement_singleResource; milestone 5 is the same statement for this paper's dynamic program and can be bridged to it. A complete development here produces machine-checked Theorem 2 and Lemma 3 for the Markov chain choice model, which are reusable for any assortment problem under this model.

Difficulty

The obvious argument chooses, for each (t,x)(t,x)(t,x), any maximizer of the stage problem. This fails: the stage problems have ties, and an arbitrary choice of maximizers need not be nested, so the theorem is an existence statement about a coordinated choice. The dependence of Pj,SP_{j,S}Pj,S​ on SSS is through the solution of a linear system, so the effect of adding or removing a product on the other purchase probabilities has no simple sign, and optimal assortments need not be nested by revenue (Table 1, p. 1328). The link between two stage problems is that their revenues differ by the same constant for all products, and what has to be shown is that such a uniform shift moves an optimal assortment in a controlled direction. Revenues in the stage problems, rj−ΔVt+1(x)r_j-\Delta V_{t+1}(x)rj​−ΔVt+1​(x), can be negative.

Formalization scope

Products are Fin n, offer sets Finset (Fin n), the paper's ⊂\subset⊂ is non-strict inclusion ⊆\subseteq⊆. The model is a structure with fields λ\lambdaλ, ρ\rhoρ and the standing assumptions λj>0\lambda_j>0λj​>0, ∑iρj,i<1\sum_i\rho_{j,i}<1∑i​ρj,i​<1; ρj,i≥0\rho_{j,i}\ge0ρj,i​≥0 is added as a field because the ρj,i\rho_{j,i}ρj,i​ are probabilities. (PS,RS)(P_S,R_S)(PS​,RS​) is a solution of (Balance) chosen by Classical.epsilon; milestone 1 states that it is the unique nonnegative one. The value function is defined by recursion on the number of periods to go, and VtV_tVt​ is used only for 1≤t≤T+11\le t\le T+11≤t≤T+1. Maxima over S⊆NS\subseteq NS⊆N are taken over the finite family of all subsets, so they are attained; dual optimality is "feasible and no worse than every feasible point".

One hypothesis is added to Theorem 5 and milestone 5, and disclosed: ∑jλj≤1\sum_j\lambda_j\le1∑j​λj​≤1, implied by the paper's description of at most one arrival per period (it makes 1−∑jPj,S1-\sum_jP_{j,S}1−∑j​Pj,S​ a probability). No sign condition is placed on the revenues, as on the page. Theorem 2 and Lemma 3 are stated for arbitrary real revenues, as the proof of Theorem 5 applies them to rj−ΔVt+1(x)r_j-\Delta V_{t+1}(x)rj​−ΔVt+1​(x). The first inclusion of Theorem 5 is stated for x≥2x\ge2x≥2: with x−1=0x-1=0x−1=0 units there is no decision, and the paper reads S^t(0)\hat S_t(0)S^t​(0) as ∅\emptyset∅. The existential requires every S^t(x)\hat S_t(x)S^t​(x) in range to be optimal for its stage problem; a statement without that clause would be satisfied by the empty sets and is ruled out.

Theorem 2 and Lemma 3 need LP duality for (Dual) and the (Balance) system; milestone 5 needs the standard induction on the dynamic program, or a bridge to the published choice_marginal_values (with arrival probabilities 111 and choice model Pj,SP_{j,S}Pj,S​). Proofs of any milestone, and bridges to published Mathlib or platform LP duality results, are welcome.

Selected references

  • J. B. Feldman, H. Topaloglu, Revenue Management Under the Markov Chain Choice Model, Operations Research 65(5):1322–1342, 2017. https://doi.org/10.1287/opre.2017.1628
  • J. Blanchet, G. Gallego, V. Goyal, A Markov Chain Approximation to Choice Modeling, Operations Research 64(4):886–905, 2016. https://doi.org/10.1287/opre.2016.1505
  • K. Talluri, G. van Ryzin, Revenue Management Under a General Discrete Choice Model of Consumer Behavior, Management Science 50(1):15–33, 2004. https://doi.org/10.1287/mnsc.1030.0147
  • K. Talluri, G. van Ryzin, The Theory and Practice of Revenue Management, Springer, 2004. https://doi.org/10.1007/b139000
10 thms2 active usersReviewed
🏆Completed
Linear OptimizationOperations ResearchOptimization·Captain: mikedeng1

Revenue Management Under the Markov Chain Choice Model III: The Reduced Linear Program Is Equivalent to the Choice-Based Linear ProgramResearch Paper

Motivation

Network revenue management decides which products to make available to arriving customers when products share scarce resources: an airline sells itineraries (products) that consume seats on flight legs (resources), and each itinerary has a fare. When customers choose among the offered products, rather than asking for one fixed product, the standard planning tool is a deterministic linear program that replaces random choices by their expected values. Gallego, Iyengar, Phillips and Dubey (2004, Columbia CORC technical report TR-2004-01) and Liu and van Ryzin (2008) formulated this choice-based linear program; its solution drives bid-price and offer-set policies used in practice.

The difficulty is size. The choice-based program has one variable for each subset of products, 2n2^n2n in all, and is solved by column generation, whose pricing subproblem is itself an assortment problem. Feldman and Topaloglu (Oper. Res. 65(5), 2017) show that when customers choose under the Markov chain choice model of Blanchet, Gallego and Goyal (2016), the choice-based program is equivalent to a linear program with only 2n2n2n variables and m+nm+nm+n constraints. This mission formalizes that equivalence, Theorem 7 of the paper.

Setting

There are nnn products, N={1,…,n}N=\{1,\dots,n\}N={1,…,n}. Under the Markov chain choice model a customer first visits product jjj with probability λj\lambda_jλj​. If the product she visits is offered, she buys it. Otherwise she moves from product jjj to product iii with probability ρj,i\rho_{j,i}ρj,i​, or leaves without buying with probability 1−∑i∈Nρj,i1-\sum_{i\in N}\rho_{j,i}1−∑i∈N​ρj,i​. The paper assumes throughout that

λj>0and∑i∈Nρj,i<1for all j∈N.\lambda_j>0\quad\text{and}\quad\sum_{i\in N}\rho_{j,i}<1\qquad\text{for all } j\in N.λj​>0andi∈N∑​ρj,i​<1for all j∈N.

For an offer set S⊆NS\subseteq NS⊆N, Pj,SP_{j,S}Pj,S​ is the expected number of visits to product jjj while it is offered, which is its purchase probability, and Rj,SR_{j,S}Rj,S​ is the expected number of visits to jjj while it is not offered. The pair (PS,RS)(P_S,R_S)(PS​,RS​) is the solution of the (Balance) equations

Pj,S+Rj,S=λj+∑i∈Nρi,jRi,S  ∀j∈N,Pj,S=0  ∀j∉S,Rj,S=0  ∀j∈S.P_{j,S}+R_{j,S}=\lambda_j+\sum_{i\in N}\rho_{i,j}R_{i,S}\ \ \forall j\in N,\qquad P_{j,S}=0\ \ \forall j\notin S,\qquad R_{j,S}=0\ \ \forall j\in S .Pj,S​+Rj,S​=λj​+i∈N∑​ρi,j​Ri,S​  ∀j∈N,Pj,S​=0  ∀j∈/S,Rj,S​=0  ∀j∈S.

Dropping the constraints tied to SSS gives the polyhedron

H={(x,z)∈R+2n:xj+zj=λj+∑i∈Nρi,jzi  ∀j∈N}.\mathcal H=\Big\{(x,z)\in\mathbb R^{2n}_+ : x_j+z_j=\lambda_j+\sum_{i\in N}\rho_{i,j}z_i\ \ \forall j\in N\Big\}.H={(x,z)∈R+2n​:xj​+zj​=λj​+i∈N∑​ρi,j​zi​  ∀j∈N}.

The network has mmm resources, M={1,…,m}M=\{1,\dots,m\}M={1,…,m}, with capacities cqc_qcq​; the selling horizon has TTT periods; product jjj earns rjr_jrj​ and consumes aq,ja_{q,j}aq,j​ units of resource qqq. With uSu_SuS​ the probability of offering SSS in a period, the (Choice Based) linear program is

max⁡u∈R+2n{∑S⊆N∑j∈NTrjPj,SuS: ∑S⊆N∑j∈NTaq,jPj,SuS≤cq ∀q∈M, ∑S⊆NuS=1},\max_{u\in\mathbb R^{2^n}_+}\Big\{\sum_{S\subseteq N}\sum_{j\in N}T r_jP_{j,S}u_S:\ \sum_{S\subseteq N}\sum_{j\in N}Ta_{q,j}P_{j,S}u_S\le c_q\ \forall q\in M,\ \sum_{S\subseteq N}u_S=1\Big\},u∈R+2n​max​{S⊆N∑​j∈N∑​Trj​Pj,S​uS​: S⊆N∑​j∈N∑​Taq,j​Pj,S​uS​≤cq​ ∀q∈M, S⊆N∑​uS​=1},

and the (Reduced) linear program is

max⁡(x,z)∈R+2n{∑j∈NTrjxj: ∑j∈NTaq,jxj≤cq ∀q∈M, xj+zj=λj+∑i∈Nρi,jzi ∀j∈N}.\max_{(x,z)\in\mathbb R^{2n}_+}\Big\{\sum_{j\in N}T r_jx_j:\ \sum_{j\in N}Ta_{q,j}x_j\le c_q\ \forall q\in M,\ x_j+z_j=\lambda_j+\sum_{i\in N}\rho_{i,j}z_i\ \forall j\in N\Big\}.(x,z)∈R+2n​max​{j∈N∑​Trj​xj​: j∈N∑​Taq,j​xj​≤cq​ ∀q∈M, xj​+zj​=λj​+i∈N∑​ρi,j​zi​ ∀j∈N}.

In (Reduced), xjx_jxj​ is the expected number of visits to product jjj while it is available and zjz_jzj​ the expected number of visits while it is not.

Formalization targets

Goal: Theorem 7

Let (x^,z^)(\hat x,\hat z)(x^,z^) be an optimal solution of (Reduced). Then there are subsets S1,…,SK⊆NS^1,\dots,S^K\subseteq NS1,…,SK⊆N and positive scalars γ1,…,γK\gamma^1,\dots,\gamma^Kγ1,…,γK with ∑kγk=1\sum_k\gamma^k=1∑k​γk=1 such that

x^=∑k=1KγkPSk,z^=∑k=1KγkRSk,\hat x=\sum_{k=1}^K\gamma^kP_{S^k},\qquad \hat z=\sum_{k=1}^K\gamma^kR_{S^k},x^=k=1∑K​γkPSk​,z^=k=1∑K​γkRSk​,

and for any such subsets and scalars the vector u^\hat uu^ with u^Sk=γk\hat u_{S^k}=\gamma^ku^Sk​=γk and u^S=0\hat u_S=0u^S​=0 for S∉{S1,…,SK}S\notin\{S^1,\dots,S^K\}S∈/{S1,…,SK} is optimal for (Choice Based), with objective value equal to that of (x^,z^)(\hat x,\hat z)(x^,z^) in (Reduced). In particular the two programs have the same optimal value.

Milestones

  1. (Balance) has a unique and nonnegative solution for every offer set (§2, p. 1325).
  2. Lemma 1: an extreme point (x^,z^)(\hat x,\hat z)(x^,z^) of H\mathcal HH equals (PS,RS)(P_{S},R_{S})(PS​,RS​) for S={j:x^j>0}S=\{j:\hat x_j>0\}S={j:x^j​>0} (p. 1326).
  3. Lemma 10: H\mathcal HH is bounded (quoted on p. 1331; proved in the online appendix).
  4. Every point of H\mathcal HH is a positive convex combination of finitely many extreme points of H\mathcal HH (proof of Theorem 7, p. 1331).
  5. Every feasible uuu of (Choice Based) yields the feasible point x~j=∑SPj,SuS\tilde x_j=\sum_S P_{j,S}u_Sx~j​=∑S​Pj,S​uS​, z~j=∑SRj,SuS\tilde z_j=\sum_S R_{j,S}u_Sz~j​=∑S​Rj,S​uS​ of (Reduced), with the same objective value (proof of Theorem 7, p. 1332).

Significance

The result. Theorem 7 replaces a program with 2n2^n2n columns by one with 2n2n2n variables and m+nm+nm+n constraints, solvable directly by any LP solver, and it returns an optimal solution of the original program, not only its value. The optimal value is the standard upper bound on the optimal expected revenue of a network revenue management policy, and the dual variables of the capacity constraints are the bid prices used to control sales. The decomposition of part 1 is what turns the small program's solution back into offer-set frequencies that a policy can implement; Section 7 of the paper makes that decomposition algorithmic (a separate mission of this series).

Formalizing it. The theorem is proved in the paper; to the best of available knowledge no machine-checked version exists. A formal proof checks the link between polyhedral geometry (extreme points of H\mathcal HH and the solutions of (Balance)) and linear-programming optimality, and records exactly which properties of the Markov chain choice model are used: the standing assumptions enter through uniqueness and nonnegativity of (PS,RS)(P_S,R_S)(PS​,RS​) and through boundedness of H\mathcal HH.

Difficulty

The inequality "(Reduced) ≥\ge≥ (Choice Based)" is a direct computation: averaging the (Balance) equations with weights uSu_SuS​ lands in H\mathcal HH. The reverse direction is where the obvious argument fails. A point of H\mathcal HH has no offer set attached to it, and a general polyhedron need not be the convex hull of its extreme points: it can contain lines or rays. The argument requires that H\mathcal HH is bounded, which depends on the substochasticity ∑iρj,i<1\sum_i\rho_{j,i}<1∑i​ρj,i​<1, and that each extreme point is exactly some (PS,RS)(P_S,R_S)(PS​,RS​), which uses the structure of the balance equations and the uniqueness of their solution. Neither follows from general linear-programming facts. In Lean, the finite vertex representation of a bounded polyhedron is also not a one-line consequence of Mathlib's Krein–Milman theorem, which gives only the closure of the convex hull.

Formalization scope

Products are Fin n, offer sets Finset (Fin n), resources Fin m. The model is a structure Model n holding λ\lambdaλ, ρ\rhoρ (rho j i =ρj,i=\rho_{j,i}=ρj,i​, the transition from jjj to iii) and the standing assumptions λj>0\lambda_j>0λj​>0 and ∑iρj,i<1\sum_i\rho_{j,i}<1∑i​ρj,i​<1; the nonnegativity ρj,i≥0\rho_{j,i}\ge0ρj,i​≥0, implicit in the paper because the ρj,i\rho_{j,i}ρj,i​ are probabilities, is an added field. (PS,RS)(P_S,R_S)(PS​,RS​) is a solution of (Balance) chosen by Classical.epsilon, and milestone 1 is what identifies it with the paper's unique solution. H\mathcal HH is a subset of (Fin n → ℝ) × (Fin n → ℝ), extreme points are Mathlib's Set.extremePoints ℝ, and boundedness is Bornology.IsBounded. TTT is a natural number entering as a real factor exactly where the paper writes it; ccc, aaa, rrr carry no sign conditions, as in the paper. "Optimal solution" means feasible and at least as good as every feasible point; no supremum is used.

Two conventions are disclosed. The paper defines u^\hat uu^ by u^Sk=γk\hat u_{S^k}=\gamma^ku^Sk​=γk; the formalization sets u^S=∑k:Sk=Sγk\hat u_S=\sum_{k:S^k=S}\gamma^ku^S​=∑k:Sk=S​γk, which agrees when the SkS^kSk are distinct and is the only consistent reading otherwise. Milestone 5 is stated for every feasible uuu of (Choice Based), while the paper applies it to an optimal one; its argument uses only feasibility. No printed statement needed correction.

The goal is not only the decomposition of part 1, which is milestones 2 and 4 combined: it also asserts optimality of u^\hat uu^ and equality of the optimal values, and a formalization that drops part 2 does not state Theorem 7. Part 2 is required for every decomposition, and part 1 guarantees one exists, so part 2 is not vacuous.

A complete development needs the vertex representation of polytopes (reusable beyond this mission; LinearOptimization.polyhedron_resolution on the platform proves the resolution theorem in another encoding), the theory of substochastic matrices behind (Balance) (invertibility of I−QˉI-\bar QI−Qˉ​ with a nonnegative inverse), and finite-sum manipulations over Finset (Fin n). Contributions to any milestone, and bridges to existing polyhedral results, are welcome.

Selected references

  • J. B. Feldman, H. Topaloglu, Revenue Management Under the Markov Chain Choice Model, Operations Research 65(5):1322–1342, 2017. https://doi.org/10.1287/opre.2017.1628
  • J. Blanchet, G. Gallego, V. Goyal, A Markov Chain Approximation to Choice Modeling, Operations Research 64(4):886–905, 2016. https://doi.org/10.1287/opre.2016.1505
  • Q. Liu, G. van Ryzin, On the Choice-Based Linear Programming Model for Network Revenue Management, Manufacturing & Service Operations Management 10(2):288–310, 2008. https://doi.org/10.1287/msom.1070.0172
  • G. Gallego, G. Iyengar, R. Phillips, A. Dubey, Managing Flexible Products on a Network, Computational Optimization Research Center Technical Report TR-2004-01, Columbia University, 2004 (technical report; no DOI).
  • M. L. Puterman, Markov Decision Processes: Discrete Stochastic Dynamic Programming, Wiley, 1994. https://doi.org/10.1002/9780470316887
8 thms2 active usersReviewed
🏆Completed
Linear OptimizationOperations ResearchOptimization·Captain: mikedeng1

Revenue Management Under the Markov Chain Choice Model IV: Dimension Reduction Recovers a Choice-Based Solution with at Most n+1 Nested Offer SetsResearch Paper

Motivation

In network revenue management a firm sells nnn products that consume mmm capacitated resources over TTT periods, and in each period it chooses which offer set of products to make available. When customers choose among the offered products according to a choice model, the standard planning tool is the choice-based linear program, which assigns a frequency uS≥0u_S\ge 0uS​≥0 to every offer set S⊆NS\subseteq NS⊆N (Liu and van Ryzin 2008). It has 2n2^n2n variables. Its optimal solution is what the firm actually implements: a randomization over offer sets.

For the Markov chain choice model of Blanchet, Gallego and Goyal (2016), Feldman and Topaloglu (2017) show that the choice-based program is equivalent to a reduced linear program with only 2n2n2n variables. The reduced program can be solved directly, but its solution is a pair of vectors (x^,z^)(\hat x,\hat z)(x^,z^), not a randomization over offer sets. This mission formalizes the paper's Section 7, which converts the reduced solution back into a randomization over offer sets: the Dimension Reduction algorithm, and Theorem 9, which says that the algorithm produces at most n+1n+1n+1 offer sets and that these sets are nested.

Setting

Products are indexed by N={1,…,n}N=\{1,\dots,n\}N={1,…,n}. A customer arrives wanting product jjj with probability λj\lambda_jλj​. If the product she visits is offered, she buys it. Otherwise she moves to product iii with probability ρj,i\rho_{j,i}ρj,i​, and leaves without buying with probability 1−∑i∈Nρj,i1-\sum_{i\in N}\rho_{j,i}1−∑i∈N​ρj,i​. Throughout, λj>0\lambda_j>0λj​>0, ρj,i≥0\rho_{j,i}\ge0ρj,i​≥0 and ∑i∈Nρj,i<1\sum_{i\in N}\rho_{j,i}<1∑i∈N​ρj,i​<1 for all j,ij,ij,i.

For an offer set S⊆NS\subseteq NS⊆N, the (Balance) equations

Pj,S+Rj,S=λj+∑i∈Nρi,jRi,S  ∀j∈N,Pj,S=0  ∀j∉S,Rj,S=0  ∀j∈SP_{j,S}+R_{j,S}=\lambda_j+\sum_{i\in N}\rho_{i,j}R_{i,S}\ \ \forall j\in N,\qquad P_{j,S}=0\ \ \forall j\notin S,\qquad R_{j,S}=0\ \ \forall j\in SPj,S​+Rj,S​=λj​+i∈N∑​ρi,j​Ri,S​  ∀j∈N,Pj,S​=0  ∀j∈/S,Rj,S​=0  ∀j∈S

have a unique solution. Pj,SP_{j,S}Pj,S​ is the probability that product jjj is purchased, and Rj,SR_{j,S}Rj,S​ is the expected number of visits to jjj while it is not offered. Write PS=(Pj,S)jP_S=(P_{j,S})_jPS​=(Pj,S​)j​ and RS=(Rj,S)jR_S=(R_{j,S})_jRS​=(Rj,S​)j​.

The set H\mathcal HH consists of all pairs (x,z)∈Rn×Rn(x,z)\in\mathbb R^n\times\mathbb R^n(x,z)∈Rn×Rn with x,z≥0x,z\ge0x,z≥0 and xj+zj=λj+∑iρi,jzix_j+z_j=\lambda_j+\sum_{i}\rho_{i,j}z_ixj​+zj​=λj​+∑i​ρi,j​zi​ for all jjj. The (Reduced) linear program maximizes ∑jTrjxj\sum_j T r_j x_j∑j​Trj​xj​ over (x,z)∈H(x,z)\in\mathcal H(x,z)∈H subject to the capacity rows ∑jTaq,jxj≤cq\sum_j T a_{q,j}x_j\le c_q∑j​Taq,j​xj​≤cq​ for each resource qqq.

The Dimension Reduction algorithm starts from the optimal solution (x^,z^)(\hat x,\hat z)(x^,z^) of (Reduced) and sets (x1,z1)=(x^,z^)(x^1,z^1)=(\hat x,\hat z)(x1,z1)=(x^,z^). At iteration kkk it forms Sk={j:xjk>0}S^k=\{j: x^k_j>0\}Sk={j:xjk​>0}. It stops with αk=1\alpha^k=1αk=1 if Sk=∅S^k=\emptysetSk=∅. Otherwise it sets αk=min⁡{xjk/Pj,Sk:j∈Sk}\alpha^k=\min\{x^k_j/P_{j,S^k}: j\in S^k\}αk=min{xjk​/Pj,Sk​:j∈Sk}, stops if αk=1\alpha^k=1αk=1, and if not updates

xjk+1=xjk−αkPj,Sk1−αk,zjk+1=zjk−αkRj,Sk1−αk.x^{k+1}_j=\frac{x^k_j-\alpha^kP_{j,S^k}}{1-\alpha^k},\qquad z^{k+1}_j=\frac{z^k_j-\alpha^kR_{j,S^k}}{1-\alpha^k}.xjk+1​=1−αkxjk​−αkPj,Sk​​,zjk+1​=1−αkzjk​−αkRj,Sk​​.

The weights are γk=(1−α1)⋯(1−αk−1)αk\gamma^k=(1-\alpha^1)\cdots(1-\alpha^{k-1})\alpha^kγk=(1−α1)⋯(1−αk−1)αk.

Formalization targets

Goal: Theorem 9

Let (x^,z^)(\hat x,\hat z)(x^,z^) be optimal for (Reduced). The algorithm stops at some first iteration KKK with 1≤K≤n+11\le K\le n+11≤K≤n+1. At that KKK,

x^=∑k=1KγkPSk,z^=∑k=1KγkRSk,∑k=1Kγk=1,S1⊇S2⊇⋯⊇SK.\hat x=\sum_{k=1}^K\gamma^kP_{S^k},\qquad \hat z=\sum_{k=1}^K\gamma^kR_{S^k},\qquad \sum_{k=1}^K\gamma^k=1,\qquad S^1\supseteq S^2\supseteq\cdots\supseteq S^K .x^=k=1∑K​γkPSk​,z^=k=1∑K​γkRSk​,k=1∑K​γk=1,S1⊇S2⊇⋯⊇SK.

Milestones

In the order the proof uses them:

  1. (Balance) has a unique, nonnegative solution for every SSS (p. 1325).
  2. Pj,S>0P_{j,S}>0Pj,S​>0 for j∈Sj\in Sj∈S (p. 1325).
  3. For (x^,z^)∈H(\hat x,\hat z)\in\mathcal H(x^,z^)∈H: z^j≥Rj,Sx^\hat z_j\ge R_{j,S_{\hat x}}z^j​≥Rj,Sx^​​. This is Lemma 11 of the online appendix, as quoted on p. 1332.
  4. For (x^,z^)∈H(\hat x,\hat z)\in\mathcal H(x^,z^)∈H with nonempty support: αx^≤1\alpha_{\hat x}\le1αx^​≤1 (p. 1332).
  5. The two stopping configurations are (Balance) solutions: Sx^=∅S_{\hat x}=\emptysetSx^​=∅ gives (P∅,R∅)(P_\emptyset,R_\emptyset)(P∅​,R∅​), and αx^=1\alpha_{\hat x}=1αx^​=1 gives (PSx^,RSx^)(P_{S_{\hat x}},R_{S_{\hat x}})(PSx^​​,RSx^​​) (p. 1332).
  6. Lemma 8: one step of the algorithm maps H\mathcal HH into H\mathcal HH and removes a minimizing index from the support (p. 1333).
  7. The iterates stay in H\mathcal HH, and Sk+1⊆Sk∖{jk}S^{k+1}\subseteq S^k\setminus\{j^k\}Sk+1⊆Sk∖{jk} (pp. 1333–1334).

Significance

Theorem 9 makes the reduced program usable. Together with the equivalence of (Reduced) and the choice-based program, it produces an optimal choice-based solution that offers at most n+1n+1n+1 sets. By a basic-solution count the choice-based program has an optimal solution with at most m+1m+1m+1 nonzero offer sets. Theorem 9 bounds the number by the number of products instead, and adds a structural property: the offered sets form a chain. Each set is contained in the previous one, so the implemented policy offers a single nested sequence of assortments. The algorithm needs at most n+1n+1n+1 solutions of linear systems of size nnn. It never enumerates offer sets.

The result is proved in the paper. As far as a search of the platform and of Mathlib shows, it has not been formalized. Mathlib's Carathéodory theorem (Analysis/Convex/Caratheodory) bounds the number of points in a convex combination. It does not construct this decomposition, and it does not give nestedness. Formalizing Theorem 9 therefore requires the paper's argument: a verified, terminating peeling procedure on the polyhedron H\mathcal HH whose output is a convex combination of (Balance) solutions indexed by a chain of sets.

Difficulty

The algorithm divides twice. Step 2 divides by Pj,SkP_{j,S^k}Pj,Sk​, which must be positive on SkS^kSk. Step 3 divides by 1−αk1-\alpha^k1−αk, which must be nonzero whenever the algorithm continues. Both rest on nontrivial facts: the positivity of purchase probabilities, and the bound α≤1\alpha\le1α≤1 on H\mathcal HH. That bound in turn requires the comparison z^≥RSx^\hat z\ge R_{S_{\hat x}}z^≥RSx^​​. The paper proves this comparison only in its online appendix; it is a monotonicity property of the inverse of I−ρ⊤I-\rho^\topI−ρ⊤ restricted to the unoffered products. The update of Step 3 can make coordinates of zzz negative unless this comparison holds, so the naive observation that "(xk,zk)(x^k,z^k)(xk,zk) is a combination of points of H\mathcal HH" does not by itself keep the iterates in H\mathcal HH.

Termination is not a consequence of αk<1\alpha^k<1αk<1 alone. It needs the strict decrease of the support, which comes from choosing a minimizing index. The weights γk\gamma^kγk are defined by products of the (1−αl)(1-\alpha^l)(1−αl). Their sum telescopes to 111 only because the last step has αK=1\alpha^K=1αK=1.

Formalization scope

Everything lives in the namespace MarkovChainChoice.DimReduction. Products are Fin n and offer sets are Finset (Fin n). The model is a structure carrying λj>0\lambda_j>0λj​>0, ∑iρj,i<1\sum_i\rho_{j,i}<1∑i​ρj,i​<1 and, as a disclosed implicit hypothesis, ρj,i≥0\rho_{j,i}\ge0ρj,i​≥0. The pair (PS,RS)(P_S,R_S)(PS​,RS​) is chosen by Classical.epsilon among the solutions of (Balance); milestone 1 is what identifies it with the paper's unique solution. Revenues, capacities and consumptions carry no sign conditions. "Optimal for (Reduced)" means feasible and at least as good as every feasible point.

The iterates are indexed from k=1k=1k=1 as in the paper. The value at index 000 is an unused copy of the input. "The algorithm stops at iteration KKK" is the first K≥1K\ge1K≥1 with SK=∅S^K=\emptysetSK=∅ or αK=1\alpha^K=1αK=1. Iterates after the first stop involve a division by zero, which Lean evaluates to 000, and no statement refers to them. The termination bound K≤n+1K\le n+1K≤n+1 is a conclusion of the goal, not a hypothesis. The goal is about the algorithm's actual iterates, not about an arbitrary sequence that satisfies the invariants. The paper's ⊂\subset⊂ and ⊃\supset⊃ denote non-strict inclusion and are rendered as ⊆\subseteq⊆ and ⊇\supseteq⊇.

Two milestones are stated more generally than the page. The paper states the iteration invariants for the (Reduced) optimum; the Lean requires only (x^,z^)∈H(\hat x,\hat z)\in\mathcal H(x^,z^)∈H, which every feasible (Reduced) solution satisfies. Lemma 8's arg⁡min⁡\arg\minargmin may be a set; the statement holds for every minimizer. The page contains printed slips: on p. 1332, v^j\hat v_jv^j​ has numerator x^j−αRj,S\hat x_j-\alpha R_{j,S}x^j​−αRj,S​, and the paper writes Rx^R_{\hat x}Rx^​ for RSx^R_{S_{\hat x}}RSx^​​. The Lean uses the correct forms, as in Lemma 8 and Step 3. None of these slips appears in a milestone quotation.

A complete development needs a small theory of M-matrices, in the form of nonnegativity of (I−Q)−1(I-Q)^{-1}(I−Q)−1 for substochastic QQQ, together with finite induction on supports. The comparison lemma (milestone 3) and the (Balance) existence and uniqueness result (milestone 1) are reusable across the other missions on this paper. Proofs of any milestone are welcome, as are alternative proofs of the goal.

Selected references

  • J. B. Feldman, H. Topaloglu, Revenue Management Under the Markov Chain Choice Model, Operations Research 65(5):1322–1342, 2017. https://doi.org/10.1287/opre.2017.1628
  • J. Blanchet, G. Gallego, V. Goyal, A Markov Chain Approximation to Choice Modeling, Operations Research 64(4):886–905, 2016. https://doi.org/10.1287/opre.2016.1505
  • Q. Liu, G. van Ryzin, On the Choice-Based Linear Programming Model for Network Revenue Management, Manufacturing & Service Operations Management 10(2):288–310, 2008. https://doi.org/10.1287/msom.1070.0169
14 thms2 active usersReviewed
🏆Completed
Algorithmic Game TheoryLinear OptimizationOperations Research+1·Captain: mikedeng1

Online Primal-Dual Algorithms for Maximizing Ad-Auctions Revenue: The Competitive Ratio of the Primal-Dual Allocation AlgorithmResearch Paper

Motivation

Search engines sell advertisement slots next to their results through ad-auctions. Advertisers bid on keywords, and each advertiser also sets a daily budget: the most it is willing to pay in a day. Queries arrive one at a time and each must be assigned to an advertiser at once, with no knowledge of the queries still to come. The seller's revenue from an advertiser is capped by its budget, so an allocation rule that ignores budgets can exhaust a high bidder early and forgo revenue that a more even allocation would have collected. The question is how much of the offline optimum an online rule can guarantee against every arrival sequence.

Mehta, Saberi, Vazirani and Vazirani (FOCS 2005 / J. ACM 2007) gave a deterministic algorithm whose competitive ratio tends to 1−1/e1 - 1/e1−1/e when bids are small compared with budgets, and showed that no deterministic algorithm does better. Their algorithm builds on online bipartite matching (Karp, Vazirani and Vazirani, STOC 1990) and online bbb-matching (Kalyanasundaram and Pruhs, 2000). Buchbinder, Jain and Naor (ESA 2007) rederived the 1−1/e1 - 1/e1−1/e bound with an online primal-dual algorithm, which gives the ratio in closed form for every value of the bid-to-budget ratio and extends to multiple slots, stochastic information, bounded degree and budget flexibility. This mission formalizes the basic algorithm of that paper and its Theorem 1.

Setting

There is a finite nonempty set III of buyers. Buyer iii has a known budget B(i)>0B(i) > 0B(i)>0. Products j=1,…,mj = 1, \dots, mj=1,…,m arrive one by one; when product jjj arrives, every buyer's bid b(i,j)≥0b(i,j) \ge 0b(i,j)≥0 on it is revealed. The bid-to-budget ratio is

Rmax⁡=max⁡i∈I, jb(i,j)B(i).R_{\max} = \max_{i \in I,\, j} \frac{b(i,j)}{B(i)} .Rmax​=i∈I,jmax​B(i)b(i,j)​.

A fractional allocation y(i,j)≥0y(i,j) \ge 0y(i,j)≥0 assigns fractions of products to buyers; the revenue from buyer iii is the minimum of ∑jb(i,j) y(i,j)\sum_j b(i,j)\,y(i,j)∑j​b(i,j)y(i,j) and B(i)B(i)B(i).

The offline fractional problem is the packing LP, which the paper calls the dual:

max⁡∑j∑ib(i,j) y(i,j)s.t.∑iy(i,j)≤1  ∀j,∑jb(i,j) y(i,j)≤B(i)  ∀i,y≥0.\max \sum_{j}\sum_{i} b(i,j)\,y(i,j) \quad\text{s.t.}\quad \sum_i y(i,j) \le 1 \ \ \forall j,\qquad \sum_j b(i,j)\,y(i,j) \le B(i)\ \ \forall i,\qquad y \ge 0 .maxj∑​i∑​b(i,j)y(i,j)s.t.i∑​y(i,j)≤1  ∀j,j∑​b(i,j)y(i,j)≤B(i)  ∀i,y≥0.

Its LP dual, the paper's primal, is the covering LP:

min⁡∑iB(i) x(i)+∑jz(j)s.t.b(i,j) x(i)+z(j)≥b(i,j)  ∀i,j,x,z≥0.\min \sum_i B(i)\,x(i) + \sum_j z(j) \quad\text{s.t.}\quad b(i,j)\,x(i) + z(j) \ge b(i,j)\ \ \forall i,j,\qquad x, z \ge 0 .mini∑​B(i)x(i)+j∑​z(j)s.t.b(i,j)x(i)+z(j)≥b(i,j)  ∀i,j,x,z≥0.

The Allocation Algorithm has a parameter c>1c > 1c>1 and starts from x≡0x \equiv 0x≡0. When product jjj arrives it takes a buyer iii maximizing b(i,j)(1−x(i))b(i,j)(1 - x(i))b(i,j)(1−x(i)). If x(i)≥1x(i) \ge 1x(i)≥1, the product is not sold. Otherwise it charges iii the minimum of b(i,j)b(i,j)b(i,j) and iii's remaining budget, sets y(i,j)←1y(i,j) \leftarrow 1y(i,j)←1 and z(j)←b(i,j)(1−x(i))z(j) \leftarrow b(i,j)(1 - x(i))z(j)←b(i,j)(1−x(i)), and updates

x(i)←x(i)(1+b(i,j)B(i))+b(i,j)(c−1) B(i).x(i) \leftarrow x(i)\Big(1 + \frac{b(i,j)}{B(i)}\Big) + \frac{b(i,j)}{(c-1)\,B(i)} .x(i)←x(i)(1+B(i)b(i,j)​)+(c−1)B(i)b(i,j)​.

Its revenue is the total amount charged.

Formalization targets

Goal: Theorem 1

For every instance and every bound R>0R > 0R>0 with b(i,j)≤R B(i)b(i,j) \le R\,B(i)b(i,j)≤RB(i) for all i,ji, ji,j, the Allocation Algorithm run with c=(1+R)1/Rc = (1+R)^{1/R}c=(1+R)1/R, under any tie-breaking of the maximum, satisfies for every feasible y′y'y′ of the packing LP

Revenue  ≥  (1−1c)(1−R)∑j∑ib(i,j) y′(i,j).\mathrm{Revenue} \;\ge\; \Big(1 - \frac1c\Big)(1 - R)\sum_{j}\sum_{i} b(i,j)\,y'(i,j).Revenue≥(1−c1​)(1−R)j∑​i∑​b(i,j)y′(i,j).

With R=Rmax⁡R = R_{\max}R=Rmax​ this is the paper's statement that the algorithm is (1−1/c)(1−Rmax⁡)(1 - 1/c)(1 - R_{\max})(1−1/c)(1−Rmax​)-competitive; the fractional optimum bounds every integral offline allocation.

Milestones

The proof of Theorem 1 rests on three claims and three auxiliary facts, each a milestone:

  1. the inequality ln⁡(1+x)/x≥ln⁡(1+y)/y\ln(1+x)/x \ge \ln(1+y)/yln(1+x)/x≥ln(1+y)/y for 0<x≤y≤10 < x \le y \le 10<x≤y≤1;
  2. Claim (1): the final (x,z)(x, z)(x,z) is feasible for the covering LP;
  3. Claim (2): the covering cost of the run equals (1+1/(c−1))(1 + 1/(c-1))(1+1/(c−1)) times the packing value of the run's own yyy;
  4. Inequality (1): x(i)≥1c−1(c∑jb(i,j)y(i,j)/B(i)−1)x(i) \ge \frac{1}{c-1}\big(c^{\sum_j b(i,j) y(i,j)/B(i)} - 1\big)x(i)≥c−11​(c∑j​b(i,j)y(i,j)/B(i)−1) at every stage of the run;
  5. Claim (3): ∑jb(i,j) y(i,j)≤B(i)+max⁡jb(i,j)\sum_j b(i,j)\,y(i,j) \le B(i) + \max_j b(i,j)∑j​b(i,j)y(i,j)≤B(i)+maxj​b(i,j), and the amount charged to iii is at least (1−R)∑jb(i,j) y(i,j)(1 - R)\sum_j b(i,j)\,y(i,j)(1−R)∑j​b(i,j)y(i,j);
  6. weak duality for the LP pair above;

and, separately, the second sentence of Theorem 1,

lim⁡R→0+(1−1(1+R)1/R)(1−R)=1−1e.\lim_{R\to 0^+}\Big(1 - \frac{1}{(1+R)^{1/R}}\Big)(1-R) = 1 - \frac1e .R→0+lim​(1−(1+R)1/R1​)(1−R)=1−e1​.

Significance

Theorem 1 gives an explicit ratio for every value of Rmax⁡R_{\max}Rmax​, not only in the limit. It tends to the optimal deterministic ratio 1−1/e1 - 1/e1−1/e as bids become small, and it quantifies how the guarantee degrades as single bids become a larger share of a budget. The primal-dual analysis is the template for the paper's later sections and for a line of work on online packing and covering problems, surveyed in Buchbinder and Naor's monograph The Design of Competitive Online Algorithms via a Primal-Dual Approach (Foundations and Trends in TCS, 2009).

The result is proved in the paper, and the proof is short. What this mission adds is a machine-checked proof about an algorithm that is defined, not described: the run is computed by recursion from the instance, and the guarantee is proved for that run and every tie-breaking. A related private mission on the platform, The Design of Competitive Online Algorithms via a Primal-Dual Approach VI: Maximizing Ad-Auctions Revenue, states the monograph's Theorem 10.1, which is this theorem, in a form that takes the analysis's intermediate inequalities as hypotheses over arbitrary lists of won bids; the present mission states it for the algorithm itself. No machine-checked proof of Theorem 1 is known to this mission.

Difficulty

Each step of the proof is elementary; the difficulty is the bookkeeping of an online process. Claims (1) and (2) are statements about a single iteration that must be lifted to the whole run: Claim (1) uses that xxx only increases, and Claim (2) that each product is processed once. Inequality (1) is an induction over the iterations that allocate to one buyer, interleaved with iterations that allocate to others and is the only place where the value of ccc matters. Claim (3) needs a further invariant: the amount charged equals the minimum of the allocated bids and the budget.

A tempting shortcut is to take Inequality (1) and the "at most one undercharge" fact as hypotheses about some list of bids. That does not describe the algorithm and is not the theorem; here the only hypotheses are on the instance and on the tie-breaking rule.

Formalization scope

Buyers are a type I with [Fintype I] and [Nonempty I]; products are Fin m, whose order is the arrival order. Bids and budgets are real, with B(i)>0B(i) > 0B(i)>0 and b(i,j)≥0b(i,j) \ge 0b(i,j)≥0. The state of the algorithm records xxx, the amounts charged, yyy and zzz; one iteration is step, the run after kkk products is runPrefix, and revenue sums the charges of the final state. The tie-breaking rule is a function sel of the current xxx and the product, required to return a maximizer of b(i,j)(1−x(i))b(i,j)(1-x(i))b(i,j)(1−x(i)); the theorem holds for every such rule. The constant c=(1+R)1/Rc = (1+R)^{1/R}c=(1+R)1/R is a real power and requires R>0R > 0R>0. The theorem is stated for any bound RRR on the ratios, of which the exact maximum is one instance. Claims (1) and (2) are stated for every c>1c > 1c>1, which covers the paper's choice. The paper's inequality for ln⁡(1+x)/x\ln(1+x)/xln(1+x)/x allows x=0x = 0x=0, read as a limit; the Lean statement requires x>0x > 0x>0.

A statement over an unconstrained allocation, or one conditioned on the proof's own intermediate inequalities, would be trivially true or false; the targets here concern only the run the definitions compute.

The development needs finite sums, real powers and logarithms from Mathlib and an induction principle for the run. The LP pair and weak duality are reusable for the paper's extensions, and the run invariants for any primal-dual online algorithm with multiplicative updates. Proofs of any milestone are welcome, as are sharper variants, such as the exact-Rmax⁡R_{\max}Rmax​ form or the bound against integral allocations.

Selected references

  • N. Buchbinder, K. Jain, J. Naor, Online Primal-Dual Algorithms for Maximizing Ad-Auctions Revenue, Algorithms – ESA 2007, LNCS 4698, 2007. https://doi.org/10.1007/978-3-540-75520-3_24
  • A. Mehta, A. Saberi, U. Vazirani, V. Vazirani, AdWords and Generalized Online Matching, Journal of the ACM 54(5), 2007. https://doi.org/10.1145/1284320.1284321
  • R. M. Karp, U. V. Vazirani, V. V. Vazirani, An Optimal Algorithm for On-line Bipartite Matching, STOC 1990. https://doi.org/10.1145/100216.100262
  • B. Kalyanasundaram, K. R. Pruhs, An Optimal Deterministic Algorithm for Online b-Matching, Theoretical Computer Science 233(1–2), 2000. https://doi.org/10.1016/S0304-3975(99)00140-1
  • N. Buchbinder, J. Naor, The Design of Competitive Online Algorithms via a Primal-Dual Approach, Foundations and Trends in Theoretical Computer Science 3(2–3), 2009. https://doi.org/10.1561/0400000024
10 thms3 active usersReviewed
🏆Completed
Linear OptimizationOperations ResearchOptimization+1·Captain: mikedeng1

Air Travel Demand and Airline Seat Inventory Management I: Marginal Seat Allocation Among Distinct Fare ClassesTextbook

Why airlines allocate seats by fare class

An airline sells the seats of one flight leg at several prices. Low fares fill seats that would otherwise fly empty; high fares are bought by passengers who book late and cannot be predicted exactly. Seat inventory control decides how many seats each fare class may sell. Peter Belobaba's 1987 MIT dissertation (MIT Flight Transportation Laboratory Report R87-7) gave the probabilistic treatment of this problem that became the expected marginal seat revenue (EMSR) method, which was in use across the airline industry for decades.

This mission formalizes the first, simplest model of the thesis: distinct (non-nested) fare-class inventories on a single leg, where a seat assigned to a class may be sold only in that class or not at all. The thesis surveys this model in Sect. 4.2 (pp. 84–94), where an integer-programming formulation from McDonnell-Douglas and its solution by ranking marginal values are described (p. 90), and develops it probabilistically in Sect. 5.1 (pp. 102–107).

Setting

A leg has capacity nnn seats (the thesis also writes CCC). There are finitely many fare classes iii. Class iii has an average fare fi≥0f_i \ge 0fi​≥0 and receives a random number of requests ri∈{0,1,2,… }r_i \in \{0,1,2,\dots\}ri​∈{0,1,2,…}, with law pip_ipi​. The seats are split into allocations Si∈NS_i \in \mathbb{N}Si​∈N, one per class.

With SSS seats, a class books requests until its seats run out, so its bookings and spill (refused requests) are

b=min⁡(r,S),l=(r−S)+(Eq. (5.3)).b = \min(r, S), \qquad l = (r - S)^+ \qquad \text{(Eq. (5.3))}.b=min(r,S),l=(r−S)+(Eq. (5.3)).

The expected revenue of class iii is Rˉi(Si)=fi⋅bˉi(Si)\bar R_i(S_i) = f_i \cdot \bar b_i(S_i)Rˉi​(Si​)=fi​⋅bˉi​(Si​) with bˉi(Si)=E[min⁡(ri,Si)]\bar b_i(S_i) = E[\min(r_i,S_i)]bˉi​(Si​)=E[min(ri​,Si​)], and the leg's expected revenue is Rˉ=∑iRˉi(Si)\bar R = \sum_i \bar R_i(S_i)Rˉ=∑i​Rˉi​(Si​) (Eq. (5.9)). Write

Pˉi(S)=P[ri≥S],EMSRi(S)=fi⋅Pˉi(S)(Eqs. (5.11), (6.1), (6.2)),\bar P_i(S) = P[r_i \ge S], \qquad \mathrm{EMSR}_i(S) = f_i \cdot \bar P_i(S) \qquad \text{(Eqs. (5.11), (6.1), (6.2))},Pˉi​(S)=P[ri​≥S],EMSRi​(S)=fi​⋅Pˉi​(S)(Eqs. (5.11), (6.1), (6.2)),

the expected marginal seat revenue of the SSS-th seat of class iii. The value of the kkk-th seat of class iii in the integer program of p. 90 is mi(k)=EMSRi(k)m_i(k) = \mathrm{EMSR}_i(k)mi​(k)=EMSRi​(k), k=1,…,nk = 1, \dots, nk=1,…,n.

Formalization targets

Goal: the nnn largest marginal values give the optimal booking limits (p. 90)

Let TTT be any set of nnn pairs (i,k)(i,k)(i,k), 1≤k≤n1 \le k \le n1≤k≤n, such that every mi(k)m_i(k)mi​(k) with (i,k)∈T(i,k) \in T(i,k)∈T is at least every mj(l)m_j(l)mj​(l) with (j,l)∉T(j,l) \notin T(j,l)∈/T, and let SiT=#{k:(i,k)∈T}S^T_i = \#\{k : (i,k) \in T\}SiT​=#{k:(i,k)∈T}. Then ∑iSiT=n\sum_i S^T_i = n∑i​SiT​=n and

∑ifi E[min⁡(ri,Si)]  ≤  ∑ifi E[min⁡(ri,SiT)]for every S with ∑iSi≤n.\sum_i f_i\, E[\min(r_i, S_i)] \;\le\; \sum_i f_i\, E[\min(r_i, S^T_i)] \qquad \text{for every } S \text{ with } \textstyle\sum_i S_i \le n .i∑​fi​E[min(ri​,Si​)]≤i∑​fi​E[min(ri​,SiT​)]for every S with ∑i​Si​≤n.

The goal fixes no distribution, number of classes or fare ordering, and it holds for every tie-breaking among equal marginal values.

Milestones

  1. Eq. (5.6): bˉi(S)+lˉi(S)=rˉi\bar b_i(S) + \bar l_i(S) = \bar r_ibˉi​(S)+lˉi​(S)=rˉi​ for requests of finite mean.
  2. Eq. (5.11): Rˉi(S)−Rˉi(S−1)=fi⋅P[ri≥S]\bar R_i(S) - \bar R_i(S-1) = f_i \cdot P[r_i \ge S]Rˉi​(S)−Rˉi​(S−1)=fi​⋅P[ri​≥S] for S≥1S \ge 1S≥1.
  3. Eqs. (6.1)–(6.2): Pˉi\bar P_iPˉi​ and EMSRi\mathrm{EMSR}_iEMSRi​ are non-increasing in SSS.
  4. Eq. (4.5): the 0–1 vector equal to 111 on a set of nnn largest mi(k)m_i(k)mi​(k) is an optimal solution of the linear program max⁡∑i,kXikmi(k)\max \sum_{i,k} X_{ik} m_i(k)max∑i,k​Xik​mi​(k) subject to ∑Xik≤n\sum X_{ik} \le n∑Xik​≤n, 0≤Xik≤10 \le X_{ik} \le 10≤Xik​≤1.
  5. Eq. (5.13), discrete form: an allocation of exactly CCC seats maximises Rˉ\bar RRˉ among such allocations if and only if some λ\lambdaλ satisfies EMSRi(Si)≥λ\mathrm{EMSR}_i(S_i) \ge \lambdaEMSRi​(Si​)≥λ whenever Si≥1S_i \ge 1Si​≥1 and EMSRi(Si+1)≤λ\mathrm{EMSR}_i(S_i + 1) \le \lambdaEMSRi​(Si​+1)≤λ, for all iii.

Significance

The goal is the reason distinct-inventory allocation is computationally easy: a revenue-maximising allocation is obtained by sorting n×(number of classes)n \times (\text{number of classes})n×(number of classes) numbers, with no search over allocations. The same marginal-value principle underlies the EMSR rules for nested classes in the rest of the thesis, and the identity (5.11) is the link between an expected-revenue function and its marginal seat values used throughout revenue management. Milestone 5 is the integer form of the Lagrangian condition of Eq. (5.13); the thesis states it only for a continuous relaxation, and its equality form is generally unattainable with integer seats.

These results are classical and their proofs are elementary; to our knowledge none of them has been machine-checked. The mission produces a reusable formal model of a single-leg, distinct-inventory allocation problem with integer demand (bookings, spill, revenue and marginal values of a PMF ℕ), and checked statements of the marginal-allocation principle for it.

Difficulty

The thesis argues with continuous densities and derivatives, setting ∂Rˉ/∂Si\partial \bar R / \partial S_i∂Rˉ/∂Si​ equal across classes. That argument does not transfer to integer seats: the derivative of a step-shaped expected-revenue function does not exist, equality of marginal values across classes generally fails at every integer allocation, and the tail probability must be P[r≥S]P[r \ge S]P[r≥S] rather than the P[r>S]P[r > S]P[r>S] of Eq. (5.2) for the marginal identity to hold. The integer statements need their own exchange argument. A second subtlety is ties: "the nnn largest values" is not unique, and the goal must hold for every admissible choice, including choices in which a class's selected seat numbers are not an initial segment {1,…,Si}\{1, \dots, S_i\}{1,…,Si​}.

Formalization scope

All declarations live in the namespace SeatInventory.Distinct. Conventions:

  • Integer demand. The law of class iii's requests is d i : PMF ℕ; expectations are series over N\mathbb{N}N. The thesis's continuous densities are replaced by this discrete model, which the thesis itself requires for seat allocations (p. 103).
  • Tail convention. Pˉ(S)=P[r≥S]\bar P(S) = P[r \ge S]Pˉ(S)=P[r≥S], as in Eq. (6.2) and the prose of Eq. (5.11) ("the probability of selling SiS_iSi​ or more seats"), not the P[r>S]P[r > S]P[r>S] of Eq. (5.2).
  • Seat numbers start at 1, and the pairs (i,k)(i,k)(i,k) range over k∈{1,…,n}k \in \{1, \dots, n\}k∈{1,…,n}, as the 600 variables of a 150-seat, four-class problem on p. 90 indicate.
  • Nonnegative fares fi≥0f_i \ge 0fi​≥0 are assumed in every statement that needs them; with a negative fare the capacity constraint ∑Si≤n\sum S_i \le n∑Si​≤n would not bind and the claims fail.
  • Finite mean of the requests is assumed explicitly for Eq. (5.6); the thesis assumes it silently. Expected bookings are bounded and need no assumption.
  • Capacity. The goal and the LP compare against allocations with ∑iSi≤n\sum_i S_i \le n∑i​Si​≤n (the LP's constraint); milestone 5 compares allocations of exactly CCC seats (Eq. (5.8)).
  • No independence assumption. Expected revenue of distinct inventories depends only on each class's marginal law, so the statements take one law per class.
  • "Decreasing" is non-increasing. The thesis's justification in Sect. 6.1.1 gives only monotonicity; strict decrease fails for bounded demand.
  • LP integrality. "The solution will be integer" is stated as: the indicator of every set of nnn largest values is optimal. With ties the LP also has fractional optima.

The expected revenue in the goal is computed from the booking rule min⁡(ri,Si)\min(r_i, S_i)min(ri​,Si​); it is not defined as a sum of marginal values, and the optimal allocation is not defined as an argmax of Rˉ\bar RRˉ. Either shortcut would make the goal a tautology and is ruled out.

Needed infrastructure: tail sums of a PMF ℕ, telescoping of E[min⁡(r,S)]E[\min(r, S)]E[min(r,S)], and a finite exchange argument for sums of the nnn largest values of a function on a finite set; the last two are reusable for any separable concave resource-allocation problem. Contributions of proofs of any milestone, and of a verified sorting routine that produces a set of nnn largest values, are welcome.

Selected references

  • P. P. Belobaba, Air Travel Demand and Airline Seat Inventory Management, PhD thesis, MIT Flight Transportation Laboratory Report R87-7, 1987 (no DOI).
  • P. P. Belobaba, Airline yield management: an overview of seat inventory control, Transportation Science 21(2), 63–73, 1987. https://doi.org/10.1287/trsc.21.2.63
  • P. P. Belobaba, Application of a probabilistic decision model to airline seat inventory control, Operations Research 37(2), 183–197, 1989. https://doi.org/10.1287/opre.37.2.183
  • K. Littlewood, Forecasting and control of passenger bookings, AGIFORS Symposium Proceedings 12, 1972; reprinted in Journal of Revenue and Pricing Management 4(2), 111–123, 2005. https://doi.org/10.1057/palgrave.rpm.5170134
8 thms4 active usersReviewed
🏆Completed
Operations ResearchOptimizationProbability·Captain: mikedeng1

Air Travel Demand and Airline Seat Inventory Management II: The EMSR Protection Level for Two Nested Fare ClassesTextbook

Why airlines protect seats

An airline sells the seats of one flight in several fare classes at different prices, all drawn from one shared cabin. Discount fares are bought early, under advance-purchase restrictions, while most high-fare requests arrive close to departure. Accepting every early low-fare request fills the aircraft with cheap passengers and turns away late high-fare passengers; refusing too many leaves seats empty. Seat inventory control decides how many seats to keep away from the low fare. Peter Belobaba's 1987 MIT thesis (Flight Transportation Laboratory Report R87-7) introduced the expected marginal seat revenue (EMSR) model for this decision, and EMSR-type rules remain the basis of the booking-limit logic in airline revenue management systems.

Timeline. Littlewood (1972, AGIFORS Symposium Proceedings; reprinted 2005) proposed accepting a low-fare request as long as its fare is at least the high fare times the probability of selling all remaining seats to high-fare passengers. Analysts at Trans World Airlines (1973) and Richter at Lufthansa (1982) gave equivalent formulations for the dynamic case. Belobaba (1987, Ch. 5) restated the two-class rule as a static protection level for nested inventories and extended it heuristically to many classes. Brumelle and McGill (Operations Research 41, 1993) and Curry (Transportation Science 24, 1990) later proved optimality of nested protection levels for any number of classes under low-to-high arrivals, and showed that Belobaba's multi-class EMSR levels are not optimal for three or more classes. This mission concerns only the two-class result, which is correct.

Setting

A single flight leg has capacity C∈NC \in \mathbb NC∈N. Class 1 has fare f1f_1f1​, class 2 has fare f2f_2f2​, with 0≤f2≤f10 \le f_2 \le f_10≤f2​≤f1​. The numbers of requests for the two classes are random variables r1,r2r_1, r_2r1​,r2​ with values in N\mathbb NN, defined on a probability space (Ω,μ)(\Omega, \mu)(Ω,μ) and independent. There are no cancellations, no no-shows, and a refused request is lost.

The inventory is nested: a class-1 request is accepted as long as any seat is unsold. A protection level S∈{0,…,C}S \in \{0, \dots, C\}S∈{0,…,C} is the number of seats reserved for class 1; it sets the class-2 booking limit BL2=C−SBL_2 = C - SBL2​=C−S. All class-2 requests arrive before any class-1 request. Class 2 therefore books min⁡(r2,C−S)\min(r_2, C - S)min(r2​,C−S) seats and class 1 books min⁡(r1,C−min⁡(r2,C−S))\min(r_1, C - \min(r_2, C-S))min(r1​,C−min(r2​,C−S)), and the realised revenue is

RS=f2min⁡(r2,C−S)+f1min⁡(r1, C−min⁡(r2,C−S)).R_S = f_2 \min(r_2, C - S) + f_1 \min\bigl(r_1,\, C - \min(r_2, C - S)\bigr).RS​=f2​min(r2​,C−S)+f1​min(r1​,C−min(r2​,C−S)).

The expected revenue is Rˉ(S)=E[RS]\bar R(S) = \mathbb E[R_S]Rˉ(S)=E[RS​].

The tail probability of class 1 is Pˉ1(S)=P[r1≥S]\bar P_1(S) = P[r_1 \ge S]Pˉ1​(S)=P[r1​≥S], the probability of receiving SSS or more class-1 requests, and the expected marginal seat revenue of the SSS-th class-1 seat is

EMSR1(S)=f1⋅Pˉ1(S).\mathrm{EMSR}_1(S) = f_1 \cdot \bar P_1(S).EMSR1​(S)=f1​⋅Pˉ1​(S).

For a single class with SSS seats the expected revenue is f1 E[min⁡(r1,S)]f_1\,\mathbb E[\min(r_1, S)]f1​E[min(r1​,S)], and EMSR1(S)\mathrm{EMSR}_1(S)EMSR1​(S) is its increment from S−1S-1S−1 to SSS seats. The EMSR protection level S21S_2^1S21​ is the largest integer S∈{0,…,C}S \in \{0, \dots, C\}S∈{0,…,C} with

EMSR1(S)≥f2.\mathrm{EMSR}_1(S) \ge f_2 .EMSR1​(S)≥f2​.

In Lean these objects are nestedRevenue, expectedNestedRevenue, tailProb, classRevenue, emsr and emsrProtectionLevel in SeatInventory.Nested.

Formalization targets

Goal: Eqs. (5.15)–(5.16), optimality of the EMSR protection level

Rˉ(S)≤Rˉ(S21)for all S∈{0,…,C}.\bar R(S) \le \bar R(S_2^1) \qquad \text{for all } S \in \{0, \dots, C\}.Rˉ(S)≤Rˉ(S21​)for all S∈{0,…,C}.

The goal fixes no distribution: it holds for every pair of independent N\mathbb NN-valued demands, and S21S_2^1S21​ depends only on f2/f1f_2/f_1f2​/f1​ and the law of r1r_1r1​.

Milestones

  1. Eq. (5.11). f1E[min⁡(r1,S)]−f1E[min⁡(r1,S−1)]=f1P[r1≥S]f_1\mathbb E[\min(r_1,S)] - f_1\mathbb E[\min(r_1,S-1)] = f_1 P[r_1 \ge S]f1​E[min(r1​,S)]−f1​E[min(r1​,S−1)]=f1​P[r1​≥S] for S≥1S \ge 1S≥1.
  2. Eqs. (6.1)–(6.2). Pˉ1\bar P_1Pˉ1​ and, for f1≥0f_1 \ge 0f1​≥0, EMSR1\mathrm{EMSR}_1EMSR1​ are non-increasing in SSS.
  3. Eq. (4.8), Littlewood's rule, already on the platform as RevenueManagement.littlewood_marginal_value (Talluri and van Ryzin's Eq. (2.1), proved).
  4. Sect. 5.2, p. 112. Rˉ(S)≤Rˉ(S21)\bar R(S) \le \bar R(S_2^1)Rˉ(S)≤Rˉ(S21​) for S21≤S≤CS_2^1 \le S \le CS21​≤S≤C: a smaller booking limit for class 2 cannot raise expected revenue.
  5. Sect. 5.2, p. 114. With the same class-2 limit C−SC - SC−S, the expected nested revenue is at least the expected revenue of two distinct inventories with SSS and C−SC - SC−S seats, strictly if f1>0f_1 > 0f1​>0 and P[r2<C−S, r1>S]>0P[r_2 < C - S,\ r_1 > S] > 0P[r2​<C−S, r1​>S]>0.

Significance

The two-class result says that, for a static booking limit set once before sales open and low-fare demand arriving first, the airline needs only the high-fare demand distribution and the fare ratio to set the optimal limit; the low-fare forecast is irrelevant. This is the rule that the thesis then applies class by class in multi-class nested systems, and it is the base case against which the later exact multi-class theory (Brumelle–McGill, Curry) is checked. Milestone 5 makes precise why nested inventories dominate the distinct-inventory allocation of the thesis's Sect. 5.1 with the same class-2 limit.

The result is classical and proved, in the sense that the optimality of a two-class threshold policy follows from Littlewood's argument and from the dynamic-programming treatment in Talluri and van Ryzin's The Theory and Practice of Revenue Management (2004, Ch. 2). On Prove2Me, Littlewood's marginal rule and the dynamic-programming optimality of nested protection levels (RevenueManagement.static_optimal_controls) are formalized, but in Bellman form: there the protection level is defined through the value function of a dynamic program. What is not formalized is the statement in Belobaba's form, where the protection level is the explicit threshold of f1P[r1≥S]f_1 P[r_1 \ge S]f1​P[r1​≥S] against f2f_2f2​ and the objective is the explicit expected revenue of a booking limit. Connecting the two forms, and the comparison with distinct inventories, is the work of this mission.

Difficulty

The expected revenue couples the two demands through the capacity left by class 2, so Rˉ\bar RRˉ is not a sum of single-class revenues and is not separately concave in an obvious way. The step that requires care is the increment Rˉ(S)−Rˉ(S−1)\bar R(S) - \bar R(S-1)Rˉ(S)−Rˉ(S−1): it is not EMSR1(S)−f2\mathrm{EMSR}_1(S) - f_2EMSR1​(S)−f2​, as the thesis's sentence after the milestone on p. 112 suggests, because the extra protected seat matters only on the event that class 2 would have reached its limit. Independence of r1r_1r1​ and r2r_2r2​ is what makes that event's probability factor out; without independence the threshold rule is not optimal. The discrete reading matters too: with P[r1>S]P[r_1 > S]P[r1​>S] in place of P[r1≥S]P[r_1 \ge S]P[r1​≥S] the rule is off by one seat and the claim fails.

Formalization scope

Conventions the Lean statements commit to:

  • Demands are N\mathbb NN-valued measurable random variables r₁ r₂ : Ω → ℕ on a probability space μ; the goal and milestone 4 assume IndepFun r₁ r₂ μ. The thesis writes continuous densities (Eqs. (5.1)–(5.5)) but requires integer seat counts; the discrete model is used throughout.
  • Pˉ1(S)=P[r1≥S]\bar P_1(S) = P[r_1 \ge S]Pˉ1​(S)=P[r1​≥S], as in Eq. (6.2) and the prose of Eq. (5.11), not P[r1>S]P[r_1 > S]P[r1​>S] as in Eq. (5.2).
  • The EMSR protection level is the largest S∈{0,…,C}S \in \{0,\dots,C\}S∈{0,…,C} with f1P[r1≥S]≥f2f_1 P[r_1 \ge S] \ge f_2f1​P[r1​≥S]≥f2​ (Eq. (5.15)); Eq. (5.16)'s equality is the continuous idealisation and is not stated.
  • Booking order: all class-2 requests precede all class-1 requests (pp. 108, 112). This order is built into the revenue formula, not assumed separately.
  • Fares satisfy 0≤f2≤f10 \le f_2 \le f_10≤f2​≤f1​; the thesis has f1>f2f_1 > f_2f1​>f2​, and the statements also cover equality.
  • Expectations are Bochner integrals of bounded revenues, probabilities are μ.real; seat counts use truncated subtraction only where S≤CS \le CS≤C.

A trivializing formalization is ruled out: S21S_2^1S21​ is defined by the threshold of (5.15), never as an argmax of expected revenue, and the expected revenue is computed from the realised revenue of the booking process, not postulated as a sum of marginal terms.

The multi-class EMSR levels of Eqs. (5.19)–(5.29) and the dynamic revision of Eqs. (5.31)–(5.32) are out of scope. Proofs need the discrete expectation identity E[min⁡(r,S)]−E[min⁡(r,S−1)]=P[r≥S]\mathbb E[\min(r,S)] - \mathbb E[\min(r,S-1)] = P[r \ge S]E[min(r,S)]−E[min(r,S−1)]=P[r≥S] and expectation of products of independent bounded functions, both in Mathlib's reach and reusable for other single-leg revenue models. Proofs of any milestone, and a proof of the goal from milestones 1, 2 and 4 plus the matching lower-half argument, are welcome.

Selected references

  • P. P. Belobaba, Air Travel Demand and Airline Seat Inventory Management, PhD thesis, MIT, Flight Transportation Laboratory Report R87-7, 1987 (no DOI).
  • K. Littlewood, Forecasting and control of passenger bookings, AGIFORS Symposium Proceedings 12, 1972; reprinted in Journal of Revenue and Pricing Management 4, 2005. https://doi.org/10.1057/palgrave.rpm.5170134
  • S. L. Brumelle and J. I. McGill, Airline seat allocation with multiple nested fare classes, Operations Research 41, 1993. https://doi.org/10.1287/opre.41.1.127
  • R. E. Curry, Optimal airline seat allocation with fare classes nested by origins and destinations, Transportation Science 24, 1990. https://doi.org/10.1287/trsc.24.3.193
  • K. T. Talluri and G. J. van Ryzin, The Theory and Practice of Revenue Management, Springer, 2004. https://doi.org/10.1007/b139000
8 thms4 active usersReviewed
🏆Completed
Operations ResearchOptimizationProbability·Captain: mikedeng1

Air Travel Demand and Airline Seat Inventory Management III: Gaussian EMSR Protection Levels and Their SensitivityTextbook

Why protection levels and their inputs matter

An airline sells the seats of one flight leg in several fare classes at different prices. Low-fare passengers usually book first, so the airline must decide how many seats to keep back, or protect, for later high-fare passengers. Peter Belobaba's 1987 MIT dissertation introduced the expected marginal seat revenue (EMSR) rule for this decision, and EMSR-type rules became a standard of airline revenue management practice (Talluri and van Ryzin 2004). A protection level is computed from a demand forecast, and forecasts are uncertain. Section 6.2 of the dissertation asks how the protection level moves when its inputs move: the mean of forecast demand, its standard deviation, and the ratio of the two fares. That question decides where forecasting effort pays off, and this mission formalizes the answers the dissertation gives for Gaussian demand.

This is the third mission in a series on the dissertation. The first treats marginal allocation among distinct fare classes, and the second the two-class nested protection level in the discrete model, including its revenue optimality. This mission takes the continuous Gaussian model of Chapter 6 on its own terms.

Setting

Let rrr be the number of requests for a fare class, a real random variable with law μ\muμ. For a seat level S∈RS \in \mathbb RS∈R the tail probability is

Pˉ(S)=P[r≥S],\bar P(S) = P[r \ge S],Pˉ(S)=P[r≥S],

and for the fare fff of the class the expected marginal seat revenue is EMSR(S)=Pˉ(S)⋅f\mathrm{EMSR}(S) = \bar P(S)\cdot fEMSR(S)=Pˉ(S)⋅f (Eqs. (6.1)–(6.2)).

There are two classes: class 1 with fare f1f_1f1​ and class 2 with fare f2f_2f2​, where 0<f2<f10 < f_2 < f_10<f2​<f1​. Requests for class 1 are Gaussian with estimated mean rˉ\bar rrˉ and estimated standard deviation σ^>0\hat\sigma > 0σ^>0, written r1∼N(rˉ,σ^2)r_1 \sim N(\bar r, \hat\sigma^2)r1​∼N(rˉ,σ^2). A real number SSS is an EMSR protection level for class 1 against class 2 when

Pˉ1(S)=P[r1≥S]=f2f1(Eq. (6.10)).\bar P_1(S) = P[r_1 \ge S] = \frac{f_2}{f_1} \qquad \text{(Eq. (6.10))}.Pˉ1​(S)=P[r1​≥S]=f1​f2​​(Eq. (6.10)).

The standardized level ZZZ is the value "which has a probability of f2/f1f_2/f_1f2​/f1​ of being exceeded" by a standard normal variable:

P[N(0,1)≥Z]=f2f1.P[N(0,1) \ge Z] = \frac{f_2}{f_1}.P[N(0,1)≥Z]=f1​f2​​.

In the Lean development these are tailProb, emsr, gaussianLaw rbar σ, stdNormal, IsProtectionLevel rbar σ f₁ f₂ S and IsStdNormalLevel f₁ f₂ Z, all in the namespace SeatInventory.Gaussian.

Formalization targets

Goal: the Gaussian protection level and its sensitivity to σ^\hat\sigmaσ^

For σ^>0\hat\sigma > 0σ^>0 and 0<f2<f10 < f_2 < f_10<f2​<f1​:

  1. Eq. (6.10) has exactly one solution SSS, and the standard normal equation has exactly one solution ZZZ;
  2. they satisfy
S=rˉ+Zσ^(Eq. (6.12));S = \bar r + Z\hat\sigma \qquad \text{(Eq. (6.12))};S=rˉ+Zσ^(Eq. (6.12));
  1. Z<0Z < 0Z<0 if f2/f1>1/2f_2/f_1 > 1/2f2​/f1​>1/2, Z>0Z > 0Z>0 if f2/f1<1/2f_2/f_1 < 1/2f2​/f1​<1/2, Z=0Z = 0Z=0 if f2/f1=1/2f_2/f_1 = 1/2f2​/f1​=1/2 (Eq. (6.14)), and S=rˉS = \bar rS=rˉ in the last case;
  2. if σ^′>σ^\hat\sigma' > \hat\sigmaσ^′>σ^ and S′S'S′ solves (6.10) for N(rˉ,σ^′2)N(\bar r, \hat\sigma'^2)N(rˉ,σ^′2), then S′<SS' < SS′<S, S′>SS' > SS′>S or S′=SS' = SS′=S according as f2/f1f_2/f_1f2​/f1​ is above, below or equal to 1/21/21/2.

The goal states no numerical constant and no particular fare ratio; it fixes only the shape of the dependence.

Milestones, in attack order

  • Eq. (6.1)–(6.2): for any request law, Pˉ\bar PPˉ and EMSR\mathrm{EMSR}EMSR are non-increasing in SSS.
  • Eq. (6.10): the Gaussian protection level exists and is unique.
  • Eq. (6.11)–(6.12): S=rˉ+Zσ^S = \bar r + Z\hat\sigmaS=rˉ+Zσ^.
  • p. 154: with σ^\hat\sigmaσ^ and the fares fixed, replacing rˉ\bar rrˉ by rˉ+c\bar r + crˉ+c replaces SSS by S+cS + cS+c.
  • Eq. (6.14): the sign of ZZZ, and S=rˉS = \bar rS=rˉ at fare ratio 1/21/21/2 for every σ^\hat\sigmaσ^.
  • p. 154: the effect of σ^\hat\sigmaσ^ on SSS (part 4 of the goal on its own).
  • p. 157: ZZZ and SSS decrease strictly as the fare ratio f2/f1f_2/f_1f2​/f1​ increases.

The dissertation's constant-coefficient-of-variation form, Eq. (6.13), S=rˉ(1+Zk)S = \bar r(1 + Zk)S=rˉ(1+Zk) with k=σ^/rˉk = \hat\sigma/\bar rk=σ^/rˉ, follows from (6.12) by substitution and is not stated separately.

Significance

The result gives every Gaussian protection level as a closed form in one standard normal quantile. From it come the three sensitivities that Sect. 6.2 uses to argue for better forecasts. The protection level moves one-for-one with mean demand. The standard deviation moves it in a direction fixed only by whether the discount fare is above or below half the full fare. A higher fare ratio always lowers it. The dissertation uses these facts, and its Figures 6.1 and 6.2, to argue that reducing the estimated standard deviation of demand narrows the range of protection levels a forecast can produce. The same quantile structure is behind Littlewood's rule and the newsvendor critical fractile, so the statements here are the Gaussian specialization of a pattern that recurs throughout revenue management and inventory theory.

All the statements are classical and easy to believe. None of them, to our knowledge, has a machine-checked proof. Mathlib provides the Gaussian law and its affine images, but no standard normal quantile and no statement that a Gaussian tail is a strictly decreasing bijection onto (0,1)(0,1)(0,1). Formalizing this mission produces both, in a form that can be used again wherever a normal critical fractile appears.

Difficulty

Most of the work is in the existence and uniqueness of the two tail solutions. The tail S↦P[r1≥S]S \mapsto P[r_1 \ge S]S↦P[r1​≥S] must be shown continuous, strictly decreasing, and to take every value in (0,1)(0,1)(0,1). Strictness needs the Gaussian density to be positive everywhere, and existence needs a limit argument at both ends. Monotonicity alone, which holds for every law (Eqs. (6.1)–(6.2)), gives neither, because a general law can have flat stretches and jumps in its tail. The relation S=rˉ+Zσ^S = \bar r + Z\hat\sigmaS=rˉ+Zσ^ then requires transporting the tail of N(rˉ,σ^2)N(\bar r,\hat\sigma^2)N(rˉ,σ^2) to that of N(0,1)N(0,1)N(0,1) through the affine map x↦(x−rˉ)/σ^x \mapsto (x - \bar r)/\hat\sigmax↦(x−rˉ)/σ^, and the sign of ZZZ requires the symmetry of N(0,1)N(0,1)N(0,1), namely P[N(0,1)≥0]=1/2P[N(0,1) \ge 0] = 1/2P[N(0,1)≥0]=1/2. Once uniqueness is available, each sensitivity statement follows from these facts. The tempting shortcut of reading S=rˉ+Zσ^S = \bar r + Z\hat\sigmaS=rˉ+Zσ^ as a definition is ruled out below.

Formalization scope

  • Continuous seats. Protection levels and ZZZ are real numbers, as in the dissertation's own Gaussian example (Z=−0.675Z = -0.675Z=−0.675 at fare ratio 0.750.750.75). This differs from the first two missions of the series, which count seats in N\mathbb NN. For a continuous law P[r≥S]=P[r>S]P[r \ge S] = P[r > S]P[r≥S]=P[r>S], so the two definitions of Pˉ\bar PPˉ the dissertation uses (Eq. (5.2) and Eq. (6.2)) coincide here.
  • Gaussian law. N(rˉ,σ^2)N(\bar r, \hat\sigma^2)N(rˉ,σ^2) is Mathlib's gaussianReal rbar (σ^2), parameterised by the variance. Every theorem assumes σ^>0\hat\sigma > 0σ^>0; at σ^=0\hat\sigma = 0σ^=0 the law is a Dirac mass and (6.10) has no solution.
  • Fares. 0<f2<f10 < f_2 < f_10<f2​<f1​, so f2/f1∈(0,1)f_2/f_1 \in (0,1)f2​/f1​∈(0,1). This is the dissertation's "f2<f1f_2 < f_1f2​<f1​" together with positive fares.
  • Relational sensitivity. The sensitivity statements compare any two solutions of (6.10) under the two input values. Together with uniqueness, this is the same as monotonicity of the solution map. No function is defined by a choice operator.
  • Tail as a real number. Pˉ(S)\bar P(S)Pˉ(S) is the measure of [S,∞)[S,\infty)[S,∞) as a real number. The law is a probability measure, so nothing is truncated.
  • No trivialization. SSS is defined only by the tail equation (6.10) for N(rˉ,σ^2)N(\bar r, \hat\sigma^2)N(rˉ,σ^2), and ZZZ only by the tail equation for N(0,1)N(0,1)N(0,1). Neither is defined by the formula S=rˉ+Zσ^S = \bar r + Z\hat\sigmaS=rˉ+Zσ^, which would make Eq. (6.12) true by definition.
  • Not covered. The revenue optimality of the level defined by (6.10) belongs to the second mission. The multi-class EMSR rules (5.19)–(5.29) are not optimal for three or more classes and are not stated. The empirical analysis of Sect. 6.1 is out of scope.

Useful infrastructure, all reusable: the strict monotonicity, continuity and range of Gaussian tails; the standard normal quantile; and tail transport under affine maps. Contributions of these as separate lemmas are welcome.

Selected references

  • P. P. Belobaba, Air Travel Demand and Airline Seat Inventory Management, PhD thesis, MIT Flight Transportation Laboratory Report R87-7, 1987. (no DOI; the source PDF of this mission).
  • P. P. Belobaba, Application of a probabilistic decision model to airline seat inventory control, Operations Research 37(2):183–197, 1989. https://doi.org/10.1287/opre.37.2.183
  • K. Littlewood, Forecasting and control of passenger bookings, AGIFORS Symposium Proceedings 12, 1972; reprinted in Journal of Revenue and Pricing Management 4:111–123, 2005. https://doi.org/10.1057/palgrave.rpm.5170134
  • K. T. Talluri and G. J. van Ryzin, The Theory and Practice of Revenue Management, Springer, 2004. https://doi.org/10.1007/b139000
9 thms2 active usersReviewed
🏆Completed
Operations ResearchProbability·Captain: mikedeng1

Conditional Logit Analysis of Qualitative Choice Behavior 1: Independence of Irrelevant Alternatives with a Universal Benchmark Yields Logit Selection ProbabilitiesResearch Paper

Motivation

The conditional logit model is the workhorse of discrete choice analysis: it is used to forecast travel mode shares, to estimate demand for differentiated products, and, in operations research, as the multinomial logit (MNL) choice model behind assortment optimization and revenue management. Its selection probabilities have the form P(x∣s,B)=ev(s,x)/∑y∈Bev(s,y)P(x\mid s,B) = e^{v(s,x)}/\sum_{y\in B} e^{v(s,y)}P(x∣s,B)=ev(s,x)/∑y∈B​ev(s,y). Daniel McFadden's 1974 chapter Conditional Logit Analysis of Qualitative Choice Behavior gave the model two behavioural foundations, one of which is the subject of this mission: the logit form is a consequence of a single axiom on how choice probabilities change when the set of available alternatives changes.

That axiom is Luce's choice axiom, which McFadden calls Independence of Irrelevant Alternatives (IIA): the relative odds of choosing one alternative over another do not depend on which other alternatives are present. Luce (1959) introduced it; McFadden (1974, §I) showed how, together with positivity and a mild condition on which alternative sets can occur, it yields the conditional logit form with a "utility indicator" v(s,x)v(s,x)v(s,x) shared by all alternative sets.

Timeline. Luce, Individual Choice Behavior (1959): the choice axiom and its ratio-scale representation. McFadden (1974, pp. 109–110): the derivation in the econometric setting with measured attributes sss, the binary-odds identities (5)–(10), and footnote 3, which removes an extra axiom (Axiom 3) by a universal benchmark alternative. McFadden (1974, pp. 111–112): the companion random-utility characterization by extreme-value shocks, treated in mission 2 of this series.

Setting

Let XXX be the universe of objects of choice and SSS the universe of vectors of measured attributes of decision-makers. An alternative set is a finite set B⊆XB\subseteq XB⊆X; a designated family of finite sets is the family of possible alternative sets. The selection probability P(x∣s,B)P(x\mid s,B)P(x∣s,B) is the probability that an individual drawn at random from the population, with attributes sss and facing BBB, chooses x∈Bx\in Bx∈B. For every sss and possible BBB, x↦P(x∣s,B)x\mapsto P(x\mid s,B)x↦P(x∣s,B) is a probability vector on BBB. Whenever x≠yx\neq yx=y belong to a possible set, the pair {x,y}\{x,y\}{x,y} is possible too, so binary choices are defined.

  • Axiom 1 (IIA). For all possible BBB, all sss and all x,y∈Bx,y\in Bx,y∈B: P(x∣s,{x,y})P(y∣s,B)=P(y∣s,{x,y})P(x∣s,B)P(x\mid s,\{x,y\})P(y\mid s,B) = P(y\mid s,\{x,y\})P(x\mid s,B)P(x∣s,{x,y})P(y∣s,B)=P(y∣s,{x,y})P(x∣s,B).
  • Axiom 2 (Positivity). P(x∣s,B)>0P(x\mid s,B)>0P(x∣s,B)>0 for all possible BBB, all sss, all x∈Bx\in Bx∈B.
  • Binary probabilities. pxy=P(x∣s,{x,y})p_{xy}=P(x\mid s,\{x,y\})pxy​=P(x∣s,{x,y}) for x≠yx\neq yx=y, and pxx=12p_{xx}=\tfrac12pxx​=21​ by definition.
  • The function VVV. V(s,x,z)=log⁡(pxz/pzx)V(s,x,z)=\log(p_{xz}/p_{zx})V(s,x,z)=log(pxz​/pzx​).
  • Universal benchmark. An alternative zzz such that B∪{z}B\cup\{z\}B∪{z} is possible whenever BBB is.

In Lean these are IsSelectionProb, PairsPossible, Axiom1, Axiom2, binProb, altSetV and IsUniversalBenchmark in the namespace McFadden1974.IIA.

Formalization targets

Goal: footnote 3 with Equation (12)

Under Axioms 1 and 2 and a universal benchmark zzz, with v(s,x)=V(s,x,z)v(s,x)=V(s,x,z)v(s,x)=V(s,x,z), for every sss, every possible BBB (containing zzz or not) and every x∈Bx\in Bx∈B:

P(x∣s,B)=ev(s,x)∑y∈Bev(s,y).P(x\mid s,B) = \frac{e^{v(s,x)}}{\sum_{y\in B} e^{v(s,y)}}.P(x∣s,B)=∑y∈B​ev(s,y)ev(s,x)​.

The function vvv is the same for all alternative sets; this is what distinguishes the goal from Equation (10).

Milestones, in the paper's order

  1. Equation (5): for x≠yx\neq yx=y in BBB with P(x∣s,B)>0P(x\mid s,B)>0P(x∣s,B)>0, Axiom 1 gives P(x∣s,{x,y})>0P(x\mid s,\{x,y\})>0P(x∣s,{x,y})>0 and P(y∣s,{x,y})P(x∣s,{x,y})=P(y∣s,B)P(x∣s,B)\dfrac{P(y\mid s,\{x,y\})}{P(x\mid s,\{x,y\})}=\dfrac{P(y\mid s,B)}{P(x\mid s,B)}P(x∣s,{x,y})P(y∣s,{x,y})​=P(x∣s,B)P(y∣s,B)​.
  2. Equations (6)–(7): P(y∣s,B)=pyxpxyP(x∣s,B)P(y\mid s,B)=\dfrac{p_{yx}}{p_{xy}}P(x\mid s,B)P(y∣s,B)=pxy​pyx​​P(x∣s,B) and 1=(∑y∈Bpyxpxy)P(x∣s,B)1=\Big(\sum_{y\in B}\dfrac{p_{yx}}{p_{xy}}\Big)P(x\mid s,B)1=(∑y∈B​pxy​pyx​​)P(x∣s,B).
  3. Equation (8): P(x∣s,B)=1/∑y∈B(pyx/pxy)P(x\mid s,B)=1\big/\sum_{y\in B}(p_{yx}/p_{xy})P(x∣s,B)=1/∑y∈B​(pyx​/pxy​).
  4. Equation (9): pyxpxy=pyz/pzypxz/pzx\dfrac{p_{yx}}{p_{xy}}=\dfrac{p_{yz}/p_{zy}}{p_{xz}/p_{zx}}pxy​pyx​​=pxz​/pzx​pyz​/pzy​​ for x,y,zx,y,zx,y,z in a possible set.
  5. Equation (10): for a benchmark z∈Bz\in Bz∈B, P(x∣s,B)=eV(s,x,z)/∑y∈BeV(s,y,z)P(x\mid s,B)=e^{V(s,x,z)}\big/\sum_{y\in B}e^{V(s,y,z)}P(x∣s,B)=eV(s,x,z)/∑y∈B​eV(s,y,z).

Significance

The result. The goal identifies a testable axiom on choice probabilities, IIA, with a parametric functional form, the conditional logit model. It is what licenses the econometric specification v(s,x)=θ′z(s,x)v(s,x)=\theta'z(s,x)v(s,x)=θ′z(s,x) estimated in the rest of McFadden's chapter, and it is the reason the MNL model is the default in assortment and pricing problems in operations research. It also makes the model's limitations precise: any population whose choices violate IIA (the auto/red-bus/blue-bus example on p. 113 of the chapter) cannot be logit.

Formalizing it. The result is classical and proved on paper. No machine-checked statement of it exists on the platform, which has the logit form only as a definition (soft-max, MNL revenue) and IIA only in Arrow's social-choice sense, a different axiom about preference aggregation. This mission produces a formal statement of the derivation with every standing assumption explicit, including two the paper leaves implicit: that selection probabilities are normalized on binary sets, and that binary subsets of possible sets are possible.

Difficulty

The algebra is elementary; the difficulty is bookkeeping of where each axiom may be applied. Axioms 1 and 2 are assumed only on possible alternative sets. Equation (10) needs the benchmark to lie in the alternative set, and the naive argument "pick z∈Bz\in Bz∈B as benchmark" produces a function V(s,x,z)V(s,x,z)V(s,x,z) that depends on the set through the choice of zzz. The goal requires a single vvv for all sets, including sets that do not contain zzz, where neither Equation (10) nor the axioms on BBB alone say anything about zzz. A second subtlety is the diagonal: {x,x}={x}\{x,x\}=\{x\}{x,x}={x}, so pxxp_{xx}pxx​ is set to 12\tfrac1221​ by definition rather than read off a singleton choice.

Formalization scope

Alternatives form a type X with decidable equality, alternative sets are Finset X, possible sets are a Set (Finset X), and selection probabilities are a real-valued function P : S → Finset X → X → ℝ. Only values P s B x with x ∈ B and B possible are constrained; no statement depends on the others. binProb sets the diagonal to 1/2. altSetV uses Real.log, which is 0 on non-positive arguments; under Axiom 2 on the binary sets its argument is always positive where it is used.

The probability-vector hypothesis on every possible set, binary sets included, is part of every statement: without it the zero function satisfies Axiom 1 vacuously and Equations (7)–(8) fail. The goal is stated with the explicit v(s,x)=V(s,x,z)v(s,x)=V(s,x,z)v(s,x)=V(s,x,z), never as "for each BBB there is a vvv", which would only restate (10).

Nothing beyond Mathlib's finite sums, Real.exp and Real.log is needed. Proofs of the milestones and of the goal are welcome, as is a formal statement of the auto/bus example or of the converse (logit selection probabilities satisfy Axioms 1 and 2).

Selected references

  • D. McFadden, Conditional logit analysis of qualitative choice behavior, in P. Zarembka (ed.), Frontiers in Econometrics, Academic Press, New York, 1974, pp. 105–142. https://eml.berkeley.edu/reprints/mcfadden/zarembka.pdf
  • R. D. Luce, Individual Choice Behavior: A Theoretical Analysis, Wiley, New York, 1959. https://doi.org/10.1037/14396-000
7 thms2 active usersReviewed
Next

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