c++用数组int gcd辗转相除法_python(四)：辗转相除

最新推荐文章于 2023-11-11 00:34:34 发布
最新推荐文章于 2023-11-11 00:34:34 发布
阅读量269 收藏
点赞数
文章标签： c++用数组int gcd辗转相除法
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_39672572/article/details/113626565
版权

摘要：

学习leetcode_365: Water and Jug Problem的解法：辗转相除相关内容（最大公约数、裴蜀定理、欧几里得算法和扩展欧几里得算法）。

正文：

1、问题描述

You are given two jugs with capacities x and y litres. There is an infinite amount of water supply available. You need to determine whether it is possible to measure exactly z litres using these two jugs.

If z liters of water is measurable, you must have z liters of water contained within one or both buckets by the end.

Operations allowed:

  • Fill any of the jugs completely with water.
  • Empty any of the jugs.
  • Pour water from one jug into another till the other jug is completely full or the first jug itself is empty.

Example 1: 

Input: x = 3, y = 5, z = 4
Output: True

（**分析**）：

首先，将三个操作进行“量化”，一次倒满已知体积的壶的操作，相当于有一个(+1)*Volume，而一次将满壶水倒掉的操作，则相当于(-1)*Volume。这意味着结果可以表示为两个已知体积壶的整数线性组合：

equation?tex=z+%3D+mx%2Bny

所以，任意给定（x, y, z）且z<(x + y)后，那么就看上述方程是否有整数解（m, n），如果整数解存在的话结果就是True，否则就是False。

2、整数解存在的判据：裴蜀定理

前述方程是否有整数解的解法有一个已经存在的定理（裴蜀定理），其结论是该方程有整数解（m, n）的充要条件是z为x与y的最大公约数gcd(x, y)的倍数。也就是说，如果这个方程有整数解，那么z一定是

equation?tex=W%2Agcd%28x%2C+y%29%28W+%5Csubseteq+Z%29 ，反过来看，也只有当z%gcd(x,y)=0时，整数解才存在。于是，
这个问题转化为求x和y的最大公约数，有辗转相除法（欧几里得算法）：

（2.1）核心思想，是如下等式，对于正整数p和q（若p>q），p除以q的余数为r，则有：

equation?tex=gcd%28p%2Cq%29%3Dgcd%28q%2Cr%29

由于显然有p>q>r，这个等式将两个较大数的公约数求解转化为两个较小数的公约数；

（2.2）这个最大公约数等式可以往下继续写：

equation?tex=gcd%28p%2Cq%29%3Dgcd%28q%2Cr_1%29%3Dgcd%28r_1%2Cr_2%29%3D...gcd%28c%2C0%29

其中从第一个等式开始，后式括号里的第一项为前式括号里的较小数、第二项为前式括号里较大数和较小数相除的余数。最后总能到一个余数为0的情况，此时的第一项即为p和q的最大公约数；

（2.3）这个算法可以用递归的思想很好地表达出来，python代码如下：

d802c82224b09b23d7baad6c9d6dbb09.png

test of gcd：

4b76eba7624e153d3a728dccbdfbcb76.png
这两组数的最大公约数都是6

那么，leetcode_365就迎刃而解了，代码如下：

dc0dc850c7d035445068a7d220a237e4.png

test of canMeasureWater：

c365e1b8f4fd718cdd00d48369462da0.png
题目的example

2e8e600a05cc33252b3c9cbe2e2ec6cc.png
随机输入3个正整数进行测试

3、整数解的求解：扩展欧几里得算法

判断了能否测量，进一步就会思考如何进行测量呢？这个实际上就是在求前述方程的整数解。这里的解法类似于辗转相除法，也即扩展欧几里得算法。需要注意的是，由于px+qy=c的整数解从几何意义上来说是检查一条直线过整数格点的情况，所以当然可以不止一个解，而扩展欧几里得算法得到的只是其中一个。

首先从一个例子来感受一下该算法：

equation?tex=20x_0%2B12y_0%3D4 。

（a）方程系数(20, 12)变为(12, 8)：

equation?tex=12x_0%2B8x_0%2B12y_0%3D4+%5CRightarrow+12%28x_0%2By_0%29%2B8x_0+%3D4

假设

equation?tex=x_1%3Dx_0%2By_0%EF%BC%8Cy_1%3Dx_0%28x_0%3Dy_1%2C+y_0%3Dx_1-y_1%29 ，则变为：

equation?tex=12x_1%2B8y_1%3D4

（b）继续将(12, 8)变为(8, 4)：

equation?tex=8x_1%2B4x_1%2B8y_1%3D4+%5CRightarrow+8%28x_1%2By_1%29%2B4x_1+%3D4

假设

equation?tex=x_2%3Dx_1%2By_1%EF%BC%8Cy_2%3Dx_1%28x_1%3Dy_2%2C+y_1%3Dx_2-y_2%29 ，则变为：

equation?tex=8x_2%2B4y_2%3D4

（c）继续将(8, 4)变为(2, 0)：

equation?tex=8x_2%2B4y_2%3D4+%5CRightarrow+4%282x_2%2By_2%29%3D4

假设

equation?tex=x_3%3D2x_2%2By_2%EF%BC%8Cy_3%3Dx_2%3D0%28x_2%3Dy_3%3D0%2C+y_2%3Dx_3-2y_3%3Dx_3%29 ，则变为：

equation?tex=4x_3%3D4

（d）最后一个解(x3, y3)=(1, 0)，那么从括号里的递推公式，依次有：(x2, y2)=(0, 1)，(x1, y1)=(1, -1)，(x0, y0)=(-1, 2)，最后一项即为原方程的一个解。

形如

equation?tex=px%2Bqy%3Dc 的方程总是可以用（a）-（d）来求解。整个过程和辗转相除类似（方程系数就是在做辗转相除），而且如果方程有解（c%gcd(p,q)=0），那最后的解肯定是(c/gcd(p,q), 0)，因为系数辗转相除总是会到(gcd(p,q), 0)。其中变量转换的递推公式比较关键，但也很好理解，推广一下就是（//和python的整数除法定义一致）：

equation?tex=x_n%3Dy_%7Bn%2B1%7D%2C+y_n%3Dx_%7Bn%2B1%7D-%28p_n%2F%2Fq_n%29y_%7Bn%2B1%7D

python实现代码如下：

16a4914f8fe1cb3c5c57d45e96d5f858.png
这里先求解px+qy=gcd(p,q)，再转换为px+qy=c的解；这样操作只是为了逻辑更清晰，实际上可以整合在一起（递归结束时返回值从(1, 0)变为(c/gcd(p,q), 0)）

2e7aa7a452cde1a1668f40a897bb6802.png
px+qy=gcd(p,q)的解，求解过程不依赖于方程左边的c；以(p,q)为递归过程，再反向递归求解(x,y)

test of integer_solution：

bbec8c680027d9cfaea5287ca53bb10c.png
实际上还有x=-2, y=2的解，多解的原因应该是前述(c)步骤中，最后一次变量转换实际上可以不遵从递推公式

从这个结果可知，如果想用3升和5升的水壶生成4升容积，解（8，-4）意味着总共会倒入8次3升水壶、倒掉4次5升水壶，而解（-2, 2）意味着总共会倒掉2次3升水壶、倒满2次5升水壶。当然，解只是总体概括，至于具体步骤还需要进一步地分解了，不过问题就已经变得很简单了。

4、本篇内容代码的github链接

gcd_euclidean_algorithm.py​github.com
weixin_39672572
关注
c++用数组int gcd辗转相除法_算法竞赛专题解析（20）：数论GCD和LCM
weixin_39909366的博客
01-17 560
本系列文章将于2021年整理出版。前驱教材：《算法竞赛入门到进阶》 清华大学出版社有建议请加QQ 群：567554289csdn同步：点击本文末尾“阅读原文”。目录1. GCD定义2. GCD性质3. GCD编码3.1 欧几里得算法3.2 更相减损术3.3 Stein算...
c++用数组int gcd辗转相除法_从辗转相除法到求逆元，数论算法初体验
weixin_32263265的博客
02-03 641
今天是算法和数据结构专题的第22篇文章，我们一起来聊聊辗转相除法。辗转相除法又名欧几里得算法，是求最大公约数的一种算法，英文缩写是gcd。所以如果你在大牛的代码或者是书上看到gcd，要注意，这不是某某党，而是指的辗转相除法。在介绍这个算法之前，我们先来看下最大公约数问题。暴力解法这个问题应该很明确了，我们之前数学课上都有讲过。给我们纸笔让我们求都没有问题，分解因数找下共同的部分，很快就算出来了。但...
有关辗转相除法原理和C++函数的实现
darkmoonbird的博客
11-11 324
辗转相除法的数形结合与代码实现
c#数组求余，辗转相除递归
weixin_44202936的博客
09-07 226
public int Gcd(int a,int b) { if(b%a==0) { return a; } else { return Gcd(b%a); } } main() { int gcd = 0; int[] num = new int[]{15,39,45,102}; for(int I =0;i<num.Length;i++) { gcd = Gcd(num[i],num...
c++求 两数 GCD
lijialinapple的专栏
05-12 1716
#include using namespace std; int GCD(int num1, int num2) { int mod=0; if ((mod = (num1%num2)) !=0) { return GCD(num2, mod); } else { return num2; } } int main() { int num1, num2; cou
python辗转相除法求最大公约数最小公倍数_辗转相除法求最大公约数/最小公倍数...
weixin_36205072的博客
01-14 3425
http://blog.csdn.net/jtujtujtu/article/details/44071712009辗转相除法求最大公约数：辗转相除法, 又名欧几里德算法(Euclidean algorithm)乃求两个正整数之最大公因子的算法。它是已知最古老的算法, 其可追溯至前300年。它首次出现于欧几里德的《几何原本》(第VII卷，命题i和ii)中，而在中国则可以追溯至东汉出现的《九章算术》...
辗转相除法
qq_29301113的博客
01-12 2428
辗转相除法， 又名欧几里德算法（Euclidean algorithm），是求最大公约数的一种方法。它的具体做法是：用较大数除以较小数，再用出现的余数（第一余数）去除除数，再用出现的余数（第二余数）去除第一余数，如此反复，直到最后余数是0为止。如果是求两个数的最大公约数，那么最后的除数就是这两个数的最大公约数。 另一种求两数的最大公约数的方法是更相减损法。 中文名 辗转相除法 外文名 Euclid...
辗转相除法（欧几里得算法）
upc_miao的博客
03-20 4868
辗转相除法，又叫欧几里德（Euclidean）算法， 是求最大公约数的算法。辗转相除法首次出现于欧几里得的《几何原本》（第VII卷，命题i和ii）中，而在中国则可以追溯至东汉出现的《九章算术》。 计算公式gcd(a,b) = gcd(b,a mod b)。 算法简介： 欧几里德算法是用来求两个正整数最大公约数的算法。是由古希腊数学家欧几里德在其著作《The Elements》中最早描述了这种算法,...
蓝桥杯 2016省赛 最大比例 从直觉到辗转相除法的变化
Lejeune的博客
03-16 350
蓝桥杯 2016省赛 最大比例 从直觉到辗转相除法的变化 题目链接 题面请看链接，分析一下题目，题目就是要求从一等比数列的某几项中得知，这个等比数列的最大公比。 看到题目显然是和gcd有关的，这里介绍两种方法，一种是从直觉而来的，一种是参照网上题解，自己推导而出的。 1.直觉题解         等比数列必然有一个首项...
gcd.rar_最大公约数
09-22
1. 辗转相除法（欧几里得算法）： 欧几里得算法是计算两个正整数最大公约数的最古老且最常用的方法。其基本思想是：对于任意两个正整数a和b（a>b），它们的最大公约数等于a除以b的余数c和b之间的最大公约数。反复...
GCD
a673786103的博客
03-18 205
#include&lt;iostream&gt; using namespace std; int count = 0; //ÆäÖÐ a&gt;b int gcd( int a, int b ){ count ++; if( b == 0 ) return a; return gcd(b,a%b); } int main() { int a = 27; int b = 81...
int gcd()求最大公约数
m0_73237155的博客
11-12 622
代码如下: #pragma warning(disable:4996) #include <stdio.h> #include <math.h> int gcd(int m, int n) { if (n == 0) return m; else gcd(n, m % n); } int main() { int m, n, r; do { printf("Input m,n"); scanf("%d%d", &a
辗转相除法（gcd、exgcd）
默的博客
11-24 350
用于求两个数的最大公约数。复杂度在O(log(max(a,b)))O(log(max(a,b)))O(log(max(a,b)))以内。 #include <cstdio> #include <iostream> using namespace std; int gcd(int a, int b) { if (b == 0) return a; return gcd(...
GCD算法
看了比懂王还懂
02-19 9655
网上代码多的是，我只是记下来而已 基础 int gcd(int a,int b) { int r; while(b>0) { r=a%b; a=b; b=r; } return a; } 位运算 int gcd(int a,int b) { while(b^=a^=b^
gcd公式
hpuw1234的博客
08-05 1334
求最大公约数： 1.迭代实现 int gcd(int a,intb) { while(b!=0) { int r=b; b=a%b; a=r; } return a; } 2.递归实现 int gcd(int a,int b) { if(b==0) return a; ret
GCD（求最大公约数）
qq_43731172的博客
11-23 456
#GCD 原理：类似辗转求余 两种方法： 1.循环 int main() { int a,b; scanf("%d%d",&amp;a,&amp;b); while(b) { int mid=a; a=b; b=mid%b; } printf("%d\n",a); return 0; } 2....
数组辗转相除法c++
最新发布
08-20
下面是使用C++实现数组中所有元素辗转相除法的一个示例代码： ```cpp #include using namespace std; // 函数声明 int gcd(int a, int b); int findGCDofArray(int arr[], int n); int main() { int arr[] = {...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值