[論文レビュー] Faster exponentials of power series

David Harvey|ArXiv.org|Nov 16, 2009

Coding theory and cryptography参考文献 1被引用数 4

ひとこと要約

本稿では、複素数形式の形式的べき級数環 $\mathbf{C}\llbracket x\rrbracket$ におけるべき級数 $f$ の指数関数 $\exp f$ を計算するための新しいアルゴリズムを提示する。時間計算量は $(13/6 + o(1))M(n)$ であり、以前の最良の境界である $(7/3 + o(1))M(n)$ よりも改善されている。本手法は、最初の半分の級数に対してヴァン・デル・ホーヴェンのブロックベースFFTアプローチを組み合わせ、精度を効率的に拡張するニュートン反復を用いることで、重複するFFT計算を削減し、計算量の定数因子を低減している。

ABSTRACT

We describe a new algorithm for computing exp(f) where f is a power series in C[[x]]. If M(n) denotes the cost of multiplying polynomials of degree n, the new algorithm costs (2.1666... + o(1)) M(n) to compute exp(f) to order n. This improves on the previous best result, namely (2.333... + o(1)) M(n).

研究の動機と目的

• 形式的べき級数 $f \in \mathbf{C}\llbracket x\rrbracket$ に対する $\exp f$ の計算の計算量を低減すること。
• この問題に対する以前の最良の時間計算量 $(7/3 + o(1))M(n)$ を改善すること。
• ヴァン・デル・ホーヴェンのアルゴリズムのように、すべてのブロックに対してフルFFTを実行しないことで、重複するFFT計算を回避すること。
• ブロックベース計算とニュートン反復による精度拡張を組み合わせることで、よりタイトな計算量の上限を達成すること。

提案手法

• アルゴリズムはヴァン・デル・ホーヴェンのブロック分解とFFTに基づく乗算を用い、$\exp f$ の最初の $s$ ブロックを $X^s$ を法として計算する。
• 初期の部分 $g_{[0]}$ の逆数を用いて対数微分の誤差を補正するニュートン反復を適用し、結果を $X^s$ から $X^{2s}$ へと拡張する。
• 主な構成要素として、$\delta f = x f'$ およびそのブロック適応版 $\delta_k$ である微分作用素が含まれる。これらは周波数領域における級数の効率的取り扱いを可能にする。
• FFTは級数の最初の半分でのみ実行され、2番目の半分に対しては冗長な変換を回避する。
• FFTと逆FFTを用いて、[Har09, Lemma 1] を活用し、1ブロックあたり $T(2m) + O(m(k+1))$ 時間でブロック積を効率的に計算する。
• 最終的な拡張ステップでは、誤差項 $\varepsilon_{[k]} = \delta_{k+s}^{-1} q_{[k+s]} - f_{[k+s]}$ を用いて $g_{[k+s]}$ を計算し、高精度な再構成を実現する。

実験結果

リサーチクエスチョン

• RQ1 ヴァン・デル・ホーヴェンの $(7/3 + o(1))M(n)$ の境界を上回る、べき級数の指数関数の計算がより効率的に行えるか？
• RQ2 ブロックベースのアルゴリズムにおける重複するFFT計算を、正しさを損なわずに排除できるか？
• RQ3 ニュートン反復をブロックベースFFT手法と効果的に組み合わせることで、全体の計算量を低減できるか？
• RQ4 現代のFFTおよび級数操作技術を用いた場合、$\exp f$ の計算量における最小の定数因子は何か？

主な発見

• 提案されたアルゴリズムは、$\exp f \bmod x^n$ を時間 $(13/6 + o(1))M(n)$ で計算でき、以前の最良の $(7/3 + o(1))M(n)$ よりも改善されている。
• この改善は、ブロックの2番目の半分に対してFFTを実行しないことで生じる。これは、サイズ $m$ の $s$ 個のブロックのFFTに相当するコストの削減に寄与する。
• アルゴリズムはヴァン・デル・ホーヴェンの方法で最初の $s$ ブロックを正しく計算し、その後対数微分のニュートン反復により精度を拡張する。
• 計算量は $(13s - 4)T(2m) + O(s^2 m)$ と解析され、標準的な仮定の下で、$(13/6 + o(1))M(n)$ の境界が得られる。
• 微分作用素 $\delta$ 及びその逆数の性質に依存した、ブロックインデックスに関する帰納法により、アルゴリズムの正しさが証明されている。
• 本手法は定数因子 $13/6 \approx 2.1667$ を達成し、以前の $7/3 \approx 2.3333$ よりも顕著な漸近的改善を実現している。