/****************************************************/
// 模取幂运算 计算a^b mod c
// 利用公式
// (a*b)mod(c) = ((a mod c )*b)mod c
/****************************************************/
int a_b_Mod_c(int a, int b, int c)
{//前提 a b c 都是正数
int digit[32];
int i, k, resualt = 1;
i = 0;
while(b)//把b化成2进制
{
digit[i++] = b%2;
b >>= 1;
}
//计算(a^b) mod c
for(k = i-1; k >= 0; k--)
{
resualt = (resualt * resualt) % c;
if(digit[k] == 1)
{
resualt = (resualt * a) % c;
}
}
return resualt;
}
转载自:http://blog.csdn.net/dremi/article/details/1568221
这个算法的思想我是从一本书上看到的,对合法的输入能很快的计算出结果来,其思想是利用
数学公式: (a * b ) mod c = (( a mod c) * b) mod c;
首先把 b 转化成二进制如: b0 b1 b2 b3..... b31 即 b = b0*231 + b1*230+......+ b31;
也就是把 ab = a ^ (b0*231 + b1*230+......+ b31) = [a(b0*2^31)] * [a(b1*2^30)] *..... * [ab31*2^0];
所以如果b的化成二进制的某位为0时,可以直接用这样算 resualt = (resualt * resualt)%c;
若为1 则为 resualt = ( resualt * a) % c;
例如:35 mod 4 ; 5转化成二进制为101,即35 = 3(2^2) * 3(2^0) = 34 * 31
所以35 mod 4 = 34 * 3 mod 4 = ((34)mod 4 * 3)mod 4;
31 mod 4 = 3, 32 mod 4 = (31 * 3 1) mod 4 ; 34 mod 4 = (32 *32)mod 4;