2

我已经开始了一个用于计算 pi 的第 n 位数字的 java 项目，并决定使用 BBP 算法。在我的输出中（在另一个类中），我遇到了一些奇怪的数学错误，我不知道它来自哪里。所以，我认为我没有正确地将算法放入代码中。

我从http://en.wikipedia.org/wiki/Bailey%E2%80%93Borwein%E2%80%93Plouffe_formula得到了算法

这是我的代码：

import java.lang.Math;
import java.lang.Math.*;
public class Pi
{
public static double getDigit(int n, int infinity)
{   int pow = 0;
    double[] sums = new double[4];
    int tot = 0;
    int result = 0;
    double fraction = 0;
    for(int x = 0; x < 4; x++)
    {
        for(int k = 0; k < n; k++)
        {
            tot = 8 * k + 1;
            if(x == 1)
                tot += 3;
            else if(x > 1)
                tot++;
            pow = n-k;
            result = modular_pow(16, pow, tot);
            sums[x] += (double)result / (double)tot;
        }
        for(int i = n + 1; i < infinity; i++)
        {
            tot = 8 * i + 1;
            if(x == 1)
                tot += 3;
            else if(x > 1)
                tot++;
            fraction = Math.pow(16.0, (double)pow);
            sums[x] += fraction / (double)tot;
        }
    }
    return 4 * sums[0] - 2 * sums[1] - sums[2] - sums[3];
}
public static int modular_pow(int base, int exponent, int modulus)
    {
    int result = 1;
    while(exponent > 0)
    {
        if (exponent % 2 == 1)
            result = (result * base) % modulus;
        exponent--;
        base = (base * base) % modulus;
    }
    return result;
}

提前致谢。

4

1 回答 1

1

首先，对删除旧帖子表示歉意，但是对于有意义地应用 BBP 算法的解释严重缺乏，所以我认为这对一些想要研究它的人来说可能仍然有用。

根据 Wikipedia 文章，您返回的结果需要去掉其整数部分（留下小数部分），然后乘以 16。这应该留下整数部分作为 pi 的第 n 个十六进制数字的表示。明天我会测试一下，看看是否有帮助。否则，很好的实施，易于理解和有效地完成。

于 2016-08-09T05:51:07.437 回答