Pour simuler avec une pièce une épreuve de Bernoulli de proba. de réussite
absolument quelconque la méthode théorique est on ne peut plus simple : on jette une infinité de fois la pièce, on note avec des 0 et des 1 la suite des résultats et on considère cette suite comme celle de l'écriture en base
d'un réel
. La loi de
est évidement la loi uniforme sur
et, concernant la Bernoulli à simuler, il suffit bien sûr de dire que "Réussite" <=>
pour avoir la bonne simulation.
Bon, évidement, on peut se dire que ça va être un peu chaud de jeter une infinité de fois la pièce, mais en fait, non vu qu'on peut s’arrêter dés que les digits de
(en base 2) déjà obtenues sont suffisant pour être sûr que
ou bien que
(par contre, le petit problème, c'est qu'il faut avoir l'écriture de
en base 2 et que, si
n'est pas rationnel, ça va pas être périodique)
Par exemple, pour simuler une épreuve de Bernoulli de proba de 1/3, on écrit 1/3 en base 2 :
1/3 = 0,010101010101... et donc on jette la pièce jusqu'au premier lancé ou la série de résultats obtenus diffère du début de la série 01010101... et on dit que la Bernoulli a réussi si la série obtenue correspond à l'écriture en base 2 d'un réel <1/3 donc ici, à condition d'avoir tiré 00 ou bien 01010...100 avec deux zéros à la fin. Par contre, si la série commence par 1 ou bien est de la forme 01010...1011 avec deux 1 à la fin, c'est "echec".
Et c'est avec cette méthode que j'ai "pondu" ma méthode consistant à attendre la première apparition de deux tirages identiques pour simuler un dés à 6 faces.