python　シミュレーション　基礎的な問題です。

学校の課題で衝突確率を計算し、グラフに表示するプログラムを作成しています。
解析計算プログラムとシミュレーションプログラム、グラフ作成プログラムを作成。
このうち、解析プログラムは衝突する確率を計算（これは完成できました）
シミュレーションプログラムは、プログラムの機能をM回実行、実行するごとに重複の有無を判定、最終的に「M回中何回重複したか」の回数を「M」で割れば、重複確率が求められる。（ここがわかりません）

この、シミュレーションプログラムの作成方法がわかりません。
どのように考えたらいいのでしょうか。
よろしくお願いします。

（解析プログラム）
import datetime
from collections import Counter

"""
birthday paradox
"""


def rand_date():
    """
    ランダムな日付を返す
    戻り値: (2010/1/1 + d)を計算し(月, 日)を返す
        - base_date(基準日): 2010/1/1 + d
        - d(進める日数): (0〜364)
    """
    # 2010は日付を生成するための仮の年. うるう年でなければいつでもよい
    base_date = datetime.datetime(2010, 1, 1)
    d = random.randint(0, 365)
    ret_date = base_date + datetime.timedelta(days=d)
    return ret_date.month, ret_date.day

if __name__ == "__main__":
    # 生成する誕生日の個数をセット
    if len(sys.argv) != 2:
        print("usage: rand.py [nOfSample]")
        exit(0)

    n = int(sys.argv[1])

    # 誕生日をn個生成
    birthday_list = [rand_date() for i in range(n)]

    # 結果を検証

    # 重複する誕生日を検索(重複なしの場合は空リスト)
    non_unique_items = [k for k, v in Counter(birthday_list).items() if v > 1]

    num = 1
    for i in range(1,n):
      num *= (365-i)/365

    probability = (1-num)
    # 結果を表示
    print("# non_unique_items")
    print(non_unique_items)
    print("# probability(unique)")
    print(1-num)
    print("# all items")
    print(birthday_list)

（プログラム）
#! /usr/ bin /env python
# −*− coding : utf−8 −*−
import sys
import random
import datetime
from collections import Counter
” ” ”
birthday paradox
” ” ”
def rand_date ( ) :
” ” ”
ランダムな日付を返す
戻り値: (2010/1/1 + d) を計算し( 月, 日) を返す
− base_date ( 基準日) : 2010/1/1 + d
− d( 進める日数) : (0～364)
” ” ”
# 2010は日付を生成するための仮の年. うるう年でなければいつでもよい
base_date = datetime . datetime(2010, 1 , 1)

d = random . randint (0 , 365)
ret_date = base_date + datetime . timedelta ( days=d)
return ret_date . month, ret_date . day
i f __name__ == ” __main__ ” :
# 生成する誕生日の個数をセット
i f len ( sys . argv ) != 2:
print ( ”usage : rand . py [nOfSample] ” )
exit (0)
n = i n t ( sys . argv [ 1 ] )
# 誕生日をn個生成
birthday_list = [ rand_date ( ) for i in range (n ) ]
# 結果を検証
# 重複する誕生日を検索( 重複なしの場合は空リスト)
non_unique_items = [ k for k , v in Counter ( birthday_list ) . iteritems ( ) i f v > 1]
# 結果を表示
print ( ”# non_unique_items ” )
print ( non_unique_items )
print ( ”# a l l items ” )
print ( birthday_list )