Prove2Me
Navigate
DiscoverCollectionsFormalpediaBlogsUsersMomentumMy Missions+
Prove2Me
⌕
Log in

Get started

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

Revenue Management and Pricing

6 missions · 2 completed

Missions

Open4Completed2All6
Bandit AlgorithmsOperations ResearchProbability·Captain: iris

Continuous Assortment Optimization with Logit Choice Probabilities under Incomplete Information: The Stochastic Approximation Policy SAP(α, β) Has Worst-Case Regret O(log T)Research Paper

Motivation

A retailer choosing which products to display, a platform choosing which items to recommend, and a broadcaster choosing which content to stream all face an assortment problem: offer a set of products, observe what customers buy, and earn the revenue of the purchased item. The standard demand model is the multinomial logit (MNL) model, under which the optimal assortment without capacity constraints consists of the most profitable products (Talluri and Van Ryzin, 2004). When the customers' preferences are unknown, the seller must learn them from sales, and the quality of a learning policy is measured by its regret: the cumulative revenue lost relative to an omniscient seller.

When the product catalogue is very large, or products are described by a continuous attribute such as quality, the discrete model becomes a continuum. Peeters, den Boer and Mandjes (arXiv:1807.06497, published in Operations Research, 2022) introduced a continuous analogue of the MNL model and studied learning in it. Their first main result, the goal of this mission, is that a simple stochastic approximation policy has worst-case regret of order log⁡T\log TlogT over TTT customers.

  • 1951: Robbins and Monro introduce stochastic approximation for finding the root of a function observed with noise (doi:10.1214/aoms/1177729586).
  • 2010: Rusmevichientong, Shen and Shmoys give adaptive policies for capacitated MNL assortment (doi:10.1287/opre.1100.0866).
  • 2011: Broadie, Cicek and Zeevi give mean-square rates for Kiefer–Wolfowitz type stochastic approximation (doi:10.1287/opre.1110.0970).
  • 2019: Agrawal, Avadhanula, Goyal and Zeevi give an O~(NT)\tilde O(\sqrt{NT})O~(NT​)-regret UCB policy for discrete MNL bandits (doi:10.1287/opre.2018.1832).
  • 2021/2022: Peeters, den Boer and Mandjes prove O(log⁡T)O(\log T)O(logT) worst-case regret for the continuum model, with a matching Ω(log⁡T)\Omega(\log T)Ω(logT) lower bound.

Setting

The products are the points x∈[0,1]x\in[0,1]x∈[0,1], and the revenue function w:[0,1]→[0,1]w:[0,1]\to[0,1]w:[0,1]→[0,1] is continuously differentiable with w′≥k>0w'\ge k>0w′≥k>0. An assortment is a Borel set S⊆[0,1]S\subseteq[0,1]S⊆[0,1]; S\mathcal SS is the set of all of them. A customer offered SSS buys a product in a Borel set A⊆SA\subseteq SA⊆S, or buys nothing (∅\emptyset∅), with probabilities

P(XS∈A)=∫Av(x) dx1+∫Sv(x) dx,P(XS=∅)=11+∫Sv(x) dx,\mathbb P(X^S\in A)=\frac{\int_A v(x)\,dx}{1+\int_S v(x)\,dx},\qquad \mathbb P(X^S=\emptyset)=\frac{1}{1+\int_S v(x)\,dx},P(XS∈A)=1+∫S​v(x)dx∫A​v(x)dx​,P(XS=∅)=1+∫S​v(x)dx1​,

where the preference function vvv is unknown to the seller. The expected revenue of SSS is r(S,v)=∫Svw/(1+∫Sv)r(S,v)=\int_S vw\big/(1+\int_S v)r(S,v)=∫S​vw/(1+∫S​v), and ϱ∗=max⁡S∈Sr(S,v)\varrho^*=\max_{S\in\mathcal S}r(S,v)ϱ∗=maxS∈S​r(S,v).

The class V(v‾,v‾,D)\mathcal V(\underline v,\overline v,D)V(v​,v,D) (Assumption 1) contains the integrable vvv with v‾≤v≤v‾\underline v\le v\le\overline vv​≤v≤v on [0,1][0,1][0,1] that are differentiable on (0,1)(0,1)(0,1) with ∣v′∣≤D|v'|\le D∣v′∣≤D, where v‾>v‾>0\overline v>\underline v>0v>v​>0 and v‾>w(0)/∫01(w−w(0))\overline v>w(0)/\int_0^1(w-w(0))v>w(0)/∫01​(w−w(0)).

The generalized inverse is w−1(ϱ)=min⁡{x∈[0,1]:w(x)≥ϱ}w^{-1}(\varrho)=\min\{x\in[0,1]:w(x)\ge\varrho\}w−1(ϱ)=min{x∈[0,1]:w(x)≥ϱ}, and h(ϱ)=r([w−1(ϱ),1],v)h(\varrho)=r([w^{-1}(\varrho),1],v)h(ϱ)=r([w−1(ϱ),1],v). The policy SAP(α,β)(\alpha,\beta)(α,β) starts at ϱ1∈[0,1]\varrho_1\in[0,1]ϱ1​∈[0,1], offers St=[w−1(ϱt),1]S_t=[w^{-1}(\varrho_t),1]St​=[w−1(ϱt​),1] in period ttt, observes the revenue Rt=w(Xt)1{Xt∈St}R_t=w(X_t)\mathbf 1\{X_t\in S_t\}Rt​=w(Xt​)1{Xt​∈St​}, and updates

ϱt+1=ϱt+αt+β (Rt−ϱt).\varrho_{t+1}=\varrho_t+\frac{\alpha}{t+\beta}\,(R_t-\varrho_t).ϱt+1​=ϱt​+t+βα​(Rt​−ϱt​).

Given the history, the purchase XtX_tXt​ has the law of XStX^{S_t}XSt​. The regret over TTT periods is Δπ(T,v)=∑t=1TEπ[ϱ∗−r(St,v)]\Delta_\pi(T,v)=\sum_{t=1}^T\mathbb E_\pi[\varrho^*-r(S_t,v)]Δπ​(T,v)=∑t=1T​Eπ​[ϱ∗−r(St​,v)].

Formalization targets

Goal: Theorem 1 (p.12)

For α≥v‾+1\alpha\ge\overline v+1α≥v+1 and β≥max⁡{0,α−1}\beta\ge\max\{0,\alpha-1\}β≥max{0,α−1} there is C‾>0\overline C>0C>0 such that for all T≥2T\ge2T≥2

sup⁡v∈VΔπ(T,v)≤C‾log⁡T.\sup_{v\in\mathcal V}\Delta_\pi(T,v)\le\overline C\log T.v∈Vsup​Δπ​(T,v)≤ClogT.

The constant depends on www, the class and the policy, but not on TTT or vvv.

Milestones, in proof order

  1. Well-posedness: an SAP purchase process exists for every horizon.
  2. §4.1: the optimum over S\mathcal SS is attained by an interval [w−1(ϱ),1][w^{-1}(\varrho),1][w−1(ϱ),1].
  3. Lemma 4: h(ϱ∗)=ϱ∗h(\varrho^*)=\varrho^*h(ϱ∗)=ϱ∗; the drift inequality (23) (ϱ−ϱ∗)(h(ϱ)−ϱ)≤−(ϱ−ϱ∗)2/(1+v‾)(\varrho-\varrho^*)(h(\varrho)-\varrho)\le-(\varrho-\varrho^*)^2/(1+\overline v)(ϱ−ϱ∗)(h(ϱ)−ϱ)≤−(ϱ−ϱ∗)2/(1+v); the derivative g′=(g−w)ξg'=(g-w)\xig′=(g−w)ξ and the curvature bound (27); the quadratic gap (24) h(ϱ∗)−h(ϱ)≤C(ϱ−ϱ∗)2h(\varrho^*)-h(\varrho)\le C(\varrho-\varrho^*)^2h(ϱ∗)−h(ϱ)≤C(ϱ−ϱ∗)2 with CCC uniform over V\mathcal VV.
  4. The iterates stay in [0,1][0,1][0,1], and E[Rt∣history]=h(ϱt)\mathbb E[R_t\mid\text{history}]=h(\varrho_t)E[Rt​∣history]=h(ϱt​).
  5. The recursion (25) δt+1≤δt(1−21+v‾αt+β)+α2(t+β)2\delta_{t+1}\le\delta_t\big(1-\frac{2}{1+\overline v}\frac{\alpha}{t+\beta}\big)+\frac{\alpha^2}{(t+\beta)^2}δt+1​≤δt​(1−1+v2​t+βα​)+(t+β)2α2​ for δt=E[(ϱt−ϱ∗)2]\delta_t=\mathbb E[(\varrho_t-\varrho^*)^2]δt​=E[(ϱt​−ϱ∗)2], and Lemma 5, δt≤κ/(t+β)\delta_t\le\kappa/(t+\beta)δt​≤κ/(t+β).
  6. The regret identity Δπ(T,v)=∑tE[h(ϱ∗)−h(ϱt)]\Delta_\pi(T,v)=\sum_t\mathbb E[h(\varrho^*)-h(\varrho_t)]Δπ​(T,v)=∑t​E[h(ϱ∗)−h(ϱt​)] and ∑t=1T(t+β)−1≤3log⁡T\sum_{t=1}^T(t+\beta)^{-1}\le3\log T∑t=1T​(t+β)−1≤3logT.

Significance

The theorem shows that in the continuum logit model the seller does not need to estimate the preference function: a one-dimensional Robbins–Monro recursion on the observed revenue finds the optimal revenue level, and the loss is logarithmic in the horizon. This is faster than the T\sqrt TT​ rates typical of continuum-armed bandits, and the paper's Theorem 2 shows it is optimal. The bound is worst-case over a class of preference functions, which matters when nothing is known about the customers beyond the class.

The result is proved on paper. No part of it is formalized on the platform. The mission adds the platform's first continuum-indexed choice model, a measure-theoretic model of a sequential learning policy whose decisions depend on the full purchase history, and a complete regret analysis for a stochastic approximation policy. The discrete MNL bandit of Agrawal et al. is on the platform as MNLBandit.UCB; it is a related, not a reused, model.

Difficulty

The policy itself is three lines; the difficulty is in the analysis of the map hhh. The obvious route — "hhh is a contraction, so the iterates converge" — fails: hhh is not a contraction, and ϱ↦h(ϱ)\varrho\mapsto h(\varrho)ϱ↦h(ϱ) is flat below w(0)w(0)w(0) and drops to zero above w(1)w(1)w(1). The proof instead needs the one-sided drift inequality (23), whose constant must use the upper bound v‾\overline vv, and the quadratic gap (24), which requires a second-order Taylor expansion of g(y)=r([y,1],v)g(y)=r([y,1],v)g(y)=r([y,1],v) with a curvature bound (27) that is uniform over the class. The bound (27) needs vvv differentiable, and the uniformity of CCC and κ\kappaκ over V\mathcal VV is what turns per-instance bounds into a worst-case one.

The induction in the proof of Lemma 5 multiplies by 1−2α/((1+v‾)(t+β))1-2\alpha/((1+\overline v)(t+\beta))1−2α/((1+v)(t+β)), which is negative for small ttt when v‾\overline vv is small; the printed κ=max⁡{1+β,α(1+v‾)}\kappa=\max\{1+\beta,\alpha(1+\overline v)\}κ=max{1+β,α(1+v)} does not follow from the printed step there, although the lemma's statement (some κ>0\kappa>0κ>0) remains true.

Formalization scope

Products are reals and purchases live in ℝ ⊕ Unit (Sum.inr () is no-purchase). w,v:R→Rw,v:\mathbb R\to\mathbb Rw,v:R→R carry hypotheses only on [0,1][0,1][0,1] and (0,1)(0,1)(0,1). S\mathcal SS is the measurable subsets of [0,1][0,1][0,1] (the uncapacitated case c=1c=1c=1), and ϱ∗\varrho^*ϱ∗ is a supremum (sSup), shown to be a maximum in a milestone. The paper's min⁡{x:w(x)≥ϱ}\min\{x:w(x)\ge\varrho\}min{x:w(x)≥ϱ} is empty for ϱ>w(1)\varrho>w(1)ϱ>w(1), which SAP can reach since ϱ1∈[0,1]\varrho_1\in[0,1]ϱ1​∈[0,1] is arbitrary; there w−1(ϱ):=1w^{-1}(\varrho):=1w−1(ϱ):=1, so the offered assortment {1}\{1\}{1} is a null set. Lean's sInf ∅ = 0 would instead offer [0,1][0,1][0,1] and break (23). The class has a fixed derivative bound DDD, since Assumption 1(ii)'s supremum over all admissible vvv is otherwise infinite.

Time is 0-based (n=t−1n=t-1n=t−1). A horizon-TTT process is a probability space with measurable purchases X0,…,XT−1X_0,\dots,X_{T-1}X0​,…,XT−1​ whose conditional law given the history is the choice law of the offered assortment, written as an identity of measures on history events. A Lean proof that the SAP iterate is a measurable function of the history exists, so the expectations in the regret are genuine integrals. The goal quantifies over all such processes, and a milestone asserts that one exists, so the goal is not vacuous.

A complete development needs: second-order calculus for parameter integrals y↦∫y1y\mapsto\int_y^1y↦∫y1​; the construction of a process from a history-dependent kernel (Ionescu–Tulcea, ProbabilityTheory.Kernel.traj, or an iterated Measure.bind); and conditional expectations given a history σ-algebra. The mean-square recursion and the harmonic-sum bound are reusable for other Robbins–Monro analyses. Contributions to any milestone are welcome, in any order; milestones 2–3 and 6 are independent of the process construction.

Selected references

  • Y. Peeters, A. V. den Boer, M. Mandjes, Continuous Assortment Optimization with Logit Choice Probabilities under Incomplete Information, arXiv:1807.06497v4, 2021; Operations Research, 2022. https://arxiv.org/abs/1807.06497
  • H. Robbins, S. Monro, A Stochastic Approximation Method, Annals of Mathematical Statistics 22(3), 1951. https://doi.org/10.1214/aoms/1177729586
  • 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
  • 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), 2010. https://doi.org/10.1287/opre.1100.0866
  • M. Broadie, D. Cicek, A. Zeevi, General Bounds and Finite-Time Improvement for the Kiefer-Wolfowitz Stochastic Approximation Algorithm, Operations Research 59(5), 2011. https://doi.org/10.1287/opre.1110.0970
  • S. Agrawal, V. Avadhanula, V. Goyal, A. Zeevi, MNL-Bandit: A Dynamic Learning Approach to Assortment Selection, Operations Research 67(5), 2019. https://doi.org/10.1287/opre.2018.1832
17 thms2 active usersReviewed
Operations ResearchOptimizationTheoretical Computer Science·Captain: iris

Assortment Optimization with Multi-Item Basket Purchase under Multivariate MNL Model: Without Interaction Terms, Algorithm U1 Finds a (1 − ε)-Optimal Assortment in Polynomial TimeResearch Paper

Motivation

Retailers rarely sell one item at a time. A shopper fills a basket with, say, one book, one film and one game, and the decision to buy in one category is linked to the others: if nothing in the store is attractive enough, the shopper leaves without buying anything. Assortment optimization asks which products a retailer should offer so as to maximize expected revenue under a model of customer choice. For single-item purchases under the multinomial logit (MNL) model the answer is classical: some revenue-ordered set is optimal (Talluri and van Ryzin 2004). For basket purchases the natural choice model in the marketing literature is the multivariate MNL (MVMNL) model (Russell and Petersen 2000), and its assortment problem was open.

Jasin, Lyu, Najafi and Zhang (SSRN 3818886; Manufacturing & Service Operations Management, 2023) study this problem. Their first result is for the case without interaction terms across groups. The problem is NP-hard there (by a reduction from PARTITION), but it has a fully polynomial-time approximation scheme (FPTAS), Algorithm U1. This mission formalizes that FPTAS and its analysis.

Setting

There are JJJ groups of products and III products in each group; (i,j)(i,j)(i,j) is product iii of group jjj. A customer buys a bundle y=(yj)j=1Jy=(y_j)_{j=1}^Jy=(yj​)j=1J​, at most one product per group: yj=iy_j=iyj​=i buys product iii of group jjj, and yj=0y_j=0yj​=0 buys nothing there. An assortment x∈{0,1}JIx\in\{0,1\}^{JI}x∈{0,1}JI offers product (i,j)(i,j)(i,j) iff xi,j=1x_{i,j}=1xi,j​=1, and a bundle is permissible when it buys only offered products.

Each product has a weight vi,j=eθi,j>0v_{i,j}=e^{\theta_{i,j}}>0vi,j​=eθi,j​>0 and a revenue ri,j≥0r_{i,j}\ge0ri,j​≥0. The empty bundle has weight v0=eμ0>1v_0=e^{\mu_0}>1v0​=eμ0​>1, the value of the global outside option. Without interaction terms, a nonempty bundle has weight ∏j:yj≠0vyj,j\prod_{j:y_j\neq0}v_{y_j,j}∏j:yj​=0​vyj​,j​. A bundle is chosen with probability proportional to its weight among the permissible bundles, so the expected total revenue is

R(x)=∑y permissibler(y) eμ(y)∑y permissibleeμ(y),r(y)=∑j:yj≠0ryj,j.R(x)=\frac{\sum_{y\ \text{permissible}}r(y)\,e^{\mu(y)}}{\sum_{y\ \text{permissible}}e^{\mu(y)}},\qquad r(y)=\sum_{j:y_j\neq0}r_{y_j,j}.R(x)=∑y permissible​eμ(y)∑y permissible​r(y)eμ(y)​,r(y)=j:yj​=0∑​ryj​,j​.

Within each group, the products are ordered by non-increasing revenue. A revenue-ordered assortment offers a prefix {1,…,i(j)}\{1,\dots,i(j)\}{1,…,i(j)} of every group. Write γj(xj)=1+∑ivi,jxi,j\gamma_j(x_j)=1+\sum_i v_{i,j}x_{i,j}γj​(xj​)=1+∑i​vi,j​xi,j​, Γ(x)=∏jγj(xj)\Gamma(x)=\prod_j\gamma_j(x_j)Γ(x)=∏j​γj​(xj​) and β=max⁡j(1+∑ivi,j)\beta=\max_j(1+\sum_i v_{i,j})β=maxj​(1+∑i​vi,j​).

Formalization targets

Goal: Theorem 1 (p.14)

There is a universal constant CCC such that for all J,I≥1J,I\ge1J,I≥1, every instance and every ϵ∈(0,1)\epsilon\in(0,1)ϵ∈(0,1), Algorithm U1 run with ϵ′=ϵ/(2J)\epsilon'=\epsilon/(2J)ϵ′=ϵ/(2J) returns xˉ∗\bar x^*xˉ∗ with

R(xˉ∗)≥(1−ϵ) R(x)  for every x,#operations≤C J5 I (1+log⁡β)2ϵ2.R(\bar x^*)\ge(1-\epsilon)\,R(x)\ \text{ for every } x,\qquad \#\text{operations}\le C\,\frac{J^5\,I\,(1+\log\beta)^2}{\epsilon^2}.R(xˉ∗)≥(1−ϵ)R(x)  for every x,#operations≤Cϵ2J5I(1+logβ)2​.

The paper prints log⁡2β\log^2\betalog2β. That bound fails as β→1\beta\to1β→1, because the algorithm always performs at least one operation. The goal uses (1+log⁡β)2(1+\log\beta)^2(1+logβ)2, which agrees with the printed bound up to a factor 4 when β≥e\beta\ge eβ≥e.

Milestones, in proof order

  1. Lemma 1 and (5): Pi,j(x)=vi,jxi,j∏j′≠jγj′/(v0−1+Γ(x))P_{i,j}(x)=v_{i,j}x_{i,j}\prod_{j'\neq j}\gamma_{j'}/(v_0-1+\Gamma(x))Pi,j​(x)=vi,j​xi,j​∏j′=j​γj′​/(v0​−1+Γ(x)), hence a closed form of R(x)R(x)R(x).
  2. (7): R(x)=∑jRj(x)R(x)=\sum_j R_j(x)R(x)=∑j​Rj​(x) with Rj(x)=∑iri,jvi,jxi,j/([v0−1Γ(x)+1]γj(xj))R_j(x)=\sum_i r_{i,j}v_{i,j}x_{i,j}\big/\big([\tfrac{v_0-1}{\Gamma(x)}+1]\gamma_j(x_j)\big)Rj​(x)=∑i​ri,j​vi,j​xi,j​/([Γ(x)v0​−1​+1]γj​(xj​)).
  3. (EC.1) and Lemma 2: with the other groups fixed, RRR is linear-fractional in xjx_jxj​, and some revenue-ordered assortment is optimal.
  4. (EC.5): rounding each γj\gamma_jγj​ up to a power of 1+ϵ′1+\epsilon'1+ϵ′ changes each RjR_jRj​ by at most (1+ϵ′)J−1(1+\epsilon')^{J-1}(1+ϵ′)J−1 one way and (1+ϵ′)J(1+\epsilon')^J(1+ϵ′)J the other.
  5. Algorithm U1, step 2 and EC.3 feasibility: the output maximizes the rounded objective Rˉ\bar RRˉ over revenue-ordered assortments, and K=JK=JK=J always gives a candidate.
  6. (EC.6) and the Taylor step: R(xˉ∗)≥R(x)/(1+ϵ′)2J≥(1−2Jϵ′)R(x)R(\bar x^*)\ge R(x)/(1+\epsilon')^{2J}\ge(1-2J\epsilon')R(x)R(xˉ∗)≥R(x)/(1+ϵ′)2J≥(1−2Jϵ′)R(x).
  7. The running time: O(J3M2I)O(J^3M^2I)O(J3M2I) operations, and M≤1+4Jlog⁡β/ϵM\le1+4J\log\beta/\epsilonM≤1+4Jlogβ/ϵ (corrected from O(2Jϵ−1log⁡β)O(2J\epsilon^{-1}\log\beta)O(2Jϵ−1logβ)).

Significance

The theorem shows that the problem, NP-hard even without interactions, can still be approximated to any accuracy in time polynomial in the input and 1/ϵ1/\epsilon1/ϵ. Its structural facts carry over to the rest of the paper. Lemma 1 turns a sum over (I+1)J(I+1)^J(I+1)J bundles into a product over groups, and Lemma 2 cuts the search to (I+1)J(I+1)^J(I+1)J revenue-ordered candidates. The rounding-plus-DP scheme is the template for the algorithms with pairwise group interactions and with a capacity constraint (Theorems 2–4).

The result is proved in the paper; nothing here is open mathematics. No machine-checked version is known. The formalization adds three things:

  • a reusable basket-choice model with a global outside option, the first multivariate choice model on the platform;
  • a fully specified version of Algorithm U1 with an explicit cost model;
  • corrected running-time bounds with explicit constants, and an explicit sign condition r≥0r\ge0r≥0 on revenues, which the paper leaves implicit but which Theorem 1 needs.

Difficulty

Without interactions, the revenue still does not separate across groups. The denominator v0−1+Γ(x)v_0-1+\Gamma(x)v0​−1+Γ(x) couples every group, so the natural idea of optimizing each group's threshold separately fails; the paper's Table 2 shows the loss can be large. A DP that guesses Γ(x∗)\Gamma(x^*)Γ(x∗) exactly has exponentially many states. Algorithm U1 rounds every γj\gamma_jγj​ to a power (1+ϵ′)mj(1+\epsilon')^{m_j}(1+ϵ′)mj​, so the DP state becomes the integer exponent k∈{0,…,JM}k\in\{0,\dots,JM\}k∈{0,…,JM}.

Formally, the hard parts are the following.

  • The DP correctness of step 2. One must show that the backward table over integer states, followed by the traceback, attains the maximum of Rˉ\bar RRˉ over revenue-ordered assortments with Γˉ=(1+ϵ′)K\bar\Gamma=(1+\epsilon')^KΓˉ=(1+ϵ′)K. One must also show that every revenue-ordered assortment falls under some enumerated K∈{J,…,JM}K\in\{J,\dots,JM\}K∈{J,…,JM}.
  • The operation count of an explicit program.
  • The finite-sum identity of Lemma 1, whose proof the paper omits.

Formalization scope

Groups are Fin J and products Fin I. Assortments are Fin J → Fin I → Bool, and bundles are Fin J → Option (Fin I), with none meaning nothing is bought from that group. The data are vvv, v0v_0v0​ and rrr themselves (θ=log⁡v\theta=\log vθ=logv, μ0=log⁡v0\mu_0=\log v_0μ0​=logv0​). The empty bundle weighs v0v_0v0​, and RRR is defined by the bundle formula (4), never by its closed forms (5) or (7), which are milestones. Interaction terms ψ\psiψ are part of the instance, and "no interaction terms" is a hypothesis. The standing assumptions are J,I≥1J,I\ge1J,I≥1, v>0v>0v>0, v0>1v_0>1v0​>1, r≥0r\ge0r≥0, and rrr non-increasing within each group.

Algorithm U1 is an explicit program algU1 in the unit-cost model CapacitatedAssort.General.Work, reused from the platform. The following operations are charged where they execute: +,−,×,/+,-,\times,/+,−,×,/, comparisons, ⌈⋅⌉+\lceil\cdot\rceil_+⌈⋅⌉+​, array reads and writes, and loop dispatch. The program evaluates no logarithm or exponential. It finds each rounding exponent by a multiplication loop, indexes DP states by integers, and breaks ties by the first maximizer. The universal constants CCC and C1C_1C1​ are quantified before JJJ, III, the instance and ϵ\epsilonϵ. This rules out a trivial reading in which an existential "some algorithm" could enumerate all 2JI2^{JI}2JI assortments.

Contributions are welcome on any milestone. The bundle-sum factorization (Lemma 1) and the DP-correctness lemma are the largest pieces. The model module is reusable for other basket-choice results.

Selected references

  • S. Jasin, C. Lyu, S. Najafi, H. Zhang, Assortment Optimization with Multi-Item Basket Purchase under Multivariate MNL Model, SSRN 3818886 (2023); Manufacturing & Service Operations Management, 2023. https://ssrn.com/abstract=3818886, https://doi.org/10.1287/msom.2021.0526
  • G. J. Russell, A. Petersen, Analysis of cross category dependence in market basket selection, Journal of Retailing 76(3), 2000. https://doi.org/10.1016/S0022-4359(00)00030-0
  • 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
17 thms2 active usersReviewed
Dynamic ProgrammingOperations ResearchOptimization·Captain: iris

Approximation Schemes for Capacity-Constrained Assortment Optimization Under the Nested Logit Model: A Dynamic Program Earns 1 − ε of Optimal in Time Polynomial in Input Size, 1/ε and γ⁺/γ⁻Research Paper

Motivation

The Nested Logit model is one of the most widely used discrete choice models in revenue management. It captures substitution that the Multinomial Logit model cannot: products are grouped into nests, a customer first chooses a nest and then a product in it (Davis, Gallego, Topaloglu, 2014; Gallego, Topaloglu, 2014). A retailer that decides which products to offer under a limit on shelf space or budget faces capacitated assortment optimization. Earlier constant-factor approximations depend on restricted forms of the model, such as dissimilarity parameters in [0,1][0,1][0,1] or no nest-level no-purchase option, and Segev (SSRN #3553264) removes those restrictions with an approximation scheme.

Setting

There are nnn products partitioned into nests N1,…,NMN_1,\dots,N_MN1​,…,NM​. Product iii has a price ri>0r_i>0ri​>0, a preference weight wi>0w_i>0wi​>0 and a cost cic_ici​. Nest mmm has a dissimilarity parameter γm\gamma_mγm​ and a no-purchase weight w0mw_{0m}w0m​, and w0w_0w0​ is the outside no-purchase weight. An assortment is S=(S1,…,SM)S=(S_1,\dots,S_M)S=(S1​,…,SM​) with Sm⊆NmS_m\subseteq N_mSm​⊆Nm​. For T⊆NmT\subseteq N_mT⊆Nm​ write w+(T)=I[T≠∅](w0m+w(T))w_+(T)=\mathbb I[T\ne\emptyset](w_{0m}+w(T))w+​(T)=I[T=∅](w0m​+w(T)) and Rm(T)=∑i∈Triwi/w+(T)\mathcal R_m(T)=\sum_{i\in T}r_iw_i/w_+(T)Rm​(T)=∑i∈T​ri​wi​/w+​(T). The expected revenue is

R(S)=∑mRm(Sm) w+(Sm)γmw0+∑mw+(Sm)γm.\mathcal R(S)=\frac{\sum_{m}\mathcal R_m(S_m)\,w_+(S_m)^{\gamma_m}}{w_0+\sum_{m}w_+(S_m)^{\gamma_m}} .R(S)=w0​+∑m​w+​(Sm​)γm​∑m​Rm​(Sm​)w+​(Sm​)γm​​.

An assortment is C\mathcal CC-feasible when ∑mc(Sm)≤C\sum_m c(S_m)\le\mathcal C∑m​c(Sm​)≤C. The objective is a C\mathcal CC-feasible assortment of maximum revenue.

Formalization targets

Goal

The goal is Theorem 1.1, stated about one explicit instrumented program fptas. There is a universal KKK such that, for every ϵ>0\epsilon>0ϵ>0: the program returns an assortment whenever a C\mathcal CC-feasible one exists; every assortment S~\tilde SS~ it returns is C\mathcal CC-feasible; and

R(S~)≥(1−ϵ) R(S)for every C-feasible S,\mathcal R(\tilde S)\ge(1-\epsilon)\,\mathcal R(S)\quad\text{for every }\mathcal C\text{-feasible }S,R(S~)≥(1−ϵ)R(S)for every C-feasible S,

and on rational instances it uses at most

K⋅inputSize20⋅max⁡(1,1/ϵ)8⋅(γ+/γ−)10K\cdot\mathrm{inputSize}^{20}\cdot\max(1,1/\epsilon)^{8}\cdot(\gamma^+/\gamma^-)^{10}K⋅inputSize20⋅max(1,1/ϵ)8⋅(γ+/γ−)10

operations, with γ+=max⁡{1,γ1,…,γM}\gamma^+=\max\{1,\gamma_1,\dots,\gamma_M\}γ+=max{1,γ1​,…,γM​} and γ−=min⁡{1,γ1,…,γM}\gamma^-=\min\{1,\gamma_1,\dots,\gamma_M\}γ−=min{1,γ1​,…,γM​}.

Milestones

In proof order: the representation (1); Lemma 2.1 on the dual; the single-nest results, Lemmas 4.1–4.3 and Theorem 3.1; the ranges of §3.3; the recursion (4) and the final selection, as properties of the program; Lemmas 3.3, 3.5, 3.6, 3.4 and 3.2; and the running time of §3.4.

Significance

The theorem shows that the revenue of the capacitated problem can be approached to any accuracy by combinatorial means, in the full generality of the model: parameters γm>1\gamma_m>1γm​>1, nest-level no-purchase weights, and arbitrary costs. It is the first fully polynomial-time approximation scheme for this problem when γ+/γ−\gamma^+/\gamma^-γ+/γ− is polynomial in the input size, a dependence the paper explains in §5.

The result is proved on paper and not machine-checked. This mission adds the first approximation scheme for the Nested Logit model on the platform, an instrumented dual dynamic program with a rounding-path analysis, and a formal count of its operations. It also records three defects of the printed construction with their counterexamples: the grids stop at the range of a single nest, the grid for ρ\rhoρ starts too high, and the single-nest search can return the empty set. The headline theorem survives with minimal changes to the algorithm, at the cost of log⁡M\log MlogM and constant factors. The Operations Research version was not consulted, so it is unknown whether it already contains these repairs.

Difficulty

The obvious dynamic program, tracking the sums ω\omegaω and ρ\rhoρ across nests, fails in two places. First, the single-nest problem with exact targets is NP-hard, so the targets must be relaxed to a neighbourhood; the difficulty is that w+(T)γw_+(T)^{\gamma}w+​(T)γ and Rk(T)w+(T)γ\mathcal R_k(T)w_+(T)^{\gamma}Rk​(T)w+​(T)γ are not additive over products, and their sensitivity changes direction at γ=1\gamma=1γ=1. Secondly, the grids are not closed under subtraction, so the recursion cannot split ω\omegaω exactly into ωk\omega_kωk​ and ω−ωk\omega-\omega_kω−ωk​. The recursion splits it superadditively and ρ\rhoρ subadditively, and a rounded path must be shown admissible at every nest, with distortion only (1+ξω)M(1+\xi_\omega)^{M}(1+ξω​)M and (1−ξρ)M(1-\xi_\rho)^{M}(1−ξρ​)M over MMM levels. Rounding one nest at a time, the natural first idea, accumulates errors in the wrong direction unless the factor (1+ξω)M−k(1+\xi_\omega)^{M-k}(1+ξω​)M−k decreases along the path.

Formalization scope

Products are Fin n, the partition is a map nest : Fin n → Fin M, and an assortment is a Finset (Fin n). Standing assumptions are ri,wi>0r_i,w_i>0ri​,wi​>0, M≥1M\ge1M≥1, nonempty nests, γm>0\gamma_m>0γm​>0, and w0,w0m≥0w_0,w_{0m}\ge0w0​,w0m​≥0; costs and C\mathcal CC are real of any sign. ϵ^=min⁡(ϵ,1)\hat\epsilon=\min(\epsilon,1)ϵ^=min(ϵ,1) enters the parameters, and ϵ\epsilonϵ the guarantees.

Without the program the accuracy claim is trivial, because the optimum itself exists. So the goal is about fptas, in the cost monad Work of CapacitatedAssort.General: it uses unit-cost arithmetic, comparisons, ⌊⋅⌋\lfloor\cdot\rfloor⌊⋅⌋, ⌈⋅⌉\lceil\cdot\rceil⌈⋅⌉, Real.rpow, Real.log and array accesses, each charged where it executes, and it enumerates no subset of products. This is the accounting of §3.4 and §4.2 of the paper, and exact powers are irrational for rational data, so no bit model is claimed for them.

The neighbourhood of (3) is stated multiplied out, with ∅\emptyset∅ only at ωk=ρk=0\omega_k=\rho_k=0ωk​=ρk​=0. The grids are corrected as described in the docstrings of Pomega and Lrho, the table of the single-nest step distinguishes nonempty subsets, and each difference is declared on its item.

A complete development needs estimates for real powers on both sides of γ=1\gamma=1γ=1, correctness and cost lemmas for the imperative tables, and bit-length bounds for the logarithmic factors (logs_le_inputSize). The cost lemmas for tabulate and for nested loops are reusable. Contributions to any milestone are welcome, and the analysis of the single-nest step can be done independently of the rest.

Selected references

  • D. Segev, Approximation Schemes for Capacity-Constrained Assortment Optimization under the Nested Logit Model, SSRN #3553264, March 2020. https://ssrn.com/abstract=3553264 (published as a Technical Note in Operations Research, 2022, https://doi.org/10.1287/opre.2022.2336; the numbering of this mission follows the SSRN version).
  • J. 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
  • G. Gallego, H. Topaloglu, Constrained assortment optimization for the nested logit model, Management Science 60(10), 2014. https://doi.org/10.1287/mnsc.2014.1931
24 thms2 active usersReviewed
EconomicsOperations Research·Captain: iris

Bounds and Heuristics for Multi-Product Personalized Pricing: Pricing Along a Positive Vector Earns at Least 1/(1+ln⁡(qmax⁡/qmin⁡))1/(1+\ln(q_{\max}/q_{\min}))1/(1+ln(qmax​/qmin​)) of the Optimal Personalized Profit, and This Is TightResearch Paper

Why pricing along one vector

Charging every customer segment its own optimal prices is the most a seller can extract with a posted price per segment, and solving each segment's problem is often easy. Finding one optimal non-personalized price vector is harder, because aggregate demand may be much more complex than the segments' demands, so single-factor heuristics are used instead. A natural question is how much is lost.

For one product, Bergemann, Castro and Weintraub bound the ratio of personalized to uniform profit by two when each type's profit is concave and demand is positive on a common compact set; Malueg and Snyder bound it by the number of types under milder conditions; and Elmachtoub, Gupta and Hamilton give bounds in terms of summary statistics of the demand distribution. For several products, Berbeglia and Joret analyse uniform pricing in the special case of envy-free pricing. Gallego and Berbeglia state the first tight bound of this kind for general multi-product demand, against the personalized optimum rather than the non-personalized one.

The pricing model

A firm sells a finite nonempty set NNN of products to customer types t∈Tt\in Tt∈T, distributed according to a probability measure μ\muμ. In the paper TTT is finite with weights θj>0\theta_j>0θj​>0; its remark on p.5 extends the result to continuous types. Write dt,i(p)≥0d_{t,i}(p)\ge0dt,i​(p)≥0 for the demand of type ttt for product iii at the price vector p≥0p\ge0p≥0. The type profit is Rt(p)=∑i∈Npidt,i(p)R_t(p)=\sum_{i\in N}p_id_{t,i}(p)Rt​(p)=∑i∈N​pi​dt,i​(p) and the aggregate profit is R(p)=∫Rt(p) dμ(t)R(p)=\int R_t(p)\,d\mu(t)R(p)=∫Rt​(p)dμ(t).

Let pˉt\bar p_tpˉ​t​ be an optimal price vector of type ttt (Assumption A0: optimal prices exist and are positive and finite). The personalized profit is Rˉ=∫Rt(pˉt) dμ(t)\bar{\mathcal R}=\int R_t(\bar p_t)\,d\mu(t)Rˉ=∫Rt​(pˉ​t​)dμ(t). Pricing along a positive vector fff earns Rf=sup⁡q>0R(qf)\mathcal R^f=\sup_{q>0}R(qf)Rf=supq>0​R(qf); f=ef=ef=e is uniform pricing. For q≥0q\ge0q≥0, the filtered demand G(q)G(q)G(q) is the fff-weighted demand at the personalized optima, keeping the pairs (i,t)(i,t)(i,t) with q≤pˉt,i/fiq\le\bar p_{t,i}/f_iq≤pˉ​t,i​/fi​, and H(q)=∫∑i∈Nfidt,i(qf) dμ(t)H(q)=\int\sum_{i\in N}f_id_{t,i}(qf)\,d\mu(t)H(q)=∫∑i∈N​fi​dt,i​(qf)dμ(t). Assumption A1 is G(q)≤H(q)G(q)\le H(q)G(q)≤H(q) for all q≥0q\ge0q≥0. Finally, qmin⁡q_{\min}qmin​ and qmax⁡q_{\max}qmax​ are the infimum and supremum of the ratios pˉt,i/fi\bar p_{t,i}/f_ipˉ​t,i​/fi​.

Formalization targets

The bound

Rˉ≤(1+ln⁡qmax⁡qmin⁡)Rf\bar{\mathcal R}\le\Bigl(1+\ln\frac{q_{\max}}{q_{\min}}\Bigr)\mathcal R^fRˉ≤(1+lnqmin​qmax​​)Rf

for every finite nonempty product set, every probability space of types, every positive fff and every selection pˉ\bar ppˉ​ of optimal prices, under A0 and A1, with qmin⁡>0q_{\min}>0qmin​>0 and integrable demands.

Exact tightness

For every ρ>1\rho>1ρ>1 there is an instance with one product, f=1f=1f=1, qmin⁡=1q_{\min}=1qmin​=1, qmax⁡=ρq_{\max}=\rhoqmax​=ρ, Rf>0\mathcal R^f>0Rf>0 and

Rˉ=(1+ln⁡ρ) Rf.\bar{\mathcal R}=(1+\ln\rho)\,\mathcal R^f .Rˉ=(1+lnρ)Rf.

Best constant over finite instances

For every ρ>1\rho>1ρ>1 and c<1+ln⁡ρc<1+\ln\rhoc<1+lnρ there is an instance with one product, f=1f=1f=1, finitely many types, qmin⁡=1q_{\min}=1qmin​=1, qmax⁡=ρq_{\max}=\rhoqmax​=ρ, Rf>0\mathcal R^f>0Rf>0 and Rˉ>c Rf\bar{\mathcal R}>c\,\mathcal R^fRˉ>cRf.

All three are clauses of the goal theorem, SingleFactorPricing.LogBound.theorem_1. The milestones are the displayed steps of the paper's proof (the integral forms of the profits, the split at qmin⁡q_{\min}qmin​, the identity R(qf)=qH(q)R(qf)=qH(q)R(qf)=qH(q), the pointwise bound on qG(q)qG(q)qG(q)) and its tight instance.

What the result gives

Since Rf≤R∗≤Rˉ\mathcal R^f\le\mathcal R^*\le\bar{\mathcal R}Rf≤R∗≤Rˉ, where R∗\mathcal R^*R∗ is the optimal non-personalized profit, one inequality yields a guarantee for single-factor heuristics against non-personalized pricing and a bound on the value of personalization. The constant depends only on the spread of the optimal prices measured against fff, not on the number of products or types. For example, pˉij∈[h,7.38h]\bar p_{ij}\in[h,7.38h]pˉ​ij​∈[h,7.38h] gives Rˉ≤3Re\bar{\mathcal R}\le3\mathcal R^eRˉ≤3Re. The paper also gives sufficient conditions for A1 that hold in linear, multinomial logit and envy-free models, and uses the bound to compare linear with non-linear pricing.

The result is proved in the paper. This mission adds a machine-checked version: the bound over a general type space, where the paper only remarks that the proof carries over; both meanings of tightness, of which the paper states neither separately; and an explicit probability law for the instance. No formalization of this bound for general demand is known. The platform has narrower statements of the same shape: RevenueOrdered.UDPmin.uniform_pricing_log_ratio_bound for envy-free pricing and ValuePersPricing.Deviation.theorem_2 for one product.

Where the difficulty lies

The obvious comparison fails. Personalized optimal prices need not lie on a single ray {qf}\{qf\}{qf}, so no single qqq can be compared with them, and A1 compares two sums at the same qqq rather than two profits. The bound therefore has no term-by-term proof, and the factor β\betaβ is not visible from any one type.

Three points need care in Lean. Over a general type space, Rˉ\bar{\mathcal R}Rˉ involves an exchange of integrals for functions that are only assumed integrable in the type. The second line of the printed chain contains ∫R(qf) dq/q\int R(qf)\,dq/q∫R(qf)dq/q, and nothing makes q↦R(qf)q\mapsto R(qf)q↦R(qf) measurable. And "tight" cannot mean equality with finitely many types: when qmin⁡<qmax⁡q_{\min}<q_{\max}qmin​<qmax​ and Rf>0\mathcal R^f>0Rf>0, GGG is a step function and equality is impossible, so the exact statement and the finite statement are independent, and the instance for each needs its own computation.

Formalization scope

The type space is T with a probability measure μ, and demand is d : T → (N → ℝ) → N → ℝ with N a nonempty Fintype. Non-negativity and optimality are required for p≥0p\ge0p≥0 only. KallMayerSLP.ICC.discreteLaw θ id, a published finite law, recovers the paper's finite model, and regularity is then automatic. qmin⁡q_{\min}qmin​, qmax⁡q_{\max}qmax​ and Rf\mathcal R^fRf are fixed by IsGLB, IsLUB and IsLUB, so a real ⨆ is not used; Rf\mathcal R^fRf need not be attained. The filter uses the paper's non-strict inequality. The optimality of pˉ\bar ppˉ​ stays in AssumptionA0 although the bound does not use it, so that Rˉ\bar{\mathcal R}Rˉ is the personalized optimum and not the profit of arbitrary prices. There is no cost vector (footnote 1 of the paper).

Trivial tightness is excluded: the clauses assert Rf>0\mathcal R^f>0Rf>0, and LogBoundInstance fixes qmin⁡=1q_{\min}=1qmin​=1 and qmax⁡=ρq_{\max}=\rhoqmax​=ρ as the true infimum and supremum of the ratios, not bounds on them.

A complete development needs integrals over the product of the type space and (0,∞)(0,\infty)(0,∞), the logarithmic integral, the moments and tail of the truncated law tailLaw, and an explicit family of finite instances for the last clause. The definitions are general enough to reuse for the other results of the paper. Proofs of the milestones and of the goal are welcome, as are the omitted results (Corollary 1, Proposition 1, Theorem 2), which would be separate statements.

Selected references

  • G. Gallego, G. Berbeglia, Bounds and Heuristics for Multi-Product Personalized Pricing, arXiv preprint, 2021. https://arxiv.org/abs/2102.03038
  • G. Berbeglia, G. Joret, Assortment optimisation under a general discrete choice model: A tight analysis of revenue-ordered assortments, Algorithmica 82(4), 2020. https://arxiv.org/abs/1606.01371
  • D. Bergemann, F. Castro, G. Y. Weintraub, Uniform pricing versus third-degree price discrimination, Cowles Foundation Discussion Paper, 2020. https://cowles.yale.edu/sites/default/files/2022-08/d2213-r.pdf
  • A. N. Elmachtoub, V. Gupta, M. Hamilton, The value of personalized pricing, Management Science, 2021. https://doi.org/10.1287/mnsc.2020.3821
  • D. A. Malueg, C. M. Snyder, Bounding the relative profitability of price discrimination, International Journal of Industrial Organization 24(5), 2006. https://doi.org/10.1016/j.ijindorg.2005.11.006
24 thms1 active userReviewed
🏆Completed
Algorithmic Game TheoryEconomicsOperations Research+1·Captain: mikeh

Unit-Demand Pricing for Min-Buying CustomersResearch Paper

Motivation

A seller offering several substitutable products must decide whether to charge one common price or a different price for each product. Customers may compare every product's net utility, or they may buy the cheapest product that they can afford. The second behavior is min-buying: a customer checks whether each product's valuation covers its price and selects the cheapest acceptable product. The distinction changes how inexpensive products divert purchases from expensive ones. Cui, Hamilton, Lamperski, and Ganjizadeh Tari study the resulting value of uniform pricing in symmetric and asymmetric valuation environments. The source for this mission is the September 27, 2026 manuscript, including its electronic companion.

Setting

There are n≥1n\ge1n≥1 products. A customer has a nonnegative random valuation vector V=(V1,…,Vn)V=(V_1,\ldots,V_n)V=(V1​,…,Vn​). A feasible price vector p=(p1,…,pn)p=(p_1,\ldots,p_n)p=(p1​,…,pn​) has nonnegative coordinates. The customer's payment is the smallest pip_ipi​ among products satisfying Vi≥piV_i\ge p_iVi​≥pi​, or zero if there is no acceptable product. Equal prices create no payment ambiguity. This defines expected differentiated-pricing revenue RDP(p)R_{\mathrm{DP}}(p)RDP​(p) and uniform-pricing revenue RUP(q)=RDP(q,…,q)R_{\mathrm{UP}}(q)=R_{\mathrm{DP}}(q,\ldots,q)RUP​(q)=RDP​(q,…,q).

The optimal revenues are RDPR_{\mathrm{DP}}RDP​ and RUPR_{\mathrm{UP}}RUP​. The source explicitly assumes that these optima are finite, attained, and that RUP>0R_{\mathrm{UP}}>0RUP​>0. The formal model records those conditions. It does not impose finite expected valuations: fixed finite prices have bounded payments even for valuation laws with infinite first moments. General joint probability measures describe dependent valuations; product measures describe independent valuations.

A law is exchangeable when permuting its product coordinates leaves it unchanged. It is affiliated when it has a joint density or probability mass function fff satisfying f(v)f(w)≤f(v∨w)f(v∧w)f(v)f(w)\le f(v\vee w)f(v\wedge w)f(v)f(w)≤f(v∨w)f(v∧w) for every pair of valuation vectors. A nonnegative marginal law is regular when its revenue-quantile curve is concave. Regularity uses interior quantiles, with endpoint conventions stated explicitly; it is not replaced by the stronger monotone-hazard-rate condition.

Formalization targets

Goal: Theorem 4, the exact tight regular-distribution guarantee

For independent, absolutely continuous regular valuation distributions, the target is

RUP≥x⋆RDP.R_{\mathrm{UP}}\ge x^\star R_{\mathrm{DP}}.RUP​≥x⋆RDP​.

The constant is defined exactly. Let A(0)=0A(0)=0A(0)=0 and A(x)=−log⁡(1−x2)/xA(x)=-\log(1-x^2)/xA(x)=−log(1−x2)/x for 0<x<10<x<10<x<1, and let ccc solve

c′(x)=A′(x)(1−xc(x)),c(0)=1.c'(x)=A'(x)(1-xc(x)),\qquad c(0)=1.c′(x)=A′(x)(1−xc(x)),c(0)=1.

Then x⋆x^\starx⋆ is the unique point in (0,1)(0,1)(0,1) with x⋆c(x⋆)=1x^\star c(x^\star)=1x⋆c(x⋆)=1. The main target includes this nondegenerate root characterization, the universal guarantee, and sharpness: for every ε>0\varepsilon>0ε>0 there is a finite independent continuous regular instance with RDP/RUP>1/x⋆−εR_{\mathrm{DP}}/R_{\mathrm{UP}}>1/x^\star-\varepsilonRDP​/RUP​>1/x⋆−ε. The number of products can vary with ε\varepsilonε. The manuscript's decimal x⋆≈0.6986x^\star\approx0.6986x⋆≈0.6986 is explanatory; no decimal is substituted for the exact constant. These are the two conclusions of Theorem 4, p.17.

Symmetric-market results

Theorems 1 and 2, pp.10 and 13, assert RDP=RUPR_{\mathrm{DP}}=R_{\mathrm{UP}}RDP​=RUP​ for IID valuations and for exchangeable affiliated valuations, respectively. Supporting targets state the stronger policy comparison proved in the manuscript: every price vector is dominated by a uniform price selected from that vector. Theorem 3, p.15, gives an exchangeable nonnegative valuation law for every positive nnn such that

RDP≥332log⁡(n)RUP.R_{\mathrm{DP}}\ge\frac{3}{32}\log(n)R_{\mathrm{UP}}.RDP​≥323​log(n)RUP​.

This last construction is discrete and does not carry a density assumption.

Supporting targets

The milestone list includes the independent-revenue identities, the affiliated conditional-threshold comparison, Lemmas 1–4 on triangular and canonical distributions, Lemmas EC.1–EC.2 on rejection probabilities and continuous approximation, and the analytic properties used to specify the exact constant. Additional interface targets identify the concrete triangular probability laws and connect the canonical continuation formulas to actual min-buying revenues. All main results and numbered supporting lemmas of the source are represented in this one mission.

Significance

The symmetric results specify settings where charging different prices gives no revenue advantage. The asymmetric result gives an exact worst-case guarantee and establishes that a better universal fraction cannot hold throughout the stated regular class. The exchangeable lower bound distinguishes symmetry alone from symmetry combined with affiliation. These interpretations follow from the paper's statements and comparisons in Sections 3–4.

The formalization makes the customer payment rule, probability classes, optimization assumptions, and sharpness quantifiers explicit. It also separates an ODE-defined constant from its numerical approximation. The current deliverable consists of locally compiled definitions and open theorem statements, with independent blind read-backs. Statement compilation is not a proof of a theorem; completed Lean proofs are the mission's remaining work.

Mathematical difficulties

Atomic cap distributions require inclusive purchasing and strict rejection probabilities. Substituting a right-continuous CDF for strict rejection at a cap would change the revenue. Affiliation must be used in its density-or-mass-function meaning; the conditional-tail comparison is a theorem to establish. Discrete order statistics can contain ties and need the corresponding multiplicity treatment. The tightness claim must ultimately hold for continuous regular distributions, although several intermediate extremal distributions have atoms. The continuous approximation lemma is therefore part of the target list.

The exact constant requires existence, uniqueness, and maximum properties of its concrete integral benchmark. Those properties are separate conclusions, not hypotheses of the market model. The canonical continuation value is revenue from the designated cap prices; it is not defined to equal the optimal differentiated revenue.

Formalization scope

Valuation and price vectors are real functions on a finite index type. Joint laws are probability measures; no finite support or moment restriction is inserted into the general pricing theorems. The original continuous-market conventions and positive attained optimum conditions are explicit where used. Atomic triangular laws and the exchangeable lower-bound witness have their appropriate separate scopes. Prices are optimized over all nonnegative real vectors, with no fixed ordering imposed on the feasible set; sorted order is required only in formulas that use it.

Regularity is concavity of qF−1(1−q)qF^{-1}(1-q)qF−1(1−q) for 0<q<10<q<10<q<1, avoiding undefined endpoint quantiles. For a nonnegative curve this is equivalent to concavity of its zero-endpoint extension. Affiliation covers both a Lebesgue density and a probability mass function on real-valued vectors. The ODE benchmark is an explicit integrating-factor integral; the main goal proves that its root is unique and lies strictly between zero and one. No desired revenue comparison or extremal bound is assumed inside a definition.

Reusable components include the min-buying payment model, strict-rejection identities, quantile regularity, exchangeability and affiliation, triangular laws, and finite-to-continuous comparison. Proofs of these components and of the paper's complete targets are welcome. All code targets Lean 4.33.1 and the pinned Mathlib environment.

Selected references

  • Titing Cui, Michael L. Hamilton, Jourdain Lamperski, and Kasra Ganjizadeh Tari, Unit-Demand Pricing for Min-Buying Customers, September 27, 2026 manuscript. SSRN / DOI.
26 thms1 active userReviewed
🏆Completed
EconomicsOperations ResearchOptimization·Captain: mikeh

Optimal Feature-Based Market Segmentation and PricingResearch Paper

Motivation

A seller can often predict a customer's willingness to pay from observed features, yet may be limited to a small number of prices. The decisions are then coupled: the seller must choose both the customer groups and the price assigned to each group. Grouping customers first and optimizing prices afterward does not generally solve the joint problem. Cui and Hamilton study this joint optimization in Optimal Feature-Based Market Segmentation and Pricing, using the structure of a regression model with a common prediction error.

The central question is when that statistical structure makes optimal segmentation tractable. The main positive result is Theorem 2. The companion results describe hardness without the additional shape condition, approximation by greedy price-menu selection, and the value of allowing more segments. The source for the page-specific formalization is the October 22, 2025, 60-page manuscript available under SSRN 4151103.

Setting

Let M=μ(X)≥0M=\mu(X)\ge0M=μ(X)≥0 denote a predicted valuation and let EEE denote a mean-zero error independent of the observed features. Realized valuation is V=M+EV=M+EV=M+E. The general model uses separate probability laws ν\nuν for MMM and η\etaη for EEE, with finite first moments. Independence is represented by the shared location family: every customer with predicted value xxx has valuation law x+Ex+Ex+E.

The inclusive survival function is S(t)=Pr⁡(E≥t)S(t)=\Pr(E\ge t)S(t)=Pr(E≥t). At nonnegative price ppp, expected customer revenue is r(x,p)=pS(p−x)r(x,p)=pS(p-x)r(x,p)=pS(p−x). Inclusive survival is essential when the error distribution has atoms. A finite price menu PPP earns

F(P)=∫max⁡({0}∪{r(x,p):p∈P}) dν(x).F(P)=\int \max\bigl(\{0\}\cup\{r(x,p):p\in P\}\bigr)\,d\nu(x).F(P)=∫max({0}∪{r(x,p):p∈P})dν(x).

For a fixed menu, assigning every customer to a revenue-maximizing menu price optimizes its segmentation. Write Rk=sup⁡∣P∣≤kF(P)R_k=\sup_{|P|\le k}F(P)Rk​=sup∣P∣≤k​F(P), R∞=∫sup⁡p≥0r(x,p) dν(x)R_\infty=\int\sup_{p\ge0}r(x,p)\,d\nu(x)R∞​=∫supp≥0​r(x,p)dν(x), and Rk0R_k^0Rk0​ for the same finite-menu problem with zero error. The MHR condition is encoded for density-bearing error laws by log-concavity of survival on its positive domain, the formulation used in the proof of Lemma 2(b). This permits bounded support and zero purchase probabilities.

Formalization targets

Goal: Theorem 2, optimal computation

For nnn distinct, sorted predicted values with positive probability masses, independent mean-zero MHR error, and a primitive oracle that returns an optimal block price and its attained revenue, construct a globally optimal policy with at most kkk contiguous index blocks. Its revenue must equal the specified Bellman recurrence, and the reconstructed policy must dominate every assignment of types to kkk arbitrary nonnegative prices. If each nonempty block-oracle call costs at most mmm, the specified memoized schedule has

T(n,k,m)=O(n2(k+m)).T(n,k,m)=O\bigl(n^2(k+m)\bigr).T(n,k,m)=O(n2(k+m)).

The draft states an explicit upper bound 5n2(k+m)5n^2(k+m)5n2(k+m) for oracle charges, Bellman candidate evaluations, and backtracking. This is an oracle-operation model for already sorted data; it is not a bit-complexity assertion about evaluating arbitrary real probabilities.

Companion results

Theorem 1 asserts NP-hardness even with finitely supported independent mean-zero noise of polynomial support size. Its target uses an explicit binary encoding of rational input data and polynomial-time finite Turing-machine reductions; kkk is part of the input. Remark 2 targets normalized, monotone, submodular menu revenue and the greedy guarantee

F(Gk)≥(1−e−1)max⁡A⊆P, ∣A∣≤kF(A).F(G_k)\ge(1-e^{-1})\max_{A\subseteq P,\ |A|\le k}F(A).F(Gk​)≥(1−e−1)A⊆P, ∣A∣≤kmax​F(A).

A separate finite-support target shows that nonnegative realized valuations, together with zero, suffice as candidate prices. Theorem 4 states

R∞−Rk≤U−LkR_\infty-R_k\le\frac{U-L}{k}R∞​−Rk​≤kU−L​

for any finite almost-sure prediction bounds L≤M≤UL\le M\le UL≤M≤U. Theorem 5 states diminishing increments,

Rk+2−Rk+1≤Rk+1−Rk,k≥1,R_{k+2}-R_{k+1}\le R_{k+1}-R_k,\qquad k\ge1,Rk+2​−Rk+1​≤Rk+1​−Rk​,k≥1,

under MHR error. These two targets retain general prediction distributions. Theorem 3 includes almost-sure nonnegative realized valuations. Corollary 2 uses this same setting and includes the factor two confirmed by the author during review.

Significance

Theorem 2 replaces a search over arbitrary segmentations with an optimal interval policy and a quadratic-size block table. The performance results explain how a small menu approaches feature-based personalized pricing and when the incremental value of another segment diminishes. The hardness and approximation targets delineate what remains possible when the MHR restriction is removed.

The draft makes ties, atoms, support endpoints, optimizer attainment, and computational costs explicit. Generic proofs of monotone finite-menu assignment under a TP2 revenue condition and Bellman optimality have been checked locally in the pinned Lean environment. The full paper theorems remain open formal proof obligations; well-formed statements are not completed proofs.

Mathematical difficulties and author-confirmed corrections

The interval argument needs an attaining assignment that is monotone despite ties and all-zero revenue rows. Division by survival probabilities would exclude legitimate bounded-support cases. The complexity claim also needs the block oracle's attained revenue, not merely its price. The hardness target requires an actual polynomial-time reduction and careful centering of the constructed error while preserving nonnegative predicted means.

The author confirmed two corrections on October 30, 2026. Theorem 3 now explicitly assumes V=M+E≥0V=M+E\ge0V=M+E≥0 almost surely. Corollary 2 now includes the factor two:

R∞−Rk≤22σ(E[M]+σ)k.R_\infty-R_k\le2\sqrt{\frac{2\sigma(\mathbb E[M]+\sigma)}{k}}.R∞​−Rk​≤2k2σ(E[M]+σ)​​.

Corollary 2 retains the same nonnegative realized-valuation setting used in the corrected Theorem 3. Both results are definitive milestone targets in this draft. Two earlier draft records are retained as clearly labeled superseded copies, with their mathematical statements corrected as well; they are not milestones. The nonnegative-valuation assumption supports the proof's bound Rϵ(x)≤xR_\epsilon(x)\le xRϵ​(x)≤x. For Corollary 2, minimizing 2jσ/k+(E[M]+σ)/j2j\sigma/k+(\mathbb E[M]+\sigma)/j2jσ/k+(E[M]+σ)/j gives the corrected factor two. Formal proofs must still handle conditional truncation, boundary atoms, and zero variance.

Formalization scope

General revenue targets use probability measures on the real line, nonnegative integrable predictions, and integrable mean-zero error. The finite computational model uses strictly increasing predicted values and positive weights summing to one. At most kkk price slots are allowed; repeated cuts represent unused slots. Index intervals include every type, including the largest support point. Least optimal prices implement a deterministic tie convention; no interior first-order condition or uniqueness assumption is imposed.

The definition layer contains no proof placeholders or new axioms. Real suprema and infima are defined from the feasible economic objects; their attainment and boundedness must be proved where needed. The MHR-to-TP2 bridge is a theorem target, not an assumption that deletes the analytic content. The Bellman recurrence, its partition relation, and its charged evaluation schedule are separate objects. The goal connects that recurrence to the original unrestricted assignment problem and an actual feasible reconstruction.

The general model, the rational binary-input hardness model, and the discrete candidate-price model are explicit interfaces. Their connecting lemmas, optimizer existence, the analytic MHR argument, interval-to-DP reconstruction, and the source's companion results remain available proof obligations. The optional derivative refinement of Lemma 1 and appendix extension/necessity theorems are outside this first proposal. Every uploaded item includes an independent read-back for comparison with its intended mathematical content.

Selected references

  • Titing Cui and Michael L. Hamilton, Optimal Feature-Based Market Segmentation and Pricing, Operations Research 74(3), 1354–1370 (2026). DOI. Page citations refer to the October 22, 2025 manuscript, SSRN 4151103.
  • Adrien Saumard and Jon A. Wellner, Log-concavity and strong log-concavity: a review, Statistics Surveys 8, 45–114 (2014), Section 11, cited in the source's proof of Lemma 2(b). DOI.
24 thms1 active userReviewed

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