pythonのrandom - rougeref’s diary

高体連派遣の合間に1から50までの整数から20個の整数をランダムに選択するスクリプトを作成。なぜにエントリのカテゴリがrefereeになっているかというと、それは11月14日、15日に行われる某会で用意するなにかに関係があるからです。

最初に作成したスクリプトはこんなの。

#!/usr/bin/python

import random
b=[]
while 1:
    a=int(50*random.random())+1
    if a not in b:
        b.append(a)
    if(len(b) == 20):
        break
b.sort()
print  b

実際に動かしてみるとこんな結果。

$ ./mmm.py 
[4, 5, 9, 11, 15, 16, 19, 20, 21, 22, 25, 27, 28, 30, 32, 37, 38, 39, 41, 42]
$ ./mmm.py 
[4, 7, 9, 12, 13, 14, 16, 20, 21, 23, 24, 25, 26, 30, 31, 38, 43, 44, 45, 47]
$ ./mmm.py 
[1, 2, 5, 6, 8, 10, 11, 15, 17, 20, 25, 26, 27, 30, 32, 33, 37, 38, 46, 47]

ん〜、いいのかな。なんだか偏りがあるような気がしませんか？
ということでさっそく検証してみる。上記のスクリプトを100万回ループさせて、各数字毎に選択された回数を数える。100万回×2/5=40万回くらいが平均になればいいはず。

検証用スクリプトはこんなの。

#!/usr/bin/python

import random
c={}
while i < 1000000:
    b=[]
    while 1:
        a=int(50*random.random())+1
        if a not in b:
            b.append(a)
            if c.has_key(a):
                c[a]=c[a]+1
            else:
                c[a]=1
        if(len(b) == 20):
            break
    i=i+1

tmp_list=c.items()
tmp_list.sort(key=lambda x:x[1])

i=0
while i < 50:
    print "%02d:%d" % (tmp_list[i][0],tmp_list[i][1])
    i=i+1

技術的には別に難しくありません。ramdomで1から50までの整数を作成する度に、ディクショナリのカウントアップをして、valueでsortして出力。で結果はこんな。

もっとも多く選ばれた"38"ともっとも少なく選ばれた"47"の差は1.0045倍。

もう一回。

こんどは1.0057倍。

ふむ、少なくとも日本の衆議院選挙における一票の格差よりは差分が少ないようなので、偏りはないようです。ということで、これで20問選択して作成します。