扩欧理解ex_gcd
来源:互联网 发布:淘宝店铺装模板 编辑:程序博客网 时间:2024/05/24 05:04
ex_gcd
一直以来都不是很懂扩欧是啥,只知道怎么用,现在来小小的写一波理解吧,毕竟完全靠记模板也是不行的。
首先扩欧是用来解决 ax+by=c 这种问题的。
然后 有这么一个定理 ax+by=gcd(a,b) 一定有解。
由欧几里得定理我们知道 gcd(a,b)=gcd(b,a%b)——>bx+(a%b)y=ax+by
下面就涉及到取模的本质了。
from Baidu
稍微想想就知道了。
带入:
bx1+(a-a/b*b)y1=ax+by
—>bx1+ay1-(a/b)*by1=ax+by
—>ay1+b(x1-(a/b)*by1)=ax+by
我们发现 x=y1; y=(x1-(a/b)*y1)
这样我们就得到了一组解。那么ex_gcd 的c++模板就出来了。
int ex_gcd(int a,int b,int &x,int &y){ if(b==0) { x=1; y=0; return a; } int d=ex_gcd(b,a%b,y,x); y-=(a/b)*x; return d; }void ex_gcd(int a,int b,int &d,int &x,int &y){ if(b==0) { x=1; y=0; d=a; } else { ex_gcd(b,a%b,d,y,x); y-=(a/b)*x; }}
两种大概差不多吧,就有时候选择用会方便一点。
对于ax+by=c ;我们求出了ax+by=gcd(a,b)
那么原式的解就是:
x=x'*(c/gcd(a,b)) y=y'*(c/gcd(a,b))
我们可以看出,当 c%gcd(a,b)!=0 时,没有整数解。
这时候我们得出一组特解,那么我们如何找到通解呢:
其实这就相当于你在特解的基础上加上或减去某个数,和不变,我们需要找到最小的情况,即范围最大的情况.
那么就有:
x=x1-b/gcd(a,b) *T (T为未知数)
y=y1-a/(gcd(a,b) *T
这时我们就求出x,y 的通解了。
阅读全文
1 0
- 扩欧理解ex_gcd
- voj1781 ex_gcd
- hdu1576 ex_gcd
- hdu 1576 ex_gcd
- ex_gcd 乘法逆元
- gcd以及ex_gcd的总结
- zoj 2609 Modular Inverse (ex_gcd)
- gcd以及ex_gcd的总结
- gcd以及ex_gcd的总结
- poj 1061 青蛙的约会 ex_gcd
- BZOJ 3122: [Sdoi2013]随机数生成器 ex_gcd+BSGS
- uva 12169 ex_gcd lrj—P316
- 紫书 例题 10-2 不爽的裁判 UVa12169 ex_gcd
- poj_1061_终于理ex_gcd正解的正确求法_结果要mod(b%gcd)
- 理解
- 理解
- 理解
- 理解
- Java学习笔记(八)——多态
- ZK Hadoop Hive 集群搭建
- (CodeForces
- 看不懂的微信小程序 正在布局移动电商之路
- Kafka剖析(一):Kafka背景及架构介绍
- 扩欧理解ex_gcd
- jQuery选择器(1)
- Android CheckBox让文字在选择框的左边
- [P1134]阶乘问题
- dp基础6--[HDU-2859] 最大对称矩阵
- sprintf用法
- 文章标题
- 记一次JVM调优-后续跟踪优化
- VS win32命令控制台程序与单片机串口数据传输