SOFTIST 编程方法筆記 →目录

计算圆周率(算法)

圆周率π=3.1415926....，有时候需要用到小数的很多位。一般利用级数来求出圆周率π的近似値。本文纪录用下面的級数计算小数点后１０００位的π的算法。



１．先使用浮点数double型，理解算法。 
void CalculatePiWithDouble()
{
    double x = 2;
    double z = 2;
    int a=1, b=3;
    while ( z > 1e-15)
    {
        z = z * a / b;
        x += z;
        a ++;
        b += 2;
    }
    printf("π=%.13f\r\n", x);
}

２．使用数组代替浮点数double型，以增加小数位数，提高π的精度。 

void CalculateP() 
{
    #define ARRSIZE 1010
    #define DISPCNT 1000
    int i;
    char x[ARRSIZE];
    char z[ARRSIZE]; //x[0] x[1] . x[2] x[3] x[4] .... x[ARRSIZE-1]
    int a=1, b=3, c, d, Run=1, Cnt=0;

    memset(x, 0, ARRSIZE);
    memset(z, 0, ARRSIZE);
    x[1] = 2;
    z[1] = 2;
    while (Run && (++Cnt < 200000000))
    {
        // z*=a;
        d = 0;
        for (i = ARRSIZE - 1; i > 0; i --)
        {
            c = z[i] * a + d;
            z[i] = c % 10;
            d = c / 10;
        }

        // z/=b;
        d = 0;
        for(i = 0; i < ARRSIZE; i ++)
        {
            c = z[i] + d * 10;
            z[i] = c / b;
            d = c % b;
        }

        // x+=z;
        Run = 0;
        for(i = ARRSIZE - 1; i > 0; i --)
        {
            c = x[i] + z[i];
            x[i] = c % 10;
            x[i-1] += c / 10;
            Run |= z[i];
        }
        a ++;
        b +=2;
    }
    printf("次数 = %d 次\r\n", Cnt);
    printf("π=%d%d.\r\n", x[0], x[1]);
    for(i = 0; i < DISPCNT; i ++)
    {
        if (i && (i % 50 == 0))
            printf("\r\n");
        printf("%c", (char)(x[i+2] + '0'));
    }
    printf("\r\n");
}