此方法解决这样一个问题,就是a^b mod m 的余数是多少。
如果直接计算a^b,方次很大的时候,会溢出,而且时间很长。
当然指数很小的时候直接用自带的Math函数就行,如果指数很大的时候,可以用以下的方法。
原理是利用模运算的分配率:
举例 56 mod 3=(7x8) mod3=(7 mod 3)x(8 mod 3)=1x2=2
public static ulong ModPow(ulong baseValue, ulong exponent, ulong modulus) { // 初始化结果为1,因为任何数的0次方都是1 ulong result = 1; // 底数对模数取模,以防底数过大 baseValue %= modulus; // 当指数不为0时,继续循环 while (exponent > 0) { // 如果当前指数是奇数,将底数乘到结果上,并取模 if ((exponent & 1) == 1) result = (result * baseValue) % modulus; // 底数平方并取模 baseValue = (baseValue * baseValue) % modulus; // 指数右移一位(等同于除以2) exponent >>= 1; } // 返回结果 return result; }