知识储备
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

扩展欧几里得定理基础讲解 代码及证明相关推荐

  1. 专题·扩展欧几里得定理【including 求解二元一次方程,线性同余方程

    初见安~这里是基础数论专题(3)~[详见数论专栏] p.s:本文章假设你已经掌握了欧几里得算法--辗转相除法求最大公约数(gcd) 一.二元一次方程 形如的含有两个未知数且最高次数为1的方程我们称之为 ...

  2. 欧几里得定理与扩展欧几里得定理

    欧几里的定理(辗转相除法): gcd(a,b) = gcd(b,a%b) gcd(a,b)表示a,b的最大公约数 证明: 设 a > b c = a%b a = k * b + c (k为某个整 ...

  3. 数学基础知识(扩展欧几里得定理)

    一.欧几里得定理(辗转相除法) 性质:如果 i 能整除a,也能整除b,则 i 能整除a+b,也能整除ax+by(x,y为整数),同时a mod b = a-(a/b)*b 用以上的性质我们知道如果一个 ...

  4. 扩展欧几里得定理求ax + by = c 的通解

    扩展欧几里得定理求ax + by = c 的通解: 前置条件: ax + by = c , gcd(a, b) = d 计算: a d \frac{a}{d} da​x + b d \frac{b}{ ...

  5. 扩展欧几里得定理详解和运用(就不信你看不懂!)

    1 :扩展欧几里得内容: 扩展欧几里德算法是用来在已知a, b求解一组x,y使得ax+by =c.(若 c%gcd(a,b)!=0)则无解 所以 我们求ax+by=c是不是可以转化为求 ax+by=k ...

  6. 扩展欧几里得定理的证明和代码

    1.欧几里得算法,gcd(a,b)为a b(a>b)的最大公约数,则gcd(a,b) = acd(b, a%b) 利用这个定理我们可以反复对ab模下去求得a和b的最大公约数 代码如下 int G ...

  7. 浅谈扩展欧几里得定理(附裴蜀定理)

    关于扩展欧几里得定理 众所周知,扩展欧几里得定理是用来求形如(a,b,c皆为整数)这样的方程的一组解[注,仅是一组解]的定理 它的原理比较复杂,本人学了挺久才懂了一点,这里就不谈了,扩欧的核心是它的思 ...

  8. 数论基础(1)扩展欧几里得定理

    一.引言 扩欧在朴素欧几里得定理中扩展得到,主要用于解决什么问题? 1.求两个数的最大公约数(朴素欧也可以解决这个问题) 2.ax+by=gcd(a,b),求解这个线性不定方程的一组特解. (补充:贝 ...

  9. 扩展欧几里得定理 新手向证明及代码

    知识储备 1 . 朴素欧几里得原理:gcd(a,b) == gcd(b,a % b) 2 . 负数取模:忽略符号返回绝对值就好了 3 . 模数原理:对于整数a,b必然存在整数k使得a % b == a ...

最新文章

  1. linux-awk的简单应用
  2. 获取文件名和路径函数
  3. Elasticsearch 7.0中引入的新集群协调子系统如何使用?
  4. python常用面试题_史上最全Python工程师常见面试题集锦,有这一份就够了
  5. Windows Phone开发(13):如何规范用户的输入行为 转:http://blog.csdn.net/tcjiaan/article/details/7341513...
  6. CodeForces - 1339C Powered Addition(思维+贪心)
  7. Oracle结构化查询语言(Structured Query Language)
  8. [RN] React Native 使用 Redux 比较详细和深刻的教程
  9. CCF推荐国际学术会议与学术期刊
  10. CMD控制台光标无法显示
  11. Spring Boot 学习[四] web项目实战训练(增删改查,分页,排序)
  12. win10打开蓝牙_用了N年也未必知!Win10竟藏有这么多实用功能
  13. Python爬虫报错HTTPSConnectionPool(host=‘heat.qq.com‘, port=443)处理
  14. java quartz 动态执行,浅谈SpringBoot集成Quartz动态定时任务
  15. 我们一起学linux之V4L2摄像头应用流程
  16. java 英文日期转数字_DatePickerDialog的月份为英文转成数字
  17. 手机支持电影格式一览
  18. sqlite:微信数据库
  19. 如何在七牛云测试域名失效情况下下载文件[windows]
  20. 正点原子Linux移植Qt,正点原子I.MX6U Qt综合例程源码

热门文章

  1. LocalDB安装和连接
  2. 5G无线帧结构与空口资源
  3. Android四大组件之Activity
  4. 好开心!Google收录了我的BLOG!
  5. 大端模式小端模式详解(不懂你打我,略略~~)
  6. 微软承认被黑客入侵,37GB 的源代码被泄露
  7. OSM地图本地发布(五)-----Geoserver发布海洋图层
  8. java读书心得_范文精选-读书心得体会-java夜未眠读书心得
  9. (PC+WAP)带手机端pbootcms模板房屋建造类网站
  10. 太有才了!街头创意涂鸦手绘图片欣赏【上篇】