扩展欧几里得定理基础讲解 代码及证明
1 . 朴素欧几里得原理:gcd(a,b) == gcd(b,a % b)
2 . 负数取模:忽略符号返回绝对值就好了
3 . 模数原理:对于整数a,b必然存在整数k使得a % b == a - k * b,
且此时k == a / b向下取整
对于正整数a,b,必然存在整数(不一定是正数)x,y,
使得ax+by==gcd(x,y)
证明
(来自SDFZ-SPLI的援助)
把两边同时除以gcd(x,y),由朴素欧几里得定理可以得到恒等式,说明一定存在至少一组解使得$ax+by==gcd(x,y)
求解过程
a*x1+b*y1=gcd(a, b);
b*x2+gcd(a%b)*y2=gcd(b, a%b);
∵gcd(a, b)==gcd(b, a%b)
∴a*x1+b*y1==b*x2+(a%b)*y2
∵对于a和b存在整数k使得a % b == a - k * b
∴a%b可以拆成a-k*b,k==(a / b 向下取整)
∴a * x1 + b * y1 == b * x2 + ( a - k * b ) * y2 【把k==(a / b 向下取整)代入】
∴a * x1 + b * y1 == b * x2 + ( a - ( a / b ) * b ) * y2
即a * x1 + b * y1 == b * x2 + a * y2 - b * ( a / b ) * y2
所以有a * x1 + b * y1 == a * y2 + b * ( x2 - ( a / b ) * y2 )
观察上式可知 x1 = y2 , y1 = x2 - a / b * y2
可能上面五行是小学数学?不是很难吧
所以得出结论:x1,y1是由x2,y2得来的
那么x2,y2也一定是由x3,y3得来的,所以我们来递归吧
但递归需要终止条件:
记得在辗转相除法求最大公约数的时候,我们的递归终止条件是一方为0,那我们现在也可以参考一下,先递归gcd(a , b),当b==0时,此时a的值就是我们要的最大公约数,再根据扩展欧几里得定律a * x + b * y == gcd(a , b),将a = a,b = 0带入,可以求出x和y的值,此时x == 1,y == 0,在这里终止递归。
代码实现
清楚了原理以后代码应该很简单吧
int exGcd(int a,int b,int &x,int &y) {if(b==0){x=1;y=0;return a;}int r=exGcd(b,a%b,x,y);int t=x;x=y;y=t-a/b*y;return r; }
转载于:https://www.cnblogs.com/floatiy/p/9472356.html
扩展欧几里得定理基础讲解 代码及证明相关推荐
- 专题·扩展欧几里得定理【including 求解二元一次方程,线性同余方程
初见安~这里是基础数论专题(3)~[详见数论专栏] p.s:本文章假设你已经掌握了欧几里得算法--辗转相除法求最大公约数(gcd) 一.二元一次方程 形如的含有两个未知数且最高次数为1的方程我们称之为 ...
- 欧几里得定理与扩展欧几里得定理
欧几里的定理(辗转相除法): gcd(a,b) = gcd(b,a%b) gcd(a,b)表示a,b的最大公约数 证明: 设 a > b c = a%b a = k * b + c (k为某个整 ...
- 数学基础知识(扩展欧几里得定理)
一.欧几里得定理(辗转相除法) 性质:如果 i 能整除a,也能整除b,则 i 能整除a+b,也能整除ax+by(x,y为整数),同时a mod b = a-(a/b)*b 用以上的性质我们知道如果一个 ...
- 扩展欧几里得定理求ax + by = c 的通解
扩展欧几里得定理求ax + by = c 的通解: 前置条件: ax + by = c , gcd(a, b) = d 计算: a d \frac{a}{d} dax + b d \frac{b}{ ...
- 扩展欧几里得定理详解和运用(就不信你看不懂!)
1 :扩展欧几里得内容: 扩展欧几里德算法是用来在已知a, b求解一组x,y使得ax+by =c.(若 c%gcd(a,b)!=0)则无解 所以 我们求ax+by=c是不是可以转化为求 ax+by=k ...
- 扩展欧几里得定理的证明和代码
1.欧几里得算法,gcd(a,b)为a b(a>b)的最大公约数,则gcd(a,b) = acd(b, a%b) 利用这个定理我们可以反复对ab模下去求得a和b的最大公约数 代码如下 int G ...
- 浅谈扩展欧几里得定理(附裴蜀定理)
关于扩展欧几里得定理 众所周知,扩展欧几里得定理是用来求形如(a,b,c皆为整数)这样的方程的一组解[注,仅是一组解]的定理 它的原理比较复杂,本人学了挺久才懂了一点,这里就不谈了,扩欧的核心是它的思 ...
- 数论基础(1)扩展欧几里得定理
一.引言 扩欧在朴素欧几里得定理中扩展得到,主要用于解决什么问题? 1.求两个数的最大公约数(朴素欧也可以解决这个问题) 2.ax+by=gcd(a,b),求解这个线性不定方程的一组特解. (补充:贝 ...
- 扩展欧几里得定理 新手向证明及代码
知识储备 1 . 朴素欧几里得原理:gcd(a,b) == gcd(b,a % b) 2 . 负数取模:忽略符号返回绝对值就好了 3 . 模数原理:对于整数a,b必然存在整数k使得a % b == a ...
最新文章
- linux-awk的简单应用
- 获取文件名和路径函数
- Elasticsearch 7.0中引入的新集群协调子系统如何使用?
- python常用面试题_史上最全Python工程师常见面试题集锦,有这一份就够了
- Windows Phone开发(13):如何规范用户的输入行为 转:http://blog.csdn.net/tcjiaan/article/details/7341513...
- CodeForces - 1339C Powered Addition(思维+贪心)
- Oracle结构化查询语言(Structured Query Language)
- [RN] React Native 使用 Redux 比较详细和深刻的教程
- CCF推荐国际学术会议与学术期刊
- CMD控制台光标无法显示
- Spring Boot 学习[四] web项目实战训练(增删改查,分页,排序)
- win10打开蓝牙_用了N年也未必知!Win10竟藏有这么多实用功能
- Python爬虫报错HTTPSConnectionPool(host=‘heat.qq.com‘, port=443)处理
- java quartz 动态执行,浅谈SpringBoot集成Quartz动态定时任务
- 我们一起学linux之V4L2摄像头应用流程
- java 英文日期转数字_DatePickerDialog的月份为英文转成数字
- 手机支持电影格式一览
- sqlite:微信数据库
- 如何在七牛云测试域名失效情况下下载文件[windows]
- 正点原子Linux移植Qt,正点原子I.MX6U Qt综合例程源码