约数相关算法

1. 最大公约数

$a$ 和 $b$ 的最大公约数用 $gcd(a,b)$ 表示

最大公约数基本上有这几种算法，分别是：

1. 暴力，复杂度 $O(min(a,b))$
2. 利用 九章算法 · 更相减损术 $gcd(a,b)=gcd(a-b,b)$ ，不断缩小数据，边界条件是 $gcd(a,0)=agcd(0,b)=b$
3. 利用欧几里得发现的 $gcd(a,b)=gcd(b,a\mathrm{m}\mathrm{o}\mathrm{d}b)$ 边界条件也是一样

因为第三种数据缩小的最快，所以使用第三个，其利用第二个多来几次就可以，这里给另一种证明：

当 $a<b$ ， $a\mathrm{m}\mathrm{o}\mathrm{d}b=a$ ，所以 $gcd(a,b)=gcd(b,a)=gcd(b,a\mathrm{m}\mathrm{o}\mathrm{d}b)$

当 $a\ge b$ ，设 $a=qb+r$ ，其中 $0\le r<b$ （其实就是 $a/b$ 的商和余数），设 $g=gcd(a,b)$

$r=a\mathrm{m}\mathrm{o}\mathrm{d}b$

$\because g\mid a,g\mid b$

$\therefore g\mid qb+r,g\mid qb$

$\therefore g\mid r$

$\therefore g\mid a\mathrm{m}\mathrm{o}\mathrm{d}b$

$\because a\mathrm{m}\mathrm{o}\mathrm{d}b<b$

$\therefore g=gcd(b,a\mathrm{m}\mathrm{o}\mathrm{d}b)$

于是我们可以这样求最大公约数;

int gcd(int a, int b)
{
   
	if(b==0)return a;
	return gcd(b, a%b);
}

我们也可以用循环实现：

int gcd(int a, int b)
{
   
	if(a<b)swap(a, b);
	while(b)
	{
   
		a%=b;
		swap(a, b);
	}
	return a;
}

这是最常用的方法，当然还有一个针对高精度数据的方法

二进制算法

对于两个数 $a,b$ ，当 $a=b$ 时 ， $\gcd (a,b)=a$

否则，当 $a,b$ 都是偶数时 ，$\gcd (a,b)=2$