[论文解读] 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

一句话总结

本文提出了两种高效算法，用于计算 Newman 数字和，定义为模 3 同余类中二进制数位和为偶数与奇数的整数个数之差。文章给出了 Coquet 对该和的精确界限的初等证明，证明了 $ 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.

研究动机与目标

• 开发用于精确计算 Newman 数字和 $ S_{3,0}(x) $ 的高效算法，其定义为在模 3 同余类 0 中，小于等于 $ x $ 的整数中二进制数位和为偶数与奇数的个数之差。
• 提供 Coquet 关于 $ S_{3,0}(x) $ 的渐近估计的全新初等证明，纠正原始表述中的笔误。
• 为 $ S_{3,0}(3x) $ 建立紧致的显式上下界，证明其增长为 $ \Theta(x^\lambda) $，其中 $ \lambda = \frac{\ln 3}{\ln 4} $，并通过计算算法验证这些界限。
• 将该方法推广至计算 $ S_{3,1}(x) $、$ S_{3,2}(x) $ 和 $ S_{3\cdot 2^k,r}(x) $，以实现对算术级数中数字和分布的更广泛应用。

提出的方法

• 本文提出一种基于二进制展开和模运算的递归算法，利用函数 $ t(y) = \sum (-1)^{k_i} $（其中 $ 2^{k_i} $ 为二进制位），该函数满足 $ t(y) \equiv y \pmod{3} $，用于对区间进行分类，并将求值归约为已知的基底情形。
• 推导出 $ S_{3,0}(2^m) $ 和 $ S_{3,0}([2^n, 2^n + 2^m)) $ 的闭式表达式，作为递归计算 $ S_{3,0}(x) $ 的基本构建模块。
• 基于 $ N $ 的二进制表示，构建一种快速算法以精确计算 $ S_{3,0}(N) $，从而实现高效的数值实验。
• 作者定义了辅助函数 $ \nu(x) $、$ G(x) $ 和 $ H(x) $，并证明其具有周期性（周期分别为 24、2 和 16），这是证明 Coquet 结果及推导精确界限的基础。
• 精确界限的证明依赖于一个函数方程，涉及 $ 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 是否存在高效的算法，可精确计算 Newman 数字和 $ S_{3,0}(x) $，其定义为在模 3 同余类 0 中，小于等于 $ x $ 的整数中二进制数位和为偶数与奇数的个数之差？
• RQ2 是否能使用一种初等方法重新证明 Coquet 对 $ 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} $，确认并改进了 Coquet 的原始估计。
• 给出了 Coquet 结果的初等证明，纠正了原始表述中校正项 $ \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 的周期函数，为 Newman 数字和的振荡行为提供了结构性解释。
• 开发了一种用于精确计算 $ 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} $，当 $ m $ 为奇数时 $ S_{3,0}(2^m) = 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)} $，其余情况为零，从而实现精确的递归计算。