Параметрические комбинаторные задачи и методы их исследования Текст научной статьи по специальности «Математика»

УДК 519.1

Л. Н. Бондаренко, М. Л. Шарапова

ПАРАМЕТРИЧЕСКИЕ КОМБИНАТОРНЫЕ ЗАДАЧИ И МЕТОДЫ ИХ ИССЛЕДОВАНИЯ1

Аннотация. Рассматриваются четыре комбинаторные задачи, параметризованные кратностью r элемента базового мультимножества: распределение индексов vp-монотонных перестановок, обобщенные перестановки Гесселя - Стенли и обобщенные частично упорядоченные множества Баклавского - Эдель-мана, обобщенные числа Стирлинга и обобщенные частично упорядоченные множества разбиений, обобщенные статистики и обобщенные многочлены Эйлера. Для исследования этих задач привлекаются различные методы.

Ключевые слова: мультимножество, статистика, производящая функция, vp-монотонные перестановки, перестановки Гесселя - Стенли, посеты Баклав-ского - Эдельмана, числа Стирлинга, посеты разбиений, многочлены Эйлера.

Abstract. Four combinatorial problems parametrized by multiplicity r of an element of base multiset are considered: distribution of indexes vp-monotonous permutations, the generalized Gessel - Stanley permutations and the generalized Baclawski -Edelman partially ordered sets, the generalized Stirling numbers and the generalized partially ordered sets of the partitions, the generalized statistics and generalized Eul-erian polynomials. For research of these problems various methods are involved.

Keywords: multiset, statistic, generating function, vp-monotonous permutations, Gessel - Stanley permutations, Baclawski - Edelman posets, Stirling numbers, posets partitions, Eulerian polynomials.

Введение

Разнообразие комбинаторных задач заставляет унифицировать методы их описания и исследования. Прогресс в этой области хорошо прослеживается по монографиям [1, 2], в которых анализируется большое число перечислительных задач. Одним из простейших подходов к унификации является исследование семейств комбинаторных задач, зависящих от некоторого параметра. В качестве такого параметра будем рассматривать кратность r > 1 каждого элемента мультимножества {1r ,2r,..., nr }, которое служит базой для постановки ряда комбинаторных задач. Это мультимножество для краткости будем обозначать [nr ], где n является целым положительным числом.

Все перестановки множества [nr ] образуют множество SPnr мощности #{o: oe SPn r} = (r !)_n(rn)!. Каждую перестановку oe SPn r удобно рассматривать как слово 0 = 01... orn длины | o | = rn, причем SPn r содержит и зеркальный образ (mirror image) mi(o) e SPnr слова o, т.е. o, записанное в обратном порядке. Например, в тождественной перестановке ее SPnr,

трактуемой как слово е = 1r ...nr, запись ir, i е [n], означает степень символа

i, полученную конкатенацией r символов i, а mi(e) = nr ...1r.

1 Работа выполнена при финансовой поддержке РГНФ (проект № 09-0628615 а/В).

На подмножествах перестановок из SPn r в комбинаторном анализе

часто определяют числовые функции и рассматривают их распределения относительно равномерной меры на этих подмножествах. Такие функции в работах [2, 3] называются статистиками.

При рассмотрении упорядоченных множеств вместо словосочетания «локально конечное частично упорядоченное множество» будем использовать, как и в работе [4], термин «посет», соответствующий английскому сокращению «poset» (partial ordered set).

1. Распределение индексов vp-монотонных перестановок

В криптографии из перестановки ае SPn r с помощью фиксированного ключа K.eSPnr образуют новое слово T = Tj ..лт, ieTnr, | т | = rn, над алфавитом [n] следующим образом: т = а ©к, где т- = а- + к- (modn),

i = 1,..rn, а т является наименьшим положительным вычетом.

Это преобразование, определяющее биекцию SPn r на множество слов Tnr, назовем vp-отображением, т.е. vp : а ^ т (vp - сокращение словосочетания «vector permutation»). Так как нетривиальные перестановки рассматриваются при n > 2 , то в качестве ключа будем фиксировать перестановку s или mi(e).

Определение 1. Назовем слово те Tnr монотонным, если все его последовательные символы образуют неубывающую (невозрастающую) последовательность, а соответствующую перестановку а = vp 1(т), ае SPnr , назовем vp-монотонной перестановкой.

Соответствие между словами т = а © е и mi(т) = mi(a) © mi(e) позволяет ограничиться рассмотрением только множества mon(SPn r ) vp-монотонно неубывающих перестановок, причем в отображении vp : mon(SPn r ) ^ ^ mon(Tn r ) используется ключ к = mi(e).

Определение 2. Статистику ivp(a) = (nr)-1 ^^т-, где ае mon(SPn r),

а т = vp(a), назовем индексом vp-монотонно неубывающей перестановки.

Эта статистика при r = 1 была введена в работе [3]. Она находит среднее значение символа в слове т = vp(a) и, несмотря на нерегулярность vp-

отображения на множестве mon(SPn r ), обладает рядом закономерностей.

Лемма 1. а) Индексы ivp(a) перестановок aemon(SPnr ) принимают только целые значения k = 1,...,n.

б) Если ае mon(SPnr ) и ivp(a) = k, то в слове а = ‘tínr^ префикс £ и суффикс n длины | ^ | = r (k -1) являются любыми монотонно неубывающими словами такими, что все эти а могут быть лексикографически упорядочены следующим образом: 1r .(n - k )r nr (n - k + 1)r... (n - 1)r, ..., kr ...nr 1r ...(k - 1)r .

То есть

Доказательство. С помощью равенства т = а ©к, к = ш1(е), для перестановок а = кг ...пг 1г... (к — 1)г находим т = кгп , поэтому ^р(а) принимает значения к = 1,., п. Пусть а е шоп(БРп г), а в слове т = т. тпг символы т = 1 при а + кг- > п и т = 0 при аг- + кг- < п . Тогда статистика егур(а) = г 1 ^ *= тг- дополнительна к статистике гур(а), т.е. справедливо равенство гур(а) + егур(а) = п + 1. Так как слово т = 1г (п-к+1)0г (к-1) для всех перестановок, описанных в лемме 1,б), и эта форма слова т сохраняется для любой ае шоп(8Рп г), то лемма 1 доказана.

Теорема 1. Пусть Упгк = #{а: ае шоп^Рп г), гур(а) = к} , ке [п]. Тогда имеет место следующее рекуррентное соотношение:

У0, г,к = 0, У1, г,к = 81к,

Уп,г,к = Уп-1,г,к + Уп-1,г,к-1 + (г — 1)Уп-2,г,к-1, п — 2, ке ^ , (1)

где 8у - символ Кронекера (8гу = 1 при 1 = у и 8гу = 0 при 1 Ф у ).

Доказательство. Применим метод математической индукции. При г — 1 равенства У0гк = 0, У1гк = 81 к находятся непосредственно. Если для некоторого п — 2 соотношение (1) выполнено, то оно выполняется и при п +1. Действительно, с помощью леммы 1 все слова аешоп(5Рп+1г), гур(а) = к, можно разбить на следующие блоки: 1) получается вставкой слова (п + 1)г перед суффиксом длины г (к — 1) в каждое слово лешоп(£Рпг), гур(я) = к ; 2) получается вставкой (п + 1)г перед суффиксом длины г (к — 1) в каждое яе шоп(5Рп г), гур(я) = к — 1 и перемещением подслова пг в конец слова; 3) получается вставкой (п + 1)г перед суффиксом длины г (к — 2) в каждое яе шоп(£Рп_ 1 г), гур(я) = к — 1, вставкой пг в конец слова и обменом

1 = 1,.., г — 1 символом между подсловами (п — 1)г, пг так, чтобы выполнялись условия леммы 1.

Для производящих многочленов Упг (() = ^ п=1 Уп г к^к соотношение (1) позволяет легко найти рекуррентную формулу:

У0, г (0 = 0, У1, г (0 = г, Уп, г (0 = (1 + 0 Уп—1, г и) + (г — 1) tУn—2, г (0, п — 2, (2)

которая влечет равенство Уп г ^) = ^+1 Уп г (t 1), показывающее, что t—1Уп г ^) - возвратный многочлен, иначе Упгк =Упг п—к+1.

п

п_0)Упг^г определяется стандартным образом с помощью выражения (2) и оказывается рациональной

[7

Уг ([, 7) =------------------------- , (3)

1 - (1 + [) 7 - (Г - 1) [7

а с помощью разностного уравнения (2) или функции (3) несложно найти явное выражение для многочленов Уп г ([):

У„ г (ґ) =

( ґ I----------------\П ґ I--------------\п ^

1 + ґ + УІ (1 - ґ)2 + 4гґ 1 + ґ-д/ (1 - ґ)2 + 4гґ

7(1 - ґ )2 + 4 гґ

, (4)

уи,г(1) = (^) ((1 + ()п -(1 -4~г)п)

что дает Упг(1) = (2^г) ((1 + л/г) -(1 -Vг) j - число ур-монотонных пе-

рестановок. В частности, Уп 1 (1) = 2п-1, а Уп 5(1) = 2п-1 ¥п , где Еп - числа Фибоначчи.

Разложение функции (3) по степеням [ приводит к выражению

^ п!у 7П )[фг (7)]* * >0 () 1 7 + (г - 1)72

^ 77Уп,г,*+1 —т = ¥(7)---- ----, 7 >0 ,где у(7) = -, Фг(7) =--- -------,

пк ! п! 7! 1 - 7 1 - 7

п=0

а ^(0) Ф 0, фг (0) = 0, ф'г (0) Ф 0 , т.е., по терминологии Л. М. Коганова [4], последовательность чисел {7! 1 п!Упг *+1} является псевдопорождаемой с порождающей функцией фг(г) и образующей ^(г).

Псевдопорождаемые последовательности и производящие функции являются мощными инструментами решения ряда задач перечислительной комбинаторики [4, 5]. Таким образом, соотношения (2)-(4) дают исчерпывающее описание многочленов Уп г ([) и их коэффициентов Упг *, определяющих распределение индексов ур-монотонных перестановок.

Следующее утверждение позволяет использовать вероятностные методы работы [3] для получения асимптотики чисел Упг * при п ^^ .

Теорема 2. Многочлены [-1Уп-1 г (О и [-1Уп г ([) (п > 2) образуют положительную пару, иначе все их нули - различные отрицательные числа, взаимно разделяющие друг друга, т.е. между любыми двумя нулями многочлена [-1Уп г ([) имеется один нуль многочлена [-1Уп-1 г ([).

Доказательство. В промежутке (~тс,0] последовательность

{[-1У* г ([)}7п=1 образует ряд многочленов Штурма, что следует из соотношения (2). Разность числа перемен знаков в значениях многочленов этого ряда, вычисленных при [ = -°° и [ = 0, равна п -1, что и доказывает теорему 2.

2. Обобщенные ГС-перестановки и обобщенные БЭ-посеты

При изучении полиномиальных последовательностей Стирлинга обоих родов И. Гесселем и Р. Стенли в [6] были определены специальные перестановки множества [п2], которые в [7] названы ГС-перестановками.

Определение 3. Обобщенной ГС-перестановкой будем называть слово ає БРп г , г > 1, обладающее ГС-свойством: все буквы слова а, стоящие между любыми двумя вхождениями символа і є [п], не меньше этого і.

Под определение 3 при г = 1 подходят обычные перестановки, а при г = 2 - перестановки Гесселя - Стенли. Множество всех обобщенных ГС-перестановок а є БРп г обозначим ОБп г. В работе [8] рассмотрены некоторые свойства обобщенных эйлеровых статистик на множестве ОБ п г .

По определению 3 легко строится итерационный алгоритм генерации ОБп+1 г, основанный на нахождении (гп + 1)-го слова из ОБп+іг путем

вставки (п + 1)г в выбранное слово кє ОБпг, а #{а: ає ОБп+іг} = = 1 • (г +1) •... • (гп +1).

Определение 4. Фиксируя г > 1 и п > 1, сформируем множество Рп г всех подмножеств {І1, І2,..., ігк } с Z таких, что 0 < І1 < І2 < ... < ігк < гп +1, к є [п], а числа І1, І2 - І1,.., ігк - ігк-1 имеют вид гу +1, у = 0,1,... Упорядоченное по включению множество Рп г будем называть обобщенным БЭ-посетом.

Под определение 4 при г = 1 подходит булева решетка Вп, а при г = 2 - посет К. Баклавского и П. Эдельмана, описанный в [2].

Диаграмма Хассе посета Рп г строится на базе множества [гп] - максимальный (единичный) элемент. На первом шаге все покрываемые им множества находятся по определению 4 вычеркиванием любых г рядом стоящих элементов [гп], что можно выполнить г(п -1) +1 способом. Если на к-м шаге все вершины у диаграммы Хассе построены, то покрываемые ими вершины х < у также находятся вычеркиванием из каждого множества у любых г рядом стоящих его элементов. Пустое множество 0 служит минимальным (нулевым) элементом.

Таким образом, посет Рп г имеет единственную ранговую функцию

р: Рпг ^ [п], для которой р(0) = 0, р([гп]) = п . Так как все максимальные

цепи посета Рп г имеют одинаковую длину п, то множество Рп г градуированное, а число максимальных цепей равно (г(п -1) +1) •... • (г +1) • 1.

Определим посет Рп г, диаграмма Хассе которого строится так же, как

и для посета Рпг, но ее вершинами являются слова, а не множества, причем слово £ = 1г...пг заменяет базовое множество [гп] = {1,2,.,гп}.

Лемма 2. Посеты Рпг и Рпг изоморфны.

Доказательство. Построение диаграмм Хассе посетов Рпг и Рпг показывает, что достаточно задать биекцию между словом е = Є1... £т и мно-

жеством [rn] = {/'і,...,im} . Она определяется формулами ik = r(є^ -1) + m , k = 1,.., rn , где m - наименьший положительный вычет числа k по mod r . Отметим, что эта биекция позволяет преобразовать каждое слово а є GS n r

в слово а є GSn,r.

По лемме 2 число максимальных цепей в посетах Pnr и Pnr одинаково и совпадает с

___ rin + r-11

card GS nr = card GS nr = rn —^,

, ' r(r-1)

где Г(z) - гамма-функция.

Покажем, что существует алгоритм маркировки каждой максимальной

цепи посетов Pnr и Pnr соответствующими словами ає GS nr и ає GSn,r.

Переход от вершины у є Pnr ранга (n - k) диаграммы Хассе к вершине x є Pnr, x < у, ранга n - k -1, k = 0,., n -1, состоит в вычеркивании r рядом стоящих элементов множества у. Считая, что вычеркнутые числа задают номера символов в слове а = а1,..., ат, вставим на эти r мест символы

(n - k). Тогда при переходе от единицы [rn] к нулю 0 посета Pn r найдем

слово ає GSnr, которое легко преобразуется в слово ає GSn,г. Так как этим методом находятся все перестановки ає GSn r , то построен второй алгоритм генерации GSn r.

В рангово-производящей функции U n r (t) = ^ n=oUn r ktk коэффициенты Unr k называются числами Уитни второго рода [2].

Теорема 3. Для посетов Pnr и Pnr числа Уитни второго рода

Ґ n + (r - 1)k Л

U

n,r,k

(5)

гк

Ч У

Доказательство. На шаге к = 0,1,.,п вершины диаграммы Хассе посета Рп г находятся вычеркиванием гк элементов из множества [гп] так,

чтобы в результате было не более к пробелов. Поэтому на к-м шаге имеем ' г (п - к) + к Л

вершин, а замена к на (п - к) окончательно дает выражение (5).

^ г (п - к) )

Для вычисления общего числа элементов посетов Рпг и Рпг введем

п

¡Риг-(ґ)г , полагая (ґ) = 1.

Теорема 4. Функция иг (ґ,г) рациональна и имеет вид

(1 - г)г -1

и г (ґ, г) =-(-------------------------------------г-. (6)

(1 - г) - ґг

Доказательство. Подставляя в иг (ґ, г) выражение (5), изменяя порядок суммирования и заменяя п на п + к , получим

иг (ґ, г) = £ (ґг)к £ ( ” + * 1 г”

к=0 п=0 V гк

Затем, применяя тождество

—1__=£ [п+т 1 гп,

(1 - 2)т+1 П=01 т )

окончательно приходим к формуле (6).

Обобщенные многочлены Фибоначчи ¥п г (^) можно определить рекуррентным соотношением:

^ (() = ... = ^-2,г ^) = 0, ¥г-1,г (^) = 1, ¥п г (^) = (*) + ^п-г г (^), п > г .

п=0 ¥п г (^) 2 легко вычисляется

¥г (¿, 2) = 2г-1(1 - & - 2Г )-1, а разложение ¥г (1,2) = 2-1 £ ¿=1(гГ (1 - 2)-1)к показывает, что ¥п-1 г (1) равно числу композиций п, все части которых не меньше г [9].

Так как для производящей функции (6) справедливо равенство

иг (Г, 2) = г-1 £ (1 - (е2га'^)17 г - 2)-1, I = Т-Г, к =1

связывающее ее с функцией (1 -& - 2Г)-1, то ип г (^) = ¥г(п+1)-1 г (г1 г), а общее число элементов посетов Рпг и Рпг равно ип г (1) = ¥г(п+1)-1 г (1).

Посет Рпг однороден в смысле [2], и матрицы (ип г к УН к=0,

(ип г п-к )тк=0 взаимно обратны (ип г к называется к-м числом Уитни первого рода).

Разложение функции (6) по степеням ^ приводит к псевдопорождаемой

последовательности {к!-1 пШпг к} с порождающей функцией

- г -1

фг (2) = 2(1 - 2) и образующей ^(2) = (1 - 2) . Если Фг =Фг (^),

Фг (0) = 0 - решение уравнения Фг = ^(1 -Фг)г, то паре (^(2), фг (2)) отвечает пара ([^(Фг (^))]-1, Фг (^)), задающая псевдопорождаемую последова-

тельность, обратную данной [4]. В частности, при к = 0 получим последовательность {п!ипгп}, позволяющую вычислить функцию Мебиуса

цп г = ц(0,1) = ип,г,п посета Рп г .

В рассматриваемом случае по теореме Лагранжа [5] находим ]п-1,

n

Фг (W = £ W-

л n! n=1

1(1 - z У

п-1

а [у(Ф г (w))] 1 = 1 -Ф г (w),

J Z=0

что устанавливает связь функции цпг с числами Фусса - Каталана Cn

М1n,r ( 1) Cn,г , Cn,г

1

(г - 1)n + 1

V n /

комбинаторный смысл которых подробно рассматривается в [10]. Число плоских деревьев с висячим корнем и с (r - 1)n +1 некорневыми вершинами,

имеющими степень, сравнимую с 1 по mod (г -1), также равно Cn r [5].

3. Обобщенные числа Стирлинга и обобщенные посеты разбиений

Известно, что общие задачи размещения и занятости для случая одинаковых ячеек весьма сложны [1]. Поставим задачу о нахождении числа способов Sn r k размещения n объектов с номерами 1,...,n по k одинаковым ячейкам с упорядочением уровня r при условии, что ни одна из них не остается пустой.

При r = 0 упорядочение отсутствует и Sn о k = Snk, где Snk - числа Стирлинга второго рода; при r = 1 объекты в ячейках упорядочены и n! (n -1^

- числа Лаха без знака [1]. При упорядочении уровня r > 1

Sn,1,k J I

k!

k-1

каждый объект, помещенный в ячейку, копируется г раз, а слова, составленные из номеров объектов каждой ячейки, обладают ГС-свойством.

Используя индукцию по п, приходим к рекуррентному соотношению

Б 0, г, к = §0 к, Б п, г, к = (г (п 1) + к )Бп-1, г, к + Бп-1, г, к -1, п > 1, к е^ , (7)

поэтому будем называть Бп г к обобщенными числами Стирлинга второго рода.

Каждому рассмотренному размещению при г = 0 отвечает обычное разбиение множества [п], а при г > 1 разбиение множества [пг] на упорядоченные блоки, обладающие ГС-свойством, причем Бп г 1 = #{а: ае ОБп г} .

Определение 5. Пусть Пп г - множество всех разбиений [пг] типа, определенного значением г > 0 . Обобщенным посетом разбиений Пп г назовем частично упорядоченное по измельчению множество Пп г, для которого

я<а, если каждый упорядоченный блок п содержится в упорядоченном блоке а.

Например, по этому определению в посете П3 2 упорядоченное разбиение {122331} покрывает три разбиения {11},{2233}; {22},{1331}; {33}, {1221} .

Из определения 5 следует, что посет Пп г - градуированное множество ранга (п -1), а р(п) = п -1 п |, где | п | - число блоков разбиения п, причем при г > 1 этот посет не содержит единицы. С помощью (7) несложно проверить,

что многочлен Бпг (X) = £п=1 Бпгк(к , связанный с рангово-производящей

функцией посета Пп г формулой XпБп г(X 1), удовлетворяет рекуррентному соотношению

б0,г (0 = 1, Бп,г (0 = (г (п -1) + X)Бп-1,г (0 + ХБ'п-1,г (X), п > 1, (8)

т.е. Бп г п-к - к-е число Уитни второго рода посета Пп г .

Теорема 5. Для многочленов Бпг (X) справедлива формула

Бпг (X) = е-Х (СгпИпг ) еХ, Нг = Xг+1 М, (9)

аХ

_1 п

п=0 Бпг (X) п! г имеет вид

Б0(Х,г) = ехр(Х(ег -1)), Бг(X,г) = ехр(Х((1 - гг)-1/г -1)) при г > 1, (10)

-1/г

где используется главное значение функции (1 - гг) .

Доказательство. Формула Бп г (X) = Х-гпе-г (Хг+1М /Мх)^11-1^^п-1 г (X))

следует из (8), а ее итерирование дает (9). Сравнение левых и правых частей в выражениях (10) с применением (9), тождества

(X-гп Нпг) Хк = к (к + г)... (к + г (п -1)) Хк

и разложений в степенные ряды функций ег и (1 - гг) 1/г приводит к (10).

Из доказательства теоремы 5 также вытекает, что последовательности {Бп 0 к} и {Бпг к}, г > 1, порождаемы с порождающими функциями

Ф0(г) = ег -1 и фг(г) = (1 -гг)-1/г -1 соответственно (образующие равны единице).

-1 - г

Порождающие функции Ф0^) = 1п(1 + w) и Фг (w) = г (1 - (1 + w) ),

г > 1, обратны к ф0(г) и фг(г) соответственно. Согласно [4] они порождают последовательность чисел {япг к}, обратную последовательности {Бп г к}. Числа ,5пг к являются обобщением чисел Стирлинга первого рода, так как 5п 0 к = яп к, где 5пк - обычные числа Стирлинга первого рода [1].

Многочлен 5п г (у) = £п 5п гкук связан с характеристическим многочленом [2] посета Пп г формулой

(v ^ = Zиеп ^(°’п)уП Р(Л)

при г > 1. (11)

т.е. sn r n_k - k-е число Уитни первого рода этого посета.

С помощью порождающих функций Oq(w) и Фr (w) несложно найти

производящие функции s r (v, w) = z: =0sn,r(v)n! 1 w , r > °.

s °(v, w) = z (Ф°(w)) kk !_1 vk = (1 + w)v, k=0

^ _ sr (v, w) = z (Фr (w))kk!_1 vk = exp r_1v(-(1 + w)_r j

k=0

При r > 0 формулы (11) влекут рекуррентное соотношение

s0,r,k = §0k, sn,r,k = sn_1,r,k_1 _ (n + rk _ 1)sn_1,r,k, n > 1, kе Z, с помощью которого легко вычисляется

^n_1,r = Sn,r,1 = (_1)n_V !)_1(n + r _1)! -

функция Мебиуса посета Пn r .

Отметим, что числа Snrk были введены Л. Комте для действительного

параметра r при рассмотрении действия оператора H^ из (9) на функцию f (t) [4], где также доказано соотношение Snr k = Z j.k (_rr~Js»jSj.k .

Сравнение выражений (10) и (11) для действительного параметра r Ф 0 дает соотношение sn r (v) = (_r)n S _1 (_r_1v), а из первой формулы (11) получаем известное выражение sn 0 (v) = v(v _ 1)... (v _ n +1). С помощью теоремы 5 по индукции можно показать, что многочлены t_"1Sn_1 r (t) и

t 1Sn r (t), (n > 2) также образуют положительную пару аналогично результату теоремы 2.

4. Обобщенные статистики и обобщенные многочлены Эйлера

В работе [11] исследован ряд статистик на группе перестановок. Расширим действие некоторых из них на множество GSn r. Число подъемов

перестановки является одной из простейших ее характеристик. Для перестановки ое GSnr число подъемов rise(o) = #{i: ог- <ог+1,0 < i < rn _1,00 = 0},

причем статистика crise(o) = #{i: ог- >ог+1,1 < i < rn, orn+1 = 0} , вычисляющая число спусков перестановки ое GSnr , дополнительна к rise(o), т.е. rise(o) + crise(o) = n +1 [8].

Определение 6. Код Лемера £ = l(а) перестановки ає GSnr зададим словом £ = £, ..., £т с буквами £ = °, £ = #{j: аj < а, ,1 < j < i -1},

i = 2,..., rn, а статистику imal(a) = ima(£), где ima(£) - число различных букв (integer make) в слове £, будем называть по аналогии с [11] обобщенной статистикой Дюмона.

Лемма 3. Отображение l: а ^ £ биективно.

Доказательство. По определению 6 перестановке ає GSnr соответствует слово £ с числом вхождений каждого символа, кратным г, причем позиции возрастаний в словах а и £ совпадают. В соответствии с определением 6 для восстановления по £ слова а є GSn r будем заполнять пустые (пронумерованные) позиции искомого слова а по следующему алгоритму.

Если £ содержит символы, отличные от 0, то на k-м шаге, k = 1,2,.,

в слове £ находится подслово mr с таким первым символом m = £ , что при m = i-rj — 1, j = 0,1,., величина n - k - j + 1 максимальна. Затем на г не заполненных ранее позициях слова а записываются символы n — k — j + 1, начиная с i-й позиции. Вычеркивая из £ найденное подслово mr, получаем новое

слово £ длины | £ |=| £ | — r . Затем полагаем £=£ и переходим к шагу k +1.

Если £ содержит все символы, равные 0, то на оставшиеся пустые позиции слова а записываются неиспользованные ранее символы алфавита [n] кратностью r в порядке убывания их величины.

Итак, по слову £ сформирована единственная перестановка а є GSn r .

Определение 7. Числа Anrk = #{а: ає GSnr,гі8е(а) = k}, kє[п], назовем обобщенными числа Эйлера.

При r = 1 по определению 7 получаем обычные числа Эйлера Ank .

Теорема 6. а) Числа Anrk удовлетворяют рекуррентному соотношению

A°,r,k =^°k, An,r,k ~kAn—1,r,k (r(n 1) k ^ 2)An —1,r,k —1, n > 1, k . (12)

б) Anr k = #{а: ає GSn r ,imal(а) = k}, kє [n].

Доказательство. В соответствии с определением 7 и теоремой 6,б rise и imal являются обобщенными эйлеровыми функциями (E-статистиками), определяемыми выражением Anr k = #{а: ає GSn r,Е(а) = k}, kє [n]. Поэтому

для доказательства теоремы 6 следует с помощью статистик rise и imal получить формулу (12). При n = 1 эта формула проверяется непосредственно. Пусть по индукционному предположению для перестановок X є GSn—1 r и

статистик rise, imal числа An—\rk могут быть вычислены с помощью (12).

Тогда для доказательства (12) для чисел Anr k рассмотрим два случая:

1) Если An—1rk = #{х: хє GSn—1 r,E(x) = k} , то вставить слово nr в слова {х: хє GSn—1 r ,rise(x) = k} без изменения числа подъемов rise(x) = k

можно ровно к способами, вставляя его в те ОБп_ г между буквами ту,тг+і, для которых Ту < Ту+1,0 < і <г(п _ 1) _ 1. Аналогично к способами можно вставить слово пг в слова {т: те ОБп_іг ,ішаі(т) = к} без изменения числа іша(£) = к, вставляя его в те ОБп_ г на первом месте или между двумя буквами ту, т+і, 1 < і < г (п _ 1) _ 1, для которых в слове £ = I (т) имеем

£ у ф£у+1,1< ] < у .

2) Если Ап_1гк = #{т: те ОБп_1г,Е(т) = к _1}, то вставить слово пг в слова {т: те ОБп_ г ,гі8е(т) = к _ 1} с изменением числа подъемов гІ8е(т) = к _ 1 на к можно г(п _ 1) _ к + 2 способами, вставляя его в те ОБп_ г на последнем месте или между двумя такими буквами Оу, оу+, что ту > ту+1,0 < у < г (п _ 1) _ 1. Аналогично, г (п _ 1) _ к + 2 способами можно вставить слово пг в слова {т: те ОБп_уг ,іша1(т) = к _1} с изменением числа іша(£) = к _ 1 на к, вставляя его в слово т е ОБп_ г на последнем месте или между двумя буквами ту, ту+1,1 < у < г(п — 1) —1, для которых в слове £ = / (т) имеем £ у = £+1,1< у < у .

Таким образом, если соотношение (12) выполняется для перестановок т е ОБп_1 г, то оно выполняется и для перестановок о е ОБп г.

С помощью формулы (12) для обобщенных многочленов Эйлера Ап г (і) = ^ п=1 Апг к їк проверяется рекуррентное соотношение

Ао,г (() = 1, Ап,г (і) = (г(п _ 1) +1) і Ап_1,г (і) + і (1 _ і) А'п_1,г (і), п > 1, (13)

которое при г = 1 определяет обычные многочлены Эйлера Ап (і), а также легко позволяет вычислить Апг(1) = #{о : ое ОБп г} .

В работе [2] многочлен Q(t) называется /-эйлеровым, если последовательность {/(к)}^ значений многочлена/при т > 0 имеет производящую

функцию 2 Г=о /(к) ік = Q(t)(1 _і) т 1, ¿ее Q < т, причем ¿ее / < т, а ¿ее / = т тогда и только тогда, когда Q(1) Ф 0 .

Обобщенные многочлены Эйлера Ап г (і) подходят под это определение, так как из (13) имеем

Ап,г (і) = (1 _ і) гп+1 (і (1 _ і )^г ё / Л )((1 _ і)_г(п _1)Ч Ап_ 1 (і ,

и итерацией, аналогично доказательству теоремы (5), получаем представление

Ап г (і) = (1 _ іГ+1 Нпг (1 _ і)_1, Нг =-*—гг ё. (14)

(1 _і) г_1 ё

С помощью (14) несложно получить производящую функцию только при г = 1, т.е. для обычных многочленов Эйлера Ап (і)

Использование формулы (14) позволяет по индукции показать, что

пару аналогично результату теоремы 2.

В заключение отметим, что используемый параметр г может быть продолжен на множество целых или множество действительных чисел, но при этом теряется комбинаторный смысл рассматриваемых двухиндексных по-

1. Риордан, Дж. Введение в комбинаторный анализ / Дж. Риордан. - М. : Изд-во иностр. литер., 1963. - 288 с.

2. Стенли, Р. Перечислительная комбинаторика / Р. Стенли. - М. : Мир, 1990. -440 с.

3. Бондаренко, Л. Н. Статистики на классах отображений / Л. Н. Бондаренко, М. Л. Шарапова // Дискретные модели в теории управляющих систем : VIII Меж-дунар. конф. (Москва, 6-9 апреля 2009 г.). - М. : Издательский отдел факультета ВМиК МГУ им. М. В. Ломоносова; МАКС Пресс, 2009. - С. 33-39.

4. Коганов, Л. М. Псевдопорождаемые двухиндексные последовательности / Л. М. Коганов. М. : Недра, 1989. - 86 с.

5. Гульден, Я. Перечислительная комбинаторика / Я. Гульден, Д. Джексон. - М. : Наука, 1990. - 504 с.

6. Gessel, I. Stirling polynomials / I. Gessel, R. P. Stanly // Journal of combinatorial theory. Series A. - 1978. - V 24. - № 1. - P. 24-33.

7. Коганов, Л. М. Универсальная биекция между перестановками Гесселя -Стенли и диаграммами связей соответствующих рангов / Л. М. Коганов // Успехи математических наук. - 1996. - Т. 51. - Вып. 2. - С. 165-166.

8. Бондаренко, Л. Н. Два типа r-перестановок и r-многочлены Эйлера /

Л. Н. Бондаренко, М. Л. Шарапова // Дискретная математика и ее приложения : материалы X Междунар. семинара (Москва, МГУ, 1-6 февраля 2010 г.). - М. : Изд-во механико-математического факультета МГУ, 2010. - С. 217 220.

9. Эндрюс, Г. Теория разбиений / Г. Эндрюс. - М. : Наука, 1982. - 256 с.

10. Грэхем, Р. Конкретная математика. Основание информатики / Р. Грэхем,

Д. Кнут, О. Паташник. - М. : Мир, 1998.

11. Фоата, Д. Распределения типа Эйлера и Макмагона на группе перестановок / Д. Фоата // Проблемы комбинаторного анализа : сб. статей. - М. : Мир, 1980. -

многочлены t An_i r (t) и t An r (t) (n > 2) также образуют положительную

следовательностей Vn,r,k, Un,r,k, Sn,r,k, An,r,k •

Список литературы

С. 120-141.

Бондаренко Леонид Николаевич

кандидат технических наук, доцент, кафедра дискретной математики, Пензенский государственный университет

Bondarenko Leonid Nikolaevich Candidate of engineering sciences, associate professor, sub-department of discrete mathematics,

Penza State University

E-mail: dm@pnzgu.ru

Шарапова Марина Леонидовна

ассистент, кафедра общей топологии и геометрии, Московский государственный университет им. М. В. Ломоносова

Е-шаП: dm@pnzgu.ru

УДК 519.1 Бондаренко, Л. Н.

Параметрические комбинаторные задачи и методы их исследования / Л. Н. Бондаренко, М. Л. Шарапова // Известия высших учебных заведений. Поволжский регион. Физико-математические науки. - 2010. - № 4 (16). -С.50-63.

Sharapova Marina Leonidovna Assistant, sub-department of general topology and geometry, Moscow State University named after M. V. Lomonosov