a的b次方对p取模 a乘b对p取模 快速幂

思路：对b用二进制表示，分开操作防止溢出。
对b用二进制表示：
$b=b_0\ast 2^0+b_1\ast 2^1+b_2\ast 2^2+...+b_n\ast 2^n,b_i=0or1,i=0,..,n$

$a^b=a^{b_0\ast 2^0+b_1\ast 2^1+b_2\ast 2^2+...+b_n\ast 2^n}=a^{b_0\ast 2^0}\ast a^{b_1\ast 2^1}\ast ...\ast a^{b_n\ast 2^n}$
a的b次方对p取模，每次累乘$a^{b_i\ast 2^i}$后立即取模，这样不会溢出。

$a\ast b=a\ast (b_0\ast 2^0+b_1\ast 2^1+b_2\ast 2^2+...+b_n\ast 2^n)=a\ast b_0\ast 2^0+a\ast b_1\ast 2^1+...+a\ast b_n\ast 2^n$
a乘b对p取模，每次累加$a\ast b_i\ast 2^i$后立即取模，这样不会溢出。

时间复杂度，两者都是b拆开后的项数，$O(\lfloor lo{g}_{2}b\rfloor +1)\approx O(logb)$

a = 4
b = 3
p = 3

def quick_power(a, b, p):
    tmp,sum = 1, 1
    while b:
        tmp = tmp*a%p 
        if b&1:
            sum = sum*tmp%p
        b >>= 1
    return sum
print(quick_power(a, b, p))

def quick_mul(a, b, p):
    sum, tmp = 0, 1
    flag = False
    while b:
        if not flag:
            tmp *= 1
            flag = True
        else:
            tmp = tmp*2%p
        if b&1:
            sum = (sum + a*tmp)%p
        b >>= 1
    return sum
print(quick_mul(a, b, p))