.ipynb .pdf

Binder

Contents

import random

N = 1000000
wins = 0
for i in range(N):
    A = int(6*random.random()+1)
    B = int(6*random.random()+1)
    C = int(6*random.random()+1)
    if A == B == C:
        wins = wins+A
    elif A+B+C <= 6:
        wins = wins-(A+B+C)
print(wins, N, end=' ')
print(1.0*wins/N)