print [ len(filter(lambda x: x > 0.6 *per_day, [ sum([ randint(0,1) for birth in range(0, per_day) ]) for day in range(0, 365) ])) for per_day in (15, 45) ]
wedrifid comments on 2012 Less Wrong Census/​Survey