誕生日のパラドックス

birthday_paradox.py

from random import randint

NUM = 1000


def main():

    days = 365

    for n in range(1, days + 1):
        ratio = sum(1 if f(days, n) else 0 for _ in range(0, NUM)) / NUM
        print("{}: {:1.4f}".format(n, ratio))


def f(days, n):
    values = []
    for i in range(0, n):
        v = randint(0, days)
        if v in values:
            return True
        values.append(v)
    return False


if __name__ == "__main__":
    main()