gcd及exgcd详解

前言

由于复习CSP，我的blog咕咕咕了一个多月。大有洛谷的风范 今天，我来写一下自己学习gcd及exgcd的感受。

详解gcd及exgcd

gcd

什么是gcd

gcd指的是Greatest Common Divisor，就是两个数的最大公因数。

通俗来讲，$\gcd (a,b)$就是最大化一个$Ans$，使得$$\begin{gathered} Ans|a \\ ,Ans|b \end{gathered}$$。

举个例子：

$20$的因数有$1,2,4,5,10,20$。

$15$的因数有$1,3,5,15$。

所以$\gcd (20,15)=5$。

如何求gcd

想必各位在高中都学过辗转相除法和更相减损之术，这里只讲辗转相除法。

辗转相除法

首先不妨设 $x\le y$,则$\gcd (x,y)=\gcd (x,x+y)=\gcd (x,y-x)$，所以$\gcd (x,y)=gcd(y\mathrm{m}\mathrm{o}\mathrm{d}x,x)$，因此可以递归求解。
复杂度证明：因为$y\mathrm{m}\mathrm{o}\mathrm{d}x\le x,x\le y$，所以$y\mathrm{m}\mathrm{o}\mathrm{d}x<y/2$。因此在最坏情况下为$O(nlogn)$。

另外说一句：最坏的情况恰好是斐波那契额数列的相邻两个（细思极恐）

代码

int gcd(int a,int b)
{
	if(!b)
	{
		return a;
	}
	return gcd(b,a%b);
}

exgcd

啥是exgcd

exgcd可以概括成求形如

$$ax+by=c$$

的方程。

当存在$$\gcd (a,b)|c$$时，存在解。

也就是说，这个exgcd可以用来求解方程$ax+by=\gcd (a,b)$

exgcd的证明

我们命$a=b,b=a\mathrm{m}\mathrm{o}\mathrm{d}b$，则有方程$b\times x1+(a\mathrm{m}\mathrm{o}\mathrm{d}b)\times y1=gcd(b,a\mathrm{m}\mathrm{o}\mathrm{d}b)$

又因为$\gcd (a,b)=\gcd (a\mathrm{m}\mathrm{o}\mathrm{d}b)$，且$a\mathrm{m}\mathrm{o}\mathrm{d}b=a-b\times \lfloor a/b\rfloor$

则$b\times x1+(a-b\ast \lfloor a/b\rfloor )\times y1=\gcd (a,b)$

整理得：$a\times y1+b\ast \times (x1-\lfloor a/b\rfloor \times y1)=\gcd (a,b)$

所以原方程中：$x=y1,y=x1-\lfloor a/b\rfloor \times y1$。于是我们只要递归求出$x1,y1$就能求出$x,y$。

代码

void exgcd(ll a, ll b, ll& x, ll& y, ll& c)
{
	if(!b)
	{
		y=0;
		x=1;
		c=a;
		return;
	}
	exgcd(b,a%b,y,x);
	y-=a/b*x;
}