快速幂模n运算

模运算里的求幂运算，比如 5^596 mod 1234, 当然，直接使用暴力循环也未尝不可，在书上看到一个快速模幂算法

大概思路是，a^b mod n ,先将b转换成二进制，然后从最高位开始（最高位一定为1），如果遇到一个b[i]=0，则那么此时的结果就是b[i+1]时的结果的平方，若果b[i]=1,则结果是b[i+1]时的结果的平方再乘一个a

从b的角度理解，比如，二进制为 100 ，此时b=4，当下一位为0时，也就是 1000，即b=8，则此时的a^8=(a^4)^2 ,若果下一位为1，即二进制为 1001，b=9，则结果为 a^9=((a^4)^2)*a

代码如下：

#快速幂模n运算
 
def bits(b):
    k=[]
    while b:
        if b%2!=0:
            k.append(1)
        else:
            k.append(0)
        b>>=1
    k.reverse()      #对于List等Sequence等类型的变量，比如此处的List变量，其内置函数reverse，是直接操作变量本身，调用reverse后，变量本身的值就是reverse后的值了，所以不能出现：kk=k.reverse()这样的操作
    return k
    
def quickmod(a,b,n):      #a^b mod n
    f=1
    k=bits(b)
    for i in range(len(k)):
        f=(f*f)%n
        if k[i]:
            f=(f*a)%n
    return f