————数学————|扩展欧几里得定理 新手向证明及代码

知识储备
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)a x + b y == g c d ( 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; }

    推荐阅读