朴素欧几里得:
gcd(a,b)=gcd(b,a mod b) 。
这个是比较好证明的:
假设 a=k∗b+r ,有 r=a mod b 。不妨设 d 为 a 和 b 的一个任意一个公约数,则有 a≡b≡0(mod d) 。
由于同余的性质 a−kb ≡ r ≡ 0(mod d) 因此 d 是 b 和 a mod b 的公约数。
然后 ∀ d|gcd(a,b) 都满足这个性质,所以这个定理成立。
对于方程ax+by=c,当且仅当gcd(a,b)是c的因数时有解。
ax+by=c
递归先求bu+(a%b)v=c的解
由上式可得b(u-v*(a/b))+(a%b+(a/b)*b)v=c
即av+b(u-v*(a/b))=c
所以x=v, y=u-v*(a/b)时成立
代码:
void Exgcd(int a,int b,int &x,int &y)
{
if(!b)x=1,y=0;//边界,此处是按c=1来写的
else Exgcd(b,a%b,y,x),y-=(a/b)*x;
}