Pi的计算程序

一直对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