#扩展欧几里得算法

2023-09-25 12 0

扩展欧几里得算法的解释:
找到一对整数(x,y)使得ax+by=gcd(a,b)

证明及过程:当a为0时,y为1,(b一样)

当a、b都不为0时,因为gcd(a,b)=gcd(b,a%b)

所以 ax+by = bx1+a%by1

= b*x1+(a-a/b * b)*y1

= ay1+b(x1-a/b*y1)

然后我们需要模拟欧几里得算法,交换x、y,并改变x

#include<stdio.h>int a,b,c,x,y;int exgcd(int a,int b,int &x,int &y){if(!b){x=1;y=0;return a;}int e=exgcd(b,a%b,x,y);int temp=x;x=y;y=temp-a/b*y;return e;}int main(){scanf("%d%d%d",&a,&b,&c);int k=exgcd(a,b,x,y);if(c%k)printf("Impossible\n");else{k=c/k;x*=k;y*=k;printf("x=%d,y=%d\n",x,y);}return 0;}

个人理解:即一个等于c的二元方程的整数解x,y,c必定是gcd(a,b)的整数倍。
如果不是,则无法做到整数x,y输出
如果是,针对解出来x,y先求出gcd(a,b)的值,然后乘c/gcd(a,b)即可

代码编程
赞赏

相关文章

pku1325 Machine Schedule
pku1915 Knight Moves
POJ的另一种登录方法http://162.105.81.212
pku1323 Game Prediction
pku1326 – Mileage Bank
zju1358 Moving Object Recognition