Coefficient Binomial Central : un aperçu

Ce qui suit suppose connues les propriétés de base des coefficients binomiaux. On pourra au besoin se reporter à cet article.

Par définition, le n-ème coefficient binomial central est l’entier :

    \[C_{n}=\binom{2n}{n}=\frac{\left(2n\right)!}{\left(n!\right)^{2}}\]

Il est évidemment dénommé ainsi en raison de sa position dans le triangle de Pascal (les premières valeurs de C_{n} apparaissent en rouge dans le tableau ci-dessous) :

Nous allons examiner quelques unes des nombreuses propriétés de cette suite d’entiers, en commençant par ce qu’il y a de plus simple : la parité.

1 – Parité des coefficients binomiaux centraux

Rappelons la formule de Pascal. Pour tout couple d’entiers \left(n,k\right) tel que 1\leqslant k<n :

    \[\binom{n}{k}=\binom{n-1}{k}+\binom{n-1}{k-1}\]

Il en résulte, en appliquant ceci au couple \left(2n,n\right) que :

    \[\binom{2n}{n}=\binom{2n-1}{n}+\binom{2n-1}{n-1}\]

Mais d’après la formule de symétrie :

    \[\binom{2n-1}{n-1}=\binom{2n-1}{\left(2n-1\right)-\left(n-1\right)}=\binom{2n-1}{n}\]

et donc :

    \[\fcolorbox{black}{myBlue}{$\displaystyle{\forall n\in\mathbb{N}^{\star},\quad\binom{2n}{n}=2\thinspace\binom{2n-1}{n}}$}\]

ce qui prouve la parité de {\displaystyle \binom{2n}{n}.} Noter que pour n=0, cet entier est impair (il vaut 1).

Remarque

Il y a plus simple … D’après la formule du pion, pour tout entier n\geqslant1 :

    \[n\binom{2n}n=2n\binom{2n-1}{n-1}\]

donc, après simplification :

    \[\binom{2n}n=2\binom{2n-1}{n-1}\]

2 – Un résultat plus précis

Intéressons-nous, pour tout n\geqslant1, à l’exposant de 2 dans la décomposition en facteurs premiers de C_{n}. Notons V_{n} cet entier. Nous venons de voir que V_{n}\geqslant1 mais il est naturel de se demander s’il est possible de calculer explicitement V_{n}.

Pour cela, nous allons utiliser la formule de Legendre (dont une preuve détaillée est disponible dans l’annexe de cet article). Rappelons de quoi il s’agit :

Formule de Legendre / Version 1

Pour tout entier n\geqslant1 et tout nombre premier p, l’exposant de p dans la décomposition en facteurs premiers de n! est donné par :

    \[v_{p}\left(n!\right)=\sum_{i=1}^{n}\left\lfloor \frac{n}{p^{i}}\right\rfloor\]

Dans cette formule, \left\lfloor X\right\rfloor désigne la partie entière du réel X. C’est ainsi que, par exemple, l’exposant de 2 dans la décomposition en facteurs premiers de 100! est :

    \[v_{2}\left(100!\right)=\left\lfloor 50\right\rfloor +\left\lfloor 25\right\rfloor +\left\lfloor 12,5\right\rfloor +\left\lfloor 6,25\right\rfloor +\left\lfloor 3,125\right\rfloor +\left\lfloor 1,5625\right\rfloor\]

(les termes suivants sont nuls), c’est-à-dire :

    \[v_{2}\left(100!\right)=50+25+12+6+3+1=97\]

Autrement dit, l’entier 100! est de la forme 2^{97}q avec q impair. Mais revenons à notre sujet …

Comme la valuation p-adique est additive (ce qui veut simplement dire que v_{p}\left(ab\right)=v_{p}\left(a\right)+v_{p}\left(b\right) pour tout couple \left(a,b\right) d’entiers naturels non nuls), on voit que :

    \[V_{n}=v_{2}\left(\left(2n\right)!\right)-2\thinspace v_{2}\left(n!\right)\]

Afin de poursuivre notre calcul, nous avons besoin d’une autre version de la formule de Legendre, à savoir :

Formule de Legendre / Version 2

Pour tout entier n\geqslant1 et tout nombre premier p :

    \[v_{p}\left(n!\right)=\frac{n-S_{p}\left(n\right)}{p-1}\]

S_{p}\left(n\right) désigne la somme des chiffres de l’écriture de n en base p.

Si la notion d’écriture en base p n’est pas claire pour vous, surtout pas de panique !… Tout est expliqué dans ces deux vidéos :

Quant à la preuve de cette seconde version de la formule de Legendre, les curieux la trouveront à l’annexe 1.

Avec ce nouvel outil, on constate que :

    \begin{eqnarray*}V_{n} & = & \left(2n-S_{2}\left(2n\right)\right)-2\left(n-S_{2}\left(n\right)\right)\\& = & 2\thinspace S_{2}\left(n\right)-S_{2}\left(2n\right)\end{eqnarray*}

Enfin, lorsqu’on écrit un entier en base 2, il est très simple de le multiplier par 2 : il suffit d’ajouter un zéro à la fin, ce qui ne modifie évidemment pas la somme des chiffres binaires ! Par conséquent :

    \[S_{2}\left(2n\right)=S_{2}\left(n\right)\]

et finalement :

    \[\boxed{V_{n}=S_{2}\left(n\right)}\]

En résumé, nous avons établi le :

Théorème

Pour tout entier n\geqslant1, l’exposant de 2 dans la décomposition en facteurs premiers de {\displaystyle \binom{2n}{n}} est donné par le nombre de 1 dans l’écriture binaire de n.

On notera que le nombre de 1 dans l’écriture binaire de n est évidemment égal à la somme des chiffres binaires de n.

Un exemple, juste pour voir … Choisissons n=15. On calcule :

    \[C_{15}=\binom{30}{15}=155\thinspace117\thinspace520\]

puis on décompose 15 en binaire :

    \[15=1111_{2}\]

et donc, l’exposant de 2 dans la décomposition en facteurs premiers
de C_{15} est :

    \[V_{15}=4\]

On peut confirmer cela en observant que :

    \[C_{15}=2^{4}\times3^{2}\times5\times17\times19\times23\times29\]

Dans cette décomposition en facteurs premiers, on voit bien que l’exposant de 2 vaut 4.

Corollaire

Une condition nécessaire et suffisante pour que {\displaystyle \binom{2n}{n}} soit de la forme 2q, avec q impair, est que n soit une puissance de 2.

En effet, cette condition signifie que l’exposant de 2 dans la décomposition en facteurs premiers de {\displaystyle \binom{2n}{n}} est égal à 1, ou encore que l’écriture binaire de n comporte exactement une fois le chiffre 1. Les entiers n vérfifiant cette dernière condition sont évidemment les puissances de 2.

3 – Quatre formules pour le n-ème coefficient binomial central

Nous allons établir les formules suivantes :

(F_1)   \[\binom{2n}{n}=\sum_{k=0}^{n}\binom{n}{k}^{2}\]

(F_2)   \[\binom{2n}{n}=\sum_{k=0}^{n}\left(-1\right)^{n-k}\binom{2n+1}{k}\]

(F_3)   \[\binom{2n}{n}=\sum_{k=0}^{2n}\left(-1\right)^{n-k}\binom{2n}{k}^{2}\]

(F_4)   \[\binom{2n}{n}=\sum_{j=0}^{\left\lfloor n/2\right\rfloor }\binom{n}{j}\binom{n-j}{n-2j}2^{n-2j}\]

Cette courte liste n’a évidemment rien d’exhaustif ! Il existe une foule de formules donnant une expression plus ou moins exotique du n-ème coefficient binomial central.

Preuve de la formule n° 1

D’après la formule du binôme, l’entier {\displaystyle \binom{2n}{n}} s’interprète naturellement comme le coefficient de X^{n} dans l’expression développée du polynôme \left(1+X\right)^{2n}. Par ailleurs :

    \begin{eqnarray*}\left(1+X\right)^{2n} & = & \left(1+X\right)^{n}\left(1+X\right)^{n}\\& = & \left(\sum_{k=0}^{n}\binom{n}{k}X^{k}\right)\left(\sum_{\ell=0}^{n}\binom{n}{\ell}X^{\ell}\right)\end{eqnarray*}

et donc, ce coefficient peut aussi s’écrire :

    \[\binom{2n}{n}=\sum_{k=0}^{n}\binom{n}{k}\binom{n}{n-k}\]

On invoque alors la symétrie des coefficients binomiaux, pour conclure que :

    \[\fcolorbox{black}{myBlue}{$\displaystyle{\binom{2n}{n}=\sum_{k=0}^{n}\binom{n}{k}^{2}}$}\]

Preuve de la formule n° 2

Cette question a fait l’objet du challenge 75. Le lecteur est prié de s’y reporter.

Preuve de la formule n° 3

Cette question est traitée dans cette vidéo (en anglais 🇺🇸🇬🇧) :

Mais donnons tout de même une preuve détaillée (en français 🙂)…

Pour tout n\in\mathbb{N}, notons \lambda_{n} le coefficient de X^{n} dans l’expression développée du polynôme :

    \[P_{n}=\left(1-X^{2}\right)^{n}\]

D’une part (formule du binôme) :

    \[P_{n}=\sum_{k=0}^{2n}\left(-1\right)^{k}\binom{2n}{k}X^{2k}\]

et donc :

    \[\lambda_{n}=\left\{\begin{array}{cc}0 & \text{si }n\text{ est impair}\\\left(-1\right)^{n/2}{\displaystyle \binom{n}{n/2}} & \text{sinon}\end{array}\right.\]

Et d’autre part :

    \[P_{n}=\left(1-X\right)^{n}\left(1+X\right)^{n}=\left(\sum_{k=0}^{n}\left(-1\right)^{k}\binom{n}{k}X^{k}\right)\left(\sum_{\ell=0}^{n}\binom{n}{\ell}X^{\ell}\right)\]

d’où (d’après la symétrie des coefficients binomiaux) :

    \[\lambda_{n}=\sum_{k=0}^{n}\left(-1\right)^{k}\binom{n}{k}\binom{n}{n-k}=\sum_{k=0}^{n}\left(-1\right)^{k}\binom{n}{k}^{2}\]

En remplaçant n par 2n, on conclut que :

    \[\left(-1\right)^{n}\binom{2n}{n}=\sum_{k=0}^{2n}\left(-1\right)^{k}\binom{2n}{k}^{2}\]

Il reste à multiplier chaque membre par \left(-1\right)^{n} pour obtenir :

    \[\fcolorbox{black}{myBlue}{$\displaystyle{\binom{2n}{n}=\sum_{k=0}^{2n}\left(-1\right)^{n-k}\binom{2n}{k}^{2}}$}\]

Preuve de la formule n° 4

Cette fois, considérons plutôt le polynôme :

    \[Q_{n}=\left(1+X\right)^{2n}\]

D’une part :

    \[Q_{n}=\sum_{k=0}^{2n}\binom{2n}{k}X^{k}\]

et d’autre part :

    \begin{eqnarray*}Q_{n} & = & \left(1+2X+X^{2}\right)^{n}\\& = & \sum_{k=0}^{n}\binom{n}{k}X^{k}\left(X+2\right)^{k}\\& = & \sum_{k=0}^{n}\binom{n}{k}X^{k}\left(\sum_{j=0}^{k}\binom{k}{j}2^{k-j}X^{j}\right)\\& = & \sum_{k=0}^{n}\sum_{j=0}^{k}\binom{n}{k}\binom{k}{j}2^{k-j}X^{j+k}\end{eqnarray*}

Or, on sait que :

    \[\binom{n}{k}\binom{k}{j}=\binom{n}{j}\binom{n-j}{k-j}\]

et donc, après interversion des sommes :

    \[Q_{n}=\sum_{j=0}^{n}\sum_{k=j}^{n}\binom{n}{j}\binom{n-j}{k-j}2^{k-j}X^{j+k}\]

Pour finir, en effectue un petit changement d’indice en posant \ell=j+k dans la somme interne, ce qui donne :

    \[Q_{n}=\sum_{j=0}^{n}\sum_{\ell=2j}^{n+j}\binom{n}{j}\binom{n-j}{\ell-2j}2^{\ell-2j}X^{\ell}\]

Il ne reste plus qu’à extraire le coefficient de X^{n} pour obtenir :

    \[\fcolorbox{black}{myBlue}{$\displaystyle{\binom{2n}{n}=\sum_{j=0}^{\left\lfloor n/2\right\rfloor }\binom{n}{j}\binom{n-j}{n-2j}2^{n-2j}}$}\]

Remarque

Cette dernière formule redonne (mais d’une manière bien compliquée !) la parité de {\displaystyle \binom{2n}{n}} qui a été établie au tout début du présent article.

En effet, si n est impair, alors tous les termes de la somme sont pairs et donc {\displaystyle \binom{2n}{n}} aussi. Et si n est pair, alors tous les termes sont pairs sauf peut-être celui d’indice j=n/2; par conséquent, {\displaystyle \binom{2n}{n}} est de la même parité que {\displaystyle \binom{n}{n/2}.} En posant n=2^{r}q avec q impair, une récurrence immédiate montre que
{\displaystyle \binom{2n}{n}} est de la même parité que {\displaystyle \binom{2q}{q},} donc est pair.

4 – Fonction génératrice de la suite des coefficients binomiaux centraux

Il est connu que pour tout \alpha\in\mathbb{R} et pour tout x\in\left]-1,1\right[ :

    \[\left(1+x\right)^{\alpha}=\sum_{n=0}^{\infty}\binom{\alpha}{n}x^{n}\]

{\displaystyle \binom{\alpha}{n}} désigne la coefficient binomial généralisé, défini par :

    \[\binom{\alpha}{n}=\frac{\alpha\left(\alpha-1\right)\cdots\left(\alpha-n+1\right)}{n!}\]

En choisissant \alpha=-\frac{1}{2} et en remplaçant x par -4x, on voit que pour tout x\in\left]-\frac{1}{4},\frac{1}{4}\right[ :

    \[\frac{1}{\sqrt{1-4x}}=\sum_{n=0}^{\infty}\binom{-1/2}{n}\left(-4x\right)^{n}\]

Or :

    \[\binom{-1/2}{n}=\frac{1}{n!}\prod_{k=0}^{n-1}\left(-\frac{1}{2}-k\right)=\frac{\left(-1\right)^{n}}{2^{n}\thinspace n!}\prod_{k=0}^{n-1}\left(2k+1\right)\]

et ce dernier produit peut s’écrire (en bouchant les trous, c’est-à-dire en multipliant et en redivisant aussitôt par le produit des entiers pairs de 2 à 2n) :

    \[\prod_{k=0}^{n-1}\left(2k+1\right)=\frac{\left(2n\right)!}{2\times4\times\cdots\times\left(2n\right)}=\frac{\left(2n\right)!}{2^{n}n!}\]

Par conséquent :

    \[\binom{-1/2}{n}=\frac{\left(-1\right)^{n}\left(2n\right)!}{4^{n}\left(n!\right)^{2}}\]

et finalement :

    \[\fcolorbox{black}{myBlue}{$\displaystyle{\frac{1}{\sqrt{1-4x}}=\sum_{n=0}^{\infty}\binom{2n}{n}x^{n}}$}\]

On dit que {\displaystyle x\mapsto\frac{1}{\sqrt{1-4x}}} est la fonction génératrice de la suite des coefficients binomiaux centraux.

A partir de cette formule, on peut trouver plusieurs sommes de séries numériques, en donnant une valeur particulière à x. Par exemple, pour x=\frac{1}{8} :

    \[\boxed{\sum_{n=0}^{\infty}\frac{1}{8^{n}}\binom{2n}{n}=\sqrt{2}}\]

5 – Une formule remarquable pour \zeta\left(2\right)

Dans sa tentative de résolution du problème de Bâle, le mathématicien écossais James Stirling découvrit une formule inédite pour le nombre :

    \[\sum_{n=1}^{\infty}\frac{1}{n^{2}}\]

que l’on note aujourd’hui \zeta\left(2\right). Il démontra que :

    \[\zeta\left(2\right)=\sum_{n=1}^{\infty}\frac{3}{n^{2}\binom{2n}{n}}\]

ce qui lui permit d’obtenir une bonne valeur approchée de ce nombre (voir plus bas).

Cependant, Stirling ne parvint pas à reconnaître qu’il s’agissait de {\displaystyle \frac{\pi^{2}}{6}} … ce que fit par la suite Euler. A ce sujet, on pourra consulter cette vidéo :

Donnons une preuve détaillée de la formule :

    \[\fcolorbox{black}{myBlue}{$\displaystyle{\frac{\pi^2}{6}=\sum_{n=1}^{\infty}\frac{3}{n^{2}\binom{2n}{n}}}$}\]

Admettons temporairement que pour tout x\in\left]-1,1\right[ :

    \[\fcolorbox{black}{myBlue}{$\displaystyle{\frac{\arcsin\left(x\right)}{\sqrt{1-x^{2}}}=\sum_{n=1}^{\infty}W_{2n-1}x^{2n-1}}$}\]

W_{q} désigne, pour tout q\in\mathbb{N}, la q-ème intégrale de Wallis. On rappelle que, par définition :

    \[W_{q}=\int_{0}^{\pi/2}\cos^{q}\left(t\right)\thinspace dt\]

Cette propriété sera établie à l’annexe 2.

Etant donné que :

    \[\frac{d}{dx}\left(\frac{1}{2}\arcsin^{2}\left(x\right)\right)=\frac{\arcsin\left(x\right)}{\sqrt{1-x^{2}}}\]

il en résulte que, pour tout x\in\left]-1,1\right[ :

    \[\arcsin^{2}\left(x\right)=\sum_{n=1}^{\infty}\frac{W_{2n-1}}{n}x^{2n}\]

Or, d’après un calcul classique :

    \[\frac{W_{2n-1}}{n}=\frac{1}{n}\thinspace\frac{4^{n-1}\left[\left(n-1\right)!\right]^{2}}{\left(2n-1\right)!}\]

et donc, en multipliant numérateur et dénominateur par 4n^{2} :

    \[\frac{W_{2n-1}}{n}=\frac{4^{n}}{2n^{2}\binom{2n}{n}}\]

Ainsi, pour tout x\in\left]-1,1\right[ :

    \[\arcsin^{2}\left(x\right)=\sum_{n=1}^{\infty}\frac{\left(2x\right)^{2n}}{2n^{2}\binom{2n}{n}}\]

ou, de manière équivalente :

    \[\boxed{\forall x\in\left]-2,2\right[,\thinspace2\arcsin^{2}\left(\frac{x}{2}\right)=\sum_{n=1}^{\infty}\frac{x^{2n}}{n^{2}\binom{2n}{n}}}\]

Il ne reste plus qu’à choisir x=1, ce qui donne (compte tenu de l’égalité \arcsin\left(\frac{1}{2}\right)=\frac{\pi}{6}) :

    \[\frac{\pi^{2}}{18}=\sum_{n=1}^{\infty}\frac{1}{n^{2}\binom{2n}{n}}\]

La formule annoncée est donc établie.

En termes d’accélération de la convergence, le gain est spectaculaire ! Si l’on note :

    \[A_{n}=\sum_{k=1}^{n}\frac{1}{k^{2}}\qquad\text{et}\qquad B_{n}=\sum_{k=1}^{n}\frac{3}{k^2\binom{2k}{k}}\]

on constate que :

    \[\frac{\pi^{2}}{6}-A_{100}\simeq0,01\]

tandis que :

    \[\frac{\pi^{2}}{6}-B_{100}\simeq10^{-63}\]

Annexe 1 – Preuve de la 2ème version de la formule de Legendre

Nous démontrons ici la :

Proposition

Pour tout entier n\geqslant1 et tout nombre premier p :

    \[v_{p}\left(n!\right)=\frac{n-S_{p}\left(n\right)}{p-1}\]

S_{p}\left(n\right) désigne la somme des chiffres de l’écriture de n en base p.

Commençons par exprimer n en base p :

    \[n=\sum_{k=0}^{d}\,a_{k}p^{k}\]

avec \forall k\in\left\llbracket 0,d\right\rrbracket ,\thinspace a_{k}\in\left\llbracket 0,p-1\right\rrbracket et a_{d}\neq0.

On observe que pour tout i\in\left\llbracket 1,d\right\rrbracket :

(\star)   \[\left\lfloor \frac{n}{p^{i}}\right\rfloor =\left\lfloor \sum_{k=i}^{d}\,a_{k}p^{k-i}+\sum_{k=0}^{i-1}\,a_{k}p^{k-i}\right\rfloor =\sum_{k=i}^{d}\,a_{k}p^{k-i}\]

car :

    \[\sum_{k=i}^{d}\,a_{k}p^{k-i}\in\mathbb{N}\]

et

    \[0\leqslant\sum_{k=0}^{i-1}\,a_{k}p^{k-i}\leqslant\left(p-1\right)\sum_{k=0}^{i-1}p^{k-i}<\left(p-1\right)\sum_{j=1}^{\infty}\frac{1}{p^{j}}=1\]

Donc, d’après la formule de Legendre et la formule \left(\star\right) :

    \begin{eqnarray*}v_{p}\left(n!\right) & = & \sum_{i=1}^{d}\left\lfloor \frac{n}{p^{i}}\right\rfloor\\& = & \sum_{i=1}^{d}\left(\sum_{k=i}^{d}\,a_{k}p^{k-i}\right)\\& = & \sum_{k=1}^{d}\left(a_{k}\sum_{i=1}^{k}p^{k-i}\right)\\& = & \sum_{k=1}^{d}a_{k}\frac{p^{k}-1}{p-1}\end{eqnarray*}

Or :

    \[S_{p}\left(n\right)=\sum_{k=0}^{d}\,a_{k}\]

et donc :

    \[\boxed{v_{p}\left(n!\right)=\frac{n-S_{p}\left(n\right)}{p-1}}\]

comme souhaité.

Annexe 2 – Un développement en série entière peu classique

Lemme

Pour tout x\in\left]-1,1\right[ :

    \[\frac{\arcsin\left(x\right)}{\sqrt{1-x^{2}}}=\sum_{n=0}^{\infty}W_{2n+1}x^{2n+1}\]

Partons de la série en utilisant la définition des intégrales de Wallis.

Pour tout x\in\left]-1,1\right[ :

    \begin{eqnarray*}\sum_{n=0}^{\infty}W_{2n+1}x^{2n+1} & = & \sum_{n=0}^{\infty}x^{2n+1}\int_{0}^{\pi/2}\cos^{2n+1}\left(t\right)\thinspace dt\\& = & \int_{0}^{\pi/2}\frac{x\cos\left(t\right)}{1-x^{2}\cos^{2}\left(t\right)}\thinspace dt\end{eqnarray*}

L’intégration terme à terme est justifiée car la série d’applications {\displaystyle \sum_{n\geqslant0}\left[t\mapsto x^{2n+1}\cos^{2n+1}\left(t\right)\right]} converge normalement (et donc uniformément) sur le segment \left[0,\frac{\pi}{2}\right].

En remplaçant \cos^{2}\left(t\right) par 1-\sin^{2}\left(t\right) au dénominateur puis en posant u=x\sin\left(t\right), il vient :

    \begin{eqnarray*}\sum_{n=0}^{\infty}W_{2n+1}x^{2n+1} & = & \int_{0}^{\pi/2}\frac{x\cos\left(t\right)}{1-x^{2}+x^{2}\sin^{2}\left(t\right)}\thinspace dt\\& = & \int_{0}^{x}\frac{du}{1-x^{2}+u^{2}}\\& = & \frac{1}{\sqrt{1-x^{2}}}\arctan\left(\frac{x}{\sqrt{1-x^{2}}}\right)\end{eqnarray*}

Pour finir, il est facile de voir que pour tout x\in\left]-1,1\right[ :

    \[\arctan\left(\frac{x}{\sqrt{1-x^{2}}}\right)=\arcsin\left(x\right)\]

(si vous ne voyez pas pourquoi, n’hésitez pas à demander !) Ainsi :

    \[\boxed{\forall x\in\left]-1,1\right[,\thinspace\sum_{n=0}^{\infty}W_{2n+1}x^{2n+1}=\frac{\arcsin\left(x\right)}{\sqrt{1-x^{2}}}}\]

Vos questions ou remarques sont les bienvenues. Vous pouvez laisser un commentaire ci-dessous ou bien passer par le formulaire de contact.

Partager cet article

La publication a un commentaire

  1. Oncle Junior

    Bonjour,
    Merci pour cet article instructif, notamment les deux formules de Legendre (ou plutôt la formule de Legendre et son autre formulation équivalente de l’annexe 1).

    Concernant les 4 formules de la partie 3, il me semble que la première a des applications en Dénombrement / Probabilité. En revanche, je ne connais pas d’application des formules 2-3-4 (bien que leur esthétique suffise, une application n’étant pas du tout nécessaire !).

    Bien à vous.

Laisser un commentaire