扩展欧几里德算法
内容导读
互联网集市收集整理的这篇技术教程文章主要介绍了扩展欧几里德算法,小编现在分享给大家,供广大互联网技能从业者学习和参考。文章包含796字,纯文字阅读大概需要2分钟。
内容图文
![扩展欧几里德算法](/upload/InfoBanner/zyjiaocheng/1138/f080db8bdd3243aba88fce8aa705b81e.jpg)
gcd算法:
通过辗转相除求最大公约数
#include<stdio.h> int gcd(int a,int b){ return a%b==0?b:gcd(b,a%b); } int main(){ printf("%d",gcd(15,18)); return 0; }
扩展gcd算法:
对于不完全为 0 的非负整数 a,b,若gcd(a,b)表示 a,b 的最大公约数,必然存在整数对x,y ,使得 ax+by = gcd(a,b)。
设:
ax1+by1=gcd(a,b),
bx2+a%by2=gcd(b,a%b)
因为gcd(a,b)=gcd(b,a%b),所以
ax1+by1
=bx2+a%by2
=bx2+(a-a/b*b)y2
=ay2+(x2-a/b*y2)b
所以x1=y2,y1=x2-a/b*y2
且if(b==0)不定方程 的一组解为x=1,y=0
因此扩展gcd代码为:
#include<stdio.h> #define ll long long ll exgcd(ll a,ll b,ll &x,ll &y) { if(b==0){x=1;y=0;return a;} ll r=exgcd(b,a%b,x,y); ll tmp=x;x=y;y=tmp-a/b*y; return r; } int main(){ ll a,b,ans,x,y; scanf("%lld%lld",&a,&b); ans=exgcd(a,b,x,y); printf("(%lld) * (%lld)+(%lld) * (%lld)=%lld",a,x,b,y,ans); return 0; }
原文:http://www.cnblogs.com/flipped/p/5184105.html
内容总结
以上是互联网集市为您收集整理的扩展欧几里德算法全部内容,希望文章能够帮你解决扩展欧几里德算法所遇到的程序开发问题。 如果觉得互联网集市技术教程内容还不错,欢迎将互联网集市网站推荐给程序员好友。
内容备注
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 gblab@vip.qq.com 举报,一经查实,本站将立刻删除。
内容手机端
扫描二维码推送至手机访问。