1.来源
设两数为a、b(a>b),求a和b最大公约数(a,b)的步骤如下:用a除以b,得a÷b=q......r1(0≤r1)。若r1=0,则(a,b)=b;若r1≠0,则再用b除以r1,得b÷r1=q......r2 (0≤r2).若r2=0,则(a,b)=r1,若r2≠0,则继续用r1除以r2,……如此下去,直到能整除为止。其最后一个为被除数的余数的除数即为(a, b)。
例如:a=25,b=15,a/b=1......10,b/10=1......5,10/5=2.......0,最后一个为被除数余数的除数就是5,5就是所求最大公约数。
2.原理
3.算法
自然语言描述
用辗转相除法确定两个正整数 a 和 b(a≥b) 的最大公因数gcd(a,b):
当a mod b=0 时gcd(a,b)=b,否则
gcd(a,b) = gcd(b,a mod b)
递归或循环运算得出结果
/** * * @return int * @tags @param m * @tags @param n * @tags @return * @todo 【方法二】利用辗除法 */ public static int gcd(int m, int n) { while (true) { if ((m = m % n) == 0) return n; if ((n = n % m) == 0) return m; } }
4.时间复杂度
辗转相除法的运算速度为 O(n),其中 n 为输入数值的位数。
辗转相除法处理大数时非常高效,它需要的步骤不会超过较小数的位数(十进制下)的五倍。加百利·拉梅(GabrielLamé)于1844年证明了这点,开创了计算复杂性理论。
5.应用
求不定方程的一组整数解方法
[注:以下出现的qi,ri括号中的是下标,gcd(a,b)为a,b的最大公约数]
辗转相除法可以求出特定条件的不定方程的一组整数解。
设不定方程为ax+by=c,其中a,b,c为整数,且 gcd(a,b) | c
a,b辗转相除的算式为
b=q1 a+r2
a=q2 r2+r3
r2=q3 r3+r4
...
rn-2=qn-1rn-1+rn
rn-1=qnrn
其中rn=gcd(a,b),不妨令b=r0,a=r1,rn+1=0
第i个算式为
ri-1= qi×ri+ ri+1
所以ri+1= ri-1 - qi×r(i)......(1)
用公式(1)可以得到rn=gcd(a,b)关于a,b的线性组合sa+tb=gcd(a,b)
举例说明
例如不定方程为326x+78y=4,求出一组整数解x,y
求(326,78)的算式为:
326=4*78+14
14=326-4*78
78=5*14+8
8=78-5*14
14=1*8+6
6=14-1*8
8=1*6+2
2=8-1*6
6=3*2
所以
2=8-6=8-(14-8)
=2*8-14=2*(78-5*14)-14
=2*78-11*14=2*78-11*(326-4*78)
=46*78-11*326
即2=(-11)*326+46*78
所以4=(-22)*326+92*78
所以x = - 22, y = 92是不定方程326x+78y=4的一组解。
6.相关原理
两个整数的最大公约数是能够同时整除它们的最大的正整数。
辗转相除法基于如下原理:两个整数的最大公约数等于其中较小的数和两数的相除余数的最大公约数。
例如,252和105的最大公约数是21(252 = 21 × 12;105 = 21 × 5);
因为252 ÷105 = 2......42,所以(105,42)是21。在这个过程中,较大的数缩小了,所以继续进行同样的计算可以不断缩小这两个数直至余数变为零。这时的除数就是所求的两个数的最大公约数。由辗转相除法也可以推出,两数的最大公约数可以用两数的整数倍相加来表示,如21 = 5 × 105 + (−2) × 252。这个重要的等式叫做贝祖等式(又称“裴蜀定理”)。