裏表(Phinloda のもう裏だか表だか分からないページ)

コンピュータ・プログラミング系の話がメインのそれなりにごちゃごちゃしたネタばかり出てくるサイトです。多分。
<< スマホのブラウザがまた重くなってきたのだが | Top | 「Phinloda のいつか聴いた曲」にヘンな文字列が入っていたのだが >>
45問の四択のテストで適当に選んで30問以上正解する確率は

計算したら約 0.000000475% になった、多分合っていると思うのだけど、これがうまくプログラムで書けない。

import scipy.misc as scm
import numpy as np

p = 0
n = 45
for i in range(30, n + 1):
    p += scm.comb(n, i, 1) * np.power(1/4, i) * np.power(3/4, n - i)

print(p)
4.75006548988e-09

numpy 使うのだから、ループさせないで一撃で計算できそうな気がするけど、vectorize とか使って定義した関数を要素毎に呼び出して sum 取ったりしているとあまり見た目は変わらない(笑)。もうちょっと真面目に numpy 勉強しないとダメっぽい。

JUGEMテーマ:日記・一般

| Python | 18:12 | comments(0) | trackbacks(0)
スポンサーサイト
| - | 18:12 | - | -
コメント
コメントする









この記事のトラックバックURL
http://phinloda.jugem.cc/trackback/4389
トラックバック
Powered by "JUGEM"
▲このページの先頭へ
CALENDAR
S M T W T F S
1234567
891011121314
15161718192021
22232425262728
2930     
<< September 2019 >>
NEW ENTRIES
CATEGORIES
ARCHIVES
NEW COMMENTS
NEW TRACKBACKS
LINKS
PROFILE