C程序实现更相减损术

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

        《九章算术》是中国古代的数学专著，其中的“更相减损术”可以用来求两个数的最大公约数，原文是： 可半者半之，不可半者，副置分母、子之数，以少减多，更相减损，求其等也。以等数约之。 白话文译文： （如果需要对分数进行约分，那么）可以折半的话，就折半（也就是用2来约分）。如果不可以折半的话，那么就比较分母和分子的大小，用大数减去小数，互相减来减去，一直到减数与差相等为止，用这个相等的数字来约分。

实现步骤：

1. 任意给定两个正整数；判断它们是否都是偶数。若是，则用2约简；若不是则执行第二步。
2. 以较大的数减较小的数，接着把所得的差与较小的数比较，并以大数减小数。继续这个操作，直到 所得的减数和差相等为止。
3. 则第一步中约掉的若干个2的积与第二步中等数的乘积就是所求的最大公约数。
4. 其中所说的“等数”，就是公约数。求“等数”的办法是“更相减损”法。

比较：

        更相减损术和辗转相除法的主要区别在于前者所使用的运算是“减”，后者是“除”。从算法思想上看，两者并没有本质上的区别，但是在计算过程中，如果遇到一个数很大，另一个数比较小的情况，可能要进行很多次减法才能达到一次除法的效果，从而使得算法的时间复杂度退化为O(N)，其中N是原先的两个数中较大的一个。相比之下，辗转相除法的时间复杂度稳定于O(logN)。

实现代码：

#include <stdio.h>

int gcd(int a, int b) // 递归实现
{
    int tmp;

    if (a < b) // 如果a < b,则交换a,b的值
    {
        a = a + b;
        b = a - b;
        a = a - b;
    }

    if (a - b == b) // 以较大的数减较小的数，如果减数与差相等返回
    {
        return b; // 约掉的若干个2的积与等数(公约数)的乘积就是所求的最大公约数
    }

    if (a % 2 == 0 && b % 2 == 0)
    {
        return gcd(a / 2, b / 2) * 2; // 约掉的若干个2的积
    }
    else
    {
        return gcd(a - b, b); // 所得的差与较小的数继续计算
    }
}
int main()
{
    int a;
    int b;

    printf("请输入两个整数：");
    scanf("%d%d", &a, &b);

    printf("它们的最大公约数是:%d\n", gcd(a, b));
    printf("它们的最小公倍数是:%d\n", a * b / gcd(a, b));

    return 0;
}

#include <stdio.h>

int gcd(int a, int b) // 迭代实现
{
    int count = 1;

    if (a < b) // 如果a < b,则交换a,b的值
    {
        a = a + b;
        b = a - b;
        a = a - b;
    }

    while (a % 2 == 0 && b % 2 == 0)
    {
        a /= 2;
        b /= 2;

        count *= 2; // 约掉的若干个2的积
    }

    while (a - b != b && a - b != 0) // 直到所得的减数和差相等为止。
    {
        if (a < b)
        {
            a = a + b;
            b = a - b;
            a = a - b;
        }

        a = a - b; // 以较大的数减较小的数，接着把所得的差与较小的数比较，并以大数减小数,
    }

    return b * count; // 约掉的若干个2的积与等数(公约数)的乘积就是所求的最大公约数
}
int main()
{
    int a;
    int b;

    printf("请输入两个整数：");
    scanf("%d%d", &a, &b);

    printf("它们的最大公约数是:%d\n", gcd(a, b));
    printf("它们的最小公倍数是:%d\n", a * b / gcd(a, b));

    return 0;
}