c语言 圆周率代码

01:  #include <stdio.h>
02:  
03:  const int DIGITS = 2400;  // must: DIGITS % LEN == 0
04:  const int BASE   = 10000; // BASE == 10 ** LEN
05:  const int LEN    = 4;
06:  const int TIMES  = 14;
07:  
08:  int get_item(int pi[], int n)
09:  {
10:    int item = 0;
11:    for (int k = n - 1; k >= 0; k--) {
12:      item += pi[k] * BASE;
13:      pi[k] = item % (k * 2 + 1);
14:      item /= (k * 2 + 1);
15:      if (k > 0) item *= k;
16:    }
17:    return item;
18:  }
19:  
20:  int main()
21:  {
22:    int pi[DIGITS / LEN * TIMES], n = sizeof(pi) / sizeof(pi[0]);
23:    for (int i = 0; i < n; i++) pi[i] = BASE / 5; // BASE / 5 == 2.000
24:    for (int remainder = 0; n > 0; n -= TIMES) {
25:      int item = get_item(pi, n);
26:      printf("%.*d", LEN, remainder + item / BASE);
27:      remainder = item % BASE;
28:    }
29:  }

//环境：VC6.0,Console Application
//原理：π=2+1/3*(2+2/5*(2+3/7*(2+...
//特点：内嵌汇编提速并扩大了计算范围
//限制：位数ws原则上没有限制但因为本
//算法的时间正比于ws的平方所以将位数
//控制在二、三十万以内较好。本人曾用
//奔Ⅳ2.6GHz算20万位，耗时10分钟左右
#include<stdio.h>
#include<stdlib.h>
#include<malloc.h>
#include<time.h>
long a=100000L;
void main()
{ FILE *fp;
long t1,t2;
char filename[40];
unsigned long c,d,e,i,j,ws;
unsigned long *f,*bb;
printf("位数=?");
scanf("%ld",&ws);
if(ws<1)return;
c=(ws+4)/05*17;
bb=f=(long*)malloc(04*c);
if(f==NULL)abort();
printf("将Pi存为: ");//提示输入数据文件名
scanf("%s",filename);//若打入NUL,则不存盘
fp=fopen(filename,"w");
if(fp==NULL)abort();
t1=time(NULL);
for(i=0;i<c-1;i++,bb++)
*bb=a/5;*bb=a/5;//这里并不错
for(e=0;c;c-=17,bb-=17,e=d%a)
{ static long group;
d=0;i=c;
j=c+c-1;
_asm mov eDI,bb
loopi:
_asm mov eax,[eDI]
_asm mul Dword ptr a
_asm mov ecx,edx
_asm mov ebx,eax
_asm mov eax,d
_asm mul dword ptr i
_asm add eax,ebx
_asm adc edx,ecx
_asm div dword ptr j
_asm mov d,eax
_asm mov [eDI],edx
_asm sub eDI,04
_asm dec dword ptr j
_asm dec dword ptr j
_asm dec dword ptr i
_asm jnz loopi
e += d/a;//这是为了避免重复计算
if(e>=100000)//这是一种糟糕情况
abort();//希望它最好不要发生
if(++group%200==0) //逢千位附近
printf( "%05lu\t",e);//显示五位
fprintf(fp,"%05lu",e);//写入磁盘
}
t2=time(NULL);
printf("%ld秒\n",t2-t1);
free(f);
fclose(fp);
}

 
 


 
 
两个都很厉害，收藏下来