Nombre de Stirling

En mathématiques, les nombres de Stirling apparaissent dans plusieurs problèmes combinatoires. Ils tirent leur nom de James Stirling, qui les a introduits au XVIIIe siècle. Il en existe trois sortes, nommés les nombres de Stirling de première espèce signés et non signés, et les nombres de Stirling de seconde espèce.

Notations

Diverses notations sont utilisées pour les nombres de Stirling, parmi lesquelles[1] :

  • nombres de Stirling de première espèce « signés » : s ( n , k ) {\displaystyle s(n,k)}  ;
  • nombres de Stirling de première espèce « non signés » : | s ( n , k ) | = c ( n , k ) = [ n k ] {\displaystyle |s(n,k)|=c(n,k)=\left[{n \atop k}\right]}  ;
  • nombres de Stirling de seconde espèce : S ( n , k ) = S n ( k ) = { n k } {\displaystyle S(n,k)=S_{n}^{(k)}=\left\{{n \atop k}\right\}} .

La notation avec crochets, analogue à celle utilisée pour les coefficients binomiaux, est due à Jovan Karamata, qui l'a proposée en 1935. Son usage a été encouragé par Donald Knuth[2] mais, outre son ergonomie discutable[3], elle comporte un risque de confusion avec les coefficients binomiaux de Gauss (présentés dans l'article « q-analogue »). Nous nous limiterons donc, pour chacun des trois types de nombres, à la première notation correspondante ci-dessus.

Nombre de Stirling de première espèce

Les nombres de Stirling de première espèce signés s(n, k) ,pour n, k entiers naturels, sont les coefficients du développement de la factorielle décroissante (x)n, c'est-à-dire que

( x ) n = x ( x 1 ) ( x 2 ) ( x n + 1 ) = k = 0 n s ( n , k ) x k {\displaystyle (x)_{n}=x(x-1)(x-2)\ldots (x-n+1)=\sum _{k=0}^{n}s(n,k)x^{k}}

((x)0 = 1 car il s'agit d'un produit vide).

Le nombre s(n, k) a même signe que (–1)n – k.

Les nombres de Stirling de première espèce non signés |s(n, k)| (valeurs absolues des précédents) sont les coefficients du développement de la factorielle croissante (x)n, c'est-à-dire que

( x ) n = x ( x + 1 ) ( x + 2 ) ( x + n 1 ) = k = 0 n | s ( n , k ) | x k {\displaystyle (x)^{n}=x(x+1)(x+2)\ldots (x+n-1)=\sum _{k=0}^{n}|s(n,k)|x^{k}} .

Ils ont aussi une définition combinatoire : cf. § #Interprétation combinatoire ci-dessous.

Table de valeurs

Voici une table donnant quelques valeurs des s(n, k) (suites OEIS A008275 et OEIS A008276 de l'OEIS), que l'on peut calculer ligne par ligne grâce à la relation de récurrence du § suivant, de même que le triangle de Pascal :

n \ k 0 1 2 3 4 5 6 7 8 9
0 1
1 0 1
2 0 –1 1
3 0 2 –3 1
4 0 –6 11 –6 1
5 0 24 –50 35 –10 1
6 0 –120 274 –225 85 –15 1
7 0 720 –1764 1624 –735 175 –21 1
8 0 –5040 13068 –13132 6769 –1960 322 –28 1
9 0 40320 –109584 118124 –67284 22449 –4536 546 –36 1

Expression sommatoire

En partant de la relation ( X 1 ) ( X 2 ) . . ( X n ) = k = 0 n s ( n + 1 , k + 1 ) X k {\displaystyle (X-1)(X-2)..(X-n)=\sum _{k=0}^{n}s(n+1,k+1)X^{k}} , et en utilisant les relations entre coefficients et racines, on obtient [4] :

σ k ( 1 , 2 , , n ) = 1 i 1 < i 2 < . . < i k n i 1 . . i k = ( 1 ) k s ( n + 1 , n + 1 k ) {\displaystyle \sigma _{k}(1,2,\cdots ,n)=\sum _{1\leqslant i_{1}<i_{2}<..<i_{k}\leqslant n}i_{1}..i_{k}=(-1)^{k}s(n+1,n+1-k)} , pour n k 1 {\displaystyle n\geqslant k\geqslant 1} .
Exemple

Pour tout entier n supérieur ou égal à 1 :

s ( n + 1 , n 1 ) = 1 i < j n i j = ( n 1 ) n ( n + 1 ) ( 3 n + 2 ) 24 {\displaystyle s(n+1,n-1)=\sum _{1\leqslant i<j\leqslant n}ij={\frac {(n-1)n(n+1)(3n+2)}{24}}} , voir la suite A000914 de l'OEIS.

Formule de récurrence

Les nombres de Stirling de première espèce signés vérifient la relation de récurrence

k 1 s ( n + 1 , k ) = s ( n , k 1 ) n s ( n , k ) {\displaystyle \forall k\geqslant 1\quad s(n+1,k)=s(n,k-1)-ns(n,k)}

avec les conditions initiales

s ( 0 , 0 ) = 1 e t n 1 s ( n , 0 ) = s ( 0 , n ) = 0 {\displaystyle s(0,0)=1\quad {\rm {et}}\quad \forall n\geqslant 1\quad s(n,0)=s(0,n)=0} .

Leurs valeurs absolues vérifient (avec mêmes conditions initiales) la relation de récurrence

k 1 | s ( n + 1 , k ) | = | s ( n , k 1 ) | + n | s ( n , k ) | {\displaystyle \forall k\geqslant 1\quad |s(n+1,k)|=|s(n,k-1)|+n|s(n,k)|} .

Chacune des deux relations de récurrence peut se déduire de l'autre. De plus, la première découle de la relation de récurrence des factorielles décroissantes :

( x ) n + 1 = x ( x ) n n ( x ) n {\displaystyle (x)_{n+1}=x(x)_{n}-n(x)_{n}}

et la seconde, d'un raisonnement combinatoire ou de la relation de récurrence des factorielles croissantes :

( x ) n + 1 = x ( x ) n + n ( x ) n {\displaystyle (x)^{n+1}=x(x)^{n}+n(x)^{n}} .

Identités simples

Remarquons que

k > n s ( n , k ) = 0 {\displaystyle \forall k>n\quad s(n,k)=0} ,
s ( n , n ) = 1 , s ( n , 1 ) = ( 1 ) n 1 ( n 1 ) ! , s ( n , n 1 ) = ( n 2 ) {\displaystyle s(n,n)=1,\quad s(n,1)=(-1)^{n-1}(n-1)!,\quad s(n,n-1)=-{n \choose 2}} ,
s ( n , n 2 ) = 3 n 1 4 ( n 3 ) , s ( n , n 3 ) = ( n 2 ) ( n 4 ) {\displaystyle s(n,n-2)={\frac {3n-1}{4}}{n \choose 3},\quad s(n,n-3)=-{n \choose 2}{n \choose 4}} .

Il existe d'autres identités, comme

s ( n , 2 ) = ( 1 ) n ( n 1 ) ! H n 1 {\displaystyle s(n,2)=(-1)^{n}(n-1)!\;H_{n-1}}

Hn est un nombre harmonique et

s ( n , 3 ) = ( 1 ) n 1 ( n 1 ) ! 2 [ ( H n 1 ) 2 H n 1 ( 2 ) ] {\displaystyle s(n,3)=(-1)^{n-1}{\frac {(n-1)!}{2}}\left[(H_{n-1})^{2}-H_{n-1}^{(2)}\right]}

Hn(m) est un nombre harmonique généralisé.

Des relations similaires lient les nombres de Stirling de première espèce aux polynômes de Bernoulli. Un grand nombre de relations liées aux nombres de Stirling cachent des relations similaires liées aux coefficients binomiaux. L'étude des relations entre ces deux nombres est le calcul ombral et est un domaine important de la théorie des suites de Sheffer.

Formules explicites

On peut montrer[5] la relation suivante entre nombres de Stirling de première et seconde espèce :

s ( n , m ) = k = 0 n m ( 1 ) k ( n 1 + k n m + k ) ( 2 n m n m k ) S ( n m + k , k ) {\displaystyle s(n,m)=\sum _{k=0}^{n-m}(-1)^{k}{\binom {n-1+k}{n-m+k}}{\binom {2n-m}{n-m-k}}S(n-m+k,k)}

d'où, utilisant la formule pour ces derniers qui sera donnée plus bas :

s ( n , m ) = k = 0 n m j = 0 k ( 1 ) 2 k j k ! ( n 1 + k n m + k ) ( 2 n m n m k ) ( k j ) j n m + k {\displaystyle s(n,m)=\sum _{k=0}^{n-m}\sum _{j=0}^{k}{\frac {(-1)^{2k-j}}{k!}}{\binom {n-1+k}{n-m+k}}{\binom {2n-m}{n-m-k}}{k \choose j}j^{n-m+k}}

ou encore, après simplifications :

s ( n , m ) = ( 2 n m ) ! ( m 1 ) ! k = 0 n m 1 ( n + k ) ( n m k ) ! ( n m + k ) ! j = 0 k ( 1 ) j j n m + k j ! ( k j ) ! {\displaystyle s(n,m)={\frac {(2n-m)!}{(m-1)!}}\sum _{k=0}^{n-m}{\frac {1}{(n+k)(n-m-k)!(n-m+k)!}}\sum _{j=0}^{k}{\frac {(-1)^{j}j^{n-m+k}}{j!(k-j)!}}} .

Fonction génératrice

On peut démontrer plusieurs identités en manipulant la fonction génératrice :

( 1 + t ) x = n = 0 ( x n ) t n = n = 0 t n n ! k = 0 n s ( n , k ) x k = k = 0 x k n = k t n n ! s ( n , k ) = e x ln ( 1 + t ) {\displaystyle (1+t)^{x}=\sum _{n=0}^{\infty }{x \choose n}t^{n}=\sum _{n=0}^{\infty }{\frac {t^{n}}{n!}}\sum _{k=0}^{n}s(n,k)x^{k}=\sum _{k=0}^{\infty }x^{k}\sum _{n=k}^{\infty }{\frac {t^{n}}{n!}}s(n,k)=\operatorname {e} ^{x\ln(1+t)}} .

En particulier, on peut inverser l'ordre de la sommation et prendre des dérivées, puis fixer t ou x.

Sommes finies

Si on prend x = -1, qu'on développe ( 1 + t ) 1 {\displaystyle (1+t)^{-1}} et qu'on identifie les puissances de t dans les deux développements, on obtient :

k = 0 n ( 1 ) k s ( n , k ) = ( 1 ) n n ! {\displaystyle \sum _{k=0}^{n}(-1)^{k}s(n,k)=(-1)^{n}n!}

Sommes infinies

Si on développe l'exponentielle e x ln ( 1 + t ) {\displaystyle e^{x\ln(1+t)}} et qu'on identifie les puissances de x dans les deux développements, on obtient :

n = k s ( n , k ) t n n ! = ( ln ( 1 + t ) ) k k ! {\displaystyle \sum _{n=k}^{\infty }s(n,k){\frac {t^{n}}{n!}}={\frac {\left(\ln(1+t)\right)^{k}}{k!}}} ,

valide pour | t | < 1 {\displaystyle |t|<1} .

Interprétation combinatoire

La valeur absolue du nombre de Stirling de première espèce compte le nombre de permutations de n objets composés d'exactement k cycles disjoints. Par exemple, s ( 4 , 2 ) = 11 {\displaystyle s(4,2)=11} correspond au fait que le groupe symétrique S 4 {\displaystyle S_{4}} possède trois permutations de la forme

( ) ( ) {\displaystyle (**)(**)} — 2 cycles de longueur 2

et huit permutations de la forme

( ) {\displaystyle (***)} — 1 cycle de longueur 3 et 1 cycle de longueur 1.

La valeur absolue du nombre de Stirling de première espèce compte aussi le nombre de permutations de n objets ayant exactement k records. Cette identité entre records et cycles résulte de la correspondance fondamentale de Foata. La forme produit de la série génératrice des nombres de Stirling de première espèce résulte de l'indépendance des termes du code de Lehmer d'une permutation, code très lié aux records d'une permutation. L'interprétation des nombres de Stirling en fonction du nombre de records explique l'apparition des nombres de Stirling dans l'analyse de l'algorithme de recherche du maximum, qui est la première analyse d'algorithme traitée dans le livre fondateur de Knuth, The Art of Computer Programming.

Nombre de Stirling de seconde espèce

Définition combinatoire

Les nombres de Stirling de seconde espèce S(n, k) sont définis combinatoirement de trois façons équivalentes :

  • S(n, k) est le nombre de relations d'équivalence ayant k classes d'équivalence définies sur un ensemble de n éléments ;
  • S(n, k) est le nombre de partitions d'un ensemble à n éléments en k sous-ensembles ;
  • k! × S(n, k) est le nombre de surjections d'un ensemble à n éléments sur un ensemble à k éléments. Attention, ce dernier nombre est lui aussi souvent noté S(n, k) ou s(n, k).

Formule explicite

Les nombres de Stirling de seconde espèce sont donnés par la formule explicite

S ( n , k ) = 1 k ! j = 0 k ( 1 ) k j ( k j ) j n {\displaystyle S(n,k)={\frac {1}{k!}}\sum _{j=0}^{k}(-1)^{k-j}{\binom {k}{j}}j^{n}} ,

laquelle s'obtient par exemple[6],[7] en remarquant que le nombre de surjections (d'un ensemble à n éléments vers un ensemble à k éléments) peut se compter par la formule d'inclusion-exclusion : on compte toutes les applications moins celles n'atteignant pas un certain élément, plus celles n'atteignant pas deux éléments, moins...

Relation de récurrence

À partir de la définition combinatoire, on peut également démontrer[8],[7] que ces nombres vérifient la relation de récurrence

n , k 1 S ( n , k ) = S ( n 1 , k 1 ) + k S ( n 1 , k ) {\displaystyle \forall n,k\geqslant 1\quad S(n,k)=S(n-1,k-1)+kS(n-1,k)}

avec les conditions initiales

S ( 0 , 0 ) = 1 e t n 1 S ( n , 0 ) = S ( 0 , n ) = 0 {\displaystyle S(0,0)=1\quad {\rm {et}}\quad \forall n\geqslant 1\quad S(n,0)=S(0,n)=0} .

Caractérisation algébrique

On déduit de la relation de récurrence ci-dessus[9],[7] que

k = 0 n S ( n , k ) ( X ) k = X n {\displaystyle \sum _{k=0}^{n}S(n,k)(X)_{k}=X^{n}} ,

( X ) k = X ( X 1 ) . . . ( X k + 1 ) {\displaystyle (X)_{k}=X(X-1)...(X-k+1)} (symbole de Pochhammer),

ce qui fournit une définition algébrique des nombres de Stirling de deuxième espèce, équivalente à la définition combinatoire initiale.

Cette formule est utilisée par exemple dans l'expression des sommes k = 1 n k p {\displaystyle \sum _{k=1}^{n}k^{p}} .

Table de valeurs

Voici quelques valeurs des nombres de Stirling de seconde espèce (suites OEIS A008277 et OEIS A008278 de l'OEIS) :

n \ k 1 2 3 4 5 6 7 8 9
0 1
1 0 1
2 0 1 1
3 0 1 3 1
4 0 1 7 6 1
5 0 1 15 25 10 1
6 0 1 31 90 65 15 1
7 0 1 63 301 350 140 21 1
8 0 1 127 966 1701 1050 266 28 1
9 0 1 255 3025 7770 6951 2646 462 36 1

Identités simples

On a par exemple S(n, n) = 1,

S ( n , n 1 ) = ( n 2 ) {\displaystyle S(n,n-1)={\binom {n}{2}}}

et

S ( n , 2 ) = 2 n 1 1 {\displaystyle S(n,2)=2^{n-1}-1} .

Le nombre total de partitions d'un ensemble à n éléments,

B n = k = 1 n S ( n , k ) {\displaystyle B_{n}=\sum _{k=1}^{n}S(n,k)} ,

est le n-ième nombre de Bell.

Rapport avec la distribution de Poisson

Si X est une variable aléatoire suivant une distribution de Poisson avec une moyenne λ, alors son n-ième moment est

E ( X n ) = k = 1 n S ( n , k ) λ k {\displaystyle E(X^{n})=\sum _{k=1}^{n}S(n,k)\lambda ^{k}} .

En particulier, le n-ième moment d'une distribution de Poisson de moyenne 1 est précisément le nombre de partitions d'un ensemble de taille n, qui est le n-ième nombre de Bell (formule de Dobinski).

Relation de réciprocité

D'après leur caractérisation algébrique, les nombres de Stirling de première et seconde espèce, disposés, comme dans les tables de valeurs ci-dessus, en deux matrices triangulaires infinies, constituent les deux matrices de passage (dans un sens et dans l'autre) entre deux bases de l'espace des polynômes[10] : la base canonique des monômes Xk et la base des symboles de Pochhammer X(X – 1)(X – 2)…(Xk + 1). Le fait que ces deux matrices sont inverses l'une de l'autre se traduit par[11] :

m k n s ( k , m ) S ( n , k ) = m k n S ( k , m ) s ( n , k ) = δ m n {\displaystyle \sum _{m\leq k\leq n}s(k,m)S(n,k)=\sum _{m\leq k\leq n}S(k,m)s(n,k)=\delta _{mn}}

δ m n {\displaystyle \delta _{mn}} est le symbole de Kronecker.

Notes et références

(en) Cet article est partiellement ou en totalité issu de l’article de Wikipédia en anglais intitulé « Stirling number » (voir la liste des auteurs).
  1. D'autres encore sont utilisées par Abramowitz et Stegun.
  2. (en) D. E. Knuth, Two notes on notation (source TeX).
  3. Louis Comtet, Analyse combinatoire élémentaire, Techniques Ingénieur, coll. « Mathématiques concrètes », , 687 p. (ISBN 978-2-84180-981-3, lire en ligne), p. 24.
  4. Comtet 1970, p. 48
  5. Pour une démonstration, voir par exemple Louis Comtet, Analyse combinatoire, t. 2, PUF, , p. 51.
  6. Comtet 1970, p. 38.
  7. a b et c Voir aussi cet exercice corrigé sur Wikiversité.
  8. Voir par exemple Louis Comtet, Analyse combinatoire approfondie, Techniques Ingénieur, coll. « Mathématiques concrètes », (lire en ligne), p. 3-4.
  9. Comtet 2003, p. 5.
  10. Espace vectoriel K[X] des polynômes en une indéterminée sur un corps commutatif K, ou même module libre A[X] des polynômes à coefficients dans un anneau commutatif A.
  11. Abramowitz et Stegun, p. 825.

Voir aussi

Bibliographie

  • (en) Milton Abramowitz et Irene Stegun, Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables [détail de l’édition] (lire en ligne), « Stirling Numbers », § 24.1.3-4, p. 824-825
  • (en) Louis Comtet, Advanced Combinatorics : The Art of Finite and Infinite Expansions, D. Reidel, , 343 p. (ISBN 978-90-277-0441-2, lire en ligne)
  • (en) Victor Adamchik, « On Stirling Numbers and Euler Sums », Journal of Computational and Applied Mathematics, vol. 79,‎ , p. 119-130 (lire en ligne)
  • (en) Arthur T. Benjamin, Gregory O. Preston et Jennifer J. Quinn, « A Stirling Encounter with Harmonic Numbers », Mathematics Magazine, vol. 75, no 2,‎ , p. 95-103 (lire en ligne)
  • (en) J. M. Sixdeniers, K. A. Penson et A. I. Solomon, « Extended Bell and Stirling Numbers From Hypergeometric Exponentiation », Journal of Integer Sequences, vol. 4, no 01.1.4,‎ (lire en ligne)
  • (en) Hsien-Kuei Hwang, « Asymptotic Expansions for the Stirling numbers of the first kind », Journal of Combinatorial Theory, a, vol. 71, no 2,‎ , p. 343-351 (DOI 10.1016/0097-3165(95)90010-1, lire en ligne)

Articles connexes

Liens externes

  • (en) Eric W. Weisstein, « Stirling Number of the First Kind », sur MathWorld
  • (en) E. Weisstein, « Stirling Number of the Second Kind », sur MathWorld
  • (en) rmilson (146), « Stirling numbers of the first kind », sur PlanetMath
  • (en) rmilson (146), « Stirling numbers of the second kind », sur PlanetMath
  • icône décorative Portail des mathématiques