一.求最大公约数
辗转相除法
int gcd(int a,int b){
if(a%b==0){
return b;
}else{
return gcd(b,a%b);
}
}
二.裴蜀定理与扩展欧几里得算法
裴蜀定理:
总存在整数x,y使 ax+by=gcd(a,b) 成立。
证明方法就是辗转相除法逆过来。
裴蜀定理的推论:
1.a,b最大公约数的另一种理解方法是:找到一组整数x,y使 ax+by 得到最小正数解。
2.假如方程 ax+by=c 有整数解,则c一定是gcd(a,b)的倍数。
3.可扩展到多项式。
扩展欧几里得算法:
寻找 ax+by=gcd(a,b) 的一个解。
原理(公式推导):
设g为最大公约数,gcd调用时,第k次输入的参数为a[k],b[k];
在最后一次调用中,g=x*a[n]+y*b[n]。
(实际上最后一次x==1,y==0)
则可推导上一次调用时的x,y:g=x*b[n-1]+y*((a[n-1]-(a[n-1]/b[n-1])*b[n-1])
(a[n]对应b[n-1],b[n]对应(a[n-1]/b[n-1])*b[n-1],即a%b)
整理得:
g = y*a[n-1] + (x-y*(a[n-1]/b[n-1]))*b[n-1];
这样一直推到初始的两个a,b,就得出一个解x,y。
代码:
int a,b,x,y,g;
void exgcd(int a,int b){
if(b==0){
g=a;
x=1;
y=0;
}else{
exgcd(b,a%b);
int px=x,py=y;
x=py;
y=px-py*(a/b);
}
}