平方乘积法求幂的模运算

我们想计算 $a^k\mathrm{m}\mathrm{o}\mathrm{d}n$ 的数值。想法是先将 $k$ 用二进制表达 $(k_0{k}_1\cdots k_d)$ 然后逐次归纳计算 $a^{2^s}\mathrm{m}\mathrm{o}\mathrm{d}n$ 。

首先这里给出一个正整数用二进制表达的C++代码

int k,d=0;
cin>>k;
int adic[100];
while(k!=0)
{
   if(k&1){adic[d]=1;k=(k-1)/2}
   else {adic[d=0];k=k/2}
   d++;
}//这时候adic[d] 就是a 的二进制表达了（从左到右，十进制是从右到左的位数）。

现在使用平方乘积算法

//沿用以上的代码块的记号，现在假设 adic[d] 就是 k 的二进制表达。
int b;
if(k==0) return b;
int A=a;
if(adic[0]==1)b=a;
for(int i=1;i<d;i++)
{
A=(A*A)%n;
if(k[i]==1)b=A*b%n;
}
return b;