En mathématiques, et plus précisément en combinatoire, les polynômes de Bell, nommés ainsi d'après le mathématicien Eric Temple Bell, sont des polynômes multivariés définis par :

(n,k)∈ℕ2,Bn,k(x1,x2,…)=defn!∑(m1,m2,…)∈ℕℕ∗{∑i∈ℕ∗mi=k∑i∈ℕ∗i⋅mi=n∏i∈ℕ∗ximimi!(i!)mi

Sachant que mi est forcément nul pour i > n − k + 1, on peut expliciter la borne supérieure des indices i :

Bn,k(x1,…,xn−k+1)=n!∑(m1,…,mn−k+1)∈ℕn−k+1{∑i=1n−k+1mi=k∑i=1n−k+1i⋅mi=n∏i=1n−k+1ximimi!(i!)mi

Interprétation combinatoire

modifier

Soit un ensemble de n éléments partitionné en k sous-ensembles non vides, dont m1 sous-ensembles de cardinalité 1, m2 sous-ensembles de cardinalité 2, etc. Le nombre de telles partitions est le coefficient du monôme unitaire xm1
1
xm2
2
…
dans le polynôme de Bell Bn,k(x1, x2, …). On notera que :

  • par construction, on a ∑ mi = k (nombre de sous-ensembles) et ∑ i·mi = n (nombre total d'éléments), avec chaque mi positif ou nul (nombre de sous-ensembles de cardinalité i) ;
  • par conséquent, les cardinalités des sous-ensembles forment une partition de l'entier n en k parties, avec mi la multiplicité de l'entier i dans cette partition.

Exemples

modifier

On a :

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

car il y a :

  • 6 partitions d'un ensemble à 6 éléments de la forme 5 + 1 ;
  • 15 partitions de la forme 4 + 2 ;
  • 10 partitions de la forme 3 + 3.

De même :

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

car il y a :

  • 15 partitions d'un ensemble à 6 éléments de la forme 4 + 1 + 1 ;
  • 60 partitions de la forme 3 + 2 + 1 ;
  • 15 partitions de la forme 2 + 2 + 2.

Polynômes de Bell complets

modifier

La somme

Bn(x1,…,xn)=∑k=0nBn,k(x1,…,xn−k+1)

est parfois appelée n-ème polynôme de Bell complet, et alors les polynômes Bn, k définis ci-dessus sont appelés des polynômes de Bell « partiels ». Les polynômes de Bell complets Bn peuvent être exprimés par le déterminant d’une matrice :

Bn(x1,…,xn)=|(00)x1(10)x2(20)x3(30)x4⋯(n−20)xn−1(n−10)xn−1(11)x1(21)x2(31)x3⋯(n−21)xn−2(n−11)xn−10−1(22)x1(32)x2⋯(n−22)xn−3(n−12)xn−200−1(33)x1⋯(n−23)xn−4(n−13)xn−3⋮⋮⋮⋱⋱⋮⋮0000⋯(n−2n−2)x1(n−1n−2)x20000⋯−1(n−1n−1)x1|=|(j−1i−1)xj−i+1−δj−i+1|

avec δk le symbole de Kronecker. La matrice dont Bn est le déterminant est une matrice de Hessenberg.

Table de valeurs

modifier

Le tableau suivant regroupe les premières valeurs de Bn,k(x1,…,xn−k+1) :

k
n
0 1 2 3 4 5 6
0 1
1 0 x1
2 0 x2 x12
3 0 x3 3x2x1 x13
4 0 x4 4x3x1+3x22 6x2x12 x14
5 0 x5 5x4x1+10x3x2 10x3x12+15x22x1 10x2x13 x15
6 0 x6 6x5x1+15x4x2+10x32 15x4x12+60x3x2x1+15x23 20x3x13+45x22x12 15x2x14 x16

Propriétés

modifier

Cas limites

modifier

Séries formelles exponentielles

modifier
Puissance
1k!(∑n=1∞xntnn!)k=∑n=k∞Bn,k(x1,…,xn−k+1)tnn!
Formule exponentielle
exp⁡(∑n=1∞xntnn!)=∑n=0∞Bn(x1,…,xn)tnn!
Composition

Si

f(t)=∑n=0∞fntnn! (avec f0=f(0) quelconque)

et

g(t)=∑n=1∞gntnn! (avec donc g0=g(0)=0)

alors

(f∘g)(t)=∑n=0∞[∑k=0nfkBn,k(g1,…,gn−k+1)]tnn!

Formules de récurrence

modifier
Bn+1,k+1(x1,…,xn−k+1)=∑i=0n−k(ni)xi+1Bn−i,k(x1,…,xn−i−k+1)=∑i=kn(ni)xn−i+1Bi,k(x1,…,xi−k+1)

avec Bn,0 = δn.

Bn+1(x1,…,xn+1)=∑i=0n(ni)xi+1Bn−i(x1,…,xn−i)=∑i=0n(ni)xn−i+1Bi(x1,…,xi)

avec B0 = B0,0 = 1.

Valeurs particulières

modifier
  • Bn,k(1,1,…,1)={nk} (nombre de Stirling de seconde espèce non signé)
  • Bn,k(1,−1,1,−1,…,(−1)n−k)=(−1)n−k{nk} (nombre de Stirling de seconde espèce signé)
  • Bn(1,1,…,1)=Bn (nombre de Bell)
  • Bn,k(0!,1!,…,(n−k)!)=[nk] (nombre de Stirling de première espèce non signé)
  • Bn,k(0!,−1!,2!,−3!,…,(−1)n−k(n−k)!)=(−1)n−k[nk] (nombre de Stirling de première espèce signé)
  • Bn(0!,1!,…,(n−1)!)=n!
  • Bn(0!,−1!,2!,−3!,…,(−1)n−1(n−1)!)=δn,0+δn,1
  • Bn(−0!,−1!,…,−(n−1)!)=δn,0−δn,1
  • Bn,k(1!,2!,…,(n−k+1)!)=⌊nk⌋ (nombre de Lah non signé)
  • Bn,k(−1!,2!,−3!,4!,…,(−1)n−k+1(n−k+1)!)=(−1)n⌊nk⌋ (nombre de Lah signé)

Type binomial

modifier
Bn(x1+y1,…,xn+yn)=∑k=0n(nk)Bn−k(x1,…,xn−k)Bk(y1,…,yk)

avec B0 = 1.

Réciproque

modifier

Soit f une fonction infiniment dérivable en un point a et de réciproque f -1, alors[1] :

yn=∑k=1nf(k)(a)Bn,k(x1,…,xn−k+1)⇔xn=∑k=1n[f−1](k)(f(a))Bn,k(y1,…,yn−k+1)

Cas particuliers

modifier

En prenant f (x) = ex (soit f –1(x) = ln(x)) infiniment dérivable en 0, on a :

  • f(k)(0)=1
  • [f−1](k)(f(0))=(−1)k−1(k−1)!

d’où :

yn=∑k=1nBn,k(x1,…,xn−k+1)⇔xn=∑k=1n(−1)k−1(k−1)!Bn,k(y1,…,yn−k+1)

soit :

xn=∑k=1n(−1)k−1(k−1)!Bn,k[B1(x1),…,Bn−k+1(x1,…,xn−k+1)]

En prenant f (x) = xα avec α ≠ 0 (soit f –1(x) = x1/α) infiniment dérivable en 1, on a :

  • f(k)(1)=αk_
  • [f−1](k)(f(1))=(1α)k_

avec .k la factorielle décroissante, d’où :

yn=∑k=1nαk_Bn,k(x1,…,xn−k+1)⇔xn=∑k=1n(1α)k_Bn,k(y1,…,yn−k+1)

Composition

modifier

Soient :

f(t)=∑n=0∞fntnn!
g(t)=∑n=1∞gntnn! (on note que g0=g(0)=0)
(f∘g)(t)=h(t)=∑n=0∞hntnn!

Alors :

Bn,k(h1,…,hn−k+1)=∑i=knBi,k(f1,…,fi−k+1)Bn,i(g1,…,gn−i+1)

En posant les matrices F=(Bj−1,i−1(f1,…,fj−i+1)) (triangulaire supérieure) ainsi que G et H de manière similaire, on a alors :

H=FG

Cas particuliers

modifier

En prenant f(t)=exp⁡(t) et g(t)=ln⁡(1+t), on obtient :

∑i=kn(−1)n−i{ik}[ni]=δn,k

En prenant f(t)=ln⁡(1+t) et g(t)=exp⁡(t)−1, on obtient :

∑i=kn(−1)i−k[ik]{ni}=δn,k

En prenant f(t)=exp⁡(t) et g(t)=−ln⁡(1−t), on obtient :

∑i=kn{ik}[ni]=⌊nk⌋

En prenant f(t)=11−t et g(t)=t1+t, on obtient :

∑i=kn(−1)n−i⌊ik⌋⌊ni⌋=δn,k

Factorielle décroissante

modifier
∀(a,b)∈ℕ2,∑k=1nak_Bn,k(b1_,…,bn_)=(ab)k_[2]

avec .k la factorielle décroissante.

Comportement d’échelle

modifier

Polynômes de Bell partiels

modifier
Cas général
Bn,k(αβx1,…,αβn−k+1xn−k+1)=αkβnBn,k(x1,…,xn−k+1)
Cas particuliers
Bn,k(αx1,…,αxn−k+1)=αkBn,k(x1,…,xn−k+1)
Bn,k(αx1,…,αn−k+1xn−k+1)=αnBn,k(x1,…,xn−k+1)

Polynômes de Bell complets

modifier
Cas général
Bn(αβx1,…,αβnxn)=βn∑k=0nαkBn,k(x1,…,xn−k+1)
Cas particuliers
Bn(αx1,…,αxn)=∑k=0nαkBn,k(x1,…,xn−k+1)
Bn(αx1,…,αnxn)=αnBn(x1,…,xn)
Autre expression
Bn(αx1,…,αxn)=∑k=0nαk_Bn,k[B1(x1),…,Bn−k+1(x1,x2,…,xn−k+1)]

avec .k la factorielle décroissante.

Identité de convolution

modifier

Pour des suites xn, yn, n = 1, 2, …, on peut définir un produit de convolution par :

(x♢y)n=∑j=1n−1(nj)xjyn−j

(les bornes de sommation étant 1 et n − 1, et non 0 et n).

Soit xnk♢ le n-ème terme de la suite x♢⋯♢x⏟k facteurs

Alors :

Bn,k(x1,…,xn−k+1)=xnk♢k!

Applications

modifier

Formule de Faà di Bruno

modifier

La formule de Faà di Bruno peut être énoncée à l'aide des polynômes de Bell de la manière suivante :

(f∘g)(n)(x)=∑k=0n(f(k)∘g)(x)Bn,k(g′(x),…,g(n−k+1)(x))

Moments et cumulants

modifier

Pour une variable aléatoire réelle dont le moment ordinaire mr d’ordre r existe, on a :

mr=Br(κ1,…,κr)

avec κi les cumulants.

Représentations de suites polynomiales

modifier

Pour toute suite a1, a2, … de scalaires, soit :

pn(x)=∑k=1nBn,k(a1,…,an−k+1)xk

Cette suite de polynômes est de type binomial, c'est-à-dire qu'elle satisfait l'identité binomiale suivante :

pn(x+y)=∑k=0n(nk)pk(x)pn−k(y)

pour n ≥ 0.

En fait, on a également la réciproque :

Théorème — Toutes les suites de polynômes de type binomial peuvent s’exprimer sous la forme faisant intervenir les polynômes de Bell.

Si nous posons

h(x)=∑n=1∞anxnn!

en considérant cette série comme une série formelle, alors pour tout n :

h−1(ddx)pn(x)=npn−1(x)

Notes et références

modifier
  1. ↑ (en) W.-S. Chaou, Leetsch C. Hsu, Peter J.-S. Shiue, “Application of Faà di Bruno’s formula in characterization of inverse relations”, dans Journal of Computational and Applied Mathematics, vol. 190, 2006, p. 151–169
  2. ↑ (en) Andrzej Korzeniowski, “Binomial Tails Domination for Random Graphs via Bell Polynomials”, dans JPSS, vol. 4, no 1, 2006, p. 99-105

Articles connexes

modifier