快速幂求模幂

本文实现使用快速幂算法求模幂

定义有限域$GF(p)$, 要求域内元素的幂，即计算$c=a^{b}modp$。

快速幂算法

要求$c=a^{b}modp$，首先我们将 $b$ 用二进制表示 $b=(b_n...b_1{b}_0{)}_2$。
那么，$b=b_n{2}^n+...+b_1{2}^1+b_0{2}^0$。

$c=a^{b}modp=a^{b_n{2}^n+...+b_1{2}^1+b_0{2}^0}modp$
当$b_i=0$时，$a^{b_i{2}^i}$项为无效项。
当$b_i=1$时，$a^{b_i{2}^i}$项才为有效项。
因此，只要判断b的二进制位的每一位是否为有效位即可求出幂。

实现代码

int power(int a, int b, int p)
{
    int c = 1;
    a = a % p;
    while(b > 0){
        if(b%2 == 1) 
            c = (c*a) % p;
        b = b / 2;
        a = (a*a) % p;
    }
    return c;
}