快速幂算法-Pow(x,n)

一、递归

原理:

        举个例子求解 x^{77} ,不需要对 x 进行 77 次乘法计算,而可以按照   x\rightarrow x^{2}\rightarrow x^{4}\rightarrow x^{9}\rightarrow x^{19}\rightarrow x^{38}\rightarrow x^{77} 的顺序,在 x\rightarrow x^{2},x^{2}\rightarrow x^{4},x^{19}\rightarrow x^{38} 这些步骤中,可以直接把上次的结果进行平方,而在 x^{4}\rightarrow x^{9},x^{9}\rightarrow x^{19},x^{19}\rightarrow x^{38} 这些步骤中,在对上次结果平方后,还需额外乘一个 x 。

        从右往左推导上述过程,利用分治思想,可以得到快速幂算法递归步骤:

        (1)求解 x^{n} ,可以先递归的计算 y= x^{n/2} ;

        (2)根据递归的计算结果,若 n 为偶数,则 x^{n}=y^{2} ;若 n 为奇数,则 x=x*y^{2}

        (3)递归出口为 n = 0,任意数的 0 次方均为 1 。

注意:当 n 为负数时,可以计算 x^{-n} 再取倒数得到结果。

double quickMul(double x,long long N) {
    if(N==0) return 1.0;
    double y=quickMul(x,N/2);
    return N%2==0?y*y : y*y*x;//讨论 N 的奇偶
}

double myPow(double x, int n) {
    long long N=n;//INT_MIN > INT_MAX 故要用long long
    return N>0? quickMul(x,N):1.0/quickMul(x,-N);//讨论 N 的正负
}

二、迭代

原理:

        由于递归需要使用额外的栈空间,故试着将递归转写为迭代。遇到一个问题,就是不知道什么时候需要额外乘一个 x 。需借助 整数的二进制拆分 :

x^{n}=x^{2^{k}}*x^{2^{k-1}}*...*x^{2^{1}}*x^{2^{0}}

从 x 开始不断进行平方,如果如果 n 的第 k 个(从右往左,从 0 开始计数)二进制位为 1,那么我们就将对应的贡献 x^{2^{k}} 计入答案。具体快速幂迭代算法参考下方代码。

double myPow(double x, int n) {
    long long N=n;
    if(N<0) N=-N;
    double res=1.0;
    double x_contribution=x;//贡献值
    while(N>0)
    {
        if(N%2==1) res*=x_contribution;//二进制最低位为 1 ,需把贡献记入结果
        x_contribution*=x_contribution;
        N=N/2;//不断右移, 舍弃 N 二进制的最低位
    }
    if(n<0) return 1.0/res;            
    else return res;
}

### 回答1: RSA-CRT 算法是用于加速 RSA 密钥的模幂运算的一种方法。这种方法通过使用 Chinese Remainder Theorem(中国剩余定理)来减少计算量,从而提高效率。 下面是 RSA-CRT 算法的 Python 代码实现: ``` def rsa_crt(m: int, d: int, p: int, q: int, dp: int, dq: int) -> int: qinv = pow(q, p - 2, p) m1 = pow(m, dp, p) m2 = pow(m, dq, q) h = (qinv * (m1 - m2)) % p return m2 + h * q ``` 在这段代码中,`m` 是要解密的信息,`d` 是 RSA 私钥,`p` 和 `q` 是 RSA 密钥的两个质因数,`dp` 和 `dq` 是 RSA 私钥的两个 CRT 参数。 快速幂算法是一种用于快速计算幂的算法,它通过递归的方式来实现。递归的关键在于,对于一个数的幂,如果幂是奇数,则可以用其平方的一半再乘上它本身来得到结果;如果幂是偶数,则可以直接将幂的一半平方来得到结果。 Robin-Miller 算法是一种用于快速计算模的算法。它的基本思想是,对于一个数的模运算,如果模数是奇数,则可以将其分解成两个较小的奇数之和,然后对这两个数分别进行模运算,再将结果相加得 ### 回答2: RSA-CRT算法是一种加密算法,其中使用了两个较慢的算法快速幂算法和Robin-Miller算法。下面是关于这两种算法的精简的Python代码实现。 1. 快速幂算法: ```python # 计算快速幂 def fast_power(base, exponent, modulus): result = 1 base = base % modulus while exponent > 0: if exponent % 2 == 1: result = (result * base) % modulus exponent = exponent // 2 base = (base * base) % modulus return result ``` 2. Robin-Miller算法: ```python # Robin-Miller素性测试 def is_prime(n, k): # 排除小于2的数和偶数 if n < 2 or n % 2 == 0: return False # 定义辅助函数 def cal_s_d(n): d = n - 1 s = 0 while d % 2 == 0: d /= 2 s += 1 return s, d # 进行k次测试 for _ in range(k): a = random.randint(2, n - 1) s, d = cal_s_d(n - 1) x = fast_power(a, d, n) if x == 1 or x == n - 1: continue for _ in range(s - 1): x = (x * x) % n if x == 1: return False if x == n - 1: break else: return False return True ``` 希望以上代码能对您有所帮助,如有任何问题,请随时告诉我。 ### 回答3: RSA-CRT算法中使用了快速幂算法和Robin-Miller算法来提高计算效率。下面是一个使用Python编写的简化代码: # 快速幂算法函数 def fast_power(base, exponent, modulus): result = 1 while exponent > 0: if exponent % 2 == 1: result = (result * base) % modulus base = (base * base) % modulus exponent = exponent // 2 return result # Robin-Miller算法函数 def robin_miller(n): # 根据Robin-Miller算法找到一个与n互质的数a a = 2 while True: if gcd(a, n) == 1: break a += 1 # 使用Robin-Miller算法进行快速判断n是否为素数 x = fast_power(a, n-1, n) if x != 1: return False else: return True # 主程序 def main(): p = 61 # 选择素数p q = 53 # 选择素数q n = p * q # 计算n euler_n = (p - 1) * (q - 1) # 计算n的欧拉函数值 # 使用Robin-Miller算法判断n是否为素数 if robin_miller(n) == False: print("n不是素数") return # 选择一个与euler_n互质的数e e = 17 # 计算e的模反元素d d = 0 while (e * d) % euler_n != 1: d += 1 # 加密 plaintext = 123 # 待加密的明文 ciphertext = fast_power(plaintext, e, n) # 加密后的密文 # 解密 decryptedtext = fast_power(ciphertext, d, n) # 解密后的明文 print("解密后的明文为:", decryptedtext) main()
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

棱角码农

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值