解决ax+by=c,不定方程(扩展欧几里得)

2024-09-05 21:38

本文主要是介绍解决ax+by=c,不定方程(扩展欧几里得),希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

首先有几个定理我们需要知道,在这里我也会一一证明。

——————————————————————————————————————

定理1:gcd(a,b)==gcd(b,a%b);这个是欧几里得提出并证明的。 (%是取余的意思,在数学中

可用mod表示);

以下是证明过程

——————————————————————————————————————

令a = k * b + r; (k为整数);=> r = a%b;

设d是a,b的任意一个公约数。=> d|a, d|b(d|a的意思是d能被a整除);

又 r = k * b - a    =>d|(k*b - a);

以上可得d|b, d|(a%b),   =>d是b, a % b的公约数;

以上得 d既是a,b的公约数,又是b, a % b的公约数;

以上可得 gcd(a, b) == gcd(b, a%b);

证毕

——————————————————————————————————————


定理2:a*x + b*y ==gcd(a,b)一定存在解。这个定理又叫裴蜀定理,或贝祖定理。

以下会证明过程

——————————————————————————————————————

现在还不会............

——————————————————————————————————————



以下是求解 a*x +b*y  == gcd(a,b)的过程。

——————————————————————————————————————

当 b = 0时,a * x == gcd(a,0) == a;  =>x=1,y=0;

当a>b>0时:

由定理1得gcd(a,b) == gcd(b,a%b);

易得 a*X1 + b*Y1 == b * X2 + (a % b) * Y2;

=> a*X1 + b*Y1 == b * X2 +(a-[a / b] * b) * Y2;  (此处的/是不带余除法,也就是c++中的/);

=> a*X1 + b*Y1 == a * Y2 + b * (X2 - [a / b] * y2);

以上可得 1. X1 == Y2;

               2.Y1 == X2 - [a / b] * y2;

显然 以上两个方程式可以一直递归下去;

我们只要递归到b == 0的时候,就能求出Xn = 1, Yn = 0。然后我们一直往前回溯就能求出

X1,Y1;

代码如下:

#include<stdio.h>
int exgcd(int a, int b, int &x, int &y);
int main()
{int a, b, x = 0, y = 0;scanf("%d %d",&a, &b);int gcd = exgcd(a,b,x,y);printf("%d %d %d\n", gcd, x, y);return 0;
} int exgcd(int a, int b, int &x, int &y)
{if(b>a)return exgcd( b,  a,  y,  x);if(b==0){x = 1, y = 0;return a;}int r=exgcd(b, a%b, x, y);int temp = x;x = y;y = temp - (a/b) * y;return r;}

以上是求 a*x +b*y == gcd(a,b)某一组特解X1,Y1的过程

所以a*x + b*y == gcd(a,b)的通解为 X = X1 - b/gcd(a,b)*t

                                                        Y = Y1 + a/gcd(a,b)*t        t为任意整数。

——————————————————————————————————————



以下就是求ax + by = c的过程。为了好表示,我们将上一步的ax + by == gcd(a,b) 等价为

 am + bn == gcd(a,b).

——————————————————————————————————————

当c % gcd(a,b) == 0 时有解,令 k * gcd(a,b) == c;

=> k*a*m + k*b*n == k*gcd(a,b);

=>x == k*m == c*m/gcd(a,b)  ,  y == k*n == c*n/gcd(a,b) ;


设 X0 ,Y0 是 a*x + by 的某一特解。则 该不定方程的通解为

X = X0 - b/gcd(a,b)*t;

Y = Y0 + a/gcd(a,b)*t;            t为任意一个整数

X = (c*M0 - b*t)/gcd(a,b);

Y = (c*N0 + a*t)/gcd(a,b);

——————————————————————————————————————



综上求解不定方程a*x + b*y == c的步骤为

1.  先用扩展欧几里得求出 a*m + b* y ==gcd(a,b)的一组特解 M0,N0;

2.  求出a*x + b*y ==c 的通解为    X = (c*M0 - b*t)/gcd(a,b);

                                                     Y = (c*N0 + a*t)/gcd(a,b);




这篇关于解决ax+by=c,不定方程(扩展欧几里得)的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/1140143

相关文章

关于MongoDB图片URL存储异常问题以及解决

《关于MongoDB图片URL存储异常问题以及解决》:本文主要介绍关于MongoDB图片URL存储异常问题以及解决方案,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐... 目录MongoDB图片URL存储异常问题项目场景问题描述原因分析解决方案预防措施js总结MongoDB图

SpringBoot项目中报错The field screenShot exceeds its maximum permitted size of 1048576 bytes.的问题及解决

《SpringBoot项目中报错ThefieldscreenShotexceedsitsmaximumpermittedsizeof1048576bytes.的问题及解决》这篇文章... 目录项目场景问题描述原因分析解决方案总结项目场景javascript提示:项目相关背景:项目场景:基于Spring

解决Maven项目idea找不到本地仓库jar包问题以及使用mvn install:install-file

《解决Maven项目idea找不到本地仓库jar包问题以及使用mvninstall:install-file》:本文主要介绍解决Maven项目idea找不到本地仓库jar包问题以及使用mvnin... 目录Maven项目idea找不到本地仓库jar包以及使用mvn install:install-file基

最详细安装 PostgreSQL方法及常见问题解决

《最详细安装PostgreSQL方法及常见问题解决》:本文主要介绍最详细安装PostgreSQL方法及常见问题解决,介绍了在Windows系统上安装PostgreSQL及Linux系统上安装Po... 目录一、在 Windows 系统上安装 PostgreSQL1. 下载 PostgreSQL 安装包2.

Mysql如何解决死锁问题

《Mysql如何解决死锁问题》:本文主要介绍Mysql如何解决死锁问题,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录【一】mysql中锁分类和加锁情况【1】按锁的粒度分类全局锁表级锁行级锁【2】按锁的模式分类【二】加锁方式的影响因素【三】Mysql的死锁情况【1

SpringBoot内嵌Tomcat临时目录问题及解决

《SpringBoot内嵌Tomcat临时目录问题及解决》:本文主要介绍SpringBoot内嵌Tomcat临时目录问题及解决,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,... 目录SprinjavascriptgBoot内嵌Tomcat临时目录问题1.背景2.方案3.代码中配置t

如何解决idea的Module:‘:app‘platform‘android-32‘not found.问题

《如何解决idea的Module:‘:app‘platform‘android-32‘notfound.问题》:本文主要介绍如何解决idea的Module:‘:app‘platform‘andr... 目录idea的Module:‘:app‘pwww.chinasem.cnlatform‘android-32

kali linux 无法登录root的问题及解决方法

《kalilinux无法登录root的问题及解决方法》:本文主要介绍kalilinux无法登录root的问题及解决方法,本文给大家介绍的非常详细,对大家的学习或工作具有一定的参考借鉴价值,... 目录kali linux 无法登录root1、问题描述1.1、本地登录root1.2、ssh远程登录root2、

SpringBoot应用中出现的Full GC问题的场景与解决

《SpringBoot应用中出现的FullGC问题的场景与解决》这篇文章主要为大家详细介绍了SpringBoot应用中出现的FullGC问题的场景与解决方法,文中的示例代码讲解详细,感兴趣的小伙伴可... 目录Full GC的原理与触发条件原理触发条件对Spring Boot应用的影响示例代码优化建议结论F

Pyserial设置缓冲区大小失败的问题解决

《Pyserial设置缓冲区大小失败的问题解决》本文主要介绍了Pyserial设置缓冲区大小失败的问题解决,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面... 目录问题描述原因分析解决方案问题描述使用set_buffer_size()设置缓冲区大小后,buf