扩展欧几里得算法的解释:
找到一对整数(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)即可