基于情境与可分需求模型的动态定价:
数理建模、在线学习算法与统计复杂度全景深度解析

Context-Based Dynamic Pricing with Separable Demand Models
Authors: Jinzhi Bu, David Simchi-Levi, Chonghuan Wang · Management Science (2026)

【文献出版与分类档案元数据】

【英文原题】Context-Based Dynamic Pricing with Separable Demand Models

【论文作者】Jinzhi Bu (The Hong Kong Polytechnic University), David Simchi-Levi (Massachusetts Institute of Technology), Chonghuan Wang (The University of Texas at Dallas)

【发表期刊】Management Science, Vol. 72, No. 8 (August 2026), pp. 7279–7298

【数字对象标识符 (DOI)】https://doi.org/10.1287/mnsc.2022.02260

【JEL 分类号】C44 (Operations Research; Statistical Decision Theory), C14 (Semiparametric and Nonparametric Methods), D42 (Monopoly Pricing), D83 (Search; Learning; Information and Knowledge), L81 (Retail and Wholesale Trades; e-Commerce), M31 (Marketing: Pricing and Revenue Management)

【核心关键词】可分模型 (Separable Model)、动态定价 (Dynamic Pricing)、情境信息 (Contextual Information)、在线学习 (Online Learning)、后悔界分析 (Regret Bounds)、统计复杂度 (Statistical Complexity)

一、研究背景、产业实践动机与核心学术贡献

1.1 现实商业背景:中东头部消费电子零售商的数据特征与长尾挑战

产业实践痛点:在现代电子商务与零售平台中,企业拥有海量的外生环境情境数据(Contextual Information,即协变量 x∈ℝdx \in {\mathbb{R}}^{d}),如商品物理属性、时间季节性、促销节日、宏观经济指数以及竞争对手的实时价格。企业面临的核心挑战是如何利用这些高维情境信息进行动态定价,以最大化全时域的期望累积收益。然而,真实商业场景中存在严重的长尾效应与低销量困境。

本文的研究直接源于作者团队与中东地区一家领先的大型消费电子零售商的深度产业合作(Li et al. 2023)。该企业每两天为一个销售周期,针对手机、笔记本电脑等数码产品依据历史销量与外生特征动态调整价格。研究团队对该企业 2000 个手机 SKU(Stock Keeping Units)跨越三年的真实交易数据进行了系统分析,发现高达 83% 的 SKU 在每个销售周期内的平均销量不足 5 件!这种极度普遍的低销量长尾特征使得单个 SKU 能够提供的历史需求样本极其稀疏,导致企业难以从单一产品的有限数据中准确辨识真实的价格弹性(Price Elasticity)。

Figure 1: CatBoost 非参数树模型与可分模型在真实手机 SKU 上的价格需求弹性曲线对比

1.2 传统机器学习与非参数树模型的失效:价格弹性平坦化困境

传统黑箱模型的缺陷:为了克服单一产品样本稀疏问题,工业界的标准预处理做法是采用 K-Means 聚类技术,将具有相似历史属性(如折扣幅度、竞品价格、发布时间、颜色、存储容量、品牌及平均销量)的 SKU 划分为同一产品簇(Cluster),并假设同一簇内的产品共享相同的需求预测模型。然而,当采用业界最流行的非参数回归树算法(如 Random Forest、LightGBM、CatBoost 等)进行需求拟合时,却出现了严重的“价格弹性失真”:

从上图 Figure 1 可以清晰看出,非参数树模型(CatBoost)拟合出的价格-需求曲线几乎是一条平坦的水平线,计算出的价格弹性趋近于 0。这是因为传统的非参数回归树将商家的内生定价决策 pp 与外生环境特征 xx 混为一谈,视为完全同质的输入特征,从而忽视了内生决策与外生环境在经济学因果机制上的根本差异。这种平坦的弹性曲线使得算法误以为提高价格不会损失销量,进而产生灾难性的定价决策。

相比之下,作者提出的可分需求模型 f(p)+g(x)f(p) + g(x)(采用分段线性函数逼近价格效应 f(p)f(p),采用 CatBoost 估计情境效应 g(x)g(x)),能够恢复出清晰单调递减的价格需求曲线,完美契合经济学常识。下表 Table 1 与下图 Figure 2 进一步展示了在包含 678 个 SKU 的代表性聚类簇上的预测精度对比:可分模型在平均均方误差、均方误差标准差以及时序平均决定系数 r2r^{2} 上均显著优于传统线性模型和纯 CatBoost 模型。

评估指标 (Evaluation Metric) 可分需求模型 (f(p)+g(x)f(p) + g(x)) CatBoost 非参数树模型 传统线性模型 (bp+a⊤xbp + a^{\top}x)
均方误差均值 (Average Squared Errors) 24.520 25.232 31.699
均方误差标准差 (Std. Dev. of Squared Errors) 225.192 238.185 269.583
跨周期平均决定系数 (Average of r2r^{2} Scores) 0.644 0.642 0.544

Figure 2: 三种预测模型(可分模型、CatBoost、线性模型)在测试周期上的预测平均销量与真实销量走势对比

1.3 可分需求模型 f(p) + g(x) 的理论定位与学术优势

兼具解释性与鲁棒性的中道理论:在动态定价理论中,现有文献主要存在两大建模流派:一是参数化模型(Parametric Models,如经典的线性模型 bp+a⊤xbp + a^{\top}x),虽然估计简便且易于设计在线算法,但极易遭遇模型误设(Model Misspecification)风险;二是完全非参数模型(Fully Nonparametric Models,如 d(p,x)d(p,x)),虽然理论稳健性强,但在高维情境下饱受“维度诅咒(Curse of Dimensionality)”的困扰,收敛速度极其缓慢(如 Slivkins 2011 的后悔界为 Θ∼(T(d+2)/(d+3))\overset{\sim}{\Theta}(T^{(d + 2)/(d + 3)}),Chen and Gallego 2021 为 Θ∼(T(d+2)/(d+4))\overset{\sim}{\Theta}(T^{(d + 2)/(d + 4)}))。

可分需求模型 Dt(p)=f(p)+g(x)+εtD_{t}(p) = f(p) + g(x) + \varepsilon_{t} 完美继承了参数模型的结构优势与非参数模型的鲁棒性:它不仅保持了高度的经济学可解释性,有效避免了参数模型的误设风险,更重要的是通过结构可分性,极大地降低了在线学习的统计复杂度(Statistical Complexity)。

1.4 本文核心研究问题与主要学术贡献梳理

核心科学问题:虽然 Li et al. (2023) 在离线(Offline)设定下验证了可分模型的经验效果,但在必须一边定价探索一边收集数据的在线学习(Online Learning)设定下,学术界对可分需求模型的统计学习复杂度尚缺乏理论刻画。本文的核心研究问题是:在可分需求模型下,在线动态定价问题的统计复杂度究竟由何决定?内生价格效应 f(p)f(p) 与外生情境效应 g(x)g(x) 的不同函数结构如何从本质上改变最优学习速率?

【本文核心学术贡献提要】

1. 【系统构建理论框架与三大配置】:首次系统研究了三种具有代表性的可分需求模型配置:SMLPE(线性价格效应 + 高阶光滑非参数情境效应)、SMLCE(高阶光滑非参数价格效应 + 线性情境效应)、SMNPE(双重高阶光滑非参数效应)。

2. 【创新在线学习算法设计】:针对聚合需求中 f(p)f(p) 与 g(x)g(x) 贡献无法单独观测且估计误差相互耦合的核心挑战,针对性地提出了三种高效在线学习算法:ASMLPE(基于随机价格冲击与桥接优化误差控制)、ASMLCE(基于局部多项式近似与有偏线性 Bandit 的“乐观之上二次乐观”OFU 原则)、ASMNPE(基于网格离散化与探索利用严格分离)。

3. 【证明紧致的后悔极小化极大上下界】:推导了所有算法的严谨后悔上界,并通过创新性地构造自适应于非整数阶光滑参数 kg,kfk_{g},k_{f} 与光滑度常数 L,δL,\delta 的困难实例族,证明了匹配的极小化极大下界,确立了算法在核心参数区间的速率最优性(Rate-optimality)。

4. 【揭示内生与外生变量的统计非对称性】:深刻揭示了在在线学习中,内生决策变量 pp(算法可主动控制并施加实验冲击)与外生随机变量 xx(被动由自然界独立生成)在决定学习复杂度上的本质非对称性。

文献来源 需求模型结构 后悔上界 (Upper Bound) 后悔下界 (Lower Bound) 核心假设与条件
Keskin & Zeevi (2014) bp+abp + a O∼(T)\overset{\sim}{O}(\sqrt{T}) Ω(T)\Omega(\sqrt{T}) 无情境,纯线性需求
Wang et al. (2021a) f(p)f(p) O∼(Tk+12k+1)\overset{\sim}{O}(T^{\frac{k + 1}{2k + 1}}) Ω(Tk+12k+1)\Omega(T^{\frac{k + 1}{2k + 1}}) 无情境,kk 阶光滑需求
Nambiar et al. (2019) bp+g(x)bp + g(x) O(T)O(\sqrt{T}) Ω(T)\Omega(\sqrt{T}) 基准为最优线性近似
Slivkins (2011) f(p,x)f(p,x) O(Td+2d+3)O(T^{\frac{d + 2}{d + 3}}) Ω(Td+2d+3)\Omega(T^{\frac{d + 2}{d + 3}}) 完全非参数 Lipschitz
Chen & Gallego (2021) f(p,x)f(p,x) O∼(Td+2d+4)\overset{\sim}{O}(T^{\frac{d + 2}{d + 4}}) Ω(Td+2d+4)\Omega(T^{\frac{d + 2}{d + 4}}) 非参数 Lipschitz + 局部凹
本文 (Section 2: SMLPE) bp+g(x)bp + g(x) O∼(T∨(L2T)dd+2kg)\overset{\sim}{O}\left( \sqrt{T} \vee (L^{2}T)^{\frac{d}{d + 2k_{g}}} \right) Ω(T∨(L2T)dd+2kg)\Omega\left( \sqrt{T} \vee (L^{2}T)^{\frac{d}{d + 2k_{g}}} \right) gg 为 kgk_{g} 阶光滑 (0<kg≤10 < k_{g} \leq 1)
本文 (Section 3: SMLCE) f(p)+a⊤xf(p) + a^{\top}x O∼(d(T∨(δTkf+1)12kf+1))\overset{\sim}{O}\left( d\left( \sqrt{T} \vee (\delta T^{k_{f} + 1})^{\frac{1}{2k_{f} + 1}} \right) \right) Ω(T∨(δTkf+1)12kf+1)\Omega\left( \sqrt{T} \vee (\delta T^{k_{f} + 1})^{\frac{1}{2k_{f} + 1}} \right) ff 为 kfk_{f} 阶光滑 (kf>0k_{f} > 0)
本文 (Section 4: SMNPE) f(p)+g(x)f(p) + g(x) O∼(T3/4∨(LdTd+kg)1d+2kg)\overset{\sim}{O}\left( T^{3/4} \vee (L^{d}T^{d + k_{g}})^{\frac{1}{d + 2k_{g}}} \right) Ω((LdTd+kg)1d+2kg)\Omega\left( (L^{d}T^{d + k_{g}})^{\frac{1}{d + 2k_{g}}} \right) 双重非参数光滑 (0<kg≤10 < k_{g} \leq 1)

二、动态定价与可分需求数理模型构建

2.1 决策环境、情境生成与需求生成机制

系统运行协议:考虑一个垄断商家在有限离散时域 TT 期内销售某种单品。决策时序与信息结构如下:

1. 【情境揭示】:在每一期 t∈{1,2,…,T}t \in \{ 1,2,\ldots,T\} 开始时,商家首先观察到当期的外生情境向量 xt∈[0,1]d⊂ℝdx_{t} \in \lbrack 0,1\rbrack^{d} \subset {\mathbb{R}}^{d}。假设情境序列 {xt}t≥1\{ x_{t}\}_{t \geq 1} 是从定义在 [0,1]d\lbrack 0,1\rbrack^{d} 上的未知联合概率分布 𝒫\mathcal{P} 中独立同分布(i.i.d.)抽取的随机向量。

2. 【定价决策】:商家基于当期情境 xtx_{t} 以及截止到 t−1t - 1 期的全部历史观测数据,从闭区间可行价格集合 [p_,p¯]⊂(0,∞)\lbrack\underline{p},\overline{p}\rbrack \subset (0,\infty) 中选定并展示价格 ptp_{t}。

3. 【需求实现】:市场顾客观察到价格 ptp_{t} 后产生随机需求 Dt(pt)D_{t}(p_{t})。需求生成过程遵循未知但具有可分结构的函数形式:

Dt(pt)=f(pt)+g(xt)+εtD_{t}(p_{t}) = f(p_{t}) + g(x_{t}) + \varepsilon_{t}

其中 f:[p_,p¯]→ℝf:\lbrack\underline{p},\overline{p}\rbrack \rightarrow {\mathbb{R}} 为未知的内生价格需求效应函数,g:[0,1]d→ℝg:\lbrack 0,1\rbrack^{d} \rightarrow {\mathbb{R}} 为未知的外生情境需求效应函数。随机扰动项 {εt}t≥1\{\varepsilon_{t}\}_{t \geq 1} 为相互独立同分布的零均值 σ2\sigma^{2}-亚高斯(σ2\sigma^{2}-sub-Gaussian)随机变量,即满足对任意 λ∈ℝ\lambda \in {\mathbb{R}},有 𝔼[exp(λεt)]≤exp(λ2σ2/2){\mathbb{E}}\lbrack\exp(\lambda\varepsilon_{t})\rbrack \leq \exp(\lambda^{2}\sigma^{2}/2)。

2.2 期望收益函数与全知最优定价基准

给定价格 pp 和情境 xtx_{t},商家的单期期望收益函数定义为:

r(p,xt)≜p⋅𝔼[Dt(p)|xt]=p(f(p)+g(xt))r(p,x_{t}) \triangleq p \cdot {\mathbb{E}}\lbrack D_{t}(p)|x_{t}\rbrack = p(f(p) + g(x_{t}))

如果商家预先完全知晓真实的函数形式 f(⋅)f( \cdot ) 和 g(⋅)g( \cdot ),则在面对情境 xtx_{t} 时,全知最优定价策略(Clairvoyant Optimal Policy)将选择使期望收益最大化的最优价格 pt*p_{t}^{*}:

pt*≜argmaxp∈[p_,p¯]r(p,xt)=argmaxp∈[p_,p¯]p(f(p)+g(xt))p_{t}^{*} \triangleq \arg\max_{p \in \lbrack\underline{p},\overline{p}\rbrack}r(p,x_{t}) = \arg\max_{p \in \lbrack\underline{p},\overline{p}\rbrack}p(f(p) + g(x_{t}))

2.3 容许定价策略与累积后悔(Regret)测度体系

一个容许定价策略(Admissible Pricing Policy)π\pi 由一系列决策映射序列 {πt}t≥1\{\pi_{t}\}_{t \geq 1} 构成。在第 tt 期,πt\pi_{t} 将截止到当期开始时的历史信息集 Ht=(x1,p1,D1,…,xt−1,pt−1,Dt−1,xt)H_{t} = (x_{1},p_{1},D_{1},\ldots,x_{t - 1},p_{t - 1},D_{t - 1},x_{t})(以及可能的外部随机化种子)映射为可行价格 pt∈[p_,p¯]p_{t} \in \lbrack\underline{p},\overline{p}\rbrack。策略 π\pi 在时域 TT 内的性能由累积期望后悔值(Cumulative Expected Regret)衡量:

Rπ(T)≜∑t=1T𝔼[r(pt*,xt)−r(pt,xt)]R^{\pi}(T) \triangleq \sum_{t = 1}^{T}{}{\mathbb{E}}\lbrack r(p_{t}^{*},x_{t}) - r(p_{t},x_{t})\rbrack

其中期望 𝔼[⋅]{\mathbb{E}}\lbrack \cdot \rbrack 是关于情境分布 𝒫\mathcal{P}、噪声分布以及策略 π\pi 可能采用的随机化机制的联合期望。算法设计的核心目标是寻找容许策略 π\pi,使得 Rπ(T)R^{\pi}(T) 随着时域 TT 的增长实现次线性(Sublinear)增长,即 limT→∞Rπ(T)/T=0\lim_{T \rightarrow \infty}R^{\pi}(T)/T = 0,并使后悔阶数尽可能达到信息论意义上的理论下界。

三、模型一 —— 线性价格效应可分模型 (SMLPE)

3.1 模型设定与高阶光滑函数类假设

在 SMLPE(Separable Model with Linear Pricing Effect)中,价格效应被假定为线性,而情境效应为高阶光滑非参数函数。需求模型具体为:

Dt(p)=bp+g(xt)+εt,∀p∈[p_,p¯]D_{t}(p) = bp + g(x_{t}) + \varepsilon_{t},\forall p \in \lbrack\underline{p},\overline{p}\rbrack

其中价格敏感度系数 bb 未知,但已知其落在先验区间 [b_,b¯]⊂(−∞,0)\lbrack\underline{b},\overline{b}\rbrack \subset ( - \infty,0) 内。由于 b<0b < 0,单期期望收益函数 r(p,xt)=p(bp+g(xt))r(p,x_{t}) = p(bp + g(x_{t})) 是关于价格 pp 的严格强凹(Strongly Concave)二次函数。全知最优无约束价格为 pt*=−g(xt)/(2b)p_{t}^{*} = - g(x_{t})/(2b),遵循文献通用设定(Ban and Keskin 2021),假设对所有 xt∈[0,1]dx_{t} \in \lbrack 0,1\rbrack^{d} 和 b∈[b_,b¯]b \in \lbrack\underline{b},\overline{b}\rbrack,均有 pt*∈[p_,p¯]p_{t}^{*} \in \lbrack\underline{p},\overline{p}\rbrack。

【Assumption 1:情境需求函数的高阶光滑性假设】

【定义:最大严格下界整数】:对任意实数 k>0k > 0,定义 ⌊k⌋*≜sup{i∈ℕ:i<k}\lfloor k\rfloor_{*} \triangleq \sup\{ i \in {\mathbb{N}}:i < k\} 为严格小于 kk 的最大非负整数。

【定义:多元偏导数与多重指标】:对多重指标 κ=(κ1,…,κd)∈ℕd\kappa = (\kappa_{1},\ldots,\kappa_{d}) \in {\mathbb{N}}^{d},记 |κ|≜∑i=1dκi|\kappa| \triangleq \sum_{i = 1}^{d}{}\kappa_{i},偏导算子为 ∂κg≜∂|κ|g∂x1κ1…∂xdκd\partial^{\kappa}g \triangleq \frac{\partial^{|\kappa|}g}{\partial x_{1}^{\kappa_{1}}\ldots\partial x_{d}^{\kappa_{d}}}。

【Assumption 1 (高阶光滑情境函数类)】:称函数 g:[0,1]d→ℝg:\lbrack 0,1\rbrack^{d} \rightarrow {\mathbb{R}} 属于 kgk_{g} 阶光滑函数类 𝒢d(kg,L)\mathcal{G}_{d}(k_{g},L)(光滑常数 L>0L > 0),若 g(⋅)g( \cdot ) 在 [0,1]d\lbrack 0,1\rbrack^{d} 上 ⌊kg⌋*\lfloor k_{g}\rfloor_{*} 次连续可微,且对任意满足 |κ|=⌊kg⌋*|\kappa| = \lfloor k_{g}\rfloor_{*} 的多重指标 κ\kappa 及任意 x1,x2∈[0,1]dx_{1},x_{2} \in \lbrack 0,1\rbrack^{d},均满足 Hölder 连续性:

|∂κg(x1)−∂κg(x2)|≤L∥x1−x2∥2kg−⌊kg⌋*|\partial^{\kappa}g(x_{1}) - \partial^{\kappa}g(x_{2})| \leq L\| x_{1} - x_{2}\|_{2}^{k_{g} - \lfloor k_{g}\rfloor_{*}}

【性质】:由该假设可直接导出,存在常数 g¯≥0\overline{g} \geq 0 使得 supx∈[0,1]d|g(x)|≤g¯\sup_{x \in \lbrack 0,1\rbrack^{d}}|g(x)| \leq \overline{g}。

3.2 ASMLPE 在线学习算法架构与三大创新机制

算法核心设计思想:针对 SMLPE,作者设计了算法 ASMLPE(Algorithm 1),巧妙融合了以下三大核心机制:

(1)随机价格冲击(Pricing with Random Shock)实现斜率 bb 的无偏解耦估计:传统方法若直接对 dtd_{t} 和 (pt,xt)(p_{t},x_{t}) 进行线性回归,由于内生价格 ptp_{t} 依赖于历史情境估计,会引发严重的价格内生性偏误(Price Endogeneity Bias)。ASMLPE 在贪婪投影价格 ptgp_{t}^{g} 的基础上叠加外生零均值随机冲击 Δt=±δt\Delta_{t} = \pm \delta_{t}(其中 δt=t−1/4\delta_{t} = t^{- 1/4},以各 1/2 概率取正负)。在第 tt 期末,通过仅对随机冲击 Δs\Delta_{s} 进行回归来估计斜率 bb:

∑s=1tΔsds∑s=1tΔs2=∑s=1tΔs(bps+g(xs)+εs)∑s=1tΔs2=b+∑s=1tΔs(bpsg+g(xs)+εs)∑s=1tΔs2\frac{\sum_{s = 1}^{t}\Delta_{s}d_{s}}{\sum_{s = 1}^{t}\Delta_{s}^{2}} = \frac{\sum_{s = 1}^{t}\Delta_{s}(bp_{s} + g(x_{s}) + \varepsilon_{s})}{\sum_{s = 1}^{t}\Delta_{s}^{2}} = b + \frac{\sum_{s = 1}^{t}\Delta_{s}(bp_{s}^{g} + g(x_{s}) + \varepsilon_{s})}{\sum_{s = 1}^{t}\Delta_{s}^{2}}

由于 Δs\Delta_{s} 独立于 psg,g(xs)p_{s}^{g},g(x_{s}) 和 εs\varepsilon_{s},上式第二项的期望精确为 0!这使得对价格敏感度 bb 的估计彻底解耦于对复杂未知函数 g(x)g(x) 的估计,且能够无障碍利用全时域所有历史样本(跨越所有情境区域),实现了极高的数据利用效率。

(2)情境空间网格化(Context Space Binning)与局部泰勒多项式逼近:将 dd 维情境超立方体 [0,1]d\lbrack 0,1\rbrack^{d} 沿每个维度等分为 MM 份,构成 MdM^{d} 个小立方体网格 MjM_{j}。在每个网格 MjM_{j} 内,选取中心点 mjm_{j},将 g(x)g(x) 用其 ⌊kg⌋*\lfloor k_{g}\rfloor_{*} 阶多元泰勒多项式 PMj(x)=∑i=0⌊kg⌋*∑|κ|=i∂κg(mj)κ!(x−mj)κP_{M_{j}}(x) = \sum_{i = 0}^{\lfloor k_{g}\rfloor_{*}}{}\sum_{|\kappa| = i}^{}{}\frac{\partial^{\kappa}g(m_{j})}{\kappa!}(x - m_{j})^{\kappa} 进行局部逼近。根据多元泰勒展开定理与 Assumption 1,多项式逼近误差被严格控制在 O(L/Mkg)O(L/M^{k_{g}})。

(3)残差局部多项式岭回归与桥接优化误差控制:在第 tt 期,当观察到 xt∈Mjx_{t} \in M_{j} 时,利用落入该网格的历史残差 ds−psb̂td_{s} - p_{s}{\hat{b}}_{t} 进行局部多项式岭回归得到多项式系数估计 θ̂t,j{\hat{\theta}}_{t,j},进而计算无约束贪婪价格 ptu=−θ̂t,j⊤φMj(xt)/(2b̂t)p_{t}^{u} = - {\hat{\theta}}_{t,j}^{\top}\varphi_{M_{j}}(x_{t})/(2{\hat{b}}_{t}) 并投影到安全区间 [p_+δt,p¯−δt]\lbrack\underline{p} + \delta_{t},\overline{p} - \delta_{t}\rbrack。为了克服残差中含有 b̂t−b{\hat{b}}_{t} - b 估计偏误的困难,作者引入了创新的桥接优化问题(Bridging Optimization Problem)与清洁高概率事件分析,成功控制了误差的级联传播。

3.3 定理 1(ASMLPE 后悔上界)数理推导与深入解析

【Theorem 1:SMLPE 模型的后悔上界定理】

【Theorem 1 (ASMLPE 累积后悔上界)】:假设需求函数满足 Assumption 1,令算法 ASMLPE 设置网格参数 M=⌈(L2T)1/(d+2kg)⌉M = \lceil(L^{2}T)^{1/(d + 2k_{g})}\rceil,则算法在时域 TT 内的期望累积后悔满足:

RASMLPE(T)=O∼(T∨(L2T)dd+2kg)R^{ASMLPE}(T) = \overset{\sim}{O}\left( \sqrt{T} \vee (L^{2}T)^{\frac{d}{d + 2k_{g}}} \right)

其中 O∼(⋅)\overset{\sim}{O}( \cdot ) 隐藏了常数因子以及 logT\log T 的对数多项式阶。

数理推导详解:【定理 1 核心数理推导过程拆解】:

Step 1:单期后悔分解与强凹性展开。由于收益函数 r(p,x)=p(bp+g(x))r(p,x) = p(bp + g(x)) 关于 pp 具有强凹模数 2|b|2|b|,对于任意价格 pt=ptg+Δtp_{t} = p_{t}^{g} + \Delta_{t},利用二阶泰勒展开可将单期期望收益差精确转化为价格偏误的二次项:

r(pt*,xt)−r(pt,xt)=|b|(pt−pt*)2≤2|b|(ptg−pt*)2+2|b|Δt2r(p_{t}^{*},x_{t}) - r(p_{t},x_{t}) = |b|(p_{t} - p_{t}^{*})^{2} \leq 2|b|(p_{t}^{g} - p_{t}^{*})^{2} + 2|b|\Delta_{t}^{2}

其中随机冲击部分带来的累积后悔为 ∑t=1T2|b|δt2=2|b|∑t=1Tt−1/2=O(T)\sum_{t = 1}^{T}{}2|b|\delta_{t}^{2} = 2|b|\sum_{t = 1}^{T}{}t^{- 1/2} = O(\sqrt{T})。

Step 2:贪婪价格偏误界定。利用贪婪价格计算公式 ptg≈−ĝt(xt)/(2b̂t)p_{t}^{g} \approx - {\hat{g}}_{t}(x_{t})/(2{\hat{b}}_{t}) 与最优价格 pt*=−g(xt)/(2b)p_{t}^{*} = - g(x_{t})/(2b),经过代数放缩可得:

|ptg−pt*|≤|ĝt(xt)−g(xt)|2|b̂t|+|g(xt)|⋅|b̂t−b|2|bb̂t|≤C1|ĝt(xt)−g(xt)|+C2|b̂t−b||p_{t}^{g} - p_{t}^{*}| \leq \frac{|{\hat{g}}_{t}(x_{t}) - g(x_{t})|}{2|{\hat{b}}_{t}|} + \frac{|g(x_{t})| \cdot |{\hat{b}}_{t} - b|}{2|b{\hat{b}}_{t}|} \leq C_{1}|{\hat{g}}_{t}(x_{t}) - g(x_{t})| + C_{2}|{\hat{b}}_{t} - b|

Step 3:斜率估计误差界。利用鞅差序列的大偏差集中不等式(Azuma-Hoeffding / Freedman 不等式),由于分母 ∑s=1tΔs2≈∑s=1ts−1/2≈2t\sum_{s = 1}^{t}{}\Delta_{s}^{2} \approx \sum_{s = 1}^{t}{}s^{- 1/2} \approx 2\sqrt{t},分子为均值为 0 的鞅,因此以极高概率有 |b̂t−b|≤O∼(t−1/4)|{\hat{b}}_{t} - b| \leq \overset{\sim}{O}(t^{- 1/4})。对其平方求和即可贡献 O∼(∑t=1Tt−1/2)=O∼(T)\overset{\sim}{O}(\sum_{t = 1}^{T}{}t^{- 1/2}) = \overset{\sim}{O}(\sqrt{T})。

Step 4:情境函数估计误差界与网格平衡。对于落入网格 MjM_{j} 的情境,函数估计误差由两部分构成:一是多项式逼近偏差(Bias)O(L/Mkg)O(L/M^{k_{g}}),二是局部最小二乘估计的统计方差(Variance)O∼(Md/Tt,j)\overset{\sim}{O}(\sqrt{M^{d}/T_{t,j}})。对所有网格求和并在时域累加,得到非参数学习的累积后悔为 O∼(T⋅(L/Mkg)2+Md)\overset{\sim}{O}(T \cdot (L/M^{k_{g}})^{2} + M^{d})。取 M=⌈(L2T)1/(d+2kg)⌉M = \lceil(L^{2}T)^{1/(d + 2k_{g})}\rceil 平衡偏误与方差,立即可得 O∼((L2T)d/(d+2kg))\overset{\sim}{O}((L^{2}T)^{d/(d + 2k_{g})})。结合两部分,定理 1 得证。

3.4 引理 1 与定理 2(SMLPE 后悔极小化极大下界)数理证明与构造分析

【Lemma 1 & Theorem 2:SMLPE 模型困难实例与后悔下界】

【Lemma 1 (困难实例函数的光滑性保证)】:对 0<kg≤10 < k_{g} \leq 1 和 L>0L > 0,令 M=⌈(L2T)1/(d+2kg)⌉M = \lceil(L^{2}T)^{1/(d + 2k_{g})}\rceil。对任意二值向量 w∈{0,1}Mdw \in \{ 0,1\}^{M^{d}},定义在网格 MjM_{j} 内的函数:

{|b|(p_+p¯),ifwj=0|b|(p_+p¯)+L2(D(x,∂Mj))kg,ifwj=1andD(x,∂Mj)≤141/kgM|b|(p_+p¯)+L8Mkg,ifwj=1andD(x,∂Mj)>141/kgM\left\{ \begin{matrix} |b|(\underline{p} + \overline{p}), & if\ w_{j} = 0 \\ |b|(\underline{p} + \overline{p}) + \frac{L}{2}(D(x,\partial M_{j}))^{k_{g}}, & if\ w_{j} = 1\ and\ D(x,\partial M_{j}) \leq \frac{1}{4^{1/k_{g}}M} \\ |b|(\underline{p} + \overline{p}) + \frac{L}{8M^{k_{g}}}, & if\ w_{j} = 1\ and\ D(x,\partial M_{j}) > \frac{1}{4^{1/k_{g}}M} \end{matrix} \right.\

其中 ∂Mj\partial M_{j} 为网格边界,D(x,∂Mj)≜infy∈∂Mj∥x−y∥2D(x,\partial M_{j}) \triangleq \inf_{y \in \partial M_{j}}\| x - y\|_{2}。则 gw(x)g_{w}(x) 严格属于 𝒢d(kg,L)\mathcal{G}_{d}(k_{g},L)。

【Theorem 2 (SMLPE 极小化极大后悔下界)】:对于 0<kg≤1,L>0,[b_,b¯]⊂(−∞,0)0 < k_{g} \leq 1,L > 0,\lbrack\underline{b},\overline{b}\rbrack \subset ( - \infty,0) 及 σ≥0\sigma \geq 0,存在与 TT 和 LL 无关的常数 K1>0K_{1} > 0,使得对任意容许策略 π\pi:

supg∈𝒢d(kg,L),b∈[b_,b¯],𝒫,𝒬∈ℰ(σ)Rg,b,𝒫,𝒬π(T)≥K1⋅(T∨(L2T)dd+2kg)\sup_{g \in \mathcal{G}_{d}(k_{g},L),b \in \lbrack\underline{b},\overline{b}\rbrack,\mathcal{P},\mathcal{Q} \in \mathcal{E}(\sigma)}R_{g,b,\mathcal{P},\mathcal{Q}}^{\pi}(T) \geq K_{1} \cdot \left( \sqrt{T} \vee (L^{2}T)^{\frac{d}{d + 2k_{g}}} \right)

数理推导详解:【定理 2 证明与信息论下界推导精要】:

1. 【Ω(T)\Omega(\sqrt{T}) 下界来源】:当 g(x)g(x) 为常数函数时,模型退化为纯线性需求模型。根据 Keskin and Zeevi (2014) 的 Theorem 1,由于无法预知真实的斜率参数 bb,任何算法在估计未知单参数时都必须承受至少 Ω(T)\Omega(\sqrt{T}) 的信息论下界。

2. 【Ω((L2T)d/(d+2kg))\Omega((L^{2}T)^{d/(d + 2k_{g})}) 下界来源】:基于 Lemma 1 构造的 2Md2^{M^{d}} 个假设函数族。考虑两个仅在第 jj 个分量上不同的索引 w=(w−j,0)w = (w_{- j},0) 和 w′=(w−j,1)w\prime = (w_{- j},1)。在网格 MjM_{j} 内,若 wj=0w_{j} = 0,真实最优价格为 p*=(p_+p¯)/2p^{*} = (\underline{p} + \overline{p})/2;若 wj=1w_{j} = 1,真实最优价格为 p*=(p_+p¯)/2+L/(16Mkg|b|)p^{*} = (\underline{p} + \overline{p})/2 + L/(16M^{k_{g}}|b|)。若算法出价偏离中点超过临界阈值 L/(32Mkg|b|)L/(32M^{k_{g}}|b|),则在至少一个真实函数下将产生至少 Ω((L/Mkg)2)\Omega((L/M^{k_{g}})^{2}) 的单期损失。

3. 【Bretagnolle-Huber 不等式与 KL 散度约束】:通过 Bretagnolle-Huber 不等式,判别两个环境的最小两类错误概率之和被两概率测度的 Kullback-Leibler (KL) 散度所下界约束:ℙ0(error)+ℙ1(error)≥12exp(−KL(ℙ0T∥ℙ1T)){\mathbb{P}}_{0}(error) + {\mathbb{P}}_{1}(error) \geq \frac{1}{2}\exp( - KL({\mathbb{P}}_{0}^{T} \parallel {\mathbb{P}}_{1}^{T}))。由于在单个网格内单期观测的 KL 散度为 (g1(x)−g0(x))22σ2≤L2128M2kgσ2\frac{(g_{1}(x) - g_{0}(x))^{2}}{2\sigma^{2}} \leq \frac{L^{2}}{128M^{2k_{g}}\sigma^{2}},在 MdM^{d} 个网格上积分并对先验平均,导出算法的总后悔必达到 Ω(Md⋅TMd⋅L2M2kg⋅exp(−CTL2Md+2kg))=Ω((L2T)d/(d+2kg))\Omega(M^{d} \cdot \frac{T}{M^{d}} \cdot \frac{L^{2}}{M^{2k_{g}}} \cdot \exp( - C\frac{TL^{2}}{M^{d + 2k_{g}}})) = \Omega((L^{2}T)^{d/(d + 2k_{g})})。

【结论与理论意义】:Theorem 2 证明了 Algorithm 1 (ASMLPE) 所达到的后悔上界在时域 TT 和光滑常数 LL 的依赖关系上均是严格紧致且不可改进的(Rate-optimal),建立了 0<kg≤10 < k_{g} \leq 1 常见 Hölder 连续函数类上的理论极限!

四、模型二 —— 线性上下文效应可分模型 (SMLCE)

4.1 模型设定与任意实数阶光滑价格函数类

在 SMLCE(Separable Model with Linear Contextual Effect)中,情境效应为线性参数模型,而价格效应为任意实数阶光滑的一般非参数函数。需求模型具体为:

Dt(p)=f(p)+a⊤xt+εt,∀p∈[p_,p¯]D_{t}(p) = f(p) + a^{\top}x_{t} + \varepsilon_{t},\forall p \in \lbrack\underline{p},\overline{p}\rbrack

其中 a∈ℝda \in {\mathbb{R}}^{d} 为未知的线性情境权重向量,满足范数有界条件 ∥a∥2≤a¯\| a\|_{2} \leq \overline{a}。函数 f:[p_,p¯]→ℝ+f:\lbrack\underline{p},\overline{p}\rbrack \rightarrow {\mathbb{R}}_{+} 为未知的价格需求函数,其满足以下定义:

【Assumption 2:价格需求函数的高阶光滑性假设】

【Assumption 2 (kfk_{f} 阶光滑价格函数类)】:称函数 f:[p_,p¯]→ℝ+f:\lbrack\underline{p},\overline{p}\rbrack \rightarrow {\mathbb{R}}_{+} 属于 kfk_{f} 阶光滑函数类 ℱ(kf,δ)\mathcal{F}(k_{f},\delta)(光滑常数 δ>0\delta > 0),若 f(⋅)f( \cdot ) 在 [p_,p¯]\lbrack\underline{p},\overline{p}\rbrack 上 ⌊kf⌋*\lfloor k_{f}\rfloor_{*} 次连续可微,且对任意 p,p′∈[p_,p¯]p,p\prime \in \lbrack\underline{p},\overline{p}\rbrack,满足:

|f(⌊kf⌋*)(p)−f(⌊kf⌋*)(p′)|≤δ|p−p′|kf−⌊kf⌋*|f^{(\lfloor k_{f}\rfloor_{*})}(p) - f^{(\lfloor k_{f}\rfloor_{*})}(p\prime)| \leq \delta|p - p\prime|^{k_{f} - \lfloor k_{f}\rfloor_{*}}

【重要理论特性】:与传统文献普遍要求的收益函数凹性假设(Concavity)不同,本文在 SMLCE 下完全不要求收益函数具有单峰或凹性,允许价格响应函数存在多峰(Multimodal)等极其复杂的形态。

4.2 ASMLCE 算法设计与“乐观之上的二次乐观”(Optimism over OFU)

算法机制与核心突破:针对 SMLCE 模型,作者提出了算法 ASMLCE(Algorithm 2)及回归核心 RLC(Algorithm 3):

(1)价格空间等距分段与局部多项式参数化:将价格区间 [p_,p¯]\lbrack\underline{p},\overline{p}\rbrack 划分为 NN 个等长小区间 Ij=[aj,bj)I_{j} = \lbrack a_{j},b_{j})。在每个区间 IjI_{j} 内,将 f(p)f(p) 近似为其 ⌊kf⌋*\lfloor k_{f}\rfloor_{*} 阶泰勒展开式 PIj(p)=∑i=0⌊kf⌋*f(i)(aj)i!(p−aj)iP_{I_{j}}(p) = \sum_{i = 0}^{\lfloor k_{f}\rfloor_{*}}{}\frac{f^{(i)}(a_{j})}{i!}(p - a_{j})^{i}。逼近偏误严格有界:|f(p)−PIj(p)|≤δ(p¯−p_)kf/Nkf≜Δ|f(p) - P_{I_{j}}(p)| \leq \delta(\overline{p} - \underline{p})^{k_{f}}/N^{k_{f}} \triangleq \Delta。

(2)特征联合映射与有偏线性上下文岭回归:在区间 IiI_{i} 内,将需求模型改写为增广线性形式:Dt(p)=⟨ui,φ(p)⟩+⟨a,xt⟩+βt=⟨θi,φ(p,xt)⟩+βtD_{t}(p) = \langle u_{i},\varphi(p)\rangle + \langle a,x_{t}\rangle + \beta_{t} = \langle\theta_{i},\varphi(p,x_{t})\rangle + \beta_{t},其中 φ(p,xt)=(φ(p)⊤,xt⊤)⊤\varphi(p,x_{t}) = (\varphi(p)^{\top},x_{t}^{\top})^{\top},θi=(ui⊤,a⊤)⊤\theta_{i} = (u_{i}^{\top},a^{\top})^{\top}。扰动项 βt=f(pt)−PIi(pt)+εt\beta_{t} = f(p_{t}) - P_{I_{i}}(p_{t}) + \varepsilon_{t} 含有受控有界偏误 |βt−εt|≤Δ|\beta_{t} - \varepsilon_{t}| \leq \Delta。通过调用 RLC 求解带正则化的岭回归估计量 (ût,i,ât,i)({\hat{u}}_{t,i},{\hat{a}}_{t,i}) 与设计矩阵 Vt,iV_{t,i}。

(3)“在 OFU 之上的二次乐观”(Optimism over OFU):在无情境的多模态动态定价中(如 Wang et al. 2021a),最优价格是静态常数,因此可将每个价格区间视为多臂老虎机的一个臂;但在情境定价中,最优价格随随机情境 xtx_{t} 每期剧烈变动。为此,算法在每个价格区间 IiI_{i} 内,首先应用面对不确定性时的乐观原则(OFU),计算区间内的乐观价格 p̂t,i{\hat{p}}_{t,i} 与乐观收益上界 r̂t,i{\hat{r}}_{t,i}(引入置信椭球半径 γt,i\gamma_{t,i} 与偏误补偿项 Δ\Delta):

r̂t,i≜maxp∈Iip×(⟨ût,i,φ(p)⟩+⟨ât,i,xt⟩+γt,iφ(p,xt)⊤Vt,i−1φ(p,xt)+Δ){\hat{r}}_{t,i} \triangleq \max_{p \in I_{i}}p \times (\langle{\hat{u}}_{t,i},\varphi(p)\rangle + \langle{\hat{a}}_{t,i},x_{t}\rangle + \gamma_{t,i}\sqrt{\varphi(p,x_{t})^{\top}V_{t,i}^{- 1}\varphi(p,x_{t})} + \Delta)

随后,算法在全部 NN 个价格区间的乐观收益中进行“二次全局挑优”:选取使 r̂t,i{\hat{r}}_{t,i} 达到全局最大的区间 it=argmax1≤i≤Nr̂t,ii_{t} = \arg\max_{1 \leq i \leq N}{\hat{r}}_{t,i},并执行价格 pt=p̂t,itp_{t} = {\hat{p}}_{t,i_{t}}。这一机制确保了即使最优价格随情境动态跳跃,算法也能以高概率追踪到全知最优收益!

4.3 定理 3(ASMLCE 后悔上界)数理推导与深入解析

【Theorem 3:SMLCE 模型的后悔上界定理】

【Theorem 3 (ASMLCE 累积后悔上界)】:假设需求函数满足 Assumption 2,令算法 ASMLCE 设置区间数 N=⌈δ2/(2kf+1)T1/(2kf+1)⌉+1N = \lceil\delta^{2/(2k_{f} + 1)}T^{1/(2k_{f} + 1)}\rceil + 1 且 Δ=δ(p¯−p_)kf/Nkf\Delta = \delta(\overline{p} - \underline{p})^{k_{f}}/N^{k_{f}},则算法的累积期望后悔满足:

RASMLCE(T)=O∼(d(T∨(δTkf+1)12kf+1))R^{ASMLCE}(T) = \overset{\sim}{O}\left( d\left( \sqrt{T} \vee (\delta T^{k_{f} + 1})^{\frac{1}{2k_{f} + 1}} \right) \right)

当 δ=Θ(1)\delta = \Theta(1) 时,后悔界简化为 O∼(dT(kf+1)/(2kf+1))\overset{\sim}{O}(dT^{(k_{f} + 1)/(2k_{f} + 1)})。

数理推导详解:【定理 3 数理推导核心逻辑】:

1. 【自正则化鞅集中不等式】:利用 Abbasi-Yadkori et al. (2011) 的自正则化鞅大偏差理论,在存在有界扰动偏误 Δ\Delta 的条件下,建立参数误差的置信椭球界:以至少 1−1/T21 - 1/T^{2} 的高概率,对所有 t≥1t \geq 1 和区间 ii,均有 ∥θ̂t,i−θi*∥Vt,i≤γt,i\|{\hat{\theta}}_{t,i} - \theta_{i}^{*}\|_{V_{t,i}} \leq \gamma_{t,i}。

2. 【乐观收益控制】:由“乐观之上的乐观”决策机制,令包含最优价格 pt*p_{t}^{*} 的区间为 i*i^{*},则必然有 r̂t,it≥r̂t,i*≥r(pt*,xt){\hat{r}}_{t,i_{t}} \geq {\hat{r}}_{t,i^{*}} \geq r(p_{t}^{*},x_{t})。因此单期后悔直接被选定区间的置信宽度所控制:

r(pt*,xt)−r(pt,xt)≤r̂t,it−r(pt,xt)≤2p¯(γt,it∥φ(pt,xt)∥Vt,it−1+Δ)r(p_{t}^{*},x_{t}) - r(p_{t},x_{t}) \leq {\hat{r}}_{t,i_{t}} - r(p_{t},x_{t}) \leq 2\overline{p}\left( \gamma_{t,i_{t}}\|\varphi(p_{t},x_{t})\|_{V_{t,i_{t}}^{- 1}} + \Delta \right)

3. 【椭球势能求和引理(Ellipsoidal Potential Lemma)】:利用经典的矩阵行列式引理,对每个区间 ii 内被选中的样本点求和,有 ∑t:it=i∥φ(pt,xt)∥Vt,i−1≤|Di|⋅2(b(kf)+d)log(1+|Di|/λ)\sum_{t:i_{t} = i}^{}{}\|\varphi(p_{t},x_{t})\|_{V_{t,i}^{- 1}} \leq \sqrt{|D_{i}| \cdot 2(b(k_{f}) + d)\log(1 + |D_{i}|/\lambda)}。对全部 NN 个区间应用柯西-施瓦茨不等式求和,总不确定性贡献为 O∼(dNT)\overset{\sim}{O}(d\sqrt{NT})。加上偏误累积项 2p¯TΔ=O(Tδ/Nkf)2\overline{p}T\Delta = O(T\delta/N^{k_{f}}),总后悔界为 O∼(dNT+Tδ/Nkf)\overset{\sim}{O}(d\sqrt{NT} + T\delta/N^{k_{f}})。

4. 【最优段数选取】:令 dNT≈Tδ/Nkfd\sqrt{NT} \approx T\delta/N^{k_{f}},解得最优段数 N=Θ((δ2T)1/(2kf+1))N = \Theta((\delta^{2}T)^{1/(2k_{f} + 1)}),代入即得紧致后悔上界 O∼(d(δTkf+1)1/(2kf+1))\overset{\sim}{O}(d(\delta T^{k_{f} + 1})^{1/(2k_{f} + 1)})。

4.4 引理 2 与定理 4(SMLCE 后悔下界)数理证明与构造分析

【Lemma 2 & Theorem 4:SMLCE 模型困难实例与后悔下界】

【Lemma 2 (无限光滑阶跃 Bump 函数与候选需求族)】:设 u(x)=exp(−1/(x(1−x)))𝟷{x∈[0,1]}u(x) = \exp( - 1/(x(1 - x)))\mathtt{1}_{\{ x \in \lbrack 0,1\rbrack\}},定义归一化平滑阶跃函数 gkf(x)=1Zkf(S(x)𝟷{x≤1}+S(2−x)𝟷{x>1})g_{k_{f}}(x) = \frac{1}{Z_{k_{f}}}(S(x)\mathtt{1}_{\{ x \leq 1\}} + S(2 - x)\mathtt{1}_{\{ x > 1\}})。将价格区间划分为 JJ 个等长段 Ij=[aj,bj)I_{j} = \lbrack a_{j},b_{j}),定义需求函数族:

{δ̂2p,ifp∉Ijδ̂+η⋅gkf(2J(p−aj))2p,ifp∈Ij\left\{ \begin{matrix} \frac{\hat{\delta}}{2p}, & if\ p \notin I_{j} \\ \frac{\hat{\delta} + \eta \cdot g_{k_{f}}(2J(p - a_{j}))}{2p}, & if\ p \in I_{j} \end{matrix} \right.\

适当选取缩放系数 δ̂,J\hat{\delta},J 和扰动幅度 η\eta,可保证每一个 fj(p)f_{j}(p) 均严格属于 ℱ(kf,δ)\mathcal{F}(k_{f},\delta)。

【Theorem 4 (SMLCE 极小化极大后悔下界)】:对于 kf>0,δ>0,a¯≥0k_{f} > 0,\delta > 0,\overline{a} \geq 0 及 σ≥0\sigma \geq 0,存在与 TT 和 δ\delta 无关的常数 K2>0K_{2} > 0,使得对任意容许策略 π\pi:

supf∈ℱ(kf,δ),∥a∥2≤a¯,𝒫,𝒬∈ℰ(σ)Rf,a,𝒫,𝒬π(T)≥K2⋅(T∨(δTkf+1)12kf+1)\sup_{f \in \mathcal{F}(k_{f},\delta),\| a\|_{2} \leq \overline{a},\mathcal{P},\mathcal{Q} \in \mathcal{E}(\sigma)}R_{f,a,\mathcal{P},\mathcal{Q}}^{\pi}(T) \geq K_{2} \cdot \left( \sqrt{T} \vee (\delta T^{k_{f} + 1})^{\frac{1}{2k_{f} + 1}} \right)

数理推导详解:【定理 4 证明精要与信息论测度变换】:

设基准环境为 f0(p)=δ̂/(2p)f_{0}(p) = \hat{\delta}/(2p),此时任意价格下的期望收益恒为 δ̂/2\hat{\delta}/2;而在扰动环境 fj(p)f_{j}(p) 下,仅在区间 IjI_{j} 内收益隆起高度为 η/2\eta/2。根据 Pinsker 不等式,算法在基准测度 ℙ0{\mathbb{P}}_{0} 与扰动测度 ℙj{\mathbb{P}}_{j} 下对区间 IjI_{j} 的访问次数期望之差满足:

|𝔼jπ[Tj]−𝔼0π[Tj]|≤12TKL(ℙ0π∥ℙjπ)≤14σ𝔼0π[Tj]⋅Tη|{\mathbb{E}}_{j}^{\pi}\lbrack T_{j}\rbrack - {\mathbb{E}}_{0}^{\pi}\lbrack T_{j}\rbrack| \leq \frac{1}{2}T\sqrt{KL({\mathbb{P}}_{0}^{\pi} \parallel {\mathbb{P}}_{j}^{\pi})} \leq \frac{1}{4\sigma}\sqrt{{\mathbb{E}}_{0}^{\pi}\lbrack T_{j}\rbrack} \cdot T\eta

选取在基准策略下访问次数最少的区间 j*=argmin1≤j≤J𝔼0π[Tj]j^{*} = \arg\min_{1 \leq j \leq J}{\mathbb{E}}_{0}^{\pi}\lbrack T_{j}\rbrack,必然有 𝔼0π[Tj*]≤T/J{\mathbb{E}}_{0}^{\pi}\lbrack T_{j^{*}}\rbrack \leq T/J。通过设置参数 η=Θ((δT−kf)1/(2kf+1))\eta = \Theta((\delta T^{- k_{f}})^{1/(2k_{f} + 1)}) 和 J=Θ((δ2T)1/(2kf+1))J = \Theta((\delta^{2}T)^{1/(2k_{f} + 1)}),保证了在真实环境为 fj*f_{j^{*}} 时,算法在区间 Ij*I_{j^{*}} 的探索次数不超过 T/2T/2,从而在剩余至少 T/2T/2 期中错失最优峰值收益,导致累积后悔至少为 Ω(Tη)=Ω((δTkf+1)1/(2kf+1))\Omega(T\eta) = \Omega((\delta T^{k_{f} + 1})^{1/(2k_{f} + 1)})。定理 4 得证。

五、模型三 —— 全非参数可分模型 (SMNPE)

5.1 模型设定与双重非参数光滑性假设

SMNPE(Separable Model with Nonparametric Effects)是可分需求模型的最一般形式,价格效应与情境效应同时具备未知的非参数高阶光滑结构:

Dt(p)=f(p)+g(xt)+εt,∀p∈[p_,p¯]D_{t}(p) = f(p) + g(x_{t}) + \varepsilon_{t},\forall p \in \lbrack\underline{p},\overline{p}\rbrack

【Assumption 3:双重高阶非参数光滑性】

【Assumption 3 (双重非参数光滑性假设)】:价格函数 f(⋅)∈ℱ(kf,δ)f( \cdot ) \in \mathcal{F}(k_{f},\delta) 满足 Assumption 2;情境函数 g(⋅)∈𝒢d(kg,L)g( \cdot ) \in \mathcal{G}_{d}(k_{g},L) 满足 Assumption 1。为保证模型的可识别性(Identifiability),不失一般性地假设 𝔼[g(x)]=0{\mathbb{E}}\lbrack g(x)\rbrack = 0(若均值非零,可将其常数部分直接吸收到 f(p)f(p) 中)。

5.2 ASMNPE 在线学习算法:探索与利用两阶段架构

两阶段算法运行逻辑:当 f(p)f(p) 与 g(x)g(x) 均为一般非参数函数时,SMLPE 的随机冲击法与 SMLCE 的联合岭回归法均无法直接套用。为此,作者提出了基于“探索-利用严格分离”(Explore-then-Exploit)的 ASMNPE(Algorithm 4):

(1)价格空间离散化与网格构建:将价格区间 [p_,p¯]\lbrack\underline{p},\overline{p}\rbrack 等分为 NN 份,取各段中点构成离散候选价格集 𝒫N={P[1],P[2],…,P[N]}\mathcal{P}_{N} = \{ P\lbrack 1\rbrack,P\lbrack 2\rbrack,\ldots,P\lbrack N\rbrack\},其中 P[i]=p_+p¯−p_N(i−1/2)P\lbrack i\rbrack = \underline{p} + \frac{\overline{p} - \underline{p}}{N}(i - 1/2)。将情境空间 [0,1]d\lbrack 0,1\rbrack^{d} 划分为 MdM^{d} 个小网格 MjM_{j}。

(2)纯探索阶段(前 n0Nn_{0}N 期):算法在全部 NN 个离散价格上开展完全均匀的轮换定价实验,每个候选价格精确试验 n0n_{0} 次。利用样本均值法纯净估计每个价格点上的基准需求 f(P[i])f(P\lbrack i\rbrack)(记为 sis_{i}):

si≜1n0∑t:pt=P[i],t≤n0NDt(pt)=f(P[i])+1n0∑g(xt)+1n0∑εts_{i} \triangleq \frac{1}{n_{0}}\sum_{t:p_{t} = P\lbrack i\rbrack,t \leq n_{0}N}^{}{}D_{t}(p_{t}) = f(P\lbrack i\rbrack) + \frac{1}{n_{0}}\sum_{}^{}{}g(x_{t}) + \frac{1}{n_{0}}\sum_{}^{}{}\varepsilon_{t}

由于在探索阶段定价完全独立于情境 xtx_{t},且 𝔼[g(xt)]=0{\mathbb{E}}\lbrack g(x_{t})\rbrack = 0,大数定律与亚高斯集中性保证了 sis_{i} 是 f(P[i])f(P\lbrack i\rbrack) 的高精度无偏估计,估计误差为 O∼(1/n0)\overset{\sim}{O}(1/\sqrt{n_{0}}),彻底解耦了未知情境效应 g(x)g(x) 的干扰!

(3)纯利用阶段(第 n0N+1∼Tn_{0}N + 1 \sim T 期):固定价格估计值 sis_{i} 不再更新。在每一期 tt,观察到 xt∈Mjx_{t} \in M_{j},利用利用阶段落入该网格的历史残差 ds−∑i=1Nsi𝟷{ps=P[i]}d_{s} - \sum_{i = 1}^{N}{}s_{i}\mathtt{1}_{\{ p_{s} = P\lbrack i\rbrack\}} 进行局部多项式回归,得到情境效应估计 ĝt(xt)=ût,j⊤φMj(xt){\hat{g}}_{t}(x_{t}) = {\hat{u}}_{t,j}^{\top}\varphi_{M_{j}}(x_{t})。随后在 NN 个候选价格中选取最大化代理期望收益的价格:

pt=argmaxP[i]∈𝒫NP[i]⋅(si+ût,j⊤φMj(xt))p_{t} = \arg\max_{P\lbrack i\rbrack \in \mathcal{P}_{N}}P\lbrack i\rbrack \cdot (s_{i} + {\hat{u}}_{t,j}^{\top}\varphi_{M_{j}}(x_{t}))

5.3 定理 5 与定理 6(SMNPE 后悔界与紧致性分析)

【Theorem 5 & Theorem 6:SMNPE 模型后悔上界与下界】

【Theorem 5 (ASMNPE 累积后悔上界)】:在 Assumption 3 下,令算法参数取 N=⌈T1/4⌉,n0=⌈T1/2⌉,M=⌈(L2T)1/(d+2kg)⌉N = \lceil T^{1/4}\rceil,n_{0} = \lceil T^{1/2}\rceil,M = \lceil(L^{2}T)^{1/(d + 2k_{g})}\rceil,则 ASMNPE 的期望累积后悔满足:

RASMNPE(T)=O∼(T3/4∨(LdTd+kg)1d+2kg)R^{ASMNPE}(T) = \overset{\sim}{O}\left( T^{3/4} \vee (L^{d}T^{d + k_{g}})^{\frac{1}{d + 2k_{g}}} \right)

【Theorem 6 (SMNPE 极小化极大后悔下界)】:对于 0<kg≤1,L>0,kf>0,δ>00 < k_{g} \leq 1,L > 0,k_{f} > 0,\delta > 0 及 σ≥0\sigma \geq 0,存在常数 K3>0K_{3} > 0,使得对任意容许策略 π\pi:

supf∈ℱ(kf,δ),g∈𝒢d(kg,L),𝒫,𝒬∈ℰ(σ)Rf,g,𝒫,𝒬π(T)≥K3⋅(LdTd+kg)1d+2kg\sup_{f \in \mathcal{F}(k_{f},\delta),g \in \mathcal{G}_{d}(k_{g},L),\mathcal{P},\mathcal{Q} \in \mathcal{E}(\sigma)}R_{f,g,\mathcal{P},\mathcal{Q}}^{\pi}(T) \geq K_{3} \cdot (L^{d}T^{d + k_{g}})^{\frac{1}{d + 2k_{g}}}

数理推导详解:【定理 5 与定理 6 核心推导与紧致性讨论】:

1. 【探索与利用的后悔平衡】:在探索阶段,由于未进行收益优化,共进行 n0N=T3/4n_{0}N = T^{3/4} 期探索,产生直接探索后悔 O(n0N)=O(T3/4)O(n_{0}N) = O(T^{3/4});在利用阶段,离散化价格网格导致逼近损失 O(T/N)=O(T3/4)O(T/N) = O(T^{3/4}),对 f(P[i])f(P\lbrack i\rbrack) 的估计方差导致损失 O(T/n0)=O(T3/4)O(T/\sqrt{n_{0}}) = O(T^{3/4}),而学习高维非参数情境函数 g(x)g(x) 产生的非参数估计后悔为 O∼((LdTd+kg)1/(d+2kg))\overset{\sim}{O}((L^{d}T^{d + k_{g}})^{1/(d + 2k_{g})})。四者权衡,导出 Theorem 5。

2. 【速率最优性区间】:当 L=Θ(1)L = \Theta(1) 且满足 d≥2d \geq 2 或 kg≤1/2k_{g} \leq 1/2 时,有 (d+kg)/(d+2kg)≥3/4(d + k_{g})/(d + 2k_{g}) \geq 3/4。此时定理 5 的上界与定理 6 的下界完全重合,最优后悔阶精确为 Θ∼(T(d+kg)/(d+2kg))\overset{\sim}{\Theta}(T^{(d + k_{g})/(d + 2k_{g})})!这表明在双重非参数设定下,高维情境函数的学习复杂度构成了整个在线定价系统的统计瓶颈。

六、三大模型跨维度系统对比与理论洞察

6.1 探索-利用权衡、解耦机制与数据利用率综合对比

下表 Table 3 系统总结了本文提出的三大算法在如何平衡探索-利用权衡、如何分离估计 f(p)f(p) 与 g(x)g(x) 以及历史数据利用率层面的核心异同:

模型配置与算法 如何平衡探索与利用权衡? 如何解耦估计 f(p)f(p) 与 g(x)g(x)? 历史数据利用效率与范围
SMLPE (Algorithm 1:
ASMLPE)
随机价格冲击 (Random Shock: Δt=±t−1/4\Delta_{t} = \pm t^{- 1/4}) 随机冲击独立估计 bb;残差回归估计 g(x)g(x) 全时域数据估计 bb;各网格数据局部估计 g(x)g(x)
SMLCE (Algorithm 2:
ASMLCE)
面对不确定性时的乐观原则 (OFU) + 二次选优 局部多项式逼近 + 增广特征联合岭回归 价格区间内部数据局部估计 (f,a)(f,a)
SMNPE (Algorithm 4:
ASMNPE)
探索与利用两阶段分离 (Explore-then-Exploit) 探索期纯随机实验估计 f(p)f(p);利用期残差估计 g(x)g(x) 探索期数据估计 f(P[i])f(P\lbrack i\rbrack);利用期网格数据估计 g(x)g(x)

6.2 可分结构的降维威力:对比全非参数与参数基准

结构价值体现:将本文结果与现有代表性文献对比,可分结构的降维优势显露无遗:

1. 【超越全非参数模型】:在 Lipschitz 条件下(kg=1,kf=1k_{g} = 1,k_{f} = 1 且 d≥2d \geq 2),完全非参数定价模型(Slivkins 2011)的最优后悔界为 Θ∼(T(d+2)/(d+3))\overset{\sim}{\Theta}(T^{(d + 2)/(d + 3)})。而在本文的可分模型中,SMLPE 的后悔界大幅缩减至 Θ∼(Td/(d+2))\overset{\sim}{\Theta}(T^{d/(d + 2)}),SMNPE 的后悔界为 Θ∼(T(d+1)/(d+2))\overset{\sim}{\Theta}(T^{(d + 1)/(d + 2)})。由于 dd+2<d+1d+2<d+2d+3\frac{d}{d + 2} < \frac{d + 1}{d + 2} < \frac{d + 2}{d + 3},可分结构在统计上消除了价格与情境交叉项带来的高阶复杂性,实现了指数级的样本效率提升!

2. 【克服参数模型误设】:相比 Nambiar et al. (2019) 仅能收敛到最佳线性近似模型的局限,本文算法直接收敛到真实的非参数最优定价策略,完全消除了稳态渐近偏误。

6.3 内生决策变量 p 与外生环境协变量 x 的统计非对称性

对比一维情境下的 SMLPE(d=1d = 1 时后悔为 Θ∼(T1/(1+2kg))\overset{\sim}{\Theta}(T^{1/(1 + 2k_{g})}))与 SMLCE(后悔为 Θ∼(T(kf+1)/(2kf+1))\overset{\sim}{\Theta}(T^{(k_{f} + 1)/(2k_{f} + 1)})),可以发现参数 kgk_{g} 与 kfk_{f} 扮演了截然不同的数学角色。这一深刻差异根源于变量生成机制的非对称性(Asymmetry):价格 ptp_{t} 是算法主动选取的内生决策变量,算法可以通过主动试探与设计置信区间进行定向信息搜集;而情境 xtx_{t} 是由大自然被动抽样的外生协变量,算法无法操控其取值分布,只能被动接收,因而外生非参数函数的学习复杂度显著高于内生价格函数!

七、全文数值仿真图谱深度解析

7.1 SMLPE 算法性能、敏感性与理论收敛率验证 (Figures 3, 4, 5)

在 SMLPE 数值实验中,作者评估了 ASMLPE、RPS 算法(Nambiar et al. 2019)与 Contextual Zooming 算法(Slivkins 2011)在两种设定下的表现:

Figure 3: SMLPE 模型下三种定价算法在设定 (19) 与设定 (20) 下的经验相对后悔对比

【Figure 3 解析】:在公式 (19) 的高非线性情境设定下(图 a),ASMLPE 的相对后悔随时间迅速下降至 1% 左右;而 RPS 算法由于只能拟合最佳线性模型,陷入了约 15% 的持续模型误设后悔;Contextual Zooming 算法由于未利用可分结构且受局部大导数干扰,表现极不稳定。在公式 (20) 设定下(图 b),ASMLPE 同样展现出远超基准算法的极速收敛优势。

Figure 4: SMLPE 模型下情境维度 d 与光滑度参数 kg 对 ASMLPE 经验相对后悔的影响

【Figure 4 解析】:图 (a) 显示,固定 kg=1k_{g} = 1,随着特征维度 dd 从 2 增加到 4 和 10,后悔曲线衰减变慢,验证了高维情境带来的学习复杂度上升;图 (b) 显示,固定 d=1d = 1,光滑度 kgk_{g} 从 0.2 提升至 0.6 时收敛速度显著跃升,但从 0.6 进一步提升至 1.0 时边际改善趋于平缓,完全符合定理 1 的理论预测。

Figure 5: SMLPE 模型下 ASMLPE 缩放经验后悔 (Scaled Regret) 的渐近常数收敛验证

【Figure 5 解析】:通过将经验后悔分别除以理论理论阶 t5/9t^{5/9}(针对 d=2,kg=0.8d = 2,k_{g} = 0.8)和 t5/6t^{5/6}(针对 d=10,kg=1d = 10,k_{g} = 1),缩放后悔曲线在时域增长时精确收敛为一条水平常数线,严格验证了 Theorem 1 理论上界的极度精确性!

7.2 SMLCE 算法对比、光滑度参数影响与收敛率 (Figures 6, 7)

Figure 6: SMLCE 模型下三种算法经验相对后悔对比及光滑常数 δ 敏感性分析

【Figure 6 解析】:在需求模型 Dt(pt)=−415δpt2.5+30+1d𝟷⊤xt+εtD_{t}(p_{t}) = - \frac{4}{15}\delta p_{t}^{2.5} + 30 + \frac{1}{d}\mathtt{1}^{\top}x_{t} + \varepsilon_{t} 下(kf=2.5k_{f} = 2.5),图 (a) 表明 ASMLCE 显著优于未利用可分结构的 Contextual Zooming 算法和简单的线性贪婪算法(Linear Greedy);图 (b) 表明随着光滑常数 δ\delta 从 3.75 降至 1.00,价格函数波动减小,算法后悔显著下降。

Figure 7: SMLCE 模型下 ASMLCE 缩放经验后悔 R(t)/t^(7/12) 的渐近收敛验证

【Figure 7 解析】:当 kf=2.5k_{f} = 2.5 时,理论最优后悔阶为 t(2.5+1)/(2×2.5+1)=t7/12t^{(2.5 + 1)/(2 \times 2.5 + 1)} = t^{7/12}。Figure 7 显示将经验后悔除以 t7/12t^{7/12} 后高度平稳地收敛于常数 0.18,实证确立了 Theorem 3 的严密性。

7.3 SMNPE 算法表现、理论阶收敛与高维衰减 (Figures 8, 9, 10)

Figure 8: SMNPE 双重非参数模型下三种算法的经验相对后悔对比

【Figure 8 解析】:在分段非参数价格效应与复杂情境效应下,ASMNPE 算法凭借探索与利用分离机制稳步收敛,而基准算法几乎完全无法收敛,相对后悔停滞在 20%~25% 的高位。

Figure 9: SMNPE 模型下 ASMNPE 缩放经验后悔的渐近常数收敛验证 (t^(7/9) 与 t^(11/12))

Figure 10: SMNPE 模型下维度 d 与光滑度 kg 对 ASMNPE 相对后悔的影响

【Figure 9 & 10 解析】:Figure 9 中缩放后悔分别对 t7/9t^{7/9}(d=2,kg=0.8d = 2,k_{g} = 0.8)和 t11/12t^{11/12}(d=10,kg=1d = 10,k_{g} = 1)的收敛,以及 Figure 10 中高维导致收敛减缓、光滑度提升加速收敛的规律,均再次完美印证了 Theorem 5 和 Theorem 6 的理论结论。

7.4 模型误设(不可分需求结构)下的鲁棒性表现 (Figure 11)

Figure 11: 在交互不可分需求模型 Dt(pt) = -(xt + 0.5)pt + g(xt) + εt 下 ASMLPE 的鲁棒性检验

【Figure 11 解析】:为了测试算法在现实世界中遭遇模型误设时的稳健性,作者构造了价格与情境强交互的不可分模型 Dt(pt)=−(xt+0.5)pt+g(xt)+εtD_{t}(p_{t}) = - (x_{t} + 0.5)p_{t} + g(x_{t}) + \varepsilon_{t}(此时价格敏感度随情境 xtx_{t} 动态变化)。实验结果表明,强行应用假设可分结构的 ASMLPE 算法依然表现出优异的鲁棒性:其相对后悔在 t>10000t > 10000 时被稳定压制在 4% 以下。虽然由于不可避免的模型误设无法无限收敛至 0,但 4% 以内的低误差证明了可分模型在近似复杂现实环境时具有极高的工程实用价值。

八、全文结论、管理启示与未来研究方向

8.1 全文主要结论总结

本文针对现代零售与电商平台海量低销量长尾商品的动态定价难题,从真实产业实践出发,系统构建了基于情境的可分需求模型在线学习理论体系。论文完整刻画了 SMLPE、SMLCE 与 SMNPE 三种结构配置下的统计复杂度,提出了三种创新算法,并建立了严谨匹配的后悔上下界,填补了收益管理与在线学习领域的重要理论空白。

8.2 企业动态定价系统与算法设计的管理启示

管理启示一:1. 【重视内生决策与外生特征的因果解耦】:零售平台在应用机器学习进行定价时,切忌盲目套用纯黑箱非参数树模型(如直接把价格作为普通特征输入 GBDT)。必须采用结构可分模型,将内生价格弹性与外生环境效应解耦,以恢复单调合理的价格弹性曲线,避免“价格不敏感”导致的定价灾难。

管理启示二:2. 【巧用小幅价格随机冲击实现主动学习】:对于长尾和冷启动商品,企业无需进行大幅度的破坏性降价实验,只需在日常定价中叠加小幅的随机价格扰动(如 t−1/4t^{- 1/4} 阶的微小波动),即可实现对价格敏感度的高效无偏估计,且不影响日常收益。

管理启示三:3. 【高维特征的精细化治理】:情境维度 dd 是决定非参数学习速度的核心瓶颈。企业在收集海量特征时应进行必要的降维与特征筛选,优先保留对需求有显著解释力的核心协变量,以降低在线学习的样本复杂度。

8.3 算法参数选取的工程实践策略

工程落地策略:在算法实施中,理论算法需要预先知晓光滑度参数 kgk_{g} 和 kfk_{f}。在实际业务中,企业可以通过以下三种工程策略进行参数定阶:

(1)历史成熟品类迁移:利用具有丰富历史数据的成熟同类商品,拟合离线模型并推断该品类的光滑度先验范围;

(2)信息准则模型选择:采用贝叶斯信息准则(BIC)或赤池信息准则(AIC)在不同多项式阶数之间进行自适应超参数选优;

(3)前向逐步回归检验(Forward Selection):逐阶增加多项式阶数,通过对最高阶系数的统计显著性 t 检验决定最佳截断阶数。

8.4 未来研究方向展望

1. 【可分性结构的严谨统计检验】:目前文献中尚缺乏针对高维非参数/半参数可分结构的严格假设检验方法(即需要多少样本才能以统计置信度判定真实需求函数是否满足可分性),这是极具价值的计量经济学与统计学前沿方向;

2. 【多产品联合定价与品类选品优化(Assortment Optimization)】:将单品可分模型拓展至具有替代与互补效应的多产品网络收益管理场景;

3. 【带库存约束的联合库存与动态定价】:在有限库存约束下,研究可分需求模型在非平稳环境下的双重学习与控制问题。

九、参考文献 (References)

Abbasi-Yadkori Y, Pál D, Szepesvári C (2011) Improved algorithms for linear stochastic bandits. *Advances in Neural Information Processing Systems*, 24: 2312–2320.

Auer P, Cesa-Bianchi N, Fischer P (2002) Finite-time analysis of the multiarmed bandit problem. *Machine Learning*, 47(2–3): 235–256.

Ban GY, Keskin NB (2021) Personalized dynamic pricing with machine learning: High-dimensional features and heterogeneous elasticity. *Management Science*, 67(9): 5549–5568.

Besbes O, Zeevi A (2009) Dynamic pricing without knowing the demand function: Risk bounds and near-optimal algorithms. *Operations Research*, 57(6): 1407–1420.

Besbes O, Zeevi A (2012) Blind network revenue management. *Operations Research*, 60(6): 1537–1550.

Besbes O, Zeevi A (2015) On the (surprising) sufficiency of linear models for dynamic pricing with demand learning. *Management Science*, 61(4): 723–739.

Bretagnolle J, Huber C (1979) Estimation des densités: Risque minimax. *Zeitschrift für Wahrscheinlichkeitstheorie und Verwandte Gebiete*, 47(2): 119–137.

Broder J, Rusmevichientong P (2012) Dynamic pricing under a general parametric choice model. *Operations Research*, 60(4): 965–980.

Chen N, Gallego G (2021) Nonparametric pricing analytics with customer covariates. *Operations Research*, 69(3): 974–984.

Chen N, Hu M (2023) Data-driven revenue management: The interplay of data, model, and decisions. *Service Science*, 15(2): 79–91.

Chen X, Jasin S, Shi C (2022) *The Elements of Joint Learning and Optimization in Operations Management*, Vol. 18 (Springer Nature, New York).

Chen Y, Shi C (2023) Network revenue management with online inverse batch gradient descent method. *Production and Operations Management*, 32(7): 2123–2137.

Chu W, Li L, Reyzin L, Schapire R (2011) Contextual bandits with linear payoff functions. *Proceedings of the 14th International Conference on Artificial Intelligence and Statistics*, PMLR, 15: 208–214.

Cohen MC, Lobel I, Paes Leme R (2020) Feature-based dynamic pricing. *Management Science*, 66(11): 4921–4943.

den Boer AV (2015) Dynamic pricing and learning: Historical origins, current research, and new directions. *Surveys in Operations Research and Management Science*, 20(1): 1–18.

Fan J, Guo Y, Yu M (2024) Policy optimization using semiparametric models for dynamic pricing. *Journal of the American Statistical Association*, 119(545): 552–564.

Gur Y, Momeni A, Wager S (2022) Smoothness-adaptive contextual bandits. *Operations Research*, 70(6): 3198–3216.

Hu Y, Kallus N, Mao X (2022) Smooth contextual bandits: Bridging the parametric and nondifferentiable regret regimes. *Operations Research*, 70(6): 3261–3281.

Javanmard A, Nazerzadeh H (2019) Dynamic pricing in high-dimensions. *Journal of Machine Learning Research*, 20(1): 315–363.

Ke G, Meng Q, Finley T, Wang T, Chen W, Ma W, Ye Q, Liu T (2017) LightGBM: A highly efficient gradient boosting decision tree. *Advances in Neural Information Processing Systems*, 30: 3149–3157.

Keskin NB, Zeevi A (2014) Dynamic pricing with an unknown demand model: Asymptotically optimal semi-myopic policies. *Operations Research*, 62(5): 1142–1167.

Keskin NB, Li Y, Song JS (2022) Data-driven dynamic pricing and ordering with perishable inventory in a changing environment. *Management Science*, 68(3): 1938–1958.

Keskin NB, Li Y, Sunar N (2025) Data-driven clustering and feature-based retail electricity pricing with smart meters. *Operations Research*, 73(5): 2636–2660.

Lattimore T, Szepesvári C (2020) *Bandit Algorithms* (Cambridge University Press, Cambridge, UK).

Li M, Simchi-Levi D, Tan R, Wang C, Wu MX (2023) Contextual offline demand learning and pricing with separable models. *SSRN Electronic Journal*, https://doi.org/10.2139/ssrn.4619018.

Luo Y, Sun WW, Liu Y (2024) Distribution-free contextual dynamic pricing. *Mathematics of Operations Research*, 49(1): 599–618.

Miao S, Chen X, Chao X, Liu J, Zhang Y (2022) Context-based dynamic pricing with online clustering. *Production and Operations Management*, 31(9): 3559–3575.

Nambiar M, Simchi-Levi D, Wang H (2019) Dynamic learning and pricing with model misspecification. *Management Science*, 65(11): 4980–5000.

Perchet V, Rigollet P (2013) The multi-armed bandit problem with covariates. *The Annals of Statistics*, 41(2): 693–721.

Prokhorenkova L, Gusev G, Vorobev A, Dorogush AV, Gulin A (2018) CatBoost: Unbiased boosting with categorical features. *Advances in Neural Information Processing Systems*, 31: 6639–6649.

Qiang S, Bayati M (2016) Dynamic pricing with demand covariates. *arXiv preprint*, arXiv:1604.07469.

Rigollet P, Zeevi A (2010) Nonparametric bandits with covariates. *arXiv preprint*, arXiv:1003.1630.

Rusmevichientong P, Tsitsiklis JN (2010) Linearly parameterized bandits. *Mathematics of Operations Research*, 35(2): 395–411.

Slivkins A (2011) Contextual bandits with similarity information. *Proceedings of the 24th Annual Conference on Learning Theory (COLT)*, 679–702.

Wang Y, Chen B, Simchi-Levi D (2021a) Multimodal dynamic pricing. *Management Science*, 67(10): 6136–6152.

Wang Z, Deng S, Ye Y (2014) Close the gaps: A learning-while-doing algorithm for single-product revenue management problems. *Operations Research*, 62(2): 318–331.

Xu J, Wang YX (2021) Logarithmic regret in feature-based dynamic pricing. *Advances in Neural Information Processing Systems*, 34: 13898–13910.

论文信息

Bu, Jinzhi, David Simchi-Levi, and Chonghuan Wang. 2026. “Context-Based Dynamic Pricing with Separable Demand Models.” Management Science 72(8): 7279–7298.

原文 DOI ↗Google Scholar ↗