УДК 519.114
ИНТЕРПРЕТАЦИИ ЭЛЕМЕНТОВ ДВУХ ИНТЕРПОЛЯЦИОННЫХ КЛАССОВ ЦЕЛОЧИСЛЕННЫХ ПОСЛЕДОВАТЕЛЬНОСТЕй
Бондаренко Леонид Николаевич1,
канд. техн. наук, доцент, e-mail: [email protected],
1Московский университет им. С.Ю. Витте, филиал в г. Сергиевом Посаде, г. Сергиев Посад, Россия
По разработанной автором методике в статье вводятся и исследуются два однопараметрических класса целочисленных последовательностей. Первый класс обеспечивает интерполяцию между последовательностью степеней двойки и последовательностью чисел Каталана, а второй - между последовательностью степеней двойки и последовательностью чисел Белла. Для получения разнообразных интерпретаций чисел рассматриваемых классов используются T-модели и коды Лемера перестановок. Эти интерпретации базируются на рекурсивном построении последовательностей числовых таблиц специального вида, определяющих T-модели, и свойствах кодов Лемера перестановок. Используемая методика приводит к простым алгоритмам построения двух классов множеств перестановок, отвечающих введенным классам числовых последовательностей. На полученных классах множеств перестановок удается также задать распределения вероятностей. Представление элементов второго класса числовых последовательностей с помощью чисел Стирлинга второго рода позволяет сопоставить соответствующему классу множеств перестановок класс упорядоченных разбиений множеств на определенное число блоков. Для чисел изучаемых классов последовательностей получены соотношения для их q-аналогов.
Ключевые слова: интерполяция последовательностей, T-модели, коды лемера, числа каталана, Белла и стирлинга второго рода, q-аналоги
INTERPRETATIONS OF ELEMENTS OF TWO INTERPOLATION CLASSES OF INTEGER SEQUENCES
Bondarenko L.N.1,
candidate of technical sciences, Associate Professor, e-mail: [email protected], 'Moscow Witte University, branch in Sergiev Posad, Sergiev Posad, Russia
According to the methodology developed by the author, two one-parameter classes of integer sequences are introduced and investigated in the paper. The first class provides an interpolation between a sequence ofpowers of two and a sequence of Catalan numbers, and the second class provides an interpolation between a sequence ofpowers of two and a sequence of Bell numbers. In order to obtain various interpretations of the numbers of the considered classes, T-models and permutations by Lehmer codes are used. These interpretations are based on the recursive construction ofsequences ofnumerical tables of a special kind defining T-models and the properties of permutations by Lehmer codes. The method used leads to simple algorithms for constructing two classes of permutation sets corresponding to the introduced classes of numerical sequences. On the basis of the obtained classes of permutation sets, it is also possible to set probability distributions. The representation of elements of the second class of numerical sequences using Stirling numbers of the second kind allows us to match the class of ordered partitions of sets into a certain number of blocks to the corresponding class ofpermutation sets. For the numbers of the studied classes of sequences, the relations for their q-analogues are obtained. Keywords: sequence interpolation, T-models, Lehmer codes, Catalan, Bell and Stirling numbers of the second kind, q-analogues
DOI 10.21777/2500-2112-2021-4-88-96
Введение
Методика построения однопараметрических классов возрастающих целочисленных последовательностей, осуществляющих интерполяцию между двумя заданными фиксированными последовательностями, и отвечающих им однопараметрических классов множеств перестановок степени п е N = {1,2,...} была разработана в статье [2]. Эта методика базируется на применении Г-моделей, определяемых рекурсивно и представляющих собой последовательности таблиц Т1, Т2, ... специального вида с элементами из N а также свойств кодов Лемера перестановок, позволяющих устанавливать частичный порядок на рассматриваемых множествах перестановок [1, 2]. По этой методике в статье [2] были построены классы числовых последовательностей и аналоги отвечающих им классов множеств перестановок, рассматриваемых в работах [5, 7]. Это позволило получить простые ^-аналоги соответствующих чисел, а также существенно упростить алгоритмы генерации классов используемых множеств перестановок.
Данную статью можно рассматривать как продолжение работы [2]. В ней с помощью Г-моделей строятся два однопараметрических класса возрастающих числовых последовательностей.
Первый класс осуществляет интерполяцию между последовательностью степеней 2п-1 и после-
1 (2п\
довательностью чисел Каталана Сп =- при п е N.
п +1 ^ п )
В этой записи чисел каталана используется биномиальный коэффициент, обозначаемый стандартным образом, а многие свойства чисел каталана и их многочисленные интерпретации имеются в [3, 11]. Для получения дополнительных сведений о числах Каталана можно также обратиться к обширной статье с номером А000108 из [10].
Этот класс дополняет класс, реализующий интерполяцию между последовательностью чисел Каталана Сп и последовательностью факториалов п! при п е N который рассматривался в [2, 5].
Второй класс осуществляет интерполяцию между последовательностью степеней 2п-1 и последовательностью чисел Белла Вп при п е N задаваемых как число способов всех разбиений множества {1,2,..., п} на блоки.
п
Число Белла Вп также можно определить как сумму всех чисел Стирлинга второго рода < .
описывающих число способов разбиений множества {1,2,..., п} ровно на k блоков и имеющих, как и биномиальные коэффициенты, многочисленные свойства [3]. Дополнительные сведения о числах Белла
В = $1 и
имеются в статье с номером А000110 из [10].
Этот класс дополняет класс, реализующий интерполяцию между последовательностью чисел Белла Вп и последовательностью факториалов п! при п е N который рассматривался в [2, 7].
Рассмотрение классов числовых последовательностей позволяет стандартизировать задание различных комбинаторных последовательностей, возникающих при решении соответствующих задач в информатике, статистике, химии, экономике и т.п. [1].
С помощью Г-моделей, таблицы которых представляют одну из многочисленных интерпретаций элементов изучаемых числовых последовательностей, часто удается за счет рекурсивности, заложенной в определении Г-модели, просто получать другие важные интерпретации этих элементов. В частности, такие важные комбинаторные интерпретации как множества перестановок и разбиений находят соответствующие применения в различных прикладных областях.
Для унификации обозначений для используемых в статье Г-моделей соответствующие таблицы, многочлены и т.п. имеют в двух разделах статьи одинаковые обозначения, а их определение осуществляется надлежащим образом в каждом разделе.
1. Интерполяция числовых последовательностей между 2я"1 и Сп
Для исследования класса целочисленных последовательностей, осуществляющего интерполяцию между последовательностью степеней 2п-1 и последовательностью чисел Каталана Сп при п е N5 несложно построить соответствующий однопараметрический класс Г-моделей на базе методики, разработанной в [2]. Этот класс задается тройкой (Б(т), 0т, Т(т)), в которой £(т) = {2,..., т + 1} - алфавит используемых символов, 0т : 5 ^ 2 3 ... 5 (шт(т,5) + 1) - отображение символа 5 е Б<т) в строку символов 0т (5), а Т(т) = (2) - начальная таблица. При фиксированном значении параметра
т е N каждая таблица Г-модели определяется рекурсивно соотношением Т(+1) = 0т(ТП )) при п е N.
Для преобразования последовательности таблиц Г-модели в числовую последовательность используем производящие многочлены ипт)^) = I еТ<т) , которые при значении t = 1 определяют искомые числа.
В рассматриваемом случае рекурсивность, заложенная в определение Г-модели, позволяет легко найти следующие рекуррентные соотношения:
и(т)^) = 1, (1 - t) ипт'^) + t2 и^«) = ипт1)(1), 2 < п < т;
Г т 1+1 , ч
Г 21 (т - к + 2]
unm)(t) = I (-1)к-1 Г к 2 ), п > т. (1)
к=1 V к
В формуле (1) [а] - целая часть числа а, а последовательности чисел, вычисляемые по этой формуле при т = 1,2,3,4,5,... и t = 1, принимают вид:
1, 2, 4, 8, 16, 32, 64, 128, 256, 512, 1024, 2048, 4096, ... [10, А000079];
1, 2, 5, 13, 34, 89, 233, 610, 1597, 4181, 10946, 28657, ... [10, А001519];
1, 2, 5, 14, 41, 122, 365, 1094, 3281, 9842, 29525, 88574, ... [10, А007051]; (2)
1, 2, 5, 14, 42, 131, 417, 1341, 4334, 14041, 45542, 147798, ... [10, А080937];
1, 2, 5, 14, 42, 132, 428, 1416, 4744, 16016, 54320, 184736, ... [10, А024175];
.............................................,
причем для значений т = 1,2,..., 8 дополнительные сведения об этих последовательностях можно также найти в соответствующих статьях [10].
В частности, первая строка в (2) определяет при п е N последовательность чисел 2п-1 [10, А000079], вторая строка - последовательность чисел Фибоначчи с нечетными номерами [10, А001519].
Неявное задание рекуррентным соотношением (1) многочленов ипт)^) при 2 < п < т существенно затрудняет вычисление с его помощью чисел ипт)(1).
Так как эти многочлены не зависят от параметра т е N то вычисление чисел ипт)(1) можно значительно упростить, используя предельный переход при т ^ да и обозначая Сп = и^ (1).
В этом случае для производящей функции ¥ (и) = I Сп ип при условии С0 = 1 на основе результатов из [2] легко получить уравнение
и¥ (и)2 - ¥ (и) + 1 = 0. (3)
Решение уравнения (3) приводит к известной производящей функции для чисел Каталана [3]
I Сип = 1 ,
п 2и
п=0 2и
а также из уравнения (3) следует рекуррентное соотношение [3]
Со = С = 1, Сп+1 = $С,С^,, п е N (4)
I=о
позволяющее просто вычислять числа Сп = и((п)(Х) в (1) при 2 < п < т.
Таким образом, при п е N класс последовательностей (2) осуществляет интерполяцию между последовательностью 2п-1 и последовательностью Сп.
При фиксированном значении параметра т е N определим аналогично статьям [1, 2] множество номеров таблицы Т^ рассматриваемой Г-модели рекурсивным способом.
Определение 1. Для начальной таблицы Т[т) = (2) с одним элементом 5 = 2 полагаем его номер V = 1, а Дт) = {1}; если V = v1... vn е Х(П° номер элемента 5 е Т[п), то номер 5' е Т(П) полагаем равным ^п+1 е где Vn+1 порядковый номер элемента 5' в строке вп (5) таблицы Тп+1.
Определение 1 дает возможность сопоставить каждому числу Ц;т) (1) рассматриваемого класса целочисленных последовательностей множество перестановок Р^ символов {1,2,..., п} такое, что |р (т)| = |д(т)| = и (т)(1)
Действительно, из этого определения следует, что номер V = v1 ...vn е Х(П) является кодом Леме-ра некоторой перестановки п = п1 ... пп е Р^, так как для этой перестановки при / = 1,2,..., п имеем н{ = {у : п. < п.,0 < у < , - 1,по = 0}.
Также в [2] доказано, что нахождение п е Р^ по ее коду Лемера V е и ключу к = 1... п длины | к | =п реализуется следующим алгоритмом: на ^м его шаге, где к = 1,2,..., п, символу пп-к+1 присваивается буква ключа к с номером vn-k+1, а затем она удаляется из к, и новый ключ к имеет длину | к | = п - к.
Для рассматриваемого класса Г-моделей (£(п), вп, Т(пу) существует простой алгоритм генерации множеств номеров Х<П) элементов таблиц Т^.
Теорема 1. Пусть V = v1 ...vn е ^ и ц = Vп. Тогда множество Х-пл при п е N формируется по следующему правилу: каждому V е отвечают только коды Лемера
V1, V 2, ..., V (Шп(т, ц) + 1) е Х«.
Доказательство этой теоремы несложно проводится методом математической индукции на основе отмеченных выше результатов, а описанный в это теореме алгоритм позволяет также легко генерировать множества перестановок Р^ с использованием алгоритма преобразования кодов Лемера в перестановки .
Отметим, что множество номеров Х(™) совпадает с множеством кодов Лемера 213-избегающих перестановок степени п [2].
При фиксированном т е N с помощью теоремы 1 несложно получить распределение вероятностей на множествах перестановок Р^.
Теорема 2. Пусть перестановка п = п1 ... пп е Р^ имеет код Лемера V = v1 ... vn е Х(П), полученный по теореме 1. Тогда распределение вероятностей на множестве перестановок Р^ задается следующим выражением
п-1 1
Рг(п) = П ■ / 1 Ч 1, п е Р(пт), п е N. (5)
и шт(т, V,.) +1
Доказательство. Формула (5) находится через условные вероятности с учетом равновозмож-ности получения по алгоритму теоремы 1 из предыдущего последующих кодов Лемера перестановок, генерируемых на основе этого алгоритма. Так как справедливо равенство $ (т) Рг(п) = 1, то
выражение (5) задает распределение вероятностей на множестве перестановок Р^. Доказательство закончено.
Из теоремы 2 вытекает, что только при т = 1 распределение вероятностей (5) является равномерным. Также равномерное распределение вероятностей несложно определить на множестве всех перестановок степени п, которые не описываются теоремой 2.
Также для рассматриваемого случая рассмотрим Т([) -модель, состоящую из таблиц Тпт)(?), где 0 < [ < 1, а элементами таблиц Т{пт)(с[) служат [р(у), где р(у) = Iп=1 (у/ - 1), V е ^<п™). Производящие многочлены таблиц Тпт\[) задаются выражением Н(пт)([) = I (т) [рМ, причем
н (пт)(1) = и (пт)(1).
Простую взаимосвязь между ^-многочленами, вводимыми с помощью инверсий перестановок
( п ]
п е Рпт), и многочленами Н(пт)([) устанавливает равенство ту(п) + р(у) = , в котором V е Х<т)
код Лемера перестановки п = п1... пп е Рп(т), а ту(п) = #{(/,у) : п1 > Пу, 0 < I < у < п} - число ее инверсий [2].
Для этой Т ([) -модели комбинация ^-аналогов соотношений (1) и (4) приводит к рекуррентной формуле для многочленов Н (пт)([)
Н0т)([) = Н 1т)([) = 1, Нпт)([) = I д'Н(т)([)Нпт)-1([), 2 < п < т;
I т 1+1 (6)
Г 21 ( т - к + 2 ] Нпт)([) (-1)к-1[к(к-1) Г к + 2] Н«([), п > т.
к=1 V к )[
В выражении (6) для ^-аналогов биноминальных коэффициентов (многочленов Гаусса), подробно рассматриваемых в [4], используется следующее обозначение:
т
т] П <' - [)
к)
V У[
т-к
П (1 - [) П (1 - [)
Теперь перейдем к заданию и изучению второго интерполяционного класса числовых последовательностей.
I=0
2. Интерполяция числовых последовательностей между 2я"1 и Вп
Аналогично вышеприведенному исследованию класс целочисленных последовательностей, осуществляющий при п е N интерполяцию между последовательностью степеней 2п-1 и последовательностью чисел Белла Вп, строится на базе соответствующего однопараметрического класса Г-моделей по методике статьи [2].
Этот класс определяется тройкой (Б(т), ет, Т(т)) и при фиксированном значении параметра т е N каждая таблица Г-модели с номером п е N также определяется рекурсивным выражением
Т„т1 = ет(Тп ), причем для этого используется алфавит Б(т) = {2,..., т +1}, отображение
ет : 5 ^ 55-1 (шт(т, 5) + 1) символа 5 е Б(т) в строку символов ет (5) и начальная таблица Т/т) = (2).
В этом случае для преобразования последовательности таблиц Т^ в числовую последовательность используем производящие многочлены
) = I г-1 = ! ит ^,
5еТпт)
(7)
к=1
а значения ипт)(1) определяют искомые числа.
По описанной Г-модели несложно установить рекуррентные соотношения для вычисления коэффициентов и^Ск многочленов (7)
и<Т) = 1? ит = кит к + и« к_151 < к < п, 2 < п < ш;
и<3 = кипш)к + ипт)к-1, 1 < к < Ш -1, иПШШШ = (Ш + 1)иПШ1,Ш + ипт)Ш_1, п > ш; (8)
ипт = о, ш + 1 < к < п, п > ш. Применение формул (7) и (8) позволяет записать числа ип(ш)(1) в компактной форме:
ип")(1) = IЦ<9)
так как числа Стирлинга второго рода удовлетворяют простому рекуррентному соотношению [3]
ЙММ-» -
что позволяет интерпретировать числа (9) как количество разбиений множества {1,2,..., п} не более чем на Ш + 1 блоков.
При Ш = 1,2,3,4,5,... последовательности чисел, вычисляемые по формуле (9), принимают вид:
1, 2, 4, 8, 16, 32, 64, 128, 256, 512, 1024, 2048, 4096, ... [10, А000079];
1, 2, 5, 14, 41, 122, 365, 1094, 3281, 9842, 29525, 88574, ... [10, А007051];
1, 2, 5, 15, 51, 187, 715, 2795, 11051, 43947, 175275, ... [10, А007581]; > (11) 1, 2, 5, 15, 52, 202, 855, 3845, 18002, 86472, 422005, ... [10, А056272];
1, 2, 5, 15, 52, 203, 876, 4111, 20648, 109299, 601492, ... [10, А056273];
^..........................................,
а для Ш = 1,2,..., 15 различные сведения о них содержатся в статьях [10].
Так как при Ш ^ да и п е N формула (9) задает последовательность чисел Белла Вп = и[м)(1), то класс последовательностей (11) осуществляет интерполяцию между последовательностью 2п-1 и последовательностью Вп.
Определение 1 позволяет сопоставить каждому числу (9) множество перестановок Р^ символов {1,2,...,п} такое, что |рпш)| = = и(пт)(1), а для рассматриваемого в этом разделе класса Г-моделей (Б(ш), 0Ш, Т(т)) имеется простой алгоритм генерации множеств номеров Х(Щг) элементов таблиц ^
Теорема 3. Пусть V = v1 ...vn е Х(Щ!) и ц = шах^,...,vn). Тогда множество Х(Щ)1 при п е N аналогично теореме 1 формируется по правилу: каждому V е Х<Ш) отвечают только коды Лемера V1, V 2, ..., V (Шп(т, ц) + 1) е X«.
Доказательство теоремы 3 также проводится методом математической индукции, а описанный в
это теореме алгоритм позволяет легко генерировать множества перестановок Рп(т) с использованием алгоритма преобразования кодов Лемера Х(Щ!) в перестановки Р^.
При фиксированном т е N теорема 3 позволяет с помощью кодов Лемера Х(Ш) найти соответствующее распределение вероятностей на множествах перестановок .
Теорема 4. Пусть перестановка п = п1 ... пп е Р^ имеет код Лемера V = v1 ... vn е Х(Щг), полученный по теореме 3, а ц, = шах^,...VI), , = 1,2,...,п - 1. Тогда распределение вероятностей на _ .. _ _ п можно аналогично теореме 2 определить выражением:
я-1 1
Рг(п) = П—7-ГТ' п е р(пт), п е N. (12)
^ шт( т, ц(.) +1
Доказательство. Формула (12) находится через условные вероятности с учетом равновозмож-ности получения по алгоритму теоремы 3 из предыдущего последующих кодов Лемера перестановок, генерируемых на основе этого алгоритма. Так как справедливо равенство X Р(т) Рг(п) = 1, то
пеРп
выражение (12) задает распределение вероятностей на множестве перестановок Р^. Доказательство закончено.
По теореме 3 при т ^ да определяется множество номеров Х<™), совпадающее с множеством слов ограниченного роста (КС-слов) [8, 9].
Поэтому множества при т е N можно трактовать как га-обобщения множеств КС-слов. Пусть П(т обозначает множество всех разбиений множества {1,2,..., п} не более чем на т + 1 блоков. Так как по формуле (9) имеем равенство |П ^ = )|, то аналогично статье [6] устанавливается биекция между П(т и .
Для этого каждому номеру V = v1 ... vn е ставится в соответствие упорядоченное разбиение из П(пт) по следующему правилу: в первый блок этого разбиения включаются в порядке возрастания все индексы символов v¡ = 1, во второй - все индексы символов V у = 2 и т.д.
Например, слову V = 1121323 е ставится в соответствие разбиение {1,2,4}{3,6}{5,7} е П(4), а слову V = 1234254 е Х(74) - {1}{2,5}{3}{4,7}{6} е П(74).
Алгоритм преобразования кода Лемера V = v1 ... vn е ^ в перестановку п = ... пп е Р^ также определяет биекцию между П(т и Рп(т), что позволяет использовать соотношение (12) для задания распределения вероятностей также и на множестве разбиений П^.
Так как д-аналогом рекуррентного соотношения (10) для д-чисел Стирлинга второго рода служит выражение [6, 9]:
Гп] (кЛ \п -1] Гп -1 ]
V \, Нк 1 1 -I,+{п - 1|. ■
то использование для рассматриваемого случая Т(,) -модели приводит к формуле для многочленов
н (:\д)
т+1 (к\ | м
нпт)(,) = Х,(2 )|п
к I
,
являющейся д-аналогом соотношения (9).
следует отметить, что ряд доказательств в работах [6, 8, 9] могут быть улучшены за счет применения алгоритмического подхода, основанного на использовании 7-модели, кодов Лемера перестановок и их упорядочения как номеров элементов Т-модели [1, 2], а также соответствующей Т(,) -модели.
к=1
Заключение
Применение методики интерполяции числовых последовательностей, разработанной в [2] и основанной на использовании подходящих Г-моделей, позволило построить два простых однопараметри-ческих интерполяционных класса таких последовательностей.
Основной особенностью этих классов является наличие простых алгоритмов, описанных в теоремах 1 и 3 и служащих для построения при фиксированном значении параметра соответствующих множеств перестановок, интерпретирующих элементы этих классов. Эти алгоритмы и их аналоги, в частности, могут быть использованы в информатике для генерации разнообразных массивов слов.
На базе этих алгоритмов в теоремах 2 и 4 также на соответствующих множествах перестановок получены распределения вероятностей, которые могут найти применение для исследования различных статистик на этих множествах перестановок.
В разделе 2 для элементов рассматриваемого однопараметрического класса последовательностей при фиксированном значении параметра при n е N построена интерпретация его элементов в виде соответствующих множеств упорядоченных разбиений множества {1,2,..., n}. Для этих множеств также найдены распределения вероятностей.
Полученные в статье результаты могут найти применение в различных прикладных областях.
Список литературы
1. Бондаренко Л.Н. Систематизация комбинаторных последовательностей с использованием Г-моделей и Г-диаграмм // Образовательные ресурсы и технологии. - 2020. - №2 (31). - С. 58-68.
2. Бондаренко Л.Н. Применение Г-моделей к интерполированию целочисленных последовательностей // Образовательные ресурсы и технологии. - 2021. - №3 (36). - С. 97-105.
3. Грэхем Р., Кнут Д., Паташник О. Конкретная математика. Основание информатики / перевод с английского. - Москва: Мир, 1998. - 703 с.
4. Эндрюс Г. Теория разбиений / перевод с английского. - Москва: Наука. Гл. ред. физ.-мат. лит., 1982 -256 с.
5. Barcucci E., DelLungo A., Pergola E., Pinzani R. Permutations avoiding an increasing number of length -increasing forbidden subsequences // Discrete mathematics and theoretical computer science. - 2000. - No. 4. -P. 31-44.
6. Cai Y., ReaddyM.A. ^-Stirling numbers: A new view // Advances in applied mathematics. - 2017. - Vol. 86. -P. 50-80.
7. Labette G., Leroux P., Pergola E., Pinzani R. Stirling numbers interpolation using permutations with forbidden subsequences // Discrete mathematics. - 2002. - No. 1-3. - Vol. 246. - P. 177-195.
8. Milne S. Restricted growth functions and incidence relations of the lattice of partitions of an n-set // Journal advances in mathematics. - 1977. - Vol. 26. - P. 290-305.
9. Milne S. Restricted growth functions, rank row matchings of partition lattices, and ^-Stirling numbers // Journal advances in mathematics. - 1982. - Vol. 43. - P. 173-196.
10. Sloane N.J.A. The on-line encyclopedia of integer sequences. - 2021. - http://oeis.org.
11. Stanley R.P. Catalan numbers. - New York: Cambridge university press. - 2015. - 215 p.
References
1. Bondarenko L.N. Sistematizaciya kombinatornyh posledovatel'nostej s ispol'zovaniem Г-modelej i Г-diagramm // Obrazovatel'nye resursy i tekhnologii. (Elektronnyj nauchnyj zhurnal). - 2020. - №2 (31). -S.58-68.
2. Bondarenko L.N. Primenenie T-modelej k interpolirovaniyu celochislennyh posledovatel'nostej // Obrazovatel'nye resursy i tekhnologii. (Elektronnyj nauchnyj zhurnal). - 2021. - №3 (36). - S. 97-105.
3. Grekhem R., Knut D., Patashnik O. Konkretnaya matematika. Osnovanie informatiki / perevod s anglijskogo. - Moskva: Mir, 1998. - 703 s.
4. Endryus G. Teoriya razbienij / perevod s anglijskogo. - Moskva: Nauka. Gl. red. fiz.-mat. lit., 1982 - 256 s.
5. Barcucci E., Del Lungo A., Pergola E., Pinzani R. Permutations avoiding an increasing number of length -increasing forbidden subsequences // Discrete mathematics and theoretical computer science. - 2000. - No. 4. -P. 31-44.
6. Cai Y., Readdy M.A. ^-Stirling numbers: A new view // Advances in applied mathematics. - 2017. - Vol. 86. -P. 50-80.
7. Labelle G., Leroux P., Pergola E., Pinzani R. Stirling numbers interpolation using permutations with forbidden subsequences // Discrete mathematics. - 2002. - No. 1-3. - Vol. 246. - P. 177-195.
8. Milne S. Restricted growth functions and incidence relations of the lattice of partitions of an n-set // Journal advances in mathematics. - 1977. - Vol. 26. - P. 290-305.
9. Milne S. Restricted growth functions, rank row matchings of partition lattices, and ^-Stirling numbers // Journal advances in mathematics. - 1982. - Vol. 43. - P. 173-196.
10. Sloane N.J.A. The on-line encyclopedia of integer sequences. - 2021. - http://oeis.org.
11. Stanley R.P. Catalan numbers. - New York: Cambridge university press. - 2015. - 215 p.