数论系列之欧几里得

最新推荐文章于 2025-03-20 10:42:24 发布

置顶 leader_win

最新推荐文章于 2025-03-20 10:42:24 发布

阅读量2.6k

点赞数 1

分类专栏： North--数论 North--整除北门的智慧——数论文章标签：数学算法

本文链接：https://blog.csdn.net/qq_27599517/article/details/50888092

版权

North--数论同时被 3 个专栏收录

41 篇文章

订阅专栏

北门的智慧——数论

39 篇文章

订阅专栏

North--整除

3 篇文章

订阅专栏

本文介绍了欧几里得算法，包括基本的辗转相除法及其扩展应用。通过算法原理和代码实现，阐述如何求解ax+by=c的问题，并讨论了当c与a、b互质时的特殊情况。文章提供了详细的解题思路和迭代过程，帮助读者理解这一重要的数论算法。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

欧几里得这部分分为欧几里得算法和扩展欧几里得算法。

欧几里得算法即我们经常说的辗转相除法，采用的基本原理是设a-qb+r，其中所有数都是整数，就有gcd(a,b)=gcd(b,r)，只用a，b表示的话就是gcd(a,b0=gcd(b,a%b);

代码如下：

<span style="font-size:14px;">long long gcd(long long a,long long b){
    if(b==0)return a;
    return gcd(b,a%b);
}</span>

而扩展欧几里得算法则是求 ax+by=c；

其中c%gcd(a,b)==0;不满足就无解

一般的，我们都让这些常量互质，即让a,b,c都除以gcd(a,b);

这时，我们可以得出两个方程：

ax+by=1;

bx1+a%by1=1;

为什么等于1？因为等于1好求，并且最后得出x乘c就好了。

采用第二个式子的原因是，如果一直迭代下去，就会得到一个值，y的系数为1；

此时令当前方程的x=1,y=1-a(当前的a)再递归回去就能解出最开始的x

方法大致是这样的

令 x=y1;

解出当前y与x1之间的关系然后就一步一步迭代，直到y的系数为1时，再回溯回来就是解了；

当然解不是唯一的 x的解系为x + b/gcd(a,b)*k k为整数。

附代码：

<span style="font-size:14px;">void _gcd(long long a,long long b,long long &x,long long &y){
    if(b==1){
        x=1;
        y=1-a;
        return;
    }
    else{
        long long x1,y1;
        _gcd(b,a%b,x1,y1);
        x=y1;
        y=x1-(a/b)*x;
    }
}
</span>