公钥算法的基本数论知识——欧几里得算法、扩展的欧几里得算法、 欧拉函数、费马小定理、欧拉定理

公钥算法的基本数论知识

包含内容

欧几里得算法、扩展的欧几里得算法、

欧拉函数、费马小定理、欧拉定理
http://www.huangjihao.com/index.php/archives/625

一、欧几里得算法(Euclidean Algorithm)

1、简介

欧几里德算法又称辗转相除法,是指用于计算两个正整数 a,b 的最大公约数。应用领域有数学和计算机两个方面。

计算公式 𝑔𝑐𝑑(𝑎,𝑏)=𝑔𝑐𝑑(𝑏,𝑎𝑚𝑜𝑑𝑏)


二、例子

𝑟0=973 , 𝑟1=301,计算它们的最大公约数

三、代码实现

Python实现:

# 欧几里得算法
def gcd(a, b):
    if b == 0:
        return a
    else:
        return gcd(b, a % b)


r0 = int(input("请输入第一个正整数:"))
r1 = int(input("请输入第二个正整数:"))
print("最大公约数是:", gcd(r0, r1))

C语言实现:

#include <stdio.h>
int gcd(int a,int b)
{
   
    if (b==0) {
   
        return a;
    }
    else{
   
        return gcd(b, a%b);
    }
}

int main()
{
   
    int r0,r1;
    int result;
    printf("请输入两个正整数:\n");
    scanf("%d%d",&r0,&r1);
    result = gcd(r0, r1);
    printf("%d\n",result);
    

}

二、扩展的欧几里得算法(Extended Euclidean algorithm)

1、简介

扩展的欧几里得算法是欧几里得算法的一个扩展。

通过扩展的欧几里得算法,我们不仅可以求出 𝑎 和 𝑏 的最大公约数,还可用找到整数 𝑥 和 𝑦,使得 𝑔𝑐𝑑(𝑎,𝑏)=𝑎𝑥+𝑏𝑦


2、原理

有两个数 𝑎,𝑏,对它们进行辗转相除法,可得它们的最大公约数。

此时欧几里得算法已经递归到了终点,𝑔𝑐𝑑(𝑎,𝑏)中的 𝑏 此时为 0。

也就是 𝑎∗1+𝑏∗0=𝑎

此时的 𝑥=1,𝑦=0

那么我们可以反推上去,求得要解的 𝑥 和 𝑦


3、代码实现

Python实现:

def ext_gcd(a, b):

    if b == 0:
        x1 = 1
        y1 = 0
        x = x1
        y = y1
        r = a
        return r, x, y
    else:
        r, x1, y1 = ext_gcd(b, a % b)
        x = y1
        y = x1 - a // b * y1
        return r, x, y


print("请输入两个正整数,且第一个数大于第二个数"<
  • 0
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
快速模幂运算: 在实现快速模幂运算时,需要注意到一个简单的公式: a ^ b % m = (a % m) ^ b % m 这个公式可以极大地减少计算量,加快计算速度。具体实现时,可以通过不断将指数b除以2,并将底数a平方,直到指数b为0,最后将结果对模数m取余即可得到快速幂结果。代码实现如下: ```python def fast_power(base, exponent, mod): result = 1 base = base % mod while exponent > 0: if exponent % 2 == 1: result = (result * base) % mod exponent = exponent // 2 base = (base * base) % mod return result ``` 扩展欧几里得: 扩展欧几里得算法是用于求解两个整数a和b的最大公约数gcd以及一组满足 ax+by=gcd 的整数解x和y。在实现该算法时,需要注意到以下性质: - gcd(a, b) = gcd(b, a % b) - 若ax+by=gcd(a,b),则bx'+(a%b)y'=gcd(b,a%b),其中 x'=y, y'=x-(a//b)y 从而可以递归调用该算法,直到b=0为止。最后返回gcd以及最后一次迭代中求得的x和y即可。代码实现如下: ```python def extended_euclidean_algorithm(a, b): if b == 0: return a, 1, 0 else: gcd, x, y = extended_euclidean_algorithm(b, a % b) return gcd, y, x - (a // b) * y ``` 中国剩余定理: 中国剩余定理是用于求解一组同余方程组的解。在实现该算法时,需要注意到以下性质: - 若m1,m2互质,则同余方程组x≡a1(mod m1),x≡a2(mod m2)的解为x=(a1m2M2+a2m1M1)(mod M),其中M=m1m2,M1=M/m1,M2=M/m2,且M1x1≡1(mod m1),M2x2≡1(mod m2) - 若m1,m2不互质,则同余方程组x≡a1(mod m1),x≡a2(mod m2)无解 可以逐个求解每个同余方程,最后通过合并得到一个更小的同余方程组来求解。具体实现时,可以使用扩展欧几里得算法求解M1和M2的乘法逆元,从而计算出x1和x2,再代入上述公式即可得到最终解。代码实现如下: ```python def chinese_remainder_theorem(a, m): M = 1 for mi in m: M *= mi result = 0 for i in range(len(a)): Mi = M // m[i] gcd, x, y = extended_euclidean_algorithm(Mi, m[i]) result += a[i] * Mi * x result %= M return result ```

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值