计算圆周率
#
使用FFT+AGM代码计算
#
apt install gcc make -y
wget https://www.kurims.kyoto-u.ac.jp/~ooura/pi_fftc6_src.tgz
tar -zxvf pi_fftc6_src.tgz
cd pi_fftc6/
sed -i "s|-fast|-Ofast|g" Makefile_64bit
sed -i 's|scanf("%lg", &x)|x=100000|g' pi_fftcw.c
make -f Makefile_64bit
nohup ./pi_cw &
cat pi.dat
四行代码求圆周率
#
#include <stdio.h>
int a=10000,b,c=2800,d,e,f[2801],g;
void main(){
for(;b-c;)f[b++]=a/5;
for(;d=0,g=c*2;c-=14,printf("%.4d",e+d/a),e=d%a)for(b=c;d+=f[b]*a,
f[b]=d%--g,d/=g--,--b;d*=b);}
下载
#
wget https://files.pilookup.com/pi/1-100000000.zip
wget https://files.pilookup.com/pi/100000001-200000000.zip
wget https://files.pilookup.com/pi/200000001-300000000.zip
wget https://files.pilookup.com/pi/300000001-400000000.zip
wget https://files.pilookup.com/pi/400000001-500000000.zip
wget https://files.pilookup.com/pi/500000001-600000000.zip
wget https://files.pilookup.com/pi/600000001-700000000.zip
wget https://files.pilookup.com/pi/700000001-800000000.zip
wget https://files.pilookup.com/pi/800000001-900000000.zip
wget https://files.pilookup.com/pi/900000001-1000000000.zip