欧几里得定理 及 拓展欧几里得 小结

【欧几里得定理】

定义：

• $\gcd (a,b)=\gcd (b,a\mathrm{m}\mathrm{o}\mathrm{d}b)$
  
  

应用：

故得到 求解gcd的辗转相除法

模板：

int gcd(int a,int b)
{
	return b?gcd(b,a%b):a;
}

int lcm(int a,int b)
{
    return a/gcd(a,b)*b;      //这样写防止溢出
}

---

【拓展欧几里得】

定义：

• 已知$a,b$求解一组$x,y$，使它们满足贝祖等式： $ax+by=\gcd (a,b)=d$（解一定存在）

推导：

已知：$a{x}_1+b{y}_1=\gcd (a,b)=d$ —— ①
　　　$b{x}_2+(a\mathrm{m}\mathrm{o}\mathrm{d}b)y_2=\gcd (b,a\mathrm{m}\mathrm{o}\mathrm{d}b)=d$ —— ②

由欧几里得定理$\gcd (a,b)=\gcd (b,a\mathrm{m}\mathrm{o}\mathrm{d}b)$，
可得：$a{x}_1+b{y}_1=b{x}_2+(a\mathrm{m}\mathrm{o}\mathrm{d}b)y_2$
$⟹a{x}_1+b{y}_1=b{x}_2+(a-[\frac{a}{b}]b)y_2$
$⟹a{x}_1+b{y}_1=a{y}_2+b(x_2-[\frac{a}{b}]y_2)$

所以解得：$x_1=y_2{y}_1=x_2-[\frac{a}{b}]y_2$

当我们已知$x_2,y_2$时，便可以往回推导得到$x_1,y_1$。可以发现，当$b=0$（即辗转相除结束）时，有$d=a,x=1$，所以可以得到一组$x,y$：$x=1,y=0$（$y$可以为任意值，通常取0），然后以此解回推，即可得到一组特解$x_0,y_0$。

模板：

void exgcd(int a,int b,int &d,int &x,int &y)
{
	if(b==0)
	{
		x=1;
		y=0;
		d=a;
	}
	else
	{
		exgcd(b,a%b,d,y,x);
		y-=x*(a/b)
	}
}

从特解到通解：

已解得一组特解$x_0,y_0$，满足： $a{x}_0+b{y}_0=\gcd (a,b)=d$

1. 对于：$ax+by=\gcd (a,b)=d$

   其通解为：$X=x_0+k\frac{b}{d}Y=y_0-k\frac{a}{d}(k\in Z)$
   
   

2. 对于：$ax+by=c$

   当且仅当$d|c$（ 即$\gcd (a,b)$整除$c$，或者说$c$被$\gcd (a,b)$整除 ）时，有解，

   其通解为：$X=\frac{c}{d}{x}_0+k\frac{b}{d}Y=\frac{c}{d}{y}_0-k\frac{a}{d}(k\in Z)$

应用：

①求解方程（形如$ax+by=c$）的特解/通解

例如要求求出平面内，直线$ax+by=c$上在某些范围内的整数点。

②※ 求解线性同余方程：$ax=c(\mathrm{m}\mathrm{o}\mathrm{d}b)$

详见：线性同余方程的求解 及 模板题

③求单个逆元

求$a$在$\mathrm{m}\mathrm{o}\mathrm{d}p$下的逆元，设$a$的逆元为$x$

即求解线性同余方程：$ax=1(\mathrm{m}\mathrm{o}\mathrm{d}p)$