欧几里得算法

欧几里得算法

介绍

概念

  • 欧几里得算法又称辗转相除法,是指用于计算两个非负整数a,b的最大公约数。应用领域有数学和计算机两个方面。

公式

  • 计算公式gcd(a,b) = gcd(b,a mod b)

GCD、LCM

  • 注:最大公因数求法还有辗转相减法等方法

GCD(最大公约数)

//求a,b最大公约数
int find_gcd(int a, int b)
{
	return b==0?a:gcd(b,a%b);
}

LCM(最小公倍数)

//求a,b的最小公倍数
int find_lcm(int a, int b)
{
	return a/gcd(a,b)*b;
}

扩展欧几里得算法

求方程ax+by=gcd(a,b)特解

  • 当方程满足ax+by=gcd(a,b)时,可用扩展欧几里得算法求特解(x0,y0)
/*求方程ax+by=gcd(a,b)的一组特解*/
void extend_fun(int a, int b, int &x, int &y)
{
	if(b==0)
	{
		x=1;
		y=0;
		return;
	}
	extend_fun(b,a%b,x,y);
	int temp = x;
	x = y;
	y = temp - (a/b)*y;
}

求任意方程ax+by=n的一个整数解

  • 判断方程是否有整数解,有解的条件是n可以被gcd(a,b)整除
  • 求处ax+by=gcd(a,b)的一个特解(x0,y0)
  • ax0+by0=gcd(a,b)两边同时乘n/gcd(a,b)
  • 则得到ax+by=n的一个整数解如下

void extend_fun2(int a, int b, int n, int&x, int &y)
{
	int gcd_num = find_gcd(a,b);
	if(n%gcd_num!=0) return;	//无整数解,返回
	
	int s=0,t=0;
	extend_fun(a,b,s,t);
	
	x = s*n/gcd_num;
	y = t*n/gcd_num;
}
posted @ 2022-04-06 14:33  CairBin  阅读(24)  评论(0编辑  收藏  举报