Sezione di VALCON.IT
Funzioni per generare / manipolare numeri pseudocasuali random.choice(z) Genera un elemento casuale scegliendolo tra quelli appartenenti alla struttura dati z random.choices(z, k=n) Genera n elementi casuali (con reinserimento) scegliendoli tra quelli appartenenti alla struttura dati z random.randint(a, b) Genera un numero intero nell’intervallo [a, b] random.random() Genera un numero reale nell’intervallo [0, 1) random.randrange(stop) Come … Leggi tutto
Soluzione 1 La funzione shuffle() della libreria random risolve il problema
|
1 2 3 4 5 6 |
import random # shuffle() NUMERI=[1,9,2,8] random.shuffle(NUMERI) # [8, 2, 9, 1] random.shuffle(NUMERI) # [9, 8, 1, 2] |
Soluzione 2 Osserva n volte scambia l’elemento i-esimo con un elemento a caso
|
1 2 3 4 5 6 7 8 9 10 11 12 |
import random # randint() NUMERI = [...] # Inserisci i numeri def mescola(lista): n=len(lista) # Quanti scambi? for i in range(n): p=random.randint(0,n-1) # posizione temp = lista[i] # scambia lista[i] = lista[p] lista[p] = temp |
Soluzione 3 Osserva n volte scambia una coppia di elementi a caso
|
1 2 3 4 5 6 7 8 9 10 11 12 13 |
import random # randint() NUMERI = [...] # Inserisci i numeri def mescola(lista): n=len(lista) # Quanti scambi? for i in range(n): p1=random.randint(0, n-1) # 1° posizione casuale p2=random.randint(0, n-1) # 2° posizione casuale temp = lista[p1] # scambia lista[p1] = lista[p2] lista[p2] = temp |
La scambia può essere scritta più compatta
|
1 2 3 |
... lista[p1], lista[p2] = lista[p2], lista[p1] ... |
Soluzione 1 La lista di numeri esiste già
|
1 2 3 4 5 6 7 8 |
import random # random() NUMERI=[...] # Inserisci i numeri n=len(NUMERI) for i in range(n): x=random.random() # Numero casuale in [0, 1) NUMERI[i]=x |
Soluzione 2 La lista di numeri è vuota
|
1 2 3 4 5 6 7 8 |
import random # random() NUMERI=[] # Lista vuota n =... # Quanti numeri inserire for i in range(n): x=random.random() NUMERI.append(x) |
Soluzione 3 I numeri devono essere in un certo intervallo
|
1 2 3 4 5 6 7 8 9 10 |
import random # uniform() NUMERI =[] # Lista vuota minimo =... # Valore minimo massimo=... # Valore massimo n =... # Quanti numeri inserire for i in range(n): x=random.uniform(minimo, massimo) NUMERI.append(x) |
Soluzione 4 I numeri devono essere interi in un certo intervallo
|
1 2 3 |
... x=random.randint(minimo, massimo) ... |
(Vedi la discussione) Simulazione del gioco con 3 dadi “Carnival game”
|
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 |
import random def dado(): return random.randint(1,6) VOLTE=10 bilancio=0 x=6 # Numero frotunato! for i in range(VOLTE): dado1=dado() dado2=dado() dado3=dado() vincita=0 if(dado1 == x): vincita +=1 if(dado2 == x): vincita +=1 if(dado3 == x): vincita +=1 if(vincita == 0): bilancio -= 1 else : bilancio += vincita print(dado1, dado2, dado3, "-", vincita, "-", bilancio) spe_mat=bilancio/VOLTE print(spe_mat, 100*spe_mat, "%") |
Per avere risultati significativi aumenta il valore di VOLTE elimina la print() nel for matplotlib Osserva il bilancio finale all’aumentare del numero di giocate…
|
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 |
import random import matplotlib.pyplot as plt def dado(): return random.randint(1,6) VOLTE=10000 # 100 1000 10000 ... bilancio=0 BILANCIO=[] x=6 # Numero fortunato! for i in range(VOLTE): dado1=dado() dado2=dado() dado3=dado() vincita=0 if(dado1 == x): vincita +=1 if(dado2 == x): vincita +=1 if(dado3 == x): vincita +=1 if(vincita == 0): bilancio -= 1 else : bilancio += vincita BILANCIO.append(bilancio) plt.grid (which="major") plt.plot (BILANCIO) plt.title("Carnival Game") plt.show () |
Olimpiadi Italiane di Statistica 2011 E.M. – 2 Qual è la probabilità che in una classe di 25 alunni, almeno due compiano gli anni lo stesso giorno? Con formula Vedi la discussione La probabilità con 25 persone è superiore al 50%
|
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 |
import matplotlib.pyplot as plt PERSONE=25 ANNO=365 pNo =(PERSONE+1)*[0] pSì =(PERSONE+1)*[0] pNo[0]=1 pSì[0]=0 pNo[1]=1 pSì[1]=0 for n in range(2, PERSONE+1): pNo[n] = pNo[n-1]*(ANNO-(n-1))/ANNO pSì[n] = 1-pNo[n] for n in range(0, PERSONE+1): print("%2i %f %f" %(n, pNo[n],pSì[n])) N=range(0, PERSONE+1) plt.grid() plt.plot(N,pNo) plt.plot(N,pSì) plt.legend(["No", "Sì"]) plt.title("Compleanni coincidenti") plt.xlabel("Numero persone") plt.ylabel("Probabilità") #--------------------------------------- plt.show() |
Aumentando il numero di persone a 60 si può notare che la probabilità … Leggi tutto