来自百度文库 PPT作者 张鹏
只能说下载分数太贵了。。
基本概念及思想
¡对形如a^b mod m 的运算(b一般较大)
¡但a,b,m都在long型范围内
¡算法的主要思想是分治,分而治之。将大的问题分成若干个相似的较小的问题!
¡具体实现是用递归的方法!
举例
¡2^100mod 3
¡像这种运算如果先算出2^100 的值,然后再模上3,相信比较困难!
¡我们可以将100变小点
¡2^100=(2^50)^2 =((2^25)^2)^2=((((2^1)^2)^2)…)^2
¡若我们已经得出250 mod 3的值,我们便很简单地得出2^100mod 3的值。
¡按照上述的方法继续分下去…
¡最终肯定会得到2^1mod 3 这种情况!这样就好办了!这便是递归的边界,到此我们就可以返回2 mod 3的值了!
¡另一个问题,再分时有两种情况!
¡2^100=(2^50)2 , 100是偶数
¡2^99=(2^49)2 *2 , 99是奇数
¡第二种情况需要在第一种情况上乘上一次基数。
主要代码
long mod(long a,long b,long m)
{
if(!b) return 1; //边界处理
if(b==1) return a%m; //边界处理
long ans=mod(a,b/2,m); //进入下一层
ans=ans*ans%m; //返回值ans代表ab/2
if(b&1) ans=ans*a%m; //奇数情况处理
return ans; //返回ans代表abmodm
}
调用
¡上述函数是有三个参数,且有返回值!
¡它的功能是计算ab mod m的值!
¡如何在主函数里调用呢?
¡例如我们要计算2^100mod3的结果!我们可以首先定义一个变量ans,在主函数读入工作完成后,我们可以写ans=mod(2,100,3);
¡ans 里的数值便是2^100mod3的结果了。
递推法
¡用递推法首先要将b转换为二进制数,(b)10=(b)2。
¡举例:(45)10=(101101)2
¡545mod 3 =2
¡二进制数由0和1组成,它们分别代表着偶数和奇数。
¡递推时在0和1处的处理是有差别的。
参考代码
¡二进制转换:
¡long s=0;
¡while(b>0)
¡{
¡ t[++s] = b%2;
¡ b/=2;
¡}
¡递推过程:
¡long ans=1;
¡for(i=1;i<=s;i++)
¡if(t[i])
¡ ans=ans*ans*a%m;
¡else
¡ ans=ans*ans%m;
----------------------------------------------------------------
int modular_exponent(INT a,INT b,INT n){
int ret=1;
for (;b;b>>=1,a=a*a%n)
if (b&1)
ret=ret*a%n;
return ret;
}