Le terme

est obtenu par

, puis

, puis

, etc., autrement dit

avec

venant du développement de (1+x)^n et

venant du développement de
^m)
Or le coefficient binomial en facteur de

est

dans le développement de
^n)
Donc le coefficient qui multiplie

est

, autrement dit pour l'indice k on trouve dans le développement le terme suivant

Donc pour r fixé, en mettant

en facteur on trouve bien

Maintenant si on appelle

, la somme de k=0 à r est égale à la somme de

ce qui fait qu'on a bien pour l'indice r le terme suivant

Il reste à déterminer l'ensemble des r possibles : or
^n)
comporte tous les exposants de x compris entre 0 et n, et
^m)
, entre 0 et m. Les exposants s'ajoutant quand on multiplie, cela fait dans le développements des exposants allant de 0+0 à m+n : ce sont donc bien les bornes de variations de r et on obtient ainsi
^n(1+x)^m=\sum_{r=0}^{m+n}\Big(\sum_{k+p=r} \begin{pmatrix} n \\ k \end{pmatrix}\begin{pmatrix} m \\ p \end{pmatrix}x^r\Big))
Il n'y a que 10 types de personne au monde : ceux qui comprennent le binaire et ceux qui ne le comprennent pas.