我们想计算 a k m o d n a^k\mod n akmodn 的数值。想法是先将 k k k 用二进制表达 ( k 0 k 1 ⋯ k d ) (k_0k_1\cdots k_d) (k0k1⋯kd) 然后逐次归纳计算 a 2 s m o d n a^{2^s}\mod n a2smodn 。
首先这里给出一个正整数用二进制表达的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;