モンゴメリ乗算
モンゴメリ乗算(モンゴメリじょうざん)とは、特に時間のかかる除算を実質的に行うことなく、乗算・加減算・シフト演算のみで、高速に整数の積の剰余を求めることのできるアルゴリズムである。 数100ビットを超える法による冪剰余演算は、暗号理論の分野で重要な応用を持ち、モンゴメリ乗算を用いればこれを効率的に計算することができる。 Peter Montgomeryにより提案された。モンゴメリ法ともいう。
Contents
概要
モンゴメリ乗算のアイデアは、[math]N\gt 0[/math] を法とした合同算術に関して、演算したい値を、ある定数 [math]R[/math] を掛けた表現(ここではモンゴメリ表現と呼んでおく)に変換し、この表現によってすべての計算を行った後、最後に元の領域での表現に逆変換することである。
モンゴメリ表現での加減算はそのまま実行した後、負または [math]N[/math] 以上のときのみ [math]N[/math] の加減をするだけでよい。 しかし乗算では [math]R[/math] が余分に残るので、[math]R^{-1}[/math] を掛けて [math]N[/math] による剰余を求める処理を行う必要がある。 この処理をモンゴメリリダクションといい、[math]R[/math] をうまく選ぶことにより効率的に計算することができる。
アルゴリズム
モンゴメリリダクション
上述したように、モンゴメリリダクションは、モンゴメリ乗算の基本となる演算である。 [math]N[/math] を法とする [math]R[/math] に関する [math]T[/math] ( [math]0\le T \lt NR[/math] )のモンゴメリリダクションは
- [math]\mbox{MR}(T) = TR^{-1} \bmod N[/math]
と定義される。 ここで、[math]R[/math] は [math]R\gt N[/math] および [math]\gcd(N,R)=1[/math] (つまり [math]N[/math] と互いに素)なる任意の整数であり、[math]R^{-1}[/math] は [math]RR^{-1}\equiv 1 \pmod{N}[/math] なるモジュラー逆数である。
モンゴメリリダクションは次の手続きで計算できるが、[math]R[/math] として、[math]\bmod R[/math] と [math]/ R[/math] の計算が簡単になるような値(例えば2進数であれば2の冪)を選ぶことにより、除算を実質的に行う必要がなくなり効率的に計算できる。 ここで、[math]N'[/math] は [math]NN'\equiv-1 \pmod{R}[/math] なる値であり、[math]xR-yN=1[/math] を満たす [math]0\lt y\lt R[/math] として、拡張されたユークリッドの互除法や#Rが2の冪の時のN'の効率的な求め方などであらかじめ求めておく。 同時に得られる [math]0\lt x\lt N[/math] は上述の [math]R^{-1}[/math] である。
- [math]t \leftarrow (T+(TN' \bmod R)N) / R[/math]
- if [math]t \ge N[/math] then return [math]t-N[/math] else return [math]t[/math]
モンゴメリリダクション計算手続きの正当性
この手続きの正当性は次のように示される。 まず、
- [math]T+(TN' \bmod R)N \equiv T+TN'N \equiv T-T \equiv 0 \pmod{R}[/math]
より、[math]/ R[/math] が割り切れて [math]t[/math] が整数であることがわかる。 次に、
- [math]tR \equiv T+(TN' \bmod R)N \equiv T \pmod{N}[/math]
より、[math]t \equiv TR^{-1} \pmod{N}[/math] である。 最後に、
- [math]T \lt RN[/math], [math](TN' \bmod R)N \lt RN[/math]
より、[math]0 \le t = (T+(TN' \bmod R)N) / R \lt 2N[/math] であるため、手続きが返す値は [math]N[/math] より小さい。
モンゴメリ表現への変換と逆変換
整数 [math]0\le a\lt N[/math] をモンゴメリ表現 [math]A[/math] に変換するためには、[math]R[/math] を掛けて [math]N[/math] による剰余を求めればよいので、あらかじめ [math]R_2 = R^2 \bmod N[/math] を用意して、[math]A \leftarrow \mbox{MR}(aR_2)[/math] を求めればよい。
逆変換はモンゴメリリダクションそのものであり、[math]a \leftarrow \mbox{MR}(A)[/math] である。
乗算剰余演算
被乗数 [math]0\le a\lt N[/math] と乗数 [math]0\le b\lt N[/math] の乗算剰余 [math]c = ab \bmod N[/math] は、上述の変換と逆変換を用いて、
- [math]A \leftarrow \mbox{MR}(aR_2)[/math], [math]B \leftarrow \mbox{MR}(bR_2)[/math]
- [math]C \leftarrow \mbox{MR}(AB)[/math]
- [math]c \leftarrow \mbox{MR}(C)[/math]
により求められる。
しかし、単純に乗算剰余を1回だけ求めたい場合には、
- [math]c \leftarrow \mbox{MR}(\mbox{MR}(ab)R_2)[/math]
とする方が効率的である。
冪剰余演算
冪剰余 [math]a^k \bmod N[/math] を求めたい場合、まず [math]a[/math] をモンゴメリ表現に変換しておき、[math]A^k[/math] の計算に出現する乗算のたびに積にモンゴメリリダクションを行っていけば、最後の結果に対して逆変換することによって結果を得ることができる。
通常の冪の計算ではバイナリ法などが効率的であるが、冪剰余の計算においても、同様の効率化がそのまま利用できる。
Rが2の冪の時のN'の効率的な求め方
Rが2の冪である場合には、計算機向けの効率的な求め方が存在する。
[math]NN'\equiv-1 \pmod{R}[/math]
は
[math]NN'\equiv R-1 \pmod{R}[/math]
であり、Rが2の冪であるためR-1は二進数表記で全てのビットが1になる。また2の冪での剰余を求めることは即ち二進数表記での下位ビットを取り出すことである。なのでこれは実質的に、NN'の二進数表記の下位ビット全てが1になるN'を求めることに相当する。
int result = 0;
int t = 0;
int r = R;
int i = 1;
while( r > 1 ) /* Rのトップビットを除いたビット数分繰り返す */
{
if( !( t % 2 ) ) /* ゼロになっているビットがあったら、N'のその部分を1にする(NはRと互いに素なので必ず奇数) */
{
t += N; /* 掛け算だが、二進数一桁の掛け算なので実質は足し算 */
result += i; /* N'のその部分を1にする */
}
t /= 2; /* 必ず端数が出るが切り捨てる */
r /= 2; /* Rは2の冪なので、絶対端数は出ない */
i *= 2;
}
/* return result; この時点で N' == result */
参考文献
- Peter Montgomery, "Modular Multiplication Without Trial Division," Math. Computation, vol. 44, pp. 519–521, 1985.