拓展欧几里德

最新推荐文章于 2019-01-12 14:38:05 发布

Z_sea

最新推荐文章于 2019-01-12 14:38:05 发布

阅读量938

点赞数 1

分类专栏：扩展欧几里德数论

本文链接：https://blog.csdn.net/Z_sea/article/details/79610269

版权

本文介绍了拓展欧几里德算法，通过引例展示了如何求解两个数的最大公约数，详细阐述了解题过程，包括模拟辗转相除法、求特解的过程，以及如何找到不定方程的解。此外，还探讨了如何求解同余方程和模的逆元，最终通过代码实现和实际应用来巩固理解。

摘要由CSDN通过智能技术生成

各位大佬，转载必须注明一下博客，自己写的不容易。/流眼泪

一引例：

求两个数的gcd（a,b）=a和b两个数的最大公倍数?

短除法
更相减损法
辗转相除法

1、短除法：

百度百科——短除法

其实短除法的核心唯一分解定理。我们要求的最大公约数，其实也是它本身的一部分因子。

复杂度 $O( min\left \{ \sqrt{n} \right ,\sqrt{m}\}^{k}\ ,k=gcd(n,m) )$

注意：K=gcd(n,m)这个数的因子个数。

百度百科——唯一分解定理

我们演示一遍即可。

2、更相减损术：

百度百科——更相减损术

这个方法是来自我国数学：《九章算术》可以求最大公约数。

复杂度：O(N)

3、辗转相除法：（欧几里德）

百度百科——辗转相除法

复杂度： $log_{2}(min(n,m) )$

图像小部件

贴上代码：以上三种方法的使用。

#include<bits/stdc++.h>
using namespace std;
void DC(){
    int a=888,b=664,gcd=1;
    while(1){
        int flag=0;
        for(int i=2;i<=min(a,b);i++){
            if(a%i==0&&b%i==0){
                gcd*=i;
                a/=i,b/=i;
                flag=1;
                break;
            }
        }
        if(flag==0){    //两者为互质
            break;
        }
    }
    printf("最大公约数：gcd(664,888)=%d\n",gcd);
    printf("最小公倍数：lcm(664,888)=%d\n",gcd*a*b);
}
void XJ(){      //更相减损术
    int a=888,b=664,temp=0,cnt=0;
    while(a-b!=0){
        //printf("%d %d %d\n",a,b,temp);
        temp=a-b;
        a=b;
        b=temp;
        if(a<b) swap(a,b);
    }
    printf("%d\n",temp);
}
void GCD(){    //非递归写法
    int a&#