Question

所以我想用任意精度逐位输出Pi。我在这里https://rosettacode.org/wiki/Pi#C.23

找到了一些代码

public IEnumerable<long> PiDigits(long b = 10)
    {
        BigInteger
            k = 1,
            l = 3,
            n = 3,
            q = 1,
            r = 0,
            t = 1
            ;

        // skip integer part
        var nr = b * (r - t * n);
        n = b * (3 * q + r) / t - b * n;
        q *= b;
        r = nr;

        for (; ; )
        {
            var tn = t * n;
            if (4 * q + r - t < tn)
            {
                yield return (long)n;
                nr = b * (r - tn);
                n = b * (3 * q + r) / t - b * n;
                q *= b;
            }
            else
            {
                t *= l;
                nr = (2 * q + r) * l;
                var nn = (q * (7 * k) + 2 + r * l) / t;
                q *= k;
                l += 2;
                ++k;
                n = nn;
            }
            r = nr;
        }
    }

它可以正常工作，但是没有任何解释。我想了解这个算法，但是看起来太复杂了。变量的名字什么也没给我。所以，请您告诉我我可以继续进行研究的算法是什么。谢谢！

Answer 1

这是从实现spigot algorithm转换而来的by Jeremy Gibbons的实现。

他的论文解释了该算法的内部工作原理。您应该在第5和6页上找到所需的信息。

用任意精度逐位计算Pi

1 个答案: