欧几里得定理:
公式表述
gcd(a,b)=gcd(b,a mod b)
证明:a可以表示成a = kb + r,则r = a mod b
假设d是a,b的一个
公约数,则有
d|a, d|b,而r = a - kb,因此d|r
因此d是(b,a mod b)的
公约数
假设d 是(b,a mod b)的
公约数,则
d | b , d |r ,但是a = kb +r
因此d也是(a,b)的
公约数
因此(a,b)和(b,a mod b)的公约数是一样的,其最大公约数也必然相等,得证;
拓展欧几里得定理:
其实就是:gcd(a,b)=ax+by;求解x,y的方法:
例如:gcd(225,21)=3
225=10*21+15, 21=1*15+6
15=2*6+3, 6=2*3+0
反过来:
3=15-2*6
=15-2*(21-1*15)
=3*15-2*21
=3*(225-10*21)-2*21=3*225-32*21
求解 x,y的方法的理解
设 a>b。
1,显然当 b=0,gcd(a,b)=a。此时 x=1,y=0;
2,ab<>0 时
设 ax1+by1=gcd(a,b);
bx2+(a mod b)y2=gcd(b,a mod b);
根据朴素的
欧几里德原理有 gcd(a,b)=gcd(b,a mod b);
则:ax1+by1=bx2+(a mod b)y2;
即:ax1+by1=bx2+(a-[a/b]*b)y2=ay2+bx2-[a/b]*by2;
也就是ax1+by1==ay2+b(x2-[a/b]*y2);
根据恒等定理得:x1=y2; y1=x2-[a/b]*y2;
#include <iostream>
#include <cstdio>
usingnamespace std;
intx,y,q;
voidextend_Eulid(inta,intb){
if(b == 0){
x = 1;y = 0;q = a;
}else{
extend_Eulid(b,a%b);
inttemp = x;x = y;y = temp - a/b*y;
}
}
intmain(){
inta,b;
cin>>a>>b;
extend_Eulid(a,b);
printf("%d=(%d)*%d+(%d)*%d\n",q,x,a,y,b);
return0;
}
以上内容摘自百度百科;