Tuesday's child problem

TCP.py

import random, time, collections, pprint
  
sample = collections.defaultdict(float)

now = time.time()

known_child = ('M', 2)

while sum(sample.values()) < 100000:
  genders = [random.choice('MF') for dummy in '..']
  birthtimes = [random.uniform(0, now) for dummy in '..']
  birthdays = [time.localtime(t).tm_wday for t in birthtimes]
  gendertimes = zip(genders, birthdays)
  if known_child in gendertimes:
      gendertimes.remove(known_child)
      sample[gendertimes[0]] += 1
      if (sum(sample.values()) % 10000) == 0:
          print sum(sample.values())

males,females = [sum(c for (g, t), c in sample.items() if g == gender)
               for gender in 'MF']

print males / (males + females)