求最大公约数欧几里得算法原理证明

  辗转相除法， 又名欧几里德算法（Euclidean algorithm），是求最大公约数的一种方法。它的具体做法是：用较大数除以较小数，再用出现的余数（第一余数）去除除数，再用出现的余数（第二余数）去除第一余数，如此反复，直到最后余数是0为止。如果是求两个数的最大公约数，那么最后的除数就是这两个数的最大公约数。

递归gcd代码

unsigned int gcd(unsigned int a, unsigned int b)
{
	return b == 0? a : gcd(b, a% b);
}

非递归代码

unsigned int gcd(unsigned int a, unsigned int b)
{
	unsigned int rem;
	while(b != 0)
	{
		rem = a % b;
		a = b;
		b = rem;
	}
	return a;
}

数学证明

$$\begin{gathered} 证明欧几里得算法，即gcd(a,b)=gcd(b,amodb),前提条件: \\ \{\begin{array}{l}a=md\\ b=nd\\ m和n互质\end{array}\Rightarrow d是a和b的最大公约数 \end{gathered}$$
$证明：$
$$\begin{gathered} 设gcd(a,b)=d_1,则\{\begin{array}{l}a=m{d}_1\\ b=n{d}_1\\ m,n互质\end{array} \\ 令a=qb+r \\ 则r=a-qb \\ =m{d}_1-qn{d}_1 \\ =(m-qn)d_1 \\ 要证明欧几里得算法，只需证明b和r的最大公因数等于d_1即可 \\ 由上述过程可得，r可以被d_1整除，即d_1是r和b的公约数。 \\ 因此只需要证明n和m-qn互质即可。 \\ 下面用反证法来证明： \\ 假设n与m-qn不互质，设gcd(n,m-qn)=d_2(d_2>1) \\ 则有\{\begin{array}{l}n=x{d}_2\\ m-qn=y{d}_2\end{array} \\ 即m=(qx+y)d_2\Rightarrow m,n不互质，这与前提m,n互质相矛盾 \\ 因此假设不成立,n和m-qn互质 \\ 所有r和b的最大公约数为d_1 \\ gcd(a,b)=gcd(b,amodb)，命题得证 \end{gathered}$$

欧几里得算法时间复杂度分析

定理： $若M>N,则(MmodN)<\frac{M}{2}。$

证明：
$$\begin{gathered} 存在两种情况。 \\ 若N\le \frac{M}{2}，则由于余数小于N，故定理成立 \\ 若N>\frac{M}{2}。但是此时M仅含有一个N，从而余数为M-N<\frac{M}{2} \\ 定理得证。 \end{gathered}$$

  由此定理可得出，欧几里得算法的迭代次数最多是$2lo{g}_{2}N=O(lo{g}_{2}N)$。