Prove2Me
Navigate
DiscoverCollectionsFormalpediaBlogsUsersMomentumMy Missions+
Prove2Me
⌕
Log in
← Collections

Revenue Management and Choice Models

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

54 missions

Missions

41–54 of 54
OpenCompletedAll
Algorithmic Game TheoryOperations ResearchOptimization·Captain: mikedeng1

Purchasing, Pricing, and Quick Response in the Presence of Strategic Consumers: Under Condition (6), Quick Response Is More Valuable with Strategic Consumers than with Only Myopic OnesResearch Paper

Motivation

Fashion and consumer-electronics retailers sell a product at full price early in a season and mark down what is left. Consumers learn the pattern, and some of them wait for the markdown. Such strategic consumers lower the revenue of the full-price period, and the retailer's stocking decision affects how deep the markdown is expected to be. Quick response — a second, more expensive replenishment placed after demand is observed — is usually valued as a way to match supply with exogenous demand (Fisher and Raman 1996; Cachon and Terwiesch, Matching Supply with Demand, 2005). Cachon and Swinney ask how strategic waiting changes that value.

The source is the authors' working paper of April 2007, revised November 25, 2007, not the 2009 Management Science version, whose numbering and wording may differ. Its answer: with strategic consumers the retailer stocks less (Theorem 1), and under an explicit cost condition, quick response is worth more to a retailer facing strategic consumers than to one facing only myopic consumers (Theorem 3).

Setting

A retailer sells over two periods. It sells at the exogenous full price ppp in period 1 and at a markdown price s∈[0,p]s\in[0,p]s∈[0,p] chosen at the start of period 2. Leftover units are worth 000. First-period demand D≥0D\ge0D≥0 has density fff and distribution function FFF, and fff satisfies the monotone scaled likelihood ratio (MSLR) property: for every λ∈(0,1]\lambda\in(0,1]λ∈(0,1], x↦f(λx)/f(x)x\mapsto f(\lambda x)/f(x)x↦f(λx)/f(x) is monotonic on the support of fff.

The market has three segments:

  • myopic consumers, (1−α)D(1-\alpha)D(1−α)D of them, with value vMv_MvM​, who only buy in period 1;
  • strategic consumers, αD\alpha DαD of them, with value vMv_MvM​ in period 1 and second-period values uniform on [v‾,vˉ][\underline v,\bar v][v​,vˉ];
  • an unlimited pool of bargain hunters with value vBv_BvB​, who only buy on sale.

The standing assumptions are vˉ≤p\bar v\le pvˉ≤p and v‾≥vM−p+vB\underline v\ge v_M-p+v_Bv​≥vM​−p+vB​. Let Gˉ(s)\bar G(s)Gˉ(s) be the fraction of strategic values above sss.

By a threshold argument (Lemma 1), strategic consumers with value below some v^\hat vv^ buy at ppp and the rest wait. A fraction ξ=1−Gˉ(v^)α\xi=1-\bar G(\hat v)\alphaξ=1−Gˉ(v^)α of demand then buys in period 1, and the inventory left for period 2 is I=(q−ξD)+I=(q-\xi D)^+I=(q−ξD)+. The period-2 revenue R(s,I)R(s,I)R(s,I) counts the waiting strategic consumers with value at least sss and, if s≤vBs\le v_Bs≤vB​, the bargain hunters, up to the inventory III. The retailer's expected profit at unit cost ccc is

π(q,v^)=E[pmin⁡(q,ξD)−cq+sup⁡0≤s≤pR(s,I)].\pi(q,\hat v)=\mathbb E\Big[p\min(q,\xi D)-cq+\sup_{0\le s\le p}R(s,I)\Big].π(q,v^)=E[pmin(q,ξD)−cq+0≤s≤psup​R(s,I)].

With quick response, units ordered before the season cost c1c_1c1​ and units ordered after observing DDD cost c2c_2c2​, with c1≤c2≤pc_1\le c_2\le pc1​≤c2​≤p. The second order covers all first-period demand and may add stock for the sale. The resulting profit is πr(q,v^)\pi_r(q,\hat v)πr​(q,v^).

In the sale period, waiting strategic consumers are rationed: they effectively face the inventory θI\theta IθI, where θ∈[0,1]\theta\in[0,1]θ∈[0,1] measures their place in the queue. A strategic consumer with value v^\hat vv^ who waits gains, in expectation,

ψ(v^)=(v^−vB)Pr⁡(D<Dl and a unit is received),\psi(\hat v)=(\hat v-v_B)\Pr(D<D_l\text{ and a unit is received}),ψ(v^)=(v^−vB​)Pr(D<Dl​ and a unit is received),

where DlD_lDl​ is the demand level below which the retailer clears stock at sl=vBs_l=v_Bsl​=vB​. A rational expectations equilibrium (q∗,v∗)(q^*,v^*)(q∗,v∗) is a pair in which q∗q^*q∗ maximizes π(⋅,v∗)\pi(\cdot,v^*)π(⋅,v∗) and v∗v^*v∗ is a best response of consumers who correctly expect q∗q^*q∗. The superscript mmm denotes the benchmark with only myopic consumers (α=0\alpha=0α=0): πm\pi^mπm and πrm\pi^m_rπrm​ are the optimal myopic profits without and with quick response.

Formalization targets

Goal: Theorem 3

Assume MSLR and no rationing, 0<α≤10<\alpha\le10<α≤1, vB<c1<pv_B<c_1<pvB​<c1​<p, c1≤c2≤pc_1\le c_2\le pc1​≤c2​≤p, and condition (6):

vM−pvˉ−vB ≥ c2−c1c2−vB.\frac{v_M-p}{\bar v-v_B}\ \ge\ \frac{c_2-c_1}{c_2-v_B}.vˉ−vB​vM​−p​ ≥ c2​−vB​c2​−c1​​.

Let (q∗,v∗)(q^*,v^*)(q∗,v∗) be any equilibrium without quick response, (qr∗,vr∗)(q_r^*,v_r^*)(qr∗​,vr∗​) any equilibrium with it, and πm\pi^mπm, πrm\pi_r^mπrm​ the myopic optima. Then

πr(qr∗,vr∗)−π(q∗,v∗) ≥ πrm−πm.\pi_r(q_r^*,v_r^*)-\pi(q^*,v^*)\ \ge\ \pi_r^m-\pi^m .πr​(qr∗​,vr∗​)−π(q∗,v∗) ≥ πrm​−πm.

Milestones

The milestones follow the paper's path:

  • the threshold structure (Lemma 1);
  • the optimal sale price (Lemma 2) and quasi-concavity of π\piπ with first-order condition (2) (Lemma 3);
  • the fill probability and the limits of the best response (Lemma 4);
  • existence and the comparison q∗≤qmq^*\le q^mq∗≤qm, π∗≤πm\pi^*\le\pi^mπ∗≤πm (Theorem 1), with the myopic newsvendor F(qm)=(p−c)/(p−vB)F(q^m)=(p-c)/(p-v_B)F(qm)=(p−c)/(p−vB​);
  • the quick-response analogues (Lemma 5, Theorem 2 (i)), with the myopic fractile F(qrm)=(c2−c1)/(c2−vB)F(q^m_r)=(c_2-c_1)/(c_2-v_B)F(qrm​)=(c2​−c1​)/(c2​−vB​);
  • the statement that under (6) every equilibrium with quick response has vr∗=vˉv^*_r=\bar vvr∗​=vˉ (Theorem 2, last sentence).

Corollary 1 is the percentage form, Δ/π∗≥Δm/πm\Delta/\pi^*\ge\Delta_m/\pi^mΔ/π∗≥Δm​/πm.

Significance

Theorem 3 identifies a second channel through which quick response creates value. Beyond matching supply to demand, it lets the retailer keep its initial stock low enough that a deep markdown becomes unlikely, so strategic consumers buy at full price. Under (6), all of them do. Quick response thus reduces strategic waiting without withholding availability, unlike the inventory-signalling remedies in the literature, and the theorem quantifies when this effect dominates.

The results are proved in the working paper, partly in a technical appendix. No machine-checked version of them, or of the underlying markdown game, is known. A formalization pins down several statements that the paper states loosely:

  • the uniqueness claims of Lemmas 2 and 5;
  • the case condition of Lemma 4 (i), which is false as printed;
  • the sign in display (5);
  • the boundary cases of the threshold lemma.

It also produces reusable components: the newsvendor with salvage and the reactive-capacity fractile under a general density, and a rational-expectations equilibrium predicate for a retailer–consumer game.

Difficulty

The profit π(⋅,v^)\pi(\cdot,\hat v)π(⋅,v^) is not concave: with strategic consumers it is concave–convex (Figure 4 of the paper). The newsvendor argument therefore does not give a unique optimal order, and Lemma 3's quasi-concavity rests on MSLR in a short appendix step.

Existence (Theorem 1) needs a fixed point of the map q↦q\mapstoq↦ best response, but the consumer best response is a correspondence, not a function, so the printed intermediate-value argument does not apply directly. Theorem 3 needs a statement about every equilibrium with quick response, while Theorem 2's proof only exhibits one. Ruling out an equilibrium with vr∗<vˉv_r^*<\bar vvr∗​<vˉ requires comparing the derivative (4) of πr\pi_rπr​ with the myopic derivative along the whole demand distribution.

Formalization scope

Lean represents prices, quantities and valuations as reals, demand by a density f:R→Rf:\mathbb R\to\mathbb Rf:R→R on [0,∞)[0,\infty)[0,∞), and expectations as Lebesgue integrals against fff. The model carries the standing assumptions of §3 as fields, plus the following additions and corrections, each disclosed in the item notes:

  • vB>0v_B>0vB​>0: DlD_lDl​ divides by sl=vBs_l=v_Bsl​=vB​.
  • p<vMp<v_Mp<vM​, strengthening vM≥pv_M\ge pvM​≥p: Lemma 4 (ii) and Theorem 2's last claim need it.
  • A finite mean: πr\pi_rπr​ contains E[pξD]\mathbb E[p\xi D]E[pξD].
  • The paper's "θc≤θ\theta_c\le\thetaθc​≤θ" (p. 15) is replaced by the no-rationing condition slGˉ(v^)≤θsmGˉ(sm)s_l\bar G(\hat v)\le\theta s_m\bar G(s_m)sl​Gˉ(v^)≤θsm​Gˉ(sm​) for every belief, which is exactly Dl≤DθD_l\le D_\thetaDl​≤Dθ​. The printed condition agrees with it only when sm=v^s_m=\hat vsm​=v^.
  • Lemma 4 (i) is stated with the corrected case split.
  • Lemmas 2 and 5 claim uniqueness only off the tie points.
  • c<pc<pc<p is the reading of the "p−c>0p-c>0p−c>0" step in the proof of Theorem 1.
  • In the quick-response profit, ξD\xi DξD replaces the DDD that the proof of Theorem 2 prints.

Optimal revenues are suprema over all prices s∈[0,p]s\in[0,p]s∈[0,p] (and q2≥0q_2\ge0q2​≥0 with quick response). "Optimal order" means a maximizer over all q≥0q\ge0q≥0, never a stationary point, and the myopic benchmarks are the same functions at α=0\alpha=0α=0. Defining the optimal revenue by Lemma 2's closed form would make Lemma 2 and the first-order conditions definitional; it is not done.

The fill rate is min⁡{(1−ξ)x,θI}/((1−ξ)x)\min\{(1-\xi)x,\theta I\}/((1-\xi)x)min{(1−ξ)x,θI}/((1−ξ)x), set to 111 when no strategic consumer waits. With Lean's 0/0=00/0=00/0=0 instead, vˉ\bar vvˉ would be a best response to every order and Theorem 2's last claim would be trivial.

A complete development needs:

  • differentiation under the integral for piecewise-smooth integrands;
  • quasi-concavity from a single-crossing derivative;
  • a fixed-point argument for the equilibrium correspondence;
  • the newsvendor and reactive-capacity fractiles.

The last two are reusable beyond this mission. Proofs of any milestone, alternative existence arguments, and sorry-free proofs of the newsvendor items are welcome. The comparison "qr∗≤q∗q_r^*\le q^*qr∗​≤q∗, πr∗≥π∗\pi_r^*\ge\pi^*πr∗​≥π∗" of Theorem 2 and §8's numerical study are outside the scope.

Selected references

  • G. P. Cachon, R. Swinney, Purchasing, Pricing, and Quick Response in the Presence of Strategic Consumers, working paper, revised November 25, 2007; published in Management Science 55(3), 2009. https://doi.org/10.1287/mnsc.1080.0948
  • M. L. Fisher, A. Raman, Reducing the Cost of Demand Uncertainty Through Accurate Response to Early Sales, Operations Research 44(1), 1996. https://doi.org/10.1287/opre.44.1.87
  • J. F. Muth, Rational Expectations and the Theory of Price Movements, Econometrica 29(3), 1961. https://doi.org/10.2307/1909635
19 thms2 active usersReviewed
Dynamic ProgrammingOperations ResearchOptimization+1·Captain: mikedeng1

Coordinating Inventory Control and Pricing Strategies with Random Demand and Fixed Ordering Cost: The Finite Horizon Case 2: General Demand: Sym-k-Concave Profit-to-Go, Optimal (s, S, A, p) PolicyResearch Paper

Motivation

A seller with replenishable inventory must decide how much to order and what price to charge before demand is known. A price that increases current revenue can also change the inventory left for later periods; a fixed charge for placing any positive order adds a discontinuity to the decision. Chen and Simchi-Levi study this coordination problem over a finite horizon with backlogging and random, price-dependent demand. Their general-demand result characterizes an optimal order-and-price policy even when a conventional two-threshold ordering rule can fail. The paper gives explicit counterexamples to ordinary kkk-concavity and to optimality of a standard (s,S,p)(s,S,p)(s,S,p) policy in this setting (Chen and Simchi-Levi, 2004, §4, Lemmas 3–4).

Setting

There are periods t=1,…,Tt=1,\ldots,Tt=1,…,T. At the start of period ttt, the firm has inventory xxx, which may be negative because unmet demand is backlogged. It may order to any level y≥xy\ge xy≥x. If y>xy>xy>x, it pays a fixed ordering cost kkk as well as a variable cost ct(y−x)c_t(y-x)ct​(y−x). After ordering it chooses an expected demand level ddd in a closed interval [d‾t,d‾t][\underline d_t,\overline d_t][d​t​,dt​]. The corresponding price is Pt(d)P_t(d)Pt​(d), the inverse of the period's decreasing demand curve, and expected revenue is Rt(d)=dPt(d)R_t(d)=dP_t(d)Rt​(d)=dPt​(d).

Actual demand is αtd+βt\alpha_t d+\beta_tαt​d+βt​. The random pair (αt,βt)(\alpha_t,\beta_t)(αt​,βt​) has a period-dependent law μt\mu_tμt​, with E[αt]=1\mathbb E[\alpha_t]=1E[αt​]=1 and E[βt]=0\mathbb E[\beta_t]=0E[βt​]=0. No sign restriction is placed on αt\alpha_tαt​ in the paper's general model. The inventory after demand is y−αtd−βty-\alpha_t d-\beta_ty−αt​d−βt​, incurring a convex holding or backlog cost hth_tht​ and entering the next period. The paper assumes the random perturbations are independent across periods. Its Bellman equation uses their individual period laws (Chen and Simchi-Levi, 2004, §2, Assumptions 1–5).

The profit-to-go vt(x)v_t(x)vt​(x) starts with vT+1(x)=0v_{T+1}(x)=0vT+1​(x)=0. For a post-order inventory yyy and expected demand ddd, write

gt(y,d)=Rt(d)−cty+E[−ht(y−αtd−βt)+vt+1(y−αtd−βt)].g_t(y,d)=R_t(d)-c_ty+\mathbb E[-h_t(y-\alpha_t d-\beta_t)+v_{t+1}(y-\alpha_t d-\beta_t)].gt​(y,d)=Rt​(d)−ct​y+E[−ht​(y−αt​d−βt​)+vt+1​(y−αt​d−βt​)].

The firm maximizes gt(y,d)g_t(y,d)gt​(y,d) over admissible ddd and then maximizes the resulting value, less the fixed ordering charge, over y≥xy\ge xy≥x. This defines vt(x)v_t(x)vt​(x) by equation (2) of the paper. A symmetrically kkk-convex real function fff satisfies

f((1−λ)x0+λx1)≤(1−λ)f(x0)+λf(x1)+max⁡{λ,1−λ}kf((1-\lambda)x_0+\lambda x_1)\le(1-\lambda)f(x_0)+\lambda f(x_1)+\max\{\lambda,1-\lambda\}kf((1−λ)x0​+λx1​)≤(1−λ)f(x0​)+λf(x1​)+max{λ,1−λ}k

for every x0,x1∈Rx_0,x_1\in\mathbb Rx0​,x1​∈R and λ∈[0,1]\lambda\in[0,1]λ∈[0,1]. A function is symmetrically kkk-concave when its negative is symmetrically kkk-convex. This is Definition 4.1, equation (8), and preserves the symmetry between the two endpoint inventories (Chen and Simchi-Levi, 2004, p. 891).

Formalization targets

The goal is Theorem 4.1(c)–(d). If Gt(y)=max⁡d∈[d‾t,d‾t]gt(y,d)G_t(y)=\max_{d\in[\underline d_t,\overline d_t]}g_t(y,d)Gt​(y)=maxd∈[d​t​,dt​]​gt​(y,d), then for every period

Gt and vt are symmetrically k-concave.G_t\text{ and }v_t\text{ are symmetrically }k\text{-concave}.Gt​ and vt​ are symmetrically k-concave.

There are thresholds st≤Sts_t\le S_tst​≤St​ and a possibly empty set At⊆[st,(st+St)/2]A_t\subseteq[s_t,(s_t+S_t)/2]At​⊆[st​,(st​+St​)/2]. An optimal policy orders to StS_tSt​ if x<stx<s_tx<st​ or x∈Atx\in A_tx∈At​, and places no order otherwise. It chooses the maximizing expected demand for the resulting post-order level. In particular, all states that order to StS_tSt​ can use the same demand level dt(St)d_t(S_t)dt​(St​) and hence the same price Pt(dt(St))P_t(d_t(S_t))Pt​(dt​(St​)).

The milestones state the preceding polynomial growth and continuity assertions from Theorem 4.1(a)–(b), the threshold structure of Lemma 5(d), and the named steps in the proof on p. 892. The link from the paper's earlier kkk-convexity to symmetric kkk-convexity is also included. These targets preserve the paper's general demand law rather than restricting it to the additive case.

Significance

The theorem describes an optimal decision at every inventory level under random price-dependent demand. The exceptional set AtA_tAt​ records states where ordering to StS_tSt​ is optimal inside the interval where a single reorder threshold need not describe all optimal actions. Thus the result retains a concise policy description while allowing the behavior exhibited by the paper's counterexamples (Chen and Simchi-Levi, 2004, §4).

Formalizing the result requires a reusable account of symmetric kkk-convexity, a fixed-cost ordering envelope, and a finite-horizon Bellman recursion with real-valued expectations. The source proves Theorem 4.1 on paper. This mission poses its statements as Lean goals; the theorem files currently contain sorry and therefore are not machine-checked proofs. The published platform definition BertsekasKConvex supplies the earlier notion from Definition 2.1, so the bridge to Definition 4.1 can be stated without duplicating it.

Difficulty

In the additive-demand case, a maximizing expected-demand choice can be selected so that post-demand inventory changes monotonically with the initial inventory. For general demand αtd+βt\alpha_t d+\beta_tαt​d+βt​, that monotonicity can fail: the random multiplier changes how a pricing choice moves the next inventory state. The earlier ordered-chord kkk-concavity argument therefore does not carry over unchanged. A second issue is the fixed ordering charge. It creates ties between ordering and not ordering, so the policy region can include isolated or noninterval states within [st,(st+St)/2][s_t,(s_t+S_t)/2][st​,(st​+St​)/2]. The statement must assert Bellman optimality at every xxx, including those ties, rather than only the geometry of AtA_tAt​ (Chen and Simchi-Levi, 2004, pp. 891–892).

Formalization scope

The Lean model uses expected demand ddd as the decision variable and obtains price from Pt(d)P_t(d)Pt​(d), matching the paper's reformulation after Assumption 5. The admissible demand interval is nonempty. Periods are natural numbers 1,…,T1,\ldots,T1,…,T, with T>0T>0T>0 and terminal value vT+1=0v_{T+1}=0vT+1​=0. Random demand is integrated against a probability measure on R×R\mathbb R\times\mathbb RR×R. The first moments and the paper's demand moment condition are explicit; holding-cost and continuation integrability prevent a nonintegrable real integral from acquiring Lean's default zero value. The paper's O(∣y∣ρ)O(|y|^\rho)O(∣y∣ρ) statements are encoded by constants multiplying 1+∣y∣ρ1+|y|^\rho1+∣y∣ρ.

The formal assumptions include k≥0k\ge0k≥0 and ct≥0c_t\ge0ct​≥0, and set cT+1=0c_{T+1}=0cT+1​=0 where the paper's Assumption 3 uses that otherwise undefined terminal coefficient. The latter two restrictions make the printed finite-horizon claim valid under its cost interpretation; they are stated rather than silently supplied. The paper assumes temporal independence, while the formal Bellman recursion starts from the marginal laws and so needs no separate joint process. The paper's (y,p)(y,p)(y,p) in Theorem 4.1(b) is expressed as (y,d)(y,d)(y,d) through the continuous one-to-one price–expected-demand correspondence. The model places no positivity condition on αt\alpha_tαt​.

The optimized demand value and ordering value are real suprema. The statements include maximizing choices and Bellman attainment; a proof cannot rely on a default value for an empty or unbounded supremum. The goal requires an actual optimal action for every inventory state. Contributions that establish integrability, continuity, coercivity, symmetric convexity preservation, or the ordering-envelope structure are useful beyond this mission.

Selected references

  • Xin Chen and David Simchi-Levi, Coordinating Inventory Control and Pricing Strategies with Random Demand and Fixed Ordering Cost: The Finite Horizon Case, Operations Research 52(6), 887–896, 2004. DOI: 10.1287/opre.1040.0127.
12 thms3 active usersReviewed
Operations ResearchOptimization·Captain: mikedeng1

Assortment Optimization under Variants of the Nested Logit Model 3: With Fully-Captured Nests, the Nested-by-Revenue LP Optimum Scaled by the Factor (6) Is Feasible for the Full LPResearch Paper

Motivation

Assortment optimization asks a retailer which products to offer when customers substitute among them. The nested logit model groups products into nests and is one of the most used choice models in revenue management, because it relaxes the independence of irrelevant alternatives of the plain multinomial logit model while keeping choice probabilities in closed form. Davis, Gallego and Topaloglu (Oper. Res. 62(2), 2014) study how the tractability of the assortment problem under this model depends on two features: whether the dissimilarity parameters of the nests are at most one, and whether a customer who selects a nest always buys there (fully-captured nests).

When the dissimilarity parameters are at most one and the nests are fully captured, offering the jjj highest-revenue products in every nest is optimal (Theorem 4 of the paper). This mission concerns what survives when a dissimilarity parameter exceeds one, the regime the paper calls possibly synergistic products. Then the problem is NP-hard (Theorem 5), and the paper shows that the same nested-by-revenue assortments still achieve an explicit, data-dependent fraction of the optimal expected revenue.

Setting

There are nests i∈Mi \in Mi∈M and products j∈N={1,…,n}j \in N = \{1, \dots, n\}j∈N={1,…,n} in every nest. Product jjj of nest iii has a revenue rijr_{ij}rij​ and a preference weight vij>0v_{ij} > 0vij​>0, ordered so that ri1≥ri2≥⋯≥rinr_{i1} \ge r_{i2} \ge \dots \ge r_{in}ri1​≥ri2​≥⋯≥rin​. Nest iii has a dissimilarity parameter γi>0\gamma_i > 0γi​>0, and v0≥0v_0 \ge 0v0​≥0 is the weight of leaving without choosing a nest. Throughout this mission the nests are fully captured: the within-nest no-purchase weights vi0v_{i0}vi0​ are zero. For an assortment Si⊆NS_i \subseteq NSi​⊆N in nest iii,

Vi(Si)=∑j∈Sivij,Ri(Si)=∑j∈SirijvijVi(Si),Ri(∅)=0,V_i(S_i) = \sum_{j \in S_i} v_{ij}, \qquad R_i(S_i) = \frac{\sum_{j \in S_i} r_{ij} v_{ij}}{V_i(S_i)}, \quad R_i(\emptyset) = 0,Vi​(Si​)=j∈Si​∑​vij​,Ri​(Si​)=Vi​(Si​)∑j∈Si​​rij​vij​​,Ri​(∅)=0,

and the expected revenue of (S1,…,Sm)(S_1, \dots, S_m)(S1​,…,Sm​) is

Π(S1,…,Sm)=∑i∈MVi(Si)γiRi(Si)v0+∑i∈MVi(Si)γi.\Pi(S_1, \dots, S_m) = \frac{\sum_{i \in M} V_i(S_i)^{\gamma_i} R_i(S_i)}{v_0 + \sum_{i \in M} V_i(S_i)^{\gamma_i}}.Π(S1​,…,Sm​)=v0​+∑i∈M​Vi​(Si​)γi​∑i∈M​Vi​(Si​)γi​Ri​(Si​)​.

Problem (2) maximizes Π\PiΠ; its optimal value is Z∗Z^*Z∗. The nested-by-revenue assortment Nij={1,…,j}N_{ij} = \{1, \dots, j\}Nij​={1,…,j} collects the jjj highest-revenue products of nest iii, with Ni0=∅N_{i0} = \emptysetNi0​=∅ and N+={0,1,…,n}N_+ = \{0, 1, \dots, n\}N+​={0,1,…,n}.

Problem (2) is equivalent to the linear program (3): minimize xxx subject to v0x≥∑iyiv_0 x \ge \sum_i y_iv0​x≥∑i​yi​ and yi≥Vi(Si)γi(Ri(Si)−x)y_i \ge V_i(S_i)^{\gamma_i}(R_i(S_i) - x)yi​≥Vi​(Si​)γi​(Ri​(Si​)−x) for every nest iii and every Si⊆NS_i \subseteq NSi​⊆N. Problem (4) keeps the second family of constraints only for candidate assortments; here the candidates are {Nij:j∈N+}\{N_{ij} : j \in N_+\}{Nij​:j∈N+​}, which gives a linear program with 1+m1 + m1+m variables and 1+m(1+n)1 + m(1 + n)1+m(1+n) constraints. The performance factor of the nested-by-revenue assortments is

α=max⁡i∈M, j=2,…,n{Ri(Ni,j−1)Ri(Nij)∧(Ri(Nij)Ri(Ni,j−1) Vi(Nij)γiVi(Ni,j−1)γi)},a∧b=min⁡{a,b}.(6)\alpha = \max_{i \in M,\ j = 2, \dots, n} \left\{ \frac{R_i(N_{i,j-1})}{R_i(N_{ij})} \wedge \left(\frac{R_i(N_{ij})}{R_i(N_{i,j-1})}\, \frac{V_i(N_{ij})^{\gamma_i}}{V_i(N_{i,j-1})^{\gamma_i}}\right) \right\}, \qquad a \wedge b = \min\{a, b\}. \tag{6}α=i∈M, j=2,…,nmax​{Ri​(Nij​)Ri​(Ni,j−1​)​∧(Ri​(Ni,j−1​)Ri​(Nij​)​Vi​(Ni,j−1​)γi​Vi​(Nij​)γi​​)},a∧b=min{a,b}.(6)

Formalization targets

Goal: Theorem 7

Assume vi0=0v_{i0} = 0vi0​=0 for every nest, γi>1\gamma_i > 1γi​>1 for some nest, positive revenues and n≥2n \ge 2n≥2. If (x^,y^)(\hat x, \hat y)(x^,y^​) is an optimal solution of (4) over the nested-by-revenue assortments, then

(αx^,αy^) is feasible for (3).(\alpha \hat x, \alpha \hat y) \ \text{is feasible for (3)}.(αx^,αy^​) is feasible for (3).

Combined with Theorem 1 of the paper this yields Z∗≤α Π(S^)Z^* \le \alpha\, \Pi(\hat S)Z∗≤αΠ(S^) for the assortment S^\hat SS^ read off the small linear program; that consequence is a companion item.

Milestones

  1. Problem (3) is a relaxation of the fractional problem (7), in which yiy_iyi​ dominates (∑jvijzij)γi[∑jrijvijzij/∑jvijzij−x]\big(\sum_j v_{ij} z_{ij}\big)^{\gamma_i}\big[\sum_j r_{ij} v_{ij} z_{ij} / \sum_j v_{ij} z_{ij} - x\big](∑j​vij​zij​)γi​[∑j​rij​vij​zij​/∑j​vij​zij​−x] for all zi∈[0,1]nz_i \in [0,1]^nzi​∈[0,1]n.
  2. Lemma 6: the inner maximization (8) of (7) has a solution of the form zi1=⋯=zi,k−1=1z_{i1} = \dots = z_{i,k-1} = 1zi1​=⋯=zi,k−1​=1, zik∈[0,1]z_{ik} \in [0,1]zik​∈[0,1], zi,k+1=⋯=zin=0z_{i,k+1} = \dots = z_{in} = 0zi,k+1​=⋯=zin​=0.
  3. y^i≥0\hat y_i \ge 0y^​i​≥0 and x^≥0\hat x \ge 0x^≥0.
  4. The inequalities (20) and (22) with the coefficients αik1=Ri(Ni,k−1)/Ri(Nik)\alpha^1_{ik} = R_i(N_{i,k-1})/R_i(N_{ik})αik1​=Ri​(Ni,k−1​)/Ri​(Nik​) and 1∨αik21 \vee \alpha^2_{ik}1∨αik2​.
  5. Ri(Nij)≤Ri(Ni,j−1)R_i(N_{ij}) \le R_i(N_{i,j-1})Ri​(Nij​)≤Ri​(Ni,j−1​), and Lemma 14: α≥1\alpha \ge 1α≥1 when some γi>1\gamma_i > 1γi​>1.
  6. The inequalities (24) and (25): αy^i\alpha \hat y_iαy^​i​ dominates the objective of (8) at αx^\alpha \hat xαx^ for every vector of Lemma 6's shape.

Two companion statements follow the goal: the bounds α≤ρ\alpha \le \rhoα≤ρ and α≤2κ\alpha \le 2\kappaα≤2κ when revenues, respectively preference weights, within each nest differ by at most the factors ρ\rhoρ, κ\kappaκ; and the guarantee Z∗≤α Π(S^)Z^* \le \alpha\, \Pi(\hat S)Z∗≤αΠ(S^).

Significance

Because problem (2) is NP-hard once a dissimilarity parameter exceeds one (Theorem 5 of the paper), an exact polynomial algorithm is not expected, and a guarantee for a polynomial-size candidate family is the natural substitute. Theorem 7 gives one with an explicit factor computed from the data: α≤ρ\alpha \le \rhoα≤ρ when revenues within each nest are balanced, α≤2κ\alpha \le 2\kappaα≤2κ when preference weights within each nest are balanced, and the nests may differ arbitrarily from one another. The same linear-programming argument (Theorem 1) is reused for partially-captured nests in §6 of the paper.

The result is proved in the paper (Appendix A.1); no machine-checked version is known. A formalization verifies an appendix argument that is stated with little detail, and fixes the conventions that the page leaves implicit, such as the value of the objective of (8) at zi=0z_i = 0zi​=0 and the role of the standing assumption that some γi\gamma_iγi​ exceeds one.

Difficulty

Theorem 7 is a feasibility statement for the exponentially many constraints of (3), one per subset of every nest, while optimality of (x^,y^)(\hat x, \hat y)(x^,y^​) only controls the n+1n + 1n+1 nested-by-revenue constraints per nest. The obvious approach, comparing each subset directly with a nested-by-revenue assortment, fails: when γi>1\gamma_i > 1γi​>1 the nested-by-revenue assortments are not optimal within a nest, and the example of §4.1 shows that they can lose an unbounded factor. The constant α\alphaα must absorb the gap between a fractional assortment and its two neighbouring nested-by-revenue assortments, which the dual-style solution (x^,y^)(\hat x, \hat y)(x^,y^​) of the small program does not see.

Formalization scope

Lean 4 with Mathlib. Nests form a finite type ι; products are Fin n, indexed from 000, so NijN_{ij}Nij​ is nbr n j = {k | k < j} and product kkk of the page is ⟨k - 1, _⟩. Powers are real powers (Real.rpow), and division is total (x/0=0x / 0 = 0x/0=0), which gives Ri(∅)=0R_i(\emptyset) = 0Ri​(∅)=0 and value 000 for the objective of (8) at zi=0z_i = 0zi​=0.

Standing assumptions in every statement: v0≥0v_0 \ge 0v0​≥0, vi0≥0v_{i0} \ge 0vi0​≥0, vij>0v_{ij} > 0vij​>0, rij≥0r_{ij} \ge 0rij​≥0, γi>0\gamma_i > 0γi​>0, revenues ordered within each nest (§1); vi0=0v_{i0} = 0vi0​=0 for every nest and γi>1\gamma_i > 1γi​>1 for some nest (§4, p. 16). Added and disclosed: rij>0r_{ij} > 0rij​>0 wherever (6) appears, so that no denominator of (6) vanishes; n≥2n \ge 2n≥2, the nonempty range of (6); v0>0v_0 > 0v0​>0 only in the guarantee Z∗≤α Π(S^)Z^* \le \alpha\,\Pi(\hat S)Z∗≤αΠ(S^), inherited from Theorem 1. The factor α\alphaα enters as a real number with the hypothesis that it is the greatest element of the set of terms of (6). An optimal solution of (4) is a feasible pair whose xxx is minimal among feasible pairs.

A trivializing formalization is ruled out: α\alphaα is the maximum of (6), not an arbitrary upper bound, and the goal states feasibility for the full program (3) over every subset of every nest, not for (7) or for the candidate collection.

Needed infrastructure: maximization of a continuous function on [0,1]n[0,1]^n[0,1]n, the greedy solution of a continuous knapsack, monotonicity of weighted averages and elementary real-power calculus. These are reusable beyond this mission. Proofs of any milestone, of the companions, and of the §4.1 example are welcome.

Selected references

  • J. M. Davis, G. Gallego, H. Topaloglu, Assortment optimization under variants of the nested logit model, Operations Research 62(2), 2014. https://doi.org/10.1287/opre.2014.1256 (formalized from the revised manuscript of June 18, 2013).
  • K. Talluri, G. 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
14 thms1 active userReviewed
Operations ResearchOptimization·Captain: mikedeng1

Assortment Optimization under Variants of the Nested Logit Model 5: For General Nests, the Nested-by-Preference-and-Revenue LP Optimum Scaled by the Factor (12) Is Feasible for the Full LPResearch Paper

Assortment planning with nested choice

A retailer that groups its products into categories (brands, store sections, flight classes) and decides which products to display in each faces the assortment problem: offering more products attracts more customers but also diverts sales away from the most profitable products. The nested logit model is the standard description of customer choice in this setting. A customer first selects a category (a nest), then a product within it. Choice-based models of this kind are the basis of revenue management under customer choice (Talluri and van Ryzin 2004).

Davis, Gallego and Topaloglu (DGT 2014) study the assortment problem under the nested logit model with two features that earlier work excluded: dissimilarity parameters larger than one, under which products in a nest act as complements rather than substitutes, and a no-purchase option inside each nest, under which a customer may enter a nest and still leave without buying. They show that the problem is NP-hard once either feature is present. For each regime they give a small linear program whose solution yields an assortment with a provable performance guarantee. This mission formalizes the guarantee for the most general instances, where both features occur together (§6.1, Theorem 11).

Timeline. Rusmevichientong, Shmoys and Topaloglu (2010) bound nested-by-revenue assortments under a multinomial logit mixture. [DGT 2014] prove that nested-by-revenue assortments are optimal for dissimilarity parameters at most one without within-nest no-purchase options (Theorem 4). They give factor-(6) guarantees with synergistic products, a factor-two guarantee via knapsack relaxations for partially-captured nests (Theorem 10), and the general factor (12) of Theorem 11. Li, Rusmevichientong and Topaloglu (2015) extend the nested-by-revenue result to ddd-level nested logit models.

The model

There are nests i∈M={1,…,m}i\in M=\{1,\dots,m\}i∈M={1,…,m} and, in each nest, products j∈N={1,…,n}j\in N=\{1,\dots,n\}j∈N={1,…,n}. Product jjj of nest iii has revenue rij≥0r_{ij}\ge0rij​≥0 and preference weight vij>0v_{ij}>0vij​>0, with ri1≥ri2≥⋯≥rinr_{i1}\ge r_{i2}\ge\dots\ge r_{in}ri1​≥ri2​≥⋯≥rin​. Nest iii has a no-purchase weight vi0≥0v_{i0}\ge0vi0​≥0 and a dissimilarity parameter γi>0\gamma_i>0γi​>0, and v0≥0v_0\ge0v0​≥0 is the weight of choosing no nest at all. For an assortment Si⊆NS_i\subseteq NSi​⊆N,

Vi(Si)=vi0+∑j∈Sivij,Ri(Si)=∑j∈SirijvijVi(Si).V_i(S_i)=v_{i0}+\sum_{j\in S_i}v_{ij},\qquad R_i(S_i)=\frac{\sum_{j\in S_i}r_{ij}v_{ij}}{V_i(S_i)} .Vi​(Si​)=vi0​+j∈Si​∑​vij​,Ri​(Si​)=Vi​(Si​)∑j∈Si​​rij​vij​​.

A customer chooses nest iii with probability Vi(Si)γi/(v0+∑lVl(Sl)γl)V_i(S_i)^{\gamma_i}/(v_0+\sum_l V_l(S_l)^{\gamma_l})Vi​(Si​)γi​/(v0​+∑l​Vl​(Sl​)γl​), and the expected revenue is

Π(S1,…,Sm)=∑iVi(Si)γiRi(Si)v0+∑iVi(Si)γi.\Pi(S_1,\dots,S_m)=\frac{\sum_{i}V_i(S_i)^{\gamma_i}R_i(S_i)}{v_0+\sum_{i}V_i(S_i)^{\gamma_i}} .Π(S1​,…,Sm​)=v0​+∑i​Vi​(Si​)γi​∑i​Vi​(Si​)γi​Ri​(Si​)​.

The optimal value Z∗Z^*Z∗ of max⁡Π\max\PimaxΠ equals the optimal value of the linear program

(3)min⁡ xs.t.v0x≥∑iyi,yi≥Vi(Si)γi(Ri(Si)−x)  ∀Si⊆N, i∈M,\text{(3)}\qquad \min\ x\quad\text{s.t.}\quad v_0x\ge\sum_i y_i,\qquad y_i\ge V_i(S_i)^{\gamma_i}\big(R_i(S_i)-x\big)\ \ \forall S_i\subseteq N,\ i\in M,(3)min xs.t.v0​x≥i∑​yi​,yi​≥Vi​(Si​)γi​(Ri​(Si​)−x)  ∀Si​⊆N, i∈M,

which has 2n2^n2n constraints per nest. Problem (4) keeps only the constraints for a chosen candidate collection of assortments in each nest.

A nest is fully captured if vi0=0v_{i0}=0vi0​=0 (i∈Mfi\in M^fi∈Mf) and partially captured if vi0>0v_{i0}>0vi0​>0 (i∈Mpi\in M^pi∈Mp). Nij={1,…,j}N_{ij}=\{1,\dots,j\}Nij​={1,…,j} is the nested-by-revenue assortment. NijkN^k_{ij}Nijk​ is the set of the jjj highest-revenue products among the kkk products of nest iii with the smallest preference weights, with Ni0k=∅N^k_{i0}=\emptysetNi0k​=∅ and Nijn=NijN^n_{ij}=N_{ij}Nijn​=Nij​.

Formalization targets

Goal: Theorem 11

Let (x^,y^)(\hat x,\hat y)(x^,y^​) be an optimal solution of (4) when the candidate collection of every nest is {Nijk:k∈N, j=0,…,k}∪{{j}:j∈N}\{N^k_{ij}:k\in N,\ j=0,\dots,k\}\cup\{\{j\}:j\in N\}{Nijk​:k∈N, j=0,…,k}∪{{j}:j∈N}, and let

β=max⁡i∈Mf, j=2,…,n{Vi(Nij)Vi(Ni,j−1)}∨max⁡i∈Mp, j=1,…,n{Vi(Nij)Vi(Ni,j−1)}∨2(12).\beta=\max_{i\in M^f,\ j=2,\dots,n}\left\{\frac{V_i(N_{ij})}{V_i(N_{i,j-1})}\right\}\vee\max_{i\in M^p,\ j=1,\dots,n}\left\{\frac{V_i(N_{ij})}{V_i(N_{i,j-1})}\right\}\vee2 \qquad (12).β=i∈Mf, j=2,…,nmax​{Vi​(Ni,j−1​)Vi​(Nij​)​}∨i∈Mp, j=1,…,nmax​{Vi​(Ni,j−1​)Vi​(Nij​)​}∨2(12).

Then (βx^,βy^)(\beta\hat x,\beta\hat y)(βx^,βy^​) is feasible for problem (3).

The theorem assumes γˉ=max⁡iγi>1\bar\gamma=\max_i\gamma_i>1γˉ​=maxi​γi​>1, as all of §6 does. Otherwise it places no restriction on the γi\gamma_iγi​ or the vi0v_{i0}vi0​.

Milestones

  1. x^≥0\hat x\ge0x^≥0 (A.4, p. 46).
  2. Every greedy knapsack assortment S^i(ϵi)\hat S_i(\epsilon_i)S^i​(ϵi​) of §5 is one of the NijkN^k_{ij}Nijk​ (pp. 24–25).
  3. The relaxed nest problem over [0,1]n[0,1]^n[0,1]n has an optimal solution of fractional-prefix form (A.4 Case 1, p. 47).
  4. Inequality (30): for a nest with γi>1\gamma_i>1γi​>1 and y^i≥0\hat y_i\ge0y^​i​≥0, βy^i\beta\hat y_iβy^​i​ bounds the relaxed objective at every fractional prefix (p. 47).
  5. Case 1: γi>1\gamma_i>1γi​>1, y^i≥0\hat y_i\ge0y^​i​≥0 gives the constraints of (3) for nest iii (pp. 47–48).
  6. Problem (31) has a nested-by-revenue optimal solution when γi>1\gamma_i>1γi​>1 and its coefficient b=βy^ib=\beta\hat y_ib=βy^​i​ is negative (p. 48).
  7. Case 2: γi>1\gamma_i>1γi​>1, y^i<0\hat y_i<0y^​i​<0 (p. 48).
  8. Case 3: γi≤1\gamma_i\le1γi​≤1, through the factor-two argument of Theorem 10 (pp. 48–49).

Two companions follow the goal. One is the resulting guarantee β Π(S^)≥Z∗≥Π(S^)\beta\,\Pi(\hat S)\ge Z^*\ge\Pi(\hat S)βΠ(S^)≥Z∗≥Π(S^), through Theorem 1. The other is the bound β≤2κ\beta\le2\kappaβ≤2κ when the preference weights within a nest differ by at most a factor κ\kappaκ (p. 26).

Significance

Theorem 11, combined with Theorem 1 of the paper, gives a polynomial-size method for an NP-hard problem. The method solves one linear program with 1+m1+m1+m variables and 1+m(1+n+n2)1+m(1+n+n^2)1+m(1+n+n2) constraints, then reads off an assortment whose expected revenue is within the factor β\betaβ of the optimum. This holds for every nested logit instance, including nests where customers may walk away and nests whose products are complements. When the weights inside each nest are within a factor κ\kappaκ of each other, the guarantee is at most 2κ2\kappa2κ.

The theorem is proved in the paper's appendix. No part of it is machine-checked. Formalizing it checks a case analysis that reuses, by reference, arguments from two other theorems: Theorem 7 (synergistic, fully-captured nests) and Theorem 10 (competitive, partially-captured nests). It makes precise what these arguments need when the two regimes are mixed in one instance. The formalization also fixes the boundary conventions the printed proof leaves implicit: fully-captured nests with k=1k=1k=1, zero-weight denominators, and the sign of y^i\hat y_iy^​i​.

Difficulty

Each nest falls into one of three regimes, and a different argument controls each. With γi≤1\gamma_i\le1γi​≤1 the nest behaves like a knapsack problem. Its guarantee of two needs the knapsack collection of §5 to sit inside {Nijk}\{N^k_{ij}\}{Nijk​}. With γi>1\gamma_i>1γi​>1 and y^i≥0\hat y_i\ge0y^​i​≥0, the constraint must be extended from nested-by-revenue sets to every subset. This goes through a continuous relaxation whose optimum has a fractional coordinate, and it costs the ratio Vi(Nik)/Vi(Ni,k−1)V_i(N_{ik})/V_i(N_{i,k-1})Vi​(Nik​)/Vi​(Ni,k−1​), which is where (12) comes from. With γi>1\gamma_i>1γi​>1 and y^i<0\hat y_i<0y^​i​<0, the scaling argument of Case 1 fails because multiplying by a factor at most one no longer preserves the inequality. The proof switches to the different objective (31), whose convexity in one coordinate forces an integral optimum.

The first idea, bounding every assortment by a nested-by-revenue one, is false here. With γi>1\gamma_i>1γi​>1 or vi0>0v_{i0}>0vi0​>0, nested-by-revenue assortments are not optimal, and the loss is exactly the factor β\betaβ.

Formalization scope

Products are Fin n; NijN_{ij}Nij​ is nbr n j. Powers are Real.rpow, and x/0=0x/0=0x/0=0, so Ri(∅)=0R_i(\emptyset)=0Ri​(∅)=0. Problems (3) and (4) are stated in constraint form: LP4Optimal means feasible and with xxx minimal among feasible points. β\betaβ is the greatest element of the finite set betaSet I, which contains 222 and the ratios of (12). Fully-captured nests skip j=1j=1j=1, as on the page. The collection is constructed: nestedPR breaks weight ties by index and revenue ties by index.

Standing assumptions, all disclosed:

  • vij>0v_{ij}>0vij​>0, rij≥0r_{ij}\ge0rij​≥0 and γi>0\gamma_i>0γi​>0. The page allows zero-weight padding products and γi=0\gamma_i=0γi​=0, but its arguments do not cover them.
  • γˉ>1\bar\gamma>1γˉ​>1 on every statement set in Theorem 11's context.
  • n≥1n\ge1n≥1 for the collection claim and the prefix claim.
  • vi0>0v_{i0}>0vi0​>0 for the statement about (31). That is the only kind of nest where Case 2 arises. For vi0=0v_{i0}=0vi0​=0, Lean's 01−γi=00^{1-\gamma_i}=001−γi​=0 would remove the page's +∞+\infty+∞.
  • v0>0v_0>0v0​>0 for the guarantee, where Theorem 1 fails otherwise.
  • κ≥1\kappa\ge1κ≥1, and vi0v_{i0}vi0​ counted among the weights of a partially-captured nest (vij≤κvi0v_{ij}\le\kappa v_{i0}vij​≤κvi0​, vi0≤κvijv_{i0}\le\kappa v_{ij}vi0​≤κvij​), for the 2κ2\kappa2κ bound.

A trivializing formalization is ruled out. The goal states only feasibility for (3), with β\betaβ the maximum of (12), not any upper bound. The collection is the page's, not an arbitrary family containing it. The goal mentions none of the cases or the relaxations.

A complete development needs continuous knapsack solutions (greedy optimality, fractional prefixes), convexity of t↦t1−γt\mapsto t^{1-\gamma}t↦t1−γ on (0,∞)(0,\infty)(0,∞), and the factor-two argument of Theorem 10. The knapsack and fractional-prefix lemmas are reusable for the companion missions of this series. Proofs of individual cases, and proofs of milestones in greater generality, are welcome.

Selected references

  • J. M. Davis, G. Gallego, H. Topaloglu, Assortment optimization under variants of the nested logit model, Operations Research 62(2), 2014 (revised manuscript of June 18, 2013). https://doi.org/10.1287/opre.2014.1256
  • P. Rusmevichientong, D. B. Shmoys, H. Topaloglu, Assortment optimization with mixtures of logits, technical report, Cornell University, 2010. http://legacy.orie.cornell.edu/~huseyin/publications/publications.html
  • G. Li, P. Rusmevichientong, H. Topaloglu, The d-level nested logit model: assortment and price optimization problems, Operations Research 63(2), 2015.
  • 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
  • D. P. Williamson, D. B. Shmoys, The Design of Approximation Algorithms, Cambridge University Press, 2011. https://doi.org/10.1017/CBO9780511921735
14 thms1 active userReviewed
Linear OptimizationOperations ResearchProbability·Captain: mikedeng1

A Re-solving Heuristic with Uniformly Bounded Loss for Network Revenue Management 1: Infrequent Re-solving with Thresholding Has Regret O(1), Uniformly in the Horizon T and the Capacities CResearch Paper

Motivation

Network revenue management decides, in real time, which customer requests to accept when every request consumes a bundle of scarce, perishable resources: seats on several flight legs of an itinerary, room-nights of a multi-night hotel stay, bandwidth on several links. The exact dynamic program is intractable for realistic networks, so practice and theory rely on heuristics built from a deterministic linear program (DLP) that replaces random demand by its mean. The question is how much revenue such heuristics lose.

The classical answer (Gallego and van Ryzin 1994, 1997; Talluri and van Ryzin 1998) is that solving the DLP once and following its solution loses O(T)O(\sqrt T)O(T​) over a horizon of length TTT. Re-solving the DLP as capacity is consumed was long believed to help, and Jasin and Kumar (2012) proved that re-solving in every period has bounded loss if the DLP solution is nondegenerate. Bumpensanti and Wang (arXiv:1802.06192v3) showed that the nondegeneracy condition matters: frequent re-solving can lose Ω(T)\Omega(\sqrt T)Ω(T​) on degenerate instances. They then proposed a policy, Infrequent Re-solving with Thresholding (IRT), whose loss is bounded by a constant independent of the horizon and of the capacities, with no nondegeneracy assumption. This mission formalizes that result.

Setting

There are nnn customer classes j∈[n]j \in [n]j∈[n] and mmm resources l∈[m]l \in [m]l∈[m]. Class-jjj customers arrive as independent Poisson processes of rate λj>0\lambda_j > 0λj​>0 on [0,T][0, T][0,T]. Accepting a class-jjj customer earns rj≥0r_j \ge 0rj​≥0 and consumes alj≥0a_{lj} \ge 0alj​≥0 units of resource lll; A=(alj)A = (a_{lj})A=(alj​) is the bill-of-materials matrix with columns AjA_jAj​, and C∈R≥0mC \in \mathbb R^m_{\ge 0}C∈R≥0m​ is the initial capacity. A customer can be accepted only if Aj≤C′A_j \le C'Aj​≤C′ componentwise, C′C'C′ being the remaining capacity; leftover capacity is worthless at TTT.

The DLP with right-hand side bbb is

max⁡x{∑jrjxj ∣ ∑jAjxj≤b, 0≤xj≤λj},\max_x \Big\{ \sum_j r_j x_j \ \Big|\ \sum_j A_j x_j \le b,\ 0 \le x_j \le \lambda_j \Big\},xmax​{j∑​rj​xj​ ​ j∑​Aj​xj​≤b, 0≤xj​≤λj​},

and vDLP(T,C)v^{\mathrm{DLP}}(T, C)vDLP(T,C) is TTT times its value at b=C/Tb = C/Tb=C/T. The hindsight optimum is vHO(T,C)=E[VHO]v^{\mathrm{HO}}(T, C) = \mathbb E[V^{\mathrm{HO}}]vHO(T,C)=E[VHO], where VHOV^{\mathrm{HO}}VHO is the value of the same LP with capacity CCC and the realized demands Λj(T)∼Poisson(λjT)\Lambda_j(T) \sim \mathrm{Poisson}(\lambda_j T)Λj​(T)∼Poisson(λj​T) as upper bounds. It bounds the expected revenue of every non-anticipating policy, and the regret of a policy π\piπ is vHO−vπv^{\mathrm{HO}} - v^\pivHO−vπ.

A probabilistic allocation on a window accepts each class-jjj arrival with a fixed probability pjp_jpj​ when capacity allows. The IRT policy uses τu=T(5/6)u\tau_u = T^{(5/6)^u}τu​=T(5/6)u and re-solving times tu∗=T−τut^*_u = T - \tau_utu∗​=T−τu​ for u=0,…,Ku = 0, \dots, Ku=0,…,K, where K=⌈log⁡log⁡T/log⁡(6/5)⌉K = \lceil \log\log T / \log(6/5)\rceilK=⌈loglogT/log(6/5)⌉. At tu∗t^*_utu∗​ it solves the DLP with right-hand side C(tu∗)/τuC(t^*_u)/\tau_uC(tu∗​)/τu​ (remaining capacity over remaining time), obtaining xux^uxu. In the epochs u<Ku < Ku<K it accepts class jjj with probability 000 if xju<λjτu−1/4x^u_j < \lambda_j\tau_u^{-1/4}xju​<λj​τu−1/4​, else 111 if xju>λj(1−τu−1/4)x^u_j > \lambda_j(1 - \tau_u^{-1/4})xju​>λj​(1−τu−1/4​), else xju/λjx^u_j/\lambda_jxju​/λj​; in the last epoch [tK∗,T][t^*_K, T][tK∗​,T] it uses xjK/λjx^K_j/\lambda_jxjK​/λj​. The family IRTK′\mathrm{IRT}^{K'}IRTK′ re-solves K′K'K′ times on the same schedule; IRT0\mathrm{IRT}^0IRT0 is static probabilistic allocation (SPA), and HOK′\mathrm{HO}^{K'}HOK′ follows IRTK′\mathrm{IRT}^{K'}IRTK′ until tK′∗t^*_{K'}tK′∗​ and then earns the hindsight optimum of what remains.

Formalization targets

Goal: Theorem 1 (p. 16)

There is a constant M=M(λ,r,A)M = M(\lambda, r, A)M=M(λ,r,A) such that for every horizon T∈{1,2,… }T \in \{1, 2, \dots\}T∈{1,2,…} and every capacity vector C≥0C \ge 0C≥0,

vHO(T,C)−vIRT(T,C)≤M.v^{\mathrm{HO}}(T, C) - v^{\mathrm{IRT}}(T, C) \le M .vHO(T,C)−vIRT(T,C)≤M.

The constant is not fixed; the content is its independence of TTT, of CCC, and of the optimal LP solution chosen at each re-solve.

Milestones

  • vHO≤vDLPv^{\mathrm{HO}} \le v^{\mathrm{DLP}}vHO≤vDLP (Sec. 2.2.2, p. 9).
  • Lemma 4 (p. 37): P(∣X−μ∣≥x)≤2e−x2/(3μ)\mathbb P(|X - \mu| \ge x) \le 2e^{-x^2/(3\mu)}P(∣X−μ∣≥x)≤2e−x2/(3μ) for X∼Poisson(μ)X \sim \mathrm{Poisson}(\mu)X∼Poisson(μ), 0<x≤μ0 < x \le \mu0<x≤μ.
  • Proposition 5 (p. 27): vDLP−vSPA≤MTv^{\mathrm{DLP}} - v^{\mathrm{SPA}} \le M\sqrt TvDLP−vSPA≤MT​, uniformly in CCC.
  • Proposition 1 (p. 16): with one re-solve at T−T5/6T - T^{5/6}T−T5/6,
vHO−vHO1≤MTe−κT1/6,vHO−vIRT1≤MTe−κT1/6+MT5/12.v^{\mathrm{HO}} - v^{\mathrm{HO}^1} \le M T e^{-\kappa T^{1/6}},\qquad v^{\mathrm{HO}} - v^{\mathrm{IRT}^1} \le M T e^{-\kappa T^{1/6}} + M T^{5/12}.vHO−vHO1≤MTe−κT1/6,vHO−vIRT1≤MTe−κT1/6+MT5/12.
  • Eq. (13) (p. 29): for every K′K'K′,
vHO−vIRTK′≤M∑u=0K′−1T(5/6)ue−κT(5/6)u/6+MT(5/6)K′/2.v^{\mathrm{HO}} - v^{\mathrm{IRT}^{K'}} \le M\sum_{u=0}^{K'-1} T^{(5/6)^u} e^{-\kappa T^{(5/6)^u/6}} + M T^{(5/6)^{K'}/2}.vHO−vIRTK′≤Mu=0∑K′−1​T(5/6)ue−κT(5/6)u/6+MT(5/6)K′/2.
  • p. 30: T(5/6)K≤eT^{(5/6)^{K}} \le eT(5/6)K≤e, and the right-hand side of (13) at K′=K(T)K' = K(T)K′=K(T) is bounded uniformly in T≥1T \ge 1T≥1.

Significance

The result separates two design choices in re-solving heuristics: how often to re-solve and how to turn an LP solution into a control. It shows that re-solving only O(log⁡log⁡T)O(\log\log T)O(loglogT) times, combined with rounding nearly-degenerate acceptance probabilities to 000 or 111, is enough for bounded regret, and that the constant is uniform over all capacity-to-horizon ratios, so degenerate DLP solutions, which occur only at particular ratios, cause no loss of order. Since vHO≥v∗v^{\mathrm{HO}} \ge v^*vHO≥v∗, it also shows that the hindsight optimum is within a constant of the optimal policy's value in this model.

The result is proved on paper but not machine-checked. A formalization would produce a reusable Poisson-arrival revenue-management model, a precise account of the regret decomposition over re-solving epochs, and an independent check of a proof with known gaps (see Difficulty). No part of it is formalized on the platform.

Difficulty

The obvious argument compares the policy with the DLP solution and controls the deviation of Poisson demand by its standard deviation; this gives only O(T)O(\sqrt T)O(T​), because the DLP value exceeds the hindsight optimum by order T\sqrt TT​ and the lost sales accumulate. Bounded regret needs a comparison with the hindsight LP, path by path: the acceptances made before the last re-solve must remain extendable to a hindsight-optimal solution with overwhelming probability. This requires LP sensitivity of the hindsight solution to the random right-hand side, uniformly over the degenerate cases, and the printed proof's statement of it (Lemma 5) is false as printed: on a two-class, single-resource instance its interval for zˉ1\bar z_1zˉ1​ is empty. A correct argument needs a proximity bound for optimal LP solutions under a change of the demand bounds, summed over all classes, not only over Jλ={j:xj∗=λj}J_\lambda = \{j : x^*_j = \lambda_j\}Jλ​={j:xj∗​=λj​}. The analysis of the schedule (the sum over epochs in (13)) is elementary but the printed chain of inequalities on p. 30 uses a monotonicity that fails for large arguments.

Formalization scope

All objects are in the definition item ResolvingNRM.IRT.Model. Classes are Fin n, resources Fin m, A : Matrix (Fin m) (Fin n) ℝ. LP values reuse piValue from the published RLPBidPrice.Unbiased.Model; the hindsight LP is over real zzz, as in (3). Expectations are explicit sums of Poisson weights. A window of probabilistic allocation is represented by its Poisson count of arrivals, i.i.d. classes with law λj/∑iλi\lambda_j/\sum_i\lambda_iλj​/∑i​λi​ and independent Bernoulli acceptance coins, the standard representation for controls constant on the window. Policies are backward recursions over the re-solving epochs. Ties in the LP are left open: a policy takes an arbitrary optimal-solution selector, and every theorem quantifies over all selectors after its constants.

Conventions and deviations from the page, each disclosed in the item concerned:

  • Standing assumptions of Sec. 2 (p. 7), left implicit there: λj>0\lambda_j > 0λj​>0, r≥0r \ge 0r≥0, A≥0A \ge 0A≥0, C≥0C \ge 0C≥0.
  • "O(g)O(g)O(g)" is read as: there is MMM, depending only on (λ,r,A)(\lambda, r, A)(λ,r,A), with the bound ≤Mg(T)\le M g(T)≤Mg(T) for all T≥1T \ge 1T≥1 (or T>0T > 0T>0) and all C≥0C \ge 0C≥0.
  • Milestones applied to sub-horizons T(5/6)uT^{(5/6)^u}T(5/6)u take a real TTT; the goal takes T∈NT \in \mathbb NT∈N.
  • Lemma 4 is stated for 0<x≤μ0 < x \le \mu0<x≤μ; as printed (all x>0x > 0x>0) it is false, e.g. μ=1\mu = 1μ=1, x=5x = 5x=5.
  • In Proposition 1 and (13), κ>0\kappa > 0κ>0 is existential and uniform in CCC; the printed κ\kappaκ depends on JλJ_\lambdaJλ​ and is derived through Lemma 5.
  • (13) is stated for every number K′K'K′ of re-solves.
  • Algorithm 3 prints the re-solve right-hand side as C(tk∗)/τkC(t^*_k)/\tau_kC(tk∗​)/τk​; it is read with index uuu.
  • For 1≤T≤e1 \le T \le e1≤T≤e the formula for KKK is undefined or negative; the formalization takes K=0K = 0K=0, so IRT is SPA there.
  • The optimal policy value v∗v^*v∗ and the paper's constant α\alphaα are not defined; Lemmas 5 and 6 are not stated.

A statement that lets MMM depend on TTT or CCC, compares IRT with the DLP instead of the hindsight optimum, uses a fixed number of re-solves, or assumes a nondegenerate or vertex LP solution is trivial or a different theorem, and is ruled out by the quantifier order of the goal. The milestone T(5/6)K≤eT^{(5/6)^K} \le eT(5/6)K≤e has a sorry-free local check.

Welcome contributions: LP proximity results (Cook–Gerards–Schrijver–Tardos type), Poisson concentration, superposition and thinning of Poisson processes, and lemmas on expectations of LP values. The source is arXiv:1802.06192v3; its printed page numbers equal the PDF page numbers.

Selected references

  • P. Bumpensanti, H. Wang, A Re-solving Heuristic with Uniformly Bounded Loss for Network Revenue Management, arXiv:1802.06192v3, 2018; Management Science 66(7), 2020. https://arxiv.org/abs/1802.06192
  • S. Jasin, S. Kumar, A Re-Solving Heuristic with Bounded Revenue Loss for Network Revenue Management with Customer Choice, Mathematics of Operations Research 37(2), 2012. https://doi.org/10.1287/moor.1110.0530
  • G. Gallego, G. van Ryzin, A Multiproduct Dynamic Pricing Problem and Its Applications to Network Yield Management, Operations Research 45(1), 1997. https://doi.org/10.1287/opre.45.1.24
  • K. Talluri, G. van Ryzin, An Analysis of Bid-Price Controls for Network Revenue Management, Management Science 44(11), 1998. https://doi.org/10.1287/mnsc.44.11.1577
  • M. Reiman, Q. Wang, An Asymptotically Optimal Policy for a Quantity-Based Network Revenue Management Problem, Mathematics of Operations Research 33(2), 2008. https://doi.org/10.1287/moor.1070.0288
  • W. Cook, A. M. H. Gerards, A. Schrijver, É. Tardos, Sensitivity Theorems in Integer Linear Programming, Mathematical Programming 34, 1986. https://doi.org/10.1007/BF01582230
10 thms1 active userReviewed
Bandit AlgorithmsLinear algebraMachine Learning+1·Captain: mikedeng1

Feature-Based Dynamic Pricing: The EllipsoidPricing Algorithm Has Worst-Case Regret O(R d² ln(T/d))Research Paper

Motivation

Online marketplaces, ad exchanges and real-estate platforms sell items that are each described by a vector of features and are rarely seen twice. A seller who prices such items cannot learn a demand curve item by item. It must learn how features map to values, while pricing, from accept/reject feedback alone. Cohen, Lobel and Paes Leme (Management Science, 2020; SSRN 2737045) model this as feature-based dynamic pricing with adversarial features and a linear valuation. They show that a direct multi-dimensional binary search (PolytopePricing) can have regret exponential in the dimension (their Theorem 1). Their EllipsoidPricing algorithm, a pricing variant of Khachiyan's ellipsoid method (Khachiyan 1979), has worst-case regret quadratic in the dimension and logarithmic in the horizon (their Theorem 2). This mission formalizes Theorem 2 and the chain of lemmas behind it.

Setting

Fix a dimension d≥2d \ge 2d≥2, a radius R>0R > 0R>0 and a horizon TTT. Nature picks an unknown parameter θ∈Rd\theta \in \mathbb{R}^dθ∈Rd with Euclidean norm ∥θ∥≤R\|\theta\| \le R∥θ∥≤R. In each period t=1,…,Tt = 1, \dots, Tt=1,…,T a product arrives with a feature vector xt∈Rdx_t \in \mathbb{R}^dxt​∈Rd, ∥xt∥≤1\|x_t\| \le 1∥xt​∥≤1, and market value θ′xt\theta' x_tθ′xt​. The seller sees xtx_txt​, posts a price ptp_tpt​, and a sale occurs iff pt≤θ′xtp_t \le \theta' x_tpt​≤θ′xt​, earning ptp_tpt​. The regret (Eq. (1)) is

Regret=∑t=1T[θ′xt−pt I{θ′xt≥pt}],\mathrm{Regret} = \sum_{t=1}^{T} \big[\theta' x_t - p_t\,\mathbb{I}\{\theta' x_t \ge p_t\}\big],Regret=t=1∑T​[θ′xt​−pt​I{θ′xt​≥pt​}],

and the worst-case regret of a policy is its maximum over θ\thetaθ and over nature's (possibly adaptive) choice of features.

An ellipsoid with center aaa and positive definite shape matrix AAA is E(A,a)={θ:(θ−a)′A−1(θ−a)≤1}E(A,a) = \{\theta : (\theta - a)' A^{-1} (\theta - a) \le 1\}E(A,a)={θ:(θ−a)′A−1(θ−a)≤1}. EllipsoidPricing with parameter ϵ>0\epsilon > 0ϵ>0 keeps an ellipsoid Et=E(At,at)E_t = E(A_t, a_t)Et​=E(At​,at​), starting from the ball E1=B(0,R)E_1 = B(0,R)E1​=B(0,R). In period ttt it computes

b‾t=xt′at−xt′Atxt,bˉt=xt′at+xt′Atxt,\underline b_t = x_t' a_t - \sqrt{x_t' A_t x_t}, \qquad \bar b_t = x_t' a_t + \sqrt{x_t' A_t x_t},b​t​=xt′​at​−xt′​At​xt​​,bˉt​=xt′​at​+xt′​At​xt​​,

the minimum and maximum of θ^′xt\hat\theta' x_tθ^′xt​ over EtE_tEt​.

  • If bˉt−b‾t≤ϵ\bar b_t - \underline b_t \le \epsilonbˉt​−b​t​≤ϵ, it exploits: it posts pt=b‾tp_t = \underline b_tpt​=b​t​ and keeps Et+1=EtE_{t+1} = E_tEt+1​=Et​.
  • Otherwise it explores: it posts pt=12(bˉt+b‾t)p_t = \tfrac12(\bar b_t + \underline b_t)pt​=21​(bˉt​+b​t​) and replaces EtE_tEt​ by the ellipsoid of Eq. (4) that covers the half of EtE_tEt​ consistent with the feedback. With b=Atxt/xt′Atxtb = A_t x_t / \sqrt{x_t' A_t x_t}b=At​xt​/xt′​At​xt​​, the new shape is A~=d2d2−1(At−2d+1bb′)\tilde A = \frac{d^2}{d^2-1}\big(A_t - \frac{2}{d+1} b b'\big)A~=d2−1d2​(At​−d+12​bb′) and the new center is at±1d+1ba_t \pm \frac{1}{d+1} bat​±d+11​b, with +++ after a sale.

Write λd(A)\lambda_d(A)λd​(A) for the smallest eigenvalue of a symmetric matrix AAA.

Formalization targets

Goal: Theorem 2

There is a universal constant C>0C > 0C>0 such that for all d≥2d \ge 2d≥2, R>0R > 0R>0, T≥2dT \ge 2dT≥2d, ∥θ∥≤R\|\theta\| \le R∥θ∥≤R and ∥xt∥≤1\|x_t\| \le 1∥xt​∥≤1, EllipsoidPricing run with ϵ=Rd2/T\epsilon = R d^2 / Tϵ=Rd2/T satisfies

Regret≤C R d2ln⁡(T/d).\mathrm{Regret} \le C \, R \, d^2 \ln(T/d).Regret≤CRd2ln(T/d).

This is the paper's O(Rd2ln⁡(T/d))O(R d^2 \ln(T/d))O(Rd2ln(T/d)) claim, with no constant fixed.

Milestones, in the order of the paper's argument

  • the closed forms of b‾t,bˉt\underline b_t, \bar b_tb​t​,bˉt​ (§5.2);
  • the containment of each half-ellipsoid in the updated ellipsoid (Eq. (4), already proved on the platform);
  • θ∈Et\theta \in E_tθ∈Et​ and At≻0A_t \succ 0At​≻0 along the run;
  • the volume formula Vol⁡E(A,a)=Vd∏iλi(A)\operatorname{Vol} E(A,a) = V_d \sqrt{\prod_i \lambda_i(A)}VolE(A,a)=Vd​∏i​λi​(A)​ and the volume decrease Vol⁡E(A~)≤e−1/2dVol⁡E(A)\operatorname{Vol} E(\tilde A) \le e^{-1/2d} \operatorname{Vol} E(A)VolE(A~)≤e−1/2dVolE(A) (§5.1);
  • Lemma 2: if z<λd(A)z < \lambda_d(A)z<λd​(A) and det⁡(A−βbb′−zI)≥0\det(A - \beta b b' - zI) \ge 0det(A−βbb′−zI)≥0, then λd(A−βbb′)≥z\lambda_d(A - \beta b b') \ge zλd​(A−βbb′)≥z;
  • Lemma 3: λd(A~)≥d2(d+1)2λd(A)\lambda_d(\tilde A) \ge \frac{d^2}{(d+1)^2} \lambda_d(A)λd​(A~)≥(d+1)2d2​λd​(A);
  • Lemma 4: if λd(A)≤ϵ2/(400d2)\lambda_d(A) \le \epsilon^2/(400 d^2)λd​(A)≤ϵ2/(400d2) and x′Ax>ϵ2/4x' A x > \epsilon^2/4x′Ax>ϵ2/4, then λd(A~)≥λd(A)\lambda_d(\tilde A) \ge \lambda_d(A)λd​(A~)≥λd​(A);
  • the eigenvalue floor λd(At)≥ϵ2/(400(d+1)2)\lambda_d(A_t) \ge \epsilon^2 / (400 (d+1)^2)λd​(At​)≥ϵ2/(400(d+1)2);
  • Lemma 1: at most 2d2ln⁡(20R(d+1)/ϵ)2 d^2 \ln(20 R (d+1)/\epsilon)2d2ln(20R(d+1)/ϵ) exploration periods;
  • an exploitation period sells and has regret at most ϵ\epsilonϵ.

Significance

The result. Theorem 2 shows that contextual posted-price learning with adversarial features costs only polynomially many "probing" sales in the dimension. A direct cutting-plane search over the polytope of consistent parameters (PolytopePricing) cannot do this: its worst-case regret is exponential in ddd (Theorem 1). The ellipsoid analysis via the smallest eigenvalue (Lemmas 2–4) is reused in the paper's noisy-valuation extension (ShallowPricing, §6) and in later work on contextual pricing and contextual search.

Formalizing it. No part of the paper is machine-checked. The Eq. (4) containment is already proved on the platform (Bertsimas–Tsitsiklis Theorem 8.1). The volume factor there is e−1/(2(d+1))e^{-1/(2(d+1))}e−1/(2(d+1)), weaker than the e−1/2de^{-1/2d}e−1/2d used here. The eigenvalue lemmas on rank-one perturbations (Lemmas 2–4) are new to the platform and reusable for any analysis of ellipsoid-type updates. A formal proof of the goal would also settle a gap in the printed proof, noted under the formalization scope.

Difficulty

The obvious argument is the textbook ellipsoid bound: volume shrinks by e−1/2de^{-1/2d}e−1/2d per exploration, so exploration must stop. This fails because the ellipsoid need not become small in every direction. The volume can go to zero while one axis stays long, and then a feature along that axis keeps triggering exploration. The paper replaces the volume argument by a lower bound on the smallest eigenvalue, which needs a sign analysis of the characteristic polynomial of a rank-one perturbation (Lemmas 2–4) together with the exploration test xt′Atxt>ϵ2/4x_t' A_t x_t > \epsilon^2/4xt′​At​xt​>ϵ2/4. The constant k=1/(400d2)k = 1/(400 d^2)k=1/(400d2) in Lemma 4 is chosen to make a specific inequality hold for every d≥2d \ge 2d≥2.

The regret bound needs, in addition, a per-round bound for exploration periods. The printed proof (p. 16) uses "the trivial bound of regret RRR per round". That bound is not immediate when an exploration price is accepted, since the center ata_tat​ can leave B(0,R)B(0,R)B(0,R) and the round's regret is then only bounded by xt′Atxt\sqrt{x_t' A_t x_t}xt′​At​xt​​. This is why the goal is the OOO-form with an unspecified universal constant rather than the proof's explicit inequality.

Formalization scope

  • Types. Vectors are Fin d → ℝ; θ′x\theta' xθ′x is the dot product θ ⬝ᵥ x. Both norms are Euclidean, encoded as x⋅x≤1x\cdot x \le 1x⋅x≤1 and θ⋅θ≤R2\theta\cdot\theta \le R^2θ⋅θ≤R2; Mathlib's ‖·‖ on Fin d → ℝ is the sup norm and is never used. Ellipsoids, the ball, and the update (4) are the published LinearOptimization.ellipsoid, ellipsoidBall, ellipsoidUpdateCenter and ellipsoidUpdateMatrix.
  • Eigenvalues and volume. λd(A)\lambda_d(A)λd​(A) is the infimum of the real spectrum, i.e. the minimum eigenvalue of a symmetric matrix. Volumes are Lebesgue measure on Fin d → ℝ in [0,∞][0, \infty][0,∞]. φD(z)\varphi_D(z)φD​(z) is det⁡(D−zI)\det(D - zI)det(D−zI) as printed, not Mathlib's charpoly.
  • The algorithm. EllipsoidPricing is a deterministic recursion indexed from 000 (Lean period ttt is the paper's t+1t+1t+1), started from E1=B(0,R)E_1 = B(0,R)E1​=B(0,R), which the paper allows ("or in fact any ellipsoid that contains K1K_1K1​", p. 11). The paper's "smallest ellipsoid containing Ht+1H_{t+1}Ht+1​" is replaced by its closed form (4), which the paper gives on p. 14. Because the algorithm is deterministic, a closed-loop nature generates a fixed feature sequence, and the worst case is a universal quantifier over θ\thetaθ in the ball and over sequences.
  • Added hypotheses, all disclosed in the items.
    • d≥2d \ge 2d≥2 everywhere, since (4) divides by d2−1d^2 - 1d2−1.
    • T≥2dT \ge 2dT≥2d in the goal, since ln⁡(T/d)≤0\ln(T/d) \le 0ln(T/d)≤0 for T≤dT \le dT≤d.
    • ϵ≤20R(d+1)\epsilon \le 20 R (d+1)ϵ≤20R(d+1) in Lemma 1 and the eigenvalue floor: for larger ϵ\epsilonϵ the printed bound is negative.
    • ∥x∥≤1\|x\| \le 1∥x∥≤1 in Lemma 4, the §3 normalization that its proof uses.
  • Not taken from the page. E₁ is not the Löwner–John ellipsoid of a general K1K_1K1​: the proof's "E~1\tilde E_1E~1​ lies in the ball of radius RRR" fails for it. The proof's explicit inequality Regret≤NR+(T−N)ϵ\mathrm{Regret} \le NR + (T-N)\epsilonRegret≤NR+(T−N)ϵ is not a milestone, because of the per-round gap above.
  • Ruled out. The goal may not be trivialized by any of the following:
    • a constant CCC that depends on ddd, RRR, TTT or the instance;
    • regret measured against a price the algorithm never posts, or an exploit price that is not accepted;
    • features bounded in the sup norm;
    • an explore test or update different from (4);
    • a single fixed feature sequence.
  • Contributions welcome. Proofs of the volume formula and the e−1/2de^{-1/2d}e−1/2d decrease, which are reusable for any ellipsoid-method development; the eigenvalue lemmas; and the per-round bound for accepted exploration prices that the goal needs.

Selected references

  • M. C. Cohen, I. Lobel, R. Paes Leme, Feature-Based Dynamic Pricing, Management Science 66(11), 2020. https://doi.org/10.1287/mnsc.2019.3485 (authors' copy: https://ssrn.com/abstract=2737045)
  • L. G. Khachiyan, Polynomial algorithms in linear programming, USSR Comput. Math. Math. Phys. 20(1), 1980. https://doi.org/10.1016/0041-5553(80)90061-0
  • M. Grötschel, L. Lovász, A. Schrijver, Geometric Algorithms and Combinatorial Optimization, Springer, 1993. https://doi.org/10.1007/978-3-642-78240-4
  • G. H. Golub, Some modified matrix eigenvalue problems, SIAM Review 15(2), 1973. https://doi.org/10.1137/1015032
  • D. Bertsimas, J. N. Tsitsiklis, Introduction to Linear Optimization, Athena Scientific, 1997 (Theorem 8.1, the platform's ellipsoid update).
16 thms3 active usersReviewed
Operations ResearchProbability·Captain: mikedeng1

A Re-solving Heuristic with Uniformly Bounded Loss for Network Revenue Management 3: On a Degenerate Two-Class Instance, Frequent Re-solving Has Regret at Least Ω(√T)Research Paper

Motivation

Network revenue managers sell access to limited capacity over time. An airline, for example, may have several products that use the same seat inventory, with some products earning more revenue than others. The decision to accept a request must be made when it arrives, before future requests are known. A standard way to guide that decision is to solve a deterministic linear program using expected demand, then use its optimal allocation as a probability of acceptance. As inventory changes, one can solve the program again. Bumpensanti and Wang study how often to do so, and show that frequent re-solving has a real limitation when the linear program is degenerate: on a particular instance its expected loss grows at least as the square root of the selling horizon. Bumpensanti and Wang, arXiv:1802.06192v3, Sections 3.1 and 5.1.

The same paper gives an infrequent re-solving policy with uniformly bounded regret and an upper bound of order T\sqrt TT​ for the frequent re-solving policy. The lower-bound result is the counterpart that identifies a case where that order cannot be removed by the frequent policy itself. It uses a two-class example small enough to expose the issue without other network complications. Bumpensanti and Wang, arXiv:1802.06192v3, Theorem 1 and Propositions 2–3.

Setting

There is one resource, initially with TTT units of capacity, and two customer classes. Class jjj arrives through an independent Poisson process of rate one. Accepting a customer uses one unit of capacity and earns price rjr_jrj​, where 0<r2<r10<r_2<r_10<r2​<r1​. Rejected requests earn nothing, and unused capacity has no terminal value. The higher-price class has priority in the deterministic allocation, but its realized demand is random. The horizon TTT is a positive integer, divided into TTT periods of length one. Bumpensanti and Wang, arXiv:1802.06192v3, Section 2, pp. 7–9, and Appendix C.1, p. 32.

At the start of a period with kkk periods left and capacity ccc, the deterministic linear program (DLP) chooses rates x1,x2x_1,x_2x1​,x2​ that maximize r1x1+r2x2r_1x_1+r_2x_2r1​x1​+r2​x2​ subject to x1+x2≤c/kx_1+x_2\le c/kx1​+x2​≤c/k and 0≤xj≤10\le x_j\le10≤xj​≤1. In this instance its first coordinate is x1=min⁡{c/k,1}x_1=\min\{c/k,1\}x1​=min{c/k,1}. The frequent re-solving policy (FR) accepts each class-jjj request in that period with probability xjx_jxj​, provided capacity remains. It re-solves the DLP at the next period with the new capacity. This is Algorithm 2 specialized to the example. Bumpensanti and Wang, arXiv:1802.06192v3, Algorithm 2, p. 11, and Appendix D, p. 42.

The hindsight optimum knows the total demand from both classes at the end of the horizon and chooses the best feasible allocation using that information. Its expected value is vHO(T,T)v^{\mathrm{HO}}(T,T)vHO(T,T); the first argument is horizon length and the second is initial capacity. Let vFR(T,T)v^{\mathrm{FR}}(T,T)vFR(T,T) be the frequent policy's expected revenue. Since hindsight has more information, their difference is a regret benchmark. The formulation uses the paper's hindsight linear program, which is exact on this unit-consumption instance. Bumpensanti and Wang, arXiv:1802.06192v3, Eq. (3) and Definition 1, p. 9.

Formalization targets

Main result

For every pair of prices 0<r2<r10<r_2<r_10<r2​<r1​, there are M>0M>0M>0 and T0≥1T_0\ge1T0​≥1 such that, for all integer T≥T0T\ge T_0T≥T0​ and every optimal DLP selector used by FR,

vHO(T,T)−vFR(T,T)≥MT.v^{\mathrm{HO}}(T,T)-v^{\mathrm{FR}}(T,T)\ge M\sqrt T.vHO(T,T)−vFR(T,T)≥MT​.

This specializes Proposition 2's existence statement to the explicit family in its Appendix C.1 proof. The constant may depend on the prices, but is chosen before the selector and the horizon. The benchmark is the hindsight value, as in the proposition. Bumpensanti and Wang, arXiv:1802.06192v3, Proposition 2, p. 18, and Appendix C.1, pp. 32–34.

Supporting results

The milestones record the optimal DLP allocation on the example and two clauses of Lemma 7. The latter estimate the probabilities that the high-price arrival count is moderately below its mean in the first third and moderately above its mean in the last third. With T′T'T′ an integer phase length and N∼Poisson⁡(T′)N\sim\operatorname{Poisson}(T')N∼Poisson(T′), the first bound is

P(T′−4T′≤N≤T′−3T′)≥0.0013−0.9496/T′.\mathbb P(T'-4\sqrt{T'}\le N\le T'-3\sqrt{T'}) \ge 0.0013-0.9496/\sqrt{T'}.P(T′−4T′​≤N≤T′−3T′​)≥0.0013−0.9496/T′​.

The third-phase bound uses the standard normal cumulative distribution function Φ\PhiΦ and the unrounded constant Φ(7)−Φ(6)\Phi(7)-\Phi(6)Φ(7)−Φ(6). Bumpensanti and Wang, arXiv:1802.06192v3, Lemma 7 and Eqs. (49)–(53), p. 40.

Significance

This result shows that solving the DLP every period does not, by itself, give a horizon-independent expected loss. The example has only one resource and two customer classes; the gap cannot be attributed to a large network. The paper's separate bounded-regret result uses a different re-solving schedule and acceptance rule, so formalizing this lower bound helps distinguish guarantees for the two policies. Together with the upper bound for FR, it identifies the square-root order as the relevant scale for this policy on the example. Bumpensanti and Wang, arXiv:1802.06192v3, Sections 4–5.

The mathematical proposition is proved in the cited paper; the Lean goal here remains an open theorem statement. Formalizing it requires precise interfaces for a Poisson arrival model, a capacity-limited randomized policy, a hindsight LP benchmark, and asymptotic lower bounds. Those interfaces can be reused in later revenue-management results, while the two-class instance gives a concrete check on the conventions.

Difficulty

The initial DLP solution allocates the average capacity to the high-price class. A simple intuition might therefore predict that repeated re-solving preserves capacity for that class. The remaining-capacity ratio, however, responds to realized arrivals; when the ratio rises above one, the re-solved LP assigns positive acceptance probability to the lower-price class. The policy then spends capacity before later high-price arrivals are known. The lower-bound analysis needs a path event that controls arrivals through the middle phase and a joint event involving the policy's accepted requests. Marginal Poisson counts alone do not determine those admissions. Bumpensanti and Wang, arXiv:1802.06192v3, Appendix C.1, pp. 32–34.

Formalization scope

Lean uses Fin 2 for classes and Fin 1 for the resource. The general model takes positive Poisson rates, nonnegative prices and nonnegative consumption. On the lower-bound instance the rates and consumptions equal one, capacity equals TTT, and prices satisfy 0<r2<r10<r_2<r_10<r2​<r1​. These are the operational standing assumptions of the paper's Section 2 and the positive-price reading used by the example. Optimal DLP ties are represented by quantifying over every selector. The LP value is the published RLPBidPrice.Unbiased.piValue definition, evaluated on nonnegative capacity vectors; the rate and capacity conventions make its real supremum well posed.

The window law superposes the independent class Poisson processes into a Poisson total count, independent class labels, and independent Bernoulli acceptance marks. A fold over ordered arrivals applies Algorithm 2's capacity check after every acceptance. The expected hindsight value sums over independent total demand counts, and the FR value is a backward recursion over the remaining unit periods. The target's MMM is strictly positive and fixed before the horizon, so neither a zero constant nor a horizon-dependent constant can satisfy it trivially.

The two Lemma 7 milestones cover only Q1Q_1Q1​ and Q3Q_3Q3​. Its continuous-time Q2Q_2Q2​ event and the joint admission event require a path probability space and are outside this draft. Equation (52) yields Φ(7)−Φ(6)\Phi(7)-\Phi(6)Φ(7)−Φ(6); the printed 9.8531×10−109.8531\times10^{-10}9.8531×10−10 rounds that value upward. The source proof divides TTT into three exact integer-length phases, while Proposition 2 is stated for all sufficiently large TTT; this gap requires attention in a complete proof. The source here is arXiv:1802.06192v3, whose printed and PDF page numbers coincide.

Selected references

  • Bumpensanti, P. and Wang, H., A Re-solving Heuristic with Uniformly Bounded Loss for Network Revenue Management, arXiv preprint arXiv:1802.06192v3, 2018. PDF.
7 thms1 active userReviewed
Bandit AlgorithmsMachine LearningOperations Research+1·Captain: mikedeng1

MNL-Bandit: A Dynamic Learning Approach to Assortment Selection I: The Epoch-Based UCB Policy Has Regret at Most C₁√(NT log NT) + C₂N log² NT Under the No-Purchase AssumptionResearch Paper

Motivation

A retailer that decides which products to display, an online platform that decides which items to show in a recommendation slot, and an airline that decides which fare classes to open all face the same problem: the set of options offered changes what customers buy, and the substitution pattern is not known in advance. The multinomial logit (MNL) model is the standard model of this substitution in revenue management (Talluri and van Ryzin 2004), and optimizing the offered set under a known MNL model is a classical, efficiently solvable problem (e.g. Rusmevichientong, Shen and Shmoys 2010, for a capacity constraint).

The MNL-Bandit asks what happens when the MNL parameters are unknown and must be learned from the purchases themselves, while revenue is being earned. Earlier approaches (Rusmevichientong, Shen and Shmoys 2010; Sauré and Zeevi 2013) separate exploration from exploitation and need the gap between the best and the second-best assortment to be known. Agrawal, Avadhanula, Goyal and Zeevi (Oper. Res. 2019, arXiv:1706.03880) give a single upper-confidence-bound policy whose regret is of order NT\sqrt{NT}NT​ up to logarithmic factors, with no such knowledge. This mission formalizes that result, Theorem 1 of the paper.

Setting

There are NNN products with known revenues ri∈[0,1]r_i\in[0,1]ri​∈[0,1]. At each time t=1,…,Tt=1,\dots,Tt=1,…,T the seller offers an assortment StS_tSt​ from a family S\mathcal SS of feasible subsets of {1,…,N}\{1,\dots,N\}{1,…,N}, and one customer either buys a product ct∈Stc_t\in S_tct​∈St​ or buys nothing (ct=0c_t=0ct​=0). Given St=SS_t=SSt​=S, the choice follows the MNL model with attraction parameters v1,…,vN≥0v_1,\dots,v_N\ge0v1​,…,vN​≥0 and v0=1v_0=1v0​=1:

pi(S)=vi1+∑j∈Svj(i∈S∪{0}),pi(S)=0 otherwise,p_i(S)=\frac{v_i}{1+\sum_{j\in S}v_j}\quad(i\in S\cup\{0\}),\qquad p_i(S)=0\ \text{otherwise},pi​(S)=1+∑j∈S​vj​vi​​(i∈S∪{0}),pi​(S)=0 otherwise,

independently of the past. The expected revenue of SSS is R(S,v)=∑i∈Srivi/(1+∑j∈Svj)R(S,v)=\sum_{i\in S}r_iv_i\big/\big(1+\sum_{j\in S}v_j\big)R(S,v)=∑i∈S​ri​vi​/(1+∑j∈S​vj​). The family S\mathcal SS is described by totally unimodular constraints, S={S:A x(S)≤b}\mathcal S=\{S : A\,x(S)\le b\}S={S:Ax(S)≤b} with x(S)x(S)x(S) the incidence vector of SSS, AAA totally unimodular and bbb integral (cardinality constraints are the main example). A policy chooses StS_tSt​ from the past choices, and its regret is

Regπ(T,v)=T R(S∗,v)−Eπ[∑t=1TR(St,v)],R(S∗,v)=max⁡S∈SR(S,v).\mathrm{Reg}_\pi(T,v)=T\,R(S^*,v)-\mathbb E_\pi\Big[\sum_{t=1}^TR(S_t,v)\Big],\qquad R(S^*,v)=\max_{S\in\mathcal S}R(S,v).Regπ​(T,v)=TR(S∗,v)−Eπ​[t=1∑T​R(St​,v)],R(S∗,v)=S∈Smax​R(S,v).

The seller knows rrr and S\mathcal SS but not vvv.

Algorithm 1 works in epochs: epoch ℓ\ellℓ offers one assortment SℓS_\ellSℓ​ repeatedly until a customer buys nothing. Let v^i,ℓ\hat v_{i,\ell}v^i,ℓ​ be the number of purchases of iii in epoch ℓ\ellℓ, Ti(ℓ)T_i(\ell)Ti​(ℓ) the number of the first ℓ\ellℓ epochs that offered iii, and vˉi,ℓ\bar v_{i,\ell}vˉi,ℓ​ the average of v^i,τ\hat v_{i,\tau}v^i,τ​ over those epochs. The algorithm forms the upper confidence bounds

vi,ℓUCB=vˉi,ℓ+vˉi,ℓ 48log⁡(Nℓ+1)Ti(ℓ)+48log⁡(Nℓ+1)Ti(ℓ),v^{\mathrm{UCB}}_{i,\ell}=\bar v_{i,\ell}+\sqrt{\bar v_{i,\ell}\,\frac{48\log(\sqrt N\ell+1)}{T_i(\ell)}}+\frac{48\log(\sqrt N\ell+1)}{T_i(\ell)},vi,ℓUCB​=vˉi,ℓ​+vˉi,ℓ​Ti​(ℓ)48log(N​ℓ+1)​​+Ti​(ℓ)48log(N​ℓ+1)​,

starting from vi,0UCB=1v^{\mathrm{UCB}}_{i,0}=1vi,0UCB​=1, and offers next the assortment in S\mathcal SS that maximizes R(S,v⋅,ℓUCB)R(S,v^{\mathrm{UCB}}_{\cdot,\ell})R(S,v⋅,ℓUCB​).

Formalization targets

Goal: Theorem 1 (p. 11)

Under Assumption 4.1 (every vi≤v0=1v_i\le v_0=1vi​≤v0​=1, and S\mathcal SS is closed under taking subsets) there are absolute constants C1,C2C_1,C_2C1​,C2​ such that for every instance and every horizon TTT,

Regπ(T,v)≤C1NTlog⁡NT+C2Nlog⁡2NT.\mathrm{Reg}_\pi(T,v)\le C_1\sqrt{NT\log NT}+C_2N\log^2NT .Regπ​(T,v)≤C1​NTlogNT​+C2​Nlog2NT.

The constants are not fixed: the goal asserts the order of the regret, which is what the paper claims.

Milestones

The milestones follow the proof in Appendix A of the paper:

  1. The law of an epoch's purchase count: Lemma A.1, its moment generating function given the offered assortment.
  2. Concentration: Theorem 5, Chernoff bounds for geometric variables; Corollary D.1; and Lemma A.2 for the averages vˉi,ℓ\bar v_{i,\ell}vˉi,ℓ​ along Algorithm 1.
  3. Lemma 4.1: vi,ℓUCB≥viv^{\mathrm{UCB}}_{i,\ell}\ge v_ivi,ℓUCB​≥vi​ with probability at least 1−6/(Nℓ)1-6/(N\ell)1−6/(Nℓ), and its rate of convergence.
  4. Optimism of the revenue estimate: Lemma A.3 (monotonicity of the optimal revenue in vvv) and Lemma 4.2.
  5. The per-epoch error: Lemma A.4 (a Lipschitz bound) and Lemma 4.3.
  6. The epoch decomposition of the regret, (A.14).

Significance

Theorem 1 shows that learning the choice model costs only O~(NT)\tilde O(\sqrt{NT})O~(NT​) revenue, with no dependence on how separated the instance is. The paper also proves a lower bound of order NT/K\sqrt{NT/K}NT/K​ under a KKK-cardinality constraint (its Theorem 2), so for small KKK the bound is optimal up to logarithmic factors. The epoch device used here reduces the analysis of a choice model with substitution to unbiased i.i.d. estimates of each viv_ivi​.

The result is proved in the paper (Appendix A, with the concentration bounds in Appendix D). It has not been machine-checked. A formal proof would also settle several printed slips that the formalization had to repair (listed under Formalization scope). The concentration bounds for geometric variables, Theorem 5 and Corollary D.1, are useful outside this mission.

Difficulty

The estimates v^i,ℓ\hat v_{i,\ell}v^i,ℓ​ are counted in epochs whose assortments are chosen adaptively from earlier estimates, so they are not independent a priori. The paper's argument that they are i.i.d. geometric with mean viv_ivi​ has to be made rigorous for an adaptive policy. The variables are unbounded, so the standard Chernoff–Hoeffding bounds for bounded variables do not apply, and the number of samples Ti(ℓ)T_i(\ell)Ti​(ℓ) is random, which requires a union bound over all its possible values. Finally, the regret is a sum over customers while the analysis is per epoch, and the two are linked through the expected epoch length 1+∑j∈Sℓvj1+\sum_{j\in S_\ell}v_j1+∑j∈Sℓ​​vj​ under a random number of epochs and a horizon that cuts the last one.

Formalization scope

  • Model. Products are Fin N; a choice is none (no purchase) or some i. The parameter v0v_0v0​ is normalized to 111, as the paper allows. Algorithm 1 is deterministic, so a history of horizon TTT is a sequence of TTT choices with the product law ∏tpct(St)\prod_tp_{c_t}(S_t)∏t​pct​​(St​), and every expectation is a finite sum. The revenue R(S,v)R(S,v)R(S,v) is the published definition ChoiceCDLP.MNL.mnlObjective v r 1 S.
  • Feasible family. S\mathcal SS is a finite family with the TU representation (2.3), closed under subsets, and nonempty (the paper presupposes nonemptiness when it writes S∗=arg⁡max⁡S^*=\arg\maxS∗=argmax).
  • Algorithm. Algorithm 1 is a policy defined by recursion on the history. It receives rrr and an argmax selector for S\mathcal SS, never vvv. Every statement quantifies over all selectors, that is, over every tie-breaking rule. While a product has never been offered, its vUCBv^{\mathrm{UCB}}vUCB stays at the initial value 111.
  • Epoch events. The lemmas about epoch ℓ\ellℓ are stated for every finite horizon TTT, on the event that epoch ℓ\ellℓ is completed by customer TTT; uniformity in TTT is the paper's infinite-horizon statement.
  • Repairs of the page, all disclosed in the items:
    1. Lemma A.2's misprinted log⁡(ℓ+1)\log(\ell+1)log(ℓ+1) is replaced by log⁡(Nℓ+1)\log(\sqrt N\ell+1)log(N​ℓ+1).
    2. Lemmas 4.2 and 4.3 are indexed by the estimate built at the end of epoch ℓ\ellℓ.
    3. Lemma 4.3's free index iii becomes the sum over i∈Sℓ+1i\in S_{\ell+1}i∈Sℓ+1​.
    4. Lemmas 4.1 and 4.3 use the explicit constants C1=72+24C_1=\sqrt{72}+\sqrt{24}C1​=72​+24​ and C2=144C_2=144C2​=144 of the proof.
    5. Lemma A.3 and Lemma 4.2 require positive parameters on the optimal set; the printed versions are false otherwise.
    6. (A.14) is an inequality, since the horizon cuts the last epoch.
  • Not stated. Corollary A.1 (the epoch estimates are i.i.d. geometric across the epochs of the adaptive policy) is not posed as an item: a single-epoch law would be a weaker statement than the page's.
  • No trivialization. The policy cannot see vvv, the constants of the goal precede every other quantifier, and the regret is the paper's (2.6), so the goal cannot be met by a policy that offers S∗S^*S∗ or by constants that depend on the instance.
  • Welcome contributions. Infrastructure for adaptive sampling (the i.i.d. property of the epoch estimates), Chernoff bounds for geometric and sub-exponential variables, and a Wald-type identity for epochs.

The paper's own proofs are in its Appendices A and D.

Selected references

  • S. Agrawal, V. Avadhanula, V. Goyal, A. Zeevi, MNL-Bandit: A Dynamic Learning Approach to Assortment Selection, Operations Research 67(5):1453–1485, 2019. arXiv:1706.03880v2, doi:10.1287/opre.2018.1832
  • P. Rusmevichientong, Z.-J. M. Shen, D. B. Shmoys, Dynamic assortment optimization with a multinomial logit choice model and capacity constraint, Operations Research 58(6):1666–1680, 2010. doi:10.1287/opre.1100.0866
  • D. Sauré, A. Zeevi, Optimal dynamic assortment planning with demand learning, Manufacturing & Service Operations Management 15(3):387–404, 2013. doi:10.1287/msom.2013.0429
  • K. Talluri, G. van Ryzin, Revenue management under a general discrete choice model of consumer behavior, Management Science 50(1):15–33, 2004. doi:10.1287/mnsc.1030.0147
  • M. Mitzenmacher, E. Upfal, Probability and Computing, Cambridge University Press, 2005.
15 thms1 active userReviewed
Operations ResearchOptimizationTheoretical Computer Science·Captain: mikedeng1

Assortment Optimisation Under a General Discrete Choice Model: A Tight Analysis of Revenue-Ordered Assortments II: Revenue-Ordered Assortments Earn OPT/(1 + ln ν), ν the Optimum's Purchase RatioResearch Paper

Motivation

A retailer that can display only some of its products must choose an assortment: the set of products offered to an arriving customer. Customers substitute: whether a given product is bought depends on what else is on the shelf. The assortment problem asks for the offer set that maximises expected revenue under a model of this substitution behaviour. It is a core problem of revenue management (Talluri and van Ryzin, Management Science 2004), and it is NP-hard already for a mixture of two multinomial logit models (Rusmevichientong, Shmoys, Tong and Topaloglu, POMS 2014).

The heuristic used most widely in practice is revenue-ordered assortments: sort the products by price and offer, for some threshold, every product priced at least that threshold. It is optimal under the multinomial logit model (Talluri and van Ryzin 2004) but not in general. Berbeglia and Joret (arXiv:1606.01371v3, 2019; Algorithmica 2020) give a tight analysis of its approximation ratio for every regular discrete choice model, a class that includes all random utility models. They prove three incomparable guarantees. This mission formalizes the third, Theorem 3.3, whose ratio depends on the purchase behaviour of an optimal assortment rather than on the prices. Companion missions in this series cover the price-ratio bound (Theorem 3.2) and the tightness of all three bounds (Theorem 3.4).

Setting

Let C\mathcal CC be a finite nonempty set of products. A system of choice probabilities gives, for every offer set S⊆CS\subseteq\mathcal CS⊆C and product xxx, the probability P(x,S)\mathcal P(x,S)P(x,S) that a customer offered SSS buys xxx. Buying nothing is the option 000, with P(0,S)=1−∑x∈SP(x,S)\mathcal P(0,S)=1-\sum_{x\in S}\mathcal P(x,S)P(0,S)=1−∑x∈S​P(x,S). The model is regular when

  1. P(x,S)≥0\mathcal P(x,S)\ge0P(x,S)≥0 for products and for x=0x=0x=0;
  2. P(x,S)=0\mathcal P(x,S)=0P(x,S)=0 for x∉Sx\notin Sx∈/S;
  3. ∑x∈SP(x,S)≤1\sum_{x\in S}\mathcal P(x,S)\le1∑x∈S​P(x,S)≤1;
  4. P(x,S)≥P(x,S′)\mathcal P(x,S)\ge\mathcal P(x,S')P(x,S)≥P(x,S′) whenever S⊆S′S\subseteq S'S⊆S′ and x∈S∪{0}x\in S\cup\{0\}x∈S∪{0}.

Axiom 4 at x=0x=0x=0 says that enlarging the offer set never makes buying nothing more likely.

Prices are a function r:C→R>0r:\mathcal C\to\mathbb R_{>0}r:C→R>0​. Offering SSS earns rev(S)=∑x∈SP(x,S) r(x)\mathrm{rev}(S)=\sum_{x\in S}\mathcal P(x,S)\,r(x)rev(S)=∑x∈S​P(x,S)r(x), and OPT=max⁡S⊆Crev(S)\mathrm{OPT}=\max_{S\subseteq\mathcal C}\mathrm{rev}(S)OPT=maxS⊆C​rev(S). Let r1<⋯<rkr_1<\dots<r_kr1​<⋯<rk​ be the distinct values of rrr, and let Si={x:r(x)≥ri}S_i=\{x:r(x)\ge r_i\}Si​={x:r(x)≥ri​} for i∈[k]i\in[k]i∈[k]. The revenue-ordered strategy earns

RO=max⁡1≤i≤krev(Si).\mathrm{RO}=\max_{1\le i\le k}\mathrm{rev}(S_i).RO=1≤i≤kmax​rev(Si​).

For an assortment S∗S^*S∗, the purchase profile is

Ni=∑x∈S∗, r(x)≥riP(x,S∗)(i∈[k]),Nk+1:=0,N_i=\sum_{x\in S^*,\ r(x)\ge r_i}\mathcal P(x,S^*)\qquad(i\in[k]),\qquad N_{k+1}:=0,Ni​=x∈S∗, r(x)≥ri​∑​P(x,S∗)(i∈[k]),Nk+1​:=0,

the probability that a customer offered S∗S^*S∗ buys something priced at least rir_iri​. It is non-increasing in iii.

Formalization targets

Goal: Theorem 3.3 (p. 9)

Let S∗S^*S∗ be optimal, rev(S∗)=OPT\mathrm{rev}(S^*)=\mathrm{OPT}rev(S∗)=OPT, suppose N1>0N_1>0N1​>0, and let ℓ∈[k]\ell\in[k]ℓ∈[k] be maximum with Nℓ>0N_\ell>0Nℓ​>0. Then

OPT≤(∑i=1ℓNi−Ni+1Ni)ROand∑i=1ℓNi−Ni+1Ni≤1+ln⁡ν,ν=N1Nℓ.\mathrm{OPT}\le\Big(\sum_{i=1}^{\ell}\frac{N_i-N_{i+1}}{N_i}\Big)\mathrm{RO} \qquad\text{and}\qquad \sum_{i=1}^{\ell}\frac{N_i-N_{i+1}}{N_i}\le1+\ln\nu,\quad\nu=\frac{N_1}{N_\ell}.OPT≤(i=1∑ℓ​Ni​Ni​−Ni+1​​)ROandi=1∑ℓ​Ni​Ni​−Ni+1​​≤1+lnν,ν=Nℓ​N1​​.

The first inequality is the paper's sum-form factor, the bound that Theorem 3.4 shows to be tight. The second is the closed form 1/(1+ln⁡ν)1/(1+\ln\nu)1/(1+lnν).

Milestones (in proof order)

  • Lemma 2.1 (p. 6): ∑x∈SP(x,S)≤∑x∈S′P(x,S′)\sum_{x\in S}\mathcal P(x,S)\le\sum_{x\in S'}\mathcal P(x,S')∑x∈S​P(x,S)≤∑x∈S′​P(x,S′) for S⊆S′S\subseteq S'S⊆S′.
  • First observation of the proof (p. 9): Ni≤∑x∈SiP(x,Si)N_i\le\sum_{x\in S_i}\mathcal P(x,S_i)Ni​≤∑x∈Si​​P(x,Si​) and Niri≤∑x∈SiP(x,Si)riN_ir_i\le\sum_{x\in S_i}\mathcal P(x,S_i)r_iNi​ri​≤∑x∈Si​​P(x,Si​)ri​.
  • Inequality (6) (p. 9): Niri≤RON_ir_i\le\mathrm{RO}Ni​ri​≤RO for every i∈[k]i\in[k]i∈[k].
  • Revenue identity (p. 9): rev(S∗)=∑i=1ℓ(Ni−Ni+1)ri=∑i=1ℓNi−Ni+1NiNiri\mathrm{rev}(S^*)=\sum_{i=1}^{\ell}(N_i-N_{i+1})r_i=\sum_{i=1}^{\ell}\frac{N_i-N_{i+1}}{N_i}N_ir_irev(S∗)=∑i=1ℓ​(Ni​−Ni+1​)ri​=∑i=1ℓ​Ni​Ni​−Ni+1​​Ni​ri​.
  • Logarithmic step (p. 9, the comparison 1/∑≥1/(1+ln⁡ν)1/\sum\ge1/(1+\ln\nu)1/∑≥1/(1+lnν) in Theorem 3.3, used in the last inequality of the proof): for N1≥⋯≥Nℓ>0N_1\ge\dots\ge N_\ell>0N1​≥⋯≥Nℓ​>0 and Nℓ+1=0N_{\ell+1}=0Nℓ+1​=0, ∑i=1ℓ(Ni−Ni+1)/Ni≤1+ln⁡(N1/Nℓ)\sum_{i=1}^{\ell}(N_i-N_{i+1})/N_i\le1+\ln(N_1/N_\ell)∑i=1ℓ​(Ni​−Ni+1​)/Ni​≤1+ln(N1​/Nℓ​). The paper states this step without proof.

Significance

Theorem 3.3 gives a guarantee that is independent of prices. The bound 1+ln⁡(rk/r1)1+\ln(r_k/r_1)1+ln(rk​/r1​) of Theorem 3.2 grows when prices are spread out. The bound 1+ln⁡ν1+\ln\nu1+lnν is small whenever an optimal assortment sells its expensive products with probability comparable to its overall purchase probability. In §4 the paper combines it with a reduction from unit-demand pricing to derive, for example, the 1/(1+ln⁡m)1/(1+\ln m)1/(1+lnm) guarantee of uniform pricing for the unit-demand min-pricing problem (Corollary 4.8, originally due to Aggarwal et al.). Together with Theorems 3.1 and 3.2, it describes the performance of the most common assortment heuristic over the whole class of regular models, with no parametric assumption on customer behaviour.

The theorem is proved on paper. To our knowledge it has no machine-checked proof, and no regular choice model has been formalized on the platform. This mission produces the regular model as a reusable definition, a statement of Theorem 3.3 in which every hypothesis is explicit, and formal versions of the proof's identities and inequalities. The logarithmic step is asserted without proof in the source, so a formal proof of it completes the paper's argument.

Difficulty

Each step is elementary. The difficulty is bookkeeping. The revenue of S∗S^*S∗ has to be regrouped by distinct price levels rather than by products: several products may share a price, and kkk counts values. The regrouping uses an Abel-type rearrangement with the boundary convention Nk+1=0N_{k+1}=0Nk+1​=0. Comparing NiN_iNi​ with the purchase probability of SiS_iSi​ needs regularity twice. First, axiom 4 for products passes from S∗S^*S∗ to S∗∩SiS^*\cap S_iS∗∩Si​. Then the no-purchase case passes from S∗∩SiS^*\cap S_iS∗∩Si​ to SiS_iSi​. A proof that uses axiom 4 only for products fails at the second step, and the claim is false without it.

Formalization scope

  • Products form a finite type C with [Fintype C] [DecidableEq C]. The goal adds [Nonempty C], the paper's k≥1k\ge1k≥1. Offer sets are Finset C, and P\mathcal PP is P : C → Finset C → ℝ. The no-purchase option is not a product: P(0,S)\mathcal P(0,S)P(0,S) is the derived quantity noPurchase P S. IsRegular P carries axioms 1–4, including both no-purchase cases.
  • OPT\mathrm{OPT}OPT is Finset.sup' over all subsets. RO\mathrm{RO}RO is Finset.sup' over the indices 1,…,k1,\dots,k1,…,k of the threshold sets only.
  • Indices are 1-based natural numbers. level r i is rir_iri​ for 1≤i≤k1\le i\le k1≤i≤k. purchaseProfile P r S i is NiN_iNi​ for 1≤i≤k1\le i\le k1≤i≤k and 000 otherwise, which builds in Nk+1=0N_{k+1}=0Nk+1​=0.
  • Optimality is the hypothesis rev(S∗)=OPT\mathrm{rev}(S^*)=\mathrm{OPT}rev(S∗)=OPT. The index ℓ\ellℓ is a variable with hypotheses 1≤ℓ≤k1\le\ell\le k1≤ℓ≤k, Nℓ>0N_\ell>0Nℓ​>0, and Ni≯0N_i\not>0Ni​>0 for ℓ<i≤k\ell<i\le kℓ<i≤k. N1>0N_1>0N1​>0 is kept as in the paper.
  • The approximation factor is stated in product form, OPT≤D⋅RO\mathrm{OPT}\le D\cdot\mathrm{RO}OPT≤D⋅RO, never as a ratio. Both the sum form and the logarithmic form are stated. ln⁡\lnln is Real.log, applied to ν≥1\nu\ge1ν≥1.
  • The milestones other than the goal are stated for an arbitrary S∗S^*S∗, because the proof does not use optimality there.

The following statements are trivial or false and are not this mission: a maximum over all subsets in place of RO\mathrm{RO}RO; regularity without its no-purchase case; a purchase profile taken from a non-optimal set while rev(S∗)\mathrm{rev}(S^*)rev(S∗) is still called the optimum; the logarithmic form alone; a specific choice model (MNL, Markov chain) in place of an arbitrary regular P\mathcal PP.

A complete development needs finite sums regrouped by the values of a function and an elementary logarithm inequality. The regular choice model and the revenue-ordered sets are shared with the other missions of this series and are reusable for any analysis of assortment heuristics. Proofs of any milestone are welcome. So are alternative arguments for the logarithmic step and a proof that NNN is non-increasing.

Selected references

  • G. Berbeglia and G. Joret, Assortment Optimisation Under a General Discrete Choice Model: A Tight Analysis of Revenue-Ordered Assortments, arXiv:1606.01371v3, 2019; Algorithmica 82, 2020. https://arxiv.org/abs/1606.01371v3
  • K. Talluri and G. 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
  • P. Rusmevichientong, D. Shmoys, C. Tong and H. Topaloglu, Assortment Optimization under the Multinomial Logit Model with Random Choice Parameters, Production and Operations Management 23(11), 2014. https://doi.org/10.1111/poms.12191
9 thms2 active usersReviewed
CombinatoricsOperations ResearchOptimization+1·Captain: mikedeng1

Assortment Optimisation Under a General Discrete Choice Model: A Tight Analysis of Revenue-Ordered Assortments III: The Three Revenue-Ordered Approximation Bounds Are TightResearch Paper

Motivation

A retailer or an airline chooses which products to offer, and customers choose among what is offered, or buy nothing. Choosing the offer set that maximizes expected revenue is the assortment problem, central to revenue management (Talluri and van Ryzin, 2004). It is NP-hard even for simple mixtures of logit models, so practice relies on heuristics. The most common one is the revenue-ordered assortments strategy: offer the products whose price is above a threshold, and pick the best threshold.

Berbeglia and Joret (arXiv:1606.01371) analyse this heuristic under every regular choice model, the broad class in which enlarging the offer set never raises the probability of choosing any given alternative, including not buying. They prove three approximation guarantees, then show that none of the three can be improved. This mission formalizes that last result, Theorem 3.4.

Timeline:

  • Talluri and van Ryzin (2004) showed that revenue-ordered assortments are optimal under the multinomial logit model.
  • Rusmevichientong, Shmoys, Tong and Topaloglu (2014) showed the problem NP-hard under a mixture of two logit models, and proved that revenue-ordered assortments earn at least OPT/(e(1+ln⁡(rk/r1)))\mathrm{OPT}/(e(1+\ln(r_k/r_1)))OPT/(e(1+ln(rk​/r1​))) under mixed logit.
  • Aouad, Farias, Levi and Segev (2018) proved an Ω(1/ln⁡(rk/r1))\Omega(1/\ln(r_k/r_1))Ω(1/ln(rk​/r1​)) guarantee under random utility models, and that the assortment problem there is NP-hard to approximate within Ω(1/n1−ϵ)\Omega(1/n^{1-\epsilon})Ω(1/n1−ϵ) and Ω(1/log⁡1−ϵ(rk/r1))\Omega(1/\log^{1-\epsilon}(r_k/r_1))Ω(1/log1−ϵ(rk​/r1​)).
  • Berbeglia and Joret (2016–2020, Algorithmica) proved the guarantees 1/k1/k1/k, 1/∑i(ri−ri−1)/ri≥1/(1+ln⁡(rk/r1))1/\sum_{i}(r_i-r_{i-1})/r_i\ge1/(1+\ln(r_k/r_1))1/∑i​(ri​−ri−1​)/ri​≥1/(1+ln(rk​/r1​)) and a purchase-probability bound under any regular model, and an instance on which all three, in their sum forms, are attained in the limit.

Setting

The products form a finite set C\mathcal CC. A system of choice probabilities gives, for every offer set S⊆CS\subseteq\mathcal CS⊆C and product xxx, the probability P(x,S)\mathcal P(x,S)P(x,S) that a customer buys xxx. The no-purchase probability is P(0,S)=1−∑x∈SP(x,S)\mathcal P(0,S)=1-\sum_{x\in S}\mathcal P(x,S)P(0,S)=1−∑x∈S​P(x,S). The model is regular if

  1. P(x,S)≥0\mathcal P(x,S)\ge0P(x,S)≥0 for x∈C∪{0}x\in\mathcal C\cup\{0\}x∈C∪{0};
  2. P(x,S)=0\mathcal P(x,S)=0P(x,S)=0 for x∉Sx\notin Sx∈/S;
  3. ∑x∈SP(x,S)≤1\sum_{x\in S}\mathcal P(x,S)\le1∑x∈S​P(x,S)≤1;
  4. P(x,S)≥P(x,S′)\mathcal P(x,S)\ge\mathcal P(x,S')P(x,S)≥P(x,S′) for S⊆S′S\subseteq S'S⊆S′ and every x∈S∪{0}x\in S\cup\{0\}x∈S∪{0}.

Each product has a revenue r(x)>0r(x)>0r(x)>0. Offering SSS earns rev(S)=∑x∈SP(x,S) r(x)\mathrm{rev}(S)=\sum_{x\in S}\mathcal P(x,S)\,r(x)rev(S)=∑x∈S​P(x,S)r(x), and OPT=max⁡S⊆Crev(S)\mathrm{OPT}=\max_{S\subseteq\mathcal C}\mathrm{rev}(S)OPT=maxS⊆C​rev(S).

Let r1<⋯<rkr_1<\dots<r_kr1​<⋯<rk​ be the distinct revenues, with r0:=0r_0:=0r0​:=0, and Si={x:r(x)≥ri}S_i=\{x: r(x)\ge r_i\}Si​={x:r(x)≥ri​}. The heuristic earns

RO=max⁡i∈[k]rev(Si).\mathrm{RO}=\max_{i\in[k]}\mathrm{rev}(S_i).RO=i∈[k]max​rev(Si​).

For an optimal S∗S^*S∗ let Ni=∑x∈S∗, r(x)≥riP(x,S∗)N_i=\sum_{x\in S^*,\,r(x)\ge r_i}\mathcal P(x,S^*)Ni​=∑x∈S∗,r(x)≥ri​​P(x,S∗), Nk+1=0N_{k+1}=0Nk+1​=0, and ℓ\ellℓ the largest index with Nℓ>0N_\ell>0Nℓ​>0. Section 3 of the paper proves

  • (A) OPT≤k⋅RO\mathrm{OPT}\le k\cdot\mathrm{RO}OPT≤k⋅RO (Theorem 3.1);
  • (B) OPT≤Dr⋅RO\mathrm{OPT}\le D_r\cdot\mathrm{RO}OPT≤Dr​⋅RO, with Dr=∑i=1kri−ri−1ri≤1+ln⁡(rk/r1)D_r=\sum_{i=1}^{k}\frac{r_i-r_{i-1}}{r_i}\le 1+\ln(r_k/r_1)Dr​=∑i=1k​ri​ri​−ri−1​​≤1+ln(rk​/r1​) (Theorem 3.2);
  • (C) OPT≤DN(S∗)⋅RO\mathrm{OPT}\le D_N(S^*)\cdot\mathrm{RO}OPT≤DN​(S∗)⋅RO, with DN(S∗)=∑i=1ℓNi−Ni+1Ni≤1+ln⁡(N1/Nℓ)D_N(S^*)=\sum_{i=1}^{\ell}\frac{N_i-N_{i+1}}{N_i}\le 1+\ln(N_1/N_\ell)DN​(S∗)=∑i=1ℓ​Ni​Ni​−Ni+1​​≤1+ln(N1​/Nℓ​) (Theorem 3.3).

Formalization targets

Goal: Theorem 3.4

For every k≥1k\ge1k≥1 and every δ>0\delta>0δ>0 there are a finite nonempty product set, a regular P\mathcal PP, revenues r>0r>0r>0 with exactly kkk distinct values, and an optimal S∗S^*S∗ with N1>0N_1>0N1​>0, such that

k⋅RO<(1+δ) OPT,Dr⋅RO<(1+δ) OPT,DN(S∗)⋅RO<(1+δ) OPT.k\cdot\mathrm{RO}<(1+\delta)\,\mathrm{OPT},\qquad D_r\cdot\mathrm{RO}<(1+\delta)\,\mathrm{OPT},\qquad D_N(S^*)\cdot\mathrm{RO}<(1+\delta)\,\mathrm{OPT}.k⋅RO<(1+δ)OPT,Dr​⋅RO<(1+δ)OPT,DN​(S∗)⋅RO<(1+δ)OPT.

So none of the bounds (A), (B), (C) stays true when multiplied by 1+δ1+\delta1+δ, for any number kkk of distinct revenues.

The tight instance (milestones)

The paper's witness has products (i,j)(i,j)(i,j) with i∈[k]i\in[k]i∈[k] and j∈[i]j\in[i]j∈[i]. Product (i,j)(i,j)(i,j) has revenue ε−j\varepsilon^{-j}ε−j, and P((i,j),S)=εi\mathcal P((i,j),S)=\varepsilon^iP((i,j),S)=εi when (i,j)∈S(i,j)\in S(i,j)∈S and (i,1),…,(i,j−1)∉S(i,1),\dots,(i,j-1)\notin S(i,1),…,(i,j−1)∈/S, and 000 otherwise, for 0<ε≤120<\varepsilon\le\tfrac120<ε≤21​. The milestones follow the proof on pp. 10–11:

  • the axiom (iii) bound;
  • (7);
  • (8);
  • (9);
  • regularity;
  • the distinct revenues ri=ε−ir_i=\varepsilon^{-i}ri​=ε−i;
  • RO=rev(C)<1/(1−ε)\mathrm{RO}=\mathrm{rev}(\mathcal C)<1/(1-\varepsilon)RO=rev(C)<1/(1−ε);
  • OPT=k\mathrm{OPT}=kOPT=k, attained by {(i,i)}\{(i,i)\}{(i,i)};
  • the limits OPT/RO→k\mathrm{OPT}/\mathrm{RO}\to kOPT/RO→k and Dr→kD_r\to kDr​→k;
  • Ni=εi+⋯+εkN_i=\varepsilon^i+\dots+\varepsilon^kNi​=εi+⋯+εk and DN→kD_N\to kDN​→k as ε→0+\varepsilon\to0^+ε→0+.

Significance

With Theorems 3.1–3.3, Theorem 3.4 closes the analysis: in terms of the parameters kkk, DrD_rDr​ and DND_NDN​, revenue-ordered assortments are understood exactly under regular choice. It complements the hardness results of Aouad et al., which already show that no efficient strategy can do much better than (A) and (B) in general; Theorem 3.4 shows that the analysis of this particular heuristic is exact. The instance is also a concrete regular choice model in which the optimal offer set is far from every nested one.

On the formal side, the mission produces a reusable formal definition of regular discrete choice models, revenue-ordered assortments and the three bound quantities. These are shared, under other sub-namespaces, with the companion missions on Theorems 3.2 and 3.3. The result is proved in the paper. As far as is known it has not been machine-checked anywhere; the remaining work is formalizing the paper's proof, including the steps it leaves to the reader.

Difficulty

The proof is a construction, and the paper verifies most of it in a sentence each. The work lies in those sentences:

  • the regularity of the instance at the no-purchase option, (9), which needs the row decomposition (8);
  • the claim that {(i,i)}\{(i,i)\}{(i,i)} is optimal among all 2k(k+1)/22^{k(k+1)/2}2k(k+1)/2 offer sets, asserted without proof;
  • the claim that the full set is the best threshold set;
  • the evaluation of DN(S∗)D_N(S^*)DN​(S∗), which needs ℓ=k\ell=kℓ=k.

The naive witness, one instance per bound, does not help: the goal asks for one instance with exactly kkk revenues on which all three bounds are nearly attained, for every kkk.

Formalization scope

  • Products are an arbitrary finite type C. The no-purchase option is not a product: P(0,S)\mathcal P(0,S)P(0,S) is the derived quantity noPurchase P S. Offer sets are Finset C.
  • IsRegular carries axioms (i)–(iv), with (i) and (iv) stated both for products and for the no-purchase option. Regularity at x=0x=0x=0 is essential: without it the guarantees fail.
  • OPT\mathrm{OPT}OPT is Finset.sup' over all subsets, including ∅\emptyset∅. RO\mathrm{RO}RO is Finset.sup' over the kkk threshold sets only, and needs Nonempty C.
  • The distinct revenues are the sorted image of r, indexed by Fin k from 000: the Lean index iii is the paper's i+1i+1i+1, and r0=0r_0=0r0​=0 is a separate case. ℓ\ellℓ lives in WithBot (Fin k).
  • Bounds are multiplicative (OPT≤D⋅RO\mathrm{OPT}\le D\cdot\mathrm{RO}OPT≤D⋅RO); there is no ratio RO/OPT\mathrm{RO}/\mathrm{OPT}RO/OPT in the goal. Limits are along 𝓝[>] 0.
  • The tight instance keeps the paper's 1-based pairs (i,j)(i,j)(i,j) as a subtype of Fin (k+1) × Fin (k+1). Its regularity is a theorem to prove, never a field assumed.

Ruled out:

  • a fixed kkk, since "for every kkk" is the content;
  • tightness of the logarithmic forms 1/(1+ln⁡(rk/r1))1/(1+\ln(r_k/r_1))1/(1+ln(rk​/r1​)) and 1/(1+ln⁡ν)1/(1+\ln\nu)1/(1+lnν), which this instance does not show (there 1+ln⁡(rk/r1)→∞1+\ln(r_k/r_1)\to\infty1+ln(rk​/r1​)→∞ while OPT/RO→k\mathrm{OPT}/\mathrm{RO}\to kOPT/RO→k), so only the sum forms are claimed tight;
  • an instance whose regularity is assumed;
  • a heuristic maximizing over all subsets.

Contributions welcome: proofs of the milestones, especially (8), (9) and the optimality of {(i,i)}\{(i,i)\}{(i,i)}, and general lemmas on sorted distinct values of a finite function.

Selected references

  • G. Berbeglia and G. Joret, Assortment Optimisation Under a General Discrete Choice Model: A Tight Analysis of Revenue-Ordered Assortments, arXiv:1606.01371v3, 2019; published in Algorithmica, 2020. https://arxiv.org/abs/1606.01371
  • K. Talluri and G. 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
  • P. Rusmevichientong, D. Shmoys, C. Tong and H. Topaloglu, Assortment Optimization under the Multinomial Logit Model with Random Choice Parameters, Production and Operations Management 23(11), 2014. https://doi.org/10.1111/poms.12191
  • A. Aouad, V. Farias, R. Levi and D. Segev, The Approximability of Assortment Optimization Under Ranking Preferences, Operations Research 66(6), 2018. https://doi.org/10.1287/opre.2018.1724
15 thms2 active usersReviewed
Algorithmic Game TheoryCombinatoricsOperations Research+2·Captain: mikedeng1

Assortment Optimisation Under a General Discrete Choice Model: A Tight Analysis of Revenue-Ordered Assortments V: Stackelberg Matroid Pricing Is an Assortment Problem Under a Regular Choice ModelResearch Paper

Motivation

In Stackelberg network pricing a leader sets prices on some resources, and a follower then buys the cheapest structure available to them, paying the leader for the priced resources they use. The Stackelberg Minimum Spanning Tree problem, introduced by Cardinal, Demaine, Fiorini, Joret, Langerman, Newman and Weimann (Algorithmica 2011), is the version in which the follower buys a minimum spanning tree. The best known approximation factors for it are those of uniform pricing, which gives every priced edge the same price (Berbeglia–Joret, p. 21).

Independently, revenue management studies assortment optimisation: a seller chooses which products to offer, and customers choose among the offered products according to a discrete choice model. Berbeglia and Joret (arXiv:1606.01371, Algorithmica 2020) analyse the revenue-ordered heuristic (offer every product whose revenue is at least some threshold) under any regular choice model, and prove tight approximation bounds.

§4.6 of that paper shows that the two problems are the same problem: a Stackelberg Matroid pricing instance is an assortment problem under a regular choice model, and uniform pricing is the revenue-ordered heuristic on that model. The bounds of Cardinal et al. on uniform pricing thus become special cases of the general bounds on revenue-ordered assortments. This mission formalizes that correspondence (Theorem 4.16).

Setting

Choice model. Let C\mathcal CC be a finite set of products. A system of choice probabilities assigns to each offered set S⊆CS\subseteq\mathcal CS⊆C and product xxx a probability P(x,S)\mathcal P(x,S)P(x,S), with no-purchase probability P(0,S)=1−∑x∈SP(x,S)\mathcal P(0,S)=1-\sum_{x\in S}\mathcal P(x,S)P(0,S)=1−∑x∈S​P(x,S). It is regular if (i) all probabilities are nonnegative, (ii) P(x,S)=0\mathcal P(x,S)=0P(x,S)=0 for x∉Sx\notin Sx∈/S, (iii) ∑x∈SP(x,S)≤1\sum_{x\in S}\mathcal P(x,S)\le1∑x∈S​P(x,S)≤1, and (iv) P(x,S)≥P(x,S′)\mathcal P(x,S)\ge\mathcal P(x,S')P(x,S)≥P(x,S′) whenever S⊆S′S\subseteq S'S⊆S′ and x∈S∪{0}x\in S\cup\{0\}x∈S∪{0}. With revenues r:C→R>0r:\mathcal C\to\mathbb R_{>0}r:C→R>0​, rev(S)=∑x∈SP(x,S)r(x)\mathrm{rev}(S)=\sum_{x\in S}\mathcal P(x,S)r(x)rev(S)=∑x∈S​P(x,S)r(x) and OPT=max⁡Srev(S)\mathrm{OPT}=\max_S\mathrm{rev}(S)OPT=maxS​rev(S). The revenue-ordered assortment generated by yyy is {y′∈C:r(y′)≥r(y)}\{y'\in\mathcal C:r(y')\ge r(y)\}{y′∈C:r(y′)≥r(y)}.

Greedy algorithm. Given a family of independent sets, a linear ordering LLL and a set FFF, greedy(F,L)\mathrm{greedy}(F,L)greedy(F,L) scans FFF in the order induced by LLL, starting from ∅\emptyset∅, and adds each element that keeps the current set independent.

Stackelberg Matroid problem. An instance is a matroid M=(E,X)M=(E,\mathcal X)M=(E,X), a bipartition E=R⊔BE=R\sqcup BE=R⊔B into red and blue elements, and red costs c:R→R>0c:R\to\mathbb R_{>0}c:R→R>0​; some base of MMM lies inside RRR. The leader chooses prices p:B→R>0p:B\to\mathbb R_{>0}p:B→R>0​. The customer buys a minimum-weight base by running greedy on R∪BR\cup BR∪B with an ordering L∗L^*L∗ that is non-decreasing in weight (ccc on RRR, ppp on BBB) and puts blue elements first on ties. The leader earns revStack(p,L∗)=∑e∈B∩greedyM(R∪B,L∗)p(e)\mathrm{rev}_{\mathrm{Stack}}(p,L^*)=\sum_{e\in B\cap\mathrm{greedy}_M(R\cup B,L^*)}p(e)revStack​(p,L∗)=∑e∈B∩greedyM​(R∪B,L∗)​p(e).

The assortment instance. Let c1<⋯<ckc_1<\dots<c_kc1​<⋯<ck​ be the distinct red costs. Products are C=B×{c1,…,ck}\mathcal C=B\times\{c_1,\dots,c_k\}C=B×{c1​,…,ck​} with r((e,q))=∣B∣ qr((e,q))=|B|\,qr((e,q))=∣B∣q. The auxiliary matroid M′M'M′ on R∪CR\cup\mathcal CR∪C declares XXX independent when it holds at most one pair (e,q)(e,q)(e,q) per blue eee and (R∩X)∪{e:(e,q)∈X}(R\cap X)\cup\{e:(e,q)\in X\}(R∩X)∪{e:(e,q)∈X} is independent in MMM. With an ordering LLL of R∪CR\cup\mathcal CR∪C that is non-decreasing in cost and puts products before red elements on ties, P((e,q),S)=1/∣B∣\mathcal P((e,q),S)=1/|B|P((e,q),S)=1/∣B∣ if (e,q)∈greedyM′(R∪S,L)(e,q)\in\mathrm{greedy}_{M'}(R\cup S,L)(e,q)∈greedyM′​(R∪S,L) and 000 otherwise.

Formalization targets

Goal: Theorem 4.16 (p. 24)

For the instance above, with B≠∅B\ne\emptysetB=∅ and every admissible LLL:

P is regular,r>0,max⁡S⊆Crev(S)=max⁡p>0, L∗revStack(p,L∗),\mathcal P\ \text{is regular},\quad r>0,\qquad \max_{S\subseteq\mathcal C}\mathrm{rev}(S)=\max_{p>0,\ L^*}\mathrm{rev}_{\mathrm{Stack}}(p,L^*),P is regular,r>0,S⊆Cmax​rev(S)=p>0, L∗max​revStack​(p,L∗),

and uniform pricing corresponds to revenue-ordered assortments: for each cic_ici​ some revenue-ordered assortment earns the revenue of the uniform price cic_ici​, and every revenue-ordered assortment earns the revenue of some uniform price cic_ici​.

Milestones

  1. Lemma 4.14 (p. 23): for F⊆F′F\subseteq F'F⊆F′, ∣greedyM(F′,L)∣≥∣greedyM(F,L)∣|\mathrm{greedy}_M(F',L)|\ge|\mathrm{greedy}_M(F,L)|∣greedyM​(F′,L)∣≥∣greedyM​(F,L)∣ and F∩greedyM(F′,L)⊆greedyM(F,L)F\cap\mathrm{greedy}_M(F',L)\subseteq\mathrm{greedy}_M(F,L)F∩greedyM​(F′,L)⊆greedyM​(F,L).
  2. Property (14) (p. 31): orderings agreeing on a block partition make greedy pick the same number of elements per block.
  3. Lemma 4.15 (p. 23): the Stackelberg revenue does not depend on the compatible ordering.
  4. M′M'M′ is a matroid (p. 32).
  5. Regularity of P\mathcal PP (pp. 32–33).
  6. rev(S)\mathrm{rev}(S)rev(S) equals the Stackelberg revenue of pS(e)=min⁡{q:(e,q)∈S}p_S(e)=\min\{q:(e,q)\in S\}pS​(e)=min{q:(e,q)∈S} (pp. 33–34).
  7. Rounding prices up to the cost levels does not decrease revenue (p. 34).
  8. For prices in {c1,…,ck}∪{+∞}\{c_1,\dots,c_k\}\cup\{+\infty\}{c1​,…,ck​}∪{+∞}, rev(Sp)\mathrm{rev}(S_p)rev(Sp​) equals the Stackelberg revenue of ppp (pp. 34–35).

Significance

Theorem 4.16 transfers every guarantee for revenue-ordered assortments under regular models to uniform pricing in Stackelberg Matroid pricing. Specialising Theorems 3.1, 3.2 and 3.3 of the paper through it gives exactly the three bounds on uniform pricing proved by Cardinal et al. for Stackelberg Minimum Spanning Tree (Theorem 3 of their paper), which those authors showed to be tight (p. 24). It also places uniform pricing inside a general picture: it is a threshold policy for a regular choice model whose purchase probabilities come from a matroid greedy algorithm, and the paper notes that the argument extends to polymatroids.

The result is proved in the paper, with two steps left to the reader (M′M'M′ is a matroid) or justified briefly (rounding prices). No machine-checked proof exists. A formalization also produces a reusable greedy algorithm on finite matroids with its monotonicity properties (Lemma 4.14, (14)), which Mathlib at the pinned revision does not have.

Difficulty

The obvious argument identifies the customer's run of greedy on (M,R∪B)(M,R\cup B)(M,R∪B) with the run of greedy on (M′,R∪S)(M',R\cup S)(M′,R∪S) step by step. That identification fails as stated: the choice probabilities are defined with one fixed ordering LLL of R∪CR\cup\mathcal CR∪C, while the customer's ordering L∗L^*L∗ is any ordering compatible with the prices, and ties among blue elements of equal price, or between a product and a red element of equal cost, may be broken differently. Equality of revenues therefore needs the tie-independence property (14) applied to M′M'M′ as well as to MMM, which in turn needs M′M'M′ to be a matroid. A second gap is axiom (iv) for the no-purchase option: it requires that offering more products never decreases the number of products greedy selects, which is Lemma 4.14 (i)–(ii) for M′M'M′ combined, not a pointwise statement. Finally, the paper's "+∞+\infty+∞" price must be handled with the red base: without a base inside RRR, a blue element priced above all costs can still be bought and the optimum is unbounded.

Formalization scope

  • Elements and orderings. The matroid is a Mathlib Matroid α; RRR, BBB are Finset α with R∪BR\cup BR∪B equal to the ground set; costs and prices are functions α → ℝ, used only on RRR and BBB. A linear ordering is a duplicate-free List covering the set; greedy on FFF scans the whole list and skips elements outside FFF, so FFF is never re-sorted.
  • Greedy is defined for an arbitrary independence predicate on finite sets, so that it applies to M′M'M′ before M′M'M′ is shown to be a matroid.
  • Customer model. Compatibility (15a)–(15b), including blue priority on ties, is part of the definition of an admissible customer ordering. The red base is a field of every instance.
  • Assortment instance. Elements of R∪CR\cup\mathcal CR∪C live in α ⊕ (α × ℝ). The paper's conditions (1)–(2) for M′M'M′ contain two typos (a missing "∈X\in X∈X" and X\mathcal XX for XXX); the intended conditions are formalized. The price +∞+\infty+∞ is the real number 1+∑f∈Rc(f)1+\sum_{f\in R}c(f)1+∑f∈R​c(f), strictly above every red cost.
  • Goal shape. The theorem is stated for the explicit instance of the proof, for every admissible LLL. The optimum of the Stackelberg side is a maximum (IsGreatest) over positive real prices and all compatible customer orderings. Cost levels are quantified as elements of the set of red costs rather than by index.
  • Ruled out. An existential statement ("some regular instance has the same optimum") is met by a single product with revenue OPTStack\mathrm{OPT}_{\mathrm{Stack}}OPTStack​ and is not this theorem; so is a customer model without tie-breaking, a formalization without the red base, or Lemma 4.14 only for F=F′F=F'F=F′.

Contributions welcome: proofs of Lemma 4.14 and (14) for Mathlib matroids (reusable beyond this mission), the matroid property of parallel extensions such as M′M'M′, and the three revenue identities.

Selected references

  • G. Berbeglia and G. Joret, Assortment Optimisation Under a General Discrete Choice Model: A Tight Analysis of Revenue-Ordered Assortments, arXiv:1606.01371v3, 2019 (Algorithmica, 2020). https://arxiv.org/abs/1606.01371
  • J. Cardinal, E. D. Demaine, S. Fiorini, G. Joret, S. Langerman, I. Newman and O. Weimann, The Stackelberg Minimum Spanning Tree Game, Algorithmica 59(2):129–144, 2011. https://doi.org/10.1007/s00453-009-9299-y
  • A. Schrijver, Combinatorial Optimization: Polyhedra and Efficiency, Vol. B, Algorithms and Combinatorics 24, Springer, 2003 (matroids and the greedy algorithm). https://link.springer.com/book/9783540443896
14 thms1 active userReviewed
Operations ResearchOptimization·Captain: mikedeng1

Assortment Optimisation Under a General Discrete Choice Model: A Tight Analysis of Revenue-Ordered Assortments VI: Revenue-Ordered Offer Sets Grow with Capacity Left and Shrink with Time LeftResearch Paper

Motivation

In airline and hotel revenue management a firm sells a fixed stock of a perishable resource (seats on a flight leg, rooms on a night) over a finite selling horizon, and in each period decides which fare classes to open. Customers do not buy a fixed fare: they choose among the fares on offer, or leave. Talluri and van Ryzin (Management Science, 2004) formulated this single-leg, choice-based problem as a dynamic program over the time remaining and the units remaining.

Practical revenue management systems rarely offer arbitrary sets of fares. They use nested controls: fares are opened from the most expensive downwards, so the open set is always "every fare above some threshold". Whether restricting to such revenue-ordered offer sets is natural depends on how the optimal threshold moves as the state changes. For mixtures of multinomial logit models, Rusmevichientong, Shmoys, Tong and Topaloglu (Production and Operations Management, 2014, Theorem 6) showed that the best revenue-ordered threshold moves monotonically in both the remaining capacity and the remaining time.

Berbeglia and Joret (arXiv:1606.01371v3, §5, Theorem 5.1) extend these two monotonicity properties to every regular discrete choice model, the class that contains essentially every model used in revenue management, including all random utility models. This mission formalizes that result.

Setting

A finite nonempty set C\mathcal CC of products is sold. For a choice set S⊆CS\subseteq\mathcal CS⊆C and a product xxx, P(x,S)\mathcal P(x,S)P(x,S) is the probability that a customer offered SSS buys xxx; the no-purchase probability is P(0,S)=1−∑x∈SP(x,S)\mathcal P(0,S)=1-\sum_{x\in S}\mathcal P(x,S)P(0,S)=1−∑x∈S​P(x,S). The system P\mathcal PP is regular when (i) all these probabilities are nonnegative, (ii) P(x,S)=0\mathcal P(x,S)=0P(x,S)=0 for x∉Sx\notin Sx∈/S, (iii) ∑x∈SP(x,S)≤1\sum_{x\in S}\mathcal P(x,S)\le 1∑x∈S​P(x,S)≤1, and (iv) P(x,S)≥P(x,S′)\mathcal P(x,S)\ge\mathcal P(x,S')P(x,S)≥P(x,S′) whenever S⊆S′S\subseteq S'S⊆S′, for every x∈Sx\in Sx∈S and for the no-purchase option x=0x=0x=0.

Each product xxx has a revenue r(x)>0r(x)>0r(x)>0. Let r1<r2<⋯<rkr_1<r_2<\cdots<r_kr1​<r2​<⋯<rk​ be the distinct revenues. For ℓ∈[k]={1,…,k}\ell\in[k]=\{1,\dots,k\}ℓ∈[k]={1,…,k} the revenue-ordered assortment is Sℓ={x∈C:r(x)≥rℓ}S_\ell=\{x\in\mathcal C:r(x)\ge r_\ell\}Sℓ​={x∈C:r(x)≥rℓ​}; a larger index ℓ\ellℓ gives a smaller set, S1=CS_1=\mathcal CS1​=C.

In the multi-period model one customer arrives per period. With ttt periods remaining and qqq units left, the firm offers some SℓS_\ellSℓ​. The values of the revenue-ordered dynamic program are

Jt(q,ℓ)=∑x∈SℓP(x,Sℓ)(r(x)+Jt−1(q−1))+P(0,Sℓ) Jt−1(q)(t,q>0),\mathcal J_t(q,\ell)=\sum_{x\in S_\ell}\mathcal P(x,S_\ell)\bigl(r(x)+\mathcal J_{t-1}(q-1)\bigr)+\mathcal P(0,S_\ell)\,\mathcal J_{t-1}(q)\qquad(t,q>0),Jt​(q,ℓ)=x∈Sℓ​∑​P(x,Sℓ​)(r(x)+Jt−1​(q−1))+P(0,Sℓ​)Jt−1​(q)(t,q>0),

Jt(q,ℓ)=0\mathcal J_t(q,\ell)=0Jt​(q,ℓ)=0 when t=0t=0t=0 or q=0q=0q=0, and Jt(q)=max⁡ℓ∈[k]Jt(q,ℓ)\mathcal J_t(q)=\max_{\ell\in[k]}\mathcal J_t(q,\ell)Jt​(q)=maxℓ∈[k]​Jt​(q,ℓ). The optimal revenue-ordered index is the smallest maximiser,

ℓt∗(q)=min⁡{ℓ∈[k]:Jt(q,ℓ)=Jt(q)},\ell^*_t(q)=\min\{\ell\in[k]:\mathcal J_t(q,\ell)=\mathcal J_t(q)\},ℓt∗​(q)=min{ℓ∈[k]:Jt​(q,ℓ)=Jt​(q)},

and ΔJt(q)=Jt(q)−Jt(q−1)\Delta\mathcal J_t(q)=\mathcal J_t(q)-\mathcal J_t(q-1)ΔJt​(q)=Jt​(q)−Jt​(q−1) is the marginal value of capacity.

Formalization targets

Goal: Theorem 5.1

For every t≥1t\ge 1t≥1 and q≥1q\ge 1q≥1,

ℓt∗(q)≤ℓt∗(q−1)  if q≥2,ℓt∗(q)≥ℓt−1∗(q)  if t≥2.\ell^*_t(q)\le\ell^*_t(q-1)\ \text{ if } q\ge 2,\qquad \ell^*_t(q)\ge\ell^*_{t-1}(q)\ \text{ if } t\ge 2 .ℓt∗​(q)≤ℓt∗​(q−1)  if q≥2,ℓt∗​(q)≥ℓt−1∗​(q)  if t≥2.

More units left give a weakly larger optimal offer set; more periods left give a weakly smaller one. The paper states the theorem for t∈[T]t\in[T]t∈[T], q∈[Q]q\in[Q]q∈[Q]; the horizon and the capacity only bound ttt and qqq, so the goal is stated for all t,q≥1t,q\ge 1t,q≥1.

Milestones

  1. Lemma 2.1 (p. 6): ∑x∈SP(x,S)≤∑x∈S′P(x,S′)\sum_{x\in S}\mathcal P(x,S)\le\sum_{x\in S'}\mathcal P(x,S')∑x∈S​P(x,S)≤∑x∈S′​P(x,S′) for S⊆S′S\subseteq S'S⊆S′.
  2. Lemma .1 (p. 35): with L∗(δ)\mathcal L^*(\delta)L∗(δ) the set of indices ℓ\ellℓ maximising ∑x∈SℓP(x,Sℓ)(r(x)+δ)\sum_{x\in S_\ell}\mathcal P(x,S_\ell)(r(x)+\delta)∑x∈Sℓ​​P(x,Sℓ​)(r(x)+δ), if δ1+rk≥0\delta_1+r_k\ge 0δ1​+rk​≥0 and δ1≤δ2\delta_1\le\delta_2δ1​≤δ2​ then min⁡L∗(δ2)≤min⁡L∗(δ1)\min\mathcal L^*(\delta_2)\le\min\mathcal L^*(\delta_1)minL∗(δ2​)≤minL∗(δ1​).
  3. Equation (16) (p. 36): Jt(q)=max⁡ℓ∑x∈SℓP(x,Sℓ)(r(x)−ΔJt−1(q))+Jt−1(q)\mathcal J_t(q)=\max_{\ell}\sum_{x\in S_\ell}\mathcal P(x,S_\ell)(r(x)-\Delta\mathcal J_{t-1}(q))+\mathcal J_{t-1}(q)Jt​(q)=maxℓ​∑x∈Sℓ​​P(x,Sℓ​)(r(x)−ΔJt−1​(q))+Jt−1​(q).
  4. Equations (17)–(18) (p. 36): ℓt∗(q)=min⁡L∗(−ΔJt−1(q))\ell^*_t(q)=\min\mathcal L^*(-\Delta\mathcal J_{t-1}(q))ℓt∗​(q)=minL∗(−ΔJt−1​(q)).
  5. Marginal value at most rkr_krk​ (p. 36): ΔJt(q)≤rk\Delta\mathcal J_t(q)\le r_kΔJt​(q)≤rk​.
  6. Marginal value non-increasing in capacity (p. 36, citing Talluri–van Ryzin, Lemma 4): ΔJt(q+1)≤ΔJt(q)\Delta\mathcal J_t(q+1)\le\Delta\mathcal J_t(q)ΔJt​(q+1)≤ΔJt​(q).
  7. Marginal value non-decreasing in time (p. 37, citing Talluri–van Ryzin, Lemma 5): ΔJt(q)≤ΔJt+1(q)\Delta\mathcal J_t(q)\le\Delta\mathcal J_{t+1}(q)ΔJt​(q)≤ΔJt+1​(q).

Milestones 5–7 are asserted or cited in the paper, not proved there; milestones 6–7 are stated for this restricted dynamic program, which is the one the paper applies them to.

Significance

The theorem says that a firm restricted to revenue-ordered offer sets can implement its policy as a nested booking control: as seats sell out, the threshold fare can only rise, and as departure approaches with seats in hand, it can only fall. This is the structure that standard revenue management systems already assume, and Rusmevichientong et al. point out that such monotonicity can be used to implement them. The paper's contribution is that the property depends only on regularity of the choice model, not on its multinomial-logit form.

Formalizing it adds a machine-checked statement of the single-leg choice-based dynamic program over a general choice model, a precise account of the tie-breaking rule (the smallest maximiser), and machine-checked versions of the two marginal-value monotonicity facts that the paper cites from Talluri and van Ryzin rather than proves. To our knowledge none of these statements has been formalized in any proof assistant; the paper's proofs are informal.

Difficulty

The paper's proof is short, but two of its steps are not proved there. The marginal-value inequalities are imported from Talluri and van Ryzin, whose dynamic program optimises over all offer sets; here only the kkk revenue-ordered sets are available and the empty set is not, so those proofs have to be redone for this program, and the bound ΔJ≤rk\Delta\mathcal J\le r_kΔJ≤rk​ is needed to keep each stage's shifted problem well behaved. The claim ΔJ≤rk\Delta\mathcal J\le r_kΔJ≤rk​ itself is justified on the page only by an informal sentence.

The second subtlety is tie-breaking. The monotonicity is about the smallest optimal index. Several indices can be optimal at once, and the comparison of Lemma .1 transfers optimality of one index from one shift to another only through Lemma 2.1, that is, through the no-purchase case of the regularity axiom. Without that case the statement fails.

Formalization scope

  • Products form a finite nonempty type C; choice probabilities are P : C → Finset C → ℝ, defined on all pairs. The no-purchase option is not a product: its probability is the derived quantity noPurchase P S = 1 − ∑_{x∈S} P x S.
  • IsChoiceSystem P is axioms (i)–(iii); IsRegular P adds axiom (iv) for products and for the no-purchase option. Milestones 5–7 assume only axioms (i)–(iii), as the paper says suffices; Lemma .1, (16), (18) and the goal assume regularity and r>0r>0r>0.
  • Revenue levels use the paper's 1-based index: level r ℓ is rℓr_\ellrℓ​ for 1≤ℓ≤k1\le\ell\le k1≤ℓ≤k, topRevenue r is rkr_krk​, and roSet r ℓ is SℓS_\ellSℓ​ as a Finset. The paper's sorting of products (Sℓ={1,…,j(ℓ)}S_\ell=\{1,\dots,j(\ell)\}Sℓ​={1,…,j(ℓ)}) is notation only and is not used.
  • The dynamic program is defined by structural recursion on ttt; there is no stochastic process. Maxima over [k][k][k] are Finset.sup', and the minima defining ℓt∗(q)\ell^*_t(q)ℓt∗​(q) and min⁡L∗(δ)\min\mathcal L^*(\delta)minL∗(δ) are Finset.min' over sets proved nonempty.
  • ΔJt(q)\Delta\mathcal J_t(q)ΔJt​(q) is used only for q≥1q\ge 1q≥1; milestones are stated with t+1t+1t+1, q+1q+1q+1 in place of the paper's t−1t-1t−1, q−1q-1q−1 to avoid natural-number subtraction, which the goal keeps under its hypotheses q≥2q\ge 2q≥2, t≥2t\ge 2t≥2.

Formalizations that would trivialize the statement are ruled out: the dynamic program maximises over the kkk revenue-ordered sets only, never over all subsets (that is Talluri and van Ryzin's program); ℓ∗\ell^*ℓ∗ is the smallest maximiser, never the largest; and TTT, QQQ, kkk and the choice model are arbitrary, never fixed.

The definitions (regular choice model, revenue-ordered assortments, the restricted dynamic program) are reusable for other results on nested policies. Proofs of the cited Talluri–van Ryzin lemmas for this program are especially welcome, as they are the part the paper leaves to the literature.

Selected references

  • G. Berbeglia and G. Joret, Assortment Optimisation Under a General Discrete Choice Model: A Tight Analysis of Revenue-Ordered Assortments, arXiv:1606.01371v3, 2019; Algorithmica 82, 2020. https://arxiv.org/abs/1606.01371
  • K. Talluri and G. 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
  • P. Rusmevichientong, D. Shmoys, C. Tong and H. Topaloglu, Assortment Optimization under the Multinomial Logit Model with Random Choice Parameters, Production and Operations Management 23(11), 2014. https://doi.org/10.1111/poms.12191
14 thms2 active usersReviewed
CombinatoricsOperations ResearchOptimization·Captain: mikedeng1

Assortment Optimization under Variants of the Nested Logit Model 6: For δ > 1, the Powers-of-δ LP Optimum Scaled by (δ^(2γ̄+1), δ^(γ̄+1)) Is Feasible for the Full LPResearch Paper

Motivation

Assortment optimization asks which products a retailer should offer when customers choose among the offered products according to a probabilistic choice model, so as to maximize expected revenue. Under the nested logit model the products are grouped into nests: a customer first picks a nest, then a product inside it. The model is the standard relaxation of the independence of irrelevant alternatives property of the multinomial logit, and it is used throughout revenue management and transportation demand modelling.

Davis, Gallego and Topaloglu (Oper. Res. 62(2), 2014) classify the complexity of the assortment problem under four variants of the nested logit model, according to whether the dissimilarity parameters are at most one and whether a customer who picks a nest always buys there. The problem is NP-hard as soon as some dissimilarity parameter exceeds one (their Theorem 5), and also when the nests have positive no-purchase weights (Theorem 8), so for the general variant approximation is the realistic aim. §6.2 of the paper gives an approximation scheme for the most general variant: for any δ>1\delta > 1δ>1 it restricts each nest to a short list of candidate assortments indexed by the powers of δ\deltaδ, solves a small linear program, and loses at most a factor δ2γˉ+1\delta^{2\bar\gamma+1}δ2γˉ​+1 of the optimal revenue. This mission formalizes the guarantee behind that scheme, Theorem 12.

Setting

There are nests MMM and products N={1,…,n}N = \{1, \dots, n\}N={1,…,n} in each nest. Product jjj of nest iii has revenue rij≥0r_{ij} \ge 0rij​≥0 and preference weight vij>0v_{ij} > 0vij​>0; the products are ordered so that ri1≥⋯≥rinr_{i1} \ge \dots \ge r_{in}ri1​≥⋯≥rin​. Nest iii has a no-purchase weight vi0≥0v_{i0} \ge 0vi0​≥0 and a dissimilarity parameter γi>0\gamma_i > 0γi​>0, and v0≥0v_0 \ge 0v0​≥0 is the weight of leaving without choosing a nest. Offering Si⊆NS_i \subseteq NSi​⊆N in nest iii gives

Vi(Si)=vi0+∑j∈Sivij,Ri(Si)=∑j∈SirijvijVi(Si),Π(S1,…,Sm)=∑iVi(Si)γiRi(Si)v0+∑iVi(Si)γi.V_i(S_i) = v_{i0} + \sum_{j \in S_i} v_{ij}, \qquad R_i(S_i) = \frac{\sum_{j \in S_i} r_{ij} v_{ij}}{V_i(S_i)}, \qquad \Pi(S_1, \dots, S_m) = \frac{\sum_i V_i(S_i)^{\gamma_i} R_i(S_i)}{v_0 + \sum_i V_i(S_i)^{\gamma_i}}.Vi​(Si​)=vi0​+j∈Si​∑​vij​,Ri​(Si​)=Vi​(Si​)∑j∈Si​​rij​vij​​,Π(S1​,…,Sm​)=v0​+∑i​Vi​(Si​)γi​∑i​Vi​(Si​)γi​Ri​(Si​)​.

The optimal expected revenue Z∗=max⁡ΠZ^* = \max \PiZ∗=maxΠ is the optimal value of the linear program (3): minimize xxx subject to v0x≥∑iyiv_0 x \ge \sum_i y_iv0​x≥∑i​yi​ and yi≥Vi(Si)γi(Ri(Si)−x)y_i \ge V_i(S_i)^{\gamma_i}(R_i(S_i) - x)yi​≥Vi​(Si​)γi​(Ri​(Si​)−x) for every nest iii and every Si⊆NS_i \subseteq NSi​⊆N. Problem (4) keeps the second family of constraints only for a candidate collection of assortments in each nest.

Let γˉ=max⁡iγi\bar\gamma = \max_i \gamma_iγˉ​=maxi​γi​, assumed >1> 1>1 throughout §6, and fix δ>1\delta > 1δ>1. Put viL=vi0+min⁡jvijv^L_i = v_{i0} + \min_j v_{ij}viL​=vi0​+minj​vij​, viU=vi0+∑jvijv^U_i = v_{i0} + \sum_j v_{ij}viU​=vi0​+∑j​vij​, and let liLl^L_iliL​, liUl^U_iliU​ be the least integers with δl≥viL\delta^{l} \ge v^L_iδl≥viL​, δl≥viU\delta^l \ge v^U_iδl≥viU​. For each level l=liL,…,liUl = l^L_i, \dots, l^U_il=liL​,…,liU​, problem (15) maximizes ∑j∈Srijvij\sum_{j \in S} r_{ij} v_{ij}∑j∈S​rij​vij​ over the assortments with δl−1≤Vi(S)≤δl\delta^{l-1} \le V_i(S) \le \delta^lδl−1≤Vi​(S)≤δl; its value is G^il\hat G_{il}G^il​. An assortment S^il\hat S_{il}S^il​ is feasible for (15) and satisfies δ∑j∈S^ilrijvij≥G^il\delta \sum_{j \in \hat S_{il}} r_{ij} v_{ij} \ge \hat G_{il}δ∑j∈S^il​​rij​vij​≥G^il​. The candidate collection of nest iii is {S^il:l=liL,…,liU}∪{∅}\{\hat S_{il} : l = l^L_i, \dots, l^U_i\} \cup \{\emptyset\}{S^il​:l=liL​,…,liU​}∪{∅}.

Formalization targets

Goal: Theorem 12 (p. 28)

If (x^,y^)(\hat x, \hat y)(x^,y^​) is an optimal solution of problem (4) over the candidate collections {S^il}∪{∅}\{\hat S_{il}\} \cup \{\emptyset\}{S^il​}∪{∅}, then

(δ2γˉ+1x^, δγˉ+1y^) is feasible for problem (3).\big(\delta^{2\bar\gamma+1}\hat x,\ \delta^{\bar\gamma+1}\hat y\big) \text{ is feasible for problem (3).}(δ2γˉ​+1x^, δγˉ​+1y^​) is feasible for problem (3).

Milestones (Appendix A.6, pp. 53–54)

  1. x^≥0\hat x \ge 0x^≥0.
  2. Every nonempty assortment of nest iii lies in some level liL≤l≤liUl^L_i \le l \le l^U_iliL​≤l≤liU​.
  3. If δl−1≤a≤δl\delta^{l-1} \le a \le \delta^lδl−1≤a≤δl, then aγ−1≥(δl)γ−1δ−[γ−1]+a^{\gamma-1} \ge (\delta^l)^{\gamma-1}\delta^{-[\gamma-1]^+}aγ−1≥(δl)γ−1δ−[γ−1]+.
  4. Under the same hypothesis, (δl)γ−1≥δ−[1−γ]+aγ−1(\delta^l)^{\gamma-1} \ge \delta^{-[1-\gamma]^+} a^{\gamma-1}(δl)γ−1≥δ−[1−γ]+aγ−1.
  5. δγˉδ−[γi−1]+δ−[1−γi]+≥1\delta^{\bar\gamma}\delta^{-[\gamma_i-1]^+}\delta^{-[1-\gamma_i]^+} \ge 1δγˉ​δ−[γi​−1]+δ−[1−γi​]+≥1 and δγˉ+γi+1≤δ2γˉ+1\delta^{\bar\gamma+\gamma_i+1} \le \delta^{2\bar\gamma+1}δγˉ​+γi​+1≤δ2γˉ​+1.
  6. The scaled pair satisfies δγˉ+1y^i≥Vi(Si)γi(Ri(Si)−δ2γˉ+1x^)\delta^{\bar\gamma+1}\hat y_i \ge V_i(S_i)^{\gamma_i}(R_i(S_i) - \delta^{2\bar\gamma+1}\hat x)δγˉ​+1y^​i​≥Vi​(Si​)γi​(Ri​(Si​)−δ2γˉ​+1x^) for every nonempty SiS_iSi​.
  7. The same inequality for Si=∅S_i = \emptysetSi​=∅.

Companions

  • The guarantee stated after Theorem 12: with v0>0v_0 > 0v0​>0, the assortment assembled from the candidates solving problem (5) earns at least Z∗/δ2γˉ+1Z^*/\delta^{2\bar\gamma+1}Z∗/δ2γˉ​+1.
  • Proposition 15 (p. 51): when (15) is feasible, one of the explicit assortments S^(JL,JS)\hat S(J_L, J_S)S^(JL​,JS​), built from at most q=⌈δ/(δ−1)⌉q = \lceil \delta/(\delta-1)\rceilq=⌈δ/(δ−1)⌉ large and at most qqq small products by a greedy continuous knapsack, is feasible for (15) and within a factor δ\deltaδ of G^il\hat G_{il}G^il​.
  • The count liU−liL≤1+log⁡δ(viU/viL)l^U_i - l^L_i \le 1 + \log_\delta(v^U_i/v^L_i)liU​−liL​≤1+logδ​(viU​/viL​) (p. 28).

Significance

Theorem 12 is the analytical half of the approximation scheme. Together with the paper's Theorem 1 it shows that a linear program with 1+m1 + m1+m variables and at most 1+m(2+log⁡δ(viU/viL))1 + m(2 + \log_\delta(v^U_i/v^L_i))1+m(2+logδ​(viU​/viL​)) constraints yields an assortment within a factor δ2γˉ+1\delta^{2\bar\gamma+1}δ2γˉ​+1 of the optimum, for the most general variant, which is NP-hard, and Proposition 15 makes the candidate assortments computable. Letting δ↓1\delta \downarrow 1δ↓1 trades accuracy for running time, so the result is the paper's answer to how well the general problem can be approximated by this LP approach.

The result is proved in the paper; the work here is to formalize the known proof. To the best of a search of the Prove2Me library (local index and platform mirror, October 2026), no nested logit approximation result, and no knapsack lemma matching Proposition 15, has a machine-checked statement or proof. The definitions of the shared model (the instance, ViV_iVi​, RiR_iRi​, Π\PiΠ, the linear programs (3) and (4)) are common to the six missions of this series.

Difficulty

The obvious argument compares an arbitrary assortment SiS_iSi​ with the candidate of its level and uses the candidate's constraint in (4). This fails as a direct comparison because the nest weight Vi(⋅)γiV_i(\cdot)^{\gamma_i}Vi​(⋅)γi​ enters twice with different exponents, γi−1\gamma_i - 1γi​−1 in front of the revenue and γi\gamma_iγi​ in front of x^\hat xx^, and within one level ViV_iVi​ may vary by a factor δ\deltaδ. When γi>1\gamma_i > 1γi​>1 and when γi≤1\gamma_i \le 1γi​≤1 the monotonicity of t↦tγi−1t \mapsto t^{\gamma_i - 1}t↦tγi​−1 goes in opposite directions, so the losses must be tracked separately by [γi−1]+[\gamma_i - 1]^+[γi​−1]+ and [1−γi]+[1 - \gamma_i]^+[1−γi​]+, and they are absorbed by the common factor δγˉ\delta^{\bar\gamma}δγˉ​ only because γˉ>1\bar\gamma > 1γˉ​>1. The other half, Proposition 15, concerns a knapsack with both a lower and an upper bound on the total weight, where a feasible solution must be produced as well as a good objective value.

Formalization scope

Products are Fin n (0,…,n−10, \dots, n-10,…,n−1 for 1,…,n1, \dots, n1,…,n), nests a finite type, powers VγV^{\gamma}Vγ real powers, and x/0=0x/0 = 0x/0=0, which gives Ri(∅)=0R_i(\emptyset) = 0Ri​(∅)=0. The shared model carries the standing assumptions of §1 with the disclosed pins vij>0v_{ij} > 0vij​>0, rij≥0r_{ij} \ge 0rij​≥0 and γi>0\gamma_i > 0γi​>0 (the page allows γi=0\gamma_i = 0γi​=0 and zero-weight padding products, under which its convention 0γi=00^{\gamma_i} = 00γi​=0 and its proofs fail). Every statement of the mission also assumes n≥1n \ge 1n≥1 (for viLv^L_iviL​), γˉ\bar\gammaγˉ​ is the greatest of the γi\gamma_iγi​ (so at least one nest exists) with γˉ>1\bar\gamma > 1γˉ​>1, and δ>1\delta > 1δ>1. Integer powers δl\delta^lδl are zpow, and liLl^L_iliL​, liUl^U_iliU​ are written ⌈log⁡δviL⌉\lceil \log_\delta v^L_i\rceil⌈logδ​viL​⌉, ⌈log⁡δviU⌉\lceil \log_\delta v^U_i\rceil⌈logδ​viU​⌉, which equal the minima of the paper. An optimal solution of (4) is a feasible pair whose xxx is minimal among feasible pairs. The companion guarantee adds v0>0v_0 > 0v0​>0, the pin of Theorem 1. In Proposition 15 the capacity row of (33) includes v0v_0v0​, correcting a printed slip, and the running-time claim is not stated.

The assortments S^il\hat S_{il}S^il​ enter Theorem 12 as an arbitrary family with the two properties of p. 28; the theorem is stated for every such family. A level at which (15) has no feasible assortment imposes nothing on S^il\hat S_{il}S^il​, as on the page. Requiring S^il\hat S_{il}S^il​ to lie in its level unconditionally would make the hypothesis unsatisfiable for most instances and the goal vacuous; that encoding, and any goal that assumes displays (34)–(36) or the exponent bounds, is ruled out.

Contributions welcome: proofs of the real-power milestones (3)–(5), which are self-contained, of the two cases (6)–(7), and of Proposition 15, whose fractional knapsack lemma (the greedy solution of a continuous knapsack sorted by ratio is optimal) is reusable beyond this mission.

Selected references

  • J. M. Davis, G. Gallego, H. Topaloglu, Assortment optimization under variants of the nested logit model, Operations Research 62(2), 2014; cited from the revised manuscript of June 18, 2013. https://doi.org/10.1287/opre.2014.1256
  • A. M. Frieze, M. R. B. Clarke, Approximation algorithms for the m-dimensional 0–1 knapsack problem: worst-case and probabilistic analyses, European Journal of Operational Research 15(1), 1984. (cited on p. 50 of the paper for the continuous knapsack (33); link not verified here)
  • M. R. Garey, D. S. Johnson, Computers and Intractability: A Guide to the Theory of NP-Completeness, W. H. Freeman, 1979.
11 thms1 active userReviewed
Linear OptimizationOperations ResearchProbability·Captain: mikedeng1

A Re-solving Heuristic with Uniformly Bounded Loss for Network Revenue Management 2: Frequent Re-solving Loses at Most O(√T) Against the Deterministic LP, Uniformly in the CapacitiesResearch Paper

Motivation

Network revenue management is the problem of selling several limited, perishable resources (seats on flight legs, hotel room-nights, machine hours) to customers who arrive over time and each want a fixed bundle of them. A seller who accepts every request early may run out of the resources that the most valuable later customers need, and a seller who is too cautious leaves capacity unsold at the end of the horizon. Airlines, hotels and car-rental firms solve instances of this problem daily (Talluri and van Ryzin, 2004).

The exact optimal policy is a dynamic program over the vector of remaining capacities, which is intractable for realistic networks. The standard remedy replaces the random demand by its mean, which gives a linear program, the deterministic LP (DLP), and turns its solution into an admission rule. Its value is an upper bound on what any policy can earn (Gallego and van Ryzin, 1997). Solving the LP once, at time zero, loses O(T)O(\sqrt{T})O(T​) against that bound over a horizon of length TTT. A natural improvement is to re-solve the LP as capacity is consumed.

Timeline of the re-solving question, as the source surveys it (Sec. 1, pp. 3–4):

  • 1997: Gallego and van Ryzin: static policies built from the DLP lose Θ(T)\Theta(\sqrt{T})Θ(T​).
  • 2002: Cooper gives an example in which re-solving the DLP makes booking-limit control worse.
  • 2008: Reiman and Wang re-solve exactly once, at an endogenous random time, with probabilistic allocation, and obtain o(T)o(\sqrt{T})o(T​) loss.
  • 2012: Jasin and Kumar prove that re-solving after every unit of time with probabilistic allocation (the policy FR below) loses O(1)O(1)O(1), provided the DLP solution is nondegenerate. Wu et al. (2015) obtain O(1)O(1)O(1) for one resource, with a constant that blows up as the solution approaches degeneracy.
  • 2018: Bumpensanti and Wang (arXiv:1802.06192) show that FR can lose Ω(T)\Omega(\sqrt{T})Ω(T​) against the hindsight optimum on a degenerate instance, propose a modified policy with O(1)O(1)O(1) loss in all cases, and prove that FR never loses more than O(T)O(\sqrt{T})O(T​) against the DLP. This mission formalizes the last of these results.

Setting

There are nnn customer classes j∈[n]j \in [n]j∈[n] and mmm resources l∈[m]l \in [m]l∈[m]. Over the horizon [0,T][0, T][0,T], class-jjj customers arrive according to independent Poisson processes with rates λj>0\lambda_j > 0λj​>0. Accepting a class-jjj customer earns rj≥0r_j \ge 0rj​≥0 and consumes alj≥0a_{lj} \ge 0alj​≥0 units of each resource lll; A=(alj)A = (a_{lj})A=(alj​) is the bill-of-materials matrix and AjA_jAj​ its jjj-th column. The initial capacity vector is C≥0C \ge 0C≥0. A customer can be accepted only if Aj≤C′A_j \le C'Aj​≤C′ componentwise, where C′C'C′ is the capacity remaining at its arrival.

For a capacity-per-unit-time vector b≥0b \ge 0b≥0, let

v(b)=max⁡{∑j=1nrjxj : ∑j=1nAjxj≤b, 0≤xj≤λj}.v(b) = \max\Big\{ \sum_{j=1}^n r_j x_j \ :\ \sum_{j=1}^n A_j x_j \le b,\ 0 \le x_j \le \lambda_j \Big\}.v(b)=max{j=1∑n​rj​xj​ : j=1∑n​Aj​xj​≤b, 0≤xj​≤λj​}.

The DLP value is vDLP(T,C)=T v(C/T)v^{\mathrm{DLP}}(T, C) = T\, v(C/T)vDLP(T,C)=Tv(C/T).

The Frequent Re-solving policy (FR) divides the horizon into TTT unit periods [t,t+1)[t, t+1)[t,t+1). At the start of period ttt, with remaining capacity C(t)C(t)C(t), it sets b(t)=C(t)/(T−t)b(t) = C(t)/(T - t)b(t)=C(t)/(T−t), computes an optimal solution x(t)x(t)x(t) of the LP with right-hand side b(t)b(t)b(t), and during the period accepts each class-jjj arrival with probability xj(t)/λjx_j(t)/\lambda_jxj​(t)/λj​, subject to the capacity check. Its expected revenue is vFR(T,C)v^{\mathrm{FR}}(T, C)vFR(T,C). The LP may have several optimal solutions; the results hold for every rule choosing among them.

Formalization targets

Goal: Proposition 3

There is a constant MMM, depending only on λ\lambdaλ, rrr and AAA, such that for every integer T≥1T \ge 1T≥1, every capacity C≥0C \ge 0C≥0 and every choice of optimal LP solutions,

vDLP(T,C)−vFR(T,C)≤MT.v^{\mathrm{DLP}}(T, C) - v^{\mathrm{FR}}(T, C) \le M \sqrt{T}.vDLP(T,C)−vFR(T,C)≤MT​.

The content is the uniformity: MMM does not depend on CCC, so the bound holds whether capacity is scarce, abundant, or degenerate for the LP.

Milestones, in the order the proof uses them

  1. Eq. (32), p. 35. The capacity check costs FR at most ∑jrj(log⁡T+1)+∑jrjλj\sum_j r_j(\log T + 1) + \sum_j r_j \lambda_j∑j​rj​(logT+1)+∑j​rj​λj​ relative to the LP revenue it targets:
vFR≥E[∑t=0T−1∑j=1nrjxj(t)]−∑j=1nrj(log⁡T+1)−∑j=1nrjλj.v^{\mathrm{FR}} \ge \mathbb{E}\Big[\sum_{t=0}^{T-1} \sum_{j=1}^n r_j x_j(t)\Big] - \sum_{j=1}^n r_j (\log T + 1) - \sum_{j=1}^n r_j \lambda_j .vFR≥E[t=0∑T−1​j=1∑n​rj​xj​(t)]−j=1∑n​rj​(logT+1)−j=1∑n​rj​λj​.
  1. Eq. (33), p. 35. LP sensitivity in the right-hand side: v(b)−v(b′)≤∑lrmax⁡l(bl−bl′)+v(b) - v(b') \le \sum_l r^l_{\max} (b_l - b'_l)^+v(b)−v(b′)≤∑l​rmaxl​(bl​−bl′​)+ with rmax⁡l=max⁡jrjI(alj>0)/aljr^l_{\max} = \max_j r_j \mathbb{I}(a_{lj} > 0)/a_{lj}rmaxl​=maxj​rj​I(alj​>0)/alj​.
  2. Lemma 8, p. 43 (corrected). E[(bl−bl(t))+]≤Kl∑i=0t−1(T−i−1)−2\mathbb{E}[(b_l - b_l(t))^+] \le K_l \sqrt{\sum_{i=0}^{t-1} (T-i-1)^{-2}}E[(bl​−bl​(t))+]≤Kl​∑i=0t−1​(T−i−1)−2​ with Kl=∑jalj2λjK_l = \sqrt{\sum_j a_{lj}^2 \lambda_j}Kl​=∑j​alj2​λj​​.
  3. p. 36. ∑t=0T−1∑i=0t−1(T−i−1)−2≤2T+2\sum_{t=0}^{T-1} \sqrt{\sum_{i=0}^{t-1} (T-i-1)^{-2}} \le 2\sqrt{T} + \sqrt{2}∑t=0T−1​∑i=0t−1​(T−i−1)−2​≤2T​+2​.
  4. p. 36, explicit bound.
vDLP−vFR≤∑l=1mrmax⁡lKl(2T+2)+n rmax⁡(log⁡T+1)+n rmax⁡λmax⁡.v^{\mathrm{DLP}} - v^{\mathrm{FR}} \le \sum_{l=1}^m r^l_{\max} K_l (2\sqrt{T} + \sqrt{2}) + n\, r_{\max} (\log T + 1) + n\, r_{\max} \lambda_{\max} .vDLP−vFR≤l=1∑m​rmaxl​Kl​(2T​+2​)+nrmax​(logT+1)+nrmax​λmax​.

Significance

The result. Because vDLPv^{\mathrm{DLP}}vDLP bounds the expected revenue of every admissible policy, Proposition 3 says that FR loses at most O(T)O(\sqrt{T})O(T​) against the optimal policy in every instance, degenerate or not. Re-solving therefore never does worse, in order, than solving once. Together with the Ω(T)\Omega(\sqrt{T})Ω(T​) lower bound on a degenerate instance (Proposition 2 of the same paper), it shows that the order T\sqrt{T}T​ is exact for FR. The nondegeneracy assumption of the earlier O(1)O(1)O(1) analysis cannot be dropped. The explicit form (milestone 5) bounds the loss by constants computable from (λ,r,A)(\lambda, r, A)(λ,r,A).

Formalizing it. The result is proved in the source; no part of it is machine-checked. A formalization adds three things. It gives a precise model of an adaptive, randomized admission policy in a Poisson network, which other results on re-solving and bid-price policies can reuse. It gives a checked LP sensitivity bound. And it checks the paper's constants: the printed constant of Lemma 8 is wrong (see Formalization scope), and the mission states the corrected one.

Difficulty

The obvious argument compares FR with the DLP period by period: the LP that FR solves at time ttt differs from the original only in its right-hand side, so the loss should be controlled by how far b(t)b(t)b(t) drifts below b=C/Tb = C/Tb=C/T. Two things break a naive version of this. First, b(t)b(t)b(t) is a ratio whose denominator T−tT - tT−t shrinks to 111, so fluctuations late in the horizon are amplified; a crude bound on the drift at each ttt sums to more than T\sqrt{T}T​. Second, FR's realized consumption is not the LP's target: the capacity check rejects customers, and the LP solutions x(i)x(i)x(i) depend on the whole past, so the consumption in different periods is not independent.

Uniformity in CCC is the whole point. Arguments that rely on a margin between C/TC/TC/T and the degenerate points of the LP, as in the nondegenerate analysis, give constants that blow up as that margin vanishes.

Formalization scope

The source is arXiv:1802.06192v3, whose printed page numbers equal the PDF page numbers. Everything lives in the namespace ResolvingNRM.FRUpper.

  • LP value. v(b)v(b)v(b) is the published piValue A r b lam of RLPBidPrice.Unbiased.Model (a real supremum, equal to the LP maximum for b≥0b \ge 0b≥0).
  • Optimal solutions. The "arg⁡max⁡\arg\maxargmax" of Algorithm 2 is an arbitrary optimal-solution selector sel; every theorem quantifies over all selectors, with constants chosen before the selector.
  • Randomness. Within a period the acceptance probabilities are fixed, so the period's arrivals are represented exactly as a Poisson number of customers in arrival order, with i.i.d. classes of law λj/∑iλi\lambda_j / \sum_i \lambda_iλj​/∑i​λi​ and independent Bernoulli acceptance coins. Expectations are series over this law (windowExp). FR's value is a backward recursion over periods (frTail), and expectations of functions of C(t)C(t)C(t) are a forward recursion (frStateExp).
  • Horizon. TTT is a positive integer; capacities are real vectors.
  • Standing assumptions of p. 7, left implicit there and hypotheses here: λj>0\lambda_j > 0λj​>0, rj≥0r_j \ge 0rj​≥0, alj≥0a_{lj} \ge 0alj​≥0, C≥0C \ge 0C≥0.
  • O(⋅)O(\cdot)O(⋅). "=O(T)= O(\sqrt{T})=O(T​)" is ∃M, ∀ sel,T≥1,C≥0\exists M,\ \forall\, \mathrm{sel}, T \ge 1, C \ge 0∃M, ∀sel,T≥1,C≥0. The paper's threshold T1T_1T1​ is dropped, which is equivalent because 0≤vFR0 \le v^{\mathrm{FR}}0≤vFR and vDLP≤T∑jrjλjv^{\mathrm{DLP}} \le T\sum_j r_j\lambda_jvDLP≤T∑j​rj​λj​.
  • Corrected slip, Lemma 8. The printed Kl=∑jalj2λj2K_l = \sqrt{\sum_j a_{lj}^2 \lambda_j^2}Kl​=∑j​alj2​λj2​​ is false: the proof replaces the conditional variance xj(i)≤λjx_j(i) \le \lambda_jxj​(i)≤λj​ of a Poisson increment by λj2\lambda_j^2λj2​. With one class, one resource, a=1a = 1a=1, λ=0.01\lambda = 0.01λ=0.01, C=λTC = \lambda TC=λT, T=1000T = 1000T=1000 and t=2t = 2t=2, an exact computation gives E[(b−b(2))+]≈1.96⋅10−5\mathbb{E}[(b - b(2))^+] \approx 1.96 \cdot 10^{-5}E[(b−b(2))+]≈1.96⋅10−5, above the printed bound 1.42⋅10−51.42 \cdot 10^{-5}1.42⋅10−5. The mission uses Kl=∑jalj2λjK_l = \sqrt{\sum_j a_{lj}^2 \lambda_j}Kl​=∑j​alj2​λj​​ in Lemma 8 and in the explicit bound.
  • Index slip. The paper's sums ∑l=1L\sum_{l=1}^L∑l=1L​ run over its mmm resources and are read as l∈[m]l \in [m]l∈[m].

A trivializing formalization is ruled out. The constant MMM is chosen before CCC. FR is defined with every optimal LP solution, not only nondegenerate or vertex ones. The capacity check is kept arrival by arrival; without it, FR's expected revenue would equal the LP revenue it targets and the gap would be trivially small.

Contributions welcome on every milestone. The LP sensitivity bound (33) is a statement about linear programs alone and milestone 4 about real numbers alone; both are reusable outside this mission, as is the window model of a randomized admission policy.

Selected references

  • Y. Bumpensanti, H. Wang, A Re-solving Heuristic with Uniformly Bounded Loss for Network Revenue Management, arXiv:1802.06192v3, 2018 (the source of this mission; later in Management Science 66(7), 2020). https://arxiv.org/abs/1802.06192
  • G. Gallego, G. van Ryzin, A Multiproduct Dynamic Pricing Problem and Its Applications to Network Yield Management, Operations Research 45(1), 24–41, 1997.
  • W. L. Cooper, Asymptotic Behavior of an Allocation Policy for Revenue Management, Operations Research 50(4), 720–727, 2002.
  • M. I. Reiman, Q. Wang, An Asymptotically Optimal Policy for a Quantity-Based Network Revenue Management Problem, Mathematics of Operations Research 33(2), 257–282, 2008.
  • S. Jasin, S. Kumar, A Re-Solving Heuristic with Bounded Revenue Loss for Network Revenue Management with Customer Choice, Mathematics of Operations Research 37(2), 313–345, 2012.
  • K. T. Talluri, G. J. van Ryzin, The Theory and Practice of Revenue Management, Springer, 2004.

Bibliographic details of the non-arXiv entries are those of the source's reference list (pp. 24–25).

9 thms1 active userReviewed
Previous

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