利用辗转相除法求两个自然数的最大公因数 程序如下：
//利用辗转相除法求两个自然数的最大公因数
int gcd(int a, int b)
{
int r;
while(b)
{
r = a%b;
a = b;
b = r;
}
return a;
}
相关理论如下：
「辗转相除法」又叫做「欧几里得算法」，是公元前 300 年左右的希腊数学家欧几里得在他的著作《几何原本》提出的。利用这个方法,可以较快地求出两个自然数的最大公因数，即 HCF 或叫做 gcd。所谓最大公因数，是指几个数的共有的因数之中最大的一个，例如 8 和 12 的最大公因数是 4，记作 gcd(8,12)=4。
在介绍这个方法之前,先说明整除性的一些特点,注以下文的所有数都是正整数,以后不再重覆.
我们可以这样给出整除以的定义:
对於两个自然数 a 和 b，若存在正整数 q，使得 a=bq，则 b 能整除 a，记作 b | a，我们叫 b 是 a 的因数，而 a 是 b 的倍数。那麼如果 c | a，而且 c | b，则 c 是 a 和 b 的公因数。
由此,我们可以得出以下一些推论：
推论一：如果 a | b，若 k 是整数，则 a | kb。因为由 a | b 可知 ha=b，所以 (hk)a=kb，即 a | kb.
推论二：如果 a | b 以及 a | c，则 a | (b±c)。因为由 a | b 以及 a | c，可知 ha=b，ka=c，二式相加,得 (h+k)a=b+c,即 a | (b+c).同样把二式相减可得 a | (b-c)。
推论三：如果 a | b 以及 b | a，则 a=b。因为由 a | b 以及 b | a，可知 ha=b，a=kb，因此 a=k(ha)，hk=1，由於 h 和 k 都是正整数，故 h=k=1，因此 a=b。
辗转相除法是用来计算两个数的最大公因数，在数值很大时尤其有用而且应用在电脑程式上也十分简单。其理论如下：
如果 q 和 r 是 m 除以 n 的商及余数，即 m=nq+r，则 gcd(m,n)=gcd(n,r)。 证明是这样的:
设 a=gcd(m,n)，b=gcd(n,r)
则由 a | m 及 a | n，可得 a | (m-nq)(由推论一及推论二得出的)，即 a | r ，又 a | n，所以 a | b。
由 b | r 及 b | n，可得 b | (nq+r)，即 b | m，又 b | n，所以b | a。
因为 a | b 并且 b | a，所以 a=b，即 gcd(m,n)=gcd(n,r)。
举例计算 gcd(546, 429),由於 546=1(429)+117,429=3(117)+78,117=1(78)+39,78=2(39),因此
gcd(546, 429)
=gcd(429, 117)
=gcd(117, 78)
=gcd(78, 39)
=39
此处再添加一个程序例子，不过不是利用辗转相除法
求最大公约数和最小公倍数:
#include <iostream>
using namespace std;
int gec,lcm;
void process(int x,int y)
{
gcd=x<y?x:y;
lcm=x<y?y:x;
for(gcd=x<y?x:y;gcd>1;gcd--)if(x%gcd==0&&y%gcd==0)break;
for(lcm=x<y?x:y;lcm>1;lcm++)if(lcm%x==0&&lcm%y==0)break;
return;
}