On Wagstaff primes in the $k$-Lucas number sequence

以下は、Herbert Batte 氏による論文「On Wagstaff primes in the k-Lucas number sequence（k-ルカス数列におけるワグスタッフ素数について）」の技術的概要です。

1. 問題の背景と定義

本論文は、数論における線形漸化式数列と特殊な素数の交差点を研究対象としています。具体的には、**k-一般化ルカス数（k-generalized Lucas numbers）**の数列の中に、**ワグスタッフ素数（Wagstaff primes）**が含まれるかどうかを決定するディオファントス方程式を解くことを目的としています。

• k-一般化ルカス数 $L_n^{(k)}$:
  整数 $k \ge 2$ に対して、以下の漸化式と初期条件で定義されます。
  $$L_n^{(k)} = L_{n-1}^{(k)} + L_{n-2}^{(k)} + \dots + L_{n-k}^{(k)} \quad (n \ge 2)$$
  初期条件：$L_0^{(k)} = 2, L_1^{(k)} = 1$、および $2-k \le j \le -1$ の範囲で $L_j^{(k)} = 0$。
  （$k=2$ の場合は古典的なルカス数列 $2, 1, 3, 4, 7, 11, \dots$ となります。）

• ワグスタッフ素数:
  奇素数 $p$ に対して、$W = \frac{2^p + 1}{3}$ の形で表される素数です（例：3, 11, 43, 683 など）。

• 研究課題:
  以下のディオファントス方程式を満たす非負整数 $n, k$ と奇素数 $p$ の組 $(n, k, p)$ をすべて求めます。
  $$L_n^{(k)} = \frac{2^p + 1}{3} \quad \dots (1.1)$$
  先行研究（Rezaiguia ら）では、同様の方程式が k-一般化フィボナッチ数列 $F_n^{(k)}$ に対して解かれた際、解が $3$ のみであることが示されました。本論文は、初期値が異なるルカス数列において、より多くの解が存在する可能性を検証し、完全な解の集合を決定することを狙いとしています。

2. 手法とアプローチ

本論文は、**対数の線形形式（Linear forms in logarithms）**の理論と、**LLL 簡約法（LLL reduction method）**を組み合わせた数論的アプローチを採用しています。

2.1 対数の線形形式による上限評価

• ビネ公式（Binet formula）の利用:
  k-ルカス数は、特性多項式 $\Psi_k(x) = x^k - x^{k-1} - \dots - 1$ の最大実根 $\alpha$ を用いて近似できます。主要項は $f_k(\alpha)(2\alpha - 1)\alpha^{n-1}$ であり、残りの項は非常に小さいことが示されます。
• Matveev の定理の適用:
  方程式 (1.1) を変形し、$\Lambda = b_1 \log \gamma_1 + \dots + b_t \log \gamma_t$ の形の非零な線形形式を構成します。ここで、$\gamma_i$ は代数的数、$b_i$ は整数係数です。
  Matveev の定理を用いて、この線形形式 $\Lambda$ の絶対値の下限（lower bound）を導出します。一方、方程式の構造から得られる上限（upper bound）と比較することで、変数 $n$ と $k$ の間に成り立つ不等式を導き出します。
• 結果:
  この比較により、$n$ が $k$ の関数として非常に緩やかに増加する上限（$n < 2.0 \cdot 10^{15} k^4 (\log k)^3$）が得られます。

2.2 LLL 簡約法による上限の縮小

• 問題点:
  上記の Matveev による上限は理論的には有効ですが、数値的に計算するには依然として大きすぎます（$n$ が $10^{29}$ 程度になる可能性）。
• LLL 法の適用:
  得られた線形形式を格子（lattice）の近似問題として捉え、LLL 簡約法（Lenstra–Lenstra–Lovász algorithm）を適用します。
  具体的には、係数 $C$ を大きく設定した格子行列を構成し、LLL アルゴリズムによって基底を簡約化します。これにより、線形形式の最小非零値に対するより厳密な下限（$\delta$）が得られます。
• 結果:
  この手法を各 $k$ に対して適用することで、$n$ の上限を実用的な範囲（$n \le 141$）まで劇的に縮小することに成功しました。

2.3 場合分けと最終的な計算

• ケース 1: $2 \le n \le k$:
  この範囲では $L_n^{(k)} = 3 \cdot 2^{n-2}$ となるため、方程式は $9 \cdot 2^{n-2} = 2^p + 1$ に帰着します。$n=2$ のみ解を持ち（$p=3$）、$n \ge 3$ では偶数と奇数の矛盾により解なしとなります。
• ケース 2: $k > 190$:
  上記の上限評価と、$n < 2^{k/2}$ の条件下でのより精密な近似式（文献 [1] の式 (2.9)）を用いることで、$k > 190$ においては解が存在しないことを証明します。
• ケース 3: $k \le 190$ かつ $n > k$:
  LLL 法により $n \le 141$ まで縮小された範囲において、$k \in [2, 190]$、$n \in [k+1, 141]$、$p \in [3, n+3]$ の全組み合わせを SageMath 10.6 を用いて計算機的に検索します。

3. 主要な結果（定理）

論文の主要な定理（Theorem 1.1）は、方程式 (1.1) の整数解 $(n, k, p)$ は以下の集合に限定されることを示しています。

$$(n, k, p) \in \{(5, 2, 5), (6, 4, 7)\} \cup \{(2, k, 3) : k \ge 2\}$$

具体的には以下の 3 つのケースのみが解となります：

1. $L_2^{(k)} = 3$: 任意の $k \ge 2$ に対して、$3 = \frac{2^3 + 1}{3}$ となり、$p=3$ が解となります。
2. $L_5^{(2)} = 11$: 古典的なルカス数列（$k=2$）の第 5 項で、$11 = \frac{2^5 + 1}{3}$ となり、$p=5$ が解となります。
3. $L_6^{(4)} = 43$: $k=4$ のルカス数列の第 6 項で、$43 = \frac{2^7 + 1}{3}$ となり、$p=7$ が解となります。

4. 意義と貢献

• 先行研究の拡張:
  k-フィボナッチ数列におけるワグスタッフ素数の研究（Rezaiguia ら）を、初期値が異なる k-ルカス数列へと拡張しました。初期値の違いにより解の構造がどう変わるかという点で重要な知見を提供しています。
• 完全な解の決定:
  対数の線形形式と LLL 法を組み合わせることで、無限に広がる可能性のある解の空間を有限の範囲に限定し、計算機による完全な検索を可能にしました。これにより、解のリストがこれ以上存在しないことが厳密に証明されました。
• 数論的手法の適用:
  高度な数論的ツール（Matveev の定理、LLL 法）を具体的なディオファントス方程式の解決に応用するプロセスを詳細に示しており、同種の問題（線形漸化式と特殊素数の関係）を扱う際の手順として参考になります。

5. 結論

本論文は、k-一般化ルカス数列の中に現れるワグスタッフ素数を完全に分類することに成功しました。その結果、無限族の解（$n=2$ の場合）と、2 つの孤立した解（$(n,k,p)=(5,2,5)$ と $(6,4,7)$）のみが存在することが示されました。これは、線形漸化式数列における特殊素数の分布に関する理解を深める重要な成果です。