学校の課題で衝突確率を計算し、グラフに表示するプログラムを作成しています。
解析計算プログラムとシミュレーションプログラム、グラフ作成プログラムを作成。
このうち、解析プログラムは衝突する確率を計算(これは完成できました)
シミュレーションプログラムは、プログラムの機能をM回実行、実行するごとに重複の有無を判定、最終的に「M回中何回重複したか」の回数を「M」で割れば、重複確率が求められる。(ここがわかりません)
この、シミュレーションプログラムの作成方法がわかりません。
どのように考えたらいいのでしょうか。
よろしくお願いします。
python
1(解析プログラム) 2import datetime 3from collections import Counter 4 5""" 6birthday paradox 7""" 8 9 10def rand_date(): 11 """ 12 ランダムな日付を返す 13 戻り値: (2010/1/1 + d)を計算し(月, 日)を返す 14 - base_date(基準日): 2010/1/1 + d 15 - d(進める日数): (0〜364) 16 """ 17 # 2010は日付を生成するための仮の年. うるう年でなければいつでもよい 18 base_date = datetime.datetime(2010, 1, 1) 19 d = random.randint(0, 365) 20 ret_date = base_date + datetime.timedelta(days=d) 21 return ret_date.month, ret_date.day 22 23if __name__ == "__main__": 24 # 生成する誕生日の個数をセット 25 if len(sys.argv) != 2: 26 print("usage: rand.py [nOfSample]") 27 exit(0) 28 29 n = int(sys.argv[1]) 30 31 # 誕生日をn個生成 32 birthday_list = [rand_date() for i in range(n)] 33 34 # 結果を検証 35 36 # 重複する誕生日を検索(重複なしの場合は空リスト) 37 non_unique_items = [k for k, v in Counter(birthday_list).items() if v > 1] 38 39 num = 1 40 for i in range(1,n): 41 num *= (365-i)/365 42 43 probability = (1-num) 44 # 結果を表示 45 print("# non_unique_items") 46 print(non_unique_items) 47 print("# probability(unique)") 48 print(1-num) 49 print("# all items") 50 print(birthday_list) 51
python
1(プログラム) 2#! /usr/ bin /env python 3# −*− coding : utf−8 −*− 4import sys 5import random 6import datetime 7from collections import Counter 8” ” ” 9birthday paradox 10” ” ” 11def rand_date ( ) : 12” ” ” 13ランダムな日付を返す 14戻り値: (2010/1/1 + d) を計算し( 月, 日) を返す 15− base_date ( 基準日) : 2010/1/1 + d 16− d( 進める日数) : (0~364) 17” ” ” 18# 2010は日付を生成するための仮の年. うるう年でなければいつでもよい 19base_date = datetime . datetime(2010, 1 , 1) 20 21d = random . randint (0 , 365) 22ret_date = base_date + datetime . timedelta ( days=d) 23return ret_date . month, ret_date . day 24i f __name__ == ” __main__ ” : 25# 生成する誕生日の個数をセット 26i f len ( sys . argv ) != 2: 27print ( ”usage : rand . py [nOfSample] ” ) 28exit (0) 29n = i n t ( sys . argv [ 1 ] ) 30# 誕生日をn個生成 31birthday_list = [ rand_date ( ) for i in range (n ) ] 32# 結果を検証 33# 重複する誕生日を検索( 重複なしの場合は空リスト) 34non_unique_items = [ k for k , v in Counter ( birthday_list ) . iteritems ( ) i f v > 1] 35# 結果を表示 36print ( ”# non_unique_items ” ) 37print ( non_unique_items ) 38print ( ”# a l l items ” ) 39print ( birthday_list ) 40
回答2件
あなたの回答
tips
プレビュー
バッドをするには、ログインかつ
こちらの条件を満たす必要があります。