求解最大公約數和最小公倍數問題
來源:程序員人生 發布時間:2015-04-17 08:34:24 閱讀次數:5441次
思路:
最大公約數問題也是1個非常典型的遞歸算法的利用。每次遞歸使得原來求兩個大數之間的公約數轉變成求兩個略微小點的數之間的公約數,要求轉換的進程要保證不會改變公約數的值。這就要看其中轉換的原理了。
原理從《幾何本來》中得出--展轉相除。假定f(x, y) 表示x,y的最大公約數是g,而k = x/y,b= x%y,則g必能整出b。由于x = ky + b,b = x - ky,b/g = (x-ky)/g1定為整數,所以必有g整除b。
以下所示:
f(42, 30) = f(30, 12) = f(12, 6)= f(6, 0) = 6
代碼以下:
-
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;
}
書中還引申出展轉相減法,原理跟上面所述差不多。
代碼以下:
-
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;
}
最后書中提到1種相除和相減相結合的方法,保證了不做過量冗余的步驟。這就是把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)
此解法結合了上述兩種方法:第1種方法遞歸次數相當少但每次取模挺耗時(到底耗時情況如何,我沒有具體了解過,但肯定比加減法耗時些),而第2種遞歸可能會相當多,但每次運算都是減法運算, 很快。綜合二者,用2為基數,可以保證遞歸次數不那末多,同時運算也快。
代碼以下:
-
int gcd(int x , int y)
-
{
-
if(x < y)
-
return gcd(y , x) ;
-
if(y == 0)
-
return x ;
-
if(isEven(x))
-
{
-
if(isEven(y))
-
return gcd(x - y, y) ;
-
else
-
return gcd(x , y >>1) ;
-
}
-
else
-
{
-
if(isEven(y))
-
return gcd(x >> 1, y) ;
-
else
-
return 2 * gcd(x >> 1, y >> 1) ;
-
-
}
-
-
}
求解最小公倍數:
從題面上看,好像我們需要求解的是兩個題目,但其實就是1個題目。那就是求最大公約數?為何呢?我們可以假想這兩個數m和n,假定m和n的最大公約數是a。那末我們可以這樣寫:
m = b *a;
n = c * a;
所以m和n的最小公倍數就應當是a*b*c啊,那不就是m * n / a,其中m和n是已知的,而a就是那個需要求解的最大公約數。所以就有了下面的代碼,
-
int GetMinCommonMultiple(int m, int n)
-
{
-
assert(m && n);
-
-
return m * n / gcd(m, n);
-
}
生活不易,碼農辛苦
如果您覺得本網站對您的學習有所幫助,可以手機掃描二維碼進行捐贈