求教一个计算圆周率的最优算法

2002-05-10 07:05:28 137 阅读 AI 写同款 · GEO 优化AI 选题数据监测一键分发 ›

《CLR via C#》CLR基础 CLR是公共语言运行库(Common Language Runtime)，又叫公共语言运行时。CLR和Java虚拟机一样也是一个运行时环境，是一个可由多种编程语言使用的运行环境，可由面向CLR的所有语言使用。它负责资源管理（内存分配和垃圾收集等），并保证应用和底层操作系统之间必要的分离。 .NET框架是一种以CLR为基础、支持多种语言（C#、F#、VB.NET、C++、Python等）的开发的软件... 阅读详情

写谢谢

...全文

Kindle电子书DRM破解全攻略：从下载到转换的完整流程（含最新DeDRM插件配置）本文详细介绍了如何合法移除个人购买的Kindle电子书的DRM保护，实现跨设备阅读与备份。核心内容包括使用最新DeDRM插件（如noDRM分支）配合Calibre进行配置，并针对亚马逊最新的KFX加密格式提供了密钥提取与实战操作指南。文章强调该方法仅适用于个人合理使用，旨在帮助读者真正掌控自己的数字资产。阅读详情

dengcainiao

数据结构与算法 33,026 成员 35,334 内容

+ 关注 + 加入社区

切换为时间正序

Tice 2002-05-20

打赏
举报

那要算多久？

qtf 2002-05-12

打赏
举报

这么厉害！！！
10亿？

xuying 2002-05-11

打赏
举报

上面的程序可以算任意位。

xuying 2002-05-11

打赏
举报

#include<stdio.h>
#include<stdlib.h>

#define Fcalloc calloc
#define Ffree free
#define Size_T unsigned long

long kf, ks;
long *mf, *ms;
long cnt, n, temp, nd;
long i;
long col, col1;
long loc, stor[21];

void shift(long *l1, long *l2, long lp, long lmod) {
long k;

k = ((*l2) > 0 ? (*l2) / lmod: -(-(*l2) / lmod) - 1);
*l2 -= k * lmod;
*l1 += k * lp;
}

void yprint(long m) {
if (cnt<n) {
if (++col == 11) {
col = 1;
if (++col1 == 6) {
col1 = 0;
printf("\n");
printf("%4ld",m%10);
}
else printf("%3ld",m%10);
}
else printf("%ld",m);
cnt++;
}
}

void xprint(long m) {
long ii, wk, wk1;

if (m < 8) {
for (ii = 1; ii <= loc; )
yprint(stor[(int)(ii++)]);
loc = 0;
}
else {
if (m > 9) {
wk = m / 10;
m %= 10;
for (wk1 = loc; wk1 >= 1; wk1--) {
wk += stor[(int)wk1];
stor[(int)wk1] = wk % 10;
wk /= 10;
}
}
}
stor[(int)(++loc)] = m;
}

void memerr(int errno) {
printf("\a\nOut of memory error #%d\n", errno);
if (2 == errno)
Ffree(mf);
_exit(2);
}

int main(int argc, char *argv[]) {
int i=0;
char *endp;

printf("\n");
printf("PI Caculater made by Godman\n");
printf("$Id: pi.c 1.14 1998/12/03 15:04:23 GUOXM Exp $\n");
printf("==============================================\n");

stor[i++] = 0;
if (argc < 2) {
puts("\aUsage: PI <number_of_digits>");
return(1);
}
n = strtol(argv[1], &endp, 10);
if (NULL == (mf = Fcalloc((Size_T)(n + 3L), (Size_T)sizeof(long))))
memerr(1);
if (NULL == (ms = Fcalloc((Size_T)(n + 3L), (Size_T)sizeof(long))))
memerr(2);
printf("\nApproximation of PI to %ld digits\n", (long)n);

cnt = 0;
kf = 25;
ks = 57121L;
mf[1] = 1L;

for (i = 2; i <= (int)n; i += 2) {
mf[i] = -16L;
mf[i+1] = 16L;
}

for (i = 1; i <= (int)n; i += 2) {
ms[i] = -4L;
ms[i+1] = 4L;
}

printf("\n 3.");
while (cnt < n) {
for (i = 0; ++i <= (int)n - (int)cnt; ) {
mf[i] *= 10L;
ms[i] *= 10L;
}
for (i =(int)(n - cnt + 1); --i >= 2; ) {
temp = 2 * i - 1;
shift(&mf[i - 1], &mf[i], temp - 2, temp * kf);
shift(&ms[i - 1], &ms[i], temp - 2, temp * ks);
}
nd = 0;
shift((long *)&nd, &mf[1], 1L, 5L);
shift((long *)&nd, &ms[1], 1L, 239L);
xprint(nd);
}
printf("\n\n======================\n");
printf("Calculations Completed!\n");
printf("You can fine this program and source code in ftp://166.111.18.56\n");
Ffree(ms);
Ffree(mf);
return(0);
}

kingfighter_han 2002-05-11

打赏
举报

NowCan 2002-05-11

打赏
举报

10亿位的都见过。

Tice 2002-05-10

打赏
举报

我又看过了，最多可以算到3335万位

Tice 2002-05-10

打赏
举报

我知道有个叫super pi的程序，可以算到100多万位

one_add_one 2002-05-10

打赏
举报

最优的我不知道，比较好的，我到时知道几个。

atlantis13579 2002-05-10

打赏
举报

肯定不是最优算法,但可以算.

long a=10000,b,c=2800,d,e,f[2801],g;
main(){for(;b-c;)f[b++]=a/5;
for(;d=0,g=c*2;c-=14,printf("%.4d",e+d/a),e=d%a)
for(b=c;d+=f[b]*a,f[b]=d%--g,d/=g--,--b;d*=b);}

求教一个计算圆周率的最优算法

相关推荐