[论文解读] Root Refinement for Real Polynomials

Michael Kerber, Michael Sagraloff|arXiv (Cornell University)|Apr 7, 2011

Polynomial and algebraic computation参考文献 14被引用 5

一句话总结

该论文提出 Aqir，一种经过认证的算法，通过自适应的近似算术，将无平方多项式实根的隔离区间精炼至 $2^{-L}$ 精度。通过将二次区间精炼与区间算术验证及工作精度控制相结合，该算法实现了接近最优的位复杂度 $\tilde{O}(d(d\Gamma_f + \Sigma_f)^2 + dL)$，在保持理论最优边界的同时，具有简单且实用的实现特性。

ABSTRACT

We consider the problem of approximating all real roots of a square-free polynomial $f$. Given isolating intervals, our algorithm refines each of them to a width of $2^{-L}$ or less, that is, each of the roots is approximated to $L$ bits after the binary point. Our method provides a certified answer for arbitrary real polynomials, only considering finite approximations of the polynomial coefficients and choosing a suitable working precision adaptively. In this way, we get a correct algorithm that is simple to implement and practically efficient. Our algorithm uses the quadratic interval refinement method; we adapt that method to be able to cope with inaccuracies when evaluating $f$, without sacrificing its quadratic convergence behavior. We prove a bound on the bit complexity of our algorithm in terms of the degree of the polynomial, the size and the separation of the roots, that is, parameters exclusively related to the geometric location of the roots. Our bound is near optimal and significantly improves previous work on integer polynomials. Furthermore, it essentially matches the best known theoretical bounds on root approximation which are obtained by very sophisticated algorithms. We also investigate the practical behavior of the algorithm and demonstrate how closely the practical performance matches our asymptotic bounds.

研究动机与目标

• 解决任意具有实系数的无平方多项式实根精炼问题。
• 提供一种经过认证、数值稳定的算法，仅使用系数的有限精度近似值。
• 实现接近理论最优的位复杂度，且独立于系数表示的复杂度。
• 通过自适应工作精度和区间算术验证，确保实际效率。
• 证明几何根参数（分离度、根界）决定复杂度，而非系数结构。

提出的方法

• 将二次区间精炼（QIR）适配为支持近似多项式求值与逐步增加的工作精度。
• 使用区间算术验证正确性并确保输出的认证性。
• 采用自适应精度策略，可在低于最坏情况边界的情况下成功运行。
• 将精炼过程划分为线性阶段（类似二分法）与二次阶段（快速收敛），并通过归一化确保二次收敛。
• 利用马勒界与导数估计来限制所需精度与收敛速率。
• 引入双倍宽度区间 $\tilde{I}$ 以确保正规性，并实现从第一步起的立即二次收敛。

实验结果

研究问题

• RQ1 能否设计一种简单且经过认证的算法，实现实根精炼的近似最优位复杂度？
• RQ2 实根精炼的复杂度如何依赖于几何根参数而非系数复杂度？
• RQ3 能否利用自适应精度与近似算术，在保持正确性与效率的同时实现目标？
• RQ4 确保正确且高效精炼所需的最小工作精度是多少？
• RQ5 实际性能与理论渐近界之间的接近程度如何？

主要发现

• 该算法实现了 $\tilde{O}(d(d\Gamma_f + \Sigma_f)^2 + dL)$ 的位复杂度，该复杂度接近最优，并与目前已知的最佳理论边界一致。
• 所需的工作精度被限制在 $\tilde{O}(d\Gamma_f + \Sigma_f + L)$ 位，且独立于系数表示。
• 复杂度仅依赖于几何参数：根分离度 $\Sigma_f$、对数根界 $\Gamma_f$ 与期望精度 $L$。
• 该算法具有鲁棒性与实用性，实验结果表明其在速度上显著优于精确算法。
• 理论复杂度界是紧致的，因为当 $m = \Theta(d)$ 时，仅读取输出就至少需要 $\Theta(dL)$ 次操作。
• 使用双倍宽度区间 $\tilde{I}$ 可确保正规性，从而实现从第一步起的立即二次收敛。