求π算法

作者在 2011-08-12 18:38:21 发布以下内容
#include <stdlib.h>  
#include <stdio.h>  
long a=10000,b,c=2800,d,e,f[2801],g,j,k,a1,a2,a3,a4;  
int main()  
{  
    for(;b-c;)  
    f[b++]=a/5;
    for(j=0;d=0,g=c*2;j++)
    {
      for(b=c;d+=f[b]*a,f[b]=d%--g,d/=g--,--b;d*=b);
      c-=14;
      if(j==0)
      {
          k=e+d/a;
          a1=k/1000;
          a2=k%1000/100;
          a3=k%1000%100/10;
          a4=k%10;
          printf("%d.%d%d%d",a1,a2,a3,a4);
      }
      else
      printf("%.4d",e+d/a);
      e=d%a;
    }
      return 0;
}
 
默认分类 | 阅读 415 次
文章评论,共0条
游客请输入验证码
文章分类
文章归档
最新评论