【数据结构与算法】最大公约数与最小公倍数

最大公因数（英语：highest common factor，hcf）也称最大公约数（英语：greatest common divisor，gcd）是数学词汇，指能够整除多个非零整数的最大正整数。例如8和12的最大公因数为4。表示：gcd(8,12) = 4

最小公倍数（英语：least common multiple，lcm）是数论中的一个概念。若有一个数X，可以被另外两个数A,B整除，且X同时大于或等于A和B，则X为A和B的公倍数。A和B的公倍数有无限个，而所有正的公倍数中，最小的公倍数就叫做最小公倍数。例如8和12的最小公倍数为24。表示：lcm(8,12) = 24

两个整数的最小公倍数与最大公因数之间有如下的关系：

公式

根据上面的示例：|8*12|/4=24

也就是说，我们掌握了最大公约数的求法，也就能求最小公倍数

那么我们就以最大公约数为例

枚举(穷举)

根据已有的数学知识我们知道，m>n的话，m,n的最大公约数永远不可能是m（因为n/m不能整除），最大只可能是n。所以我们只需要从n开始依次向较小数找到“能同时被n,m整除的第一个整数”即为两数的最大公约数。

代码示例：

#include <stdio.h>
int main()
{int a, b;scanf("%d %d", &a, &b);int i;int gcd;gcd = 1;for(i = (a<b?a: b); i > 0; i--){if(a % i == 0 && b % i == 0){gcd = i;break;}}printf("gcd = %d\n", gcd);return 0;
}

另一种穷举原理：

求出两数的所有公因子，再把公因子累乘得到最大公约数。

代码示例：

#include <iostream>
using namespace std;
int CommFactor2(int m, int n); 
int main()
{int a, b;cin >> a >> b;cout << "这两个数的最大公约数为：" << CommFactor2(a,b)<< endl;return 0;
}
int CommFactor2(int m,int n)
{int i;int factor = 1;for (i=2;i<=m&&i<<n;i++){while(m % i == 0 && n % i == 0)  //这里不能用if语句，因为可能会有重复的公因子{factor = factor * i;m = m / i;n = n / i;}}return factor;
}

辗转相除法

辗转相除法，又称欧几里得算法。辗转相除法基于如下原理：两个整数的最大公约数等于其中较小的数和两数相除余数的最大公约数。表示：gcd(a,b) = gcd(b,a mod b) (a>b)

辗转相除法

代码示例：

#include <stdio.h>
int main(int argc, const char* argv[])
{int a,b,temp;scanf("%d %d", &a, &b);while (b != 0){temp = a % b;a = b;b = temp;}printf("gcd = %d\n", a);return 0;
}

更相减损法

更相减损法，又称辗转相减法。更相减损法出自《九章算术》：“可半者半之，不可半者，副置分母、子之数，以少减多，更相减损，求其等也。以等数约之。”(原本是为了约分而设计的)

具体方法为两个数之间大的数字减小的数字，之后将得到的差作为减数，较小的数作为被减数，再次相减，直到与所得的差相同，此时的差即为两个数之间的最大公约数。

用(a，b)表示a和b的最大公因数：有结论(a，b)=(a，k*a+b)，其中a、b、k都为自然数。

基于上面的原理，就能实现我们的迭代相减法：(78，14)=(64，14)=(50，14)=(36，14)=(22，14)=(8，14)=(8，6)=(2，6)=(2，4)=(2，2)=(0，2)=2

代码示例：

#include<iostream>
using namespace std;
int gcd(int a, int b)
{while(a != b)if(a > b) a -= b;else b -= a;return a;
}
int main(){int m,n;cin>>m>>n;cout<<"gcd:"<<gcd(m,n)<<endl;return 0;
}

Stein算法

Stein算法是针对欧几里德算法在对大整数进行运算时，需要试商导致增加运算时间的缺陷而提出的改进算法。

欧几里德算法是计算两个数最大公约数的传统算法，无论从理论还是从实际效率上都是很好的。但是却有一个致命的缺陷，这个缺陷在素数比较小的时候一般是感觉不到的，只有在大素数时才会显现出来：一般实际应用中的整数很少会超过64位（当然现在已经允许128位了），对于这样的整数，计算两个数之间的模是很简单的。对于字长为32位的平台，计算两个不超过32位的整数的模，只需要一个指令周期，而计算64位以下的整数模，也不过几个周期而已。但是对于更大的素数，这样的计算过程就不得不由用户来设计，为了计算两个超过64位的整数的模，用户也许不得不采用类似于多位数除法手算过程中的试商法，这个过程不但复杂，而且消耗了很多CPU时间。对于现代密码算法，要求计算128位以上的素数的情况比比皆是，比如说RSA加密算法至少要求500bit密钥长度，设计这样的程序迫切希望能够抛弃除法和取模。

Stein算法很好的解决了欧几里德算法中的这个缺陷，Stein算法只有整数的移位和加减法。下面就来说一下Stein算法的原理：

  1. 若a和b都是偶数，则记录下公约数2，然后都除2（即右移1位）
  2. 若其中一个数是偶数，则偶数除2，因为此时2不可能是这两个数的公约数了
  3. 若两个都是奇数，则a = |a-b|，b = min(a,b)，因为若d是a和b的公约数，那么d也是|a-b|和min(a,b)的公约数
这里对3.进行简单的证明：
不妨设奇数A>B，A和B的公约数为X，即A=jX，B=kX，其中j，k均为正整数且j>k。
A−B=(j−k)X
因为j，k均为整数，所以X也是A-B的公约数。
min(A,B)=B
所以A-B与min(A,B)公约数相同，因为A，B都是奇数，所以A-B必然是偶数，偶数又可以二除移位了。

代码示例1：

#include<iostream>
using namespace std;
int SteinGCD(int a, int b) {if (a < b) { int t = a; a = b; b = t; }if (b == 0) return a;if ((a & 1) == 0 && (b & 1) == 0)return SteinGCD(a >> 1, b >> 1) << 1;else if ((a & 1) == 0 && (b & 1) != 0)return SteinGCD(a >> 1, b);else if ((a & 1) != 0 && (b & 1) == 0)return SteinGCD(a, b >> 1);elsereturn SteinGCD(a - b, b);
}
int main()
{int m,n;cin>>m>>n;cout<<"gcd:"<<SteinGCD(m,n)<<endl;return 0;
}

代码示例2：

#include<iostream>
using namespace std;
int gcd(int u, int v)
{if (u == 0)	return v;if (v == 0)	return u;if (~u & 1){if (v & 1)return gcd(u >> 1, v);elsereturn gcd(u >> 1, v >> 1) << 1;}if (~v & 1)return gcd(u, v >> 1);if (u > v)return gcd((u - v) >> 1, v);return gcd((v - u) >> 1, u);
}
int main()
{int m,n;cin>>m>>n;cout<<"gcd:"<<gcd(m,n)<<endl;return 0;
}

Stein算法非递归方式

代码示例：

#include<iostream>
using namespace std;
int SteinGCD(int a, int b) {int acc = 0;while ((a & 1) == 0 && (b & 1) == 0) {acc++;a >>= 1;b >>= 1;}while ((a & 1) == 0) a >>= 1;while ((b & 1) == 0) b >>= 1;if (a < b) { int t = a; a = b; b = t; }while ((a = (a - b) >> 1) != 0) {while ((a & 1) == 0) a >>= 1;if (a < b) { int t = a; a = b; b = t; }}return b << acc;
}
int main()
{int m,n;cin>>m>>n;cout<<"gcd:"<<SteinGCD(m,n)<<endl;return 0;
}

库函数

__gcd(a,b)——库algorithm

注意：gcd前面有两个下划线

gcd(a,b)——库numeric

注意：仅作扩展，如无必要，不要使用

位运算

利用位运算的特性，将两数交换改成位运算。
inline可加可不加。我实际试验中，在1e8次执行后，加与不加的时间差在80ms左右，而两者本来的运行时间均在3000ms上下，即差别不大

代码示例：

inline int gcd(int a, int b)
{while(b ^= a ^= b ^= a %= b);    return a;
}

利用取模特点

很快，几乎与__gcd(a,b)的时间一致

代码示例：

int gcd(int a, int b){if(!a || !b)return max(a, b);for(int t; t = a % b; a = b, b = t);return b;
}

以上方法都为最大公约数的求法，若要求最小公倍数，利用关系计算即可。

迭乘法(求最小公倍数)

由公倍数的定义出发，如果一个数k是a和b的公倍数，那么k可以表示成 a*m或 b*n，而当(a*m) % b == 0时，k是最小公倍数

代码示例：

#include<stdio.h>
int main(){int a,b;scanf("%d %d", &a, &b);int i = 1;while((a*i) % b != 0){i++;}printf("最小公倍数是：%d\n", a * i);return 0;
}

参考博文：

https://blog.csdn.net/Hsuesh/article/details/111992593

https://blog.csdn.net/ly_6699/article/details/90719315

https://blog.csdn.net/JH13thpig/article/details/124362053

https://blog.csdn.net/weq2011/article/details/127953257

https://blog.csdn.net/Holmofy/article/details/76401074

https://blog.csdn.net/wyd_333/article/details/126111037

本文来自互联网用户投稿，该文观点仅代表作者本人，不代表本站立场。本站仅提供信息存储空间服务，不拥有所有权，不承担相关法律责任。如若转载，请注明出处：http://www.hqwc.cn/news/623479.html

如若内容造成侵权/违法违规/事实不符，请联系编程知识网进行投诉反馈email:809451989@qq.com，一经查实，立即删除！

相关文章

2.4G漂移小车电子方案 酷得智能科技

2.4G漂移小车电子方案 酷得智能科技

漂移高速遥控车是一种专门设计用于执行高速漂移动作的遥控车模型。以下是一些关于漂移高速遥控车的功能介绍&#xff1a; 1、高速性能&#xff1a;漂移车通常配备有强力的电机和电池&#xff0c;以便在保持高速的同时进行漂移动作。 2、漂移能力&#xff1a;漂移车的轮胎和悬挂…
阅读更多...
springboot同时支持jsp+vue页面启动

springboot同时支持jsp+vue页面启动

1、参考文档链接 参考上面文档边百度边改&#xff0c;现在可以了&#xff0c;分享下 2、Java项目目录结构 3、pom.xml内容 <?xml version"1.0" encoding"UTF-8"?><project xmlns"http://maven.apache.org/POM/4.0.0" xmlns:xsi&quo…
阅读更多...
dbever可视化工具使用

dbever可视化工具使用

分别遵循 1 2 3步骤 可以查看除了执行后的控制台以外的其他窗口。
阅读更多...
波分设备的功能和应用

波分设备的功能和应用

波分设备主要利用光波在光纤中传输数据&#xff0c;具有以下功能&#xff1a; 高速数:据传输。波分技术能提供比传统电信技术更高的数据传输速度。 大带宽:通过将不同频率的光波分配到不同的频段中&#xff0c;波分技术使得更多数据可以同时传输。低衰减&#xff1a;光波在光…
阅读更多...
Java——封装、访问修饰符、包

Java——封装、访问修饰符、包

目录 一.封装的概念 二.访问限定符 三.封装扩展之包 1.包的概念 2.导入包中的类 3.自定义包 4.包的访问权限控制举例 5.常见的包 一.封装的概念 面向对象程序三大特性&#xff1a;封装、继承、多态。而类和对象阶段&#xff0c;主要研究的就是封装特性。何为封装呢&am…
阅读更多...
electron项目打包慢、打包报错

electron项目打包慢、打包报错

项目使用了electron框架&#xff0c;在第一次打包或者网络条件不好的环境下进行打包时熟速度慢的出奇&#xff0c;甚至经常出现打包失败的情况&#xff08;如上面图片的报错&#xff09;。 这是因为&#xff0c;在electron打包的过程中&#xff0c;需要去官方源https://github.…
阅读更多...
如何在 VM 虚拟机中安装 OpenEuler 操作系统保姆级教程（附链接）

如何在 VM 虚拟机中安装 OpenEuler 操作系统保姆级教程（附链接）

一、VMware Workstation 虚拟机 若没有安装虚拟机的可以参考下篇文章进行安装&#xff1a; 博客链接https://eclecticism.blog.csdn.net/article/details/135713915 二、OpenEuler 镜像 点击链接前往官网 官网 选择第一个即可 三、安装 OpenEuler 打开虚拟机安装 Ctrl …
阅读更多...
客户资料不翼而飞？企业数据保护攻略

客户资料不翼而飞？企业数据保护攻略

在数字化经济时代&#xff0c;企业的客户资料等同于商业生命线&#xff0c;一旦泄露&#xff0c;后果不堪设想。例如&#xff0c;2017年Equifax的数据泄露事件&#xff0c;造成超过1.4亿用户的个人信息外泄&#xff0c;不仅给用户带来风险&#xff0c;也让公司名誉受损&#xf…
阅读更多...
IP定位技术原理详细阐述

IP定位技术原理详细阐述

IP定位技术原理主要基于IP地址与地理位置之间的关联&#xff0c;通过一系列的技术手段&#xff0c;实现对网络设备的物理位置进行精确或大致的定位。以下是对IP定位技术原理的详细阐述。 首先&#xff0c;我们需要了解IP地址的基本概念。IP地址是互联网协议地址的简称&#xff…
阅读更多...
【Linux】Socket编程接口 | 实现简单的UDP网络程序

【Linux】Socket编程接口 | 实现简单的UDP网络程序

文章目录 一、预备知识理解源IP地址和目的IP地址理解源mac地址和目的mac地址认识端口号理解源端口号和目的端口号理解“端口号&#xff08;PORT&#xff09;”和“进程ID&#xff08;PID&#xff09;” 认识TCP和UDP协议TCP协议UDP协议 网络字节序为什么网络字节序采用的是大端…
阅读更多...
linux项目部署 解决Nginx浏览器刷新出现404，但是不刷新是能够正常请求成功

linux项目部署 解决Nginx浏览器刷新出现404，但是不刷新是能够正常请求成功

文章目录 目录 文章目录 安装流程 小结 概要安装流程技术细节小结 概要 提示&#xff1a;部署成功&#xff0c;访问登录页面登录也成功&#xff0c;强制刷新浏览器报404问题 进入到系统 刷新页面 解决流程 参考如图&#xff0c;再下面添加这条配置信息 location / {try_file…
阅读更多...
中职人工智能技术与应用课程及就业方向

中职人工智能技术与应用课程及就业方向

泰迪智能科技推出中职人工智能专业一站式建设解决方案&#xff0c;从解决人工智能专业的专业建设难点出发&#xff0c;让教学实训更简单为理念&#xff0c;助力学校完成人工智能相关专业的专业建设。 人工智能技术应用专业课程&#xff1a; 人工智能应用导论、程序设计基础…
阅读更多...
推荐文章
最新文章

Copyright @ 2022~2023