import random
xwins = 0
owins = 0
draws = 0
wins = [{1,2,3},{4,5,6},{7,8,9},
{1,4,7},{2,5,8},{3,6,9},
{1,5,9},{3,5,7}]
nums = [1,2,3,4,5,6,7,8,9]
durations = [0 for i in range(9)]
iters = 1000000
for iter in range(iters):
done = False
tics = random.sample(nums, 5)
for w in wins:
if w&set(tics[0:5:2]) == w:
xwins += 1
done = True
durations[4] += 1
break
if done:
continue
remains = [n for n in nums if n not in tics]
tics.append(random.choice(remains))
for w in wins:
if w&set(tics[1:6:2]) == w:
owins += 1
done = True
durations[5] += 1
break
if done:
continue
remains = [n for n in nums if n not in tics]
tics.append(random.choice(remains))
for w in wins:
if w&set(tics[0:7:2]) == w:
xwins += 1
done = True
durations[6] += 1
break
if done:
continue
remains = [n for n in nums if n not in tics]
tics.append(random.choice(remains))
for w in wins:
if w&set(tics[1:8:2]) == w:
owins += 1
done = True
durations[7] += 1
break
if done:
continue
remains = [n for n in nums if n not in tics]
tics.append(random.choice(remains))
for w in wins:
if w&set(tics[0:9:2]) == w:
xwins += 1
done = True
durations[8] += 1
break
if done:
continue
draws += 1
mysum = sum((xwins,owins,draws))
print(xwins/mysum,owins/mysum,draws/mysum )
othersum = sum(durations)
for i in range(4,9):
print('{}: {}'.format(i+1, round(durations[i]/othersum,2)))
Simulation
