Lotto Simulator

Matematica e tecnologia

La vera matematica, la statistica, gli algoritmi e la fisica dietro Lotto Simulator — presi dal codice che gira nell’app.

Apri Lotto Simulator

Come leggere questa pagina

Ogni formula qui sotto è copiata dall’implementazione, non da una brochure. Dove un componente è una regola deterministica o un’euristica anziché una formula, descriviamo l’algoritmo esattamente come viene eseguito, senza inventare matematica che non c’è.

Lotto Simulator è un ambiente di analisi, apprendimento e simulazione. Genera combinazioni con sedici metodi, le sottopone a un Tribunale matematico composto da dodici giurati, sei avvocati della difesa e un Giudice supremo, verifica l’equità degli storici delle estrazioni reali, analizza le coincidenze di calendario e riproduce un’urna fisica in 3D. Nulla di tutto ciò modifica la probabilità di un numero in un’estrazione ufficiale — la sezione 8 spiega perché, numeri alla mano.

Notazione usata in tutta la pagina

dove
NN
dimensione dell’insieme dei numeri principali (per esempio 50 nell’Eurojackpot)
kk
quanti numeri principali vengono estratti e giocati (5 nell’Eurojackpot)
B,  bB,\; b
dimensione dell’insieme separato dei numeri bonus e quanti numeri bonus si giocano (12 e 2 nell’Eurojackpot)
DD
numero di estrazioni passate nella finestra di analisi selezionata
cnc_n
quante volte il numero nn è uscito tra i numeri principali di quelle DD estrazioni
t=0,1,…t = 0, 1, \dots
indice dell’estrazione, dalla più recente: t=0t = 0 è l’ultima estrazione
(nr)\binom{n}{r}
coefficiente binomiale, il numero di modi di scegliere rr elementi tra nn

1. Architettura del sistema

Una sola applicazione, un solo codice sorgente: la stessa web app funziona nei browser desktop e mobile e, senza modifiche, dentro le app iOS e Android (WebView Capacitor). I calcoli pesanti e a pagamento girano sul server; tutto ciò che l’utente vede viene renderizzato sul dispositivo.

  1. Dispositivo (web, iOS, Android)Interfaccia, schedina e set di lavoro, i quattro metodi FREE del generatore, grafici delle frequenze, test χ², calcoli di calendario, urna 3D. Casuale puro e l’urna 3D attingono la casualità dalla Web Crypto API.
  2. Funzioni edge (TypeScript, Deno)I dodici metodi PRO del generatore, tutti i diciannove modelli del Tribunale, i dodici giurati, i sei avvocati, il Giudice supremo, le matrici della ruota, il registro delle previsioni. Il Tribunale lavora su flussi casuali riproducibili, inizializzati da un seme.
  3. Database (PostgreSQL)Risultati ufficiali delle estrazioni con la relativa versione delle regole, tabelle dei premi, account e diritti di accesso, il registro delle previsioni a sola aggiunta e la sua valutazione.
  4. Pipeline dei datiJob pianificati raccolgono i risultati ufficiali, convalidano ogni estrazione secondo le regole della sua epoca e pubblicano il set di dati con le regole attuali.

Chi calcola cosa

  • Nel browser: i metodi FREE (Numeri caldi, Combinato, Copertura del segmento, Casuale puro) e tutte le statistiche descrittive. La build pubblica rimuove fisicamente dal client le implementazioni PRO: le loro funzioni sono sostituite da stub, quindi un risultato PRO può arrivare solo dal server.
  • Sul server: ogni metodo PRO, l’intero Tribunale matematico e il Giudice supremo. Una richiesta contiene la lotteria, il metodo, il numero di righe (1–50) e la finestra di analisi; il server carica da sé le estrazioni, così il client non può fornirgli uno storico inventato.
  • Rifiuto anziché invenzione: un modello sul server richiede almeno 5 estrazioni nella finestra di analisi selezionata. Con meno estrazioni rifiuta con una risposta esplicita «storico insufficiente» invece di produrre righe dal nulla.

Calcolo riproducibile del Tribunale

Ogni giurato, avvocato e modello del Tribunale lavora sul proprio flusso pseudocasuale, il cui seme è un hash della versione dell’algoritmo, del personaggio, della lotteria, dello storico completo delle estrazioni e della data dell’estrazione di riferimento. Gli stessi dati producono quindi sempre lo stesso verdetto, e qualsiasi modifica dello storico cambia tutti i flussi.

h1=FNV1a⁡(H; h0),h2=FNV1a⁡(g ∥ d; h1),s=FNV1a⁡(v ∥ κ; h2),ui=mulberry32⁡(s)i∈[0,1)h_1 = \operatorname{FNV1a}(H;\ h_0), \quad h_2 = \operatorname{FNV1a}(g \,\|\, d;\ h_1), \quad s = \operatorname{FNV1a}(v \,\|\, \kappa;\ h_2), \qquad u_i = \operatorname{mulberry32}(s)_i \in [0,1)
Seme di un calcolo del Tribunale e flusso di numeri uniformi che ne deriva.
dove
vv
versione dell’algoritmo del motore del Tribunale
κ\kappa
chiave del personaggio o del modello, per esempio l’identificativo di un giurato
HH
tutte le estrazioni della finestra, numeri principali e bonus in ordine
g,  dg,\; d
identificativo della lotteria e data dell’estrazione di riferimento
h0h_0
la base di offset (offset basis) di FNV-1a; FNV1a⁡(x; h)\operatorname{FNV1a}(x;\,h) calcola l’hash del testo xx partendo dallo stato hh
uiu_i
l’ii-esimo numero casuale uniforme del flusso

Persistenza e provenienza

Ogni combinazione porta con sé un record di provenienza: la sua origine (metodo del generatore, giurato, consenso della giuria, ruota, estrazione 3D, Giudice supremo, inserimento manuale), i numeri di partenza e un registro di eventi a sola aggiunta con esami, proposte, posizioni degli avvocati, verdetti, decisioni e trasformazioni (al massimo 200 eventi). Un record viene accettato solo se, riapplicando le sue trasformazioni ai numeri iniziali, si ottiene esattamente la riga; altrimenti la riga viene mostrata come di origine non disponibile.

Rm=(τm∘⋯∘τ2∘τ1)(R0)R_m = \bigl(\tau_m \circ \cdots \circ \tau_2 \circ \tau_1\bigr)(R_0)
Verifica della provenienza: la riga attuale R_m deve coincidere con la riga iniziale R_0 a cui si applicano, in ordine, tutte le sostituzioni registrate τ₁ … τₘ.
  • Discendenza: quando una riga viene modificata a mano, si crea un nuovo record che ha come genitore il precedente (si conservano fino a quattro record precedenti), così la storia anteriore alla modifica manuale resta visibile.
  • Registro delle previsioni: ogni raccomandazione (una riga generata, la sostituzione di un giurato, la posizione di un avvocato, un verdetto) viene scritta in un registro prima dell’estrazione di riferimento e valutata automaticamente dopo il risultato ufficiale.

2. Generatore principale: sedici metodi

Il generatore offre sedici opzioni. Quattro metodi FREE girano nel browser; dodici metodi PRO girano sul server. Tutti scelgono kk numeri distinti da 1..N1..N più i numeri bonus della lotteria; differiscono solo nel modo in cui i numeri vengono pesati o costruiti.

Il Generatore di serie con la quantità, il metodo di analisi e i tre filtri predefiniti
Il pannello del generatore: selettore del metodo e filtri predefiniti

Elementi comuni

I metodi pesati assegnano un peso wn≥0w_n \ge 0 a ogni numero e poi estraggono in sequenza senza reinserimento: ogni scelta rimuove il numero e rinormalizza i pesi rimanenti.

P(xj=n∣x1,…,xj−1)=wn∑m∉{x1,…,xj−1}wmP\bigl(x_j = n \mid x_1,\dots,x_{j-1}\bigr) = \frac{w_n}{\displaystyle\sum_{m \notin \{x_1,\dots,x_{j-1}\}} w_m}
Campionamento pesato sequenziale senza reinserimento, usato da tutti i metodi basati su pesi.
dove
xjx_j
il jj-esimo numero scelto per la riga, j=1,…,kj = 1,\dots,k
wnw_n
il peso del numero nn secondo il metodo

Casualità: nel browser, Casuale puro e il passaggio di unicità usano la Web Crypto API, mentre Numeri caldi, Combinato e Copertura del segmento usano il generatore standard del browser (Math.random). Sul server i metodi attingono a casualità crittografica, tranne il Modello fisico (generatore standard) e il seme di Caos, che mescola casualità crittografica e un timer ad alta risoluzione. Nel browser un intero non distorto in 0..M−10..M-1 si ottiene con il campionamento per rifiuto, così nessun resto è favorito:

r=x mod M,x∼U{0,…,232−1} ∣ x<232−(232 mod M)r = x \bmod M, \qquad x \sim \mathcal{U}\{0,\dots,2^{32}-1\} \ \bigm|\ x < 2^{32} - \bigl(2^{32} \bmod M\bigr)
Campionamento per rifiuto: un valore crittografico x a 32 bit viene estratto di nuovo finché non scende sotto il più grande multiplo di M, poi viene ridotto modulo M.

Numeri bonus

I numeri bonus (Eurojackpot 2 su 12, EuroMillions 2 su 12, Powerball 1 su 26, Mega Ball 1 su 24, numero Viking 1 su 5, Powerball AU 1 su 20) vengono estratti dal loro insieme, indipendentemente dai numeri principali. Lotto, SuperEnalotto e Lotto Max non hanno un numero bonus per il giocatore: la loro pallina aggiuntiva viene estratta dall’operatore dall’urna principale.

wj′={cj′+0.05μ=freq1μ∈{rnd,phys,chaos,world}cj′+2μ∉{freq,rnd,phys,chaos,world}w'_j = \begin{cases} c'_j + 0.05 & \mu = \mathtt{freq} \\ 1 & \mu \in \{\mathtt{rnd}, \mathtt{phys}, \mathtt{chaos}, \mathtt{world}\} \\ c'_j + 2 & \mu \notin \{\mathtt{freq}, \mathtt{rnd}, \mathtt{phys}, \mathtt{chaos}, \mathtt{world}\} \end{cases}
Peso w′ⱼ del numero bonus j con conteggio bonus c′ⱼ per il metodo μ: freq = Numeri caldi; rnd, phys, chaos, world = Casuale puro, Modello fisico, Caos, metodi mondiali; tutti gli altri metodi usano c′ⱼ + 2. Con meno di 5 estrazioni ogni metodo usa pesi uniformi.

Filtri FREE e unicità

Le righe prodotte nel browser passano tre filtri attivi per impostazione predefinita: non tutti i numeri pari o tutti dispari; la somma entro un intervallo; nessuna sequenza di più di tre numeri consecutivi. Una riga scartata viene sostituita da una riga casuale uniforme che supera i filtri (i suoi numeri bonus vengono estratti con i pesi bayesiani cj′+2c'_j + 2).

Sˉ=k(N+1)2,Δ=max⁡ ⁣(18, round⁡(0.22 Nk)),S∈[max⁡(k,round⁡(Sˉ−Δ)), min⁡(kN,round⁡(Sˉ+Δ))]\bar S = \frac{k(N+1)}{2}, \qquad \Delta = \max\!\bigl(18,\ \operatorname{round}(0.22\,Nk)\bigr), \qquad S \in \bigl[\max(k, \operatorname{round}(\bar S - \Delta)),\ \min(kN, \operatorname{round}(\bar S + \Delta))\bigr]
Intervallo di somma predefinito. Eurojackpot: S̄ = 127.5, Δ = 55, intervallo 73–183.

Per evitare di consegnare due volte la stessa combinazione sullo stesso dispositivo, ogni riga FREE è identificata dal rango lessicografico dei suoi numeri principali tra tutte le (Nk)\binom{N}{k} combinazioni e confrontata con le ultime 20 000 righe emesse per quella lotteria; un rango ripetuto viene alterato con scambi casuali e, come ultima risorsa, si prende il successivo rango inutilizzato da una permutazione affine casuale di tutti i ranghi. Una riga modificata in questo modo non viene ricontrollata dai filtri.

rank⁡(x1<⋯<xk)=∑i=1k ∑v=xi−1+1xi−1(N−vk−i),x0=0;ρj=(a j+b) mod (Nk), gcd⁡ ⁣(a,(Nk))=1\operatorname{rank}(x_1 < \dots < x_k) = \sum_{i=1}^{k} \ \sum_{v = x_{i-1}+1}^{x_i - 1} \binom{N - v}{k - i}, \quad x_0 = 0; \qquad \rho_j = (a\,j + b) \bmod \binom{N}{k},\ \gcd\!\bigl(a, \tbinom{N}{k}\bigr) = 1
Rango lessicografico di una combinazione e percorso affine a periodo completo su tutti i ranghi.

Percorsi di esecuzione FREE e PRO

PercorsoMetodiDove giraRighe per esecuzioneStorico
FREENumeri caldi, Combinato, Copertura del segmento, Casuale purobrowserfino a 5 (fino a 50 per gli utenti PRO)la finestra FREE delle estrazioni recenti
PROgli altri dodici metodiserverfino a 50la finestra di analisi selezionata, almeno 5 estrazioni

Numeri caldi

FREENel browser

Scopo: favorire i numeri usciti più spesso nella finestra di analisi.

wn=cn+0.05w_n = c_n + 0.05
Peso del numero n; con meno di 5 estrazioni tutti i pesi valgono 1.
Dettagli avanzati

Combinato

FREENel browser

Scopo: mescolare numeri caldi, medi e freddi in una proporzione fissa.

I numeri vengono ordinati per cnc_n (a parità vince il numero più piccolo) e divisi in terzi di dimensione ⌈N/3⌉\lceil N/3 \rceil: caldi, medi, freddi. La riga ii prende scelte uniformi da ciascun terzo:

h=max⁡ ⁣(1,round⁡(0.4k)),a=⌊k−h2⌋,(mi,ci)={(a, k−h−a)i≡1(mod2)(k−h−a, a)i≡0(mod2)h = \max\!\bigl(1, \operatorname{round}(0.4k)\bigr), \qquad a = \Bigl\lfloor \frac{k-h}{2} \Bigr\rfloor, \qquad (m_i, c_i) = \begin{cases} (a,\ k-h-a) & i \equiv 1 \pmod 2 \\ (k-h-a,\ a) & i \equiv 0 \pmod 2 \end{cases}
La riga i prende h numeri caldi, mᵢ medi e cᵢ freddi; le righe dispari e pari si scambiano le quote dei medi e dei freddi.
Dettagli avanzati

Copertura del segmento

FREENel browser

Scopo: distribuire una riga in modo uniforme su tutto l’intervallo dei numeri.

s=⌊Nk⌋,xp∼U{ p s+1,…,(p+1) s }  (p<k−1),xk−1∼U{(k−1) s+1,…,N}s = \Bigl\lfloor \frac{N}{k} \Bigr\rfloor, \qquad x_p \sim \mathcal{U}\{\,p\,s + 1, \dots, (p+1)\,s\,\}\ \ (p < k-1), \qquad x_{k-1} \sim \mathcal{U}\{(k-1)\,s + 1, \dots, N\}
Un numero uniforme da ciascuno dei k segmenti consecutivi; l’ultimo segmento prende il resto.
Dettagli avanzati

Casuale puro

FREENel browser

Scopo: il riferimento onesto — ogni numero ugualmente probabile.

P(x1,…,xk)=(Nk)−1P(x_1,\dots,x_k) = \binom{N}{k}^{-1}
Un rimescolamento di Fisher–Yates crittografico dà a ogni combinazione la stessa probabilità prima del filtraggio.
Dettagli avanzati

Matrice della ruota

PROServer

Scopo: costruire più righe che insieme coprano il maggior numero possibile di coppie di un piccolo insieme di numeri — un disegno di copertura greedy.

L’insieme ha min⁡(N,k+5)\min(N, k+5) numeri presi alternativamente dai terzi caldo, freddo e medio. Tutte le (∣P∣k)\binom{|P|}{k} combinazioni dell’insieme sono candidate (oltre 40 000 si usa un campione casuale); a ogni passo si aggiunge la candidata con il punteggio migliore e se ne segnano le coppie come coperte:

score⁡(R)=100 ∣(R2)∖C∣  −  6∑n∈Run  −  2 [max⁡R−min⁡R<⌊max⁡P/3⌋]\operatorname{score}(R) = 100\,\Bigl|\tbinom{R}{2} \setminus \mathcal C\Bigr| \;-\; 6 \sum_{n \in R} u_n \;-\; 2\,\bigl[\max R - \min R < \lfloor \max P / 3 \rfloor\bigr]
Punteggio greedy di copertura delle coppie; le parità vengono risolte in modo uniformemente casuale.
Dettagli avanzati
dove
(R2), C\tbinom{R}{2},\ \mathcal C
le coppie della riga candidata RR e le coppie già coperte dalle righe scelte
unu_n
quante righe già scelte contengono il numero nn
[ ⋅ ][\,\cdot\,]
vale 1 se la condizione è vera, altrimenti 0

Caldo mondiale e Combinata mondiale

PROServer

Scopo: trasferire sulla lotteria attuale il profilo di frequenza di tutte e nove le lotterie supportate.

Ogni numero estratto xx di ogni lotteria gg viene ricondotto alla sua posizione relativa q=(x−1)/(Ng−1)q = (x-1)/(N_g - 1) e distribuito con un nucleo gaussiano stretto sui numeri della lotteria attuale; ogni lotteria riceve lo stesso peso per estrazione, e la frequenza propria della lotteria attuale viene aggiunta ancora una volta in più:

S(n)=10−4+∑g1Dg∑x∈Xgexp⁡ ⁣(−(qn−qx)22σ2)+1.8 cnD,qn=n−1N−1,  σ=0.035S(n) = 10^{-4} + \sum_{g} \frac{1}{D_g} \sum_{x \in X_g} \exp\!\Bigl(-\frac{(q_n - q_x)^2}{2\sigma^2}\Bigr) + 1.8\,\frac{c_n}{D}, \qquad q_n = \frac{n-1}{N-1},\ \ \sigma = 0.035
Punteggio mondiale del numero n; X_g sono tutti i numeri principali estratti nella lotteria g (fino a 400 estrazioni recenti con le regole attuali, D_g estrazioni).
Dettagli avanzati
  • Caldo mondiale: kk scelte uniformi tra i max⁡(3k,12)\max(3k, 12) numeri con il punteggio più alto.
  • Combinata mondiale: ⌈0.45k⌉\lceil 0.45k \rceil dal terzo superiore, max⁡(1,⌊0.2k⌋)\max(1, \lfloor 0.2k \rfloor) dal terzo inferiore, il resto da quello centrale. I numeri bonus sono uniformi.

Catene di Markov

PROServer

Scopo: pesare i numeri in base a quanto spesso hanno seguito i numeri dell’estrazione precedente.

wn=1∣L∣∑a∈LT(a→n)+1T(a)+N,T(a)=∑bT(a→b)w_n = \frac{1}{|L|} \sum_{a \in L} \frac{T(a \to n) + 1}{T(a) + N}, \qquad T(a) = \sum_{b} T(a \to b)
Probabilità di transizione del primo ordine con lisciamento di Laplace, mediata sull’ultima estrazione L.
Dettagli avanzati
dove
T(a→b)T(a \to b)
quante volte il numero bb è uscito nell’estrazione immediatamente successiva a un’estrazione contenente aa
LL
i numeri principali dell’ultima estrazione

Gauss · teorema del limite centrale

PROServer

Scopo: costruire righe con una somma tipica. La somma di kk numeri estratti senza reinserimento ha

μ=k(N+1)2,σ2=k(N+1)(N−k)12\mu = \frac{k(N+1)}{2}, \qquad \sigma^2 = \frac{k(N+1)(N-k)}{12}
Media e varianza della somma di una riga nel campionamento senza reinserimento (sostituite dai valori empirici quando esistono almeno 5 somme).
Dettagli avanzati

Ogni numero viene pesato con una gaussiana centrata sul valore medio per posizione, e si estraggono fino a 40 righe candidate; vince quella con la somma più vicina a μ\mu (la ricerca si ferma in anticipo entro 0.5σ0.5\sigma):

wn=exp⁡ ⁣(−(n−μ/k)22s2)+0.02,s=max⁡ ⁣(3, 1.6 σk)w_n = \exp\!\Bigl(-\frac{(n - \mu/k)^2}{2 s^2}\Bigr) + 0.02, \qquad s = \max\!\Bigl(3,\ 1.6\,\frac{\sigma}{\sqrt{k}}\Bigr)
Peso per numero del metodo di Gauss.

Modello a intervalli Δ

PROServer

Scopo: imitare la distanza tra numeri vicini nelle estrazioni reali.

Dalle ultime 120 estrazioni raccoglie il numero più piccolo di ogni estrazione e tutte le differenze tra numeri vicini ordinati. Una riga si costruisce come un inizio casuale più differenze casuali; una riga che supera NN viene scartata (fino a 300 tentativi, poi una riga uniforme):

x1∼{min⁡Yt},xj+1=xj+Δj,  Δj∼{ yt,i+1−yt,i }x_1 \sim \{\min Y_t\}, \qquad x_{j+1} = x_j + \Delta_j,\ \ \Delta_j \sim \{\,y_{t,i+1} - y_{t,i}\,\}
Bootstrap empirico dei primi numeri e degli intervalli; Y_t = (y_{t,1} < … < y_{t,k}) è l’estrazione t ordinata.
Dettagli avanzati

Bayes · Dirichlet

PROServer

Scopo: la frequenza con una contrazione (shrinkage) verso l’uniforme.

p^n=cn+αkD+αN,α=3\hat p_n = \frac{c_n + \alpha}{kD + \alpha N}, \qquad \alpha = 3
Media a posteriori di una distribuzione a priori di Dirichlet(α) simmetrica aggiornata con i conteggi osservati; il metodo campiona con pesi c_n + 3.
Dettagli avanzati

Analisi degli scostamenti

PROServer

Scopo: favorire i numeri che non escono da molto tempo.

wn=(gap⁡ngˉ)1.6+0.15,gˉ=max⁡ ⁣(1,Nk)w_n = \Bigl(\frac{\operatorname{gap}_n}{\bar g}\Bigr)^{1.6} + 0.15, \qquad \bar g = \max\!\Bigl(1, \frac{N}{k}\Bigr)
gap_n è l’indice dell’estrazione più recente che contiene n (0 = l’ultima estrazione); un numero mai uscito riceve D + ḡ.
Dettagli avanzati

Modello fisico della macchina della lotteria

PROServer

Scopo: lasciare che una simulazione 2D semplificata decida l’ordine di uscita.

NN dischi di raggio r=clamp⁡(0.12/N, 0.006, 0.018)r = \operatorname{clamp}(0.12/\sqrt N,\ 0.006,\ 0.018) si muovono in un quadrato unitario con raffiche d’aria casuali, una deriva costante verso l’uscita, smorzamento e urti elastici tra masse uguali; viene estratto il disco che esce dall’apertura al centro della parete inferiore.

vx←0.992 (vx+ξx),vy←0.992 (vy+0.0015+ξy),ξx∼U(±0.003), ξy∼U(±0.004)v_x \leftarrow 0.992\,(v_x + \xi_x), \quad v_y \leftarrow 0.992\,(v_y + 0.0015 + \xi_y), \qquad \xi_x \sim \mathcal U(\pm 0.003),\ \xi_y \sim \mathcal U(\pm 0.004)
Un passo del modello 2D (y punta verso il basso, verso l’uscita; al massimo 2400 passi). Un disco esce quando 0.43 < x < 0.57 sulla parete inferiore; se escono meno di k dischi, la riga viene completata in modo uniforme.
Dettagli avanzati
vi′=vi−((vi−vj)⋅n^) n^,vj′=vj+((vi−vj)⋅n^) n^\mathbf v_i' = \mathbf v_i - \bigl((\mathbf v_i - \mathbf v_j)\cdot\hat{\mathbf n}\bigr)\,\hat{\mathbf n}, \qquad \mathbf v_j' = \mathbf v_j + \bigl((\mathbf v_i - \mathbf v_j)\cdot\hat{\mathbf n}\bigr)\,\hat{\mathbf n}
Urto elastico tra masse uguali: le componenti normali delle velocità si scambiano.

Caos · mappa logistica

PROServer

Scopo: mostrare che un sistema del tutto deterministico può essere imprevedibile.

xt+1=r xt (1−xt),r=3.99,nt=1+⌊xtN⌋x_{t+1} = r\,x_t\,(1 - x_t), \quad r = 3.99, \qquad n_t = 1 + \lfloor x_t N \rfloor
Dopo 100 iterazioni di riscaldamento da un seme casuale, si raccolgono valori finché non si trovano k numeri distinti.
Dettagli avanzati

Collasso quantistico

PROServerSperimentale

Scopo: un peso controcorrente che favorisce i numeri estratti raramente.

wn=exp⁡ ⁣(−1.15 cncˉ)+0.02,cˉ=max⁡ ⁣(1,1N∑mcm)w_n = \exp\!\Bigl(-1.15\,\frac{c_n}{\bar c}\Bigr) + 0.02, \qquad \bar c = \max\!\Bigl(1, \frac{1}{N}\sum_m c_m\Bigr)
Peso a «barriera di frequenza»; i numeri casuali provengono da un buffer di valori crittografici a 16 bit.
Dettagli avanzati

Sistema di paradossi

PROServerSperimentale

Scopo: righe di ricerca costruite su proprietà controintuitive delle combinazioni casuali. Si usano a turno cinque archetipi:

  • Folla: preferisce i numeri sopra il 31 negli insiemi di almeno 42 numeri (i compleanni finiscono al 31), peso 2.6 sopra la soglia e 0.28 sotto, ×1.1 per i multipli di 10, poi impone almeno ⌈0.6k⌉\lceil 0.6k \rceil numeri sopra la soglia. Meno persone giocano questi numeri, quindi un’eventuale vincita sarebbe divisa con meno vincitori — la probabilità di vincere non cambia.
  • Scontro: metà della riga dai numeri più frequenti, il resto da quelli assenti da più tempo.
  • Grappolo: contiene una coppia di numeri consecutivi — una proprietà molto più comune nelle estrazioni casuali di quanto suggerisca l’intuizione:
P(∃ i:xi+1=xi+1)=1−(N−k+1k)(Nk)P\bigl(\exists\, i : x_{i+1} = x_i + 1\bigr) = 1 - \frac{\binom{N-k+1}{k}}{\binom{N}{k}}
Eurojackpot circa 35%, Lotto 7/34 circa 78%.
Dettagli avanzati
  • Specchio: abbina nn e N+1−nN + 1 - n.
  • Regressione: i numeri la cui frequenza recente (ultime 10 estrazioni) si discosta di più da quella di lungo periodo.

3. Il Tribunale matematico

Il Tribunale esamina una combinazione che hai già. Giurati e modelli possono proporre di sostituire un numero; gli avvocati della difesa si esprimono a favore o contro ogni proposta; il Giudice supremo decide sulla base di un campo strutturale delle estrazioni passate. I numeri cambiano solo quando lo decidi tu — e ogni passaggio viene registrato.

Selezione dei dodici giurati attorno all’emblema del Tribunale
La giuria: dodici personaggi, ciascuno con la propria strategia matematica

Lo svolgimento di un caso

  • Proposta. Un giurato o un modello esamina ogni possibile scambio singolo della riga e propone al massimo una sostituzione — oppure propone di mantenere la riga.
  • Difesa. Gli avvocati prendono posizione su ogni proposta: sostegno, obiezione o protesta. Un avvocato non propone né applica mai un numero.
  • Verdetto. Il Giudice supremo accetta la proposta, offre un’alternativa o mantiene il numero, basandosi sul campo strutturale delle estrazioni passate.
  • Decisione. Solo tu cambi i numeri: accetta una proposta, segui la difesa, accetta il Giudice, mantieni l’originale o torna indietro. Ogni decisione viene salvata con il suo autore.

Da un modello a una preferenza per numero

Ogni modello viene ridotto a una preferenza relativa rel⁡(n)\operatorname{rel}(n) per ogni numero, dove 1 significa «la media del modello», insieme a un errore standard. I modelli definiti da un vettore di pesi esplicito sono esatti; gli altri vengono misurati eseguendoli:

rel⁡(n)=wn1N∑mwm,se⁡(n)=0\operatorname{rel}(n) = \frac{w_n}{\frac{1}{N}\sum_m w_m}, \qquad \operatorname{se}(n) = 0
Preferenza relativa esatta di Numeri caldi, Markov, Gauss, Bayes, Analisi degli scostamenti, Collasso quantistico, Quantistico-astrale e Campo strutturale.
p^n=1R∑r=1R[ n∈Xr ],rel⁡(n)=p^nk/N,se⁡(n)=p^n(1−p^n)/Rk/N\hat p_n = \frac{1}{R}\sum_{r=1}^{R} \bigl[\,n \in X_r\,\bigr], \quad \operatorname{rel}(n) = \frac{\hat p_n}{k/N}, \quad \operatorname{se}(n) = \frac{\sqrt{\hat p_n (1 - \hat p_n)/R}}{k/N}
Preferenza Monte Carlo dei modelli procedurali: X₁ … X_R sono R = 600 righe del modello generate con seme (120 per il modello fisico). Il Consenso invece fa la media di tre esecuzioni (vedi la sua scheda).

I diciannove modelli

Il Tribunale usa i sedici metodi del generatore della sezione 2 più tre modelli propri. Quattro di essi non usano lo storico e sono registrati come privi di preferenza per i numeri (Copertura del segmento, Casuale puro, Modello fisico, Caos); in pratica Copertura del segmento favorisce i numeri dei segmenti corti quando kk non divide NN, e Caos favorisce gli estremi dell’intervallo, come mostrano le loro schede. Due modelli sono compositi (Consenso e Campo strutturale).

ModelloTipoBase
Numeri caldipesicn+0.05c_n + 0.05
Combinatocampionatoterzi caldo / medio / freddo
Copertura del segmentocampionato, senza storicoun numero per segmento
Casuale purocampionato, uniformerimescolamento di Fisher–Yates
Catene di Markovpesitransizioni lisciate dall’ultima estrazione
Gauss · TLCpesigaussiana attorno a μ/k\mu/k
Modello a intervalli Δcampionatobootstrap degli intervalli
Bayes · Dirichletpesicn+3c_n + 3
Analisi degli scostamentipesi(gap⁡/gˉ)1.6+0.15(\operatorname{gap}/\bar g)^{1.6} + 0.15
Modello fisicocampionato, senza storicosimulazione di dischi in 2D
Caoscampionato, senza storicomappa logistica, sbilanciata verso gli estremi
Collasso quantisticopesi, sperimentalee−1.15cn/cˉ+0.02e^{-1.15 c_n/\bar c} + 0.02
Sistema di paradossicampionato, sperimentalecinque archetipi
Quantistico-astralepesi, sperimentalesegno e illuminazione della Luna
Matrice della ruotacampionatodisegno di copertura delle coppie
Caldo mondialecampionatoprofilo a nucleo tra lotterie
Combinata mondialecampionatoprofilo a nucleo tra lotterie
Consenso dei modelliensemble14 modelli membri
Campo strutturalecompositosei componenti descrittive

Campo strutturale

Server

Il campo è la misura del Giudice supremo: un composito, pesato a mano, di sei componenti descrittive. Frequenze, connettività e copertura delle coppie vengono divise per il loro massimo su tutti i numeri; ritardo e tipicità sono già in [0,1][0,1].

F(n)=100 (0.25 An+0.15 Sn+0.15 Gn+0.20 Ln+0.15 Tn+0.10 (1−Pn))F(n) = 100\,\bigl(0.25\,A_n + 0.15\,S_n + 0.15\,G_n + 0.20\,L_n + 0.15\,T_n + 0.10\,(1 - P_n)\bigr)
Punteggio del campo strutturale del numero n, da 0 a 100.
Dettagli avanzati
dove
An, SnA_n,\ S_n
frequenza nell’intera finestra e nelle ultime max⁡(3,min⁡(10,D))\max(3, \min(10, D)) estrazioni
GnG_n
ritardo, min⁡(gap⁡n/gˉ,2)/2\min(\operatorname{gap}_n/\bar g, 2)/2
LnL_n
somma dei lift di coppia ℓnm=Onm(cn/D)(cm/D) D\ell_{nm} = \dfrac{O_{nm}}{(c_n/D)(c_m/D)\,D} su tutti i partner mm, con OnmO_{nm} il numero di estrazioni che contengono entrambi
TnT_n
tipicità strutturale media delle estrazioni che contengono nn: 1−0.6⋅2 ∣e−12∣−0.4min⁡(∣zS∣,2)/21 - 0.6 \cdot 2\,|e - \tfrac12| - 0.4 \min(|z_S|, 2)/2, con ee la quota di numeri pari e zSz_S la somma standardizzata
PnP_n
copertura delle coppie, il numero totale di osservazioni di coppie che coinvolgono nn

Consenso dei modelli

Server

Quattordici modelli membri di dodici famiglie producono ciascuno 10 righe; un numero ottiene punti in base a quanti modelli e quante famiglie lo hanno scelto, a quanto spesso e a quanto in modo uniforme:

C(n)=100(0.40 Mn14+0.25 Fn12+0.20 fnmax⁡mfm+0.15 min⁡(stab⁡n,1))C(n) = 100\Bigl(0.40\,\frac{M_n}{14} + 0.25\,\frac{F_n}{12} + 0.20\,\frac{f_n}{\max_m f_m} + 0.15\,\min(\operatorname{stab}_n, 1)\Bigr)
Punteggio di consenso; il Tribunale fa la media di tre esecuzioni indipendenti e usa la loro dispersione come errore standard.
Dettagli avanzati
dove
Mn, FnM_n,\ F_n
numero di modelli membri e di famiglie di modelli le cui righe contengono nn
fnf_n
numero totale di righe che contengono nn
stab⁡n\operatorname{stab}_n
numero medio di comparse per modello membro che ha scelto nn, diviso per 10

Modello quantistico-astrale

ServerSperimentale

Un modello di intrattenimento che pesa i numeri in base alla Luna nella data dell’estrazione di riferimento, più una piccola quota del campo strutturale:

wn=1+[ n mod 12=σ ](0.10+0.20 I100)+0.15 F(n)100w_n = 1 + \bigl[\,n \bmod 12 = \sigma\,\bigr]\Bigl(0.10 + 0.20\,\frac{I}{100}\Bigr) + 0.15\,\frac{F(n)}{100}
Peso quantistico-astrale; σ ∈ {0,…,11} è il segno zodiacale della longitudine eclittica della Luna nella data dell’estrazione di riferimento, I l’illuminazione in percentuale.
Dettagli avanzati
I=100⋅1−cos⁡(2πA/P)2,P=29.530588853 dI = 100 \cdot \frac{1 - \cos(2\pi A / P)}{2}, \qquad P = 29.530588853\ \mathrm{d}
Illuminazione della Luna in percentuale a partire dall’età A nel mese sinodico di durata P (giorni).

Come un modello esamina una riga

Per ogni possibile scambio a→ba \to b — il numero aa esce dalla riga, il numero bb entra — il modello calcola il guadagno e il relativo rumore. Un guadagno conta solo se è statisticamente netto; una perdita conta sempre.

δ=rel⁡(b)−rel⁡(a),z=δSE⁡δ,effect⁡={δSE⁡δ=0 ∨ z≥z∗min⁡(δ,0)SE⁡δ>0 ∧ z<z∗\delta = \operatorname{rel}(b) - \operatorname{rel}(a), \qquad z = \frac{\delta}{\operatorname{SE}_\delta}, \qquad \operatorname{effect} = \begin{cases} \delta & \operatorname{SE}_\delta = 0 \ \vee\ z \ge z^{*} \\ \min(\delta, 0) & \operatorname{SE}_\delta > 0 \ \wedge\ z < z^{*} \end{cases}
Esame da parte di un modello. I modelli a pesi esatti hanno SE = 0, quindi i loro guadagni contano sempre; lo scambio è ammesso quando effect ≥ 0.10.
SE⁡δ=1k/N2 pˉ (1−pˉ)R,pˉ=max⁡ ⁣(rel⁡(a)+rel⁡(b)2 kN, 12R)\operatorname{SE}_\delta = \frac{1}{k/N}\sqrt{\frac{2\,\bar p\,(1 - \bar p)}{R}}, \qquad \bar p = \max\!\Bigl(\frac{\operatorname{rel}(a) + \operatorname{rel}(b)}{2}\,\frac{k}{N},\ \frac{1}{2R}\Bigr)
Errore standard di uno scambio per un modello campionato: la proporzione combinata sotto l’ipotesi che a e b siano ugualmente probabili. Per il Consenso SE_δ = √(se(a)² + se(b)²) sulle sue tre esecuzioni; i modelli a pesi esatti hanno SE_δ = 0.
z∗=Φ−1 ⁣(1−1−Φ(1.5)k(N−k))z^{*} = \Phi^{-1}\!\Bigl(1 - \frac{1 - \Phi(1.5)}{k(N-k)}\Bigr)
Soglia con correzione di Bonferroni su tutti i k(N−k) scambi possibili; Eurojackpot: 3.43.
zz di baseEurojackpot, 225 scambiLotto 7/34, 189 scambi
0.52.9952.942
1.03.1933.142
1.53.4343.387
2.03.7163.672
3.04.3784.340

I dodici giurati

Ogni giurato ha una strategia, le proprie fonti (modelli o segnali dello storico), un effetto minimo e un livello di significatività. Il verdetto di un giurato su una riga è «mantieni» oppure esattamente una sostituzione; i giurati non sollevano obiezioni né proteste — quello è il ruolo della difesa.

score⁡(n)=∑mωmrel⁡m(n),Δ=λ (score⁡(b)−score⁡(a))\operatorname{score}(n) = \sum_m \omega_m \operatorname{rel}_m(n), \qquad \Delta = \lambda\,\bigl(\operatorname{score}(b) - \operatorname{score}(a)\bigr)
Guadagno numerico di uno scambio dalle fonti pesate del giurato; λ = D/(D+40) per Eva Lind, altrimenti 1.
E={ΔSE⁡Δ=0 ∨ z≥z∗(zmin⁡)min⁡(Δ,0)SE⁡Δ>0 ∧ z<z∗(zmin⁡) + ωS⋅0.25∑cΔℓc + ωC ΔΓ,E≥Emin⁡ ⇒ a→bE = \begin{cases} \Delta & \operatorname{SE}_\Delta = 0 \ \vee\ z \ge z^{*}(z_{\min}) \\ \min(\Delta, 0) & \operatorname{SE}_\Delta > 0 \ \wedge\ z < z^{*}(z_{\min}) \end{cases} \ +\ \omega_S \cdot 0.25 \sum_c \Delta \ell_c \ +\ \omega_C\,\Delta \Gamma, \qquad E \ge E_{\min} \ \Rightarrow\ a \to b
Valore effettivo di uno scambio. Vince il miglior scambio ammesso (a parità: il nuovo numero più piccolo, poi la posizione precedente); se non ce n’è nessuno, il giurato mantiene la riga.
dove
ωm\omega_m
peso della fonte mm nella definizione del giurato
zz
Δ/SE⁡Δ\Delta / \operatorname{SE}_\Delta, con SE⁡Δ=λ∑m(ωmSE⁡m)2\operatorname{SE}_\Delta = \lambda\sqrt{\sum_m (\omega_m \operatorname{SE}_m)^2}
z∗(zmin⁡)z^{*}(z_{\min})
il livello di significatività del giurato corretto per k(N−k)k(N-k) scambi, come nella tabella sopra
Δℓc\Delta \ell_c
variazione della log-verosimiglianza della componente strutturale cc (in nat)
ΔΓ\Delta \Gamma
variazione del punteggio di copertura della schedina

Segnali calcolati dallo storico

Tre giurati usano segnali letti direttamente dalle estrazioni anziché da un generatore. Con p=k/Np = k/N:

rel⁡(x)=cxDp,SE⁡ab=2(N−k)(N−1) p D\operatorname{rel}(x) = \frac{c_x}{Dp}, \qquad \operatorname{SE}_{ab} = \sqrt{\frac{2(N-k)}{(N-1)\,p\,D}}
Test di frequenza. In estrazioni eque Var(c_a − c_b) = 2Dp(1−p)·N/(N−1), perché i numeri di una stessa estrazione sono correlati negativamente; dividendo per il conteggio atteso Dp si ottiene SE_ab.
rel⁡(x)=max⁡ ⁣(0, fˉx−sd⁡(fx)K),fx,w=cx,ws p,K=min⁡ ⁣(4,⌊D/5⌋)\operatorname{rel}(x) = \max\!\Bigl(0,\ \bar f_x - \frac{\operatorname{sd}(f_x)}{\sqrt K}\Bigr), \qquad f_{x,w} = \frac{c_{x,w}}{s\,p}, \qquad K = \min\!\bigl(4, \lfloor D/5 \rfloor\bigr)
Replicazione. Frequenza in K finestre disgiunte di s = ⌊D/K⌋ estrazioni; il vantaggio deve resistere a una penalità di un errore standard.
raw⁡(x)=1+∣cx−DpDp(1−p)∣,rel⁡=raw⁡raw⁡‾,SE⁡=1−2/πraw⁡‾\operatorname{raw}(x) = 1 + \left|\frac{c_x - Dp}{\sqrt{Dp(1-p)}}\right|, \qquad \operatorname{rel} = \frac{\operatorname{raw}}{\overline{\operatorname{raw}}}, \qquad \operatorname{SE} = \frac{\sqrt{1 - 2/\pi}}{\overline{\operatorname{raw}}}
Anomalia. Deviazione assoluta dal conteggio atteso; √(1−2/π) è la deviazione standard di |Z| per una normale standard Z.

Struttura e copertura

Le componenti strutturali confrontano una riga con il profilo delle estrazioni passate: parità, bassi/alti, grappoli e intervalli sono categoriali con lisciamento di Laplace; somma e ampiezza sono gaussiane.

ℓ1=ln⁡nj+1D+J,ℓ2=−12(x−μσ)2,ℓ3=1k−1∑iln⁡nΔi+1nΔ+N\ell_1 = \ln\frac{n_j + 1}{D + J}, \qquad \ell_2 = -\frac12\Bigl(\frac{x - \mu}{\sigma}\Bigr)^2, \qquad \ell_3 = \frac{1}{k-1}\sum_{i} \ln\frac{n_{\Delta_i} + 1}{n_\Delta + N}
Log-verosimiglianze: ℓ₁ categoriale (parità, bassi/alti, grappoli; J categorie, n_j estrazioni passate nella categoria j), ℓ₂ gaussiana (somma, ampiezza), ℓ₃ intervalli tra numeri vicini.
Γ(R)=−∣(R2)∩P′∣(k2)+0.5 ∣{⌊(n−1)/10⌋:n∈R}∣k\Gamma(R) = -\frac{\bigl|\tbinom{R}{2} \cap \mathcal P'\bigr|}{\binom{k}{2}} + 0.5\,\frac{\bigl|\{\lfloor (n-1)/10 \rfloor : n \in R\}\bigr|}{k}
Copertura della schedina: 𝒫′ è l’insieme delle coppie delle altre righe della schedina; il secondo termine conta le decine del campo numerico toccate dalla riga.

Eva Lind · Tiratrice scelta

Strategia: frequenza di lungo periodo con contrazione bayesiana. Fonti: Bayes 0.8 + Numeri caldi 0.2; λ=D/(D+40)\lambda = D/(D+40); Emin⁡=0.12E_{\min} = 0.12; nessun test di significatività (pesi esatti).

Regola: il miglior scambio è sempre il numero meno frequente della riga → il numero più frequente fuori dalla riga; è ammesso quando Δ≥0.12\Delta \ge 0.12.

Dettagli avanzati

Maya Chen · Architetta della forma

Strategia: la forma della riga. Fonti: Gauss 0.5 + Combinato 0.5; peso della struttura 2 su parità, bassi/alti, somma, ampiezza e grappoli; zmin⁡=1z_{\min} = 1 (Eurojackpot z∗=3.19z^{*} = 3.19); Emin⁡=0.12E_{\min} = 0.12.

Regola: il termine di struttura conta sempre, quindi uno scambio che avvicina la riga alla forma tipica delle estrazioni passate può essere ammesso per la sola struttura.

Dettagli avanzati

Clara Moretti · Cacciatrice di anomalie

Strategia: anomalie in entrambe le direzioni. Fonte: il segnale di anomalia; struttura degli intervalli 0.4, conteggiata solo per uno scambio significativo; zmin⁡=1.5z_{\min} = 1.5 (Eurojackpot z∗=3.43z^{*} = 3.43); Emin⁡=0.10E_{\min} = 0.10.

z=∣zb∣−∣za∣2⋅0.60281z = \frac{|z_b| - |z_a|}{\sqrt 2 \cdot 0.60281}
Per il segnale di anomalia la media si annulla: conta solo la differenza delle deviazioni assolute.
Dettagli avanzati

Leila Haddad · Osservatrice globale

Strategia: il profilo mondiale. Fonti: Caldo mondiale 0.5 + Combinata mondiale 0.5; zmin⁡=1.5z_{\min} = 1.5; Emin⁡=0.12E_{\min} = 0.12; richiede i dati delle altre lotterie.

Dettagli avanzati

Freya Olsen · Ingegnera della copertura

Strategia: copertura della schedina senza preferenza per i numeri. Nessuna fonte; peso della copertura 2; Emin⁡=0.20E_{\min} = 0.20.

Dettagli avanzati

Nora Blake · Chiaroveggente

Sperimentale

Strategia: il sistema astrologico di intrattenimento per la data dell’estrazione. Fonte: Quantistico-astrale; Emin⁡=0.05E_{\min} = 0.05; richiede una data di estrazione di riferimento.

Dettagli avanzati

Derek Cole · Maestro delle catene

Strategia: le transizioni dall’ultima estrazione. Fonte: Catene di Markov; struttura degli intervalli 0.4, sempre conteggiata; Emin⁡=0.12E_{\min} = 0.12. Poiché i pesi di Markov sommano a uno, rel⁡(n)=N wn\operatorname{rel}(n) = N\,w_n.

Dettagli avanzati

Marcus Reed · Segnale d’allarme

Strategia: l’assenza più lunga. Fonti: Analisi degli scostamenti 0.8 + Numeri caldi 0.2; Emin⁡=0.06E_{\min} = 0.06; nessun test di significatività.

Dettagli avanzati

Victor Moreau · Fisico dell’urna

Strategia: una distorsione fisica dell’urna si ripeterebbe in ogni periodo. Fonte: replicazione; zmin⁡=2z_{\min} = 2 (Eurojackpot z∗=3.72z^{*} = 3.72); Emin⁡=0.15E_{\min} = 0.15.

SE⁡ab=max⁡ ⁣(sd⁡w(fb,w−fa,w)K, 2(N−k)(N−1) p Ks)\operatorname{SE}_{ab} = \max\!\Bigl(\frac{\operatorname{sd}_w(f_{b,w} - f_{a,w})}{\sqrt K},\ \sqrt{\frac{2(N-k)}{(N-1)\,p\,K s}}\Bigr)
La soglia minima di rumore impedisce che quattro finestre fortunate sembrino significative.
Dettagli avanzati

Noah Bennett · Scettico vigile

Strategia: rifiutare il puro caso al livello più severo. Fonte: test di frequenza; zmin⁡=3z_{\min} = 3 (Eurojackpot z∗=4.38z^{*} = 4.38); Emin⁡=0.15E_{\min} = 0.15.

z=cb−ca2(N−k) Dp/(N−1)z = \frac{c_b - c_a}{\sqrt{2(N-k)\,Dp/(N-1)}}
Forma chiusa della statistica test di Noah.
Dettagli avanzati

Elias Berg · Prestigiatore controcorrente

Sperimentale

Strategia: contro la folla — numeri più rari. Fonti: Collasso quantistico 0.85 + Sistema di paradossi 0.15; zmin⁡=0.5z_{\min} = 0.5 (Eurojackpot z∗=3.00z^{*} = 3.00); Emin⁡=0.08E_{\min} = 0.08.

Dettagli avanzati

Leo Fischer · Direttore del consenso

Strategia: la stabilità del consenso dei modelli. Fonte: Consenso (tre esecuzioni indipendenti); zmin⁡=2z_{\min} = 2 (Eurojackpot z∗=3.72z^{*} = 3.72); Emin⁡=0.10E_{\min} = 0.10.

Dettagli avanzati

Generazione da parte della giuria

Quando la giuria genera una riga invece di esaminarne una, ogni giurato estrae 40 righe candidate con i propri pesi max⁡(0.02,score⁡(n))\max(0.02, \operatorname{score}(n)) e tiene la migliore secondo

rowScore⁡(R)=1k∑n∈Rscore⁡(n)+ωS⋅0.25∑cℓc(R)+ωCcov⁡(R)\operatorname{rowScore}(R) = \frac{1}{k}\sum_{n \in R} \operatorname{score}(n) + \omega_S \cdot 0.25 \sum_c \ell_c(R) + \omega_C \operatorname{cov}(R)
La scheda di voto di un giurato. La riga della giuria prende i k numeri con più voti; a parità prevale il sostegno normalizzato più alto, poi il numero più piccolo.

I sei avvocati della difesa

Gli avvocati valutano uno scambio proposto nc→npn_c \to n_p — il numero attuale e quello proposto, usciti xcx_c e xpx_p volte in DD estrazioni — con i propri segnali. I valori positivi sono a favore della proposta:

SegnaleDefinizione
frequenzaz=(xp−xc)/D2pˉ(1−pˉ)/D, pˉ=xc+xp2Dz = \dfrac{(x_p - x_c)/D}{\sqrt{2\bar p(1-\bar p)/D}},\ \bar p = \dfrac{x_c + x_p}{2D}
a posterioriz=Φ−1(P(θp>θc)), θ∼Beta⁡(1+x,1+D−x)z = \Phi^{-1}\bigl(P(\theta_p > \theta_c)\bigr),\ \theta \sim \operatorname{Beta}(1 + x, 1 + D - x)
Bayes, campo, astroz=rel⁡(p)−rel⁡(c)2 sd⁡(rel⁡)z = \dfrac{\operatorname{rel}(p) - \operatorname{rel}(c)}{\sqrt 2\,\operatorname{sd}(\operatorname{rel})}
strutturaz=sign⁡(Δℓ)2∣Δℓ∣z = \operatorname{sign}(\Delta\ell)\sqrt{2|\Delta\ell|}
stabilitàz=(W−L)/6z = (W - L)/\sqrt{6} su sei periodi consecutivi

La probabilità a posteriori viene calcolata in modo esatto, senza approssimazione normale:

P(θB>θA)=∑i=0aB−1B(aA+i, bA+bB)(bB+i) B(1+i, bB) B(aA,bA)P(\theta_B > \theta_A) = \sum_{i=0}^{a_B - 1} \frac{\mathrm B(a_A + i,\ b_A + b_B)}{(b_B + i)\,\mathrm B(1 + i,\ b_B)\,\mathrm B(a_A, b_A)}
Forma chiusa per distribuzioni a posteriori Beta con parametri interi; la funzione Beta usa un log-Gamma di Lanczos.

Regola di decisione per la sostituzione, con zˉ\bar z la media pesata dei segnali dell’avvocato:

ψ={Szˉ≥zsup ∧ χ=1Pzˉ≤−zprot ∧ χ=1Ozˉ∈(−zprot, zsup) ∨ χ=0\psi = \begin{cases} \mathsf S & \bar z \ge z_{\mathrm{sup}} \ \wedge\ \chi = 1 \\ \mathsf P & \bar z \le -z_{\mathrm{prot}} \ \wedge\ \chi = 1 \\ \mathsf O & \bar z \in (-z_{\mathrm{prot}},\ z_{\mathrm{sup}}) \ \vee\ \chi = 0 \end{cases}
Posizione ψ di un avvocato: S sostegno, P protesta, O obiezione. χ = 1 quando i controlli dell’avvocato sono superati: per il sostegno — unanimità dei segnali, concordanza di un secondo segnale e veto strutturale; per la protesta — solo unanimità e concordanza.

Thomas Keller · Avvocato d’ufficio

FREE

Segnali: frequenza + Bayes + campo strutturale, che devono concordare tutti; sostegno a 1, protesta a −1; richiede 20 estrazioni, e da 5 a 19 estrazioni funziona in modalità limitata con Bayes + a posteriori.

Dettagli avanzati

Elena Rossi · Avvocata statistica

Segnale: frequenza, corretta per i confronti multipli; veto strutturale a −1; sostegno e protesta a 1.96; almeno 30 estrazioni.

Dettagli avanzati

Amelia Grant · Contro la falsa sicurezza

Segnale: probabilità a posteriori di superiorità; sostegno a P≥0.90P \ge 0.90 (z=1.2816z = 1.2816), protesta a P≤0.20P \le 0.20 (z=−0.8416z = -0.8416); almeno 10 estrazioni.

Dettagli avanzati
DDxcx_cxpx_pP(θp>θc)P(\theta_p > \theta_c)Posizione
30170.987sostegno
30430.354obiezione
30520.128protesta
1006160.987sostegno

Julien Morel · Avvocato della struttura

Segnale: variazione totale delle sei log-verosimiglianze strutturali; sostegno a 0.6325 (Δℓ≥0.2\Delta\ell \ge 0.2 nat), protesta a 0.8944 (Δℓ≤−0.4\Delta\ell \le -0.4); almeno 20 estrazioni.

Dettagli avanzati

Ingrid Falk · Stabilità nel tempo

Segnale: vittorie meno sconfitte del numero proposto su sei periodi consecutivi; deve concordare con la frequenza complessiva; soglia 1.5, quindi almeno W−L≥4W - L \ge 4; almeno 30 estrazioni.

Dettagli avanzati

Adrian Vega · Astrologia d’intrattenimento

Sperimentale

Segnale: i pesi quantistico-astrali per la data dell’estrazione di riferimento; sostegno e protesta a 0.5.

Dettagli avanzati

Il Giudice supremo

Il Giudice decide su un numero contestato usando il campo strutturale FF, normalizzato al suo massimo: ϕ(n)=F(n)/max⁡F\phi(n) = F(n)/\max F. Un numero è debole se si trova nel 35% più basso del campo.

θ=ϕ(⌊0.35N⌋+1),m=0.08⋅[ nP>12 ncounsel ]\theta = \phi_{(\lfloor 0.35 N \rfloor + 1)}, \qquad m = 0.08 \cdot \bigl[\, n_{\mathsf P} > \tfrac12\, n_{\mathrm{counsel}} \,\bigr]
θ è una statistica d’ordine dei valori del campo ordinati in senso crescente; il numero attuale è debole quando φ(n_c) < θ. Il margine m si applica quando protesta più della metà degli avvocati che hanno preso posizione.
V={ncϕ(nc)≥θnpϕ(np)>ϕ(nc)+m ∧ ϕ(np)≥θa∗ϕ(a∗)>ϕ(nc)+mncϕ(a∗)≤ϕ(nc)+mV = \begin{cases} n_c & \phi(n_c) \ge \theta \\ n_p & \phi(n_p) > \phi(n_c) + m \ \wedge\ \phi(n_p) \ge \theta \\ a^{*} & \phi(a^{*}) > \phi(n_c) + m \\ n_c & \phi(a^{*}) \le \phi(n_c) + m \end{cases}
Il verdetto V è il numero che il Giudice mantiene nella posizione; si applica la prima riga che corrisponde. a* è il numero con il campo più alto fuori dalla riga e diverso dalla proposta.

Analisi e raccomandazione del Giudice supremo

Per righe intere il Giudice riporta un punteggio della riga e suggerisce di sostituire fino a ⌊k/2⌋\lfloor k/2 \rfloor dei numeri più deboli sotto θ\theta con i numeri più forti fuori dalla riga. Sei tu a scegliere quali scambi applicare. Senza righe, il Giudice formula una raccomandazione: bozze casuali uniformi vengono migliorate allo stesso modo e restituite con la loro provenienza.

RowScore⁡(R)=100k∑n∈Rϕ(n)\operatorname{RowScore}(R) = \frac{100}{k}\sum_{n \in R} \phi(n)
Punteggio della riga nell’analisi del Giudice supremo.

Decisioni, lotti e cronologia

Il Tribunale alla prova dei fatti

Ogni affermazione viene scritta nel registro delle previsioni prima della sua estrazione e poi valutata rispetto al caso. Per affermazioni su nin_i numeri con probabilità di base p=k/Np = k/N:

E[H]=∑ini p,Var⁡[H]=∑ini p(1−p),lift⁡=HE[H],z=H−E[H]Var⁡[H]\mathbb E[H] = \sum_i n_i\,p, \qquad \operatorname{Var}[H] = \sum_i n_i\,p(1-p), \qquad \operatorname{lift} = \frac{H}{\mathbb E[H]}, \qquad z = \frac{H - \mathbb E[H]}{\sqrt{\operatorname{Var}[H]}}
H è il numero di centri. Una fonte entra in classifica solo con almeno 30 numeri su almeno 3 estrazioni.

4. L’urna 3D: motore fisico

L’estrazione 3D è una simulazione di corpi rigidi in tempo reale: le palline sono corpi con massa, inerzia, attrito e restituzione, spinti dalla gravità, da un rotore miscelatore e da un campo d’aria, e la pallina estratta è quella che cade fisicamente attraverso l’uscita. È un modello fisico credibile di un’urna, non la registrazione di una macchina reale.

Urna del SuperEnalotto con 90 palline in miscelazione mentre viene selezionata una pallina
90 palline del SuperEnalotto nella camera simulata

Motore e coordinate

La fisica gira in Rapier 3D (versione 0.14, WebAssembly); il rendering usa Three.js (r169). Il sistema di riferimento è destrorso, con +y+y verso l’alto e l’origine al centro della camera; l’asse del rotore è l’asse orizzontale xx e +z+z punta verso l’osservatore.

m=43πr3ρ=0.03351,I=25mr2=5.36×10−4,m ∣g∣=1.005m = \tfrac43 \pi r^3 \rho = 0.03351, \qquad I = \tfrac25 m r^2 = 5.36 \times 10^{-4}, \qquad m\,|g| = 1.005
Massa, momento d’inerzia e peso di una pallina in unità del modello.

Scala e similitudine di Froude

Un’unità del modello vale 10 cm: una pallina da 40 mm misura 0.4 unità e la camera è una sfera da 600 mm. La gravità è ridotta a 3030 unità/s² invece dei 98.198.1 reali. Un moto dominato da gravità e inerzia appare identico quando il numero di Froude è lo stesso, quindi la simulazione corrisponde a una macchina in scala reale con un orologio più lento:

Fr=v2gL,γ=g′g=3098.1=0.306,t′t=1γ≈1.81\mathrm{Fr} = \frac{v^2}{g L}, \qquad \gamma = \frac{g'}{g} = \frac{30}{98.1} = 0.306, \qquad \frac{t'}{t} = \frac{1}{\sqrt\gamma} \approx 1.81
Le grandezze con l’apice appartengono al modello. A parità di numero di Froude e di lunghezza, le velocità scalano con √γ e i tempi con 1/√γ.

Forze

Ogni pallina obbedisce alla seconda legge di Newton con le forze che il codice applica davvero:

m d2xdt2=m g+FD+FW+FS+FC,I dωdt=τW+τCm\,\frac{d^2 \mathbf x}{dt^2} = m\,\mathbf g + \mathbf F_D + \mathbf F_W + \mathbf F_S + \mathbf F_C, \qquad I\,\frac{d\boldsymbol\omega}{dt} = \boldsymbol\tau_W + \boldsymbol\tau_C
F_D resistenza dell’aria, F_W parete analitica, F_S aspirazione all’uscita, F_C contatti con palline, rotore e chiusura (calcolati da Rapier); τ i momenti corrispondenti. Rapier applica anche un piccolo smorzamento lineare e angolare, e le velocità sono limitate a 22 unità/s e 105 rad/s come rete di sicurezza.

Campo d’aria e resistenza quadratica

Il soffiatore è modellato come un campo prescritto di velocità dell’aria u(x,t)\mathbf u(\mathbf x, t); ogni pallina subisce una resistenza quadratica verso la velocità locale dell’aria. La costante di resistenza è fissata da una velocità limite di 48 unità/s:

FD=m kd ∥u−v∥ (u−v),kd=∣g∣vt2=30482=0.01302\mathbf F_D = m\,k_d\,\lVert \mathbf u - \mathbf v \rVert\,(\mathbf u - \mathbf v), \qquad k_d = \frac{|g|}{v_t^2} = \frac{30}{48^2} = 0.01302
In aria ferma una pallina in caduta raggiunge v_t = √(g/k_d) = 48 unità/s, circa 8.7 m/s in scala reale.
uy↑=J(1−ρh2R2) clamp⁡ ⁣(6−y6+R, 0, 1),J=102u_y^{\uparrow} = J \Bigl(1 - \frac{\rho_h^2}{R^2}\Bigr)\, \operatorname{clamp}\!\Bigl(\frac{6 - y}{6 + R},\, 0,\, 1\Bigr), \qquad J = 102
Corrente ascensionale centrale: massima sull’asse vicino al fondo (102 a y = −3, 68 al centro, 34 in cima); ρ_h è la distanza orizzontale dall’asse.

Parete sferica analitica

Oltre alla mesh, il vetro viene imposto come sfera esatta a livello di velocità, il che impedisce il tunnelling e dà un rotolamento corretto. Con normale uscente n^\hat{\mathbf n}, velocità uscente vnv_n e penetrazione ξ=∣x∣−(R−r−0.06)\xi = |\mathbf x| - (R - r - 0.06):

v∗={−min⁡ ⁣(2, 0.25 ξ/Δt)ξ>0−ξ/Δtξ≤0,e={12(eb+0.28)vn>0.60vn≤0.6v^{*} = \begin{cases} -\min\!\bigl(2,\ 0.25\,\xi/\Delta t\bigr) & \xi > 0 \\ -\xi/\Delta t & \xi \le 0 \end{cases}, \qquad e = \begin{cases} \tfrac12 (e_b + 0.28) & v_n > 0.6 \\ 0 & v_n \le 0.6 \end{cases}
Tolleranza di contatto speculativo v* e restituzione e (e_b è la restituzione della pallina; zero sotto la velocità di tolleranza, così una pallina ferma non può vibrare).
Δvn=vn−min⁡(v∗, −e vn),Fn=−m ΔvnΔt n^\Delta v_n = v_n - \min(v^{*},\, -e\,v_n), \qquad \mathbf F_n = -\frac{m\,\Delta v_n}{\Delta t}\,\hat{\mathbf n}
Risposta normale, applicata solo se v_n > v*.
jt=min⁡ ⁣(μ m Δvn, 27m ∥s∥),Ft=−jtΔt s^,τ=(−rn^)×Ftj_t = \min\!\Bigl(\mu\, m\,\Delta v_n,\ \tfrac27 m\,\lVert \mathbf s \rVert\Bigr), \qquad \mathbf F_t = -\frac{j_t}{\Delta t}\,\hat{\mathbf s}, \qquad \boldsymbol\tau = (-r\hat{\mathbf n}) \times \mathbf F_t
Attrito di Coulomb nel punto di contatto; s è la velocità di slittamento tangenziale v + ω × (−r n̂), μ = ½(μ_b + 0.35) con μ_b l’attrito della pallina.

Contatti tra corpi

I contatti pallina–pallina, pallina–rotore e pallina–chiusura sono risolti dal risolutore di vincoli di Rapier con rilevamento continuo delle collisioni sulle palline. I coefficienti si combinano facendo la media:

CoppiaAttritoRestituzione
pallina – pallina0.160.55
pallina – vetro0.2550.415
pallina – pala0.190.525
pallina – braccio o albero0.170.35
pallina – chiusura0.230.375

Lo smorzamento è basso durante l’estrazione (lineare 0.05, angolare 0.035) e forte mentre le palline si assestano (2.4). Non esiste un parametro di resistenza al rotolamento: il rotolamento rallenta solo per lo smorzamento e l’attrito di contatto.

Passo temporale e integrazione

La simulazione avanza a passi fissi indipendenti dalla frequenza di aggiornamento dello schermo; ogni fotogramma esegue tanti passi fissi quanti ne richiede il tempo reale, al massimo cinque:

A←min⁡(A+Δtf, 5 Δt),n=⌊A/Δt⌋,A←A−n Δt,Δt=160 sA \leftarrow \min\bigl(A + \Delta t_f,\ 5\,\Delta t\bigr), \qquad n = \bigl\lfloor A / \Delta t \bigr\rfloor, \qquad A \leftarrow A - n\,\Delta t, \qquad \Delta t = \tfrac{1}{60}\ \mathrm{s}
Accumulatore a passo fisso: un fotogramma di durata Δt_f esegue n passi fissi (al massimo 5); ogni passo è diviso in 2 sottopassi di Rapier da h = 1/120 s con 6 iterazioni del risolutore.

All’interno di un sottopasso il motore aggiorna le velocità dalle forze e poi le posizioni dalle nuove velocità — lo schema di Eulero semi-implicito (simplettico) usato dai motori basati su impulsi, mostrato qui come principio generale:

vt+h=vt+h M−1(F+Fc),xt+h=xt+h vt+h\mathbf v_{t+h} = \mathbf v_t + h\,M^{-1}\bigl(\mathbf F + \mathbf F_{\mathrm c}\bigr), \qquad \mathbf x_{t+h} = \mathbf x_t + h\,\mathbf v_{t+h}
Forma generale: F sono le forze applicate, F_c le forze di vincolo (di contatto) trovate dal risolutore iterativo di Rapier.

Le forze proprie dell’applicazione (aria, parete, aspirazione) vengono calcolate una volta per passo dallo stato all’inizio del passo e mantenute per entrambi i sottopassi. Il rotore è cinematico e il suo obiettivo viene riaggiornato a ogni sottopasso:

ω←ω+(ω∗−ω)(1−e−4h),θ←θ+ω h\omega \leftarrow \omega + (\omega^{*} - \omega)\bigl(1 - e^{-4h}\bigr), \qquad \theta \leftarrow \theta + \omega\,h
Avvicinamento esponenziale della velocità del rotore al suo obiettivo; la miscelazione passa ciclicamente per 2.4, 2.9, 2.1 e 2.7 rad/s.

Come viene selezionata una pallina

Dopo 1 s di avvio e 2.6 s di miscelazione la chiusura si apre. Finché l’uscita è aperta, la corrente ascensionale sopra di essa si spegne gradualmente, il getto centrale scende al 62% e un’aspirazione di circa 2.2 pesi di pallina (circa 4 dopo 5 s) attira le palline vicine all’uscita:

FS=β (−1.2 x, −2.2, −1.2 z),β={1tc≤5 s1.8tc>5 s\mathbf F_S = \beta\,(-1.2\,x,\ -2.2,\ -1.2\,z), \qquad \beta = \begin{cases} 1 & t_c \le 5\ \mathrm{s} \\ 1.8 & t_c > 5\ \mathrm{s} \end{cases}
Aspirazione vicino all’uscita (y < −1.6 e distanza orizzontale < 0.6); t_c è il tempo trascorso dall’apertura della chiusura.

La pallina estratta è la prima il cui centro scende sotto y=−3.18y = -3.18; se nello stesso passo succede a più palline, la più bassa. In quel momento non si consulta alcun numero casuale. I numeri casuali di un’estrazione sono la disposizione iniziale delle palline e una piccola velocità e rotazione iniziali, presi dalla Web Crypto API (di nuovo quando viene caricato un insieme bonus separato). La pallina percorre poi il tragitto di uscita fino alla rastrelliera; quel moto è scriptato, non simulato.

Rendering, simulazione e casualità

5. Statistica delle estrazioni reali

Le schermate di analisi descrivono le estrazioni storiche. Rispondono a «che cosa è successo» e «è compatibile con una lotteria equa», mai a «che cosa succederà».

Frequenze, numeri caldi e freddi

Per ogni numero l’app conta le uscite nella finestra e le mostra accanto al conteggio atteso di una lotteria equa; la deviazione standard è riportata qui come riferimento:

E[cn]=kDN,sd⁡(cn)=D kN(1−kN)\mathbb E[c_n] = \frac{kD}{N}, \qquad \operatorname{sd}(c_n) = \sqrt{D\,\frac{k}{N}\Bigl(1 - \frac{k}{N}\Bigr)}
Conteggio atteso (mostrato nell’app) e sua deviazione standard (riferimento) per un numero su D estrazioni indipendenti.
Mappa di calore dell’Eurojackpot dei 50 numeri principali colorata dal blu al rosso
Mappa di calore: il colore mostra solo la frequenza passata

Test di equità: χ² con correzione

Il test χ² verifica se i conteggi sono compatibili con probabilità uguali. In un’estrazione i kk numeri sono presi senza reinserimento, il che rende i conteggi correlati negativamente; la statistica multinomiale semplice risulta allora troppo piccola, e l’app la corregge:

χP2=∑n=1N(cn−e)2e,e=kDN,χ2=χP2 N−1N−k ∼ χN−12\chi^2_P = \sum_{n=1}^{N} \frac{(c_n - e)^2}{e}, \quad e = \frac{kD}{N}, \qquad \chi^2 = \chi^2_P\,\frac{N-1}{N-k} \ \sim\ \chi^2_{N-1}
Statistica di Pearson con la correzione per l’estrazione senza reinserimento; servono almeno 10 estrazioni.
p=Q(N−12, χ22)=Γ(N−12,χ22)Γ(N−12)p = Q\Bigl(\frac{N-1}{2},\ \frac{\chi^2}{2}\Bigr) = \frac{\Gamma\bigl(\frac{N-1}{2}, \frac{\chi^2}{2}\bigr)}{\Gamma\bigl(\frac{N-1}{2}\bigr)}
Funzione gamma incompleta superiore regolarizzata, calcolata con serie e frazione continua.
Valore pEsito
p>0.05p > 0.05nessuna deviazione significativa
0.01<p≤0.050.01 < p \le 0.05deviazione al limite
p≤0.01p \le 0.01deviazione statisticamente significativa
Verifica di equità con il χ² corretto, i gradi di libertà, il valore p e il test Q posizionale
Test di equità χ² e verifica posizionale

Test Q posizionale

Con almeno 30 estrazioni complete l’app verifica anche la media di ogni posizione ordinata (Coronel-Brizio et al., arXiv:0806.4595). L’ii-esimo più piccolo di kk numeri estratti da 1..N1..N ha media e covarianza

μi=(N+1) ik+1,Vij=i (k−j+1)(N+1)(N−k)(k+1)2(k+2)  (i≤j),Q=m δTV−1δ ∼ χk2\mu_i = \frac{(N+1)\,i}{k+1}, \qquad V_{ij} = \frac{i\,(k - j + 1)(N+1)(N-k)}{(k+1)^2 (k+2)}\ \ (i \le j), \qquad Q = m\,\boldsymbol\delta^{\mathsf T} V^{-1} \boldsymbol\delta \ \sim\ \chi^2_k
δ è il vettore delle differenze tra le medie osservate per posizione e μ su m estrazioni.

Coppie, confronti e verifiche

Simulare una vita di gioco

La simulazione su 50 anni usa la probabilità di almeno un jackpot in TT giocate indipendenti:

PT=1−(1−1C)T,T=50⋅52⋅d⋅rP_T = 1 - \Bigl(1 - \frac{1}{C}\Bigr)^{T}, \qquad T = 50 \cdot 52 \cdot d \cdot r
Probabilità di almeno un jackpot in T giocate; C è il numero di combinazioni del jackpot, d le estrazioni a settimana, r le righe per estrazione.

6. Analisi per calendario

L’analisi per calendario cerca coincidenze tra le date delle estrazioni e i numeri estratti — per esempio il giorno del mese che compare tra le palline — e confronta ogni frequenza osservata con quella che produrrebbe il puro caso.

Selezione delle estrazioni

Le estrazioni si possono filtrare per anno, mese, stesso mese in tutti gli anni o intervallo di date, e per epoca delle regole. Gli elenchi di frequenza (caldi, freddi, assenti da più tempo) usano i conteggi grezzi della selezione; un numero impossibile in tutte le epoche selezionate viene escluso, e una selezione che mescola epoche diverse viene segnalata con un avviso, perché un numero valido solo nell’epoca più recente sembra falsamente «freddo». Le regole di data qui sotto tengono conto delle epoche: un’estrazione conta per una regola solo se almeno un candidato è un numero valido secondo le regole proprie di quell’estrazione.

r^=∣{ t∈E:Ct∩Xt≠∅ }∣∣E∣,E={ t:Ct∩Ωt≠∅ }\hat r = \frac{\bigl|\{\, t \in \mathcal E : C_t \cap X_t \ne \varnothing \,\}\bigr|}{|\mathcal E|}, \qquad \mathcal E = \bigl\{\, t : C_t \cap \Omega_t \ne \varnothing \,\bigr\}
Frequenza osservata di una regola di data: C_t sono i numeri candidati dell’estrazione t, X_t i suoi numeri estratti, Ω_t i numeri validi secondo le sue regole; ℰ sono le estrazioni ammissibili.

Dieci regole di data

RegolaNumeri candidati per un’estrazione del GG.MM.AAAA
giornoGG
meseMM
annoAAAA mod 100
oral’ora dell’estrazione, dove l’operatore la pubblica
cifre dell’annole cifre distinte di AAAA
somma delle cifre dell’annosomma delle cifre di AAAA
somma dell’anno a due cifresomma delle cifre di AAAA mod 100
somma delle cifre del giornosomma delle cifre di GG
somma delle cifre della datasomma di tutte le cifre di GG, MM e AAAA
radice numericala somma delle cifre della data ridotta a una sola cifra

Un candidato pari a 0 viene ignorato, e un’estrazione conta una sola volta per una regola anche quando più dei suoi numeri corrispondono.

Attesa casuale

Per un’estrazione di kk numeri da SS valori possibili e un insieme di mm numeri candidati, la probabilità che non esca nessun candidato è ipergeometrica:

P0=∏i=0k−1S−m−iS−i=(S−mk)(Sk),P≥1=1−P0P_0 = \prod_{i=0}^{k-1} \frac{S - m - i}{S - i} = \frac{\binom{S-m}{k}}{\binom{S}{k}}, \qquad P_{\ge 1} = 1 - P_0
P₀: non esce nessun candidato; P≥1: ne esce almeno uno. I candidati fuori dall’intervallo dell’epoca dell’estrazione vengono scartati.
E[r^]=1∣E∣∑t∈EP≥1(t),P∪=1−(1−PM)(1−PB)\mathbb E[\hat r] = \frac{1}{|\mathcal E|}\sum_{t \in \mathcal E} P_{\ge 1}^{(t)}, \qquad P^{\cup} = 1 - \bigl(1 - P^{\mathrm M}\bigr)\bigl(1 - P^{\mathrm B}\bigr)
Frequenza attesa di corrispondenza r̂ sulle estrazioni ammissibili ℰ, ciascuna con la propria epoca; P^M e P^B per l’insieme principale e quello bonus vengono combinati come indipendenti.
Corrispondenze di calendario con le frequenze osservate accanto all’attesa casuale
Corrispondenze di calendario confrontate con il caso

Interpretare le frequenze

Le coppie vengono elencate da 30 estrazioni e le terne da 120 (solo le combinazioni viste almeno due volte, le dieci più frequenti); sotto le 100 estrazioni compare un avviso di campione piccolo. Con 100 estrazioni, una frequenza vicina al 20% ha ancora una banda al 95% di circa ±1.960.2⋅0.8/100≈±7.8\pm 1.96\sqrt{0.2 \cdot 0.8/100} \approx \pm 7.8 punti percentuali.

7. Archivio storico ed epoche delle regole

Le lotterie cambiano le loro regole. Un’estrazione ha senso solo insieme alle regole in vigore alla sua data, per questo ogni estrazione memorizzata porta con sé la sua versione delle regole.

Epoche delle regole

LotteriaEpoche (principali · bonus)
Lotto (Norvegia)7/34 + 3 numeri supplementari → 7/34 + 1 supplementare (dal 21.03.2015)
Vikinglotto6/48 + supplementare e numero fortunato → numero Viking 1–8 (2017) → 1–5 (dal 09.06.2021)
Eurojackpot5/50 + 2/8 → 2/10 (2014) → 2/12 (dal 25.03.2022)
Powerball5/59 + 1/39 → 5/59 + 1/35 (2012) → 5/69 + 1/26 (dal 07.10.2015)
Mega Millions5/52 + 1/52 → 5/56 + 1/46 → 5/75 + 1/15 → 5/70 + 1/25 → 5/70 + 1/24 (dal 08.04.2025)
EuroMillions5/50 + 2/9 → 2/11 (2011) → 2/12 (dal 27.09.2016)
SuperEnalotto6/90 + Jolly dal 1997; il SuperStar proviene da un’urna separata ed è escluso dai modelli
Lotto Max7/49 → 7/50 (2019) → 7/52 (dal 14.04.2026)
Powerball AU5/45 + 1/45 → 6/40 + 1/20 (2013) → 7/35 + 1/20 (dal 19.04.2018)

Regola di compatibilità

Modelli e analisi usano solo estrazioni comparabili con le regole attuali. Due epoche vengono confrontate gruppo per gruppo:

(k,nmin⁡,nmax⁡)e=(k,nmin⁡,nmax⁡)now,Be=Bnow(k, n_{\min}, n_{\max})_e = (k, n_{\min}, n_{\max})_{\mathrm{now}}, \qquad \mathcal B_e = \mathcal B_{\mathrm{now}}
L’epoca e è comparabile per i numeri principali quando quantità e intervallo coincidono (a sinistra), e per i numeri bonus quando coincide ogni gruppo bonus 𝓑 (a destra). Un’estrazione con numeri principali comparabili ma regole bonus diverse conserva i numeri principali e scarta quelli bonus.
Storico dell’Eurojackpot con la scheda dell’epoca delle regole attuali e le estrazioni contrassegnate con la loro versione delle regole
L’archivio contrassegna ogni estrazione con la sua epoca delle regole

8. Probabilità e limiti

Ogni estrazione ufficiale è un esperimento casuale indipendente. Nessuno storico, modello, tribunale o simulazione cambia la probabilità di alcuna combinazione.

Probabilità del jackpot

PJ=1(Nk)(Bb),1(505)(122)=12 118 760×66=1139 838 160P_J = \frac{1}{\binom{N}{k}\binom{B}{b}}, \qquad \frac{1}{\binom{50}{5}\binom{12}{2}} = \frac{1}{2\,118\,760 \times 66} = \frac{1}{139\,838\,160}
Probabilità del jackpot di una riga con le regole attuali (a destra: Eurojackpot); le lotterie senza numero bonus per il giocatore usano solo il primo fattore.
LotteriaRegoleProbabilità del jackpot
Lotto (Norvegia)7/341 : 5 379 616
Vikinglotto6/48 + 1/51 : 61 357 560
Eurojackpot5/50 + 2/121 : 139 838 160
EuroMillions5/50 + 2/121 : 139 838 160
Powerball5/69 + 1/261 : 292 201 338
Mega Millions5/70 + 1/241 : 290 472 336
SuperEnalotto6/901 : 622 614 630
Lotto Max7/521 : 133 784 560
Powerball AU7/35 + 1/201 : 134 490 400

Ogni categoria di premio

Indovinare esattamente ii dei kk numeri principali e jj dei bb numeri bonus ha come probabilità il prodotto di due probabilità ipergeometriche:

P(i,j)=(ki)(N−kk−i)(Nk)⋅(bj)(B−bb−j)(Bb)P(i, j) = \frac{\binom{k}{i}\binom{N-k}{k-i}}{\binom{N}{k}} \cdot \frac{\binom{b}{j}\binom{B-b}{b-j}}{\binom{B}{b}}
Per una pallina aggiuntiva estratta dalla stessa urna (Lotto, SuperEnalotto, Lotto Max) il secondo fattore è (k−i)/(N−k) per «+1».
Categoria EurojackpotProbabilità
5 + 21 : 139 838 160
5 + 11 : 6 991 908
5 + 01 : 3 107 515
4 + 21 : 621 503
4 + 11 : 31 075
3 + 21 : 14 125
4 + 01 : 13 811
2 + 21 : 985
3 + 11 : 706
3 + 01 : 314
1 + 21 : 188
2 + 11 : 49

Perché nessun metodo può battere il caso

Se le estrazioni sono indipendenti e ogni combinazione è ugualmente probabile, la probabilità condizionata della prossima estrazione dato qualsiasi storico è la stessa che senza di esso, e il numero atteso di numeri principali indovinati è lo stesso per ogni riga, comunque sia stata prodotta:

P(Xt+1=A∣X1,…,Xt)=P(Xt+1=A)=(Nk)−1,E[H]=k⋅kNP\bigl(X_{t+1} = A \mid X_1, \dots, X_t\bigr) = P(X_{t+1} = A) = \binom{N}{k}^{-1}, \qquad \mathbb E[H] = k \cdot \frac{k}{N}
H è il numero di numeri principali indovinati da una riga. Eurojackpot: 5 · 5/50 = 0.5 per riga con qualsiasi metodo.

Apri Lotto Simulator