圆周率源程序1.txt

vc编写 计算圆周率程序

#include "stdio.h"
#include "math.h"                 
#define N 100000
#define M 1001
void div(long a[],unsigned int z)
{int i;
 unsigned long temp,p=0;
 for(i=0;i<M;i++)
   {temp=p*N+a[i];
    p=temp%z;
    a[i]=temp/z;
   }
}
void add(long pi[],long a[],int k)
{int i;
 for(i=0;i<M;i++)
  {if((k/2)%2==0)
     pi[i]+=a[i];
   else
     {pi[i]-=a[i];
      while(pi[i]<0)
       {pi[i]+=N;
        pi[i-1]--;
       }
     }
  }
}
void fff(long a[],long pi[],int t)
{int k=1,s,i;
 long b[M];
 a[0]=(long)t*N;
 switch(t)
  { case 8 : a[0]*=6; break;
    case 57: a[0]*=2; break;
  };
 s=5*M/log10(t);
 do
  {if(t!=239)
    div(a,t*t);
   else
    {div(a,t);
     div(a,t);
    }
   for(i=0;i<M;i++)
     b[i]=a[i];
   div(b,k);
   add(pi,b,k);
   k+=2;
  }while(k<=s);
}
void main()
{int i,k;
 static long a[M],pi[M];
 fff(a,pi,8);
 fff(a,pi,57);
 fff(a,pi,239);
 for(i=0;i<M;i++)
   pi[i]*=4;
 for(i=1;i<M;i++)
   {pi[i-1]+=pi[i]/N;
    pi[i]-=(pi[i]/N)*N;
   }
 printf("\n");
 for(i=0;i<M-1;i++)
   {if(pi[i]<N/10)
     for(k=log10((float)pi[i]);k<(int)log10(N)-1;k++)
       printf("0");
    printf("%ld",pi[i]);
   }
}