e60k.c

此程序可以计算自然对数的底到小数点后6万位。如果编译器对汉字有限制

/*这个程序是本人自编的*/ 
#include "stdio.h"
#include "time.h"
main()
{
   long a[15010],c,d,k,i,m,n;
   long t1,t2;
   FILE *fp;
   fp=fopen("6万位的e.txt","w");
   fprintf(fp,"计算自然对数的底e到6万位！\n"); 
   t1=time(NULL);
   a[0]=1L;
   d=10000L;
   for(i=1L;i<=15010L;i++)
   {
     a[i]=0L;  /*清零*/
   }
   for(n=15928L;n>=2;n--)
   {
     for(i=0L;i<=15009;i++)
     {
       c=a[i]%n;
       a[i]=a[i]/n;
       a[i+1]+=c*d;
     }
     a[0]+=1L;
   }
   fprintf(fp,"e=%d.\n", a[0]+1L);
   for(m=0L;m<=599L;m++)
   {
     k=25L*m;
     for(n=1L;n<=25L;n++)
     {
       fprintf(fp,"%.4d",a[k+n]);
     }
     fprintf(fp,"%10d\n",k*4L+100L);
   }
   t2=time(NULL);
   fprintf(fp,"%ds 计算完成！",t2-t1); 
   fclose(fp);
   printf("%ds  OK!!!\n",t2-t1);
   /*getchar();*/
}