Представления вероятностными полиномами Бернулли и вырожденными полиномами Бернулли

Основная информация

• ID статьи: 2510.21558
• Название: Representations by probabilistic Bernoulli and degenerate Bernoulli polynomials
• Авторы: Dae San Kim (Sogang University), Taekyun Kim (Kwangwoon University)
• Классификация: math.NT (Теория чисел), math.PR (Теория вероятностей)
• Дата подачи: 24 октября 2025 г.
• Ссылка на статью: https://arxiv.org/abs/2510.21558v1

Аннотация

В данной работе исследуется проблема представления произвольных полиномов с использованием вероятностных полиномов Бернулли и вырожденных полиномов Бернулли, связанных со случайными величинами Y, для которых производящая функция моментов существует в окрестности начала координат. Кроме того, рассматривается проблема представления произвольных полиномов с помощью соответствующих полиномов высшего порядка. С использованием умбрального исчисления авторы разработали явные формулы для этих представлений и продемонстрировали результаты для нескольких дискретных и непрерывных случайных величин Y.

Исследовательский контекст и мотивация

1. Основная исследовательская проблема

Центральная проблема, которую решает данная работа: как представить произвольный полином в виде линейной комбинации вероятностных полиномов Бернулли и вероятностных вырожденных полиномов Бернулли с явными формулами для коэффициентов.

2. Значимость проблемы

• Теоретическое значение: Полиномы Бернулли и их варианты занимают фундаментальное место в теории чисел и комбинаторике; теория их представлений является важной частью теории специальных функций
• Историческое наследие: Исследование вырожденных специальных полиномов восходит к пионерской работе Carlitz (1979) о вырожденных полиномах Бернулли и Эйлера
• Современное развитие: Вероятностные расширения специальных полиномов в последнее время получили широкое распространение, объединяя теорию вероятностей и теорию специальных функций

3. Ограничения существующих методов

• Вероятностные числа Стирлинга, определённые через производящие функции кумулянтов в литературе 2 и 18, не обладают свойствами ортогональности и взаимности
• Отсутствие этих свойств затрудняет решение обратной задачи (восстановление полинома из коэффициентов представления)
• Существующие методы доказательства (например, доказательство тождества Miki) часто очень сложны и требуют глубоких инструментов, таких как p-адический анализ и квантовая теория поля

4. Исследовательская мотивация

• Установить теоретическую базу вероятностных чисел Стирлинга на основе ортогональности
• Разработать лаконичные формулы представления полиномов, избегая сложных методов доказательства
• Проверить эффективность и практичность теории на конкретных примерах

Основные вклады

1. Установлена полная теоретическая база представлений: даны явные формулы для представления произвольного полинома с использованием вероятностных полиномов Бернулли $B_k^Y(x)$ и вероятностных вырожденных полиномов Бернулли $\beta_{k,\lambda}^Y(x)$ (теоремы 3.1 и 3.3)
2. Расширение на случай высшего порядка: предоставлены формулы для представления произвольных полиномов с использованием полиномов Бернулли высшего порядка $B_k^{Y,(r)}(x)$ и $\beta_{k,\lambda}^{Y,(r)}(x)$ (теоремы 4.1 и 4.2)
3. Развита ключевая теория ортогональности: доказано, что вероятностные числа Стирлинга $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 типичных случайных величин (распределения Бернулли, биномиальное, Пуассона, геометрическое, экспоненциальное, гамма-распределение) даны явные представления $x^n$
5. Упрощены доказательства известных тождеств: с использованием формулы (3.22) даны простые доказательства тождества Miki и тождества FPZ, избегая сложных исходных методов

Подробное описание методов

Определение задачи

Входные данные: произвольный полином $p(x) \in \mathbb{C}[x]$ степени n

Выходные данные: коэффициенты представления $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. Определение вероятностных чисел Стирлинга

Вероятностные числа Стирлинга второго рода определяются через производящую функцию: $\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)$ — дельта-ряд (с $a_0=0, a_1=E[Y]\neq 0$).

Вероятностные числа Стирлинга первого рода определяются через композиционную обратную: $\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. Рамки умбрального исчисления

Работа использует умбральное исчисление для построения теории. Ключевые элементы включают:

• Последовательности Шеффера: $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)$
• Представление вероятностных полиномов Бернулли в виде последовательности Шеффера: $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: Представление вероятностными полиномами Бернулли (теорема 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) На основе производных и чисел Стирлинга: $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: Представление вероятностными вырожденными полиномами Бернулли (теорема 3.3)

Структура полностью аналогична алгоритму 1, требуется только заменить $S_1^Y$ на $S_{1,\lambda}^Y$, а $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$, теория вырождается в теорию представлений классических полиномов Бернулли, формулы упрощаются до: $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. Распределение Бернулли: $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. Распределение Пуассона: параметр $\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. Гамма-распределение: параметры $\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)$: вероятностные числа Стирлинга первого рода
• $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 путём интегрирования
• Коэффициенты включают комбинации чисел Бернулли $B_l$ и чисел Стирлинга $S_2(r,l)$

Пример 2: Распределение Бернулли ($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$ — числа Бернулли второго рода, определяемые через $\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: Распределение Пуассона ($\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$)

Это наиболее сложный пример. Требуется использование чисел Фробениуса-Эйлера $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 требует квотиентов Ферма или p-адического анализа), метод данной работы требует только вычисления интегралов и разностей
2. Унификация: все примеры следуют одной и той же вычислительной схеме, различаются только конкретные $f_Y(t)$ и $S_1^Y(n,k)$
3. Сложность вычисления:
   • Дискретные распределения обычно более лаконичны (например, Бернулли, Пуассон)
   • Непрерывные распределения могут включать более сложные интегралы (например, геометрическое распределение)
   • Экспоненциальное распределение наиболее лаконично
4. Дополнительная сложность вырожденной версии: представления вырожденными полиномами Бернулли обычно включают дополнительные суммирования по числам Стирлинга

Связанные работы

1. Исторический контекст развития

Классическая теория:

• Теория представлений полиномов Бернулли — основа теории специальных функций
• Классический результат, даваемый формулой (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: пионерское исследование вырожденных чисел Стирлинга, чисел Бернулли и Эйлера
• Недавние работы Kim и др. 13,16,19,20,23: систематическое развитие теории вырожденных специальных полиномов

Вероятностные расширения:

• Adell и др. 1,2,3: введение концепции вероятностных чисел Стирлинга
• Kim и др. 18,21,22: развитие теории вероятностных вырожденных полиномов

2. Отличие данной работы от связанных работ

Отличие от Adell-Bényi 2:

• 2 определяет $S_1^Y(n,k)$ через производящую функцию кумулянтов
• Данная работа определяет через композиционную обратную, обеспечивая ортогональность
• Ключевое преимущество: ортогональность делает обратную задачу разрешимой

Отличие от Kim-Kim 18:

• 18 рассматривает вырожденный случай, но не даёт общей теории представлений
• Данная работа унифицированно рассматривает невырожденный и вырожденный случаи

Сравнение с Kim-Kim 16:

• 16 даёт представления вырожденными полиномами Бернулли $\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: использование формул квотиентов Ферма по модулю $p^2$
• Shiratani-Yokoyama 30: p-адический анализ
• Gessel 12: две формы выражения чисел Стирлинга

Метод данной работы: прямое применение формулы (3.22), требующее только вычисления производных и интегралов

Заключение и обсуждение

Основные выводы

1. Полнота теории: установлена полная теория представлений вероятностными полиномами Бернулли и вырожденной версией, включая базовый и высший порядок
2. Вычислительная эффективность: предоставлены три эквивалентные формулы для вычисления коэффициентов, подходящие для различных вычислительных сценариев
3. Широкая применимость: теория применима к любой случайной величине, производящая функция моментов которой существует в окрестности начала координат
4. Упрощение доказательств: предоставлены более лаконичные пути доказательства известных тождеств

Ограничения

1. Ограничения условий:
   • Требуется $E[Y] \neq 0$
   • Производящая функция моментов должна существовать в окрестности начала координат
   • Исключаются некоторые важные распределения (например, распределение Коши)
2. Сложность вычисления:
   • Требуется предварительное вычисление $S_1^Y(n,k)$ и $S_{1,\lambda}^Y(n,k)$
   • Для сложных распределений (например, геометрического) формулы могут быть очень громоздкими
3. Численная устойчивость:
   • Включает высокие порядки разностей и числа Стирлинга, что может привести к проблемам численной устойчивости
   • Работа не обсуждает численную реализацию
4. Глубина теории:
   • Главным образом, это выведение комбинаторных тождеств
   • Отсутствует асимптотический анализ или исследование глубоких теоретико-числовых свойств

Направления будущих исследований

Работа не указывает явно направления будущих исследований, но можно предположить:

1. Расширение на другие специальные полиномы: например, полиномы Эйлера, полиномы Генокки и т.д.
2. Многомерное обобщение: исследование многомерных вероятностных полиномов Бернулли
3. Численные алгоритмы: разработка стабильных и эффективных методов численного вычисления
4. Исследование приложений: поиск приложений в теории чисел, комбинаторике, квантовой теории поля

Глубокая оценка

Преимущества

1. Инновационность методов (★★★★☆)

• Рамки ортогональности: через обеспечение ортогональности чисел Стирлинга решена ключевая проблема литературы 2,18
• Применение умбрального исчисления: систематическое применение теории умбрального исчисления делает доказательства лаконичными и элегантными
• Унифицированная теория: невырожденный, вырожденный и высший порядок объединены в единую рамку

2. Полнота теории (★★★★★)

• Четыре основные теоремы: охватывают все важные случаи (теоремы 3.1, 3.3, 4.1, 4.2)
• Два основных предложения: устанавливают ортогональность и взаимность (предложения 1.1, 1.2)
• Систематическая подготовка: раздел 1 подробно вводит необходимый фон

3. Достаточность экспериментов (★★★★☆)

• 6 случайных величин: охватывают типичные дискретные и непрерывные распределения
• Два типа представлений: каждый пример даёт невырожденную и вырожденную версию
• Детальные вычисления: показаны ключевые промежуточные шаги (например, формулы 5.5, 5.19-5.20)

Недостатки:

• Отсутствуют численные проверки
• Не сравнивается вычислительная эффективность различных формул

4. Ясность изложения (★★★★★)

• Чёткая структура: от подготовки → умбральное исчисление → основные результаты → примеры, логика строгая
• Стандартная нотация: единообразное использование верхнего индекса Y для обозначения связи со случайной величиной
• Достаточные детали: шаги доказательства подробны, удобны для понимания читателем

Недостатки

1. Теоретические ограничения

• Жёсткие условия: $E[Y] \neq 0$ исключает симметричные распределения (например, стандартное нормальное распределение)
• Отсутствие анализа ошибок: не обсуждается ошибка усечения или численная точность

2. Недостатки экспериментов

• Отсутствие численной реализации: все результаты в символической форме, не даны численные примеры
• Отсутствие сравнения производительности: какая из трёх форм формул вычисляется быстрее?
• Отсутствие визуализации: не показаны графики полиномов или коэффициентов

3. Ограниченные приложения

• Теоретическая ориентация: главным образом математические выведения, отсутствуют сценарии практического применения
• Слабая связь с теорией вероятностей: хотя введены случайные величины, не глубоко исследуется вероятностный смысл

4. Обзор литературы

• Разрозненные связанные работы: распределены в введении и разделе 5, недостаточно сконцентрированы
• Недостаточное сравнение: техническое сравнение с 2,18 недостаточно подробно

Оценка влияния

1. Вклад в область (★★★★☆)

• Заполнение теоретического пробела: решена проблема ортогональности вероятностных чисел Стирлинга
• Методологический вклад: показана мощь умбрального исчисления в вероятностных расширениях
• Связь различных областей: объединены теория вероятностей, комбинаторика, теория специальных функций

Потенциальное влияние:

• Может стать стандартной ссылкой в теории вероятностных специальных функций
• Может вдохновить исследования вероятностных расширений других специальных полиномов

2. Практическая ценность (★★★☆☆)

• Символические вычисления: применимо к системам компьютерной алгебры (например, Mathematica, Maple)
• Теоретический инструмент: предоставляет новый инструмент для доказательства комбинаторных тождеств
• Образовательная ценность: подходит как дополнительный материал для курсов специальных функций

Ограничения:

• Сценарии прямого применения неясны
• Требуется дальнейшее развитие для применения к практическим задачам

3. Воспроизводимость (★★★☆☆)

• Явные формулы: все формулы имеют чёткие определения
• Зависимость от внешних результатов: ключевое вычисление $S_1^Y(n,k)$ зависит от литературы 14
• Отсутствие кода: не предоставлен код реализации

Рекомендации:

• Предоставить реализацию на Mathematica или Python
• Создать онлайн-калькулятор

Применимые сценарии

1. Теоретические исследования

• Доказательство комбинаторных тождеств: упрощение доказательств сложных тождеств
• Теория специальных функций: расширение теории полиномов Бернулли
• Теория чисел: возможное применение к исследованию свойств сравнимости чисел Бернулли

2. Символические вычисления

• Разложение полиномов: разложение произвольных полиномов в базис специальных функций
• Вычисление интегралов: упрощение интегралов с использованием свойств полиномов Бернулли

3. Образование

• Курсы специальных функций: демонстрация современных методов исследования
• Комбинаторика: продвинутые приложения чисел Стирлинга
• Умбральное исчисление: конкретные примеры приложений

4. Потенциальные приложения

• Квантовая теория поля: применение чисел Бернулли в вычислении диаграмм Фейнмана
• Теория Громова-Виттена: связь с тождеством FPZ
• Асимптотический анализ: возможное применение к асимптотическому разложению некоторых сумм

Комплексная оценка

Ключевые ссылки

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

• Введение концепции вероятностных чисел Стирлинга, но определение не обладает ортогональностью

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 серия

• Стандартная ссылка на умбральное исчисление

Резюме: Это высококачественная теоретическая математическая работа, вносящая существенный вклад в теорию вероятностных специальных функций. Через установление рамок ортогональности автор решает ключевые проблемы существующей литературы и разрабатывает полную теорию представлений. Основная ценность работы заключается в систематичности теории и элегантности методов, обеспечивая прочную основу для последующих исследований. Основные направления улучшения — добавление численных экспериментов и исследование практических приложений.