[論文レビュー] Two algorithms for evaluation of the Newman digit sum, and a new proof of Coquet's theorem

Vladimir Shevelev|arXiv (Cornell University)|Sep 6, 2007

Mathematical and Theoretical Analysis参考文献 1被引用数 8

ひとこと要約

本稿では、3を法とする剰余類における偶数および奇数の二進数の桁の和を持つ整数の個数の差として定義されるニューマン桁和の計算のための2つの効率的なアルゴリズムを提示する。この和に関するコーファーの鋭い上限・下限の初等的証明が行われ、$ S_{3,0}(3x) $ が $ \frac{2\sqrt{3}}{3}x^\lambda $ と $ \frac{55}{3}\left(\frac{3}{65}\right)^\lambda x^\lambda $ の間で有界であることを示し、$ \lambda = \frac{\ln 3}{\ln 4} $ である。また、実験的検証のための高速な計算手法が提供される。

ABSTRACT

We give two simple algorithms for the evaluation of difference between the numbers of multiples of 3 with even and odd binary digit sums in interval [0,x), and give an elementary proof of Coquet's sharp estimates for it.

研究の動機と目的

• 3を法とする剰余類0における、$ x $ 以下の整数で二進数の桁の和が偶数であるものの個数と奇数であるものの個数の差として定義されるニューマン桁和 $ S_{3,0}(x) $ の正確な評価のための効率的アルゴリズムの開発。
• コーファーの $ S_{3,0}(x) $ に関する鋭い漸近的推定値の、複素解析やフラクタル関数を避ける初等的な証明の提供。なお、元の式の誤植を是正する。
• $ S_{3,0}(3x) $ に対するタイトな明示的上限・下限の確立。$ \lambda = \frac{\ln 3}{\ln 4} $ を用いて、$ S_{3,0}(3x) $ が $ \Theta(x^\lambda) $ のオーダーで成長することを示し、計算アルゴリズムによるこれらの境界の検証。
• 数列の和 $ S_{3,1}(x) $、$ S_{3,2}(x) $、および $ S_{3\cdot 2^k,r}(x) $ を計算するための方法の拡張。これにより、算術級数における桁和分布への広範な応用が可能となる。

提案手法

• 本稿では、二進展開と剰余算術に基づく再帰的アルゴリズムを導入し、関数 $ t(y) = \sum (-1)^{k_i} $ を用いる。ここで $ 2^{k_i} $ は $ y $ の二進数の桁を表し、$ t(y) \equiv y \pmod{3} $ を満たす。この性質を用いて区間を分類し、既知の基本ケースに評価を還元する。
• 再帰的評価の基本的構成要素として、$ S_{3,0}(2^m) $ および $ S_{3,0}([2^n, 2^n + 2^m)) $ の閉形式表現を導出する。
• 自然数 $ N $ に対して $ S_{3,0}(N) $ の正確な評価を再帰的分解に基づき行う高速アルゴリズムを構築し、効率的な数値実験を可能にする。
• 補助関数 $ \nu(x) $、$ G(x) $、$ H(x) $ を定義し、それぞれ周期24、2、16で周期的であることを証明する。この性質がコーファーの結果の証明および鋭い境界の導出を支える。
• 鋭い境界の証明は、$ S_{3,0}(x) $、$ S_{3,0}(\lfloor x/4 \rfloor) $、および桁和 $ s(x) $ を含む関数方程式に依存し、周期的補正項が振動的挙動を捉える。
• 再帰的構造の鍵となるのは、$ x \mod 12 $ および $ \lfloor x/4 \rfloor \mod 2 $ に基づく場合分けであり、差分 $ \delta(x) = S_{3,0}(x) - S_{3,0}(4\lfloor x/4 \rfloor) $ の明示的公式が与えられる。

実験結果

リサーチクエスチョン

• RQ1 3を法とする剰余類0における、$ x $ 以下の整数で二進数の桁の和が偶数および奇数であるものの個数の差として定義されるニューマン桁和 $ S_{3,0}(x) $ の正確な評価のための効率的アルゴリズムは何か？
• RQ2 複素解析やフラクタル関数を用いずに、コーファーの $ S_{3,0}(x) $ に関する鋭い境界を初等的に再証明できるか？
• RQ3 $ S_{3,0}(3x) $ の正確な最良の上限・下限は何か？また、指数 $ \lambda = \frac{\ln 3}{\ln 4} $ とどのように関係しているか？
• RQ4 3を法とする桁和の再帰的構造をどのように活用し、$ S_{3,0}(x) $ の振動的挙動を捉える周期的補正関数を導出できるか？
• RQ5 この方法を $ S_{3,1}(x) $、$ S_{3,2}(x) $、および $ S_{3\cdot 2^k,r}(x) $ の計算に一般化し、鋭い推定値を得られるか？

主な発見

• すべての $ x \geq 2 $ に対して、$ \frac{2\sqrt{3}}{3}x^\lambda \leq S_{3,0}(3x) \leq \frac{55}{3}\left(\frac{3}{65}\right)^\lambda x^\lambda $ という鋭い境界が確立され、$ \lambda = \frac{\ln 3}{\ln 4} $ である。これはコーファーの元の推定値を確認・精錬するものである。
• コーファーの結果に対する初等的証明が与えられ、元の式における補正項 $ \eta(x) $ の誤植が是正された。元の記述では $ s(3x-3) $ とされていたが、正しくは $ s(3x-1) $ である。
• 関数 $ F(x) = (-1)^{s(x)}(S_{3,0}(x) - 3S_{3,0}(\lfloor x/4 \rfloor)) $ が周期24で周期的であることが証明され、ニューマン桁和の振動的挙動の構造的解釈が得られた。
• $ S_{3,0}(N) $ の正確な評価のための高速アルゴリズムが開発され、理論的境界の検証に適した数値実験に非常に適している。
• 本稿では、$ S_{3,0}(2^m) $ および $ S_{3,0}([2^n, 2^n + 2^m)) $ の明示的公式が導出され、$ m $ が偶数のときは $ S_{3,0}(2^m) = 2 \cdot 3^{m/2 - 1} $、奇数のときは $ 3^{(m-1)/2} $ であることが示され、再帰的評価の基盤となる。
• $ \delta(x) = S_{3,0}(x) - S_{3,0}(4\lfloor x/4 \rfloor) $ の差分が、$ x \mod 12 $ に依存する区分的関数として明示的に計算され、$ j=1,2,3 $ のとき $ (-1)^{s(x-j)} $、それ以外は0である。これにより、正確な再帰的計算が可能となる。