On Newton-Raphson iteration for multiplicative inverses modulo prime powers

Jean-Guillaume Dumas (1)

1 CASYS - Calculs Algébriques et Systèmes Dynamiques

We study algorithms for the fast computation of modular inverses. Newton-Raphson iteration over $p$-adic numbers gives a recurrence relation computing modular inverse modulo $p^m$, that is logarithmic in $m$. We solve the recurrence to obtain an explicit formula for the inverse. Then we study different implementation variants of this iteration and show that our explicit formula is interesting for small exponent values but slower or large exponent, say of more than $700$ bits. Overall we thus propose a hybrid combination of our explicit formula and the best asymptotic variants. This hybrid combination yields then a constant factor improvement, also for large exponents.

Mots clés: Hensel lifting.; prime power; Hensel lifting; modular inverse

Domaines: Calcul formel [cs.SC]; Logiciel mathématique [cs.MS]

Jean-Guillaume Dumas. On Newton-Raphson iteration for multiplicative inverses modulo prime powers. IEEE Transactions on Computers, 2014, 63 (8), pp.2106-2109. ⟨10.1109/TC.2013.94⟩. ⟨hal-00736701v5⟩