C语言扩展欧几里得算法代码

更新时间：2013年09月01日 15:03:01 作者：

这篇文章介绍了扩展欧几里得算法的实现代码，有需要的朋友可以参考一下

给定两个正整数m和n,我们计算它们的最大公因子d和两个整数a和b,使得a*m+b*n=d

算法流程
　　E1.置a'=b=1;a=b'=0;c=m,d=n;

　　E2.计算d和r,使得c=q*d+r;

　　E3.若r==0;则退出,当前已有a*m+b*n=d;

　　E4;c=d;d=r;t=a';a'=a;a=t-q*a;t=b';b'=b;b=t-q*b;返回E2.

证明

　　对于已有的m和n,假设m>n;如果刨除变量a,b,a',b';算法与欧几里得算法完全一样,为计算最大公约数的算法.

　　最终要求的为a*m+b*n=d=GCD(m,n);如果改式子成立由欧几里得算法可推出a'*n+b'*(m%n)=GCD(n,m%n);

　　因为GCD(m,n)=GCD(n,m%n);

　　所以a*m+b*n=a'*n+b'*(m%n)

　　　　　　　　=a'*n+b'*(m-(m/n)*n)

　　　　　　　　=a'*n+b'*m-b'*(m/n)*n

　　　　　　　　=b'*m+(a'-b'*(m/n))*n

　　所以a=b';b=a'-b'*(m/n);

　　可以推出根据a‘、b'可以计算a、b。

代码实现

复制代码代码如下:

void EGCD(int m,int n)
{
    int a,a1,b,b1,c,d,q,r,t;
    a1=b=1,a=b1=0,c=m,d=n;
    while(1)
    {
        q=c/d,r=c%d;
        if(r==0)
        {
            printf("(%d)*%d+(%d)*%d=%d\n",a,m,b,n,d);
            return;
        }
        c=d,d=r,t=a1,a1=a,a=t-q*a,t=b1,b1=b,b=t-q*b;
    }
}

您可能感兴趣的文章:

C语言超详细介绍与实现线性表中的带头双向循环链表
带头双向循环链表：结构最复杂，一般用在单独存储数据。实际中使用的链表数据结构，都是带头双向循环链表。另外这个结构虽然结构复杂，但是使用代码实现以后会发现结构会带来很多优势，实现反而简单
2022-03-03
Embarcadero Dev-C++输出中文乱码问题图文详解
Dev-C++(或者叫做 Dev-Cpp)是Windows环境下的一个轻量级C/C++ 集成开发环境(IDE),下面这篇文章主要给大家介绍了关于Embarcadero Dev-C++输出中文乱码问题的相关资料,需要的朋友可以参考下
2023-01-01
C语言语义陷阱超详细梳理总结
这篇文章主要介绍了C语言常见的一些语义陷阱，梳理的比较全面，对我们做开发的过程中有一定帮助，感兴趣的朋友快来看看吧
2022-03-03
C++遗传算法类文件实例分析
这篇文章主要介绍了C++遗传算法的一个类文件,是学习遗传算法的绝佳参考资料,需要的朋友可以参考下
2014-08-08
C++17之std::any的具体使用
本文主要介绍了C++17之std::any的具体使用，文中通过示例代码介绍的非常详细，具有一定的参考价值，感兴趣的小伙伴们可以参考一下
2022-02-02
C语言递归思想实现汉诺塔详解
大家好，本篇文章主要讲的是C语言递归思想实现汉诺塔详解，感兴趣的同学赶快来看一看吧，对你有帮助的话记得收藏一下
2022-01-01
c语言合并两个已排序数组的示例(c语言数组排序)
如何将两个已排序数组合并成一个排序数组,下面我们给出使用c语言合并两个已排序数组的示例,需要的朋友可以参考下
2014-03-03
详解Matlab实现动态表白图的绘制
这篇文章主要利用Matlab实现绘制独特的表白动图，文中的示例代码讲解详细，对我们学习Matlab有一定的帮助，感兴趣的小伙伴可以了解一下
2022-05-05
一文详解C++中的mutable关键字
在C++中mutable关键字正如字面意思所示,表示「可变的」之意,一般在以下两种情况中使用较多,一是修饰类中的变量,用来突破const的限制,二是在Lambda表达式中使用,用来捕获修改表达式之外的变量值,下面我们就针对这两种使用场景逐个介绍
2023-10-10
浅析C++11新特性的Lambda表达式
C++11 新增了很多特性，lambda 表达式是其中之一，本文涉及到C++11这次更新中较为重要的lambda表达式。有需要的朋友们可以参考学习。
2016-08-08