[英]Why don't we need to perform modulo operation on every operands in Fast Power algorithm?
今天,我练习了一个难题“快速力量”,其中使用了一个公式: (a * b) % p = (a % p * b % p) % p
来计算(a^n)%p
,类似这样: 2^31 % 3 = 2
但是,当我找到所使用的答案时,我感到非常困惑((temp * temp) % b * a) % b;
解决n为奇数时的情况,例如2^3
(temp是(temp * temp) % b * a
递归或(temp * temp) % b
)。
它应该不是((temp * temp) % b * a%b) % b
吗?
由于根据此公式,所有内容都应先%b
次在一起。
应该不是
((temp * temp) % b * a % b) % b
吗?
否。对于a
,如果您事先知道a
不会溢出(a小于b),则无需进行修改。
这个想法是用于加法和乘法的模块化算术工作。 通常执行(a + b) % M = (a % M + b % M) % M
和(a * b) % M = (a % M * b % M) % M
这样的操作以避免(a * b)
和(a + b)
并将值保持在一定范围内。
例:
const int Mod = 7;
int a = 13;
int b = 12;
int b = b % Mod; // b now contains 5 which is certainly smaller than Mod
int x = (a % Mod * b) % Mod; // you won't need to mod b again if you know beforehand b is smaller than Mod
C ++幂函数的实现:
#define MOD 1000000007
// assuming x and n both be positive and initially smaller than Mod
int power(int x, int n) {
if(n == 0) return x;
int half = power(x, n / 2) % Mod;
int ret = (half * half) % Mod; // you didn't need to do (half % Mod * half % Mod) % Mod because you already know half is smaller than Mod and won't overflow.
// Modulas being performed on the multiplied output, so now ret will be smaller than Mod
if(n & 1) {
ret = (ret * x) % Mod; // you didn't need to do (ret % Mod * x % Mod) % Mod
// because you already know ret and x is smaller than Mod
}
return ret;
}
Mod是一项昂贵的操作。 因此,您应尽可能避免使用它。
声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.