tmp006.py

import random

def proc(n):
    s1 = 0
    lims = [3, 3, 0]
    for _ in range(n):
        for i in range(len(lims)):
            s = 0
            while True:
                if random.choice((0, 0, 1)):
                    break
                s += 1
            if s >= lims[i]:
                s1 += s
                break
    return s1


n = 1000000
m = proc(n)
print(n, m, m / n)