代碼 ：
#include
long a=10000,b,c=2800,d,e,f[2801],g;
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);
}
輸出結果是
3141592653589793238462643383279502884197169399375105820974
9445923078164062862089986280348253421170679821480865132823
0664709384460955058223172535940812848111745028410270193852
1105559644622948954930381964428810975665933446128475648233
7867831652712019091456485669234603486104543266482133936072
6024914127372458700660631558817488152092096282925409171536
4367892590360011330530548820466521384146951941511609433057
2703657595919530921861173819326117931051185480744623799627
4956735188575272489122793818301194912983367336244065664308
6021394946395224737190702179860943702770539217176293176752
3846748184676694051320005681271452635608277857713427577896
0917363717872146844090122495343014654958537105079227968925
8923542019956112129021960864034418159813629774771309960518
7072113499999983729780499510597317328160963185
