请教：快速求得圆周率小数点后10000位的算法

[liqm]

请教：快速求得圆周率小数点后10000位的算法

[Arter]

拉马努金公式

[bravomyself]

#include <math.h>
#include <stdio.h>
void main(void)
{float s;
int b,x,n,c,i,j,d,l,a[10000];
printf("Please input the desired bit.\n");
scanf("%d",&x);
for(s=0,n=1;n<=10000;n++)
[
s=s+log(2*n+1)/n);
if(s/log(10)>x)break;
}
for(i=0;i<=x+1;i++)
a[i]=0;
for(c=1,j=n;j>=1;j--)
{
d=2*j-1;
for(i=0;i<=x;i++)
{
a[i]=c/d;
c=c(c%d)*10+a[i+1];
}
a[x+1]=c/d;
for(b=0,i=x+1;i>=0;i--)
{
a[i]=a[i]*j+b;b=a[i]/10;
a[i]%=10;
}
a[0]++;
c=a[0];
}
for(b=0,i=x+1;i>=0;i--)
{
a[i]=a[i}*2+b;
b=a[i]/10;
a[i]%=10;
}
printf("\npi=%d",a[0]);
for(l=10,i=1;i<=x;i++)
{
printf("%d",a[i]);l++;
if(l%10==0)
printf(" ");'
if(l%30==0)
printf("\n");
}