高校数学の質問スレ(医者・東大卒専用) Part438 (991レス)
前次1-
抽出解除 必死チェッカー(本家) (べ) 自ID レス栞 あぼーん

641: 132人目の素数さん [sage] 02/01(土)23:09 ID:RwePvY7v(1/2)
import random
from collections import Counter
import statistics

black_numbers = {2, 4, 6, 8, 10, 11, 13, 15, 17, 20, 22, 24, 26, 28, 29, 31, 33, 35}

def sim():
bribe = 0
a = random.randint(0, 36)
while a in black_numbers:
bribe += a
a = random.randint(0, 36)
return bribe + a

# 1000000回のシミュレーションを実行
results = [sim() for _ in range(1000000)]

# 最頻値を求める
frequencies = Counter(results)
most_frequent = frequencies.most_common(1)[0]

# 最頻値を表示
print(f"最頻値: {most_frequent[0]} ドル")

# 中央値を求める
median_value = statistics.median(results)

# 中央値を表示
print(f"中央値: {median_value} ドル")

# 期待値を求める
mean_value = statistics.mean(results)

# 期待値を表示
print(f"期待値: {mean_value} ドル")
642: 132人目の素数さん [sage] 02/01(土)23:30 ID:RwePvY7v(2/2)
>>638
シミュレーション結果と乖離
https://i.imgur.com/mxcZRmN.png


100以下の確率96.3%
95%以上の確率での釈放に必要な賄賂は92ドル
となりました。

シミュレーションのコード
black={2, 4, 6, 8, 10, 11, 13, 15, 17, 20, 22, 24, 26, 28, 29, 31, 33, 35};

sim[] :=Module[{bribe,a},
bribe=0;
a=RandomChoice[Range[0,36],1][[1]];
While[ContainsAny[black,{a}],
bribe+=a;
a=RandomChoice[Range[0,36],1][[1]]
];
bribe+a
]
前次1-
スレ情報 赤レス抽出 画像レス抽出 歴の未読スレ AAサムネイル

ぬこの手 ぬこTOP 0.038s