cpp2.cpp

vc编写 计算圆周率程序

#include "stdio.h"
#include "math.h"
#define N 100000  /*数组中每个元素储存一个五位数 */
#define M 201    /*数组含M个元素,共储存5*M位数*/
void fdiv(long a[],unsigned int z) /*超长整型数a[]除以整数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 fadd(long pi[],long a[],int k)
/*a[]除以k后,添上符号,加到pi[]中*/
{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)/*计算ARCTAN(1/t)*/
{int k=1,s,i;
 long b[M];
 a[0]=t*N;
 s=5*M/log10(t);
 do
  { fdiv(a,t*t);
    for(i=0;i<M;i++)
     b[i]=a[i];
   fdiv(b,k);
   fadd(pi,b,k);
   k+=2;
  }while(k<=s);
}
void main()
{int i,k;
 static long a[M],pi[M];
 fff(a,pi,2); /*计算ARCTAN(1/2)*/
 fff(a,pi,3); /*计算ARCTAN(1/3)*/
 for(i=0;i<M;i++)
   pi[i]*=4;
 for(i=1;i<M;i++)      /*整理进位,使每个元素的值不大于N*/
   {pi[i-1]+=pi[i]/N;
    pi[i]-=(pi[i]/N)*N;
   }
 printf("\n");
 for(i=0;i<M-1;i++)  /*控制输出格式,空格补0*/
   {if(pi[i]<N/10)
     for(k=log10((float)pi[i]);k<(int)log10(N)-1;k++)
       printf("0");
    printf("%ld",pi[i]);
   }
 }