最大公约数
正整数a和b的最大公约数是指a与b的所有公约数中最大的那个公约数。例如4和6的最大公约数为2;3和9的最大公约数为3;一般可用gcd(a, b)来表示a和b的最大公约数。求解最大公约数常用欧几里得算法(辗转相除法)。
欧几里得算法基于下面定理:
设a、b均为正整数, 则gcd(a, b) = gcd(b, a%b)
接下来进行证明:(我写的又臭又长,参数有点多,但是很容易懂(^-^),不想深究的话,跳过证明即可)
证明:设a = kb + r,其中k和r分别为a除以b得到商和余数
则有r = a - kb成立。
设u为a和b的最大公约数;
那么a = x1 * u,b = y1 * u;
代入到r = a - kb中;
得r = x1 * u - k * y1 * u
= (x - ky)* u
因此r 整除u;
再设b,r的最大公约数为v;
则 b = x2 * v,r = y2 * v;
由a = kb + r
= k * x2 * v + y2 * v;
=( k * x2 + y2)* v;
所以a整除于v;
总结: r整除(a,b)的最大公约数,
a整除(b, r)的最大公约数
a又整除r;
所以(a,b)的最大公约数 = (b,r)的公约数
最小公倍数
很简单,直接两数相乘再除以最大公约数即可
现在直接上代码:
#include <cstdio>
//求最大公约数的辗转相除法递归写法
int gcd(int a, int b){
if(b == 0) return a;
else return gcd(b, a % b);
}
int main(){
int n, m;
printf("输入两个数");
scanf("%d%d", &n, &m);
printf("最大公约数%d\n", gcd(n, m));
printf("最小公被数%d\n", n / gcd(n, m) * m ); // 防止n * m溢出,先除再乘
return 0;
}