gcd(greatest common divisor)
欧几里得算法:
模板如下:
int gcd(int a,int b)
{
if(!b)return a;
return gcd(b,a%b);
}
简单写下证明:
求证gcd(a,b)=gcd(b,a%b)
令a=k*b+mod;x=gcd(a,b),y=gcd(b,a%b)=gcd(b,mod);
x是a,b的最大公约数,显然可以整除a,b,由于mod=a-k*b,x可以整除mod。
x是b,mod公约数。
y是b,mod的最大公约数,显然可以整除b,mod,进而可以整除a。
y是a,b公约数。
x整除y,y整除x,所以x=y。
复杂度为o(logmax(a,b));
看到一篇文章,有助于更为直观的理解gcd
“ 在《几何原本》第七卷的命题 2当中, Euclid给出了一种求最大公度单位的通用算法,这就是后来所说的 Euclid算法。这种方法其实非常直观。假如我们要求线段 a和线段 b的最大公度单位,不妨假设 a比 b更长。如果 b正好能度量 a,那么考虑到 b当然也能度量它自身,因而 b就是 a和 b的一个公度单位;如果 b不能度量 a,这说明 a的长度等于 b的某个整倍数,再加上一个零头。我们不妨把这个零头的长度记作 c。如果有某条线段能够同时度量 b和 c ,那么它显然也就能度量 a。也就是说,为了找到 a和 b的公度单位,我们只需要去寻找 b和 c的公度单位即可。怎样找呢?我们故技重施,看看 c是否能正好度量 b。如果 c 正好能度量 b,c就是 b和 c的公度单位,从而也就是 a和 b的公度单位;如果 c不能度量 b,那看一看 b被 c度量之后剩余的零头,把它记作 d,然后继续用 d度量 c ,并不断这样继续下去,直到某一步没有零头了为止。”————《跨越千年的RSA算法》
PS:这篇文章写得很好,推荐大家去读^O^~
http://www.matrix67.com/blog/archives/5100
扩展欧几里得(解同模方程)
=============================== _(:з」∠)_==============我是模板========================
int exgcd(int a,int b,int &x,int &y)
{
if(!b)
{
x=1;y=0;
return ;
}
exgcd(b,a%b,y,x);
y-=(a/b)*x;
}
============================== _(:з」∠)_==============模板债见========================
主要就两个式子(递归思想)
Ax1+By1=gcd(a,b);
Bx2+A%By2=gcd(b,a%b);=>一直递归到a%b=0;x=1,y=0;
stein算法
补充gcd的性质
a,b均为偶数时,gcd(a,b)=gcd(a>>1,b>>1)<<1;
a为奇数,b为偶数时,gcd(a,b)=gcd(a,b>>1);
a为偶数,b为奇数时,gcd(a,b)=gcd(a>>1,b);
a,b均为奇数时,gcd(a,b)=(abs(a-b),min(a,b));
=============================== _(:з」∠)_============我是模板=========================
int sgcd(int a,int b)
{
if(!a)return b;
if(!b)return a;
if(!(a&1)&&!(b&1))
return sgcd(a>>1,b>>1)<<1;
if(!(a&1))
return sgcd(a>>1,b);
if(!(b&1))
return sgcd(a,b>>1);
return sgcd(abs(a-b),min(a,b));
}
=============================== _(:з」∠)_==============模板债见========================