Gcd && LCM 问题--------stien算法和辗转相除法

题目描述及问题分析

		已知两位整数m,n,求其最大公约数和最小公倍数

*(解决这类问题的方法有很多,在这里重点总结效率最高的栈转相除法)

  				以m = 12 , n = 16为例
   						m 	 n 
   						12 	16
  第一步:   ---------------------------gcd(12,16);
  	 t = 12%16 == 12;
  	 m = 16;
     n = 12;	
     
  第二步:  --------------------------gcd(16,12)  ---->gcd(16,12%16);
  	t = 16%12 == 4; 			
  	m = 12;
    n = 4;

  第三步: -------------------------gcd(12,4) ------->gcd(12,16%12);
  	t = 12%4 == 0; 
  	m = 4; 
  	n = 0;

这里演示了辗转相除法的原理,给了几个重要信息:
1.出循环的条件: n == 0;-------------取余分母为0无意义
2.最终的gcd == m;
这个也是递归调用的出递归的条件,即递归出口

循环

int gcd( int m , int n) {
	while(n) { //为分母
		int t = m%n;
		m = n;
		n = t;
	}
	return m;  //为分子
}

递归

基于辗转相除的特性,我们可以用递归来解决该问题,简单明了 每一步的处理形式都形同,通过规律可知下一步为:gcd(n,m%n)

int gcd( int m , int n ) {
	return n == 0 ? m : gcd( n , m%n );
	//当n == 0 时,递归结束

主函数:

int main() {
	int m,n;
	cin>>m>>n;
	cout<<"最大公约数: "<<gcd(m,n)<<endl;
	cout<<"最小公倍数: "<<(m*n)/gcd(m,n);
	//最大公约数 * 最小公倍数 == m*n;
	return 0;
}

c++中STL中的__gcd(x,y)函数

#include <bits/stdc++.h>
#include <algorithm>//其函数在这个库里面
using namespace std;
int main() {
	int m,n;
	cin>>m>>n;
	cout<<"最大公约数: "<<__gcd(m,n)<<endl;  //两者的数据类型要相同
	cout<<"最小公倍数: "<<(m*n)/__gcd(m,n);
	return 0;
}

位运算

inline int gcd(int a,int b) 
{    
    if(b) while((a%=b) && (b%=a));    
    return a+b;
}

stein算法(处理超过64位的整数)

代码:

int SteinGcd ( int m, int n ) {
	if(!min(n,n))//最小为0输出最大,递归出口,两种方法都是这个条件
		return max(m,n);//多种情况,这个单独写,用三目运算符反而麻烦
	if(!(m%2) && !(n%2))//m和n都为偶数
		return 2*SteinGcd (m/2,n/2);
	if(!(m%2))//m为偶数
		return SteinGcd (m/2,n);
	if(!(n%2))//n为偶数
		return SteinGcd (m,n/2);
	return SteinGcd ((m+n)/2,abs(m-n)/2);//m和n都为奇数
}

①递归出口: 最小为0输出最大
分析m,n均不为0的情况
②$m_n,n_n$为偶数
$m_{n+1}=m_n/2$
$n_{n+1}=n_n/2$
SteinGcd是返回的数据 $SteinGc{d}_{n+1}=SteinGc{d}_n\ast 2$
③$m_n,$为偶数
$m_{n+1}=m_n/2$
$n_{n+1}=n_n$
SteinGcd是返回的数据 $SteinGc{d}_{n+1}=SteinGc{d}_n$
④$n_n,$为偶数
$m_{n+1}=m_n$
$n_{n+1}=n_n/2$
SteinGcd是返回的数据 $SteinGc{d}_{n+1}=SteinGc{d}_n$
⑤②$m_n,n_n$为奇数
$m_{n+1}=|n_n-m_n|$
$n_{n+1}=min(n_n,m_n)$
SteinGcd是返回的数据 $SteinGc{d}_{n+1}=SteinGc{d}_n$