如何求出最大公约数

求最大公约数你有什么办法
第一种方法还是我们的暴力法，但是在两个数的数据比较大的时候这种方法就不不好用了，还是先写代码

#include <iostream>
using namespace std;
int greatest_common_divisor(int a,int b)
{
	int max=a>b?a:b;
	int min=a<b?a:b;
	if(max%min==0)
		return min;
	for(int i=min/2;i>1;i--)
	{
	  if(min%i==0&&max%i==0)
	    return i;
	}
 return 1;
}
int main()
{
cout<<greatest_common_divisor(27,18)<<endl
;
return 0;
}

这种方法的弊端是当 求 10000和10001的最大公约数时就要循环4999次
但是不得不佩服现在的计算机，结果还是秒出，改进的方法

方法2：辗转相除法
两个整数a 和 b, 他们 的最大公约数等于a除以b的余数c和b之间的最大公约数

#include <iostream>
using namespace std;
int greatest_common_divisor(int a,int b)
{
	int max=a>b?a:b;
	int min=a<b?a:b;
	if(max%min==0)
		return min;
	return greatest_common_divisor( max%min,min);
}
int main()
{
cout<<greatest_common_divisor(25,10)<<endl
;
return 0;
}

但是当a,b数据比较大的时候，取模运算也比较复杂，
第三种方法： 更相减损术
两个整数a 和 b, 他们 的最大公约数等于a减b的差值c 和较小数b之间的最大公约数

#include <iostream>
using namespace std;
int greatest_common_divisor(int a,int b)
{
	int max=a>b?a:b;
	int min=a<b?a:b;
	if(max%min==0)
		return min;
	return greatest_common_divisor( max-min,min);
}
int main()
{
cout<<greatest_common_divisor(25,10)<<endl
;
return 0;
}

上面的那种方法避免了大整数取模的麻烦，但是 如果两个数是 10000和1的话，那么就要运算9999次，
第四种方法 位运算 在更相减损术的基础上使用移位运算
当 a b 都为偶数的时候，greatest_common_divisor（a>>1,b>>1）
当 a为偶数 b 都为奇数的时候 greatest_common_divisor（a>>1,b）
当 a为奇数 b 都为偶数的时候 greatest_common_divisor（a,b>>1）
当 a b为奇数的时候，利用更相减损术，就可以变出偶数了，然后依照上述过程进行递归

#include <iostream>
using namespace std;
int greatest_common_divisor(int a,int b)
{
	if(a==b)
		return a;
	if(a&1==0&&b&1==0)
		return  greatest_common_divisor(a>>1,b>>1)<<1;//注意，求得最后的数还要再乘回去
	else if(a&1==0&&!(b&1==0))
         return  greatest_common_divisor(a>>1,b);
    else if(!(a&1==0)&&(b&1==0))
         return  greatest_common_divisor(a,b>>1);
	else
	{
	int max=a>b?a:b;
	int min=a<b?a:b;
	return greatest_common_divisor( max-min,min);
	}
	
}
int main()
{
cout<<greatest_common_divisor(18,10)<<endl;
return 0;
}