辗转相除法
又名欧几里德算法(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
运行结果如下: