解决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

相关文章

csu 1446 Problem J Modified LCS (扩展欧几里得算法的简单应用)

这是一道扩展欧几里得算法的简单应用题,这题是在湖南多校训练赛中队友ac的一道题,在比赛之后请教了队友,然后自己把它a掉 这也是自己独自做扩展欧几里得算法的题目 题意:把题意转变下就变成了:求d1*x - d2*y = f2 - f1的解,很明显用exgcd来解 下面介绍一下exgcd的一些知识点:求ax + by = c的解 一、首先求ax + by = gcd(a,b)的解 这个

如何解决线上平台抽佣高 线下门店客流少的痛点!

目前,许多传统零售店铺正遭遇客源下降的难题。尽管广告推广能带来一定的客流,但其费用昂贵。鉴于此,众多零售商纷纷选择加入像美团、饿了么和抖音这样的大型在线平台,但这些平台的高佣金率导致了利润的大幅缩水。在这样的市场环境下,商家之间的合作网络逐渐成为一种有效的解决方案,通过资源和客户基础的共享,实现共同的利益增长。 以最近在上海兴起的一个跨行业合作平台为例,该平台融合了环保消费积分系统,在短

科研绘图系列:R语言扩展物种堆积图(Extended Stacked Barplot)

介绍 R语言的扩展物种堆积图是一种数据可视化工具,它不仅展示了物种的堆积结果,还整合了不同样本分组之间的差异性分析结果。这种图形表示方法能够直观地比较不同物种在各个分组中的显著性差异,为研究者提供了一种有效的数据解读方式。 加载R包 knitr::opts_chunk$set(warning = F, message = F)library(tidyverse)library(phyl

【生成模型系列(初级)】嵌入(Embedding)方程——自然语言处理的数学灵魂【通俗理解】

【通俗理解】嵌入(Embedding)方程——自然语言处理的数学灵魂 关键词提炼 #嵌入方程 #自然语言处理 #词向量 #机器学习 #神经网络 #向量空间模型 #Siri #Google翻译 #AlexNet 第一节:嵌入方程的类比与核心概念【尽可能通俗】 嵌入方程可以被看作是自然语言处理中的“翻译机”,它将文本中的单词或短语转换成计算机能够理解的数学形式,即向量。 正如翻译机将一种语言

Spring框架5 - 容器的扩展功能 (ApplicationContext)

private static ApplicationContext applicationContext;static {applicationContext = new ClassPathXmlApplicationContext("bean.xml");} BeanFactory的功能扩展类ApplicationContext进行深度的分析。ApplicationConext与 BeanF

pip-tools:打造可重复、可控的 Python 开发环境,解决依赖关系,让代码更稳定

在 Python 开发中,管理依赖关系是一项繁琐且容易出错的任务。手动更新依赖版本、处理冲突、确保一致性等等,都可能让开发者感到头疼。而 pip-tools 为开发者提供了一套稳定可靠的解决方案。 什么是 pip-tools? pip-tools 是一组命令行工具,旨在简化 Python 依赖关系的管理,确保项目环境的稳定性和可重复性。它主要包含两个核心工具:pip-compile 和 pip

最大公因数:欧几里得算法

简述         求两个数字 m和n 的最大公因数,假设r是m%n的余数,只要n不等于0,就一直执行 m=n,n=r 举例 以18和12为例 m n r18 % 12 = 612 % 6 = 06 0所以最大公因数为:6 代码实现 #include<iostream>using namespace std;/

【VUE】跨域问题的概念,以及解决方法。

目录 1.跨域概念 2.解决方法 2.1 配置网络请求代理 2.2 使用@CrossOrigin 注解 2.3 通过配置文件实现跨域 2.4 添加 CorsWebFilter 来解决跨域问题 1.跨域概念 跨域问题是由于浏览器实施了同源策略,该策略要求请求的域名、协议和端口必须与提供资源的服务相同。如果不相同,则需要服务器显式地允许这种跨域请求。一般在springbo

速盾高防cdn是怎么解决网站攻击的?

速盾高防CDN是一种基于云计算技术的网络安全解决方案,可以有效地保护网站免受各种网络攻击的威胁。它通过在全球多个节点部署服务器,将网站内容缓存到这些服务器上,并通过智能路由技术将用户的请求引导到最近的服务器上,以提供更快的访问速度和更好的网络性能。 速盾高防CDN主要采用以下几种方式来解决网站攻击: 分布式拒绝服务攻击(DDoS)防护:DDoS攻击是一种常见的网络攻击手段,攻击者通过向目标网

Jenkins 插件 地址证书报错问题解决思路

问题提示摘要: SunCertPathBuilderException: unable to find valid certification path to requested target...... 网上很多的解决方式是更新站点的地址,我这里修改了一个日本的地址(清华镜像也好),其实发现是解决不了上述的报错问题的,其实,最终拉去插件的时候,会提示证书的问题,几经周折找到了其中一遍博文