super pi を作ろう - pyopyopyo - Linuxとかプログラミングの覚え書き -

円周率πを100万桁まで計算するプログラムを作ってみました。super pi もどきです。ソースコードと解説を載せます。 ソースコードは80行程度の非常に短いものですが、このような簡単な実装でも最近のPCなら100万桁の円周率が一分も掛らず計算できます。 ポイント 要点は、 計算精度が問題となるので、GMP(任意精度計算ライブラリ)を使う 計算アルゴリズムは、Square AGM *1と呼ばれるものを使う です。 GMPについては、http://d.hatena.ne.jp/pyopyopyo/20090303 をご覧下さい。 ソースコード /** * pi.cpp * * written by pyopyopyo at gmail dot com */ #include <math.h> #include <gmpxx.h> #include <iostream> // N桁の円周率を計