最大公约数，欧几里得算法（gcd）

对于整数而言：
若x能够整除y，则记x|y。称x为y的约数。
若m是x除y的余数，则记作x $mod$ y=m，或x%y=m。称x被y取模的结果为m，m为模数。或称x模y等于m。$mod$就是取模运算（取余运算）。
x|y，等价于y $mod$ x=0

最大公约数

若x|y，则x是y的约数。
若x|a,且x|b，则称x是a、b的公约数。
那么，所有a、b的公约数中最大的那一个，称为a、b的最大公约数（Greatest Common Divisor），记作gcd(a,b)。
一些时候简写为g，gcd，g(a,b)，(a,b)，…

欧几里得算法（gcd）

欧几里得算法可以快速求出两数a，b的最大公约数。
其表述如下递归式：
$gcd(a,b)=\{\begin{array}{c}gcd(b,amodb)(b\ne 0)\\ a(b=0)\end{array}$

证明一下：
假设此时a>=b，否则进行一轮递归之后，a>=b。
设s|a，且s|b。
则有：$a=k_{1}s,b=k_{2}s(k_1,k2\in \mathbb{Z})$
所以：$a-b=k_{1}s-k_{2}s=(k_1-k_2)s$
即：$s|a-b$

如果a-b仍大于等于b，把a-b看做新的a，重复上面的过程，继续用b来减它。
最终得到$s|a-k\cdot b(k\in \mathbb{Z})$，使得$a-k\cdot b<b$，此时$a-k\cdot b=amodb$。
这很显然，如果a一直减去b，直到小于b，那么必定得到除以b的余数。

也就是说$s|amodb$。
这就说明任意a、b的约数都是a $mod$ b的约数。

因为任意a、b的约数都是a $mod$ b的约数，因此a、b的最大公约数也是a $mod$ b的约数。
即$gcd(a,b)=gcd(b,amodb)$

QED.

就有了代码：

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

其中!b表示b==0.

当然也有$gcd(a,b)=gcd(a,amodb)$。可是这样如果b<a，就需要手动交换，代码写起来不方便，一般不采用这种形式。

你可以证明一下，a $mod$ b <= $\frac{a}{2}$，这表明每递归两次，数据规模至少缩小一半。
若a、b与n的数据规模相当（同阶）：欧几里得算法(gcd)的时间复杂度为O(logn)。

库中的gcd函数

在实现的时候，一般不采用手写的gcd函数。事实上，比赛的编译环境支持使用自带的gcd函数：__gcd(a,b)
__gcd(a,b)包含在algorithm库内，返回两个数的最大公约数。复杂度O(logn)。

后记

于是皆大欢喜。