Haskellで円周率1億桁を計算する、あるいは円周率計算にHaskellの多倍長整数の改良を見る - 純粋関数空間

きっかけ http://d.hatena.ne.jp/itchyny/20120304 この記事を見て、私も昔円周率を計算したことがあるのを思い出しました。 http://d.hatena.ne.jp/tanakh/20070506#p1 公式は同じくチュドノフスキーを使用、GHCのIntegerをそのまま用いて円周率を計算、 当時（5年前）のマシン（おそらくCore2あたり）で1億桁が15分ほどだったようです。 SPOJというサイトのコンテスト向けに、 $\sqrt{2}$ 200万桁計算の高速化 をCで書くとこから始めて、 ナイーブ乗算、カラツバ法、そしてFFT乗算の実装にいたり、円周率へ。 最初AGMを実装し、 古くて新しい方法、Arctan系の公式の実装を経て Binary Splitting Methodを知り、 チュドノフスキーの公式へ至りました。 さらにHaskellでの実装

[itchyny]

｢なぜこのコードで円周率が｣は, 私のエントリー http://d.hatena.ne.jp/itchyny/20120304/1330870932 と, 論文 http://bellard.org/pi/pi2700e9/pipcrecord.pdf に詳しい.

[l-_-ll]

数値計算

[y-kawaz]

何でこのコードで円周率が出てくるのか分からんｗｗ

[Jxck]

ソース短いなぁ。

[yamamucho]

コード短い

[fm315]

「Haskellで円周率1億桁を計算する、あるいは円周率計算にHaskellの多倍長整数の改良を見る」

[terazzo]

Karatsubaって日本人っぽい名前と思ったけどロシアの人なのか。唐鍔さん。