Поліноми Белла

Матеріал з testwiki
Версія від 19:18, 15 листопада 2023, створена imported>Weqwtur (growthexperiments-addlink-summary-summary:3|0|0)
(різн.) ← Попередня версія | Поточна версія (різн.) | Новіша версія → (різн.)
Перейти до навігації Перейти до пошуку

У комбінаториці поліноми Белла, що названі на честь Еріка Темпла Белла, використовуються для вивчення заданих розділів. Вони пов'язані з числами Стірлінга та Белла. Вони також зустрічаються у багатьох програмах, наприклад у формулі Фаа ді Бруно.

Поліноми Белла

Експоненційні поліноми Белла

Часткові або неповні експоненційні поліноми Белла — це трикутний масив поліномів, заданий

Bn,k(x1,x2,,xnk+1)=n!j1!j2!jnk+1!(x11!)j1(x22!)j2(xnk+1(nk+1)!)jnk+1,

де сума береться за всіма послідовностями j1, j2, j3, ..., jnk +1 невід’ємних цілих чисел, таким чином, що б виконувалися наступні дві умови :

j1+j2++jnk+1=k,
j1+2j2+3j3++(nk+1)jnk+1=n.

Сума

Bn(x1,,xn)=k=1nBn,k(x1,x2,,xnk+1)

називається n-м повним експоненційними полінонмоми Белла.

Звичайні поліноми Белла

Так само часткові звичайні поліноми Белла, на відміну від звичайного експоненціального многочлена Белла, визначений вище, задається формулою

B^n,k(x1,x2,,xnk+1)=k!j1!j2!jnk+1!x1j1x2j2xnk+1jnk+1,

де сума пробігає всі послідовності j1, j2, j3, ..., jnk +1 невід'ємних цілих чисел, таких що

j1+j2++jnk+1=k,
j1+2j2++(nk+1)jnk+1=n.

Звичайні поліноми Белла можна виразити в термінах експоненційних поліномів Белла:

B^n,k(x1,x2,,xnk+1)=k!n!Bn,k(1!x1,2!x2,,(nk+1)!xnk+1).

Як правило, під поліномами Белла маються на увазі експоненційні поліноми Белла, якщо не зазначено іншого.

Комбінаторний сенс

Експоненційні поліноми Белла безпосередньо стосується способів розбиття множин. Наприклад, якщо ми розглядаємо множину {A, B, C}, її можна розбити на два непусті, підмножини що не перетинаються, які також називають частинами або блоками, 3 різними способами:

{{A}, {B, C}}
{{B}, {A, C}}
{{C}, {B, A}}

Таким чином, ми можемо закодувати інформацію щодо цих розбиттів як

B3,2(x1,x2)=3x1x2.

Тут підписники B3,2 говорять нам, що ми розглядаємо поділ набору з 3-х елементів на 2 блоки. Підрядник кожного xi вказує на наявність блоку з елементами j (або блоку розміру i) в заданому розділі. Отже, x2 вказує на наявність блоку з двома елементами. Аналогічно, x1 вказує на наявність блоку з одним елементом. Експонент xij вказує на те, що в одному розбитті є j таких блоків розміром i. Тут, оскільки і x1, і x2 є експонент 1, це вказує, що в даному розділі є лише один такий блок. Коефіцієнт одночлена вказує, скільки таких перегородок є. У нашому випадку є 3 розділи набору з 3 елементами на 2 блоки, де в кожній секції елементи розділені на два блоки розмірами 1 і 2.

Оскільки будь-який набір може бути розділений на один блок лише одним способом, вищезгадане тлумачення означало б, що Bn,1 = xn. Аналогічно, оскільки існує лише один спосіб поділу множини з n елементами на n одиниць, Bn,n = x1n.

Як складніший приклад розглянемо

B6,2(x1,x2,x3,x4,x5)=6x5x1+15x4x2+10x32

Це говорить нам про те, що якщо набір з 6 елементами розділений на 2 блоки, то ми можемо мати 6 розділів з блоками розміру 1 і 5, 15 розділів з блоками розміру 4 і 2, і 10 розділів з 2 блоками розміром 3.

Сума підписок у монометах дорівнює загальній кількості елементів. Таким чином, кількість одночленів, що з'являються в частковому многочлені Белла, дорівнює кількості способів, що ціле число n може бути виражене як підсумок k додатних цілих чисел. Це і є розбиття числа n на k частин. Наприклад, у вище наведених прикладах ціле число 3 можна розділити на дві частини лише як 2 + 1. Таким чином, у B3,2 містить лише один одночлен. Однак ціле число 6 можна розділити на дві частини як 5 + 1, 4 + 2 і 3 + 3. Таким чином, B6,2 містить три одночлени. Дійсно, індекси змінних у мономери такі самі, як ті, які задаються цілим розділом, що вказує на розміри різних блоків. Таким чином, загальна кількість одночленів, що з'являються в повному поліномі Белла Bn, дорівнює загальній кількості розбиттів числа n.

Також ступінь кожного одночлена, що є сумою експонентів кожної змінної в мономі, дорівнює кількості блоків, на які поділяється множина. Тобто j 1 + j 2 + ... = k. Таким чином, задавши повний многочлен Белла B n, ми можемо відокремити частковий многочлен Белла Bn,k, зібравши всі ці мономи зі ступенем k.

Нарешті, якщо знехтувати розмірами блоків і поставити всі xi = x, то підсумовування коефіцієнтів часткового многочлена Белла Bn, k дасть загальну кількість способів розподілу множини з n елементів k блоки, що те саме, що числа Стірлінга другого роду. Крім того, підсумовування всіх коефіцієнтів повного многочлена Белла Bn дасть нам загальну кількість способів поділити набір з п елементами на підмножини, що не перекриваються, що те саме, що і число Белла.

Загалом, якщо ціле число n розбиттів на суму, в якій «1» з'являється j1 раз, «2» з'являється j2 рази і так далі, то кількість розбиття множини розміром n, які згортаються до цього розділу цілого числа n, коли члени множини стають невідрізними, — це відповідний коефіцієнт у многочлени.

Приклади

Нехай, ми маємо

B6,2(x1,x2,x3,x4,x5)=6x5x1+15x4x2+10x32

оскільки є

6 способів розбити множину 6 як 5   +   1,
15 способів розбити множину 6 як 4   +   2, і
10 способів розбити множину 6 як 3   +   3.

Подібно

B6,3(x1,x2,x3,x4)=15x4x12+60x3x2x1+15x23

оскільки є

15 способів розбити множину 6 на 4   +   1   +   1,
60 способів розділити множину 6 як 3   +   2   +   1, і
15 способів розбити множину 6 як 2   +   2   +   2.

Властивості

Генератриса

Експоненційні часткові поліноми Белла можна визначити подвійним розкладом у ряд його генератриси:

Φ(t,u)=exp(uj=1xjtjj!)=nk0Bn,k(x1,,xnk+1)tnn!uk=1+n=1tnn!k=1nukBn,k(x1,,xnk+1).

Інакше кажучи, що є фактично те ж саме, шляхом розкладу у ряд k-го степеню:

1k!(j=1xjtjj!)k=n=kBn,k(x1,,xnk+1)tnn!,k=0,1,2,

Повні експоненційні поліноми Белла визначається за допомогою Φ(t,1) або інакше:

Φ(t,1)=exp(j=1xjtjj!)=n=0Bn(x1,,xn)tnn!.

Таким чином, n -й повні поліноми Белла задається як

Bn(x1,,xn)=(t)nexp(j=1xjtjj!)|t=0.

Так само звичайні часткові поліноми Белла можна визначити твірною функцією

Φ^(t,u)=exp(uj=1xjtj)=nk0B^n,k(x1,,xnk+1)tnukk!.

Або, що еквівалентно, розкладом у ряд k-ї степеня:

(j=1xjtj)k=n=kB^n,k(x1,,xnk+1)tn.

Дивись також перетворення твірної функції для розкладу у ряд твірної функцій поліномів Белла що є композицією послідовностей твірних функцій та степенів, логаритмів та експоненційної функції, що послідовності твірних функцій. Кожна з цих формул можна знайти у відповідних розділах праці КомтетаШаблон:Sfn.

Рекурентні співвідношення

Повні поліноми Белла можна визначити за допомогою рекурентних співвідношень як

Bn+1(x1,,xn+1)=i=0n(ni)Bni(x1,,xni)xi+1

з початковим значенням B0=1 .

Часткові поліноми Белла також можна ефективно обчислені рекурентним співвідношенням:

Bn,k=i=1nk+1(n1i1)xiBni,k1,

де

B0,0=1;
Bn,0=0 for n1;

B0,k=0 for k1.

Повні поліноми Белла також задовольняють такій диференціальній рекурентній формуліШаблон:Sfn:

Bn(x1,,xn)=1n1[i=2nj=1i1(i1)(i2j1)xjxijBn1(x1,,xn1)xi1+i=2nj=1i1xi+1(ij)2Bn1(x1,,xn1)xjxij+i=2nxiBn1(x1,,xn1)xi1].

У формі визначника

Повні поліноми Белла можна виразити у вигляді визначника:

Bn(x1,,xn)=det[x1(n11)x2(n12)x3(n13)x4(n14)x5xn1x1(n21)x2(n22)x3(n23)x4xn101x1(n31)x2(n32)x3xn2001x1(n41)x2xn30001x1xn400001xn5000001x1].

Числа Стірлінга та Белла

Значення поліномів Белла Bn,k (x1, x2, ...) для набору факторіалів дорівнює числу Стірлінга першого роду (без знаку):

Bn,k(0!,1!,,(nk)!)=c(n,k)=|s(n,k)|=[nk].

Поліноми Белла Bn, k (x1, x2, ...) від набору одиниць дорівнює числам Стірлінга другого роду :

Bn,k(1,1,,1)=S(n,k)={nk}.

Сума цих значень дає значення повного полінома Белла від набору одиниць:

Bn(1,1,,1)=k=1nBn,k(1,1,,1)=k=1n{nk},

що є n-м числом Белла.

Зворотні співвідношення

Якщо ми визначимо

xn=k=1n(1)k1(k1)!Bn,k(y1,,ynk+1).

то ми маємо таке зворотне співвідношення

(xy)n=j=1n1(nj)xjynj.

Поліноми Тушара

Поліноми Тушара Tn(x)=k=0n{nk}xk може бути виражена як значення повного многочлена Белла від аргументів, що всі дорівнюють х:

Tn(x)=Bn(x,x,,x).

Тотожність згортки

Для послідовностей x n, y n, n = 1, 2, ..., визначте згортку по:

(xy)n=j=1n1(nj)xjynj.

Межі підсумовування — 1 і n  — 1, а не 0 і n.

Дозволяти xnkn-й член послідовності

xxk factors.

Тоді [1]

Bn,k(x1,,xnk+1)=xnkk!.

Наприклад, давайте обчислити B4,3(x1,x2). Ми маємо

x=(x1 , x2 , x3 , x4 ,)

xx=(0, 2x12 , 6x1x2 8x1x3+6x22)
xxx=(0 , 0 , 6x13 36x12x2)

і, таким чином,

B4,3(x1,x2)=(xxx)43!=6x12x2.

Інші тотожности

  • Bn,k(1!,2!,,(nk+1)!)=(n1k1)n!k!=L(n,k) що надає число Лаха.
  • Bn,k(1,2,3,,nk+1)=(nk)knk що повертає ідемпотентне число.
  • Повний многочлен Белла задовольняє відношення типу двочлена:
    Bn(x1+y1,,xn+yn)=i=0n(ni)Bni(x1,,xni)Bi(y1,,yi).
    Bn,k(xq+1(q+1q),xq+2(q+2q),)=n!(q!)k(n+qk)!Bn+qk,k(,0,0,xq+1,xq+2,)

Це виправляє опущення фактора (q!)k у книзі КонтаШаблон:Sfn.

  • Коли 1a<n,
Bn,na(x1,,xa+1)=j=a+12aj!a!(nj)(x1)njBa,ja(x22,x33,,x2(a+1)j2(a+1)j).
  • Особливі випадки часткових многочленів Белла:
Bn,1(x1,,xn)=xn[8pt]Bn,2(x1,,xn1)=12k=1n1(nk)xkxnk[8pt]Bn,n(x1)=(x1)n[8pt]Bn,n1(x1,x2)=(n2)(x1)n2x2[8pt]Bn,n2(x1,x2,x3)=(n3)(x1)n3x3+3(n4)(x1)n4(x2)2[8pt]Bn,n3(x1,x2,x3,x4)=(n4)(x1)n4x4+10(n5)(x1)n5x2x3+15(n6)(x1)n6(x2)3[8pt]Bn,n4(x1,x2,x3,x4,x5)=(n5)(x1)n5x5+5(n6)(x1)n6[3x2x4+2(x3)2]+105(n7)(x1)n7(x2)2x3+105(n8)(x1)n8(x2)4.

Приклади

Перші декілька повних многочленів Белла:

B0=1,B1(x1)=x1,B2(x1,x2)=x12+x2,B3(x1,x2,x3)=x13+3x1x2+x3,B4(x1,x2,x3,x4)=x14+6x12x2+4x1x3+3x22+x4,B5(x1,x2,x3,x4,x5)=x15+10x2x13+15x22x1+10x3x12+10x3x2+5x4x1+x5B6(x1,x2,x3,x4,x5,x6)=x16+15x2x14+20x3x13+45x22x12+15x23+60x3x2x1+15x4x12+10x32+15x4x2+6x5x1+x6,B7(x1,x2,x3,x4,x5,x6,x7)=x17+21x15x2+35x14x3+105x13x22+35x13x4+210x12x2x3+105x1x23+21x12x5+105x1x2x4+70x1x32+105x22x3+7x1x6+21x2x5+35x3x4+x7.

Застосування

Формула Фаа ді Бруно може бути викладена з точки зору поліномів Белла наступним чином:

dndxnf(g(x))=k=1nf(k)(g(x))Bn,k(g(x),g(x),,g(nk+1)(x)).

Аналогічно, версію формули Фаа ді Бруно можна подати з використанням многочленів Белла наступним чином. Припустимо


f(x)=n=1ann!xnandg(x)=n=1bnn!xn.

Тоді


g(f(x))=n=1k=1nbkBn,k(a1,,ank+1)n!xn.

Зокрема, повні поліноми Белла фігурують у розкладі експоненти формальний степеневий ряд:

exp(i=1aii!xi)=n=0Bn(a1,,an)n!xn,

що також представляє генератису для експоненти повних многочленів Белла на фіксованій послідовності аргументів a1,a2, .

Обернення ряду

Нехай дві функції f і g виражаються у формальних рядах потужностей як

f(w)=k=0fkwkk!,andg(z)=k=0gkzkk!,

такий, що g — композиційна інверсія f, визначена g(f (w)) = w або f (g(z)) = z. Якщо f0 = 0 і f1 ≠ 0, то явна форма коефіцієнтів оберненої може бути задана в терміні многочленів Белла якШаблон:Sfn

Z(Sn)=Bn(0!a1,1!a2,,(n1)!an)n!.

з f^k=fk+1(k+1)f1, і nk¯=n(n+1)(n+k1) — це фактор, що зростає, і g1=1f1.

Симетричні многочлени

Елементарний симетричний многочлен en і степеневий многочлен суми потужності pn можуть бути пов'язані один з одним за допомогою поліномів Белла як:


en=(1)nn!Bn(p1,1!p2,2!p3,3!p4,,(1)n1(n1)!pn),pn=(1)n1(n1)!k=1n(1)k1(k1)!Bn,k(e1,2!e2,3!e3,,(nk+1)!enk+1)=(1)nnk=1n1kB^n,k(e1,,enk+1).

Ці формули дозволяють виразити коефіцієнти монічних поліноми через поліноми Белла з нульовими аргументами. Наприклад, разом із теоремою Кейлі — Гамільтона вони призводять до вираження детермінантної n × n квадратної матриці A через сліди її потужностей:

det(A)=(1)nn!Bn(s1,s2,,sn),where sk=(1)k1(k1)!tr(Ak).

Індекс циклу симетричних груп

Індекс циклу симетричної групи Sn можна виразити через повні многочлени Белла так:

Z(Sn)=Bn(0!a1,1!a2,,(n1)!an)n!.

Поліноми Ерміта

Поліноми Ерміта імовірністів можна виразити через поліноми Белла як

Hen(x)=Bn(x,1,0,,0),

де xi = 0 для всіх i > 2; що передбачає комбінаторну інтерпретацію коефіцієнтів поліномів Ерміта. Це можна побачити, порівнюючи генератрису поліномів Ерміта

exp(xtt22)=n=0Hen(x)tnn!

з поліномами Белла.

Програмне забезпечення

Поліноми Белла реалізовані в:

Дивитися також

Примітки

Шаблон:Reflist

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