一.题目分析
1.辗转相除法
通过查询可知,辗转相除法(又名欧几里德法)C语言中用于计算两个正整数a,b的最大公约数和最小公倍数。它的具体做法是:用较小数除较大数,再用出现的余数(第一余数)去除除数,再用出现的余数(第二余数)去除第一余数,如此反复,直到最后余数是0为止。如果是求两个数的最大公约数,那么最后的除数就是这两个数的最大公约数。这里用函数嵌套调用法,其算法过程为:
前提:设两数为a,b设其中a 做被除数,b做除数,temp为余数
(1)大数放a中、小数放b中;
(2)求a/b的余数;
(3)若temp=0则b为最大公约数;
(4)如果temp!=0则把b的值给a、temp的值给a;
(5)返回第二步。
2.穷举法
穷举法(也叫枚举法)穷举法求两个正整数的最大公约数的解题步骤:从两个数中较小数开始由大到小列举,直到找到公约数立即中断列举,得到的公约数便是最大公约数 。对两个正整数a,b如果能在区间[a,0]或[b,0]内能找到一个整数temp能同时被a和b所整除,则temp即为最大公约数。
3.更相减损法
更相减损法(辗转相减法)是一种简便的求出两数最大公约数的方法。更相减损法(辗转相减法),即尼考曼彻斯法,其特色是做一系列减法,从而求得最大公约数。例如:两个自然数35和14,用大数减去小数(35,14)->(21,14)->(7,14),此时,7小于14,要做一次交换,把14作为被减数,即(14,7)->(7,7),再做一次相减,结果为0,这样也就求出了最大公约数7 。
更相减损法求最大公约数可分为2步:
第一步:任意给定两个正整数;判断它们是否都是偶数。若是,则用2约简;若不是则执行第二步。
第二步:以较大的数减较小的数,接着把所得的差与较小的数比较,并以大数减小数。继续这个操作,直到所得的减数和差相等为止。
则第一步中约掉的若干个2与第二步中等数的乘积就是所求的最大公约数。其中所说的“等数”,就是最大公约数。求“等数”的办法是“更相减损”法。所以更相减损法也叫等值算法。
4.Stein算法
研究最大公约数的性质,发现有 gcd( kx,ky ) = k*gcd( x,y ) 这么一个非常好的性质。试取 k=2,则有 gcd( 2x,2y ) = 2 * gcd( x,y )。很快联想到将两个偶数化小的方法,整理一下,对两个正整数 x>y :
(1)均为偶数 gcd(x,y)=2gcd(x/2,y/2);
(2)均为奇数 gcd(x,y)=gcd((x+y)/2,(x-y)/2);
(3)x奇y偶gcd( x,y )=gcd(x,y/2);
(4)x偶y奇gcd( x,y )=gcd(x/2,y)或gcd(x,y)=gcd(y,x/2);
现在已经有了递归式,还需要再找出一个退化情况。注意到gcd(x,x)=x,就用这个。
三.算法实现
程序源代码如下:
#include<stdio.h>
#include<stdlib.h>
#include<math.h> //更相减损法pow( )的头文件
#include<time.h> //计算时间所需头文件
//辗转相除法(函数嵌套调用)
int divisor1(int a,int b) //自定义函数求两数的最大公约数
{
int temp; //定义整型变量
if(a<b) //通过比较求出两个数中的最大值和最小值
{ temp=a;a=b;b=temp;} //设置中间变量进行两数交换
while(b!=0) //通过循环求两数的余数,直到余数为0
{
temp=a%b;
a=b; //变量数值交换
b=temp;
}
return (a); //返回最大公约数到调用函数处
}
//穷举法(函数嵌套调用&