[论文解读] Counting Points on Hyperelliptic Curves using Monsky-Washnitzer Cohomology

Kiran S. Kedlaya|ArXiv.org|May 3, 2001

Advanced Numerical Analysis Techniques参考文献 8被引用 204

一句话总结

本文提出了一种在奇特征有限域上对超椭圆曲线进行点计数的多项式时间算法，利用Monsky-Washnitzer上同调计算弗罗贝尼乌斯特征多项式的$p$-进近似。该方法在具有有理魏尔斯特拉斯点的$g$-亏格曲线在$\bF_{p^n}$上的渐近时间复杂度为$O(g^{4+\theta}n^{3+\theta})$，相较于以往高亏格曲线的指数时间方法实现了显著改进。

ABSTRACT

We describe an algorithm for counting points on an arbitrary hyperelliptic curve over a finite field of odd characteristic, using Monsky-Washnitzer cohomology to compute a p-adic approximation to the characteristic polynomial of Frobenius. For fixed p, the asymptotic running time for a curve of genus g over the field of p^n elements is O(g^{4+ε} n^{3+ε}).

研究动机与目标

• 开发一种高效算法，用于在奇特征有限域上对超椭圆曲线的有理点进行计数，克服先前指数时间方法的局限性。
• 解决将椭圆曲线点计数方法（如Schoof方法）扩展至高亏格时的计算不可行性，因为雅可比矩阵的显式表示变得不切实际。
• 提供一种$p$-进上同调方法，避免依赖于显式雅可比模型或形式群计算，这些方法在亏格上呈指数复杂度。
• 在具有有理魏尔斯特拉斯点的曲线上，实现关于亏格$g$和域扩张次数$n$的多项式时间复杂度，从而实现高亏格曲线的实际计算。

提出的方法

• 该方法采用Monsky-Washnitzer刀子上同调，计算曲线在特征零上的提升的德拉姆上同调，通过坐标环的弱完备化确保积分的收敛性。
• 通过坐标环上的幂级数展开和弗罗贝尼乌斯提升，计算弗罗贝尼乌斯自同态在上同调空间上作用的$p$-进近似。
• 利用重复的弗罗贝尼乌斯扭转变换和通过重复平方法进行的矩阵指数运算，计算在$\bF_{p^n}$上关于上同调基的弗罗贝尼乌斯作用。
• 通过扩展的欧几里得算法对微分形式模曲线的定义方程进行约化，以保持精度和效率。
• 通过迭代向量乘法和矩阵求逆计算弗罗贝尼乌斯矩阵的特征多项式，从而提取zeta函数。
• 该方法利用$\bZ_p/(p^N)$-扩张中的快速算术运算，并结合优化的多项式求值方法（如Paterson-Stockmeyer），以降低时间复杂度。

实验结果

研究问题

• RQ1 能否使用$p$-进上同调方法在关于亏格和域扩张次数的多项式时间内计算有限域上超椭圆曲线的zeta函数？
• RQ2 Monsky-Washnitzer上同调在多大程度上可以替代Schoof或Satoh等传统方法，用于高亏格点计数？
• RQ3 是否可能在不依赖显式雅可比模型的情况下，实现关于亏格$g$和域扩张次数$n$的多项式时间复杂度？
• RQ4 有理魏尔斯特拉斯点的存在如何影响算法的可行性与复杂度？
• RQ5 该算法能否被改进以更高效地计算弗罗贝尼乌斯的牛顿多边形，而非完整的特征多项式？

主要发现

• 该算法在具有有理魏尔斯特拉斯点的$g$-亏格超椭圆曲线在$\bF_{p^n}$上的渐近时间复杂度为$O(g^{4+\theta}n^{3+\theta})$。
• 该方法在亏格$g$和域扩张次数$n$上均为多项式时间，与Schoof或Satoh等早期方法不同，后者在亏格或域大小上呈指数时间复杂度。
• 主要计算开销来自矩阵指数运算和特征多项式计算，两者均因$2g \times 2g$矩阵运算而需要$O(g^{4+\theta})$时间。
• $n^{3+\theta}$因子源于在剩余域上对弗罗贝尼乌斯自同态$\tau = \rho^k$的重复应用，每步指数运算需$O(\theta)$次应用。
• 该算法避免依赖形式群或雅可比方程，因此适用于传统模型难以处理的高亏格曲线。
• 在步骤2中存在一种潜在优化：通过并行化基向量上的弗罗贝尼乌斯作用，但步骤3中的$g^4$瓶颈无法通过并行化缓解。