蒙哥马利约减

蒙哥马利解决的问题

求：
$$x\mathrm{mod}q$$
模运算在硬件中非常的慢，运比乘法慢的多，要进行优化，因此要想办法用乘
法及其它简单运算来替代模运算。
取$R=2^n$,则在R上的模运算和除运算都非常简单，假设n=32;
$$\frac{x}{R}=x>>32$$
$$x\mathrm{mod}R=x\&((1<<32)-1)\mathrm{mod}R$$
$$x\ast R=x<<32$$
$R$就是如上的一个数，它上面的除模乘运算都非常简单。

目标求$t\equiv x{R}^{-1}\mathrm{mod}N$,N为一个质数

$$R{R}^{-1}+N{N}^{\prime }=1;（要求R和N^{\prime }满足该关系）$$
$$\begin{array}{rl}x-(x\ast N^{\prime }\mathrm{mod}R)\ast N& \equiv x-(xN{N}^{\prime })\mathrm{mod}R\\ & \equiv x+(x(R{R}^{\prime }-1))\mathrm{mod}R\\ & \equiv 0\mathrm{mod}R(重要特性)\end{array}$$
由于上面的式子满足恰好能整除，所以可以直接使用除法而不用乘逆，并且，R特殊可用移位替代。

所有
$$\begin{array}{rl}(x-(x\ast N^{\prime }\mathrm{mod}R)\ast N)/R& \equiv (x-(x\ast N^{\prime }\mathrm{mod}R)\ast N)\ast R^{-1}\mathrm{mod}N\\ & \equiv x\ast R^{-1}\mathrm{mod}N\end{array}$$

Modular Multiplication Without Trial Division

例
$N=3457,R=2^{16}$
$3457\ast 12929-682\ast 2^{16}=1,3310=2^{16}\mathrm{mod}3457$
对于一个模3457的数$y$,有
$$\begin{gathered} y\equiv (y\ast R-(y\ast R)\ast N^{\prime }\mathrm{mod}R\ast N)/R\mathrm{mod}N \\ y\equiv (y\ast 3310-y\ast 3310\ast 12929\%2^{16}\ast 3457)/2^{16} \end{gathered}$$

int16_t montgomery_reduce(int32_t a)
{
  int32_t t;
  int16_t u;

  u = a * QINV;
  t = (int32_t)u * CTRU_Q;
  t = a - t;
  t >>= 16;
  return t;
}