powmod快速幂法

//a:底数 n:指数  m:模 
//取模运算法则:(a * b) % p = (a % p * b % p) % p 

借助上面法则,只需要在循环乘积的每一步都进行“取模”运算,而不是等到最后直接对结果“取模”,也能达到同样的效果。

法则是理解重点!

接着传统的求幂算法之所以时间复杂度非常高,就是因为当指数n非常大的时候,需要执行的循环操作次数也非常大。所以快速幂算法的核心思想就是每一步都把指数分成两半,而相应的底数做平方运算。这样不仅能把非常大的指数给不断变小,所需要执行的循环次数也变小,而最后表示的结果却一直不会变。

long long FastPower(long long a, long long n,long long m) 
{
    long long sum = 1;
    while (n > 0) 
    {   //如果指数为偶数
        if (n % 2 == 0) 
        {
            n = n / 2;//把指数缩小为一半
            a = a * a % m;//底数平方
        } 
        //如果指数为奇数
        else
        {
            n = n - 1;//把指数减去1,使其变成偶数
            sum = sum * a % m;//此时要把指数为奇数时分离出来的底数的一次方
            n = n / 2;//此时指数为偶数,可以继续
            a = a * a % m;
        }
    }
    return sum;
}

指数缩小一半,底数平方是算法重点!

可以优化!

long long FastPower(long long a, long long n,long long m) 
{
    long long sum = 1;
    while (n > 0)
    {
        if (n % 2 == 1) //指数为奇数时
        {
            sum = sum * a % m;//分离出底数的一次方运算,使得分离出得另一半指数是偶数
                              //则可直接进行下面运算
        }
        n = n / 2;
        a = (a * a) % m;
    }
    return sum;
}

最后得优化:位运算代替(因为直接在硬件上运行  速度会更快)

#include<stdio.h>
//a:底数 n:指数  m:模 
//取模运算法则:(a * b) % p = (a % p * b % p) % p 
long long quick_mup(long long a,long long n,long long m)
{
	
    long long res=1%m,t=a;
    while(n)
    {
       if(n&1) //此处等价于if(n%2==1)  if(n&0)等价于n%2==0 
       res = res * t % m;
       t= t * t % m;
       n >>= 1;//等价于n=n/2 
    }
    return res;
}
int main()
{
    long long a,n,m;
    scanf("%lld%lld%lld",&a,&n,&m);
    printf("%lld",quick_mup(a,n,m));
    return 0;
} 


 

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
好的,关于Python实现包含快速模幂乘、大整数幂取模的RSA算,我可以给你一些相关的代码示例。首先,你需要生成一对RSA公私钥: ```python import random import math def generate_keypair(p:int, q:int): n = p * q phi = (p - 1) * (q - 1) e = 65537 d = egcd(e, phi)[1] d = d % phi return ((e, n), (d, n)) def egcd(a:int, b:int): if a == 0: return (b, 0, 1) else: g, y, x = egcd(b % a, a) return (g, x - (b // a) * y, y) p = 961748941 q = 982451653 public_key, private_key = generate_keypair(p, q) ``` 接下来,你可以使用公钥加密消息,然后使用私钥解密: ```python def encrypt(public_key, message:str): e, n = public_key cipher = [pow(ord(char), e, n) for char in message] return cipher def decrypt(private_key, cipher:list): d, n = private_key message = [chr(pow(char, d, n)) for char in cipher] return ''.join(message) message = "RSA算测试" cipher = encrypt(public_key, message) decrypted_message = decrypt(private_key, cipher) print(decrypted_message) ``` 最后,你可能需要使用快速模幂乘和大整数幂取模的函数,以避免在RSA算中执行大量的指数计算: ```python def powmod(base:int, exponent:int, modulus:int): result = 1 while exponent > 0: if exponent & 1: result = (result * base) % modulus exponent = exponent >> 1 base = (base * base) % modulus return result def modexp(base:int, exponent:int, modulus:int): if exponent == 0: return 1 elif exponent == 1: return base % modulus else: temp = modexp(base, exponent // 2, modulus) if exponent % 2 == 0: return (temp * temp) % modulus else: return ((temp * temp) % modulus * base) % modulus ``` 希望这些代码对你有所帮助!

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值