Числа Эйлера I и II рода — различия между версиями
VolhovM (обсуждение | вклад) |
VolhovM (обсуждение | вклад) |
||
Строка 14: | Строка 14: | ||
Тогда рекуррентная формула имеет вид: | Тогда рекуррентная формула имеет вид: | ||
− | <tex dpi = " | + | :<tex dpi = "160">\left\langle{n\atop m}\right\rangle = (m + 1)\left\langle{n - 1\atop m}\right\rangle + (n - m)\left\langle{n - 1\atop m - 1}\right\rangle</tex> |
Примем также следующие начальные значения: | Примем также следующие начальные значения: | ||
− | <tex dpi = " | + | :<tex dpi = "160">\langle{n\atop m}\rangle = 0</tex>, если <tex dpi = "130">m < 0</tex> |
+ | |||
+ | Также для четных <tex dpi = "130">k</tex>: | ||
+ | :<tex dpi = "160">\left\langle{n\atop m}\right\rangle = [k = 0]</tex>, | ||
+ | Запись [выражение] означает нотацию Айверсона, где | ||
+ | :<tex dpi = "160"> [statement] =</tex><tex dpi = "140"> \begin{cases} 1 & \text{statement} \\ 0 & \text{!statement} \end{cases}</tex> | ||
===Пример=== | ===Пример=== | ||
Рассмотрим все перестановки порядка <tex dpi = "130">4</tex>, в которых есть ровно 2 подъема (в квадратных скобках один или больше подъемов подряд): | Рассмотрим все перестановки порядка <tex dpi = "130">4</tex>, в которых есть ровно 2 подъема (в квадратных скобках один или больше подъемов подряд): | ||
− | <tex dpi = " | + | :<tex dpi = "140"> \left\langle{4\atop 2}\right\rangle = 11: |
[124]3, | [124]3, | ||
[13][24], | [13][24], | ||
Строка 39: | Строка 44: | ||
Согласно алгоритму вывода рекуррентной формулы мы можем добавить '4' в следующие позиции всех перестановок порядка <tex dpi = "130">3</tex> с двумя подъемами, не увеличив количество подъемов: | Согласно алгоритму вывода рекуррентной формулы мы можем добавить '4' в следующие позиции всех перестановок порядка <tex dpi = "130">3</tex> с двумя подъемами, не увеличив количество подъемов: | ||
− | <tex dpi = " | + | :<tex dpi = "140"> |
\left\langle{3\atop 2}\right\rangle = 1: | \left\langle{3\atop 2}\right\rangle = 1: | ||
[123] => (4)[123], [1(4)][23], [12(4)]3 | [123] => (4)[123], [1(4)][23], [12(4)]3 | ||
Строка 46: | Строка 51: | ||
Далее рассмотрим все перестановки порядка <tex dpi = "130">3</tex> с одним подъемом, причем операцией вставки <tex dpi = "130">4</tex> мы будем увеличивать количество перестановок на 1: | Далее рассмотрим все перестановки порядка <tex dpi = "130">3</tex> с одним подъемом, причем операцией вставки <tex dpi = "130">4</tex> мы будем увеличивать количество перестановок на 1: | ||
− | <tex dpi = " | + | :<tex dpi = "140"> \left\langle{3\atop 1}\right\rangle = 4:</tex> |
− | <tex dpi = " | + | :<tex dpi = "140">[13]2 => [13(4)]2, [13][2(4)];</tex> |
− | <tex dpi = " | + | :<tex dpi = "140">2[13] => [2(4)][13], 2[13(4)];</tex> |
− | <tex dpi = " | + | :<tex dpi = "140">[23]1 => [23(4)]1, [23][1(4)];</tex> |
− | <tex dpi = " | + | :<tex dpi = "140">3[12] => [3(4)][12], 3[12(4)];</tex> |
Таким образом мы убеждаемся в верности формулы: | Таким образом мы убеждаемся в верности формулы: | ||
− | <tex dpi = "160"> \left\langle{4\atop 2}\right\rangle = (2 + 1) \left\langle{3\atop 2}\right\rangle + (4 - 2)\left\langle{3\atop 1}\right\rangle = 11;</tex> | + | :<tex dpi = "160"> \left\langle{4\atop 2}\right\rangle = (2 + 1) \left\langle{3\atop 2}\right\rangle + (4 - 2)\left\langle{3\atop 1}\right\rangle = 11;</tex> |
==Треугольник чисел Эйлера I рода и явная формула== | ==Треугольник чисел Эйлера I рода и явная формула== | ||
Строка 66: | Строка 71: | ||
Приведем также без вывода явную формулу для вычисления чисел Эйлера I рода: | Приведем также без вывода явную формулу для вычисления чисел Эйлера I рода: | ||
− | <tex dpi = "180">\left\langle{n\atop m}\right\rangle = \sum_{k=0}^{m}[ (-1)^k {n+1\choose k} (m+1-k)^n]</tex> | + | :<tex dpi = "180">\left\langle{n\atop m}\right\rangle = \sum_{k=0}^{m}[ (-1)^k {n+1\choose k} (m+1-k)^n]</tex> |
Строка 76: | Строка 81: | ||
|- align="center" | |- align="center" | ||
| style="background:white; color:black; width:50px;" | | | style="background:white; color:black; width:50px;" | | ||
− | | style="background:white; color:black; width:50px;" | ''''' | + | | style="background:white; color:black; width:50px;" | '''''m = 0''''' |
| style="background:white; color:black; width:50px;" | '''''1''''' | | style="background:white; color:black; width:50px;" | '''''1''''' | ||
| style="background:white; color:black; width:50px;" | '''''2''''' | | style="background:white; color:black; width:50px;" | '''''2''''' | ||
Строка 288: | Строка 293: | ||
Стоит отметить, что гистрограмма, построенная на значениях чисел Эйлера I рода аппроксимируется к гистограмме, построенной на биноминальных коээфициентах (оба графика, представленные '''справа''', смасштабированы; масштаб указан на гистограмме): | Стоит отметить, что гистрограмма, построенная на значениях чисел Эйлера I рода аппроксимируется к гистограмме, построенной на биноминальных коээфициентах (оба графика, представленные '''справа''', смасштабированы; масштаб указан на гистограмме): | ||
− | [[Файл:Euler_I_hist.gif|300px|thumb | + | [[Файл:Euler_I_hist.gif|300px|thumb|Числа Эйлера I рода (m < 90)]] |
− | [[Файл:Binomial_hist.gif|300px|thumb | + | [[Файл:Binomial_hist.gif|300px|thumb|Биномиальные коээфициенты (m < 60)]] |
'''Полезные факты о числах Эйлера I рода''' | '''Полезные факты о числах Эйлера I рода''' | ||
Строка 299: | Строка 304: | ||
:<tex dpi = "160">\sum_{k=0}^{n} \left\langle{n\atop m}\right\rangle = n!,\ n \ge 0, \,</tex> | :<tex dpi = "160">\sum_{k=0}^{n} \left\langle{n\atop m}\right\rangle = n!,\ n \ge 0, \,</tex> | ||
− | 3. <tex dpi = "160">\sum_{m=0}^n (-1)^m {\left\langle{n\atop m}\right\rangle}{n-1\choose m}^{-1}=0.</tex> | + | 3. Еще одно условие равенства нулю: |
+ | :<tex dpi = "160">\sum_{m=0}^n (-1)^m {\left\langle{n\atop m}\right\rangle}{n-1\choose m}^{-1}=0.</tex> | ||
+ | |||
==Числа Эйлера II рода== | ==Числа Эйлера II рода== | ||
+ | '''''Числа Эйлера II рода''''' — количество перестановок мультимножества от <tex dpi = "130">1</tex> до <tex dpi = "130">n</tex> вида <tex dpi = "130">\{1,1,2,2..n,n\}</tex>, обладающих свойством "все элементы перестановки, встречающиеся между двумя вхождниями <tex dpi = "130">z</tex> для любого <tex dpi = "130">z</tex>, больше, чем <tex dpi = "130">z</tex>", таких, что в каждой из них существует ровно <tex dpi = "130">m</tex> подъемов. Числа Эйлера II рода обозначаются как <tex dpi = "160"> \scriptstyle \left\langle \!\! \left\langle {n \atop m} \right\rangle \!\! \right\rangle </tex> | ||
+ | |||
+ | |||
+ | '''Пример''' | ||
+ | |||
+ | Рассмотрим <tex dpi = "130"> n = 3</tex>. Тогда существует 15 перестановок такого вида, среди которых одна не имеет подъемов, 8 штук имеют всего 1 подъем, и 6 перестановок имеют 2 подъема: | ||
+ | |||
+ | :<tex dpi = "140">112233,\; 122133,\; 112332,\; 123321,\; 133122,\; 122331. </tex> | ||
+ | :<tex dpi = "140"> 221133,\; 221331,\; 223311,\; 233211,\; 113322,\; 133221,\; 331122,\; 331221, </tex> | ||
+ | :<tex dpi = "140"> 332211,\; </tex> | ||
+ | |||
+ | {{Лемма | ||
+ | |statement=Количество перестановок мультимножества <tex dpi = "130">\{1,1,2,2..n,n\}</tex> со свойством "все элементы перестановки, встречающиеся между двумя вхождниями <tex dpi = "130">z</tex> для любого <tex dpi = "130">z</tex>, больше, чем | ||
+ | <tex dpi="130">z</tex>" равно двойному факториалу <tex dpi="130">(2n-1)!!</tex>. | ||
+ | |neat = 1 | ||
+ | }} | ||
+ | |||
+ | |||
+ | ===Рекурсивная формула=== | ||
+ | Числа Эйлера II рода можно выразить рекурсивно следующим образом: | ||
+ | :<tex dpi = "160"> \left\langle \!\! \left\langle {n \atop m} \right\rangle \!\! \right\rangle = (2n-m-1) \left\langle \!\! \left\langle {n-1 \atop m-1} \right\rangle \!\! \right\rangle + (m+1) \left\langle \!\! \left\langle {n-1 \atop m} \right\rangle \!\! \right\rangle, </tex> | ||
+ | |||
+ | С начальным условием для <tex dpi = "130>n = 0</tex>: | ||
+ | :<tex dpi = "160"> \left\langle \!\! \left\langle {0 \atop m} \right\rangle \!\! \right\rangle = [m=0]. </tex> | ||
+ | |||
+ | |||
+ | |||
+ | ===Треугольник чисел Эйлера II рода=== | ||
+ | Значения чисел Эйлера II рода представлены в данном массиве. Нижнедиагональная его часть называется треугольником чисел Эйлера II рода. | ||
+ | |||
+ | ::{| class="number_triangle" | ||
+ | |||
+ | |- align="center" | ||
+ | | style="background:white; color:black; width:50px;" | | ||
+ | | style="background:white; color:black; width:50px;" | '''''m = 0''''' | ||
+ | | style="background:white; color:black; width:50px;" | '''''1''''' | ||
+ | | style="background:white; color:black; width:50px;" | '''''2''''' | ||
+ | | style="background:white; color:black; width:50px;" | '''''3''''' | ||
+ | | style="background:white; color:black; width:50px;" | '''''4''''' | ||
+ | | style="background:white; color:black; width:50px;" | '''''5''''' | ||
+ | | style="background:white; color:black; width:50px;" | '''''6''''' | ||
+ | | style="background:white; color:black; width:50px;" | '''''7''''' | ||
+ | | style="background:white; color:black; width:50px;" | '''''8''''' | ||
+ | |- align="center" | ||
+ | | style="background:white; color:black;" | '''''n = 0''''' | ||
+ | | style="background:#FFDEAD; color:black;" | '''1''' | ||
+ | | style="background:#FFDEAD; color:red;" | '''0''' | ||
+ | | style="background:#FFDEAD; color:red;" | '''0''' | ||
+ | | style="background:#FFDEAD; color:red;" | '''0''' | ||
+ | | style="background:#FFDEAD; color:red;" | '''0''' | ||
+ | | style="background:#FFDEAD; color:red;" | '''0''' | ||
+ | | style="background:#FFDEAD; color:red;" | '''0''' | ||
+ | | style="background:#FFDEAD; color:red;" | '''0''' | ||
+ | | style="background:#FFDEAD; color:red;" | '''0''' | ||
+ | |||
+ | |- align="center" | ||
+ | | style="background:white; color:black;" | '''''1''''' | ||
+ | | style="background:#FFDEAD; color:black;" | '''1''' | ||
+ | | style="background:#FFDEAD; color:red;" | '''0''' | ||
+ | | style="background:#FFDEAD; color:red;" | '''0''' | ||
+ | | style="background:#FFDEAD; color:red;" | '''0''' | ||
+ | | style="background:#FFDEAD; color:red;" | '''0''' | ||
+ | | style="background:#FFDEAD; color:red;" | '''0''' | ||
+ | | style="background:#FFDEAD; color:red;" | '''0''' | ||
+ | | style="background:#FFDEAD; color:red;" | '''0''' | ||
+ | | style="background:#FFDEAD; color:red;" | '''0''' | ||
+ | |- align="center" | ||
+ | | style="background:white; color:black;" | '''''2''''' | ||
+ | | style="background:#FFDEAD; color:black;" | '''1''' | ||
+ | | style="background:#FFDEAD; color:black;" | '''2''' | ||
+ | | style="background:#FFDEAD; color:red;" | '''0''' | ||
+ | | style="background:#FFDEAD; color:red;" | '''0''' | ||
+ | | style="background:#FFDEAD; color:red;" | '''0''' | ||
+ | | style="background:#FFDEAD; color:red;" | '''0''' | ||
+ | | style="background:#FFDEAD; color:red;" | '''0''' | ||
+ | | style="background:#FFDEAD; color:red;" | '''0''' | ||
+ | | style="background:#FFDEAD; color:red;" | '''0''' | ||
+ | |- align="center" | ||
+ | | style="background:white; color:black;" | '''''3''''' | ||
+ | | style="background:#FFDEAD; color:black;" | '''1''' | ||
+ | | style="background:#FFDEAD; color:black;" | '''8''' | ||
+ | | style="background:#FFDEAD; color:black;" | '''6''' | ||
+ | | style="background:#FFDEAD; color:red;" | '''0''' | ||
+ | | style="background:#FFDEAD; color:red;" | '''0''' | ||
+ | | style="background:#FFDEAD; color:red;" | '''0''' | ||
+ | | style="background:#FFDEAD; color:red;" | '''0''' | ||
+ | | style="background:#FFDEAD; color:red;" | '''0''' | ||
+ | | style="background:#FFDEAD; color:red;" | '''0''' | ||
+ | |- align="center" | ||
+ | | style="background:white; color:black;" | '''''4''''' | ||
+ | | style="background:#FFDEAD; color:black;" | '''1''' | ||
+ | | style="background:#FFDEAD; color:black;" | '''22''' | ||
+ | | style="background:#FFDEAD; color:black;" | '''58''' | ||
+ | | style="background:#FFDEAD; color:black;" | '''24''' | ||
+ | | style="background:#FFDEAD; color:red;" | '''0''' | ||
+ | | style="background:#FFDEAD; color:red;" | '''0''' | ||
+ | | style="background:#FFDEAD; color:red;" | '''0''' | ||
+ | | style="background:#FFDEAD; color:red;" | '''0''' | ||
+ | | style="background:#FFDEAD; color:red;" | '''0''' | ||
+ | |- align="center" | ||
+ | | style="background:white; color:black;" | '''''5''''' | ||
+ | | style="background:#FFDEAD; color:black;" | '''1''' | ||
+ | | style="background:#FFDEAD; color:black;" | '''52''' | ||
+ | | style="background:#FFDEAD; color:black;" | '''328''' | ||
+ | | style="background:#FFDEAD; color:black;" | '''444''' | ||
+ | | style="background:#FFDEAD; color:black;" | '''120''' | ||
+ | | style="background:#FFDEAD; color:red;" | '''0''' | ||
+ | | style="background:#FFDEAD; color:red;" | '''0''' | ||
+ | | style="background:#FFDEAD; color:red;" | '''0''' | ||
+ | | style="background:#FFDEAD; color:red;" | '''0''' | ||
+ | |- align="center" | ||
+ | | style="background:white; color:black;" | '''''6''''' | ||
+ | | style="background:#FFDEAD; color:black;" | '''1''' | ||
+ | | style="background:#FFDEAD; color:black;" | '''114''' | ||
+ | | style="background:#FFDEAD; color:black;" | '''1452''' | ||
+ | | style="background:#FFDEAD; color:black;" | '''4400''' | ||
+ | | style="background:#FFDEAD; color:black;" | '''3708''' | ||
+ | | style="background:#FFDEAD; color:black;" | '''720''' | ||
+ | | style="background:#FFDEAD; color:red;" | '''0''' | ||
+ | | style="background:#FFDEAD; color:red;" | '''0''' | ||
+ | | style="background:#FFDEAD; color:red;" | '''0''' | ||
+ | |- align="center" | ||
+ | | style="background:white; color:black;" | '''''7''''' | ||
+ | | style="background:#FFDEAD; color:black;" | '''1''' | ||
+ | | style="background:#FFDEAD; color:black;" | '''240''' | ||
+ | | style="background:#FFDEAD; color:black;" | '''5610''' | ||
+ | | style="background:#FFDEAD; color:black;" | '''32120''' | ||
+ | | style="background:#FFDEAD; color:black;" | '''58140''' | ||
+ | | style="background:#FFDEAD; color:black;" | '''33984''' | ||
+ | | style="background:#FFDEAD; color:black;" | '''5040''' | ||
+ | | style="background:#FFDEAD; color:red;" | '''0''' | ||
+ | | style="background:#FFDEAD; color:red;" | '''0''' | ||
+ | |- align="center" | ||
+ | | style="background:white; color:black;" | '''''8''''' | ||
+ | | style="background:#FFDEAD; color:black;" | '''1''' | ||
+ | | style="background:#FFDEAD; color:black;" | '''494''' | ||
+ | | style="background:#FFDEAD; color:black;" | '''19950''' | ||
+ | | style="background:#FFDEAD; color:black;" | '''195800''' | ||
+ | | style="background:#FFDEAD; color:black;" | '''644020''' | ||
+ | | style="background:#FFDEAD; color:black;" | '''785304''' | ||
+ | | style="background:#FFDEAD; color:black;" | '''341136''' | ||
+ | | style="background:#FFDEAD; color:black;" | '''40320''' | ||
+ | | style="background:#FFDEAD; color:red;" | '''0''' | ||
+ | |- align="center" | ||
+ | | style="background:white; color:black;" | '''''9''''' | ||
+ | | style="background:#FFDEAD; color:black;" | '''1''' | ||
+ | | style="background:#FFDEAD; color:black;" | '''1004''' | ||
+ | | style="background:#FFDEAD; color:black;" | '''67260''' | ||
+ | | style="background:#FFDEAD; color:black;" | '''1062500''' | ||
+ | | style="background:#FFDEAD; color:black;" | '''5765500''' | ||
+ | | style="background:#FFDEAD; color:black;" | '''12440064''' | ||
+ | | style="background:#FFDEAD; color:black;" | '''11026296''' | ||
+ | | style="background:#FFDEAD; color:black;" | '''3733920''' | ||
+ | | style="background:#FFDEAD; color:black;" | '''362880''' | ||
+ | |} | ||
==Ссылки== | ==Ссылки== |
Версия 01:51, 19 декабря 2013
Числа Эйлера I рода (Eulerian numbers) — количество перестановок чисел от 1 до n таких, что в каждой из них существует ровно m подъемов. Числа Эйлера I рода обозначают как или же .
Определение: |
Пусть | и - элементы некоторой перестановки порядка причем . Тогда пара называется подъемом (ascent) данной перестановки.
Содержание
Вывод рекуррентной формулы
Пусть у нас есть некая перестановка
. Тогда операцией вставки элемента с номером n в какую-либо из позиций мы получим перестановок вида . Далее рассмотрим два случая:1. Количество подъемов в перестановке
равно количеству подъемов в . Этого можно добиться, вставляя элемент на самое первое место в (всего возможностей) или перед последним последним элементом каждого подъема(еще раз).2. Количество подъемов в новой перестановке на один больше предыдущего количества. Этого эффекта добиваемся вставкой элемента
в конце каждой перестановки или после элемента перестановки со значением . Таких элементов, как не трудно догадаться, будет .Тогда рекуррентная формула имеет вид:
Примем также следующие начальные значения:
- , если
Также для четных
:- ,
Запись [выражение] означает нотацию Айверсона, где
Пример
Рассмотрим все перестановки порядка
, в которых есть ровно 2 подъема (в квадратных скобках один или больше подъемов подряд):Согласно алгоритму вывода рекуррентной формулы мы можем добавить '4' в следующие позиции всех перестановок порядка
с двумя подъемами, не увеличив количество подъемов:Далее рассмотрим все перестановки порядка
с одним подъемом, причем операцией вставки мы будем увеличивать количество перестановок на 1:Таким образом мы убеждаемся в верности формулы:
Треугольник чисел Эйлера I рода и явная формула
Явная формула
Приведем также без вывода явную формулу для вычисления чисел Эйлера I рода:
Треугольник чисел Эйлера I рода
На значениях
чисел Эйлера I рода можно построить массив , нижнедиагональная часть которого названа треугольником чисел Эйлера I рода.m = 0 1 2 3 4 5 6 7 8 9 10 11 12 n = 0 1 0 0 0 0 0 0 0 0 0 0 0 0 1 1 0 0 0 0 0 0 0 0 0 0 0 0 2 1 1 0 0 0 0 0 0 0 0 0 0 0 3 1 4 1 0 0 0 0 0 0 0 0 0 0 4 1 11 11 1 0 0 0 0 0 0 0 0 0 5 1 26 66 26 1 0 0 0 0 0 0 0 0 6 1 57 302 302 57 1 0 0 0 0 0 0 0 7 1 120 1191 2416 1191 120 1 0 0 0 0 0 0 8 1 247 4293 15619 15619 4293 247 1 0 0 0 0 0 9 1 502 14608 88234 156190 88234 14608 502 1 0 0 0 0 10 1 1013 47840 455192 1310354 1310354 455192 47840 1013 1 0 0 0 11 1 2036 152637 2203488 9738114 15724248 9738114 2203488 152637 2036 1 0 0 12 1 4083 478271 10187685 66318474 162512286 162512286 66318474 10187685 478271 4083 1 0
Стоит отметить, что гистрограмма, построенная на значениях чисел Эйлера I рода аппроксимируется к гистограмме, построенной на биноминальных коээфициентах (оба графика, представленные справа, смасштабированы; масштаб указан на гистограмме):
Полезные факты о числах Эйлера I рода
1. Нетрудно увидеть, что каждый ряд ненулевых значений симметричен относительно своей середины, то есть:
2. Сумма всех значений каждого ряда равна
:3. Еще одно условие равенства нулю:
Числа Эйлера II рода
Числа Эйлера II рода — количество перестановок мультимножества от
до вида , обладающих свойством "все элементы перестановки, встречающиеся между двумя вхождниями для любого , больше, чем ", таких, что в каждой из них существует ровно подъемов. Числа Эйлера II рода обозначаются как
Пример
Рассмотрим
. Тогда существует 15 перестановок такого вида, среди которых одна не имеет подъемов, 8 штук имеют всего 1 подъем, и 6 перестановок имеют 2 подъема:Лемма: |
Количество перестановок мультимножества со свойством "все элементы перестановки, встречающиеся между двумя вхождниями для любого , больше, чем
" равно двойному факториалу . |
Рекурсивная формула
Числа Эйлера II рода можно выразить рекурсивно следующим образом:
С начальным условием для
:
Треугольник чисел Эйлера II рода
Значения чисел Эйлера II рода представлены в данном массиве. Нижнедиагональная его часть называется треугольником чисел Эйлера II рода.
m = 0 1 2 3 4 5 6 7 8 n = 0 1 0 0 0 0 0 0 0 0 1 1 0 0 0 0 0 0 0 0 2 1 2 0 0 0 0 0 0 0 3 1 8 6 0 0 0 0 0 0 4 1 22 58 24 0 0 0 0 0 5 1 52 328 444 120 0 0 0 0 6 1 114 1452 4400 3708 720 0 0 0 7 1 240 5610 32120 58140 33984 5040 0 0 8 1 494 19950 195800 644020 785304 341136 40320 0 9 1 1004 67260 1062500 5765500 12440064 11026296 3733920 362880