Question

The birthday problem is as follows: given a group of n people in a room, what is the probability that two or more of them have the same birthday? It is possible to determine the answer to this question via simulation. Using the starting template provided to you, complete the function called calc birthday probability that takes as input n and returns the probability that two or more of the n people will have the same birthday. To do this, the function should create a list of size n and generate n birthdays in the range 1 to 365 randomly, inclusive of the end-points 1 and 365. It should then check to see if any of the n birthdays are identical. The function should perform this experiment 106 times and calculate the fraction of time during which two or more people had the same birthday. The function will be called as follows from your main program:

IN PYTHON PLEASE

import random

def calc_birthday_probability (num_people):

random.seed (2020) # Don't change this value

num_trials = 1000

probability = 0

""" FIXME: Complete this function to return the probability of two or more people in the room having the same birthday. """

return probability

Answer 1

Thanks for the question. Below is the code you will be needing. Let me know if you have any doubts or if you need anything to change. 

Let me know for any help with any other questions.

Thank You !
===========================================================================

import random

def calc_birthday_probability(num_people):
    random.seed(2020)  # Don't change this value
    num_trials = 1000
    num_duplicates=0

    for attempt in range(int(num_trials)):
        birthday_list = [0]*366
        for person in range(num_people):
            birthday = random.randint(1,365)
            birthday_list[birthday]+=1

        for day in birthday_list:
            if day>1:
                num_duplicates+=1
                break

    return num_duplicates/num_trials

def main():
    num_people = 10
    p = calc_birthday_probability(num_people)
    print('Num People: {}, Probability: {:.5f}'.format(num_people,p))
main()

=================================================================