[论文解读] A Successive Approximation Algorithm for Computing the Divisor Summatory Function

Richard Sladkey|arXiv (Cornell University)|Jun 15, 2012

Numerical Methods and Algorithms参考文献 2被引用 3

一句话总结

本文提出了一种新颖的逐次逼近算法，通过基于连分数邻近点（Farey neighbors）和坐标变换的几何方法，在 $ O(n^{1/3}) $ 时间和 $ O\log n $ 空间内计算除数求和函数 $ T(n) = \sum_{k=1}^n \tau(k) $。通过将双曲区域递归分解为切线与三角形子区域，该方法实现了亚平方根时间复杂度，打破了传统方法长期存在的 $ O(n^{1/2}) $ 时间瓶颈。

ABSTRACT

An algorithm is presented to compute isolated values of the divisor summatory function in O(n^(1/3)) time and O (log n) space. The algorithm is elementary and uses a geometric approach of successive approximation combined with coordinate transformation.

研究动机与目标

• 打破计算除数求和函数的标准方法所长期存在的 $ O(n^{1/2}) $ 时间复杂度瓶颈。
• 开发一种高效且基础的算法，以极低的空间消耗计算 $ T(n) $ 的孤立值。
• 将Voronoï误差分解的几何方法扩展为一种构造性、递归的格点计数算法。
• 通过利用对称性、Farey序列以及递归区域划分，实现亚平方根时间复杂度。

提出的方法

• 该算法基于Farey邻近点对应的斜率，将双曲线 $ xy = n $ 的区域逐层分解为一系列切线。
• 每个由两条切线与一条双曲线弧围成的区域，均通过坐标变换映射到新坐标系，以简化格点计数。
• 三角形区域中的格点通过基于三角数的公式 $ \Delta(i) = i(i+1)/2 $ 进行计数，并根据边界排除情况进行调整。
• 通过引入第三条切线，将上层区域递归划分为更小的曲线三角形。
• 坐标变换将原始的 $ xy = n $ 双曲线映射为一般二次型，从而实现对面积和格点的高效计算。
• 该算法结合部分和 $ S(n, x_1, x_2) = \sum_{x=x_1}^{x_2} \lfloor n/x \rfloor $ 与递归区域处理，避免冗余计算。

实验结果

研究问题

• RQ1 能否基于Farey序列的几何递归方法，在 $ O(n^{1/3}) $ 时间内计算除数求和函数？
• RQ2 如何利用双曲线 $ xy \leq n $ 的对称性与格点结构，将时间复杂度降低至 $ O(n^{1/2}) $ 以下？
• RQ3 Farey邻近点与坐标变换在高效计算曲线区域中格点数量的过程中起到何种作用？
• RQ4 通过使用切线对双曲区域进行递归划分，能否获得一个具有证明性亚平方根复杂度且空间消耗低的算法？
• RQ5 当将该方法扩展至高阶求和函数 $ T_k(n) $ 时，其性能如何扩展？

主要发现

• 该算法以 $ O(n^{1/3}) $ 时间和 $ O(\log n) $ 空间计算 $ T(n) $，相较于标准的 $ O(n^{1/2}) $ 方法有显著改进。
• 该方法基于Farey邻近点与切线的递归区域分解，通过多边形逼近双曲区域。
• 通过坐标变换与三角形区域中的格点计数，避免了对 $ O(n) $ 项的直接求和。
• 该方法可推广至高阶求和函数，对 $ T_k(n) $ 可达 $ O(n^{1 - 4/(3k)}) $ 复杂度，其中 $ T_3(n) $ 的复杂度为 $ O(n^{5/9}) $。
• 该算法易于并行化，尤其在求和与区域处理阶段，可在多核系统上实现可扩展的性能。
• 通过细致的壳层分解与子区域的递归细化，该方法有效避免了重复计数与三重计数。