We study three questions related to Machin's type formulas. The first one gives all two terms Machin formulas where both arctangent functions are evaluated $2$-integers, that is values of the form $b/2^a$ for some integers $a$ and~$b$. These formulas are computationally useful because multiplication or division by a power of two is a very fast operation for most computers. The second one presents a method for finding infinitely many formulas with $N$ terms. In the particular case $N=2$ the method is quite useful. It recovers most known formulas, gives some new ones, and allows to prove in an easy way that there are two terms Machin formulas with Lehmer measure as small as desired. Finally, we correct an oversight from previous result and give all Machin's type formulas with two terms involving arctangents of powers of the golden section.
论文ID : 2302.00154标题 : Three essays on Machin's type formulas作者 : Armengol Gasull (Universitat Autònoma de Barcelona), Florian Luca (Wits University & UNAM), Juan L. Varona (Universidad de La Rioja)分类 : math.NT (Number Theory)发表期刊 : Indagationes Mathematicae (to appear)论文链接 : https://arxiv.org/abs/2302.00154 本文研究了与Machin型公式相关的三个问题。首先，给出了所有两项Machin公式，其中反正切函数都在2-整数（形如b / 2 a b/2^a b / 2 a N N N N = 2 N=2 N = 2
Machin公式起源于1706年John Machin发现的恒等式：
4 arctan 1 5 − arctan 1 239 = π 4 4 \arctan\frac{1}{5} - \arctan\frac{1}{239} = \frac{\pi}{4} 4 arctan 5 1 − arctan 239 1 = 4 π
结合Gregory在1671年发现的反正切级数展开：
arctan x = ∑ m = 0 ∞ ( − 1 ) m 2 m + 1 x 2 m + 1 , ∣ x ∣ < 1 \arctan x = \sum_{m=0}^{\infty} \frac{(-1)^m}{2m+1} x^{2m+1}, |x| < 1 arctan x = ∑ m = 0 ∞ 2 m + 1 ( − 1 ) m x 2 m + 1 , ∣ x ∣ < 1
Machin利用这个公式计算出了π的100位小数。
计算效率问题 ：寻找计算π更高效的Machin型公式，特别是涉及2的幂次的形式，因为计算机对2的幂次运算极其高效理论完备性 ：系统性地分类和构造各种Machin型公式Lehmer测度最小化 ：寻找Lehmer测度任意小的公式以提高计算效率纠正历史遗漏 ：完善涉及黄金分割的Machin公式分类已知的Machin型公式数量有限 缺乏系统性的构造方法 对特殊形式（如涉及2的幂次）的公式研究不足 先前关于黄金分割的研究存在疏漏 完全分类了涉及2-整数的两项Machin公式 ：给出了10个零散解和2个参数族提出了"Machin公式机器"方法 ：能够系统性地生成具有N N N 证明了存在Lehmer测度任意小的两项Machin公式 ：提供了构造性证明纠正并完善了黄金分割幂次的Machin公式分类 ：发现了16个本质不同的恒等式寻找形如以下的所有解：
x 1 arctan ( z 1 ) + x 2 arctan ( z 2 ) = π 4 x_1 \arctan(z_1) + x_2 \arctan(z_2) = \frac{\pi}{4} x 1 arctan ( z 1 ) + x 2 arctan ( z 2 ) = 4 π x 1 , x 2 ∈ Q x_1, x_2 \in \mathbb{Q} x 1 , x 2 ∈ Q z k ∈ ( 0 , 1 ) z_k \in (0,1) z k ∈ ( 0 , 1 ) z k = 2 a k / b k z_k = 2^{a_k}/b_k z k = 2 a k / b k b k / 2 a k b_k/2^{a_k} b k / 2 a k
利用复数表示，将问题转化为：
( 1 + i z 1 ) 12 u 1 ( 1 + i z 2 ) 12 ∣ u 2 ∣ = ( 1 − i z 1 ) 12 u 1 ( 1 − i z 2 ) 12 ∣ u 2 ∣ (1 + iz_1)^{12u_1}(1 + iz_2)^{12|u_2|} = (1 - iz_1)^{12u_1}(1 - iz_2)^{12|u_2|} ( 1 + i z 1 ) 12 u 1 ( 1 + i z 2 ) 12∣ u 2 ∣ = ( 1 − i z 1 ) 12 u 1 ( 1 − i z 2 ) 12∣ u 2 ∣
通过分析Gaussian整数的因式分解，结合已知的丢番图方程解：
x 2 + 1 = 2 y n x^2 + 1 = 2y^n x 2 + 1 = 2 y n x 2 + 2 a = y n x^2 + 2^a = y^n x 2 + 2 a = y n 定义有理函数R j ( n , x ) R_j(n,x) R j ( n , x ) j = 0 , 1 , 2 , 3 j = 0,1,2,3 j = 0 , 1 , 2 , 3 d d x arctan ( R j ( n , x ) ) = n 1 + x 2 \frac{d}{dx}\arctan(R_j(n,x)) = \frac{n}{1+x^2} d x d arctan ( R j ( n , x )) = 1 + x 2 n
基于De Moivre公式，定义：
R j ( n , x ) = tan ( n arctan ( x ) + j π / 4 ) R_j(n,x) = \tan(n\arctan(x) + j\pi/4) R j ( n , x ) = tan ( n arctan ( x ) + jπ /4 )
具体表达式：
R 0 ( n , x ) = numer n ( x ) denom n ( x ) R_0(n,x) = \frac{\text{numer}_n(x)}{\text{denom}_n(x)} R 0 ( n , x ) = denom n ( x ) numer n ( x ) R 1 ( n , x ) = denom n ( x ) + numer n ( x ) denom n ( x ) − numer n ( x ) R_1(n,x) = \frac{\text{denom}_n(x) + \text{numer}_n(x)}{\text{denom}_n(x) - \text{numer}_n(x)} R 1 ( n , x ) = denom n ( x ) − numer n ( x ) denom n ( x ) + numer n ( x ) R 2 ( n , x ) = − 1 R 0 ( n , x ) R_2(n,x) = -\frac{1}{R_0(n,x)} R 2 ( n , x ) = − R 0 ( n , x ) 1 R 3 ( n , x ) = − 1 R 1 ( n , x ) R_3(n,x) = -\frac{1}{R_1(n,x)} R 3 ( n , x ) = − R 1 ( n , x ) 1 对于满足∑ k = 1 N r k = 0 \sum_{k=1}^N r_k = 0 ∑ k = 1 N r k = 0 r k r_k r k ∑ k = 1 N r k n k arctan ( R j k ( n k , x ) ) = r s π \sum_{k=1}^N \frac{r_k}{n_k} \arctan(R_{j_k}(n_k, x)) = \frac{r}{s}\pi ∑ k = 1 N n k r k arctan ( R j k ( n k , x )) = s r π
分类所有形如：
x 1 arctan ( ϕ a 1 ) + x 2 arctan ( ϕ a 2 ) = π 4 x_1 \arctan(\phi^{a_1}) + x_2 \arctan(\phi^{a_2}) = \frac{\pi}{4} x 1 arctan ( ϕ a 1 ) + x 2 arctan ( ϕ a 2 ) = 4 π ϕ = 1 + 5 2 \phi = \frac{1+\sqrt{5}}{2} ϕ = 2 1 + 5
通过分析K = Q ( i , 5 ) K = \mathbb{Q}(i,\sqrt{5}) K = Q ( i , 5 )
发现了10个零散解和2个参数族：
零散解示例：
( − 1 , 1 239 , 4 , 1 5 ) (-1, \frac{1}{239}, 4, \frac{1}{5}) ( − 1 , 239 1 , 4 , 5 1 ) ( 1 , 1 41 , 2 , 2 5 ) (1, \frac{1}{41}, 2, \frac{2}{5}) ( 1 , 41 1 , 2 , 5 2 ) 参数族：
( 1 , 1 2 a 2 + 1 , 1 , 2 a 2 2 a 2 + 1 ) (1, \frac{1}{2^{a^2}+1}, 1, \frac{2^{a^2}}{2^{a^2}+1}) ( 1 , 2 a 2 + 1 1 , 1 , 2 a 2 + 1 2 a 2 ) ( 1 , 1 2 a 2 − 1 , 1 , 2 a 2 − 1 2 a 2 ) (1, \frac{1}{2^{a^2}-1}, 1, \frac{2^{a^2}-1}{2^{a^2}}) ( 1 , 2 a 2 − 1 1 , 1 , 2 a 2 2 a 2 − 1 ) 证明了对任意ε > 0 \varepsilon > 0 ε > 0 n arctan 1 b 1 − arctan a 2 b 2 = π 4 n \arctan\frac{1}{b_1} - \arctan\frac{a_2}{b_2} = \frac{\pi}{4} n arctan b 1 1 − arctan b 2 a 2 = 4 π ε \varepsilon ε
构造性证明基于π的连分数收敛子，具体结果见表1和表2。
表1展示了基于π的连分数收敛子的公式，Lehmer测度从0.901429递减到0.172125。
表2展示了基于2 m 2^m 2 m
发现了16个本质不同的恒等式，包括：
π 4 = 1 3 arctan ( ϕ 3 ) + 1 3 arctan ( ϕ ) \frac{\pi}{4} = \frac{1}{3}\arctan(\phi^3) + \frac{1}{3}\arctan(\phi) 4 π = 3 1 arctan ( ϕ 3 ) + 3 1 arctan ( ϕ ) π 4 = − 1 2 arctan ( ϕ 5 ) + 3 2 arctan ( ϕ ) \frac{\pi}{4} = -\frac{1}{2}\arctan(\phi^5) + \frac{3}{2}\arctan(\phi) 4 π = − 2 1 arctan ( ϕ 5 ) + 2 3 arctan ( ϕ ) Störmer (1895) ：证明了形如x 1 arctan ( 1 / m 1 ) + x 2 arctan ( 1 / m 2 ) = π / 4 x_1\arctan(1/m_1) + x_2\arctan(1/m_2) = \pi/4 x 1 arctan ( 1/ m 1 ) + x 2 arctan ( 1/ m 2 ) = π /4 Lehmer (1938) ：提出了Lehmer测度作为计算效率的度量现代发展 ：Chudnovsky算法等更快的π计算方法本文的"Machin公式机器"扩展了4 中的结果，后者使用Rédei多项式的递归关系获得了R 3 ( n , x ) R_3(n,x) R 3 ( n , x )
完全解决了涉及2-整数的两项Machin公式分类问题 提供了系统性构造Machin公式的通用方法 证明了可以构造Lehmer测度任意小的公式 纠正并完善了黄金分割相关公式的分类 对于N > 2 N > 2 N > 2 虽然能构造小Lehmer测度的公式，但分数的分子分母位数很大，实际计算效率需要权衡 现代π计算已有更高效的算法（如Chudnovsky算法），Machin型公式的实用价值有限 扩展到更多项的Machin公式 研究其他特殊形式的反正切参数 探索与椭圆函数、模函数的联系 理论完备性 ：系统性地解决了几个重要的分类问题方法创新 ：提出的"Machin公式机器"具有一般性和实用性构造性证明 ：不仅证明了存在性，还给出了具体的构造方法历史纠错 ：认真对待并纠正了先前研究的疏漏巧妙的复数化 ：将实数问题转化为复数域的代数问题丢番图方程的应用 ：有效利用了已知的丢番图方程解连分数理论的应用 ：利用π的连分数性质构造小Lehmer测度公式实用价值有限 ：在现代高精度π计算中，这些公式的优势不明显复杂度问题 ：构造的小Lehmer测度公式涉及的整数位数很大推广性 ：方法对更一般情况的适用性存在限制理论贡献 ：为Machin型公式理论提供了系统性的框架计算数学 ：为特殊函数恒等式的研究提供了新工具数论应用 ：展示了代数数论方法在具体问题中的威力论文引用了24篇重要文献，涵盖了从Machin (1706)、Gregory (1671)的经典工作到现代的相关研究，特别是Störmer、Lehmer等人在该领域的奠基性工作。
