gcd

gcd

gcd(a,b)=gcd(b,a mod b)
边界条件:gcd(a,0)=a

gcd(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的最大公因数
}