模取幂运算

当一个幂运算很大，而模为整型数时，通常的做法（先求幂再取模），结果很可能是溢出，因此我们可以利用以下的公式可以有效的避免这一问题

 

数学公式： (a * b ) mod C = (( a mod C) * b) mod C; (a, b, C 皆是是整数)

证明：

设：a=m+k,  m 的值为 C 的 n 倍数（n>=0，且为整数）, k < C

则:  (a * b ) mod C= ( ( m + k ) * b ) mod C  = ( mb + kb ) mod C = ( Cnb + kb ) mod C

       ( Cnb + kb ) mod C = kb mod C

因为：

      k = a mod C

所以：

     (a * b ) mod C = ( ( a mod C ) * b ) mod C

首先把 b 转化成二进制如： b₀ b₁ b₂ b₃..... b₃₁ 即 b ＝ b₀*2³¹ + b₁*2³⁰+......+ b₃₁;
也就是把 a^b  =  a ^ (b₀*2³¹ + b₁*2³⁰+......+ b₃₁)  ＝ [a^{(b₀*2^31)}]  * [a^{(b₁*2^30)}] *..... * [a^{b₃₁*2^0}];
所以如果b的化成二进制的某位为0时，可以直接用这样算 resualt ＝ （resualt * resualt）%c;
若为1 则为 resualt = ( resualt * a) % c;

例如：3⁵ mod 4 ;

         5 转化成二进制为101，即3⁵ = 3^{(2^2)} *   3^{(2^0)} = 3⁴ * 3¹
所以：3⁵ mod 4 = 3⁴ * 3 mod 4 = ((3⁴)mod 4 * 3)mod 4;
3¹ mod 4 = 3, 3² mod 4 = (3¹  * 3 ¹) mod 4 ; 3⁴ mod 4 = (3² *3²)mod 4;

 

int a_b_Mod_c(int a, int b, int c) { int digit[32]; int i, k, resualt = 1; i = 0; while(b) { digit[i++] = b&0x00000001; b >>= 1; } for(k = i-1; k >= 0; k--) { resualt = (resualt * resualt) % c; if(digit[k] == 1) { resualt = (resualt * a) % c; } } return resualt; }