[C++]数论（二）

往期回顾

数论（一）

裴蜀定理

内容

裴蜀定理是指方程$ax+by=m$当且仅当m是a和b的最大公因数d的倍数时，方程有整数解
特别的，只有a和b互质时，$ax+by=1$才有整数解

证明

设存在$x,y$，满足$ax+by$为正数最小值$d$。任选$m,n$不等于$x,y$，此时$am+bn$为$e$，则$e\ge d$。
$e=p\times d+r$（$p$为商，$r$是余数），则$r＜d$
那么$r=e-p\times d=am+bn-p\times (ax+by)=(m-px)\times a+(n-py)\times b$
则有$m-px$和$n-py$满足比d更小的值$r$，$r$只能为0
所以$e$是$d$的倍数
设$gcd(a,b)=l，a=sl，b=tl$，则$d=l(sx+ty)$，所以$d$一定是$gcd(a,b)$的倍数
综上可知，$d$为$gcd(a,b)$的倍数，$e$为$d$的倍数

构造

构造有很多，而扩展欧几里得算法的解就是其中之一。
详情看数论（一），这里不再构造

欧拉定理

欧拉函数

欧拉函数，即$\phi (n)$，指1~$n$中与$n$互质的数的个数

欧拉函数的计算

$\phi (1)=1$
如果$i$是质数，那么$\phi (i)=i-1$
如果$i$是合数，那么设$i={\prod }_{i=1}^m{p}_i^{k_i}$且$p_i$为质数，则$\phi (i)=n\times {\prod }_{i=1}^m\frac{p_i-1}{p_i}$
如果只求一个数的欧拉函数，可以在分解质因数中求
如果求线性的欧拉函数，则可以在线性筛中求

欧拉定理

若$gcd(a,m)=1$，那么$a^{\phi (m)}\equiv 1(\mathrm{mod}m)$

乘法逆元

我们知道，$a\times a^{-1}\equiv 1(\mathrm{mod}b)$
那么如果给定$a,b$，怎么求$x=a^{-1}$呢

扩展欧几里得算法求逆元

我们知道，扩展欧几里得算法求的是$ax+by=gcd(a,b)$的解
这和乘法逆元又有什么关系呢？
我们先来看看条件吧。
$ax\equiv 1(\mathrm{mod}b)$一定有解，说明了$ax+by=gcd(a,b)\equiv 1(\mathrm{mod}b)$
而$1\le gcd(a,b)\le b$，因此$gcd(a,b)=1$
因此，我们用扩展欧几里得算法可以求出$ax+by=1$的解
在这种情况下，显然$ax\equiv 1(\mathrm{mod}b)$
但别忘了现在的$x$可能不是正整数，也不一定是最小的解，因此我们要的答案是$(x\%b+b)\%b$
懂了吗？
奉上代码一份：

#include<iostream>
using namespace std;
long long x,y;
void gcd(long long a,long long b){
	if(!b){
		x=1,y=0;
		return ;
	}
	gcd(b,a%b);
	long long t=x;
	x=y,y=t-a/b*y; 
}
int main(){
	long long a,b;
	cin>>a>>b;
	gcd(a,b);
	x=(x%b+b)%b;
	cout<<x;
	return 0;
}

噢，对了，大家还可以顺便把这题给切了，刷个绿题谁不愿意呢？

欧拉定理求逆元

知道我这次为什么必须要讲欧拉定理了吧。
As we all know，$a^{\phi (b)}\equiv 1(\mathrm{mod}b)$
$\therefore a\times a^{\phi (b)-1}\equiv 1(\mathrm{mod}b)$
因此，我们要求的乘法逆元就是$a^{\phi (b)-1}\%b$
至于怎么求，快速幂嘛，对吧，反正我也懒得写了嘿嘿
当然如果$b$是质数，我们用费马小定理也是得出一样的结论的

线性求逆元

如果要求1~n的逆元了，怎么样才能最快？
我想，最缓(xun)慢(su)的做法就是暴力了。
我们可以来找找逆元之间的关系
$\lfloor \frac{p}{a}\rfloor \times a=p-(p\%a)$
$\lfloor \frac{p}{a}\rfloor \times a\equiv -(p\%a)(\mathrm{mod}p)$
$a\equiv \frac{-(p\%a)}{\lfloor \frac{p}{a}\rfloor }(\mathrm{mod}p)$
$a^{-1}\equiv \frac{\lfloor \frac{p}{a}\rfloor }{-(p\%a)}(\mathrm{mod}p)$
$a^{-1}\equiv -\lfloor \frac{p}{a}\rfloor \times (p\%a{)}^{-1}(\mathrm{mod}p)$
所以我们就可以通过这个来计算线性上的逆元了

inv[1]=1;
for(int i=2;i<=n;i++){
	inv[i]=(long long)(b-b/i)*inv[b%i]%b
}

今天要不就到这吧。