今天从算法竞赛的书上看到求这两个东西的方法,稍微记一下。
根据欧几里得算法,即著名的辗转相除法:
gcd(a,b)=gcd(b,a%b),边界条件为当b=0时,gcd(a,b)=a;
(扩展欧几里得请看另一篇博客)
int gcd(int x,int y){return y?gcd(y,x%y):x;}
那么,最小公倍数也是这么求的。
假如我们要求a,b的最小公倍数,我们先把他们拆成若干质数的积:(以下未用markdown)
a=p1e1p2e2···pnen
b=p1f1p2f2···pnfn
那么
gcd(a,b)=p1min(e1,f1) ···pnmin(en,fn)
lcm(a,b)=p1max(e1,f1) ···pnmax(en,fn)
∴gcd(a,b)*lcm(a,b)=a*b