【原创】更相减损术 stein算法 欧几里得算法 拓展欧几里得算法 扩展欧几里得算法 逆元的计算与筛法 解模线性方程

欧几里得

说在前面

数论学复习 Part 6。
然后再来一章CRT和组合数，就飞往概率，以此为跳板去向DP。
计划很美啊你。

P.S. 这么说来拉格朗日插值可以说是数论学复习的Part 0了啊。

有一个*龙门粗口*的*这就是祖安人打招呼的方式*博主在以前发了一篇*儒雅随和*的博客，也是讲这个的，虽然写的很*祖安粗口*，但证明还是可堪一看的，（反正是借鉴来的）。

为什么我还会记得多年前的我写博客的亚子呢？

更相减损术 辗转相减法

可半者半之，不可半者，副置分母、子之数，以少减多，更相减损，求其等也。以等数约之。

$gcd(a,b)=gcd(b,a-b)$。

int gcd(int a,int b)
{
	if(a<b) swap(a,b);
	return b?gcd(b,a-b):a;
}

再加上“可半者半之”。

int gcd(int a,int b)
{
	if(a<b) swap(a,b);
	if(a&1) a>>=1;
	if(b&1) b>>=1;
	return b?gcd(b,a-b):a;
}

$$\begin{gathered} P.S.更相损减术本来是用来约分的，所以“可半者半之”这句话被狭义理解为分子分母都是偶数时约掉2。 \\ 而我纯路人，看到可半者半之就直接断章取义了。 \end{gathered}$$

这样的话，stein算法和辗转相减法一模一样。$\xcancel{(中国数学领先外国1700年！)}$

原理是一个数$d$它整除$a$和$b$的话，它就整除$a,b$的任意线性组合即$d\mid pa+qb(p,q\in Z)$。

证明私以为完全不需要了。

欧几里得算法 辗转相除法

$gcd(a,b)=gcd(b,a\%b)$
模运算是什么？就是多次反复地减直到减不动。
就是多次调用上一个算法。
正确性就更不用赘述了。

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

贝祖（裴蜀）定理

对于任意整数$a,b$，一定存在整数$x,y$，使得$ax+by=gcd(a,b)$。
特别的，如果$a,b$互质，就一定存在整数$x,y$满足$ax+by=1$。

设$gcd(a,b)=d$，$a=pd,b=qd$，则$ax+by=gcd(a,b)⟺px+qy=1$。易知$gcd(p,q)=1$。
也就是说把所有一般命题转到了“特别的”命题上来。
怎么证明存在性问题？简单啊，给出一组$(x,y)$不就好了吗？
原式$⟺px\equiv 1\mathrm{m}\mathrm{o}\mathrm{d}q$
就是求p关于q的一个逆元，这个逆元的存在性我已经在之前不用贝祖证明出来了（1.1），现在就不再赘述一次。

P.S.换一个等价命题应该会证得舒服点，比如$ac\equiv bc\mathrm{m}\mathrm{o}\mathrm{d}p,(c,p)=1\to a\equiv b\mathrm{m}\mathrm{o}\mathrm{d}p$。

拓展欧几里得算法

拓展欧几里得算法能够在辗转相除的同时求出一组$x,y$使得$ax+by=gcd(a,b)$。

怎么在辗转相除的时候求呢？我们要看$gcd(b,a\%b)$是吧？
我们设已经辗转相除求出了$b{x}_0+(a\%b)y_0=gcd(b,a\%b)$，然后试图用$x_0,y_0$递推出$x,y$。

因为$gcd(a,b)=gcd(b,a\%b)$，所以有$ax+by=b{x}_0+(a\%b)y_0$。设$a=kb+r$，则$原方程⟺(kb+r)x+by=b{x}_0+r{y}_0⟺b(kx+y)+rx=b{x}_0+r{y}_0$。
然后随便多项式恒等定理一下就得到$\{\begin{array}{l}kx+y=x_0\\ x=y_0\end{array}\Rightarrow \{\begin{array}{l}x=y_0\\ y=x_0-k{y}_0=x_0-\lfloor \frac{a}{b}\rfloor y_0\end{array}$。

然后到最后$b=0,a=gcd$的时候，有$x_k=1,y_k=0$。

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

或者

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

逆元

定义

除法怎么取模？
小数分数怎么取模？整除、整数的概念都没有了，你怎么取模？
那怎么除法呢？
怎么规避小数分数呢？
$\frac{5}{2}=2.5$的话，$5=2\ast 2.5$。
所以乘法就行了。

我们如果说一个数$a$关于$p$的逆元是$a^{-1}$的话，就等价于定义$\frac{x}{a}\mathrm{m}\mathrm{o}\mathrm{d}p=x{a}^{-1}\mathrm{m}\mathrm{o}\mathrm{d}p$。
要求$a\times a^{-1}\equiv 1\mathrm{m}\mathrm{o}\mathrm{d}p$。

这个逆元存在的前提上文也有说，$a,p$要互质。

求法

法一

如果$p$是质数的话，我们就有费马小定理，$a^{p-1}\equiv 1\mathrm{m}\mathrm{o}\mathrm{d}p$。此时$a$的逆元显然就是$a\times a^{p-2}\equiv 1\mathrm{m}\mathrm{o}\mathrm{d}p$。
这个东西据说还叫蒙哥马利快速幂模。

法二

另外一种方法就是用exgcd。
exgcd(a,b,&x,&y)得到的x就是a关于b的逆元，y就是b关于a的逆元。

筛法

给一个数p，递推地求逆元，求出一个范围内所有数关于p的逆元。
那么p肯定是质数或者相对是质数了。

现在正在求$a$关于$p$的逆元，设$p=ka+r$，则有$ka+r\equiv 0(\mathrm{m}\mathrm{o}\mathrm{d}p)$，设$a$的逆元为$a^{-1}$，$r$的逆元是$r^{-1}$，在两边同时乘以$a^{-1}{r}^{-1}$，则有$k{r}^{-1}+a^{-1}\equiv 0(\mathrm{m}\mathrm{o}\mathrm{d}p)⟺a^{-1}\equiv -k{r}^{-1}\equiv -\lfloor \frac{p}{a}\rfloor (p\mathrm{m}\mathrm{o}\mathrm{d}a{)}^{-1}(\mathrm{m}\mathrm{o}\mathrm{d}p)$。

inv[0]=inv[1]=1;
for(int i=2;i<=N;i++) inv[i]=1ll*(mod-mod/i)*inv[mod%i]%mod;

解模线性方程

$ax+by=c$。

首先如果$gcd(a,b)\nmid c$的话，肯定无解，因为左边必然是$gcd(a,b)$的倍数。
如果$gcd(a,b)\mid c$的话，设$c=k\times gcd(a,b)$，那么我们跑一次扩欧，得到$a{x}_0+b{y}_0=gcd(a,b)$的一组解$(x_0,y_0)$，然后再乘上$k$就搞定了。
$\{\begin{array}{l}x=k{x}_0\\ y=k{y}_0\end{array}$

如果想要求最小正整数解的话，则就$(x\%p+p)\%p$就是了。

$ax\equiv b\mathrm{m}\mathrm{o}\mathrm{d}p$。

和上面那个差不多。