CCCPaste

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)