一直对PI的计算程序感到好奇,到网上找了找相关文章。
照抄写了一个程序。
用了第三个式子。ARRAYSIZE代表位数。:)
#include <stdlib.h>
#include <stdio.h>
#include <string.h>
#define ARRAYSIZE 1000
int main() {
char x[ARRAYSIZE], z[ARRAYSIZE];
int a=1, b=3, d, Cnt=0, Run = 1;
memset(x, 0, sizeof(x));
memset(z, 0, sizeof(z));
x[ARRAYSIZE-2]= 2;
z[ARRAYSIZE-2]= 2;
while(Run && ++Cnt<ARRAYSIZE*100) {
//z*=a
d = 0;
for (int i=0; i<ARRAYSIZE; i++) {
int c = z[i]*a + d;
z[i] = c%10;
d = c/10;
}
// z/=b
d = 0;
for (int i=ARRAYSIZE-1; i>=0; i--) {
int c = z[i] + d*10;
d = c % b;
z[i] = c / b;
}
// x+=z
d = 0;
Run = 0;
for (int i=0; i<ARRAYSIZE; i++) {
int c = x[i] + z[i] + d;
x[i] = c%10;
d = c/10;
Run |= z[i];
}
a++;
b+=2;
}
printf(" Run count: %d/n", Cnt);
printf(" pi = ");
for (int i=ARRAYSIZE-1; i>=0; i--) {
printf("%1d", x[i]);
if (i == ARRAYSIZE-2) putchar('.');
}
putchar('/n');
}
还有算的Pi后一千位。
Run count: 3310
pi = 03.
14159265358979323846264338327950288419716939937510
58209749445923078164062862089986280348253421170679
82148086513282306647093844609550582231725359408128
48111745028410270193852110555964462294895493038196
44288109756659334461284756482337867831652712019091
45648566923460348610454326648213393607260249141273
72458700660631558817488152092096282925409171536436
78925903600113305305488204665213841469519415116094
33057270365759591953092186117381932611793105118548
07446237996274956735188575272489122793818301194912
98336733624406566430860213949463952247371907021798
60943702770539217176293176752384674818467669405132
00056812714526356082778577134275778960917363717872
14684409012249534301465495853710507922796892589235
42019956112129021960864034418159813629774771309960
51870721134999999837297804995105973173281609631859
50244594553469083026425223082533446850352619311881
71010003137838752886587533208381420617177669147303
59825349042875546873115956286388235378759375195778
185778053217122680661300192787661119590921638728
原文地址: http://www.cppfans.com/articles/basecalc/c_pi_10000.asp
照抄写了一个程序。
用了第三个式子。ARRAYSIZE代表位数。:)
#include <stdlib.h>
#include <stdio.h>
#include <string.h>
#define ARRAYSIZE 1000
int main() {
char x[ARRAYSIZE], z[ARRAYSIZE];
int a=1, b=3, d, Cnt=0, Run = 1;
memset(x, 0, sizeof(x));
memset(z, 0, sizeof(z));
x[ARRAYSIZE-2]= 2;
z[ARRAYSIZE-2]= 2;
while(Run && ++Cnt<ARRAYSIZE*100) {
//z*=a
d = 0;
for (int i=0; i<ARRAYSIZE; i++) {
int c = z[i]*a + d;
z[i] = c%10;
d = c/10;
}
// z/=b
d = 0;
for (int i=ARRAYSIZE-1; i>=0; i--) {
int c = z[i] + d*10;
d = c % b;
z[i] = c / b;
}
// x+=z
d = 0;
Run = 0;
for (int i=0; i<ARRAYSIZE; i++) {
int c = x[i] + z[i] + d;
x[i] = c%10;
d = c/10;
Run |= z[i];
}
a++;
b+=2;
}
printf(" Run count: %d/n", Cnt);
printf(" pi = ");
for (int i=ARRAYSIZE-1; i>=0; i--) {
printf("%1d", x[i]);
if (i == ARRAYSIZE-2) putchar('.');
}
putchar('/n');
}
还有算的Pi后一千位。
Run count: 3310
pi = 03.
14159265358979323846264338327950288419716939937510
58209749445923078164062862089986280348253421170679
82148086513282306647093844609550582231725359408128
48111745028410270193852110555964462294895493038196
44288109756659334461284756482337867831652712019091
45648566923460348610454326648213393607260249141273
72458700660631558817488152092096282925409171536436
78925903600113305305488204665213841469519415116094
33057270365759591953092186117381932611793105118548
07446237996274956735188575272489122793818301194912
98336733624406566430860213949463952247371907021798
60943702770539217176293176752384674818467669405132
00056812714526356082778577134275778960917363717872
14684409012249534301465495853710507922796892589235
42019956112129021960864034418159813629774771309960
51870721134999999837297804995105973173281609631859
50244594553469083026425223082533446850352619311881
71010003137838752886587533208381420617177669147303
59825349042875546873115956286388235378759375195778
185778053217122680661300192787661119590921638728
原文地址: http://www.cppfans.com/articles/basecalc/c_pi_10000.asp