[论文解读] Constructive Polynomial Partitioning for Algebraic Curves in $\\mathbb{R}^3$ with Applications

Boris Aronov, Esther Ezra+1 more|arXiv (Cornell University)|Jan 1, 2019

Computational Geometry and Mesh Generation参考文献 24被引用 8

一句话总结

本文提出了三维实空间中代数曲线多项式分割的首个高效算法构造，实现将空间分解为 O(D³ log³D) 个胞腔，每个胞腔与 O(n/D²) 条曲线相交，构造时间复杂度为 O(n²)。该方法被应用于在 O(n^{3/2+ε}) 时间内消除 n 条非竖直直线之间的深度环，实现近乎最优的性能，并提供了一种构造性、算法化的替代方案，以取代 Guth 的存在性分割结果。

ABSTRACT

In 2015, Guth proved that for any set of $k$-dimensional bounded complexity varieties in $\\mathbb{R}^d$ and for any positive integer $D$, there exists a polynomial of degree at most $D$ whose zero set divides $\\mathbb{R}^d$ into open connected sets, so that only a small fraction of the given varieties intersect each of these sets. Guth's result generalized an earlier result of Guth and Katz for points. Guth's proof relies on a variant of the Borsuk-Ulam theorem, and for $k>0$, it is unknown how to obtain an explicit representation of such a partitioning polynomial and how to construct it efficiently. In particular, it is unknown how to effectively construct such a polynomial for bounded-degree algebraic curves (or even lines) in $\\mathbb{R}^3$. We present an efficient algorithmic construction for this setting. Given a set of $n$ input algebraic curves and a positive integer $D$, we efficiently construct a decomposition of space into $O(D^3\\log^3{D})$ open "cells," each of which meets $O(n/D^2)$ curves from the input. The construction time is $O(n^2)$. For the case of lines in $3$-space we present an improved implementation, whose running time is $O(n^{4/3} \\log^{O(1)} n)$. The constant of proportionality in both time bounds depends on $D$ and the maximum degree of the polynomials defining the input curves. As an application, we revisit the problem of eliminating depth cycles among non-vertical lines in $3$-space, recently studied by Aronov and Sharir (2018), and show an algorithm that cuts $n$ such lines into $O(n^{3/2+\\epsilon})$ pieces that are depth-cycle free, for any $\\epsilon > 0$. The algorithm runs in $O(n^{3/2+\\epsilon})$ time, which is a considerable improvement over the previously known algorithms.

研究动机与目标

• 开发一种针对 ℝ³ 中代数曲线的高效、构造性多项式分割算法，以克服 Guth 原始结果中非构造性的局限。
• 提供一种实用的方法，将空间分解为每个胞腔仅与少量输入曲线相交的胞腔，从而支持算法应用。
• 将构造性分割方法应用于三维空间中非竖直直线之间深度环的消除问题，实现近乎最优的切割数与运行时间界限。

提出的方法

• 构造一个次数为 O(D log D) 的分割多项式 f，使得 ℝ³⧵Z(f) 被分解为 O(D³ log³D) 个胞腔，每个胞腔与 O(n/D²) 条曲线相交。
• 采用适配于常数次数曲线的多项式分割框架变体，结合代数几何与计算代数技术。
• 应用两阶段切割过程：首先沿分割多项式的零点集进行切割；其次在胞腔边界处进行额外切割，以解决残余深度环。
• 利用高效算法求解垂直半平面上的方程组 {P=0, ∂P/∂z=0}，在每条直线处以常数时间计算二级切割，假设采用合适的计算模型。
• 采用递归方法，各子问题独立处理，递推关系控制切割次数与运行时间。
• 整合垂直分解与梯形分解技术，以管理胞腔边界，确保深度环消除的正确性。

实验结果

研究问题

• RQ1 能否在算法上构造 ℝ³ 中代数曲线的多项式分割，而非如 Guth 结果中那样仅以存在性方式给出？
• RQ2 此类分割的构造时间复杂度是多少？对于常数次数曲线，是否可在近乎二次时间复杂度内完成？
• RQ3 该构造性分割能否用于在三维空间中以近乎最优的切割数与运行时间界限，消除非竖直直线之间的深度环？
• RQ4 与依赖矩阵乘法或近似算法的先前方法相比，该算法构造在效率与性能上表现如何？
• RQ5 该方法在多大程度上可推广至非竖直的常数次数代数曲线，而不仅限于直线？

主要发现

• 本文提出一种算法，在 O(n²) 时间内为 ℝ³ 中的 n 条代数曲线构造多项式分割，生成 O(D³ log³D) 个胞腔，每个胞腔与 O(n/D²) 条曲线相交。
• 对于 ℝ³ 中的直线，构造时间进一步优化至 O(n^{4/3} polylog n)，优于一般曲线情形。
• 该算法可在 O(n^{3/2+ε}) 时间内消除 n 条两两不相交的非竖直直线之间的所有深度环，且切割数为 O(n^{3/2+ε})，其中 ε > 0 任意。
• 切割数近乎最优，与 Aronov 和 Sharir（2018）的存在性上界仅相差多对数因子。
• 该方法为 Guth 的非构造性分割提供了构造性替代方案，使深度环消除可被实际实现。
• 该算法可轻松扩展至非竖直常数次数代数曲线，仅需少量修改，且保持相同的渐近界。