P120 6.1 <谭> 输入两个正整数m和n, 求其最大公约数和最小公倍数. <1> 用辗转相除法求最大公约数算法描述:m对n求余为a, 若a不等于0则 m <- n, n <- a, 继续求余否则 n 为最大公约数<2> 最小公倍数 = 两个数的积 / 最大公约数#include <stdio.h>int main(){    int m, n;    int m_cup, n_cup, res;          /*被除数, 除数, 余数*/    printf("Enter two integer:\n");    scanf("%d %d", &m, &n);    if (m > 0 && n >0)    {        m_cup = m;        n_cup = n;        res = m_cup % n_cup;        while (res != 0)        {            m_cup = n_cup;            n_cup = res;            res = m_cup % n_cup;        }        printf("Greatest common divisor: %d\n", n_cup);        printf("Lease common multiple  : %d\n", m * n / n_cup);    }    else printf("Error!\n");    return 0;}★ 关于辗转相除法, 搜了一下, 在我国古代的《九章算术》中就有记载，现摘录如下: 约分术曰：“可半者半之，不可半者，副置分母、子之数，以少减多，更相减损，求其等也。以等数约之。” 其中所说的“等数”，就是最大公约数。求“等数”的办法是“更相减损”法，实际上就是辗转相除法。 辗转相除法求最大公约数，是一种比较好的方法，比较快。 对于52317和75569两个数，你能迅速地求出它们的最大公约数吗？一般来说你会找一找公共的使因子，这题可麻烦了，不好找，质因子大。 现在教你用辗转相除法来求最大公约数。 先用较大的75569除以52317，得商1，余数23252，再以52317除以23252，得商2，余数是5813，再用23252做被除数，5813做除数，正好除尽得商数4。这样5813就是5569和52317的最大公约数。你要是用分解使因数的办法，肯定找不到。 那么，这辗转相除法为什么能得到最大公约数呢？下面我就给大伙谈谈。 比如说有要求a、b两个整数的最大公约数，a＞b，那么我们先用a除以b，得到商8，余数r1：a÷b＝q1…r1我们当然也可以把上面这个式子改写成乘法式：a＝bq1＋r1------l） 如果r1＝0，那么b就是a、b的最大公约数3。要是r1≠0，就继续除，用b除以r1，我们也可以有和上面一样的式子： b＝r1q2＋r2-------2） 如果余数r2＝0，那么r1就是所求的最大公约数3。为什么呢？因为如果2）式变成了b＝r1q2，那么b1r1的公约数就一定是a1b的公约数。这是因为一个数能同时除尽b和r1，那么由l）式，就一定能整除a，从而也是a1b的公约数。 反过来，如果一个数d，能同时整除a1b，那么由1）式，也一定能整除r1，从而也有d是b1r1的公约数。 这样，a和b的公约数与b和r1的公约数完全一样，那么这两对的最大公约数也一定相同。那b1r1的最大公约数，在r1＝0时，不就是r1吗？所以a和b的最大公约数也是r1了。 有人会说，那r2不等于0怎么办？那当然是继续往下做，用r1除以r2，……直到余数为零为止。 在这种方法里，先做除数的，后一步就成了被除数，这就是辗转相除法名字的来历吧。 一般用辗转相除法，都列成下面的式子：                 不过，《九章》中的辗转相除法略有些不同，它叫“更相减损”，是辗转相减的方法。这也很好理解，除法就是一种连续地减去除数的一种简便运算，一直减到结果比除数小为止。比如我们用“更相减损法”来求91和49的最大公约数，可以由91减49一次，得余42；再由49减42一次，余7；更由42减7，这一回要减五次，余的还是7，再减，就是0了。那么这个7就是91和49的最大公约数。这个7就是约分术中所谓的“等数”，因为减得结果和最后一次的减数相等了，就叫等数。辗转相除法在小学中学都没教过，恐怕是有点难讲清其中的道理。不过，两千多年前的古人居然有此创造，咱们后人再学不会，可就惭愧了，何况这还是一种很实用的方法。注: 上面这些讲解均来自网络,由于忘记来源,不能标明出处,还望谅解! ######################################改进########################################经网友"私塾先生"的指点(见评论),作如下修改(还是辗转相除法):#include <stdio.h>int main(){    int m, n;    int m_temp, n_temp, res;    printf("Enter two integer:\n");    scanf("%d %d", &m, &n);    m_temp = m;    n_temp = n;    if (m > 0 && n >0)    {        do        {            res = n % m;            n = m;            m = res;        } while (m != 0);        printf("Greatest common divisor: %d\n", n);        printf("Lease common multiple  : %d\n", m_temp * n_temp / n);    }    else printf("Error!\n");    return 0;}测试结果:===================================Enter two integer:22 33↙Greatest common divisor: 11Lease common multiple  : 66=================================== ★ 换个写法,换汤不换药, 只写了最大公约数部分#include <stdio.h>int main(){    int m, n;    printf("Enter two integer:\n");    scanf("%d %d", &m, &n);     if (m > 0 && n >0)    {        int res;         while((res = m % n) != 0)        {            m = n;            n = res;        }         printf("Greatest common divisor: %d\n", n);    }    else printf("Error!\n");    return 0;}测试结果:==================================Enter two integer:22 99↙Greatest common divisor: 11==================================######################################拾遗########################################搜了一下,在一个帖中找到了另一种方法（查看该帖）,我猜这应该就是《九章》中的"更相减损" 法吧，是用减法做的，思路："如果两数不相等，那么就用大数减去小数，一直减到两数相等为止，这就是最大公约数。"#include <stdio.h>int main(){    int m, n;    printf("Enter two integer:\n");    scanf("%d %d", &m, &n);     if (m > 0 && n >0)    {        while(m != n)        {            if(m > n)                m = m - n;            else                n = n - m;        }         printf("Greatest common divisor: %d\n", m);    }    else printf("Error!\n");    return 0;}运行结果:===============================Enter two integer:22 33↙Greatest common divisor: 11===============================######################################再拾########################################又看到一个用穷举法做的( 查看出处 ), 实现如下:#include <stdio.h>int main(){    int m, n;    printf("Enter two integer:\n");    scanf("%d %d", &m, &n);     if (m > 0 && n >0)    {        int i;        int r; /* 最终为最大公约数 */         for(i = 1; i <= m || i <= n; i++)            if(m % i == 0 && n % i == 0)                r = i;         printf("Greatest common divisor: %d\n", r);    }    else printf("Error!\n");    return 0;}测试结果:==============================Enter two integer:22 99↙Greatest common divisor: 11==============================######################################三拾########################################又一种新的方法( 查看出处 ) , 至于什么原理, 是否合理, 还未仔细研究……。改写了一下，先搁这儿：#include <stdio.h>int main(){    int m, n;    printf("Enter two integer:\n");    scanf("%d %d", &m, &n);     if (m > 0 && n >0)    {        int r;              /* 最终为最大公约数 */         r = m < n ? m : n;  /* r为两个数中小者  */        while(r > 0)        {            if( !(m % r) )                if( !(n % r) )                    break;            r--;        }         printf("GCD: %d\n", r);    }    else printf("Error!\n");    return 0;}运行结果:============================Enter two integer:22 605↙GCD: 11============================其循环部分用for形式如下:        for( ; r > 0; r--)        {            if(!(m % r))                if(!(n % r))                    break;        } ★ 大概查了一下,还有种递归的算法,等看到"函数"一章时再补上吧！

评论