模板——最大公约数

一. 欧几里德算法

即辗转相除法

由 gcd( a , b ) = gcd ( b , a % b )

#include<iostream>
#include<cstring>
#include<cstdio>
using namespace std;

void swap(int &a,int &b){
	int tem=a;
	a=b;
	b=tem;
}

int Gcd(int a,int b){
	if(a<b)
		swap(a,b);
	return b==0?a:Gcd(b,a%b);
}

int main(){
	int a,b;
	while(cin>>a>>b){
		cout<<Gcd(a,b)<<endl; 
	}
	return 0;
}

二. Stein算法

Stein算法是一种计算两个数最大公约数的算法，是针对欧几里德算法在对大整数进行运算时，需要试商导致增加运算时间的缺陷而提出的改进算法。

原理:

  ① 当两个数都为偶数时，gcd( ka , kb ) = k * gcd( a , b ) ，当 k = 2 时,说明两个偶数的最大公约数必然能被 2 整除。
  ② 当一个为奇数，一个为偶数时， gcd( 2 * a , b ) = gcd( a , b )，可以先将偶数除以2。
  ③ 当两个数都为奇数时，gcd( a , b ) = gcd( a-b , b ) 。

#include<iostream>
#include<cstring>
#include<cstdio>
using namespace std;
#define ll long long

void swap(ll &a,ll &b){
	ll tem=a;
	a=b;
	b=tem;
}

ll Stein(ll a,ll b){
	if(a<b)	swap(a,b);
	if(a==b) return a;
	ll r=0;
	while(b>0){
		if(!(a&1)&&!(b&1)){//a为偶数，b为偶数 
			a>>=1;
			b>>=1;
			r++;
		}
		else if(!(a&1)&&(b&1)){//a为偶数，b为奇数 
			a>>=1;
		}
		else if((a&1)&&!(b&1)){//a为奇数，b为偶数 
			b>>=1;
		}
		else if((a&1)&&(b&1)){//a为奇数，b为奇数 
			a=(a-b)>>1;
		}
		if(a<b)	
			swap(a,b);
	}
	return a<<r;
}

int main(){
	ll a,b;
	while(cin>>a>>b){
		cout<<Stein(a,b)<<endl;
	}
	return 0;
}