# Afrikan tähti: Kairo 4 askeleen päässä
# Pitkä matematiikka kevät 2021 tehtävä 7
# Kuinka monta nopanheittoa tarvitaan Kairoon
# Juhani Kaukoranta 26.3.2021
from random import randint
max = 1000 # simulointikertojen lukumäärä
heittokertoja = 0
for simulointi in range(0,max): 
    heittoja = 0 
    summa = 0
    for k in range(0,4): # enintään 4 heittoa riittää
        n = randint(1,6)
        summa += n # noppien silmälukujen summa
        if summa >= 4: #päästiin Kairoon
            heittoja += k+1 # tarvittujen heittojen määrä
            heittokertoja += heittoja # kumulatiivinen keskiarvo varten
            break # päästiin Kairoon
print("Kairoon pääsyn heittokertojen keskiarvo = ",heittokertoja/max)
print("Teoreettinen odotusarvo 343/216 = 1.588")