Representations by probabilistic Bernoulli and degenerate Bernoulli polynomials

基本信息

• 论文ID: 2510.21558
• 标题: Representations by probabilistic Bernoulli and degenerate Bernoulli polynomials
• 作者: Dae San Kim (Sogang University), Taekyun Kim (Kwangwoon University)
• 分类: math.NT (Number Theory), math.PR (Probability)
• 提交时间: 2025年10月24日
• 论文链接: https://arxiv.org/abs/2510.21558v1

摘要

本文研究任意多项式使用概率Bernoulli多项式和退化Bernoulli多项式的表示问题，这些多项式与矩母函数在原点邻域存在的随机变量Y相关联。此外，本文探讨了用高阶对应多项式表示任意多项式的问题。借助umbral演算，作者发展了这些表示的显式公式，并针对若干离散和连续随机变量Y展示了结果。

研究背景与动机

1. 研究问题

本文要解决的核心问题是：如何将任意多项式表示为概率Bernoulli多项式和概率退化Bernoulli多项式的线性组合，并给出显式的系数公式。

2. 问题重要性

• 理论意义：Bernoulli多项式及其变体在数论、组合数学中具有基础性地位，其表示理论是特殊函数理论的重要组成部分
• 历史传承：退化特殊多项式的研究源于Carlitz在1979年关于退化Bernoulli和Euler多项式的开创性工作
• 现代发展：概率扩展的特殊多项式近年来受到广泛研究，将概率论与特殊函数理论结合

3. 现有方法的局限性

• 文献2和18中基于累积量生成函数定义的概率Stirling数不具备正交性和逆关系性质
• 缺乏这些性质使得逆问题（从表示系数反推多项式）难以解决
• 现有证明方法（如Miki恒等式的证明）往往非常复杂，涉及p-adic分析、量子场论等高深工具

4. 研究动机

• 建立基于正交性的概率Stirling数理论框架
• 发展简洁的多项式表示公式，避免复杂的证明技术
• 通过具体例子验证理论的有效性和实用性

核心贡献

1. 建立了完整的表示理论框架：给出任意多项式用概率Bernoulli多项式$B_k^Y(x)$和概率退化Bernoulli多项式$\beta_{k,\lambda}^Y(x)$表示的显式公式（定理3.1和3.3）
2. 扩展到高阶情形：提供了用高阶概率Bernoulli多项式$B_k^{Y,(r)}(x)$和$\beta_{k,\lambda}^{Y,(r)}(x)$表示任意多项式的公式（定理4.1和4.2）
3. 发展了关键的正交性理论：证明了概率Stirling数$S_1^Y(n,k)$和$S_2^Y(n,k)$，以及退化版本$S_{1,\lambda}^Y(n,k)$和$S_{2,\lambda}^Y(n,k)$满足正交关系和逆关系（命题1.1和1.2）
4. 提供了丰富的具体例子：针对6种常见随机变量（Bernoulli、二项、Poisson、几何、指数、Gamma分布）给出了$x^n$的显式表示
5. 简化了已知恒等式的证明：利用公式(3.22)给出了Miki恒等式和FPZ恒等式的简单证明，避免了原有的复杂技术

方法详解

任务定义

输入：度数为n的任意多项式$p(x) \in \mathbb{C}[x]$

输出：表示系数$a_0, a_1, \ldots, a_n$，使得 $p(x) = \sum_{k=0}^n a_k B_k^Y(x) \quad \text{或} \quad p(x) = \sum_{k=0}^n a_k \beta_{k,\lambda}^Y(x)$

约束条件：随机变量Y的矩母函数$E[e^{Yt}]$在原点邻域存在，且$E[Y] \neq 0$

理论基础

1. 概率Stirling数的定义

第二类概率Stirling数通过生成函数定义： $\frac{1}{k!}(E[e^{Yt}] - 1)^k = \sum_{n=k}^\infty S_2^Y(n,k) \frac{t^n}{n!}$

引入记号$e_Y(t) = E[e^{Yt}] - 1$，则$e_Y(t)$是delta级数（$a_0=0, a_1=E[Y]\neq 0$）。

第一类概率Stirling数通过复合逆定义： $\frac{1}{k!}(\bar{e}_Y(t))^k = \sum_{n=k}^\infty S_1^Y(n,k) \frac{t^n}{n!}$ 其中$\bar{e}_Y(t)$是$e_Y(t)$的复合逆，满足$e_Y(\bar{e}_Y(t)) = \bar{e}_Y(e_Y(t)) = t$。

2. 正交性和逆关系（命题1.1）

$\sum_{k=l}^n S_2^Y(n,k) S_1^Y(k,l) = \delta_{n,l}$

这个正交性导出重要的逆关系： $a_n = \sum_{k=0}^n S_2^Y(n,k) b_k \Leftrightarrow b_n = \sum_{k=0}^n S_1^Y(n,k) a_k$

3. Umbral演算框架

论文使用umbral演算建立理论。关键要素包括：

• Sheffer序列：$s_n(x) \sim (g(t), f(t))$当且仅当 $\frac{1}{g(\bar{f}(t))} e^{x\bar{f}(t)} = \sum_{k=0}^\infty s_k(x) \frac{t^k}{k!}$
• 微分算子性质：$f(t)s_n(x) = ns_{n-1}(x)$
• 概率Bernoulli多项式的Sheffer表示： $B_n^Y(x) \sim \left(g(t) = \frac{e^t-1}{f(t)}, f(t)\right)$ 其中$\bar{f}(t) = \log E[e^{Yt}]$

核心算法流程

算法1：概率Bernoulli多项式表示（定理3.1）

步骤1：计算$a_0$$a_0 = \int_0^1 \frac{t}{f(t)} p(x) dx$

步骤2：构造辅助函数$a(x) = p(x+1) - p(x) = \Delta p(x)$

步骤3：计算$a_{r+1}$（$r=0,1,\ldots,n-1$）

三种等价形式：

(a) 基于差分算子： $a_{r+1} = \frac{1}{r+1} \sum_{j=r}^{n-1} S_1^Y(j,r) \frac{1}{j!} \Delta^{j+1} p(0)$

(b) 基于导数和Stirling数： $a_{r+1} = \frac{1}{r+1} \sum_{k=r}^{n-1} \sum_{j=r}^k \frac{1}{k!} S_2(k,j) S_1^Y(j,r) \Delta p^{(k)}(0)$

(c) 基于直接展开： $a_{r+1} = \frac{1}{r+1} \sum_{j=r}^{n-1} \sum_{k=0}^{j+1} (-1)^{j+1-k} \frac{1}{j!} \binom{j+1}{k} S_1^Y(j,r) p(k)$

算法2：概率退化Bernoulli多项式表示（定理3.3）

结构与算法1完全类似，只需将$S_1^Y$替换为$S_{1,\lambda}^Y$，$f(t)$替换为$f(t)$（其复合逆为$\bar{f}(t) = \log E[e_\lambda^Y(t)]$）。

技术创新点

1. 正交性的关键作用

通过建立正交关系，作者巧妙地将表示问题转化为求解线性系统。具体地：

• 从$p(x) = \sum_{k=0}^n a_k B_k^Y(x)$出发
• 计算差分$\Delta p(x) = \sum_{k=1}^n k a_k \sum_{j=0}^{k-1} S_2^Y(k-1,j) (x)_j$
• 利用正交性反解出$a_k$

2. 与经典情形的统一

当$Y=1$时，理论退化为经典Bernoulli多项式的表示理论，公式简化为： $a_k = \frac{1}{k!} \int_0^1 p^{(k)}(x) dx$

3. 高阶表示的处理

对于高阶情形（定理4.1和4.2），区分$r>n$和$r\leq n$两种情况：

• $r>n$时：所有系数都涉及积分算子$I^{r-k}$
• $r\leq n$时：前$r$项涉及积分算子，后续项涉及差分算子

4. 计算技巧

通过引入线性算子$I$（积分算子）和差分算子$\Delta$的代数性质，将复杂的表达式转化为可计算的形式。

实验设置

数据集（随机变量选择）

论文选择了6种具有代表性的随机变量进行验证：

离散随机变量：

1. Bernoulli分布：$p(0)=1-p, p(1)=p$（$0<p\leq 1$）
2. 二项分布：参数$(m,p)$，$p(i) = \binom{m}{i} p^i (1-p)^{m-i}$
3. Poisson分布：参数$\alpha>0$，$p(i) = e^{-\alpha} \frac{\alpha^i}{i!}$
4. 几何分布：参数$0<p<1$，$p(i) = (1-p)^{i-1} p$

连续随机变量： 5. 指数分布：参数$\alpha>0$，$f(y) = \alpha e^{-\alpha y}$（$y\geq 0$） 6. Gamma分布：参数$\alpha,\beta>0$，$f(y) = \frac{\beta e^{-\beta y} (\beta y)^{\alpha-1}}{\Gamma(\alpha)}$

计算要素

对每个随机变量Y，需要计算（引用文献14的结果）：

• $f_Y(t)$：$\bar{f}_Y(t) = \log E[e^{Yt}]$的复合逆
• $f_{Y,\lambda}(t)$：$\bar{f}_{Y,\lambda}(t) = \log E[e_\lambda^Y(t)]$的复合逆
• $S_1^Y(n,k)$：第一类概率Stirling数
• $S_{1,\lambda}^Y(n,k)$：退化版本

验证策略

对每个随机变量，计算$x^n$的表示： $x^n = \sum_{k=0}^n a_k B_k^Y(x) \quad \text{和} \quad x^n = \sum_{k=0}^n a_k \beta_{k,\lambda}^Y(x)$

实验结果

主要结果展示

例1：指数分布（$\alpha>0$）

这是最简洁的例子。从文献14： $f_Y(t) = \alpha(1-e^{-t}), \quad S_1^Y(n,k) = (-1)^{n-k} \binom{n}{k} (n-1)^{n-k} \alpha^k$

结果： $x^n = \frac{1}{\alpha} B_0^Y(x) + \sum_{k=1}^n \left\{ \frac{1}{k} \sum_{j=k-1}^{n-1} (-1)^{j-k+1} \binom{j}{k-1} (j-1)^{j-k+1} \alpha^{k-1} \frac{1}{j!} \Delta^{j+1} 0^n \right\} B_k^Y(x)$

对于退化版本： $x^n = \frac{1}{\alpha} \sum_{r=0}^n \sum_{l=0}^r \binom{n}{r} S_2(r,l) (-1)^{l-r} (\alpha\lambda)^l B_l \beta_{0,\lambda}^Y(x) + \cdots$

分析：

• $a_0 = \frac{1}{\alpha}$非常简洁，通过引理5.1计算积分得到
• 系数涉及Bernoulli数$B_l$和Stirling数$S_2(r,l)$的组合

例2：Bernoulli分布（$0<p\leq 1$）

从文献14： $f_Y(t) = \log\left(1 + \frac{1}{p}(e^t-1)\right), \quad S_1^Y(n,k) = \frac{1}{p^n} S_1(n,k)$

结果： $x^n = \sum_{l=0}^n \frac{1}{p^{l-1}} S_2(n,l) b_l B_0^Y(x) + \sum_{k=1}^n \left\{ \frac{1}{k} \sum_{j=k-1}^{n-1} \frac{1}{p^j} S_1(j,k-1) \frac{1}{j!} \Delta^{j+1} 0^n \right\} B_k^Y(x)$

其中$b_l$是第二类Bernoulli数，由$\frac{t}{\log(1+t)} = \sum_{l=0}^\infty b_l \frac{t^l}{l!}$定义。

关键计算（公式5.5）： $\frac{t}{f_Y(t)} x^n = \sum_{r=0}^n \sum_{l=0}^r \binom{n}{r} \frac{1}{p^{l-1}} S_2(r,l) b_l B_{n-r}(x)$

例3：Poisson分布（$\alpha>0$）

$f_Y(t) = \log\left(1 + \frac{t}{\alpha}\right), \quad S_1^Y(n,k) = \sum_{l=k}^n \frac{1}{\alpha^l} S_1(l,k) S_1(n,l)$

结果： $x^n = \sum_{l=0}^n \binom{n}{l} \frac{1}{n-l+1} \frac{1}{\alpha^{l-1}} b_l B_0^Y(x) + \sum_{k=1}^n \left\{ \frac{1}{k} \sum_{j=k-1}^{n-1} \sum_{l=k-1}^j \frac{1}{\alpha^l} S_1(l,k-1) S_1(j,l) \frac{1}{j!} \Delta^{j+1} 0^n \right\} B_k^Y(x)$

例4：几何分布（$0<p<1$）

这是最复杂的例子。需要用到Frobenius-Euler数$H_j^{(r)}(u)$： $\left(\frac{1-u}{e^t-u}\right)^r = \sum_{n=0}^\infty H_n^{(r)}(u) \frac{t^n}{n!}$

$a_0$的计算（公式5.23）： $a_0 = \frac{1}{p} \sum_{j=0}^n \sum_{l=0}^\infty \sum_{r=0}^l (-1)^r \frac{1}{l!} \binom{n}{j} \binom{l}{r} \left(\frac{p}{1-p}\right)^l b_l H_j^{(r)}\left(\frac{p}{p-1}\right) (1-p(1-\delta_{n,j}))$

案例分析：引理5.1的作用

引理5.1： $\int_0^1 B_n(x) dx = \delta_{n,0}, \quad \int_0^1 B_n(-x) dx = (-1)^n$

这个引理在所有例子的$a_0$计算中起关键作用。例如在指数分布情形： $a_0 = \int_0^1 \frac{t}{\alpha(1-e^{-t})} x^n dx = \frac{1}{\alpha} (-1)^n \int_0^1 B_n(-x) dx = \frac{1}{\alpha}$

实验发现

1. 简化效果显著：相比文献中的复杂证明（如Miki恒等式需要Fermat商或p-adic分析），本文方法只需计算积分和差分
2. 统一性：所有例子都遵循相同的计算框架，只是具体的$f_Y(t)$和$S_1^Y(n,k)$不同
3. 计算复杂度：
   • 离散分布通常更简洁（如Bernoulli、Poisson）
   • 连续分布可能涉及更复杂的积分（如几何分布）
   • 指数分布最为简洁
4. 退化版本的额外复杂性：退化Bernoulli多项式的表示通常涉及额外的Stirling数求和

相关工作

1. 历史发展脉络

经典理论：

• Bernoulli多项式的表示理论是特殊函数论的基础内容
• 公式(3.22)给出的经典结果：$p(x) = \sum_{k=0}^n a_k B_k(x)$，其中$a_k = \frac{1}{k!} \int_0^1 p^{(k)}(x) dx$

退化理论：

• Carlitz (1979) 4：开创性地研究退化Stirling数、Bernoulli数和Euler数
• Kim等人近期工作13,16,19,20,23：系统发展退化特殊多项式理论

概率扩展：

• Adell等人1,2,3：引入概率Stirling数的概念
• Kim等人18,21,22：发展概率退化多项式理论

2. 本文与相关工作的区别

与Adell-Bényi 2的区别：

• 2基于累积量生成函数定义$S_1^Y(n,k)$
• 本文基于复合逆定义，确保正交性
• 关键优势：正交性使得逆问题可解

与Kim-Kim 18的区别：

• 18处理退化情形但未给出一般表示理论
• 本文统一处理非退化和退化情形

与Kim-Kim 16的比较：

• 16给出退化Bernoulli多项式$\beta_{k,\lambda}(x)$（$Y=1$）的表示
• 本文扩展到一般随机变量Y

3. 应用实例的创新

Miki恒等式（公式1.1）： $\sum_{k=1}^{n-1} \frac{B_k(x) B_{n-k}(x)}{k(n-k)} = \frac{2}{n} \sum_{k=0}^{n-2} \frac{1}{n-k} \binom{n}{k} B_{n-k} B_k(x) + \frac{2}{n} H_{n-1} B_n(x)$

传统证明方法：

• Miki 24：使用Fermat商模$p^2$的公式
• Shiratani-Yokoyama 30：p-adic分析
• Gessel 12：Stirling数的两种表达式

本文方法：直接应用公式(3.22)，只需计算导数和积分

结论与讨论

主要结论

1. 理论完整性：建立了概率Bernoulli多项式和退化版本的完整表示理论，包括基础版本和高阶版本
2. 计算有效性：提供了三种等价的系数计算公式，适用于不同的计算场景
3. 广泛适用性：理论适用于矩母函数在原点邻域存在的任意随机变量
4. 简化证明：为已知恒等式提供了更简洁的证明途径

局限性

1. 条件限制：
   • 要求$E[Y] \neq 0$
   • 矩母函数必须在原点邻域存在
   • 排除了一些重要分布（如Cauchy分布）
2. 计算复杂度：
   • 需要预先计算$S_1^Y(n,k)$和$S_{1,\lambda}^Y(n,k)$
   • 对于复杂分布（如几何分布），公式可能非常复杂
3. 数值稳定性：
   • 涉及高阶差分和Stirling数，可能存在数值稳定性问题
   • 论文未讨论数值实现
4. 理论深度：
   • 主要是组合恒等式的推导
   • 缺乏渐近分析或深层数论性质的探讨

未来方向

论文未明确提出未来方向，但可推测：

1. 扩展到其他特殊多项式：如Euler多项式、Genocchi多项式等
2. 多变量推广：研究多元概率Bernoulli多项式
3. 数值算法：发展稳定高效的数值计算方法
4. 应用探索：在数论、组合学、量子场论中寻找应用

深度评价

优点

1. 方法创新性（★★★★☆）

• 正交性框架：通过确保Stirling数的正交性，解决了文献2,18中的关键缺陷
• Umbral演算应用：系统运用umbral演算理论，使证明简洁优雅
• 统一理论：将非退化、退化、高阶情形纳入统一框架

2. 理论完整性（★★★★★）

• 四个主要定理：覆盖所有重要情形（定理3.1, 3.3, 4.1, 4.2）
• 两个基础命题：建立正交性和逆关系（命题1.1, 1.2）
• 系统的预备知识：第1节详尽介绍所需背景

3. 实验充分性（★★★★☆）

• 6种随机变量：涵盖常见的离散和连续分布
• 两种表示：每个例子都给出非退化和退化版本
• 详细计算：展示了关键中间步骤（如公式5.5, 5.19-5.20）

不足：

• 缺乏数值验证
• 未比较不同公式的计算效率

4. 写作清晰度（★★★★★）

• 结构清晰：从预备知识→umbral演算→主要结果→例子，逻辑严密
• 记号规范：统一使用上标Y表示与随机变量的关联
• 细节充分：证明步骤详细，便于读者理解

不足

1. 理论局限

• 条件苛刻：$E[Y] \neq 0$排除了对称分布（如标准正态分布）
• 缺乏误差分析：未讨论截断误差或数值精度

2. 实验不足

• 无数值实现：所有结果都是符号形式，未给出数值例子
• 无性能比较：三种公式形式哪种计算最快？
• 无可视化：未展示多项式或系数的图形

3. 应用有限

• 理论导向：主要是数学推导，缺乏实际应用场景
• 与概率论的联系薄弱：虽然引入随机变量，但未深入探讨概率意义

4. 文献综述

• 相关工作分散：分布在引言和第5节，不够集中
• 对比不充分：与2,18的技术对比不够详细

影响力评估

1. 对领域的贡献（★★★★☆）

• 填补理论空白：解决了概率Stirling数的正交性问题
• 方法论贡献：展示了umbral演算在概率扩展中的威力
• 连接不同领域：将概率论、组合学、特殊函数理论结合

潜在影响：

• 可能成为概率特殊函数理论的标准参考
• 启发其他特殊多项式的概率扩展研究

2. 实用价值（★★★☆☆）

• 符号计算：可用于计算机代数系统（如Mathematica、Maple）
• 理论工具：为证明组合恒等式提供新工具
• 教学价值：适合作为特殊函数课程的补充材料

局限：

• 直接应用场景不明确
• 需要进一步发展才能用于实际问题

3. 可复现性（★★★☆☆）

• 公式明确：所有公式都有清晰定义
• 依赖外部结果：关键的$S_1^Y(n,k)$计算依赖文献14
• 无代码：未提供实现代码

建议：

• 提供Mathematica或Python实现
• 建立在线计算器

适用场景

1. 理论研究

• 组合恒等式证明：简化复杂恒等式的证明
• 特殊函数理论：扩展Bernoulli多项式理论
• 数论：可能用于Bernoulli数的同余性质研究

2. 符号计算

• 多项式展开：在特殊函数基下展开任意多项式
• 积分计算：利用Bernoulli多项式的性质简化积分

3. 教学

• 特殊函数课程：展示现代研究方法
• 组合数学：Stirling数的高级应用
• Umbral演算：具体应用实例

4. 潜在应用

• 量子场论：Bernoulli数在Feynman图计算中的应用
• Gromov-Witten理论：与FPZ恒等式的联系
• 渐近分析：可能用于某些和式的渐近展开

综合评分

参考文献（关键文献）

2 J. A. Adell, B. Bényi, Probabilistic Stirling numbers and applications, Aequat. Math. 98 (2024), 1627-1646.

• 提出概率Stirling数概念，但定义不具备正交性

4 L. Carlitz, Degenerate Stirling, Bernoulli and Eulerian numbers, Utilitas Math. 15 (1979), 51-88.

• 退化特殊数的开创性工作

14 D. S. Kim, T. Kim, Probabilisitc Stirling and degenerate Stirling numbers, Preprint.

• 提供本文所需的$S_1^Y(n,k)$计算结果

16 D. S. Kim, T. Kim, Representing polynomials by degenerate Bernoulli polynomials, Quaest. Math. 46 (2022), no. 5, 959-980.

• $Y=1$情形的前期工作

27-28 S. Roman, The umbral calculus系列

• Umbral演算的标准参考

总结：这是一篇高质量的理论数学论文，在概率特殊函数理论方面做出了实质性贡献。通过建立正交性框架，作者解决了现有文献的关键缺陷，并发展了完整的表示理论。论文的主要价值在于理论的系统性和方法的优雅性，为后续研究奠定了坚实基础。主要改进空间在于增加数值实验和探索实际应用。