1

我在做生日悖论，想知道有多少人可以通过使用 python 满足两个人生日相同的 0.5 概率。

我尝试不使用数学公式通过在 python 中使用 random 和 randint 来查找给定人数的概率

import random
def random_birthdays():
    bdays = []
    bdays = [random.randint(1, 365) for i in range(23)]
    bdays.sort()
    for x in range(len(bdays)):
        while x < len(bdays)-1:
            print x
            if bdays[x] == bdays[x+1]:
                #print(bdays[x])
                return True
            x+=1
        return False
count  = sum(random_birthdays() for _ in range(1000))
print('In a sample of 1000 classes each with 23 pupils, there were', count, 'classes with individuals with the same birthday')

我期待一些提示或代码可以帮助我解决这个问题。

4

1 回答 1

0

好吧，你的代码有问题，你只检查连续的生日相等。最好使用集合检查它

沿线

import random

def sim(n):
    """simulate birthdays for n people"""
    a = set([random.randint(1, 365) for _ in range(n)])
    if len(a) == n:
        return False
    return True

print(sim(23))
print(sim(23))
print(sim(23))
print(sim(23))
print(sim(23))

如果人们的生日是同一天，上面的函数将返回 true，否则返回nfalse。

对于 n = 20, 21, ...., 25 调用它 1000000 次并计算有多少 True 与 False

运行代码

nt = 0
nf = 0
n = 23
for k in range(0, 1000000):
    if sim(n):
        nt += 1
    else:
        nf += 1

print((nt, nf))

对于 n = 23 和 n = 22 产生

(506245, 493755)
(475290, 524710)
于 2019-09-06T11:52:56.113 回答