GCD和LCM算法

news2024/11/25 7:00:12

目录

一 整除

定义

性质

二 GCD

1)定义

2)性质

 3)GCD编程

①暴力法

②欧几里得算法

③更相减损术

④Stein算法

三 LCM 

①暴力法

②最大公约数法

四 裴蜀定理

 例题:裴蜀定理

五 算法实践

[蓝桥杯 2019 省 B] 等差数列


一 整除

定义

a 能整除b,记为 a|b。其中,a 和 b为整数,且a不等于0, b是a 的倍数,a 是b的约数(因子)。

如:13|182,一5|35,-3|36;

6 的因子是士1、2、土3、土6。

性质

1)若a,b,c为整数,且alb、blc,则alc;

2)若a、b、m,n为整数,且cla、clb,则cl(ma+nb);

3)定理:带余除法。如果a和b为整数且b>0,则存在唯一的整数 q、r,使a=bq+r,0<r<b

二 GCD

1)定义

整数 a 和6的最大公约数是指能同时整除a 和b 的最大整数,记为 gcd(a,b)。

例如:gcd(15,81)=3,gcd(0,44)=44,gd(0,0)=0,ged(-6,-15)=3,ged(-17,289)=17。

注意:由于-a 的因子和a 的因子相同,因此 ged(a,b)=ged(lal,lbl)。编码时只需要关注正整数的最大公约数。

2)性质

1) gcd(a,b)=gcd(a,a+b)=gcd(a, ka+b)

2) gcd(ka,kb)=k*gcd(a, b)。

3) 定义多个整数的最大公约数: gcd(a,b,c)=gcd[gcd(a,b),c]。

4)若 gcd(a,b)=d,则 gcd(a/d,b/d)=1,即a/d 与b/d 互素。

5) gcd(a+cb,b)=gcd(a,b).

 3)GCD编程

①暴力法

#include<stdio.h>
#include<stdlib.h>
int main()
{
     int a,b,i;
     scanf("%d",&a);
     scanf("%d",&b);
     i=fmin(a,b);
     while(i)
     {
          if(a%i==0&&b%i==0)
           break;
      i--; 
     } 
 printf("%d",i);
 return 0;
}

②欧几里得算法

辗转相除法原理:用较小数除较大数,再用出现的余数(第一余数)去除除数,再用出现的余数(第二余数)去除第一余数,如此反复,直到最后余数是0为止。如果是求两个数的最大公约数,那么最后的除数就是这两个数的最大公约数。这个和更相减损术有着异曲同工之处。

辗转相除法求 GCD

gcd(a,b)=gcd(b,a mod b)。代码如下

int gcd(int a, int b)//一般要求a>=0,b>0.若a=b=0,代码也正确,则返回0
{
    return b? gcd(b, a% b):a;
}

这是最常用的方法,它极为高效

拉梅定理给出了复杂度分析

拉梅定理:用欧几里得算法计算两个正整数的最大公约数,需要的除法次数不会超过两个整数中较小的那个十进制数的位数的 5 倍。

推论:用欧几里得算法求 gcd(a,b),a>b,需要 O((loga)3)次位运算。

欧几里得算法的缺点是需要做取模运算,而高精度的除法取模比较耗时,此时可以使用“更相减损术

和 Stein 算法,它们只用到了减法和移位操作。

③更相减损术

定义:(如果需要对分数进行约分,那么)可以折半的话,就折半(也就是用2来约分)。如果不可以折半的话,那么就比较分母和分子的大小,用大数减去小数,互相减来减去,一直到减数与差相等为止,用这个相等的数字来约分。

算法的计算基于这一性质: gcd(a,b)=gcd(b,a-b)=gcd(a,a-b)。计算步骤:用较大的数减较小的

数,把所得的差与较小的数比较,然后继续做减法操作,直到减数与差相等为止。

int gcd(int a, int b)
{	
	while(a != b)
	{
	if(a> b)   a=a-b;
	else
		b=b-a;
	}
return a;
}

更相减损术虽然避免了欧几里得的取模计算,但是计算次数比欧几里得算法多很多,极端情况下需

要计算 O(max(a,b))次,如a=100,b=1时,需计算 100次

④Stein算法

Stein 算法是更相减损术的改进。求 gcd(a,b)时,可以分为几种情况进行优化。

(1)a 和b 都是偶数。gcd(a,b)=2gcd(a/2,b/2),计算减半。

(2)a 奇b偶。根据原理:若k和y互为质数有 gcd(kx,y)=gcd(x,b)。当k=2,b 为奇数时,有

gcd(a,b)=gcd(a/2,b),即偶数减半。表示 bb 存在2这个因子而 aa 不存在,则将 bb 除以2,,不考虑因子2;

(3)a 偶b奇。gcd(a,b)=gcd(a,b/2),

(4)a 和b 都是奇数。gcd(a,b)=gcd((a+b)/2,(a-b)/2)。

算法的结束条件仍然是 gcd(a,a)=a。

除 2 操作用移位就可以了,所以 Stein 算法只用到加减法和移位;

三 LCM 

 a 和b 的最小公倍数表示为 lcm(a,b),从算术基本定理推理得到。

算术基本定理:任何大于1的正整 数n 都可以唯一分解为有限个的乘积;

①暴力法

#include<stdio.h>
#include<stdlib.h>
int main()
{
     int a,b,i;
     scanf("%d",&a);
     scanf("%d",&b);
     i=fmax(a,b);
     while(i)
     {
      if(i%a==0&&i%b==0)
       break;
      i++; 
 } 
 printf("%d",i);
 return 0;
}

可以推出 gcd(a,b)lcm(a,b)=ab,即 lcm(a,b)=ab/gcd(a,b)=a/gcd(a,b)b。

注意,要先作除法再作乘法,如果先作乘法可能会溢出。

②最大公约数法

int lcm(int a, int b)
{
    return a / gcd(a,b) * b;
}

四 裴蜀定理

裴蜀定理是关于 GCD的一个定理。

裴蜀定理(Bezout’s Lemma):如果a 与b 均为整数,则有整数x和y使ax+by=gcd(a,b)。

这个等式称为Bezout 等式

推论:整数a 与b互素当且仅当存在整数 x和y,使ax+by=1.

裴蜀定理很容易证明。

可以这样理解裴蜀定理:对任意x和y,  d =ax+by,d 一定是gcd(a,b)的整数倍;最小的 d 是 gcd(a,b).

 例题:裴蜀定理

代码 

#include<stdio.h>
#define ll long long
int gcd(int a, int b)
{	
	return b? gcd(b,a%b):a;
}

int main()
{
	int n,ans=0,tmp,i;
	scanf("%d",&n);
	for(i=1;i<=n;i++)
	{
		scanf("%d",&tmp);
		if(tmp<0)	tmp=-tmp;
		ans=gcd(ans,tmp); 
	}
	printf("%d",ans);
	return 0;
}

五 算法实践

[蓝桥杯 2019 省 B] 等差数列

数学老师给小明出了一道等差数列求和的题目。但是粗心的小明忘记了一 部分的数列,只记得其中 NN 个整数。

现在给出这 NN 个整数,小明想知道包含这 NN 个整数的最短的等差数列有几项?

输入描述

输入的第一行包含一个整数 NN。

第二行包含 NN 个整数 A_1,A_2,··· ,A_NA1​,A2​,⋅⋅⋅,AN​。(注意 A_1A1​ ∼ A_NAN​ 并不一定是按等差数列中的顺序给出)

其中,

输出描述

输出一个整数表示答案。

输入输出样例

示例

5
2 6 4 10 20

10

样例说明: 包含 2、6、4、10、20 的最短的等差数列是 2、4、6、8、10、12、14、16、 18、20。

代码:

#include<stdio.h>
#define ll long long
ll cmp(const void *a,const void *b)
{
	return *(int *)a>*(int *)b;
}
ll gcd(ll a, ll b)
{	
	return b? gcd(b,a%b):a;
}


int main()
{
	ll n,a[100010],i,x;
	scanf("%lld",&n);
	for(i=1;i<=n;i++)
		scanf("%lld",&a[i]);
	qsort(&a[1],n,sizeof(a[1]),cmp);
	x=gcd(a[1],a[2]);
	printf("%lld",(a[n]-a[1])/x+1);
	return 0;
}

上述代码运行发现可以通过给出的样例,但其实是不对的。

代码:

#include<stdio.h>
#define ll long long
ll cmp(const void *a,const void *b)
{
	return *(int *)a>*(int *)b;
}
ll gcd(ll a, ll b)
{	
	return b? gcd(b,a%b):a;
}

/*
等差数列的项数公式:,an就是等差数列中的最大值,a1就是等差数列中的最小值

题目要求数列最短,因为最大值和最小值已经固定,所以要求的就是公差d

误区:数组两两之间的最小差值并非我们的公差!!!

例如:3,5,8;最小的差值是2,但是这三个数在内的等差数列其公差应该是1,

即3,4,5,6,7,8,所以这个地方我们需要对所有的差值求最大公约数
*/
int main()
{
	ll n,a[100010],i,x=0;
	scanf("%lld",&n);
	for(i=1;i<=n;i++)
		scanf("%lld",&a[i]);
	qsort(&a[1],n,sizeof(a[1]),cmp);
	for(i=2;i<=n;i++)
		x=gcd(x,a[i]-a[i-1]);
	if(x==0)
		printf("%lld",n);
	else
		printf("%lld",(a[n]-a[1])/x+1);
	return 0;
}

 

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

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

【AIOT】语音学习

声音以波的形式传播&#xff0c;即声波&#xff08;Sound Wave&#xff09;。当我们以波的视角来理解声音时&#xff0c;却又大繁若简起来&#xff1a;仅凭频率&#xff08;Frequency&#xff09;、幅度&#xff08;Magnitude&#xff09;、相位&#xff08;Phase&#xff09;便…

19 CPP拷贝构造函数

函数以值的方式返回对象时&#xff0c;可能会调用拷贝构造函数&#xff08;VS会调用&#xff0c;Linux不会&#xff0c;g编译器做了优化&#xff09; 代码说明 #include <iostream>using namespace std;//函数以值的方式返回对象时&#xff0c;可能会调用拷贝构造函数&a…

C++对C的语法增强

1、语法检查增强 声明int a; 定义int a 10; c语言的弱语法类型 2、严格的类型转换 枚举类型不能赋值整形 3、struct类型加强 c中定义结构体变量需要加上struct关键字&#xff0c;c不需要 c中结构体只能定义成员变量&#xff0c;不能定义函数。c可以定义成员函数 4、bool类…

聊聊优惠卷秒杀的下单功能

前言 案例来自黑马程序员视频&#xff1a;https://www.bilibili.com/video/BV1cr4y1671t/?spm_id_from333.999.0.0 案例分析 下单时需要判断两点&#xff1a; 秒杀是否开始或者结束&#xff0c;如果尚未开始或者已经结束则无法下单&#xff1b;库存是否充足&#xff0c;不足…

【图像隐写】DWT+SVD数字水印【含GUI Matlab源码 606期】

⛄一、DWTSVD数字水印简介 理论知识参考文献&#xff1a;基于DWT和SVD的彩色图像数字水印算法研究 一种基于DWT-SVD的图像数字水印算法\ ⛄二、部分源代码 function varargout main(varargin) % MAIN M-file for main.fig % MAIN, by itself, creates a new MAIN or raises…

想要提升英语阅读能力的小伙伴,走过路过,千万不要错过!

一、英语四六级阅读提升网站&#xff1a; 据统计&#xff1a;四六级阅读题目的来源主要来自纽约时报、经济学人、大西洋月刊、时代、新闻周刊、卫报等报刊&#xff0c;其分布比重如下图&#xff1a; 原文链接&#xff1a;英语四六级阅读外刊来源&#xff0c;知道事半功倍&#…

IIC信号为什么要加上拉电阻

IIC是一个两线串行通信总线&#xff0c;包含一个SCL信号和SDA信号&#xff0c;SCL是时钟信号&#xff0c;从主设备发出&#xff0c;SDA是数据信号&#xff0c;是一个双向的&#xff0c;设备发送数据和接收数据都是通过SDA信号。 在设计IIC信号电路的时候我们会在SCL和SDA上加一…

计算机毕业设计springboot+vue基本微信小程序的外卖点餐平台

项目介绍 餐饮行业是一个传统的行业。根据当前发展现状,网络信息时代的全面普及,餐饮行业也在发生着变化,单就点餐这一方面,利用手机点单正在逐步进入人们的生活。传统的点餐方式,不仅会耗费大量的人力、时间,有时候还会出错。小程序系统伴随智能手机为我们提供了新的方向。手机…

基于vue2整合tinymce富文本编辑器

文章目录1.导入依赖2.全部代码3.实际效果1.导入依赖 npm install tinymce -Snpm install tinymce/tinymce-vue -S2.全部代码 <template><editor:key"tinymceFlag"ref"editor"v-model"content":init"tinymceInit"></e…

国内使用ChatGPT的方式汇总

最近ChatGPT大火呀&#xff0c;小伙伴们是不是在网上看到各种和ChatGPT有趣聊天的截图&#xff0c;奈何自己实力不够&#xff0c;被网络拒之门外&#xff0c;只能眼馋别人的东西。看别人玩&#xff0c;肯定不如自己玩一把舒服的啊。今天小卷就给大家汇总了国内能使用ChatGPT的方…

【C/C++ 数据结构】-就这一篇博客让你玩爆二叉树的各种遍历问题!!!

作者&#xff1a;学Java的冬瓜 冬瓜的主页&#xff1a;☀冬瓜的主页&#x1f319; 专栏&#xff1a;【C/C 数据结构与算法】 分享&#xff1a;被苦难淬炼过的人&#xff0c;内心真诚。——都靓评宋濂 主要内容&#xff1a;二叉树的递归前序遍历、中序遍历、后序遍历、层序遍历。…

你还会想起这道题吗(another version)

链接&#xff1a;登录—专业IT笔试面试备考平台_牛客网 来源&#xff1a;牛客网 题目描述 众所周知&#xff0c;原神是由米哈游自主研发的一款全新开放世界冒险游戏。游戏发生在一个被称作「提瓦特」的幻想世界&#xff0c;在这里&#xff0c;被神选中的人将被授予「神之眼…

“价格+产品+服务”:京东手机开启“复合竞赛”

市场唯一不变的就是变化&#xff0c;那些最终生存下来的企业往往是最能适应变化的。 11月份&#xff0c;京东再次强调“低价”策略的重要性&#xff0c;背后无疑是其对消费环境变化做出的及时反应。“消费者对价格更为敏感&#xff0c;对高性价比商品的需求依旧是主流”&#…

YOLO-V5 算法和代码解析系列 —— 学习路线规划综述

目录标题为什么学习 YOLO-V5 &#xff1f;博客文章列表面向对象开源项目学习方法预备知识项目目录结构为什么学习 YOLO-V5 &#xff1f; 算法性能&#xff1a;与YOLO系列&#xff08;V1&#xff0c;V2&#xff0c;V3&#xff0c;V4&#xff09;相比&#xff0c;YOLO-V5效果最好…

【花雕动手做】有趣好玩的音乐可视化系列项目(33)---核酸托盘灯

偶然心血来潮&#xff0c;想要做一个音乐可视化的系列专题。这个专题的难度有点高&#xff0c;涉及面也比较广泛&#xff0c;相关的FFT和FHT等算法也相当复杂&#xff0c;不过还是打算从最简单的开始&#xff0c;实际动手做做试验&#xff0c;耐心尝试一下各种方案&#xff0c;…

RosonQt140——Qt Charts模块介绍和Qt绘制图表

Qt图表概述 Qt Charts能够创建时尚、互动、以数据为中心的用户界面。Qt Charts使用图形视图框架&#xff0c;便于集成。图表组件可以作为QWidget或QGraphicsWidget对象或QML类型使用。 QChart类管理不同类型的系列和其他图表相关对象的图形表示&#xff0c;如图例和坐标轴。QC…

33-Vue之ECharts高级-设置主题

ECharts高级-设置主题前言内置主题自定义主题前言 本篇来学习下ECharts中如何设置图表主题 内置主题 ECharts 中默认内置了两套主题: light dark var chart echarts.init(dom, light) var chart echarts.init(dom, dark)<!DOCTYPE html> <html lang"en&quo…

前缀树介绍,定义,图文详解分析——Java/Kotlin双版本代码

前缀树 前缀树&#xff0c;又称作字典树&#xff0c;用一个树状的数据结构储存字典中的所有单词。 列&#xff0c;一个包含can、cat、come、do、i、in、inn的前缀树如下图所示&#xff1a; 前缀树是一个多叉树&#xff0c;一个节点可能存在多个节点。除根节点外&#xff0c;每…

Spring5框架总结学习(从入门到进阶)

文章目录Spring51、如何创建一个Spring项目&#xff08;idea版&#xff09;2、 IOC容器1、XML解析工厂模式反射2、bean管理1、总述2、基于XML创建对象3、基于XML注入属性4、基于XML注入属性&#xff0c;属性值为空或特殊符号5、基于XML注入属性&#xff0c;外部bean6、基于XML注…

【第十二章 MVCC(多版本并发控制),隐藏字段,undolog(版本链),readview,原理分析(RC,RR)】

第十二章 MVCC&#xff08;多版本并发控制&#xff09;&#xff0c;隐藏字段&#xff0c;undolog&#xff08;版本链&#xff09;&#xff0c;readview&#xff0c;原理分析&#xff08;RC&#xff0c;RR&#xff09; 1.基本概念: &#xff08;1&#xff09; 当前读&#xff1…