[论文解读] Homotopy methods for multiplication modulo triangular sets

Alin Bostan, Muhammad F. I. Chowdhury+2 more|ArXiv.org|Jan 23, 2009

Polynomial and algebraic computation参考文献 22被引用 8

一句话总结

本论文提出了一种基于同伦的算法，用于快速计算多项式模三角集的乘法，利用变形参数的幂级数上的求值-插值技术。对于特定的三角集族（包括递归系统和柯西模），该算法实现了准线性时间复杂度 $ O(\tilde{u}) $，其中 $ u $ 为输入/输出规模，显著优于先前的方法。

ABSTRACT

We study the cost of multiplication modulo triangular families of polynomials. Following previous work by Li, Moreno Maza and Schost, we propose an algorithm that relies on homotopy and fast evaluation-interpolation techniques. We obtain a quasi-linear time complexity for substantial families of examples, for which no such result was known before. Applications are given to notably addition of algebraic numbers in small characteristic.

研究动机与目标

• 为解决在多变量和递归设置下，多项式模三角集乘法的高计算成本问题。
• 在标准方法导致单项式数量指数级膨胀的情况下，开发模三角集乘法的快速算法。
• 将基于同伦的求值-插值技术从单变量情形扩展至具有稀疏或低次结构的多变量三角集。
• 在小特征下，为代数数系统和柯西模提供乘法的准线性时间复杂度。
• 解决基于变形的插值过程中出现的精度问题。

提出的方法

• 引入变形参数 $ \eta $，构造新的三角集 $ \mathbf{S} = \eta\mathbf{T} + (1-\eta)\mathbf{U} $，其中 $ \mathbf{U} $ 的根位于基环中。
• 利用亨泽尔引理确保 $ \mathbf{S} $ 在 $ \mathsf{R}[[\eta]] $ 中有根，从而在形式幂级数上实现求值与插值。
• 在 $ \mathsf{R}[[\eta]] $ 上应用快速求值与插值算法，通过变形计算乘积 $ AB \mod \langle \mathbf{T} \rangle $。
• 利用同伦技术避免直接的多项式乘法与约化，后者在单项式数量上存在指数级膨胀。
• 通过将所需的 $ \eta $-进精度与 $ \mathbf{T} $ 的单项式支撑关联，控制幂级数计算中的精度。
• 在小特征下，利用牛顿恒等式与提升技术，从初等对称和恢复最终结果，确保准线性位复杂度。

实验结果

研究问题

• RQ1 当单项式支撑稀疏或次数较低时，是否可在准线性时间内完成多变量三角集模乘法？
• RQ2 如何将同伦与变形技术从单变量情形扩展至多变量三角集，以实现快速求值-插值？
• RQ3 在同伦方法中，为确保正确性与效率，幂级数计算需要何种精度控制？
• RQ4 所提方法是否可在小特征下对代数数系统与柯西模实现准线性复杂度？
• RQ5 同伦方法中的计算瓶颈是什么，如何加以缓解？

主要发现

• 对于形如 $ X_i^2 - 2X_{i-1} $ 或 $ X_i^2 - X_{i-1} $ 的三角集，模 $ \langle \mathbf{T} \rangle $ 的乘法可在 $ O(\tilde{\delta}) $ 时间内完成，其中 $ \delta = 2^n $ 为输入/输出规模。
• 对于由单变量多项式 $ F $（次数为 $ n $）导出的柯西模，模 $ \langle \mathbf{T} \rangle $ 的乘法同样为 $ O(\tilde{\delta}) $，其中 $ \delta = n! $。
• 该算法通过在幂级数上使用同伦与求值-插值，避免了直接乘法的指数级膨胀（代价为 $ \delta^{\log_2 3} $）。
• 通过在模 $ p^\alpha $ 下提升幂和，并利用牛顿恒等式恢复最终多项式，该方法在小特征下实现了准线性位复杂度。
• $ \eta $-进展开中的精度控制被证明与 $ \mathbf{T} $ 的单项式支撑密切相关，从而实现高效计算。
• 实验结果证实了该算法的实际效率，运行时间与准线性复杂度估计相符，并在性能上优于 Magma 的基于结式实现两个数量级。