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 (数論)掲載予定誌 : Indagationes Mathematicae論文リンク : https://arxiv.org/abs/2302.00154 本論文はマッキン型公式に関連する3つの問題を研究している。第一に、逆正接関数が2-整数（b / 2 a b/2^a b / 2 a N N N N = 2 N=2 N = 2
マッキン公式は1706年にジョン・マッキンが発見した恒等式に起源する：
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 π
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
マッキンはこの公式を用いてπの100桁を計算した。
計算効率の問題 ：より効率的なマッキン型公式を探索すること、特に2の累乗に関わる形式。コンピュータは2の累乗演算に極めて効率的だからである理論的完全性 ：様々なマッキン型公式を体系的に分類・構成することレーマー測度の最小化 ：計算効率を向上させるためにレーマー測度が任意に小さい公式を探索すること歴史的欠落の修正 ：黄金比に関わるマッキン公式の分類を完善すること既知のマッキン型公式の数が限定的である 体系的な構成方法が不足している 特殊形式（例えば2の累乗に関わる公式）の研究が不十分である 黄金比に関する先行研究に欠落がある 2-整数を含む2項マッキン公式の完全分類 ：10個の散在解と2つのパラメータ族を与える「マッキン公式機械」方法の提案 ：N N N レーマー測度が任意に小さい2項マッキン公式の存在証明 ：構成的証明を提供する黄金比の累乗に関するマッキン公式分類の修正と完善 ：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 ∣
ガウス整数の因数分解の分析と既知のディオファントス方程式解を組み合わせる：
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
ド・モアブルの公式に基づいて、以下を定義する：
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はπの連分数収束子に基づく公式を示しており、レーマー測度は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 ( ϕ ) シュテルマー (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 レーマー (1938) ：計算効率の尺度としてレーマー測度を提案した現代の発展 ：チュドノフスキー・アルゴリズムなどより高速なπ計算方法本論文の「マッキン公式機械」は4 の結果を拡張したもので、後者はレーディ多項式の再帰関係を使用してR 3 ( n , x ) R_3(n,x) R 3 ( n , x )
2-整数を含む2項マッキン公式の分類問題を完全に解決した マッキン公式を体系的に構成するための一般的な方法を提供した レーマー測度が任意に小さい公式を構成できることを証明した 黄金比に関連する公式の分類を修正・完善した N > 2 N > 2 N > 2 小レーマー測度の公式を構成できるが、分数の分子分母の桁数が非常に大きく、実際の計算効率はトレードオフが必要である 現代のπ計算にはより効率的なアルゴリズム（チュドノフスキー・アルゴリズムなど）が存在し、マッキン型公式の実用的価値は限定的である より多くの項を持つマッキン公式への拡張 他の特殊形式の逆正接パラメータの研究 楕円関数、モジュラー関数との関連の探索 理論的完全性 ：いくつかの重要な分類問題を体系的に解決している方法の革新性 ：提案された「マッキン公式機械」は一般性と実用性を持つ構成的証明 ：存在性を証明するだけでなく、具体的な構成方法を与える歴史的修正 ：先行研究の欠落に真摯に対処し修正している巧妙な複素数化 ：実数問題を複素数域の代数問題に変換するディオファントス方程式の応用 ：既知のディオファントス方程式解を効果的に利用する連分数理論の応用 ：πの連分数性質を利用して小レーマー測度公式を構成する実用的価値の限定 ：現代の高精度π計算においてこれらの公式の優位性は明確でない複雑性の問題 ：構成された小レーマー測度公式は非常に大きな整数を含む推広性 ：方法のより一般的な場合への適用可能性に制限がある理論的貢献 ：マッキン型公式理論に体系的な枠組みを提供する計算数学 ：特殊関数恒等式の研究に新しい道具を提供する数論応用 ：代数的数論の方法が具体的問題にいかに威力を発揮するかを示す論文は24篇の重要な文献を引用しており、マッキン (1706)、グレゴリー (1671)の古典的業績から現代の関連研究まで、特にシュテルマー、レーマーなどこの分野の基礎的業績を網羅している。
