最大公约数四种算法时间比较

一.题目分析
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);		//返回最大公约数到调用函数处
}
//穷举法（函数嵌套调用&