主题
Search

Pi 公式

下载 Mathematica Notebook下载 Wolfram Notebook

关于 pi 有许多公式，类型繁多。 其中包括级数、乘积、几何构造、极限、特殊值和 pi 迭代等。

pi 与圆和球体的性质密切相关。 对于半径为 r r 的圆，其周长和面积由下式给出

C=2pir
(1)
A=pir^2.
(2)

类似地，对于半径为 r r 的球体，其表面积和体积分别为

S=4pir^2
(3)
V=4/3pir^3.
(4)

关于 pi 的精确公式，以 反正切 的 单位分数 表示，是 马钦公式

1/4pi=4tan^(-1)(1/5)-tan^(-1)(1/(239)).
(5)

还有其他三个 类马钦公式，以及数千个其他类似的具有更多项的公式。

GregorySeries

格雷戈里和莱布尼茨发现

pi/4=sum_(k=1)^(infty)((-1)^(k+1))/(2k-1)
(6)
=1-1/3+1/5-...
(7)

（Wells 1986, p. 50），这被称为 格雷戈里级数，可以通过将 x=1 代入 莱布尼茨级数 tan^(-1)x 得到。 格雷戈里级数 中此级数的第 n 项之后的误差大于 (2n)^(-1)，因此该和收敛非常缓慢，以至于 300 项不足以正确计算 pi 到两位小数！ 然而，它可以转换为

pi=sum_(k=1)^infty(3^k-1)/(4^k)zeta(k+1),
(8)

其中 zeta(z) 是 黎曼 zeta 函数（Vardi 1991, pp. 157-158; Flajolet and Vardi 1996），因此 k 项之后的误差为 approx (3/4)^k。

亚伯拉罕·夏普 (ca. 1717) 的一个无穷和级数由下式给出

pi=sum_(k=0)^infty(2(-1)^k3^(1/2-k))/(2k+1)
(9)

（Smith 1953, p. 311）。 其中出现 pi 的其他简单级数为

1/4pisqrt(2)=sum_(k=1)^(infty)[((-1)^(k+1))/(4k-1)+((-1)^(k+1))/(4k-3)]
(10)
=1+1/3-1/5-1/7+1/9+1/(11)-...
(11)
1/4(pi-3)=sum_(k=1)^(infty)((-1)^(k+1))/(2k(2k+1)(2k+2))
(12)
=1/(2·3·4)-1/(4·5·6)+1/(6·7·8)-...
(13)
1/6pi^2=sum_(k=1)^(infty)1/(k^2)
(14)
=1+1/4+1/9+1/(16)+1/(25)+...
(15)
1/8pi^2=sum_(k=1)^(infty)1/((2k-1)^2)
(16)
=1+1/(3^2)+1/(5^2)+1/(7^2)+...
(17)

（Wells 1986, p. 53）。

1666 年，牛顿使用几何构造推导出公式

pi=3/4sqrt(3)+24int_0^(1/4)sqrt(x-x^2)dx
(18)
=(3sqrt(3))/4+24(1/(12)-1/(5·2^5)-1/(28·2^7)-1/(72·2^9)-...),
(19)

他用此公式计算 pi（Wells 1986, p. 50; Borwein et al. 1989; Borwein and Bailey 2003, pp. 105-106）。 系数可以从积分中找到

I(x)=intsqrt(x-x^2)dx
(20)
=1/4(2x-1)sqrt(x-x^2)-1/8sin^(-1)(1-2x)
(21)

通过取 I(x)-I(0) 在 0 附近的级数展开，得到

I(x)=2/3x^(3/2)-1/5x^(5/2)-1/(28)x^(7/2)-1/(72)x^(9/2)-5/(704)x^(11/2)+...
(22)

（OEIS A054387 和 A054388）。 使用欧拉的 收敛加速 变换得到

pi/2=1/2sum_(n=0)^(infty)((n!)^22^(n+1))/((2n+1)!)=sum_(n=0)^(infty)(n!)/((2n+1)!!)
(23)
=1+1/3+(1·2)/(3·5)+(1·2·3)/(3·5·7)+...
(24)
=1+1/3(1+2/5(1+3/7(1+4/9(1+...))))
(25)

（Beeler et al. 1972, Item 120）。

这对应于将 x=1/sqrt(2) 代入 幂级数 中的 超几何函数 _2F_1(a,b;c;x)，

(sin^(-1)x)/(sqrt(1-x^2))=sum_(i=0)^infty((2x)^(2i+1)i!^2)/(2(2i+1)!)=_2F_1(1,1;3/2;x^2)x.
(26)

尽管进行了收敛加速，级数 (◇) 的收敛速度仍然只有每项 1 位。 以 平方根 为代价，Gosper 指出 x=1/2 给出每项 2 位，

1/9sqrt(3)pi=1/2sum_(i=0)^infty((i!)^2)/((2i+1)!),
(27)

而 x=sin(pi/10) 给出几乎每项 3.39 位，

pi/(5sqrt(phi+2))=1/2sum_(i=0)^infty((i!)^2)/(phi^(2i+1)(2i+1)!),
(28)

其中 phi 是 黄金比例。 Gosper 还得到了

pi=3+1/(60)(8+(2·3)/(7·8·3)(13+(3·5)/(10·11·3)(18+(4·7)/(13·14·3)(23+...)))).
(29)

各种极限也收敛于 pi，一个简单的例子是

lim_(n->infty)ncos((pi(n-2))/(2n))=pi.
(30)

更有趣的例子由下式给出

lim_(n->infty)((2(-1)^(n+1)(2n)!)/(2^(2n)B_(2n)))^(1/(2n))=pi,
(31)

和

lim_(n->infty)((2(-1)^(n+1)(2n)!)/(2^(2n)B_(2n)(1-2^(-n))(1-3^(-n))(1-5^(-n))(1-7^(-n))))^(1/(2n))=pi,
(32)

其中 B_n 是 伯努利数（Plouffe 2022）。 这些公式可以用作 pi 数字的 数字提取算法。

Rabinowitz 和 Wagon (1995; Borwein and Bailey 2003, pp. 141-142) 给出了 pi 的 spigot 算法。

Bailey et al. (Bailey et al. 1997, Adamchik and Wagon 1997) 发现了一个闭式表达式，给出了另一种 数字提取算法，该算法生成以 16 为基数的 pi（或 pi^2）的数字，

pi=sum_(n=0)^infty(4/(8n+1)-2/(8n+4)-1/(8n+5)-1/(8n+6))(1/(16))^n.
(33)

该公式被称为 BBP 公式，是使用 PSLQ 算法 (Ferguson et al. 1999) 发现的，等价于

pi=int_0^1(16y-16)/(y^4-2y^3+4y-4)dy.
(34)

存在一系列 BBP 类型公式，用于表示 pi，以 (-1)^k 的幂表示，其中前几个独立公式为

pi=4sum_(k=0)^(infty)((-1)^k)/(2k+1)
(35)
=3sum_(k=0)^(infty)(-1)^k[1/(6k+1)+1/(6k+5)]
(36)
=4sum_(k=0)^(infty)(-1)^k[1/(10k+1)-1/(10k+3)+1/(10k+5)-1/(10k+7)+1/(10k+9)]
(37)
=sum_(k=0)^(infty)(-1)^k[3/(14k+1)-3/(14k+3)+3/(14k+5)+4/(14k+7)+4/(14k+9)-4/(14k+11)+4/(14k+13)]
(38)
=sum_(k=0)^(infty)(-1)^k[2/(18k+1)+3/(18k+3)+2/(18k+5)-2/(18k+7)-2/(18k+11)+2/(18k+13)+3/(18k+15)+2/(18k+17)]
(39)
=sum_(k=0)^(infty)(-1)^k[3/(22k+1)-3/(22k+3)+3/(22k+5)-3/(22k+7)+3/(22k+9)+8/(22k+11)+3/(22k+13)-3/(22k+15)+3/(22k+17)-3/(22k+19)+1/(22k+21)].
(40)

类似地，存在一系列 BBP 类型公式，用于表示 pi，以 2^k 的幂表示，其中前几个独立公式为

pi=sum_(k=0)^(infty)1/(16^k)[4/(8k+1)-2/(8k+4)-1/(8k+5)-1/(8k+6)]
(41)
=1/2sum_(k=0)^(infty)1/(16^k)[8/(8k+2)+4/(8k+3)+4/(8k+4)-1/(8k+7)]
(42)
=1/(16)sum_(k=0)^(infty)1/(256^k)[(64)/(16k+1)-(32)/(16k+4)-(16)/(16k+5)-(16)/(16k+6)+4/(16k+9)-2/(16k+12)-1/(16k+13)-1/(16k+14)]
(43)
=1/(32)sum_(k=0)^(infty)1/(256^k)[(128)/(16k+2)+(64)/(16k+3)+(64)/(16k+4)-(16)/(16k+7)+8/(16k+10)+4/(16k+11)+4/(16k+12)-1/(16k+15)]
(44)
=1/(32)sum_(k=0)^(infty)1/(4096^k)[(256)/(24k+2)+(192)/(24k+3)-(256)/(24k+4)-(96)/(24k+6)-(96)/(24k+8)+(16)/(24k+10)-4/(24k+12)-3/(24k+15)-6/(24k+16)-2/(24k+18)-1/(24k+20)]
(45)
=1/(64)sum_(k=0)^(infty)1/(4096^k)[(256)/(24k+1)+(256)/(24k+2)-(384)/(24k+3)-(256)/(24k+4)-(64)/(24k+5)+(96)/(24k+8)+(64)/(24k+9)+(16)/(24k+10)+8/(24k+12)-4/(24k+13)+6/(24k+15)+6/(24k+16)+1/(24k+17)+1/(24k+18)-1/(24k+20)-1/(24k+21)]
(46)
=1/(96)sum_(k=0)^(infty)1/(4096^k)[(256)/(24k+2)+(64)/(24k+3)+(128)/(24k+5)+(352)/(24k+6)+(64)/(24k+7)+(288)/(24k+8)+(128)/(24k+9)+(80)/(24k+10)+(20)/(24k+12)-(16)/(24k+14)-1/(24k+15)+6/(24k+16)-2/(24k+17)-1/(24k+19)+1/(24k+20)-2/(24k+21)]
(47)
=1/(96)sum_(k=0)^(infty)1/(4096^k)[(256)/(24k+1)+(320)/(24k+3)+(256)/(24k+4)-(192)/(24k+5)-(224)/(24k+6)-(64)/(24k+7)-(192)/(24k+8)-(64)/(24k+9)-(64)/(24k+10)-(28)/(24k+12)-4/(24k+13)-5/(24k+15)+3/(24k+17)+1/(24k+18)+1/(24k+19)+1/(24k+21)-1/(24k+22)]
(48)
=1/(96)sum_(k=0)^(infty)1/(4096^k)[(512)/(24k+1)-(256)/(24k+2)+(64)/(24k+3)-(512)/(24k+4)-(32)/(24k+6)+(64)/(24k+7)+(96)/(24k+8)+(64)/(24k+9)+(48)/(24k+10)-(12)/(24k+12)-8/(24k+13)-(16)/(24k+14)-1/(24k+15)-6/(24k+16)-2/(24k+18)-1/(24k+19)-1/(24k+20)-1/(24k+21)]
(49)
=1/(4096)sum_(k=0)^(infty)1/(65536^k)[(16384)/(32k+1)-(8192)/(32k+4)-(4096)/(32k+5)-(4096)/(32k+6)+(1024)/(32k+9)-(512)/(32k+12)-(256)/(32k+13)-(256)/(32k+14)+(64)/(32k+17)-(32)/(32k+20)-(16)/(32k+21)-(16)/(32k+22)+4/(32k+25)-2/(32k+28)-1/(32k+29)-1/(32k+30)]
(50)
=1/(4096)sum_(k=0)^(infty)1/(65536^k)[(32768)/(32k+2)+(16384)/(32k+3)+(16384)/(32k+4)-(4096)/(32k+7)+(2048)/(32k+10)+(1024)/(32k+11)+(1024)/(32k+12)-(256)/(32k+15)+(128)/(32k+18)+(64)/(32k+19)+(64)/(32k+20)-(16)/(32k+23)+8/(32k+26)+4/(32k+27)+4/(32k+28)-1/(32k+31)].
(51)

F. Bellard 发现了快速收敛的 BBP 类型公式

pi=1/(2^6)sum_(n=0)^infty((-1)^n)/(2^(10n))(-(2^5)/(4n+1)-1/(4n+3)+(2^8)/(10n+1)-(2^6)/(10n+3)-(2^2)/(10n+5)-(2^2)/(10n+7)+1/(10n+9)).
(52)

一个相关的积分是

pi=(22)/7-int_0^1(x^4(1-x)^4)/(1+x^2)dx
(53)

（Dalzell 1944, 1971; Le Lionnais 1983, p. 22; Borwein, Bailey, and Girgensohn 2004, p. 3; Boros and Moll 2004, p. 125; Lucas 2005; Borwein et al. 2007, p. 14）。 该积分在 20 世纪 60 年代中期被 K. Mahler 知道，并出现在 1960 年 11 月悉尼大学的一次考试中（Borwein, Bailey, and Girgensohn, p. 3）。 Beukers (2000) 和 Boros 和 Moll (2004, p. 126) 指出，尚不清楚是否存在合理的有理多项式选择，其在 0 到 1 之间的积分产生 pi-333/106，其中 333/106 是下一个收敛值。 然而，第四个收敛值存在一个积分，即

pi=(355)/(113)-1/(3164)int_0^1(x^8(1-x)^8(25+816x^2))/(1+x^2)dx.
(54)

（Lucas 2005; Bailey et al. 2007, p. 219）。 事实上，Lucas (2005) 给出了其他一些这样的积分。

Backhouse (1995) 使用了恒等式

I_(m,n)=int_0^1(x^m(1-x)^n)/(1+x^2)dx
(55)
=2^(-(m+n+1))sqrt(pi)Gamma(m+1)Gamma(n+1)×_3F_2(1,(m+1)/2,(m+2)/2;(m+n+2)/2,(m+n+3)/2;-1)
(56)
=a+bpi+cln2
(57)

对于正整数 m 和 n，其中 a、b 和 c 是有理常数，以生成许多 pi 的公式。 特别是，如果 2m-n=0 (mod 4)，则 c=0 (Lucas 2005)。

Ferguson 随后发现了一个类似的公式，从而产生了此类公式的二维格子，这些公式可以由以下两个公式生成

pi=sum_(k=0)^infty((4+8r)/(8k+1)-(8r)/(8k+2)-(4r)/(8k+3)-(2+8r)/(8k+4)-(1+2r)/(8k+5)-(1+2r)/(8k+6)+r/(8k+7))(1/(16))^k
(58)

对于 r 的任何复数值（Adamchik 和 Wagon），BBP 公式 作为特殊情况 r=0。

PiFormulasWagonIdentity

Wagon 提出的一个更通用的恒等式由下式给出

pi+4tan^(-1)z+2ln((1-2z-z^2)/(z^2+1))=sum_(k=0)^infty1/(16^k)[(4(z+1)^(8k+1))/(8k+1)-(2(z+1)^(8k+4))/(8k+4)-((z+1)^(8k+5))/(8k+5)-((z+1)^(8k+6))/(8k+6)]
(59)

（Borwein and Bailey 2003, p. 141），它在 复平面 的一个区域内成立，该区域排除了关于 实轴 对称放置的两个三角形部分，如上图所示。

一个可能更奇怪的通用恒等式类别由下式给出

pi=4sum_(j=1)^n((-1)^(j+1))/(2j-1)+((-1)^n(2n-1)!)/4sum_(k=0)^infty1/(16^k)[8/((8k+1)_(2n))-4/((8k+3)_(2n))-4/((8k+4)_(2n))-2/((8k+5)_(2n))+1/((8k+7)_(2n))+1/((8k+8)_(2n))]
(60)

对于任何正整数 n 成立，其中 (x)_n 是 波赫哈默尔符号（B. Cloitre，私人通讯，2005 年 1 月 23 日）。 更令人惊讶的是，对于 自然对数 2，存在一个非常类似的公式。

在发现以 16 为基数的数字 BBP 公式 和相关公式之后，研究了其他基数中的类似公式。 Borwein、Bailey 和 Girgensohn (2004) 最近表明，pi 没有非二元的马钦型 BBP 反正切公式，但这并不排除其他基数中 数字提取算法 的完全不同的方案。

S. Plouffe 设计了一种算法，用于计算任何基数中 n 位 数字 的 pi，步数为 O(n^3(logn)^3)。

Castellanos (1988ab, pp. 86-88) 给出了大量额外的恒等式，这些恒等式归功于拉马努金、卡塔兰和牛顿，包括几个涉及 斐波那契数 的和。 拉马努金发现

sum_(k=0)^infty((-1)^k(4k+1)[(2k-1)!!]^3)/([(2k)!!]^3)=sum_(k=0)^infty((-1)^k(4k+1)[Gamma(k+1/2)]^3)/(pi^(3/2)[Gamma(k+1)]^3)=2/pi
(61)

（Hardy 1923, 1924, 1999, p. 7）。

Plouffe (2006) 发现了美丽的公式

pi=72sum_(n=1)^infty1/(n(e^(npi)-1))-96sum_(n=1)^infty1/(n(e^(2npi)-1)) +24sum_(n=1)^infty1/(n(e^(4npi)-1)).
(62)
PiBlatnerProduct

欧拉提出的一个有趣的 无穷乘积 公式，它将 pi 与第 n 个 素数 p_n 联系起来，即

pi=2/(product_(n=1)^(infty)[1+(sin(1/2pip_n))/(p_n)])
(63)
=2/(product_(n=2)^(infty)[1+((-1)^((p_n-1)/2))/(p_n)])
(64)

（Blatner 1997, p. 119），如上图所示，它是乘积项数的函数。

类似于阿基米德的方法可用于估计 pi，方法是从一个 n 边形开始，然后关联后续 2n 边形的 面积。 令 beta 为一个 多边形 的线段中心的 角，

beta=1/4(n-3)pi,
(65)

那么

pi=(2sin(2beta))/((n-3)product_(k=0)^(infty)cos(2^(-k)beta))
(66)

（Beckmann 1989, pp. 92-94）。

Vieta (1593) 是第一个通过在上式中取 n=4 给出 pi 的精确表达式的人，给出

cosbeta=sinbeta=1/(sqrt(2))=1/2sqrt(2),
(67)

这导致 嵌套根式 的 无穷乘积，

2/pi=sqrt(1/2)sqrt(1/2+1/2sqrt(1/2))sqrt(1/2+1/2sqrt(1/2+1/2sqrt(1/2)))...
(68)

（Wells 1986, p. 50; Beckmann 1989, p. 95）。 然而，直到 1892 年 Rudio 才严格证明该表达式收敛。

一个相关公式由下式给出

pi=lim_(n->infty)2^(n+1)sqrt(2-sqrt(2+sqrt(2+sqrt(2+...+sqrt(2))))_()_(n)),
(69)

可以写成

pi=lim_(n->infty)2^(n+1)pi_n,
(70)

其中 pi_n 是使用迭代定义的

pi_n=sqrt((1/2pi_(n-1))^2+[1-sqrt(1-(1/2pi_(n-1))^2)]^2)
(71)

其中 pi_0=sqrt(2) (J. Munkhammar，私人通讯，2000 年 4 月 27 日)。 公式

pi=2lim_(m->infty)sum_(n=1)^msqrt([sqrt(1-((n-1)/m)^2)-sqrt(1-(n/m)^2)]^2+1/(m^2))
(72)

也密切相关。

一个漂亮的 pi 公式由下式给出

pi=(product_(n=1)^(infty)(1+1/(4n^2-1)))/(sum_(n=1)^(infty)1/(4n^2-1)),
(73)

其中分子是 pi/2 的 沃利斯公式 的一种形式，分母是总和为 1/2 的 伸缩和，因为

1/(4n^2-1)=1/2(1/(2n-1)-1/(2n+1))
(74)

（Sondow 1997）。

沃利斯公式 的一个特例给出

pi/2=product_(n=1)^infty[((2n)^2)/((2n-1)(2n+1))]=(2·2)/(1·3)(4·4)/(3·5)(6·6)/(5·7)...
(75)

（Wells 1986, p. 50）。 该公式也可以写成

lim_(n->infty)(2^(4n))/(n(2n; n)^2)=pilim_(n->infty)(n[Gamma(n)]^2)/([Gamma(1/2+n)]^2)=pi,
(76)

其中 (n; k) 表示 二项式系数，Gamma(x) 是 伽玛函数（Knopp 1990）。 欧拉得到

pi=sqrt(6(1+1/(2^2)+1/(3^2)+1/(4^2)+...)),
(77)

这来自 黎曼 zeta 函数 zeta(2)=pi^2/6 的特殊值。 对于所有 正整数 n，类似的 公式 来自 zeta(2n)。

拉马努金提出的一个无穷和是

1/pi=sum_(n=0)^infty(2n; n)^3(42n+5)/(2^(12n+4))
(78)

（Borwein et al. 1989; Borwein and Bailey 2003, p. 109; Bailey et al. 2007, p. 44）。 拉马努金 (1913-14) 给出了进一步的和，

4/pi=sum_(n=0)^infty((-1)^n(1123+21460n)(2n-1)!!(4n-1)!!)/(882^(2n+1)32^n(n!)^3)
(79)

和

1/pi=sqrt(8)sum_(n=0)^(infty)((1103+26390n)(2n-1)!!(4n-1)!!)/(99^(4n+2)32^n(n!)^3)
(80)
=(sqrt(8))/(9801)sum_(n=0)^(infty)((4n)!(1103+26390n))/((n!)^4396^(4n))
(81)

（Beeler et al. 1972, Item 139; Borwein et al. 1989; Borwein and Bailey 2003, p. 108; Bailey et al. 2007, p. 44）。 方程 (81) 源自 58 阶的模恒等式，尽管 Borwein 和 Borwein (1987) 之前没有提出第一个推导。 上述级数都给出

pi approx (9801)/(2206sqrt(2))=3.14159273001...
(82)

（Wells 1986, p. 54）作为第一个近似值，并且每项分别提供约 6 位和 8 位小数。 存在这样的级数是因为各种模不变量的合理性。

级数的一般形式为

sum_(n=0)^infty[a(t)+nb(t)]((6n)!)/((3n)!(n!)^3)1/([j(t)]^n)=(sqrt(-j(t)))/pi,
(83)

其中 t 是 二元二次型判别式，j(t) 是 j 函数，

b(t)=sqrt(t[1728-j(t)])
(84)
a(t)=(b(t))/6{1-(E_4(t))/(E_6(t))[E_2(t)-6/(pisqrt(t))]},
(85)

而 E_i 是 艾森斯坦级数。 类数 p 场涉及常数 A=a(t)、B=b(t) 和 C=c(t) 的 p 次 代数整数。 在所有仅由整数项组成的级数中，在最短时间内给出最多数字的级数对应于最大的 类数 1 判别式 d=-163，由 Chudnovsky 兄弟 (1987) 提出。 此处出现的 163 与 e^(pisqrt(163))（拉马努金常数）非常接近 整数 的事实中出现的 163 相同。 类似地，因子 640320^3 来自 j 函数 恒等式 j(1/2(1+isqrt(163)))。 该级数由下式给出

1/pi=12sum_(n=0)^(infty)((-1)^n(6n)!(13591409+545140134n))/((n!)^3(3n)!(640320^3)^(n+1/2))
(86)
=(163·8·27·7·11·19·127)/(640320^(3/2))sum_(n=0)^(infty)((13591409)/(163·2·9·7·11·19·127)+n)((6n)!)/((3n)!(n!)^3)((-1)^n)/(640320^(3n))
(87)

（Borwein and Borwein 1993; Beck and Trott; Bailey et al. 2007, p. 44）。 此级数每项精确给出 14 位数字。 Chudnovsky 兄弟 (1987) 给出了另一个形式的相同方程，Wolfram 语言 使用该方程来计算 pi（Vardi 1991; Wolfram Research），

pi=(426880sqrt(10005))/(A[_3F_2(1/6,1/2,5/6;1,1;B)-C_3F_2(7/6,3/2,(11)/6;2,2;B)]),
(88)

其中

A=13591409
(89)
B=-1/(151931373056000)
(90)
C=(30285563)/(1651969144908540723200).
(91)

类数 2 的最佳公式（最大判别式 -427）是

1/pi=12sum_(n=0)^infty((-1)^n(6n)!(A+Bn))/((n!)^3(3n)!C^(n+1/2)),
(92)

其中

A=212175710912sqrt(61)+1657145277365
(93)
B=13773980892672sqrt(61)+107578229802750
(94)
C=[5280(236674+30303sqrt(61))]^3
(95)

（Borwein and Borwein 1993）。 此级数每增加一项，就会增加约 25 位数字。 类数 3 的最快收敛级数对应于 d=-907，每项给出 37-38 位数字。 类数 4 的最快收敛级数对应于 d=-1555，即

(sqrt(-C^3))/pi=sum_(n=0)^infty((6n)!)/((3n)!(n!)^3)(A+nB)/(C^(3n)),
(96)

其中

A=63365028312971999585426220+28337702140800842046825600sqrt(5)+384sqrt(5)(10891728551171178200467436212395209160385656017+4870929086578810225077338534541688721351255040sqrt(5))^(1/2)
(97)
B=7849910453496627210289749000+3510586678260932028965606400sqrt(5)+2515968sqrt(3110)(6260208323789001636993322654444020882161+2799650273060444296577206890718825190235sqrt(5))^(1/2)
(98)
C=-214772995063512240-96049403338648032sqrt(5)-1296sqrt(5)(10985234579463550323713318473+4912746253692362754607395912sqrt(5))^(1/2).
(99)

每项给出 50 位数字。 Borwein 和 Borwein (1993) 开发了一种通用 算法，用于为任意 类数 生成此类级数。

Berndt (1994, pp. 352-354) 给出了拉马努金的第二个和第三个笔记本中发现的 1/pi 级数的完整列表，

4/pi=sum_(n=0)^(infty)((6n+1)(1/2)_n^3)/(4^n(n!)^3)
(100)
(16)/pi=sum_(n=0)^(infty)((42n+5)(1/2)_n^3)/(64^n(n!)^3)
(101)
(32)/pi=sum_(n=0)^(infty)((42sqrt(5)n+5sqrt(5)+30n-1)(1/2)_n^3)/(64^n(n!)^3)((sqrt(5)-1)/2)^(8n)
(102)
(27)/(4pi)=sum_(n=0)^(infty)((15n+2)(1/2)_n(1/3)_n(2/3)_n)/((n!)^3)(2/(27))^n
(103)
(15sqrt(3))/(2pi)=sum_(n=0)^(infty)((33n+4)(1/2)_n(1/3)_n(2/3)_n)/((n!)^3)(4/(125))^n
(104)
(5sqrt(5))/(2pisqrt(3))=sum_(n=0)^(infty)((11n+1)(1/2)_n(1/6)_n(5/6)_n)/((n!)^3)(4/(125))^n
(105)
(85sqrt(85))/(18pisqrt(3))=sum_(n=0)^(infty)((133n+8)(1/2)_n(1/6)_n(5/6)_n)/((n!)^3)(4/(85))^(3n)
(106)
4/pi=sum_(n=0)^(infty)((-1)^n(20n+3)(1/2)_n(1/4)_n(3/4)_n)/((n!)^32^(2n+1))
(107)
4/(pisqrt(3))=sum_(n=0)^(infty)((-1)^n(28n+3)(1/2)_n(1/4)_n(3/4)_n)/((n!)^33^n4^(2n+1))
(108)
4/pi=sum_(n=0)^(infty)((-1)^n(260n+23)(1/2)_n(1/4)_n(3/4)_n)/((n!)^3(18)^(2n+1))
(109)
4/(pisqrt(5))=sum_(n=0)^(infty)((-1)^n(644n+41)(1/2)_n(1/4)_n(3/4)_n)/((n!)^35^n(72)^(2n+1))
(110)
4/pi=sum_(n=0)^(infty)((-1)^n(21460n+1123)(1/2)_n(1/4)_n(3/4)_n)/((n!)^3(882)^(2n+1))
(111)
(2sqrt(3))/pi=sum_(n=0)^(infty)((8n+1)(1/2)_n(1/4)_n(3/4)_n)/((n!)^39^n)
(112)
1/(2pisqrt(2))=sum_(n=0)^(infty)((10n+1)(1/2)_n(1/4)_n(3/4)_n)/((n!)^39^(2n+1))
(113)
1/(3pisqrt(3))=sum_(n=0)^(infty)((40n+3)(1/2)_n(1/4)_n(3/4)_n)/((n!)^3(49)^(2n+1))
(114)
2/(pisqrt(11))=sum_(n=0)^(infty)((280n+19)(1/2)_n(1/4)_n(3/4)_n)/((n!)^3(99)^(2n+1))
(115)
1/(2pisqrt(2))=sum_(n=0)^(infty)((26390n+1103)(1/2)_n(1/4)_n(3/4)_n)/((n!)^3(99)^(4n+2)).
(116)

Borwein 和 Borwein (1987a, pp. 177-187) 首先证明了这些方程。 Borwein 和 Borwein (1987b, 1988, 1993) 证明了其他类型的方程，而 Chudnovsky 和 Chudnovsky (1987) 发现了其他超越常数的类似方程（Bailey et al. 2007, pp. 44-45）。

已知独立方程的完整列表如下所示

4/pi=sum_(n=0)^(infty)((6n+1)(1/2)_n^3)/(4^n(n!)^3)
(117)
(16)/pi=sum_(n=0)^(infty)((42n+5)(1/2)_n^3)/(64^n(n!)^3)
(118)
(32)/pi=sum_(n=0)^(infty)((42sqrt(5)n+5sqrt(5)+30n-1)(1/2)_n^3)/(64^n(n!)^3)((sqrt(5)-1)/2)^(8n)
(119)
(5^(1/4))/pi=sum_(n=0)^(infty)((540sqrt(5)n-1200n-525+235sqrt(5))(1/2)_n^3(sqrt(5)-2)^(8n))/((n!)^3)
(120)
(12^(1/4))/pi=sum_(n=0)^(infty)((24sqrt(3)n-36n-15+9sqrt(3))(1/2)_n^3(2-sqrt(3))^(4n))/((n!)^3)
(121)

对于 m=1，符号非交替，

2/pi=sum_(n=0)^(infty)((-1)^n(1/2)_n^3(12sqrt(2)n-12n-5+4sqrt(2))(sqrt(2)-1)^(4n))/((n!)^3)
(122)
2/pi=sum_(n=0)^(infty)((-1)^n(1/2)_n^3(60n-24sqrt(5)n+23-10sqrt(5))(sqrt(5)-2)^(4n))/((n!)^3)
(123)
2/pi=sum_(n=0)^(infty)((-1)^n(1/2)_n^3(420n-168sqrt(6)n+177-72sqrt(6)))/((n!)^3)
(124)
(2sqrt(2))/pi=sum_(n=0)^(infty)((-1)^n(1/2)_n^3(2sqrt(2))^(2n))/((n!)^3)
(125)

对于 m=1，符号交替，

(128)/(pi^2)=sum_(n=0)^(infty)((-1)^n(1/2)_n^5(820n^2+180n+13))/(32^(2n)(n!)^5)
(126)
(32)/(pi^2)=sum_(n=0)^(infty)((-1)^n(1/2)_n^5(20n^2+8n+1))/(2^(2n)(n!)^5)
(127)

对于 m=2 (Guillera 2002, 2003, 2006)，

(32)/(pi^3)=sum_(n=0)^infty((1/2)_n^7(168n^3+76n^2+14n+1))/(32^(2n)(n!)^5)
(128)

对于 m=3 (Guillera 2002, 2003, 2006)，并且对于 m>3 没有其他已知的（Bailey et al. 2007, pp. 45-48）。

Bellard 给出了奇异公式

pi=1/(740025)[sum_(n=1)^infty(3P(n))/((7n; 2n)2^(n-1))-20379280],
(129)

其中

P(n)=-885673181n^5+3125347237n^4-2942969225n^3+1031962795n^2-196882274n+10996648.
(130)

Gasper 引用了结果

pi=(16)/3[lim_(x->infty)x_1F_2(1/2;2,3;-x^2)]^(-1),
(131)

其中 _1F_2 是 广义超几何函数，并将其转换为

pi=lim_(x->infty)4x_1F_2(1/2;3/2,3/2;-x^2).
(132)

Gosper 提出的一个引人入胜的结果由下式给出

lim_(n->infty)product_(i=n)^(2n)pi/(2tan^(-1)i)=4^(1/pi)=1.554682275....
(133)

pi 满足 不等式

(1+1/pi)^(pi+1) approx 3.14097<pi.
(134)

D. Terr（私人通讯）注意到一个有趣的恒等式

(3,1,4)=(1,5,9)+(2,6,5) (mod 10)
(135)

涉及 pi 的前 9 位数字。

另请参阅

BBP 公式、数字提取算法、Pi、Pi 近似值、Pi 连分数、Pi 数字、Pi 迭代、Pi 平方、Spigot 算法

使用 Wolfram|Alpha 探索

WolframAlpha

更多尝试

参考文献

Adamchik, V. and Wagon, S. "A Simple Formula for pi." Amer. Math. Monthly 104, 852-855, 1997.Adamchik, V. and Wagon, S. "Pi: A 2000-Year Search Changes Direction." http://www-2.cs.cmu.edu/~adamchik/articles/pi.htm.Backhouse, N. "Note 79.36. Pancake Functions and Approximations to pi." Math. Gaz. 79, 371-374, 1995.Bailey, D. H. "Numerical Results on the Transcendence of Constants Involving pi, e, and Euler's Constant." Math. Comput. 50, 275-281, 1988a.Bailey, D. H.; Borwein, J. M.; Calkin, N. J.; Girgensohn, R.; Luke, D. R.; and Moll, V. H. Experimental Mathematics in Action. Wellesley, MA: A K Peters, 2007.Bailey, D. H.; Borwein, P.; and Plouffe, S. "On the Rapid Computation of Various Polylogarithmic Constants." Math. Comput. 66, 903-913, 1997.Beck, G. and Trott, M. "Calculating Pi: From Antiquity to Moderns Times." http://documents.wolfram.com/mathematica/Demos/Notebooks/CalculatingPi.html.Beckmann, P. A History of Pi, 3rd ed. New York: Dorset Press, 1989.Beeler, M. et al. Item 140 in Beeler, M.; Gosper, R. W.; and Schroeppel, R. HAKMEM. Cambridge, MA: MIT Artificial Intelligence Laboratory, Memo AIM-239, p. 69, Feb. 1972. http://www.inwap.com/pdp10/hbaker/hakmem/pi.html#item140.Berndt, B. C. Ramanujan's Notebooks, Part IV. New York: Springer-Verlag, 1994.Beukers, F. "A Rational Approximation to pi." Nieuw Arch. Wisk. 5, 372-379, 2000.Blatner, D. The Joy of Pi. New York: Walker, 1997.Boros, G. and Moll, V. Irresistible Integrals: Symbolics, Analysis and Experiments in the Evaluation of Integrals. Cambridge, England: Cambridge University Press, 2004.Borwein, J. and Bailey, D. Mathematics by Experiment: Plausible Reasoning in the 21st Century. Wellesley, MA: A K Peters, 2003.Borwein, J.; Bailey, D.; and Girgensohn, R. Experimentation in Mathematics: Computational Paths to Discovery. Wellesley, MA: A K Peters, 2004.Borwein, J. M. and Borwein, P. B. Pi & the AGM: A Study in Analytic Number Theory and Computational Complexity. New York: Wiley, 1987a.Borwein, J. M. and Borwein, P. B. "Ramanujan's Rational and Algebraic Series for 1/pi." Indian J. Math. 51, 147-160, 1987b.Borwein, J. M. and Borwein, P. B. "More Ramanujan-Type Series for 1/pi." In Ramanujan Revisited: Proceedings of the Centenary Conference, University of Illinois at Urbana-Champaign, June 1-5, 1987 (Ed. G. E. Andrews, B. C. Berndt, and R. A. Rankin). New York: Academic Press, pp. 359-374, 1988.Borwein, J. M. and Borwein, P. B. "Class Number Three Ramanujan Type Series for 1/pi." J. Comput. Appl. Math. 46, 281-290, 1993.Borwein, J. M.; Borwein, P. B.; and Bailey, D. H. "Ramanujan, Modular Equations, and Approximations to Pi, or How to Compute One Billion Digits of Pi." Amer. Math. Monthly 96, 201-219, 1989.Borwein, J. M.; Borwein, D.; and Galway, W. F. "Finding and Excluding b-ary Machin-Type Individual Digit Formulae." Canad. J. Math. 56, 897-925, 2004.Castellanos, D. "The Ubiquitous Pi. Part I." Math. Mag. 61, 67-98, 1988a.Castellanos, D. "The Ubiquitous Pi. Part II." Math. Mag. 61, 148-163, 1988b.Chudnovsky, D. V. and Chudnovsky, G. V. "Approximations and Complex Multiplication According to Ramanujan." In Ramanujan Revisited: Proceedings of the Centenary Conference, University of Illinois at Urbana-Champaign, June 1-5, 1987 (Ed. G. E. Andrews, B. C. Berndt, and R. A. Rankin). Boston, MA: Academic Press, pp. 375-472, 1987.Dalzell, D. P. "On 22/7." J. London Math. Soc. 19, 133-134, 1944.Dalzell, D. P. "On 22/7 and 355/113." Eureka 34, 10-13, 1971.Ferguson, H. R. P.; Bailey, D. H.; and Arno, S. "Analysis of PSLQ, An Integer Relation Finding Algorithm." Math. Comput. 68, 351-369, 1999.Finch, S. R. Mathematical Constants. Cambridge, England: Cambridge University Press, 2003.Flajolet, P. and Vardi, I. "Zeta Function Expansions of Classical Constants." Unpublished manuscript. 1996. http://algo.inria.fr/flajolet/Publications/landau.ps.Gasper, G. "Re: confluent pi." [email protected] posting, Aug. 18, 1996.Gosper, R. W. [email protected] posting, Sept. 1996.Gosper, R. W. "a product." [email protected] posting, Sept. 27, 1996.Gourdon, X. and Sebah, P. "Collection of Series for pi." http://numbers.computation.free.fr/Constants/Pi/piSeries.html.Guillera, J. "Some Binomial Series Obtained by the WZ-Method." Adv. Appl. Math. 29, 599-603, 2002.Guillera, J. "About a New Kind of Ramanujan-Type Series." Exp. Math. 12, 507-510, 2003.Guillera, J. "Generators of Some Ramanujan Formulas." Ramanujan J. 11, 41-48, 2006.Hardy, G. H. "A Chapter from Ramanujan's Note-Book." Proc. Cambridge Philos. Soc. 21, 492-503, 1923.Hardy, G. H. "Some Formulae of Ramanujan." Proc. London Math. Soc. (Records of Proceedings at Meetings) 22, xii-xiii, 1924.Hardy, G. H. Ramanujan: Twelve Lectures on Subjects Suggested by His Life and Work, 3rd ed. New York: Chelsea, 1999.Le Lionnais, F. Les nombres remarquables. Paris: Hermann, 1983.Lucas, S. K. "Integral Proofs that 355/113>pi." Austral. Math. Soc. Gaz. 32, 263-266, 2005.MathPages. "Rounding Up to Pi." http://www.mathpages.com/home/kmath001.htm.Plouffe, S. "Identities Inspired from Ramanujan Notebooks (Part 2)." Apr. 2006. http://www.lacim.uqam.ca/~plouffe/inspired2.pdf.Plouffe, S. "A Formula for the n'th Decimal Digit or Binary of pi and pi^n." https://arxiv.org/abs/2201.12601. 29 Jan 2022.Rabinowitz, S. and Wagon, S. "A Spigot Algorithm for the Digits of pi." Amer. Math. Monthly 102, 195-203, 1995.Ramanujan, S. "Modular Equations and Approximations to pi." Quart. J. Pure. Appl. Math. 45, 350-372, 1913-1914.Sloane, N. J. A. Sequences A054387 and A054388 in "The On-Line Encyclopedia of Integer Sequences."Smith, D. History of Mathematics, Vol. 2. New York: Dover, 1953.Sondow, J. "Problem 88." Math Horizons, pp. 32 and 34, Sept. 1997.Sondow, J. "A Faster Product for pi and a New Integral for ln(pi/2)." Amer. Math. Monthly 112, 729-734, 2005.Vardi, I. Computational Recreations in Mathematica. Reading, MA: Addison-Wesley, p. 159, 1991.Vieta, F. Uriorum de rebus mathematicis responsorum. Liber VII. 1593. Reprinted in New York: Georg Olms, pp. 398-400 and 436-446, 1970.Wells, D. The Penguin Dictionary of Curious and Interesting Numbers. Middlesex, England: Penguin Books, 1986.Wolfram Research, Inc. "Some Notes On Internal Implementation: Mathematical Constants." http://reference.wolfram.com/language/tutorial/SomeNotesOnInternalImplementation.html.

在 Wolfram|Alpha 中引用

Pi 公式

引用为

Weisstein, Eric W. "Pi 公式。" 来自 MathWorld--Wolfram Web 资源。 https://mathworld.net.cn/PiFormulas.html

学科分类