[論文レビュー] Kedlaya's algorithm in larger characteristic

David Harvey|ArXiv.org|Oct 31, 2006

Digital Image Processing Techniques参考文献 8被引用数 10

ひとこと要約

この論文は、有限体上の超楕円曲線のゼータ関数を計算するためのKedlayaのアルゴリズムの最適化版を提示する。多項式係数を持つ線形再帰を解くためのベビーステップ/ジャイアントステップ法を活用することで、実行時間の特性 $p$ への依存度を線形から $p^{1/2}$ に低減する。主な革新点は、微分の新しい級数表現と再帰の解法技術の向上であり、これにより $p$-進フロベニウス行列の計算が高速化される。

ABSTRACT

We show that the linear dependence on $p$ of the running time of Kedlaya's point-counting algorithm in characteristic $p$ may be reduced to $p^{1/2}$.

研究の動機と目的

• 超楕円曲線の $p$-進コホロジーにおけるフロベニウス行列の計算における、Kedlayaのアルゴリズムの実行時間における $p$-依存度を線形から $p^{1/2}$ に低減すること。
• $N$（精度）が $g$ や $n$ に依存しないように、楕円曲線上の点の $p$-進高さの計算をより効率的に行えるようにすること。
• 特に大きな $p$ および固定された $g$、$N$、$n$ に対して、既存のアルゴリズムを上回る漸近的複雑度を達成すること。
• 大規模な素数体、特に暗号的サイズに近いものも含む、大規模な素数体上の曲線を扱えるSageにおける実装を提供すること。

提案手法

• Kedlayaのフロベニウス作用の級数近似を、$p$ に依存しない $g$ と $N$ のみに依存する項数を持つ新しい表現に置き換える。
• ChudnovskyとChudnovskyの手法に基づくベビーステップ/ジャイアントステップ法を用いて、多項式係数を持つ線形再帰を解き、コホロジー関係のバッチ簡約を可能にする。
• Bostan, Gaudry, および Schost が提案した修正版Chudnovskyアルゴリズムを採用し、$O(\log(pN))$ の要因で実行時間を短縮する。
• コホロジー関係を用いて、フロベニウス作用の像を基底微分形式の線形結合として表現する還元アルゴリズムを採用し、フロベニウス行列を導出する。
• すべての計算を $p$-進設定で行い、$p^N$ の精度まで計算する。対数的要因を吸収するため、$\widetilde{O}$ 表記を用いる。
• GMP および NTL を最適化して使用したSageにおける実装を実施し、種数2、3、4の曲線における性能を検証する。

実験結果

リサーチクエスチョン

• RQ1 Kedlayaのアルゴリズムの $p$-依存度を、$\mathbf{F}_p$ 上の超楕円曲線に対して、$O(p)$ から $O(p^{1/2})$ に低減できるか。
• RQ2 フロベニウス作用の像を基底で表現するコホロジー還元ステップを、再帰構造を活用することでどのように高速化できるか。
• RQ3 どの代替的級数表現が、項数を $p$ に依存させず、$p$ に非線形な依存度を実現できるか。
• RQ4 Chudnovskyのベビーステップ/ジャイアントステップ法による再帰の解法を、$p$-進コホロジー計算に適応・最適化できるか。
• RQ5 この改善されたアルゴリズムは、$p \sim 2^{50}$ のような大規模な素数体上の曲線に対し、実行可能時間とメモリ使用量でゼータ関数を計算できるか。

主な発見

• 精度 $p^N$ での $p$-進フロベニウス行列の計算に要する実行時間は $\widetilde{O}(p^{1/2}N^{5/2}g^{\omega}n + N^4g^4n\log p)$ であり、$p$-依存度が線形から $p^{1/2}$ に低減された。
• 固定された $g$、$N$、$n$ に対して、実行時間は $\widetilde{O}(p^{1/2})$ となる。これにより、$p$ が大きい場合にはKedlayaのオリジナルアルゴリズムを上回る。
• 種数3の曲線に対して、$p = 2^{50} - 27$ の $\mathbf{F}_p$ 上で、40時間、16 GB RAM を使用してゼータ関数を $p$ を法として計算した。得られたヤコビアンの位数は約 $2^{150}$ であった。
• 種数4の曲線に対して、$p = 2^{44} + 7$ の $\mathbf{F}_p$ 上で、$p^2$ を法としてフロベニウス行列を45時間、34 GB RAM で計算し、ゼータ関数の候補を4つ得た。
• Sage 2.5.1以降の実装は、種数2、3、4の曲線に対してフロベニウス行列を正しく計算できており、性能は期待通りにスケーリングされている：$p$ が4倍になると実行時間は約2倍に増加する。
• 分母の取り扱いの改善とより洗練された $p$-進的アプローチにより、Bostan, Gaudry, および Schost の $\widetilde{O}(p^{1/2}g^{3/2+\omega}n)$ 法よりも $g^{3/2}$ の要因で性能が優れている。