辗转相除法和同余原理

辗转相除法和同余原理

  • 辗转相除法
  • 同余原理

辗转相除法

  辗转相除法就是用来求出两个数的最大公约数的方法,那么这个方法怎么用呢?举个例子:有两个数,a=12,b=8,要求这两个数的最大公约数,首先让a%b得到4,然后让a=b,b=a%b,即现在a=8,b=4;继续用a%b得到0,然后让a=b,b=a%b,现在a=4,b=0。当b等于0的时候,a的值就是原来两个数的最大公约数
  代码如下

int gcd(int a, int b)
{return b == 0 ? a : gcd(b, a % b);
}

证明辗转相除法就是证明以下关系
g c d ( a , b ) = g c d ( b , a m o d b ) a gcd(a,b)=gcd(b,a\bmod b )a gcd(a,b)=gcd(b,amodb)a
假设 a m o d b = r a\bmod b=r amodb=r
则需要证明 g c d ( a , b ) = g c d ( b , r ) gcd(a,b)=gcd(b,r ) gcd(a,b)=gcd(b,r)
因为 a m o d b = r a\bmod b=r amodb=r成立,所以下面连等式必定成立
a = b ∗ q + r \begin{align} a=b*q+r \end{align} a=bq+r
r = a − b ∗ q \begin{align} r=a-b*q \end{align} r=abq
其中q为自然整数
假设u是a和b的公共因子,则有:
a = s ∗ u b = t ∗ u \begin{split} a=s*u\\ b=t*u \end{split} a=sub=tu
将a和b带入到式子(1)当中得到
r = s ∗ u − t ∗ u ∗ q = ( s − t ∗ q ) ∗ u \begin{split} r &= s*u-t*u*q\\ &= (s-t*q)*u \end{split} r=sutuq=(stq)u
  这说明,u如果是a和b的公因子,那么u也是r的公因子,假设v是b和r的公因子,同理可得v也是a的公因子。
  综上,a和b的每个公因子,也是b和r的公因子,反之亦然。所以a和b的全体公因子集合 = b和r的全体公因子集合
g c d ( a , b ) = g c d ( b , a m o d b ) a gcd(a,b)=gcd(b,a\bmod b )a gcd(a,b)=gcd(b,amodb)a
证毕
  在求出a和b最大公约数r之后,最小公倍数就是
a ∗ b / r a*b/r ab/r
  转换成代码就是

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

  在了解到以上前置知识后,来做一道题目
在这里插入图片描述
神奇的数字
  就用示例2为例,两个数分别是2和3,要求第4个神奇的数字,其实我们可以给出一个范围,使得第4个神奇的数字是在这个范围里面的,因为是第4个神奇的数字,那么这个数字最大不会超过2*4也就是8,范围就可以确定为【0,8】,因为即使没有3,求得可以被2整除的第4个数字最大也就是8了,由于又加了一个数字,那么在范围【0,8】内神奇数字的数目肯定是变多了的,那么第4个神奇数字必定是在【0,8】之间
  在确定了范围之后,我们可以定义一个函数f(x),该函数的返回值就是从0到x包含有多少个神奇数字。还是以2和3为例
在这里插入图片描述
  可以看到,在【0,8】范围内6既是2的倍数也是3的倍数,所以不管是对于2来说还是对于3来说,6都是神奇数字,即6被算作了两次神奇数字,那么我们只要进行去重,得到的神奇数字个数就是函数f(x)的返回值
  在确定了范围还有f(x)后,我们就可以使用二分法解决这道问题
代码如下:

long long gcd(int a, int b) //求最大公约数{return b == 0 ? a : gcd(b, a % b);}long long lcm(int a, int b) //求最小公倍数{return a * b / gcd(a, b);}//检查x右边有多少个数是可以被a或者b整除long long check(long long x, int a, int b) {return x / a + x / b - x / lcm(a, b);}int nthMagicalNumber(int n, int a, int b) {long long mod = 1e9 + 7;long long l = 0;long long r = (n % mod) * min(a, b);while (l < r) {long long m = (l + r) / 2;if (check(m, a, b) >= n)r = m;elsel = m + 1;}return l % (mod);}

同余原理

加法的同余原理
( a + b ) m o d c = a m o d c + b m o d c (a+b)\bmod c=a\bmod c+b\bmod c (a+b)modc=amodc+bmodc
乘法的同余原理
( a ∗ b ) m o d c = ( a m o d c ) ∗ ( b m o d c ) (a*b)\bmod c=(a\bmod c)*(b\bmod c) (ab)modc=(amodc)(bmodc)

  以加法的同余原理为例,当题目要求结果对c进行取模运算时,但是在得到结果的时候,也就是a+b运算过程中,加出来的结果已经超出了当前类型的最大值,就发生了溢出的现象,这个时候再对溢出后的结果进行取模已经没有意义了,所以就需要在运算之前就分别对两个数进行取模防止溢出

减法的同余原理
( a − b ) m o d c = ( ( a m o d c ) − ( b m o d c ) + c ) m o d c (a-b)\bmod c=((a \bmod c)-(b \bmod c)+c) \bmod c (ab)modc=((amodc)(bmodc)+c)modc
  在减法的同余原理中需要对
( a m o d c ) − ( b m o d c ) (amodc)−(bmodc) (amodc)(bmodc)加上c是为了防止两数相减得到一个负数的情况

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

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

相关文章

【Linux笔记】进程间通信之管道

一、匿名管道 我们在之前学习进程的时候就知道了一个概念&#xff0c;就是进程间是互相独立的&#xff0c;所以就算是两个进程是父子关系&#xff0c;其中一个进程退出了也不会影响另一个进程。 也因为进程间是互相独立的&#xff0c;所以两个进程间就不能直接的传递信息或者…

使用一根网线,让Ubuntu和正点原子I.MX6ULL开发板互相ping通

1.硬件准备 准备一根网线即可 2. 让windows和I.MX6ULLping通 2.1 找根网线将I.MX6ULL和电脑连起来 2.2 让I.MX6ULL通电运行起来&#xff0c;我这里使用的是正点原子版本的内核、 2.3 进入电脑的网络连接后&#xff0c;按照如下步骤操作 2.4 将ip地址、子网掩码、默认网关…

【数据结构】二叉树的三种遍历

目录 一、数据结构 二、二叉树 三、如何遍历二叉树 一、数据结构 数据结构是计算机科学中用于组织和存储数据的方式。它定义了数据元素之间的关系以及对数据元素的操作。常见的数据结构包括数组、链表、栈、队列、树、图等。 数组是一种线性数据结构&#xff0c;它使用连续…

嵌入式系统的基础知识:了解嵌入式系统的构成和工作原理

&#xff08;本文为简单介绍&#xff0c;个人观点仅供参考&#xff09; 嵌入式系统是建立在微处理器基础上的计算机系统,用于对专门的功能进行控制、运算和接口。它结合了硬件和软件,可以提供实时的响应,广泛应用于工业控制、通信、医疗、交通等领域。 嵌入式系统的核心是微处理…

亿级推送,得物是怎么架构的?

说在前面 在40岁老架构师 尼恩的读者交流群(50)中&#xff0c;很多小伙伴拿到一线互联网企业如阿里、网易、有赞、希音、百度、滴滴的面试资格。 最近&#xff0c;尼恩指导一个小伙伴简历&#xff0c;需要织入亮点项目、黄金项目。 前段时间&#xff0c;指导小伙写了一个《高…

###51单片机学习(2)-----如何通过C语言运用延时函数设计LED流水灯

前言&#xff1a;感谢您的关注哦&#xff0c;我会持续更新编程相关知识&#xff0c;愿您在这里有所收获。如果有任何问题&#xff0c;欢迎沟通交流&#xff01;期待与您在学习编程的道路上共同进步。 目录 一. 延时函数的生成 1.通过延时计算器得到延时函数 2.可赋值改变…

信息学奥赛一本通1314:【例3.6】过河卒(Noip2002)

1314&#xff1a;【例3.6】过河卒(Noip2002) 时间限制: 1000 ms 内存限制: 65536 KB 提交数: 40991 通过数: 17884 【题目描述】 棋盘上A点有一个过河卒&#xff0c;需要走到目标B点。卒行走的规则&#xff1a;可以向下、或者向右。同时在棋盘上的某一点有一个对方…

【Go语言】第一个Go程序

第一个 Go 程序 1 安装 Go Go语言官网&#xff1a;Download and install - The Go Programming Language&#xff0c;提供了安装包以及引导流程。 以 Windows 为例&#xff0c;进入windows安装包下载地址&#xff1a;All releases - The Go Programming Language&#xff0c…

TeamCity创建git项目Timed out 超时的一个解决办法

问题&#xff1a; 当自己&#xff1a; ping github.com从本地推送到远程仓库浏览器浏览www.github.com ——都没有问题 但是在teamcity创建工程的时候就超时&#xff1a; 或者多试几次&#xff0c;终于成功了&#xff0c;然后构建的时候半途超时报错。。。。。 一种解决办…

Mermaid绘制UML图教程

Mermaid 是一种轻量级的图形描述语言&#xff0c;用于绘制流程图、时序图、甘特图等各种图表。它采用简单的文本语法&#xff0c;使得用户能够快速绘制各种复杂图表&#xff0c;而无需深入学习图形绘制工具。 一、安装Mermaid Mermaid 可以在浏览器中直接使用&#xff0c;也可…

Linux rp_filter、arp_filter、arp_ignore、arp_announce参数说明

Linux rp_filter、arp_filter、arp_ignore、arp_announce参数说明。我查看了参考资料&#xff0c;又去查阅了官方文档&#xff0c;凭着我的理解整理了以下文档。各位大神的文档写的很好&#xff0c;但都不喜欢断句啊&#xff0c;读的我这叫一个累。 参考 1.网络编程之网络丢包…

2024LeetCode分类刷题

一、数组 88. 合并两个有序数组 public void merge(int[] nums1, int m, int[] nums2, int n) {int p1 0, p2 0;int[] sorted new int[m n];while (p1 < m || p2 < n) {int current;if (p1 m) {current nums2[p2];} else if (p2 n) {current nums1[p1];} else i…