gcd
gcd(a,b)=gcd(b,a mod b)
边界条件:gcd(a,0)=agcd(a,b) * lcm(a,b)=a*b
于是求lcm:
a/gcd(a,b)*b;
先除保证中间结果不会溢出
//递归写法
int gcd(int a,int b)
{
return b==0?a:gcd(b,a%b);
}
//求a%b余数
int gcd(int a, int b)
{
int r;
while (b > 0)//(a%b的余数给b, b给a,最后b==0的a就是最终值)
{
r = a % b;
a = b;
b = r;
}
return a;
}
//这里的递归不会栈溢出
//gcd函数的递归层数不超过4.785lgN+1.6723
N=max(a,b);
让gcd递归最多的是gcd(Fn, F(n-1)), Fn是Fibonacci数
拓展欧几里得算法的板子
//gcd的同,求出ax+by=gcd(a,b)的一个解
__int64 extend_gcd(__int64 a, __int64 b, __int64& x, __int64& y)
{
if (a == 0 && b == 0) return -1;//无最大公约数
if (b == 0) { x = 1; y = 0; return a; }
__int64 d = extend_gcd(b, a % b, y, x);
y -= a / b * x;
return d;//得到a b的最大公因数
}