Probabilité de toucher la fève en découpant une galette des rois

En janvier, lors d’un goûter entre amis autour d’une galette des rois, quelqu’un m’a demandé quelle est la probabilité de toucher la fève lors de la découpe de la galette en parts égales.

J’ai fait une recherche rapide avec Google (et ChatGPT) et j’ai vu que la solution était donnée sous des hypothèses beaucoup trop restrictives, comme supposer connue la distance de la fève au centre de la galette.

Je donne ici une solution avec le minimum d’hypothèses permettant d’exprimer la probabilité sous une forme analytique.

Terminale scientifique ou, plutôt, classes préparatoires ou premier cycle universitaire.

Formulation mathématique du problème

Nous allons supposer que :

  • la fève et la galette sont des disques de rayons RFR_F et RGR_G,
  • la galette est découpée en n parts égales selon des demi-droites d’origine le centre OO de la galette,
  • le centre F de la fève est placé selon une loi de probabilité uniforme dans le disque de centre OO et de rayon RGRFR_G - R_F (sinon, une partie de la fève serait à l’extérieur de la galette).

Soit λ=RGRF.\lambda = \frac{R_{G}}{R_{F}}. Nous allons montrer que la probabilité cherchée est une fonction de λ et de n, soit P(λ,n)P(\lambda,n).

Il est à noter que les hypothèses précédentes (en gras) pourraient être relaxées. Néanmoins, elles ont l’avantage de permettre d’exprimer P(λ,n)P(\lambda, n) comme une fonction analytique de λ\lambda et de nn.

Titre par défaut

Il est clair que les hypothèses précédentes entraînent que P(λ,n)=1P(\lambda,n) = 1 si λ2λ \leqslant 2. Dans ce qui suit, nous allons donc supposer n2n \geqslant 2 (si n=1n=1, il n’y a pas de découpe !) et λ2λ \geqslant 2.

Nous allons considérer deux cas selon que F est à une distance R de O inférieure à RFR_F (cas 1) ou supérieure à RFR_F (cas 2).

Sans perte de généralité, on peut supposer que les angles des demi-droites sont égaux à 2kπn,k{0,1,...,n1}\frac{2k\pi}{n}, k \in \{ 0,1,...,n - 1\}.

Cas 1 :

Dans ce cas, O appartient à la fève et donc le couteau va nécessairement toucher la fève puisque la découpe débute en O. Soit p1p_1 la probabilité que ce cas arrive. Comme la fève est placée uniformément dans le disque de centre O et de rayon RGRFR_G - R_F, on a :

p1=πRF2π(RGRF)2=1(λ1)2(1)\tag{1} p_{1} = \frac{\pi R_{F}^{2}}{\pi\left( R_{G} - R_{F} \right)^{2}} = \frac{1}{(\lambda - 1)^{2}}

Cas 2 :

La probabilité d’occurrence du cas 2 est :

p2=1p1=λ(λ2)(λ1)2(2)\tag{2} p_{2} = 1 - p_{1} = \frac{\lambda(\lambda - 2)}{(\lambda - 1)^{2}}

Dans ce cas, soit α l’angle sous lequel la fève est vue depuis O, on a :

Angle α sous lequel la fève est vue de O
Figure 1 : Angle α\alpha sous lequel la fève est vue de OO.

Le triangle OFP est rectangle en P, par conséquent on a:

sin(α2)=RFRα=2arcsin(RFR)\sin\left( \frac{\alpha}{2} \right) = \frac{R_{F}}{R} \quad \Rightarrow \alpha = 2\arcsin\left( \frac{R_{F}}{R} \right)

Soient (X,Y)(X, Y) et (Θ,R)(\Theta,R) les coordonnées cartésiennes et polaires de F:

X=RcosΘ,Y=RsinΘ.X = R \cos\Theta, Y = R\sin \Theta.

Dans le cas 2, FF est uniformément réparti dans la couronne de centre OO et de rayons RFR_F et RGRFR_G - R_F. La densité de probabilité de (X,Y)(X, Y) est alors

p(X,Y)(x,y)=1aire(couronne)=1π((RGRF)2RF2)(3)\tag{3}p_{(X,Y)} (x, y) = \frac{1}{\text{aire(couronne)}} = \frac{1}{\pi\left((R_G - R_F)^2-R_F^2\right)}

Sachant qu’un élément d’aire dxdydx dy est égal à rdrdθrdrd\theta, la densité de probabilité p(Θ,R)p_{(\Theta, R)} de (Θ,R)(\Theta, R) s’écrit alors :

p(Θ,R)(θ,r)=rπ((RGRF)2RF2)RFr(RGRF)0θ2π(4)\tag{4} p_{(\Theta, R)}(\theta,r) = \frac{r}{\pi\left( \left( R_{G} - R_{F} \right)^{2} - R_{F}^{2} \right)}\quad\quad R_{F} \leq r \leq \left( R_{G} - R_{F} \right)\quad\quad 0 \leq \theta \leq 2\pi

La densité de probabilité de pRp_R de RR s’écrit alors :

pR(r)=02πp(θ,r)dθ=2r(RGRF)2RF2RFr(RGRF)(5)\tag{5} p_R(r) = \int_{0}^{2\pi} p(\theta,r) \mathrm d\theta = \frac{2r}{\left( R_{G} - R_{F} \right)^{2} - R_{F}^{2}}\quad\quad R_{F} \leq r \leq \left( R_{G} - R_{F} \right)

D’après la Figure (1), la probabilité qu’une demi-droite D d’angle φ\varphi touche la fève lorsque la fève est à une distance R=rR = r de O est égale à la probabilité que Θ\Theta appartienne, modulo 2π, au secteur angulaire [φα(r)2,φ+α(r)2]\lbrack \varphi - \frac{\alpha(r)}{2},\varphi + \frac{\alpha(r)}{2} \rbrack, soit α(r)2π\frac{\alpha(r)}{2\pi}.

Par conséquent, la fève sera touchée lors de la découpe si, et seulement si, Θ\Theta appartient, modulo 2π, à l’un des secteurs angulaires {[2kπnα(r)2,2kπn+α(r)2],k=0,1,...,n1}\{\left\lbrack \frac{2k\pi}{n} - \frac{\alpha(r)}{2},\frac{2k\pi}{n} + \frac{\alpha(r)}{2} \right\rbrack, k = 0,1,...,n - 1\}. Si R=rR = r, la probabilité d’intersection avec la fève est égale à :

Pinter(r)=min(1,nα(r)2π)=min(1,nπarcsin(RFr))(6)\tag{6} P_{\text{inter}}(r) = \min\left( 1,\frac{n\alpha(r)}{2\pi} \right) = \min\left( 1,\frac{n}{\pi}\text{arcsin}\left( \frac{R_{F}}{r} \right) \right)

La probabilité est égale à 1 si rRFsin(π/n)r \leq \frac{R_{F}}{\sin(\pi/n)}. Sachant que rRGRF,r \leq R_{G} - R_{F}, la probabilité d’intersection dans le cas 2 est donc égale à 1 si λ1+1/sin(π/n).\lambda \leq 1 + 1/\sin(\pi/ n). On a donc :

P(λ,n)=1siλ1+1sin(π/n)(7)\tag{7} P(\lambda,n) = 1 \quad \text{si} \quad \lambda \leq 1 + \frac{1}{\sin(\pi/n)}

En particulier, (7) implique que la fève sera touchée asymptotiquement avec une probabilité 1 lorsque le nombre de parts est de plus en plus grand.

La probabilité de toucher la fève lors de la découpe, s’obtient en randomisant R dans (5) :

P(intersection | Cas 2)=RFRGRFmin(1,nπarcsin(RFr))p(r)dr(8)\tag{8} P\left( \text{intersection | Cas 2} \right) = \int_{R_{F}}^{R_{G} - R_{F}}{\min\left( 1,\frac{n}{\pi}\arcsin\left( \frac{R_{F}}{r} \right) \right) p(r) \mathrm dr}

En reportant dans (5), on obtient :

P(intersection | Cas 2)=2(RGRF)2RF2RFRGRFmin(1,nπarcsin(RFr))rdr(9)\tag{9} P\left( \text{intersection | Cas 2} \right) = \frac{2}{\left( R_{G} - R_{F} \right)^{2} - R_{F}^{2}}\int_{R_{F}}^{R_{G} - R_{F}}{\min\left( 1,\frac{n}{\pi}\arcsin\left( \frac{R_{F}}{r} \right) \right)} r \mathrm dr

Faisons le changement de variable x=RFrx = \frac{R_{F}}{r}, on a alors :

RFRGRFmin(1,nπarcsin(RFr))rdr=RF21(λ1)1min(1,nπarcsin(x))x3dx(10)\tag{10} \int_{R_{F}}^{R_{G} - R_{F}}{\min\left( 1,\frac{n}{\pi}\arcsin\left( \frac{R_{F}}{r} \right) \right)}r \mathrm dr = R_{F}^{2}\int_{\frac{1}{(\lambda - 1)}}^{1}{\frac{\min\left( 1,\frac{n}{\pi}\arcsin(x) \right)}{x^{3}} \mathrm dx}

Par conséquent, on peut écrire :

P(intersection | Cas 2)=2nπλ(λ2)(1(λ1)sin(π/n)arcsin(x)x3dx)+2λ(λ2)sin(π/n)1dxx3(11)\tag{11} P\left( \text{intersection | Cas 2} \right) = \frac{2n}{\pi\lambda(\lambda - 2)}\left( \int_{\frac{1}{(\lambda - 1)}}^{\sin(\pi/n)}{\frac{\arcsin(x)}{x^{3}}} \mathrm dx \right) + \frac{2}{\lambda(\lambda - 2)}\int_{\sin(\pi/n)}^{1}{\frac{\mathrm dx}{x^{3}}}

La première intégrale se calcule analytiquement :

arcsin(x)x3dx=12arcsin(x)x2121x2x+C(12)\tag{12} \int_{}^{}{\frac{\arcsin(x)}{x^{3}} \mathrm dx =} - \frac{1}{2}\frac{\arcsin(x)}{x^{2}} - \frac{1}{2}\frac{\sqrt{1 - x^{2}}}{x} + C

Or, la probabilité cherchée s’écrit :

P(λ,n)=P(Cas 1)+P(intersection | Cas 2)P(Cas 2)(13)\tag{13} P(\lambda,n) = P(\text{Cas 1}) + P\left( \text{intersection | Cas 2} \right) P(\text{Cas 2})

Tous calculs faits, on obtient :

P(λ,n)=nπ(arcsin(1λ1)+λ(λ2)cotan(πn)(λ1)2)χ(λ)+(1χ(λ))(14)\tag{14} \boxed{P(\lambda,n) = \frac{n}{\pi}\left( \text{arcsin}\left( \frac{1}{\lambda - 1} \right) + \frac{\sqrt{\lambda(\lambda - 2)} - \text{cotan}\left( \frac{\pi}{n} \right)}{(\lambda - 1)^{2}} \right)\chi(\lambda) + \left( 1 - \chi(\lambda) \right)}

χ(λ)\chi(\lambda) est égale à 1 si λ1+cosec(π/n)\lambda \geqslant 1 + \cosec(\pi / n) et 0 sinon.

Application numérique

Les graphes de P(λ,4)P(λ, 4), P(λ,6)P(λ, 6) et P(λ,8)P(λ, 8) en fonction de λ\lambda sont montrés Figure 2.

Ainsi, si la galette a un rayon 10 fois plus grand que la fève :

P(10,4)0,267P(10,6)0,383P(10,8)0,489P(10, 4) \approx 0,267 \quad P(10, 6) \approx 0,383 \quad P(10, 8) \approx 0,489
image.png
Figure 2 : Probabilité de toucher la fève par une découpe de la galette en parts égales.

3 commentaires

Ça me disait quelque chose, alors j’ai cherché et j’ai retrouvé cet article d’il y a 5 ans (déjà ! :waw: ).

Du coup je remets mon commentaire de l’époque qui s’applique à nouveau :

J’dis ça j’dis rien, mais dans la boulangerie où je vais, les galettes pour 6 ou 8 personnes ont 2 fèves dedans… :-°

+4 -0
Connectez-vous pour pouvoir poster un message.
Connexion

Pas encore membre ?

Créez un compte en une minute pour profiter pleinement de toutes les fonctionnalités de Zeste de Savoir. Ici, tout est gratuit et sans publicité.
Créer un compte