Polynôme de Bell
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 :
Sachant que mi est forcément nul pour i > n − k + 1, on peut expliciter la borne supérieure des indices i :
Si m1 = k − 1 et mn−k+1 = 1, alors ∑ mi = k et ∑ i·mi = n, donc la borne supérieure des sommes doit valoir au minimum n − k + 1 pour couvrir correctement toutes les suites (mi)i≥1 valides.
Supposons qu'il existe un indice j > n − k + 1 tel que mj ≠ 0, soit mj ≥ 1. Les indices commençant à 1, on a également j ≥ 1, soit j − 1 ≥ 0, d'où (j − 1)·mj ≥ j − 1 En outre, les inégalités i ≥ 1 et mi ≥ 0 donnent i·mi ≥ mi. On obtient donc
soit n > n, ce qui est contradictoire. Donc mj = 0 pour tout j > n − k + 1, et la borne supérieure des sommes n'a pas besoin d'aller au-delà de n − k + 1 pour couvrir toutes les suites (mi)i≥1 valides.
Pour j > n − k + 1, on a mj = 0, d'où xjmjmj! (j!)mj = 1, donc la borne supérieure du produit n'a pas besoin d'aller au-delà de n − k + 1.
Interprétation combinatoire
modifierSoit 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
1xm2
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.
Formons une liste avec les n éléments de l'ensemble de départ, et construisons les k sous-ensembles d'arrivée en prenant 1 par 1 les m1 premiers éléments de la liste, puis 2 par 2 les 2·m2 éléments suivants, etc. Il y a donc n! listes différentes à partir desquelles construire ces sous-ensembles d'arrivée.
Pour m1 = 1, m2 = 2 et m3 = 1 (soit k = 4 et n = 8), et pour un ensemble de départ {A, B, …, H}, une des 8! = 40320 listes possibles est CHEFBADG, avec laquelle on construit les sous-ensembles {C}, {H, E}, {F, B} et {A, D, G}.
Cependant, l'ordre des éléments dans chaque sous-ensemble est non pertinent, donc, pour un sous-ensemble donné de cardinalité i, il y a i! listes différentes donnant le même résultat. Il faut donc diviser le nombre de listes m1 fois par 1!, m2 fois par 2!, etc., soit au final par ∏ (i!)mi.
En permutant le 2e et le 3e élément de la liste, soit CEHFBADG, le premier sous-ensemble de cardinalité 2 devient {E, H} lors de sa construction, identique à {H, E}. De même, en permutant le 4e et le 5e élément de la liste, soit CHEBFADG, le second sous-ensemble de cardinalité 2 devient {B, F} lors de sa construction, identique à {F, B}. Enfin, les 3! = 6 permutations des 6e au 8e éléments (par exemple, CHEFBGAD) donnent le même sous-ensemble {A, D, G} de cardinalité 3. Les trois zones de permutation (portant respectivement sur EH, FB et ADG) étant indépendantes, il y a donc 2!·2!·3! = 24 listes différentes permettant de construire les mêmes sous-ensembles {C}, {H, E}, {F, B} et {A, D, G} (attention : en tenant compte de l'ordre des sous-ensembles).
Enfin, les mi sous-ensembles de cardinalité i étant indiscernables, il faut également diviser le résultat précédent par m1!·m2!·… = ∏ mi!.
Les 2 sous-ensembles de cardinalité 2 sont indiscernables, donc la liste CFBHEADG (qui construit les sous-ensembles {C}, {F, B}, {H, E} et {A, D, G}) donne la même partition que CHEFBADG.
Il y a donc au final 24·2! = 48 listes donnant toutes la même partition que CHEFBADG.
De manière générale, chaque partition peut être construite par 48 listes différentes, il y a donc 40320/48 = 840 partitions différentes pour (m1, m2, m3) = (1, 2, 1).
En conclusion, l'un des monômes du polynôme de Bell B8,4(x1, x2, …) est donc 840x3x2
2x1 (note : les variables xi sont notées par ordre décroissant de leur indice, afin de suivre la convention de notation utilisée pour les partitions d'un entier).
Le nombre total de partitions est donc n!∏ (i!)mi · ∏ mi! = n!∏ mi!·(i!)mi, coefficient associé au monôme unitaire xm1
1xm2
2… = ∏ xmi
i.
Exemples
modifierOn a :
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 :
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
modifierLa somme
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 :
avec δk le symbole de Kronecker. La matrice dont Bn est le déterminant est une matrice de Hessenberg.
Table de valeurs
modifierLe tableau suivant regroupe les premières valeurs de :
k n |
0 | 1 | 2 | 3 | 4 | 5 | 6 |
|---|---|---|---|---|---|---|---|
| 0 | |||||||
| 1 | |||||||
| 2 | |||||||
| 3 | |||||||
| 4 | |||||||
| 5 | |||||||
| 6 | |||||||
Propriétés
modifierCas limites
modifier- avec δn le symbole delta de Kronecker
Séries formelles exponentielles
modifier- Puissance
Pour n et k deux entiers naturels, soient j ∈ ℕ* vérifiant j ≥ n − k + 1 et (y1, …, yj) une suite de scalaires. La formule du multinôme de Newton établit que
En posant yi = xi tii!, on obtient
Chaque membre est un polynôme en t, on peut donc en extraire le coefficient cn associé à tn, qui vaut
Comme démontré précédemment (voir en-tête de l'article), les conditions ∑ mi = k et ∑ i·mi = n font que mi = 0 pour tout i > n − k + 1. Comme on a choisi j ≥ n − k + 1, la somme constituant cn couvre bien toutes les suites (mi)i≥1 valides, d'où
On peut donc rendre j arbitrairement grand sans changer cn, qui est donc également le coefficient associé à tn de la série formelle obtenue en faisant tendre j vers l'infini, d'où
Comme Bn,k|n<k = 0, la somme peut commencer à k.
On notera que cette propriété ayant été démontrée en utilisant uniquement la définition des polynômes de Bell (et la formule du multinôme de Newton), elle pourra servir (directement ou indirectement) pour démontrer la plupart des autres propriétés ci-dessous sans causer de raisonnement circulaire.
- Composition
Si
- (avec quelconque)
et
- (avec donc )
alors
Formules de récurrence
modifieravec Bn,0 = δn.
D'après la formule de puissance des séries formelles exponentielles (voir plus haut), on a
d'où par changement de variable k → k + 1
En dérivant chaque membre par rapport à t, on a d'une part
et d'autre part
d'où
Enfin, par unicité de la décomposition des séries formelles exponentielles sur la base , on obtient
La seconde forme s'obtient immédiatement par le changement de variable i → n − i.
avec B0 = B0,0 = 1.
Pour n ≥ 0 on a n + 1 > 0, d'où Bn+1,0 = 0, donc
La seconde forme s'obtient immédiatement par le changement de variable i → n − i.
La matrice étant une matrice de Hessenberg, on peut développer son déterminant selon la dernière colonne, donnant la formule de récurrence.
Valeurs particulières
modifier- (nombre de Stirling de seconde espèce non signé)
- (nombre de Stirling de seconde espèce signé)
- (nombre de Bell)
- (nombre de Stirling de première espèce non signé)
- (nombre de Stirling de première espèce signé)
- (nombre de Lah non signé)
- (nombre de Lah signé)
Type binomial
modifieravec B0 = 1.
Réciproque
modifierSoit f une fonction infiniment dérivable en un point a et de réciproque f -1, alors[1] :
Posons
- (on vérifie que )
ainsi que
d'où
soit
En appliquant la formule de composition des séries formelles exponentielles (voir plus haut) à , on obtient
avec
Calculons maintenant la réciproque de g :
donne
d'où
soit
En appliquant la formule de composition des séries formelles exponentielles à , on obtient
d'où
Cas particuliers
modifierEn prenant f (x) = ex (soit f –1(x) = ln(x)) infiniment dérivable en 0, on a :
d’où :
soit :
En prenant f (x) = xα avec α ≠ 0 (soit f –1(x) = x1/α) infiniment dérivable en 1, on a :
avec .k la factorielle décroissante, d’où :
Composition
modifierSoient :
- (on note que )
Alors :
En posant les matrices (triangulaire supérieure) ainsi que et de manière similaire, on a alors :
Posons
- , avec donc
On a
d'où
- , avec donc
Posons
(Note : la notation pour cette fonction vient de la théorie des espèces combinatoires.) Par associativité de la composition de fonctions, on a , dont on évaluera chaque membre. En appliquant la formule de composition des séries formelles exponentielles (voir plus haut), on obtient
De même, en posant , on a
En appliquant une dernière fois la formule de composition, on obtient
L'égalité donne
Ensuite, par unicité de la décomposition des séries formelles exponentielles sur la base , on obtient
Enfin, comme , on peut changer la borne inférieure de la somme pour obtenir
Cas particuliers
modifierEn prenant et , on obtient :
En prenant et , on obtient :
En prenant et , on obtient :
En prenant et , on obtient :
Factorielle décroissante
modifieravec .k la factorielle décroissante.
Comportement d’échelle
modifierPolynômes de Bell partiels
modifier- Cas général
- Cas particuliers
Polynômes de Bell complets
modifier- Cas général
- Cas particuliers
- Autre expression
avec .k la factorielle décroissante.
Posons
- (on vérifie que )
- (on vérifie que )
Par associativité de la composition de fonctions, on a , dont on évaluera chaque membre.
- Membre de gauche :
Posons , d'où
En appliquant la formule de composition des séries formelles exponentielles (voir plus haut), on obtient
- Membre de droite :
Posons , d'où, par application de la formule de composition
Puis, en appliquant une dernière fois la formule de composition, on obtient
- Résultat
L'égalité donne
Enfin, par unicité de la décomposition des séries formelles exponentielles sur la base , on obtient
Identité de convolution
modifierPour des suites xn, yn, n = 1, 2, …, on peut définir un produit de convolution par :
(les bornes de sommation étant 1 et n − 1, et non 0 et n).
Soit le n-ème terme de la suite
Alors :
Applications
modifierFormule de Faà di Bruno
modifierLa formule de Faà di Bruno peut être énoncée à l'aide des polynômes de Bell de la manière suivante :
Moments et cumulants
modifierPour une variable aléatoire réelle dont le moment ordinaire mr d’ordre r existe, on a :
avec κi les cumulants.
Représentations de suites polynomiales
modifierPour toute suite a1, a2, … de scalaires, soit :
Cette suite de polynômes est de type binomial, c'est-à-dire qu'elle satisfait l'identité binomiale suivante :
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
en considérant cette série comme une série formelle, alors pour tout n :
Notes et références
modifier- ↑ (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
- ↑ (en) Andrzej Korzeniowski, “Binomial Tails Domination for Random Graphs via Bell Polynomials”, dans JPSS, vol. 4, no 1, 2006, p. 99-105
- (en) Eric Temple Bell, « Partition Polynomials », Ann. Math., vol. 29, nos 1/4, 1927-1928, p. 38-46 (DOI 10.2307/1967979)
- (en) Eric Temple Bell, « Exponential Polynomials », Ann. Math., vol. 35, no 2, , p. 258-277 (DOI 10.2307/1968431)
- (en) Louis Comtet, Advanced Combinatorics: The Art of Finite and Infinite Expansions, Reidel Publishing Company, Dordrecht-Holland/Boston-U.S., 1974
- (en) Steven Roman (en), The Umbral Calculus, Dover Publications
- (en) Cet article est partiellement ou en totalité issu de l’article de Wikipédia en anglais intitulé « Bell polynomials » (voir la liste des auteurs).