[论文解读] Algebraic osculation and factorization of sparse polynomials

Martin Weimann|ArXiv.org|Apr 1, 2009

Polynomial and algebraic computation参考文献 28被引用 3

一句话总结

本文提出了一种新颖的算法，通过在toric簇中利用代数接触（algebraic osculation）对稀疏二元多项式进行绝对因式分解，使用基于留数的判据，通过Newton多面体的几何结构检测不可约因子。主要贡献是提出了一种确定性的、具有几何感知能力的Hensel提升方法，相比经典的Galligo-Rupprecht算法，显著减少了重组情况。

ABSTRACT

We prove a theorem on algebraic osculation and we apply our result to the Computer Algebra problem of polynomial factorization. We consider X a smooth completion of the complex plane and D an effective divisor supported on the boundary of X. Our main result gives explicit conditions equivalent to that a given Cartier divisor on D extends to X. These osculation criterions are expressed with residues. We derive from this result a toric Hensel lifting which permits to compute the absolute factorization of a bivariate polynomial by taking in account the geometry of its Newton polytope. In particular, we reduce the number of possible recombinations when compared to the Galligo-Rupprecht algorithm.

研究动机与目标

• 解决现有LR算法在稀疏多项式中通用坐标变换效率低下的问题。
• 通过利用Newton多面体的几何结构，减少绝对因式分解中的可能重组情况。
• 开发一种避免概率假设和渐近近似方法的确定性绝对因式分解算法。
• 通过toric正规化，将基于留数的扩展判据推广至奇异曲线。
• 为使用边界交点数据和Cartier除子扩展提供理论基础。

提出的方法

• 使用$\mathbb{C}^2$的光滑toric紧化空间$X$，将由$f$定义的曲线$C$嵌入其中，以保持Newton多面体的信息。
• 通过条件$\langle \gamma, \Psi \rangle = 0$对所有$\Psi \in H^0(X, \Omega_X^2(D))$，判断$D = \partial X$上的Cartier除子$\gamma$是否能扩展至$X$。
• 利用Grothendieck留数和Dolbeault $\bar{\partial}$-分解，计算除子扩展的上同调障碍。
• 将因式分解问题转化为求解一个大小为$2\operatorname{Vol}(Q) + \deg(\gamma' \cdot \partial X)$的稀疏线性系统，未知数个数为$\operatorname{Card}(Q \cap \mathbb{Z}^2)$。
• 利用Minkowski和分解$N_f = P + Q$，建立线丛分解$\mathcal{O}_X(C) \simeq \mathcal{L}_P \otimes \mathcal{L}_Q$。
• 应用具有高阶隐函数计算（最高至阶数$na^2$）的toric Hensel提升过程，以确保因子的确定性检测。

实验结果

研究问题

• RQ1 能否基于留数的代数接触判据检测稀疏二元多项式的不可约因子？
• RQ2 如何利用Newton多面体的几何结构减少绝对因式分解中的重组情况？
• RQ3 在toric曲面的边界除子$D$上，Cartier除子何时能扩展至整个曲面？
• RQ4 与经典Hensel提升和Galligo-Rupprecht算法相比，所提方法在复杂度和精度上表现如何？
• RQ5 能否利用边界上的非toric奇点进一步减少因式分解中的重组情况？

主要发现

• 当且仅当对所有$\Psi \in H^0(X, \Omega_X^2(D))$都有$\langle \gamma, \Psi \rangle = 0$时，$D = \partial X$上的Cartier除子$\gamma$能扩展至$X$，从而提供一种基于留数的扩展判据。
• 与Galligo-Rupprecht算法相比，该算法通过使用多面体几何而非通用坐标，显著减少了可能的重组情况。
• 对于双次数为$(a,b)$的多项式，该方法分别计算$\gamma$的隐函数至阶数$b$和$a$，而GR算法需计算$a+b$个隐函数至阶数$a+b$。
• 新方法中隐函数所需的精度最高可达$na^2$，远超经典Hensel提升中的$2na = \deg(f)$界限。
• 当面多项式为无平方时，该算法成为基于零和的LR算法，通过留数条件实现确定性检测。
• 对于边界上存在奇点的曲线，非toric正规化$\widetilde{X}$允许使用定理1，且新增的例外除子可进一步减少重组数。