#!/usr/bin/env python

import random

# number of iterations per trial
numiter = 100000

# range and probability to test
r = 1000
p = 46

for num in range(100):
	count = 0
	for i in range(numiter):
		for s in range(num):
			n = random.randrange(r)
			if (n <= p):
				count = count + 1
				break

	print "n = " + str(num) + " prob " + str(float(count) / numiter)


