辗除法求最大公约数和最小公倍数

        最近两天因为忙于工作的事，学习进展较慢。

        今天看了一个辗除法求最大公约数的例子，上网查了一下其数学原理，算是理解了这个算法的数学依据。

        对于两个正整数num1，num2，且num1>num2，如果num1除以num2余数为resd，若resd=0则num2即是两数的最大公约数，若resd≠0，则令num1=num2，num2=resd，再重复进行以上过程，直到resd=0，则此时num2的值即是所求两数的最大公约数。

证明过程如下：

     用gcd表示两数的最大公约数，如a=gcd(m，n)表示m和n的最大公约数是a。

     num1>num2，若num1=a*num2+b，则只需要证明gcd(num1，num2)=gcd(num2,b)即可。

     设m=gcd(num1，num2)， n=gcd(num2,b)，则m可被num1和num2整除，则m可被a*num2整除，可推知m可被num1-a*num2整除，而b=num1-a*num2，则m可被b整除，所以m是num2和b的公约数(不一定是最大公约数)，则m≤n；同理可证：n是num1和num2的公约数，n≤m；由于m≤n和n≤m同时成立，所以m=n。

参考有关资料上的编程实例，编写的C语言程序如下：

#include<stdio.h>
int main()
{
   int num1, num2, a, b, temp;
   printf("Please input two numbers:\n");
   scanf("%d,%d", &num1, &num2);  //输入两个正整数
   if(num1 < num2)          //保证num1>num2，因为下面是用num1作被除数，num2作除数
    {
       temp = num1;
       num1 = num2;
       num2 = temp;
    }
   a =num1;
   b =num2;
   while(b != 0)  //辗转相除，直到余数为0
    {
       temp = a % b;
       a = b;
       b = temp;
    }
   printf("GCD:%d\n", a);
   printf("LCM:%d\n", num1 * num2 / a); //求最小公倍数，两个数的乘积等于它们最大公约数与最小公倍数的积
   return 0;
}

        今天在编译运行的时候出了点差错，输入两个数后得到的结果不正确，得到的是比较奇怪的数，我就把两个正整数不用输入改为自己赋初值，此时得到的结果正确，由此可见应该是输入的时候出错了。经过检查发现scanf("%d,%d",&num1, &num2); 这句输入程序中两个%d之间有逗号，而我从键盘输入时没有输这个逗号，导致了出错，要么把逗号去除，要么输数的时候两数之间加个逗号就正确了。这个错误可以给初学编程的人作为参考，避免犯类似的低级错误。