3 / 11

Esercizio: probabilità nel gioco del poker e paradosso dei compleanni

Calcolo delle probabilità delle mani del poker a 5 carte tramite combinatoria, e analisi del paradosso dei compleanni con approssimazione di Stirling e probabilità condizionata.

Esercizio

Calcoliamo la probabilità delle varie combinazioni nel gioco di poker a 5 carte. Il mazzo ha 52 carte, e l'estrazione di 5 carte senza rinserimento è equiprobabile, possiamo usare P(A)=AΩP(A) = \frac{|A|}{|\Omega|} dove Ω=(525)|\Omega| = \binom{52}{5}, e A|A| è il numero di combinazioni che soddisfano la condizione che stiamo cercando.

  1. A1A_1 = una coppia
  2. A2A_2 = due coppie
  3. A3A_3 = tris
  4. A4A_4 = scala (non reale)
  5. A5A_5 = colore
  6. A6A_6 = full
  7. A7A_7 = poker
  8. A8A_8 = scala reale

A1A_1, una coppia

Un esempio è se in una mano ci sono 2 carte dello stesso valore (es. un 4 di picche e un 4 di cuori) e 3 carte di valore diverso (che NON includono un altra coppia)

  • Numero di modi di scegliere 1 valore su 13 (da 1 a 13), quindi (131)\binom{13}{1}
  • Numero di modi di scegliere 2 semi diversi sui 4 disponibili, contando anche il fatto che l'ordine non conta, allora (42)\binom{4}{2}
  • Numero di carte rimanenti, dobbiamo evitare di formare un tris, un altra coppia o un full, abbiamo già una coppia quindi 12 valori rimanenti delle carte e quindi dobbiamo scegliere 3 valori diversi tra i 12 rimasti, allora (123)\binom{12}{3}
  • Per ogni valore scelto delle 3 carte, abbiamo 4 semi disponibili da assegnare ad ogni carta quindi 444=434\cdot4\cdot4 = 4^3

Riassumendo tutto otteniamo che A1=(131)(42)(123)43|A_1| = \binom{13}{1} \cdot \binom{4}{2} \cdot \binom{12}{3} \cdot 4^3, e quindi P(A1)=A1ΩP(A_1) = \frac{|A_1|}{|\Omega|}.

A2A_2, due coppie

Un esempio è se in una mano ci sono 2 carte dello stesso valore (es. un 4 di picche e un 4 di cuori), altre 2 carte dello stesso valore (es. un 7 di fiori e un 7 di quadri) e una carta di valore diverso.

  • Numero di modi di scegliere 2 valori su 13 per formare 2 coppie di valori diversi (132)\binom{13}{2}
  • Per ogni coppia, numero di modi di scegliere 2 semi diversi sui 4 disponibili, allora (42)(42)\binom{4}{2} \cdot \binom{4}{2}
  • Ultima carta rimanente escludendo il full (111)\binom{11}{1}
  • Per l'ultima carta, abbiamo 4 semi disponibili da assegnare, quindi (41)\binom{4}{1}

Riassumendo otteniamo che A2=(132)(42)(42)(111)(41)|A_2| = \binom{13}{2} \cdot \binom{4}{2} \cdot \binom{4}{2} \cdot \binom{11}{1} \cdot \binom{4}{1}, e quindi P(A2)=A2ΩP(A_2) = \frac{|A_2|}{|\Omega|}.

A3A_3, tris

Analogo al caso con la coppia ma usando 3 carte dello stesso valore, e 2 carte di valore diverso, allora A3=(131)(43)(122)42|A_3| = \binom{13}{1} \cdot \binom{4}{3} \cdot \binom{12}{2} \cdot 4^2, e quindi P(A3)=A3ΩP(A_3) = \frac{|A_3|}{|\Omega|}.

A4A_4, scala (non reale)

Per formare una scala, dobbiamo scegliere 5 valori consecutivi. Ce ne sono 10 dato che abbiamo A-2-3-4-5, 2-3-4-5-6, 3-4-5-6-7, ..., 10-J-Q-K-A

  • Numero di sequenze possibili (101)=10\binom{10}{1} = 10
  • Numero di semi per ogni carta della scala, 4 per ogni carta, 5 carte in totale quindi 454^5
  • Dobbiamo però togliere i 4 casi in cui le carte hanno tutti lo stesso seme (sennò sarebbe reale)

Riassumendo otteniamo che A4=(101)(454)|A_4| = \binom{10}{1} \cdot (4^5 - 4), e quindi P(A4)=A4ΩP(A_4) = \frac{|A_4|}{|\Omega|}.

A5A_5, colore

Ovvero 5 carte tutte dello stesso seme, ignorando il colore

  • Numero di carte da scegliere 5 su 13 (135)\binom{13}{5}
  • Numero di semi da scegliere 1 su 4 (41)\binom{4}{1}
  • Da rimuovere tutte le configurazioni che portano ad una scala reale, cioè (101)(41)\binom{10}{1} \cdot \binom{4}{1}

Riassumendo otteniamo che A5=(135)(41)(101)(41)|A_5| = \binom{13}{5} \cdot \binom{4}{1} - \binom{10}{1} \cdot \binom{4}{1}, e quindi P(A5)=A5ΩP(A_5) = \frac{|A_5|}{|\Omega|}.

A6A_6, full

Un esempio è se in una mano ci sono 3 carte dello stesso valore (es. un 4 di picche, un 4 di cuori e un 4 di fiori) e altre 2 carte dello stesso valore (es. un 7 di fiori e un 7 di quadri). Dobbiamo combinare una coppia con un tris:

  • Numero di modi di scegliere 1 valore su 13 (da 1 a 13), quindi (131)\binom{13}{1}
  • Numero di modi di scegliere 2 semi diversi sui 4 disponibili, allora (42)\binom{4}{2}
  • Numero di modi di scegliere 1 valore su 12 (12 perchè uno è già usato nella coppia) (121)\binom{12}{1}
  • Numero di modi di scegliere 3 semi diversi sui 4 disponibili, allora (43)\binom{4}{3}

Riassumendo otteniamo che A6=(131)(42)(121)(43)|A_6| = \binom{13}{1} \cdot \binom{4}{2} \cdot \binom{12}{1} \cdot \binom{4}{3}, e quindi P(A6)=A6ΩP(A_6) = \frac{|A_6|}{|\Omega|}.

A7A_7, poker

Un esempio è se in una mano ci sono 4 carte dello stesso valore (es. un 4 di picche, un 4 di cuori, un 4 di fiori e un 4 di quadri) e una carta di valore diverso.

  • Numero di modi di scegliere 1 valore su 13 (da 1 a 13), quindi (131)\binom{13}{1}
  • Numero di modi di scegliere i restanti 3 semi dello stesso valore, 1 ogniuno quindi 131^3
  • Numero di modi di scegliere 4 semi su 4, (44)=1\binom{4}{4} = 1
  • Numero di modi di scegliere 1 valore su 12 per l'ultima carta (121)\binom{12}{1}
  • Numero di modi di scegliere 1 seme su 4 per l'ultima carta (41)\binom{4}{1}

Riassumendo otteniamo che A7=(131)13(44)(121)(41)|A_7| = \binom{13}{1} \cdot 1^3 \cdot \binom{4}{4} \cdot \binom{12}{1} \cdot \binom{4}{1}, e quindi P(A7)=A7ΩP(A_7) = \frac{|A_7|}{|\Omega|}.

A8A_8, scala colore

Un esempio è se in una mano ci sono 5 carte dello stesso seme (es. tutte di picche) e con valori consecutivi (es. A-2-3-4-5). escluso il caso di scala reale, cioè A-K-Q-J-10

  • Numero di sequenze possibili, rimuovendo la reale (91)=9\binom{9}{1} = 9
  • Numero di modi di scegliere 1 seme su 4 (41)\binom{4}{1}

Riassumendo otteniamo che A8=(91)(41)|A_8| = \binom{9}{1} \cdot \binom{4}{1}, e quindi P(A8)=A8ΩP(A_8) = \frac{|A_8|}{|\Omega|}.

A9A_9, scala reale

C'è un solo modo per fare scala reale, cioè A-K-Q-J-10, e 4 modi per scegliere il seme, quindi A9=(41)|A_9| = \binom{4}{1}, e quindi P(A9)=A9ΩP(A_9) = \frac{|A_9|}{|\Omega|}.

Paradosso dei compleanni

Date kk persone, quale è la probabilità che due persone abbiano il compleanno nello stesso giorno?

Vogliamo sapere per quale valore kk la probabilità è maggiore di 0.5, cioè P(A)>0.5P(A) > 0.5.

I casi totali sono Ω=365k|\Omega| = 365^k, dato che ogni persona può avere uno dei 365 giorni come compleanno.

A questo punto possiamo calcolare il suo complemento, trovare questo risultato è equivalente a trovare sottrarre ad 1 la probabilità che tutte le persone abbiano compleanni diversi.

Che tutti abbiano i compleanni diversi, scegliamo per ogni persona uno dei giorni rimanenti:

365364363(365k+1)365 \cdot 364 \cdot 363 \cdot \dots \cdot (365 - k + 1)

Quindi otteniamo:

1365364(365k+1)(365)k1 - \frac{365 \cdot 364 \cdot \dots \cdot (365 - k + 1)}{(365)^k}

Che possiamo riscrivere usando la produttoria:

=1i=0k1365i365= 1 - \prod_{i=0}^{k-1} \frac{365-i}{365}

Separando il primo termine (i=0i=0) che è uguale a 1, e riscrivendo la frazione:

=1i=1k1(1i365)= 1 - \prod_{i=1}^{k-1} \left(1 - \frac{i}{365}\right)

Supponendo che k365k \ll 365 (cosa che verificheremo alla fine), allora i3651\frac{i}{365} \ll 1. Possiamo usare l'approssimazione di Taylor ex1+xe^x \approx 1+x (o meglio ln(1+x)x\ln(1+x) \approx x per x1x \ll 1), quindi 1i365ei3651 - \frac{i}{365} \approx e^{-\frac{i}{365}}.

1i=1k1ei365\cong 1 - \prod_{i=1}^{k-1} e^{-\frac{i}{365}}

Sfruttando le proprietà degli esponenziali (prodotto di potenze con stessa base è la somma degli esponenti):

=1ei=1k1i365=1e1365i=1k1i= 1 - e^{-\sum_{i=1}^{k-1} \frac{i}{365}} = 1 - e^{-\frac{1}{365} \sum_{i=1}^{k-1} i}

Sapendo che la somma dei primi nn interi è n(n+1)2\frac{n(n+1)}{2}, qui n=k1n=k-1:

=1ek(k1)2365= 1 - e^{-\frac{k(k-1)}{2 \cdot 365}}

Ora supponiamo anche che k(k1)2365=730k(k-1) \ll 2 \cdot 365 = 730. In tal caso possiamo usare l'approssimazione eϵ1+ϵ+o(ϵ2)e^{\epsilon} \approx 1 + \epsilon + o(\epsilon^2).

1[1+(k(k1)2365)]=k(k1)2365\approx 1 - \left[ 1 + \left( - \frac{k(k-1)}{2 \cdot 365} \right) \right] = \frac{k(k-1)}{2 \cdot 365}

Vogliamo quindi trovare kk tale che:

Pkk(k1)730>12P_k \approx \frac{k(k-1)}{730} > \frac{1}{2}

Proviamo con alcuni valori:

  • Se k=191918730=342730<12k=19 \rightarrow \frac{19 \cdot 18}{730} = \frac{342}{730} < \frac{1}{2}
  • Se k=202019730=380730>12k=20 \rightarrow \frac{20 \cdot 19}{730} = \frac{380}{730} > \frac{1}{2}

Controlliamo l'ipotesi iniziale: 380<730380 < 730 (l'approssimazione è accettabile).

Formula di stirling

La formula di Stirling è un'approssimazione per i fattoriali, che dice che per n grande:

n!nnen2πnn! \approx n^n e^{-n}\sqrt{2 \pi n}

Probabilità condizionata

Prendiamo un dado come esempio, ha 6 facce e se non è truccato, ogni faccia è equiprobabile. Allora:

  • Ω={1,2,3,4,5,6}\Omega = \{ 1,2,3,4,5,6 \}
  • P(i)=16i{1,,6}P(i) = \frac{1}{6} \forall i \in \{1, \cdot, 6\} equiprobabile

A questo punto, aggiungiamo un informazione, cioè che il dado è truccato e il dado da solo numeri dispari.

Noi dobbiamo stabilire le leggi di probabilità in base alle informazioni che sappiamo del sistema, ogni volta che aggiungo un informazione su un evento che posso considerare "certo" (come in questo caso di avere solo le facciate dispari), posso dire di avere una probabilità condizionata, assumendo che quell'evento si sia relizzato.

Allora chiamiamo la nostra nuova probabilità condizionata come P~\tilde{P}, allora il suo nuovo insieme sarà

  • B={1,3,5}B = \{1,3,5\}
  • Bc={2,4,6}B^c = \{2,4,6\}
  • BBc=ΩB \cup B^c = \Omega

Allora definiamo P~\tilde{P}:

P~(A)=P(AB)\tilde{P}(A) = P(A | B)

tale che P~(BB)=1\tilde{P}(B | B) = 1

A questo punto vediamo che:

  • P(AB)P(A \cap B) = la probabilità che entrambi gli eventi accadano insieme
  • P(B)P(B) = la probabilità che l'evento certo accada
P(AB)=P(AB)P(B)P(A|B) = \frac{P(A \cap B)}{P(B)}

Che nel caso di eventi equiprobabili, è la cardinalità.

Note:

ABBA    P(AB)P(B)1|A \cap B| \le |B| \le |A| \implies \frac{P(A \cap B)}{P(B)} \le 1

Quindi la probabilità condizionata non è altro che la probabilità aggiungendo l'informazione che accade un evento certo

Probabilità totale

La probabilità totale, è la probabilità che si realizzino almeno uno tra due o più eventi

Partizioniamo Ω\Omega in ii partizioni: Ω=i=1nEi\Omega = \cup_{i=1}^n E_i con EiEj=E_i \cap E_j = \emptyset se iji \ne j.

Sappiamo che per identità:

A=AΩA = A \cap \Omega

Allora sostituiamo Ω\Omega con le sue partizioni:

P(A)=P(AΩ)=P(Ai=1nEi)=P(i=1n(AEi))P(A) = P(A \cap \Omega) = P\left(A \cap \bigcup_{i=1}^n E_i\right) = P\left(\bigcup_{i=1}^n (A \cap E_i)\right)

Dato che le partizioni sono disgiunte, allora anche le loro intersezioni lo sono, e allora per gli assiomi dell probabilità abbiamo che:

P(A)=i=1nP(AEi)=i=1nP(AEi)P(Ei)P(A) = \sum_{i=1}^n P(A \cap E_i) = \sum_{i=1}^n P(A | E_i) P(E_i)

In parole, trovare la probabilità totale di A è uguale a considerare quanto è probabile A nello scenario EiE_i, per poi moltiplicarlo per quanto è probabile che EiE_i accada.

Esempio

Ci sono 5 ristoranti: R1,R2,R3,R4,R5R_1, R_2, R_3, R_4, R_5 e ogniuno è aperto con probabilità 15\frac{1}{5}

3 di questi ristoranti si trovano a destra di un bivio, e 2 a sinistra

Quale è la probabilità che trascorrendo la strada io trovi un ristorante aperto?

P(A)=P(A1,S)+P(A2,S)+P(A3,D)+P(A4,D)+P(A5,D)==(1512)2+(1512)3=12\begin{aligned} P(A) &= P(A_1, S) + P(A_2, S) + P(A_3, D) + P(A_4, D) + P(A_5, D) = \\ &= \left(\frac{1}{5} \cdot \frac{1}{2}\right) \cdot 2 + \left(\frac{1}{5} \cdot \frac{1}{2}\right) \cdot 3 = \frac{1}{2} \end{aligned}

Però possiamo usare la probabilità totale, partizionando in base alla scelta di andare a sinistra o a destra:

P(Avado a sinistra)=25P(Avado a destra)=35\begin{aligned} P(A | \text{vado a sinistra}) &= \frac{2}{5} \\ P(A | \text{vado a destra}) &= \frac{3}{5} \end{aligned}

Mettendo tutto insieme usando la formula otteniamo

P(A)=P(Avado a sinistra)P(S)+P(Avado a destra)P(D)=12(35+25)=12\begin{aligned} P(A) &= P(A | \text{vado a sinistra})P(S) + P(A | \text{vado a destra})P(D) \\ &= \frac{1}{2}\left(\frac{3}{5} + \frac{2}{5} \right) \\ &= \frac{1}{2} \end{aligned}

Proprietà della probabilità condizionata

  • Se EFE \subset F, allora P(EF)=P(E)P(F)P(E | F) = \frac{P(E)}{P(F)}
  • Se P(EF)=1P(E | F) = 1, allora EFE \subseteq F
  • Se EF=E \cap F = \emptyset, allora P(EF)=0P(E | F) = 0
  • P(EcF)=1P(EF)P(E^c | F) = 1 - P(E | F)
  • P(F)=P(FE)+P(FEc)P(F) = P(F \cap E) + P(F \cap E^c)
  • Per due eventi E1,E2E_1, E_2 con E1E2E_1 \cap E_2 \ne \emptyset, allora P(E1E2F)=P(E1F)+P(E2F)P(E1E2F)P(E_1 \cup E_2 |F) = P(E_1 | F) + P(E_2 | F) - P(E_1 \cap E_2 | F)