- 辗转相除法 (又称欧几里得算法):使用较大的数除以较小的数得到余数,再用较小的数除以余数,如此往复,直到最后得到的余数为0,那么此时的较大数为原始数的最大公约数。例如求24和16的最大公约数过程如下:(24,16)->(16,8)->(8,0),则8就是24和16的最大公约数。使用C++实现代码如下:
/**********************************************************
2. 函数:ComputeGCDivisor1
3. 描述:计算两个数的最大公约数
4. 参数:
5. 返回:
6. 备注:辗转相除法
**********************************************************/
int ComputeGCDivisor1(int a, int b)
{
int c;
int GCDivisor;
if (a<b)
{
c = a;
a = b;
b = a;
}
if (b==0)
{
GCDivisor = a;
}
else
{
GCDivisor = ComputeGCDivisor1(b, a % b);
}
return GCDivisor;
}
此处使用函数递归,即如果a>b,则a和b的最大公约数等于b和a%b的最大公约数。递归的终止条件是b=0,此时的a就是最大公约数。
- 更相减损术:老祖宗的智慧结晶,出自《九章算术》。译文:(如果需要对分数进行约分,那么)可以折半的话,就折半(也就是用2来约分)。如果不可以折半的话,那么就比较分母和分子的大小,用大数减去小数,互相减来减去,一直到减数与差相等为止,用这个相等的数字来约分。例如:24和16的最大公约数:(24,16)->(12,8)->(6,4)->(3,2)->(2,1)->(1,1),到最后数为1 ,但是由于中间约了三次2,所以最大公约数为1*2^3=8。如果遇到两个都是偶数不除以2,而是直接相减,则最后相等的数就是最大公约数。(24,16)->(16,8)->(8,8)。一下为不除以二的方法实现的代码:
/**********************************************************
1. 函数:ComputeGCDivisor2
2. 描述:计算两个数的最大公约数
3. 参数:
4. 返回:
5. 备注:更相减损术
**********************************************************/
int ComputeGCDivisor2(int a, int b)
{
while (a!=b)
{
if (a>b)
{
a = a - b;
}
else
{
b = b - a;
}
}
return a;
}
- 对辗转相除法的一些粗浅理解:首先是最大公约数的理解,即存在一个数c,能将a分解为x个c和将b分解为y个c,且c其中最大的一个。对于a>b,要寻找将能将a和b同时分解的数,则将a÷b=r…k,如果k为0,则b是最大公约数,如果k不为0,那么此时a已经分解为了r个b加上k,即r*b+k,此时只能将寻找b和k的最大公约数,如果b%k=0,证明b可以分解为整数个k,如果不等于0,则继续做b÷k的运算,直到余数为0。