求解最大公约数和最小公倍数问题
2015-04-10 16:50
239 查看
思路:
最大公约数问题也是一个非常典型的递归算法的应用。每次递归使得原来求两个大数之间的公约数转变成求两个稍微小点的数之间的公约数,要求转换的过程要保证不会改变公约数的值。这就要看其中转换的原理了。
原理从《几何原本》中得出--辗转相除。假设f(x, y) 表示x,y的最大公约数是g,而k = x/y,b= x%y,则g必能整出b。因为x = ky + b,b = x - ky,b/g = (x-ky)/g一定为整数,所以必有g整除b。
如下所示:
f(42, 30) = f(30, 12) = f(12, 6)= f(6, 0) = 6
代码如下:
[html] view
plaincopy
int gcd(int x , int y)
{
return (y == 0 )?x :gcd(y , x % y) ;
}
非递归代码:
书中还引伸出辗转相减法,原理跟上面所述差不多。
代码如下:
[html] view
plaincopy
int gcd(int x , int y)
{
if(x < y)
return gcd(y , x) ;
else if(y == 0)
return x ;
else
return gcd(x - y , y) ;
}
非递归代码:
最后书中提到一种相除和相减相结合的方法,保证了不做过多冗余的步骤。这就是把2当作每次转换的数量级。
若x,y均为偶数,f(x,y) = 2*f(x/2,y/2) = 2*f(x>>1,y>>1)
若x为偶数,y为奇数,f(x,y) = f(x/2,y) = f(x>>1,y)
若x为奇数,y为偶数,f(x,y) = f(x,y/2) = f(x,y>>1)
若x,y均为奇数,f(x,y) =f(y,x-y)
此解法结合了上述两种方法:第一种方法递归次数相当少但每次取模挺耗时(到底耗时情况如何,我没有具体了解过,但肯定比加减法耗时些),而第二种递归可能会相当多,但每次运算都是减法运算, 很快。综合两者,用2为基数,可以保证递归次数不那么多,同时运算也快。
代码如下:
[cpp] view
plaincopy
int gcd(int x , int y)
{
if(x < y)
return gcd(y , x) ;
if(y == 0)
return x ;
if(isEven(x)) //x为奇
{
if(isEven(y)) //y为奇
return gcd(x - y, y) ;
else //y为偶
return gcd(x , y >>1) ;
}
else //x为偶
{
if(isEven(y)) //y为奇数
return gcd(x >> 1, y) ;
else //y为偶
return 2 * gcd(x >> 1, y >> 1) ;
}
}
求解最小公倍数:
从题面上看,好像我们需要求解的是两个题目,但其实就是一个题目。那就是求最大公约数?为什么呢?我们可以假想这两个数m和n,假设m和n的最大公约数是a。那么我们可以这样写:
m = b *a;
n = c * a;
所以m和n的最小公倍数就应该是a*b*c啊,那不就是m * n / a,其中m和n是已知的,而a就是那个需要求解的最大公约数。所以就有了下面的代码,
[cpp] view
plaincopy
int GetMinCommonMultiple(int m, int n)
{
assert(m && n);
return m * n / gcd(m, n);
}
最大公约数问题也是一个非常典型的递归算法的应用。每次递归使得原来求两个大数之间的公约数转变成求两个稍微小点的数之间的公约数,要求转换的过程要保证不会改变公约数的值。这就要看其中转换的原理了。
原理从《几何原本》中得出--辗转相除。假设f(x, y) 表示x,y的最大公约数是g,而k = x/y,b= x%y,则g必能整出b。因为x = ky + b,b = x - ky,b/g = (x-ky)/g一定为整数,所以必有g整除b。
如下所示:
f(42, 30) = f(30, 12) = f(12, 6)= f(6, 0) = 6
代码如下:
[html] view
plaincopy
int gcd(int x , int y)
{
return (y == 0 )?x :gcd(y , x % y) ;
}
非递归代码:
int gcd(int x,int y) { int temp; while(y!=0) { temp=x; x=y; y=temp%y; } return x; }
书中还引伸出辗转相减法,原理跟上面所述差不多。
代码如下:
[html] view
plaincopy
int gcd(int x , int y)
{
if(x < y)
return gcd(y , x) ;
else if(y == 0)
return x ;
else
return gcd(x - y , y) ;
}
非递归代码:
int gcd(int n, int m) { if(n==0) return m; while (m!= 0) { if (n > m) n = n - m; else m = m - n; } return n; }或者:
int GetMaxCommonDivide_new(int n, int m) { if(n<m) { n^=m; m^=n; n^=m; //n=n+m; //m=n-m; //n=n-m; } while (m!= 0) { if (n > m) n = n - m; else m = m - n; } return n; }
最后书中提到一种相除和相减相结合的方法,保证了不做过多冗余的步骤。这就是把2当作每次转换的数量级。
若x,y均为偶数,f(x,y) = 2*f(x/2,y/2) = 2*f(x>>1,y>>1)
若x为偶数,y为奇数,f(x,y) = f(x/2,y) = f(x>>1,y)
若x为奇数,y为偶数,f(x,y) = f(x,y/2) = f(x,y>>1)
若x,y均为奇数,f(x,y) =f(y,x-y)
此解法结合了上述两种方法:第一种方法递归次数相当少但每次取模挺耗时(到底耗时情况如何,我没有具体了解过,但肯定比加减法耗时些),而第二种递归可能会相当多,但每次运算都是减法运算, 很快。综合两者,用2为基数,可以保证递归次数不那么多,同时运算也快。
代码如下:
[cpp] view
plaincopy
int gcd(int x , int y)
{
if(x < y)
return gcd(y , x) ;
if(y == 0)
return x ;
if(isEven(x)) //x为奇
{
if(isEven(y)) //y为奇
return gcd(x - y, y) ;
else //y为偶
return gcd(x , y >>1) ;
}
else //x为偶
{
if(isEven(y)) //y为奇数
return gcd(x >> 1, y) ;
else //y为偶
return 2 * gcd(x >> 1, y >> 1) ;
}
}
求解最小公倍数:
从题面上看,好像我们需要求解的是两个题目,但其实就是一个题目。那就是求最大公约数?为什么呢?我们可以假想这两个数m和n,假设m和n的最大公约数是a。那么我们可以这样写:
m = b *a;
n = c * a;
所以m和n的最小公倍数就应该是a*b*c啊,那不就是m * n / a,其中m和n是已知的,而a就是那个需要求解的最大公约数。所以就有了下面的代码,
[cpp] view
plaincopy
int GetMinCommonMultiple(int m, int n)
{
assert(m && n);
return m * n / gcd(m, n);
}
相关文章推荐
- 如何求解最大公约数和最小公倍数
- C++实现求解最大公约数和最小公倍数
- 最大公约数和最小公倍数问题(天梯)
- C++/C经典算法百题--(39-42)年龄几何,三色球问题,两个正整数的最大公约数和(GCD)和最小公倍数(LCM)
- 8.最大公约数和最小公倍数问题
- java中关于整数的一些操作问题--最大公约数、最小公倍数、质数等
- Java中循环问题:for,while,do-while 各种练习(打印上三角、下三角、输出100以内的素数、水仙花、最大公约、最小公倍数、分解质因数、兔子问题)
- 8.最大公约数和最小公倍数问题
- 最大公约数和最小公倍数问题
- 最大公约数和最小公倍数问题
- 最大公约数和最小公倍数问题
- 最大公约数和最小公倍数问题
- 最大公约数,最小公倍数,素数等问题
- 8.最大公约数和最小公倍数问题
- 最小公倍数和最大公约数问题
- 【数论】【最大公约数】【枚举约数】CODEVS 1012 最大公约数和最小公倍数问题 2001年NOIP全国联赛普及组
- 最小公倍数和最大公约数的求解
- 最大公约数和最小公倍数问题
- 8.最大公约数和最小公倍数问题
- 8.最大公约数和最小公倍数问题