欧几里得算法简介:
定理:两个整数的最大公约数等于其中较小的那个数和两数相除余数的最大公约数。最大公约数(Greatest Common Divisor)缩写为GCD。
gcd(a,b) = gcd(b,a mod b) (不妨设a>b 且r=a mod b ,r不为0)
证法一
a可以表示成a = kb + r(a,b,k,r皆为正整数,且r<b),则r = a mod b
假设d是a,b的一个公约数,记作d|a,d|b,即a和b都可以被d整除。
而r = a - kb,两边同时除以d,r/d=a/d-kb/d=m,由等式右边可知m为整数,因此d|r
因此d也是b,a mod b的公约数
假设d是b,a mod b的公约数, 则d|b,d|(a-k*b),k是一个整数。
进而d|a.因此d也是a,b的公约数
因此(a,b)和(b,a mod b)的公约数是一样的,其最大公约数也必然相等,得证。
证法二
第一步:令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)dc,b=nc=ycd,故a与b最大公约数≥cd,而非c,与前面结论矛盾】
从而可知gcd(b,r)=c,继而gcd(a,b)=gcd(b,r),得证
以上均来自百度百科介绍
程序设计:
设计思路:设定变量r记录amodb的值,使用循环结构判断r值是否为0,并记录r=0之前的余数(最大公约数)。因为循环次数不定所以采用while循环结构。
算法实现:
C语言:
#include <stdio.h>
int Gcd( int a,int b)
{
int res;//记录a mod b的值
while(b > 0)
{
res = a % b;
a = b;
b = res;
}
return a;//a为b=0所产生的结果,即最大公因子
}
int main(void)
{
int a,b;
printf("请输入两个整数\n");
scanf("%d %d",&a,&b);
printf("a和b的最大公因子为:%d,%d,%d\n",a,b,Gcd(a,b));
return 0;
}
运行结果:
D:\core>a.exe
500
300
a和b的最大公因子为:500,300,100
D:\core>