Math Battle [ 0336: πの近似計算 ]

[ 0336: πの近似計算 ]


[ 丘品花志先生 (funny teacher) の出題 ]

昨年 Google が発表した円周率の計算結果 (約 31.4 兆桁) に関して
解説記事 を書きました。
計算には Chudnovsky の式 (1988) が使われています。

これを変形した次の式を使ってどこまで計算できるか試してみてください。

Fun with Maths and Python を参考にすれば
計算式がもっと簡単になります。


[ 大宙乗児君のコメント ]

なーんだ。

教えていただいたサイトを見たら
Python のソースコードと実行結果が全部載っているではありませんか。

とりあえず僕の PC でそのままトライしてみます。


[ 広世正憲君のコメント ]

僕もさっそく pi.py というソースコードにしてやってみました。

$ python pi.py
31415926535897932384626433832795028841971693993751
058209749445923078164062862089986280348253421172983
('chudnovsky: digits', 10, 'time', 1.7881393432617188e-05)
('chudnovsky: digits', 100, 'time', 6.318092346191406e-05)
('chudnovsky: digits', 1000, 'time', 0.0008819103240966797)
('chudnovsky: digits', 10000, 'time', 0.08639097213745117)
('chudnovsky: digits', 100000, 'time', 8.564399003982544)

速いですね。凄いです。

[ 0341: 次の記事 ]

[ 0340: 素数日プログラム ]

[ 0339: emirp 探し ]

[ 0338: 0331の解きなおし ]

[ 0337: ubuntu 20.04 へ upgrade ]

[ 0336: πの近似計算 ]

[ 0335: 0253 の振り返り ]

[ 0334: Data Camp Python ]

[ 0333: 原始ピタゴラス数をさがす ]

[ 0332: LaTex の使いかた ]

[ 0331: 青い三角形の面積 ]

[ 0330: 長方形の幅 ]

[ 0329: 三角関数の関数の最大最小 ]

[ 0328: IMO 過去問 9 ]

[ 0327: 1+2+3+...+無限 ]

[ 0326: 2分木 sort も含めた性能比較 ]

[ 0325: 最悪ケースでの sort 比較 ]

[ 0324: 0307の回答 ④ ]

[ 0323: 0307の回答 ③ ]

[ 0322: 0307の回答 ② ]

[ 0321: 0307の回答 ① ]

[ 0320: 前の記事 ]

[ トップページへ ]