欧几里德算法及原理

欧几里德算法（辗转相除法）

用于计算两个整数a,b的最大公约数【gcd：Greatest Common Divisor】

原理：gcd(a,b) = gcd(b,a mod b) (条件：a>b 且a mod b 不为0)

证明①：a可以表示成a = kb + r，则r = a mod b
假设d是a,b的一个公约数，则有d|a,d|b，而r = a - kb，因此d|r
因此d也是（b,a mod b）的公约数
因此（a,b）和（b,a mod b）的公约数是一样的，其最大公约数也必然相等，得证

证明②：
第一步：令c=gcd(a,b），则设a=mc，b=nc
第二步：根据前提可知r =a-kb=mc-knc=(m-kn)c
第三步：根据第二步结果可知c也是r的因数
第四步：可以断定m-kn与n互素【否则，可设m-kn=xd，n=yd，（d>1），则m=kn+xd=kyd+xd=(ky+x)d，则a=mc=(ky+x)dc，b=nc=ycd，故a与b最大公约数≥cd，而非c，与前面结论矛盾】
从而可知gcd(b,r)=c，继而gcd(a,b)=gcd(b,r），得证

源代码：

//迭代版本 
int gcd1(int a,int b)
{
	int r;
	while(b>0)
	{
		r=a%b;
		a=b;
		b=r;
	}
	return a;
}

//递归实现 
int gcd2(int a,int b)
{
	return (b?gcd2(b,a%b):a);
}
//或者
//	int gcd2(int a,int b)
//	{
//		if(!b)             //b==0
//		{return a;}
//		return gcd2(b,a%b);
//	}

//按位运算
int gcd3(int a,int b)
{
	while(b^=a^=b^=a%=b); //b^=a^=b^=a实现a,b数值交换 
	return a;
}