Rmpfr

ついでに、R の任意精度演算パッケージ（Rmpfr）。たとえば、

「2ⁿ - 2 が n で割り切れたら n は素数である」
という命題を用いて、10000までの素数をすべて求める

library(Rmpfr)
two <- mpfr(2,10000)
for(i in 1:10000) if( (two^i-two)%%i ==0 ) cat( i," " )

なお、このパッケージも、某サーバ（CentOS 5）上の R へ入れようとすると、gmp（GNU MP）がない、GNU MPFR がないとおこられる。それぞれを wget して make install （GMP, GNU MPFR）の後に、R から install.packages( ) 。