見るだけです。ターボ推奨です。(シフト+スタート) モンテカルロ法でサンクトぺテルブルのパラドックスを実験します。 最高賞金が出たところを赤い✖で示してます。 対数グラフです。 ラインは下から5,10,15,20,30,50,100,200をあらわしてます。 私の期待値最高記録は5074.088218でした。 ---------サンクトぺテルブルのパラドックスとは---------- 簡単に言うと、 1/2の確率で1円の賞金 1/4の確率で2円の賞金 1/8の確率で4円の賞金 1/16の確率で8円の賞金 1/(2^n)の確率で2^(n-1)円の賞金 がもらえる賭けがあって、その期待値を計算すると 1/2 + 2/4 + 4/8 + 8/16 + 16/32 + 32/64 + 64/128 + ... となって∞になります。 理論上この賭けで平均∞円もらえるので、参加費がいくらであってもこの賭けに参加すべきなのか、というものです。 でもこの賭けで10000円獲得できる確率はたった32768分の1なので、なんか変ですよね、というのがサンクトぺテルブルのパラドックスです。 ------------------------------------------------------------------------- ちなみに私（四角形ABCD）は賞金40億を出したことがあります。80億分の1です。多分私の人生ではもう二度と良いことは起きないでしょう。 最高値は試行回数の半分程度になることが多いです。何度かやると最初に期待値70くらいが出ることがありますが、試行回数が少ないためです。 2020/10/24 19:17今やってみたんですけど126,797回目に33,554,332を出して期待値が273.891299になりました。これよりでかい最高期待値がでたらコメントお願いします！まあ一応464.158883までグラフに入るので振り切れることはないと思います。 2020/10/24 19:29全力で前言撤回します。27,349回目に67,108,864を出して期待値が2460.739735になりました。普通に振り切れましたごめんなさい。 2020/10/24 19:32 1,664回目に8,388,608 期待値5074.088218 ちなみに1,664回目に少し上のほうに書いてある40億(4,294,967,296)が出てたら期待値はなんと2,581,110.153846になります。
2020/10/24 ver1.1.4 軽微なバグの修正・グラフの改善 2020/10/24 ver1.1.3 ver1.1.2で修正したはずのバグが修正されていなかったので再修正 2020/10/24 ver1.1.2 軽微なバグの修正 2020/10/24 ver1.1.1 最高期待値の重大なバグの修正 2020/10/24 ver1.1.0 バグの修正・最高期待値の追加 2020/10/24 ver1.0.0
