目录
原理:通过不断地将 ab 拆分成 (2a)b/2 [b为偶数] 或a+2a(b/2)直到b拆分为1结束而2a与res+a永远不会溢出(中间要不断对2a取模)
本文章是快速幂的进阶如果是新手就不建议看了
众所周知快速幂是为了解决数的高幂次运算爆时间问题。
为什么用快速乘?
但某些题目的运算数值过大比如要求结果对E10+2取模这时就很容易出现爆精度问题
比如: a = E10 b=E10 (即10的10次方)
a*b =E20
但long 的最大范围大概为E19 此时就会爆精度
解决方式就是快速乘
快速乘不仅比正常乘法运行速度更快,而且可以防止数据溢出
【1】快速乘
快速乘的目的:1.比普通乘快方法同快速幂的原理有一定相似,2. 防止n*m溢出
应用背景:当数据过大且取得模也很大再求快速幂的过程中有a*res也会数据溢出的情况
public static long fastmul(long a ,long b ){
long res = 0 ;
while (b>0){
if ((b&1)==1){
res = (res + a) % Mod;
}
b = b>>1;
a = (a+a)%Mod;
}
return res;
}
原理:通过不断地将 ab 拆分成 (2a)b/2 [b为偶数] 或a+2a(b/2)直到b拆分为1结束而2a与res+a永远不会溢出(中间要不断对2a取模)
举个例子: 计算19*22 要求运算结果对100取模 且存储数的最大值为200
计算 19*22
第一步 (19*2)*22/2 等价 38*11
第二步 38+38*2*(11/2) 等价 19+72*5
第三步 38+76+(76*2)%100 *5/2 等价 (38+76)%100+52*2 即 14+52*2
第四步 14 +(52*2)*2/1 等价 14+4
结果为18
如果题目背景是数值最大精度为200 那么19*22就直接爆精度了,而用快速乘的方法不会直接爆
嵌套快速乘的快速幂
//快速乘 计算a*b
public static long fastmul(long a ,long b ){
long res = 0 ;
while (b>0){
if ((b&1)==1){
res = (res + a) % Mod;
}
b = b>>1;
a = (a+a)%Mod;
}
return res;
}
//fastPower 计算a的b次方
public static long fastPower(long a,long b){
long res = 1;
while(b > 0){
if ((b&1) == 1){
res = fastmul(res,a);//嵌套快速乘
}
b = b>>1;
a = fastmul(a,a);
}
return res;
}