Formule du crible/Exemple d'application

Nous allons traiter dans ce chapitre un exemple de problème que l’on résout avec la formule du crible. Le lecteur est invité, après avoir bien compris cet exemple, à faire l'exercice 1 qui est similaire.

Énoncé.

Une urne contient n boules numérotées de 1 à n. On tire indéfiniment dans cette urne avec remise. Soit X la variable aléatoire qui prend pour valeur le numéro du tirage où pour la première fois toutes les boules ont été obtenues au moins une fois. Donner la loi de X.

Corrigé.

Soit r ∈ ℕ, nous devons calculer p(X=r).

Si r < n, nous avons de façon évidente p(X=r) = 0. On ne peut pas obtenir toutes les boules si on fait moins de n tirages.

Supposons maintenant r ≥ n.

Soit B_i l’événement : « La boule numéro i n’est pas apparue au cours des r premiers tirages. »

B_i ∩ B_j représentera l’événement : « Les boules numérotées i et j ne sont pas sorties dans les r premiers tirages. »

B_i ∪ B_j représentera l’événement : « Au moins une des deux boules numérotées i ou j n’est pas apparue au cours des r premiers tirages. »

Nous avons donc :

$\begin{array}{cc}p(X>r)& =p\left(\bigcup _{i=1}^n{B}_i\right)\\ & =\sum _{k=1}^n((-1{)}^{k+1}\sum _{1⩽i_1<i_2<\cdots <i_k⩽n}p\left(\bigcap _{j=1}^k{B}_{i_j}\right))\end{array}$

$p\left({\bigcap }_{j=1}^k{B}_{i_j}\right)$ représente la probabilité que les k boules :

$B_{i_1},B_{i_2},\cdots ,B_{i_k}$

ne soient pas sorties au cours des r premiers tirages (les autres peuvent être sorties ou non).

Par conséquent :

$p\left(\bigcap _{j=1}^k{B}_{i_j}\right)={\left(\frac{n-k}{n}\right)}^r$

Et donc :

$p(X>r)=\sum _{k=1}^n((-1{)}^{k+1}\sum _{1⩽i_1<i_2<\cdots <i_k⩽n}{\left(\frac{n-k}{n}\right)}^r)$

$\sum _{1⩽i_1<i_2<\cdots <i_k⩽n}{\left(\frac{n-k}{n}\right)}^r$ est une somme de termes constants. Il y en a autant que ce qu’il y a de façon de choisir k nombres dans l’intervalle [1;n], c’est-à-dire $(\genfrac{}{}{0ex}{}{n}{k})$.

On a donc :

$\sum _{1⩽i_1<i_2<\cdots <i_k⩽n}{\left(\frac{n-k}{n}\right)}^r=(\genfrac{}{}{0ex}{}{n}{k}){\left(\frac{n-k}{n}\right)}^r$

Par conséquent :

$p(X>r)=\sum _{k=1}^n((-1{)}^{k+1}(\genfrac{}{}{0ex}{}{n}{k}){\left(\frac{n-k}{n}\right)}^r)$

Nous avons alors :

$\begin{array}{cc}p(X=r)& =p(X>r-1)-p(X>r)\\ & =\sum _{k=1}^n((-1{)}^{k+1}(\genfrac{}{}{0ex}{}{n}{k}){\left(\frac{n-k}{n}\right)}^{r-1})-\sum _{k=1}^n((-1{)}^{k+1}(\genfrac{}{}{0ex}{}{n}{k}){\left(\frac{n-k}{n}\right)}^r)\\ & =\sum _{k=1}^n((-1{)}^{k+1}(\genfrac{}{}{0ex}{}{n}{k})[{\left(\frac{n-k}{n}\right)}^{r-1}-{\left(\frac{n-k}{n}\right)}^r])\\ & =\sum _{k=1}^n((-1{)}^{k+1}(\genfrac{}{}{0ex}{}{n}{k}){\left(\frac{n-k}{n}\right)}^{r-1}[1-\left(\frac{n-k}{n}\right)])\\ & =\sum _{k=1}^n((-1{)}^{k+1}(\genfrac{}{}{0ex}{}{n}{k}){\left(\frac{n-k}{n}\right)}^{r-1}\left(\frac{k}{n}\right))\\ & =\sum _{k=1}^n((-1{)}^{k+1}(\genfrac{}{}{0ex}{}{n-1}{k-1}){\left(\frac{n-k}{n}\right)}^{r-1})\end{array}$