欧几里得算法 求最大公约数GCD:
辗转相除引理:
设数a,b(a>b),可得:
gcd(a,b)=gcd(b,a%b)
- 证明:
设 c=gcd(a,b)
可得 a=xc, b=yc (x,y互质)
又a%b=a-(a/b)*b
其中设q=a/b,上式变为 :a%b=a-qb=xc-q*(yc)
化简为:a%b= (x-qy) c
只需证:y与(x-qy)互质即得证
反证法:
设 y与(x-qy)不互质,即有存在 k=gcd(y,x-qy)>1
可表示为:y=km ;x-qy=kn
把y带入,可得 x-q(km)=kn
可得:x=k(n+qm)
此时我们发现 ,x,y不互质,矛盾
故得证
- 实现
int gcd(int a,int b){
return b?gcd(b,a%b):a;
}
扩展欧几里得
引理:
对a,b ,一定存在整数x,y使得 ax+by=gcd(a,b)
- 证明:
当 b=0 时,gcd(a,b)=a,此时 x=1 , y=0
当 b!=0 时,
设 ax1+by1=gcd(a,b)=gcd(b,a%b)=bx2+(a%b)y2
又因 a%b=a-a/b*b
则 ax1+by1=bx2+(a-a/b*b)y2
化简得:ax1+by1=ay2+b(x2-a/b*y2)
解得 x1=y2 , y1=x2-a/b*y2
因为当 b=0 时存在 x , y 为最后一组解
解 x , y 必然存在
得证
- 实现
void gcd(int a, int b, int &d, int &x, int &y)
{
if(!b) { d = a; x = 1; y = 0; }
else { gcd(b, a%b, d, y, x); y -= x*(a/b); }
}
x,y通解形式:
对ax+by=gcd(a,b)
设t(整数),可得:a(x+bt)+b(y-at)=gcd(a,b)不定方程通解
求x,y满足:ax+by=c
在求得 ax0+by0=gcd(a,b) 带入通解 两边乘以 c/gcd(a,b)
x=x0+bt/gcd(a,b)
y=y0-at/gcd(a,b)