标签:style blog color ar sp div c log amp
---恢复内容开始---
记a, b的最大公约数为gcd(a, b)。显然, gcd(a,b)=gcd(|a|,|b|).
计算最大公约数的Euclid算法基于下面定理:
【GCD递归定理】对于任意非负整数a和任意正整数b,gcd(a,b)=gcd(b,a%b)。
=============================================================
gcd(a,b)=gcd(b, a+kb) a,b,k为任意整数
即gcd(a,b)=gcd(b, a mod b) a≥0,b>0
Example:gcd(55,22)=gcd(22, 55mod22)=gcd(22,11)=11
证明:假定d=gcd(a,b),那么有d|a和d|b.对任何正整数b,a可表示为如下形式: a=kb+r ≡r mod b, a mod b =r , 因此,有(a mod b)= a-kb,k为某个整数。但由于d|b,b也能整除kb, 而d|a,故有d|(a mod b), 这表明d 也是b 和(amod b) 的公因子。由于这是可逆的,如果d 是b 和(a mod b) 的公因子,那么d|kb,且d|[kb+(a mod b)],这等同于d|a。这样a和b的公因子集合等同于b 和(a mod b) 的公因子集合。
==============================================================
Euclid算法最简单的递归版本如下:
1 int Euclid(int a,int b)
2 {
3 b == 0 ? return a : Enclid(b, a%b);
4 }
迭代版本:
---恢复内容结束---
最大公约数与欧几里得(Euclid)算法
标签:style blog color ar sp div c log amp
原文地址:http://www.cnblogs.com/LLGemini/p/3996588.html