裴蜀定理:
存在两个正整数a,b,使得a * x + b * y = gcd( a , b )
证明方法(扩展欧几里得算法):
1、当 b = 0 时
gac( a , b ) = a
易知 x = 1,y = 0,使得a * x + b * y = 0
2、当b != 0 时
gcd( a , b ) = gcd( b , a % b);------------------------1
假设 a * x + b * y = gcd( a , b );--------------------2
也有 b * x1 + (a % b) * y1 = gcd( b , a % b );---------3
由1 2 3 可知
a * x + b * y = bx1 + (a - a/b * b)y1;
=a * y1 + b * (x1 - a/b * y1);
解得:x = y1,y = x1 - a/b*y1
由此可以递归,有 (1,0) = (xn,yn) -> (xn-1,yn-1) -> (xn-2,yn-2) -> ... ->(x,y)
由欧几里得算法可知,递归到最后b一定为0
代码实现
//acwing代码
int exgcd(int a,int b,int &x,int &y)
{
if(b==0)
{
x = 1,y = 0;
return a;
}
int d = exgcd(b,a%b,y,x);
y -= a / b * x;
return d;
}
//算法竞赛代码
int exgcd(int a,int b,int &x,int &y)
{
if(b==0)
{
x = 1,y = 0;
return a;
}
int d = exgcd(b,a%b,x,y);
int z = x;x = y;y = z - (a/b)*y;
return d;
}