• 三种求最大公约数的方法


    一、辗转相除法

    原理:

    辗转相除法, 又名欧几里德算法(Euclidean algorithm)乃求两个正整数之最大公因子的算法。
    原理
    设两数为a、b(b<a),用gcd(a,b)表示a,b的最大公约数,r=a mod b 为a除以b以后的余数,k为a除以b的商,即a÷b=k.......r。辗转相除法即是要证明gcd(a,b)=gcd(b,r)。
    第一步:令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)。
    证毕。

    图示:

    辗转相除法:以大数除以小数,如果能整除,那么小数就是所求的最大公约数(Greatest CommonDivisor:gcd)。否则就用余数来除刚才的除数;
    再用这新除法的余数去除刚才的余数。依此类推,直到一个除法能够整除,这时作为除数的数就是所求的最大公约数。即:gcd(x,y)表示x与y的
    最大公约数,有gcd(x,y)=gcd(y,x%y),如此便可把原问题转化为求两个更小数的公约数,直到其中一个数为0,剩下的另外一个数就是两者的最
    大公约数。

    代码:

     1 #include <iostream>
     2 #include <stdio.h>
     3 #include <stdlib.h>
     4 #include <ctype.h>
     5 #include<string.h>
     6 #include <math.h>
     7 #include<algorithm>
     8 using namespace std;
     9 int gcd(int a,int b)
    10 {
    11     //不需要特判a与b的大小,如果a小于b,那么递归一次gcd两者就会交换值
    12     return b?gcd(b,a%b):a;
    13 }
    14 int main()
    15 {
    16     printf("%d
    ",gcd(2,4));
    17 }

    二、更相减损法

    原理:

    更相减损术是出自《九章算术》的一种求最大公约数的算法,它原本是为约分而设计的,但它适用于任何需要求最大公约数的场合。

    《九章算术》是中国古代的数学专著,其中的“更相减损术”可以用来求两个数的最大公约数,即“可半者半之,不可半者,副置分母、子之数,以少减多,更相减损,求其等也。以等数约之。”
    翻译成现代语言如下:
    第一步:任意给定两个正整数;判断它们是否都是偶数。若是,则用2约简;若不是则执行第二步。
    第二步:以较大的数减较小的数,接着把所得的差与较小的数比较,并以大数减小数。继续这个操作,直到所得的减数和差相等为止。
    则第一步中约掉的若干个2与第二步中等数的乘积就是所求的最大公约数。
    其中所说的“等数”,就是最大公约数。求“等数”的办法是“更相减损”法。

    图示:

    此算法简单来说是:将a和b减, 得到的差值赋值给c,并判断是否为0,如果不为0,则将b赋值给a,将c赋值给a, 然后重复上述步骤直到c为0。a-b=c,算法直观来看就是这个数学算式的左移操作,和辗转相除法实现方式很是接近,应该本源思路是一样的。

    代码:

     1 #include <iostream>
     2 #include <stdio.h>
     3 #include <stdlib.h>
     4 #include <ctype.h>
     5 #include<string.h>
     6 #include <math.h>
     7 #include<algorithm>
     8 using namespace std;
     9 int gcd(int a,int b)
    10 {
    11     int c = 0;
    12     if(a<b) swap(a,b);
    13     c = b;
    14     do
    15     {
    16         b = c;
    17         a>b ? c = a - b : c = b - a;
    18         a = b;
    19     } while (c != 0);
    20     return b;
    21 }
    22 int main()
    23 {
    24     printf("%d
    ",gcd(2,4));
    25 }

    三、穷举法

    找到两个数中那个小的,然后从它开始for循环,如果有满足的就输出然后break,否则就让这个值减去一再进行判断

    代码:

     1 #include <iostream>
     2 #include <stdio.h>
     3 #include <stdlib.h>
     4 #include <ctype.h>
     5 #include<string.h>
     6 #include <math.h>
     7 #include<algorithm>
     8 using namespace std;
     9 int gcd(int a,int b)
    10 {
    11     for(int i=min(a,b);i>=1;--i)
    12     {
    13         if(a%i==0 && b%i==0) return i;
    14     }
    15 }
    16 int main()
    17 {
    18     printf("%d
    ",gcd(2,4));
    19 }
  • 相关阅读:
    codeforces round #234B(DIV2) B Inna and New Matrix of Candies
    关于禁止ipad的home键解决方法
    ios cocos2d 画线出现闪烁问题
    ios 关于[xxx timeIntervalSinceNow]出现EXC_BAD_ACCESS错误的解决办法
    Codeforces Round #228 (Div. 2) B. Fox and Cross
    Codeforces Round #228 (Div. 2) A. Fox and Number Game
    c语言 %p
    xcode 编译opencv ios容易出现的错误
    ios编译ASIHTTPRequest时出现 'libxml/HTMLparser.h' file not found in ASIHTTPRequest
    当编译CCBReader时出现 “ CCBAnimationManager.m Use of undeclared identifier 'other‘ ” 解决方法
  • 原文地址:https://www.cnblogs.com/kongbursi-2292702937/p/11977981.html
Copyright © 2020-2023  润新知