C语言——辗转相除法求最大公约数与最小公倍数

辗转相除法

  又名欧几里德算法(Euclidean algorithm)乃求两个正整数之最大公因子的算法。它是已知最古老的算法，其可追溯至公元前300年前。

来源
  设两数为 a、b(a>b),求a和b最大公约数(a，b)的步骤如下:用a除以b，得a÷b=q…r1(0≤r1)。若r=0，则(a，b)=b;若r≠0，则再用b除以r，得 b÷r=q……r₂ (0≤r₂)若r₂=0，则(a，b)=r，若r2≠0，则继续用r除以r2,如此下去，直到能整除为止。其最后一个为被除数的余数的除数即为(a,b)。
  例如:a=25,b=15，a/b=1……10,b/10=1.5,10/5=2……0,最后一个为被除数余数的除数就是5,5就是所求最大公约数。

原理
  设两数为a、b(a>b),用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)cd，b=nc=ycd,则a与b的一个公约数 cd>c,故c非a与b的最大公约数，与前面结论矛盾)，因此c也是b与r的最大公约数。
  从而可知 gcd(b,r)=c，继而gcd(a,b)=gcd(b,r)。证毕。
  以上步骤的操作是建立在刚开始时r≠0的基础之上的。即m与n亦互质。

#include <stdio.h>
int main()
{
	int m,n;
    int max,min;
    int gy(int m,int n);//函数声明 
    printf("请输入需要被求出最大公约数与最小公倍数的两个数：");
    scanf("%d,%d",&m,&n);//输入 
    max=gy(m,n);//函数调用
    min=m*n/max;//求出最小共倍数 
    printf("最大公约数：%d\n最小公倍数：%d",max,min);//输出 
    return 0;
}
int gy(int m,int n)//定义求出最大公约数的函数 
{
	int t,r,max1;
	if(m<n)//保证除数大于被除数
	{
		t=m;m=n;n=t; 
	}
	r=m%n;
	if(r==0)
	{
		max1=n;//若余数为0，则最大公约数为被除数
	}
	else
	{
		while(r!=0)
	    {
			t=r;
			r=n%r;
			n=t;
		}
		max1=t;//当余数为0时，最大公约数为被除数 
	}
	return max1;
}

例如：输入2，3

运行结果如下：