数论基础——operator_

1.裴蜀定理

对于不定方程$ax+by=c,$有整数解的的充要条件是$\gcd (a,b)|c$。

也可写成$ax+by\equiv 0(\mathrm{mod}\gcd (a,b))$。

证明：

1.证充分：令$d=\gcd (a,b),d|a$且$d|b,$设$a=t_{1}d,b=t_{2}d,$那么$c=t_{1}xd+t_{2}yd,d|c$。

2.证必要：设$ax+by$最小正值为$S,$令$a=pS+m,m=a\mathrm{mod}s(S>m\ge 0)$。

$m=a-pS=a-p(a{x}_0+b{y}_0)=a-pa{x}_0-pb{y}_0=a(1-p{x}_0)+b(-p{y}_0)$

发现$m$是$ax+by$的某值,所以$m=0$才不矛盾

所以$S|a,$同理$S|b,$故$S|\gcd (a,b),S=\gcd (a,b)$。

几乎是模板：普及/提高-【模板】裴蜀定理

当$n=2,$显然有$S=\gcd (A_1,A_2);$

当$n=3,S=A_1{x}_1+A_2{x}_2+A_3{x}_3,$

令$A_1{x}_1+A_2{x}_2=y_1\ast \gcd (A_1,A_2)$

则$S=\gcd (A_1,A_2)y_1+A_3{x}_3$

再来一次裴蜀定理,就有$S=\gcd (A_1,A_2,A_3)$。

可以由归纳法推广至正整数域：$S=\gcd (A_1,A_2\cdots A_n)$。

2.拓展欧几里得

用于求解形如$ax+by=\gcd (a,b)$的方程问题。

证明：

设原方程为$aX+bY=\gcd (a,b)$

令$bx+(a\%b)y=\gcd (b,a\%b)$

有$a%b=a- \lfloor{{a}\over{b}}\rfloor*b $

所以$bx+(a-\lfloor \frac{a}{b}\rfloor \ast b)y=aX+bY$

$ay+b(x-\lfloor \frac{a}{b}\rfloor y)=aX+bY$

$X=y,Y=x-\lfloor \frac{a}{b}\rfloor y$

最后考虑一下边界即可。

几乎是模板：普及+/提高[NOIP2012 提高组] 同余方程

用$exgcd()$先求得$ax+by=\gcd (a,b)$的一组解$(x_0,y_0)$

设另一组解为$x_1,y_1$。

则由$a{x}_0+b{y}_0=\gcd (a,b)$和$a{x}_1+b{y}_1=\gcd (a,b)$

知$a(x_0-x_1)+b(y_0-y_1)=0$

$a(x_0-x_1)=-b(y_0-y_1)$

设$x_1=x_0+i,y_1=y_0+j$

就有$ai=-bj,i=-\lfloor \frac{a}{b}\rfloor \ast j$

因为$i$为最小整数$,a|j\ast b,j$的最小值为$a/\gcd (a,b)$

同理$i$的最小值为$b/\gcd (a,b)$

所以可得到原方程通解$\{\begin{array}{l}x=x_0+t\ast b/\gcd (a,b)\\ y=y_0+t\ast a/\gcd (a,b)\end{array}$。

不过本题只要输出$(x+b/\gcd (a,b))\%(b/\gcd (a,b))$。

$PS.$ 可以推广得到任意方程的通解

习题：

提高+/省选-青蛙的约会

3.乘法逆元

在一般情况下,若$a\ast x=1$则$b$为$a$的倒数

对应的，乘法逆元的定义：若$a\ast x\equiv 1(\mathrm{mod}p)$,则称$x$为$a$在模$p$意义下的乘法逆元$,$记作$a^{-1}(\mathrm{mod}p)$

显然$,$若$p|a$或$a=0$则没有逆元$,$同时逆元也可能不止一个。

求解方法：

1.拓展欧几里得

由定义：$ax\equiv 1(\mathrm{mod}p)$

则可设$ax+py=1,$解方程即可。

2.费马小定理

不会证qwq

若$p$为质数且$a$为整数$,$则$a$在模$p$意义下的乘法逆元为$a^{p-2}\%p,$一般快速幂求解。

3.递推

令$inv[i]$为$i$在模$p$意义下的乘法逆元$,$显然有$inv[1]=1$

设已知$1,2,\cdots ,i-1$的逆元

令$k=\lfloor \frac{p}{i}\rfloor ,r=p\%i,$则$p=k\ast i+r(1\le r<i)$

$k\ast i+r\equiv 0(\mathrm{mod}p)$

$k\ast i\ast r^{-1}+1\equiv 0(\mathrm{mod}p)$

$i^{-1}\equiv -k\ast r^{-1}(\mathrm{mod}p)$

最小正数解:$i^{-1}\equiv -\lfloor \frac{p}{i}\rfloor \ast (p\mathrm{mod}i{)}^{-1}(\mathrm{mod}p)$

4.阶乘逆推

原理是$\frac{1}{k}=\frac{(k-1)!}{k!}$和$\frac{1}{(k-1)!}=\frac{k}{k!}$

令$f[i]=(i!)\%p,d[i]=(i!{)}^{-1}$

显然$f[]$循环一遍就行$,d[i]$的话就先用$exgcd()$或费马小定理求出$d[n],$再由$d[i-1]=i\ast d[i]\%p$倒序求即可

那么$inv[k]=\frac{(k-1)!}{k!}=(k-1)!\ast \frac{1}{k!}\equiv f[k-1]\ast d[k](\mathrm{mod}p)$

相较于递推法$,$本方法可以求解任意$n$个数的逆元

习题

普及+/提高【模板】有理数取余

普及/提高-【模板】乘法逆元

虽然但是,个人感觉第一题比第二题简单

4.中国剩余定理(CRT)

有如下的一组同余方程组：

$$x\equiv a_1(\mathrm{mod}{m}_1)$$

$$x\equiv a_2(\mathrm{mod}{m}_2)$$

$$⋮$$

$$x\equiv a_n(\mathrm{mod}{m}_n)$$

满足 $\gcd (m_1,m_2,\dots ,m_n)=1$

解法：

$1.$ 令 $M=m_1\ast m_2\ast \cdots \ast m_n$ ， $M_i=M/m_i$ 。则有 $M_i\perp m_i$ 。

$2.$ 用 $exgcd$ 求出 $M_i$ 关于 $m_i$ 的逆元 $t_i$ 。

$3.$ 求得同余方程组的特解 $x_0={\sum }_{i=1}^{i\le n}{a}_i\ast M_i\ast t_i$ 。

证明：

显然，对于 $\forall i$ ， $\forall j\ne i$ ， $m_i|M_j$

$\therefore x_0\equiv a_i\ast M_i\ast t_i\equiv a_i(\mathrm{mod}{m}_i)$

$\therefore x_0$ 确实是一个特解。

$4.$ 求得同余方程组的通解 $x=x_0+k\cdot M$

证明：

使用反证法，假设除了 $x=x_0+k\cdot M$ 还有其他解 $x^{\prime }$ ，就有 $x\equiv a_i(\mathrm{mod}{m}_i)$ 和 $x^{\prime }\equiv a_i(\mathrm{mod}{m}_i)$

两式相减， $x-x^{\prime }\equiv 0(\mathrm{mod}{m}_i)$ ， $m_i|x-x^{\prime }$

$\therefore$ 对于 $\forall i$ ， 有 $m_i|x-x^{\prime }$ 。又 $m_i$ 两两互质，那么 $\mathrm{lcm}(m_1,m_2,\dots ,m_n)=M$ 。

$\therefore M|x-x^{\prime }$ ， $x^{\prime }=x+c\cdot M=x_0+k\cdot M$

得证。

习题：

P1495 【模板】中国剩余定理（CRT）/ 曹冲养猪

P3868 [TJOI2009] 猜数字

5.拓展中国剩余定理(EXCRT)

同上的一组同余方程组，但是 $m_i$ 不再两两互质。

解法：

考虑前 $2$ 个方程： $s\equiv a_1(\mathrm{mod}{m}_1)$ 和 $s\equiv a_2(\mathrm{mod}{m}_2)$

$\therefore s=m_1\ast x+a_1=m_2\ast y+a_2$

得 $m_1\ast x+m_2\ast (-y)=a_2-a_1$

之所以可以把 $y$ 变成 $-y$ 是因为我们只关心 $x$ 的取值。

用 $exgcd$ 求出最小正整数解 $x_0$ ，那么前两个方程的通解为 $x=x_0+k\cdot \frac{m_2}{\gcd (m_1,m_2)}$

$\therefore s$ 的通解 $s=m_1\ast x_0+k\cdot \frac{m_1{m}_2}{\gcd (m_1,m_2)}+a_1=m_1\cdot x_0+a_1+k\cdot \mathrm{lcm}(m_1,m_2)$

所以前 $2$ 个方程可以合并为新的同余方程 $s\equiv m_1\ast x_0+a_1(\mathrm{mod}\mathrm{lcm}(m_1,m_2))$

那么就可以不断合并，最终得到 $s\equiv A(\mathrm{mod}\mathrm{lcm}(m_1,m_2,\dots ,m_n))$ 。

习题：

P4777 【模板】扩展中国剩余定理（EXCRT）

6.欧拉函数

前置知识：

$3$ 种普通筛法（暴力、埃氏筛、欧拉筛）

积性函数：

数论函数：定义域为正整数域，值域为复数域的函数

积性函数：对于一个数论函数，如果 $\forall x,y$ ，且 $x\perp y$ ，都有 $f(xy)=f(x)\cdot f(y)$ ，那么 $f$ 是一个积性函数。

完全积性函数：对于一个数论函数，如果 $\forall x,y$ ，且 $x,y$ 不需要互质，都有 $f(xy)=f(x)\cdot f(y)$ ，那么 $f$ 是一个完全积性函数。

一些常见的积性函数：

$1.$ $1$ 函数： $f(n)=1$

$2.$ 幂函数： $i{d}_k(n)=n^k$

$3.$ 狄利克雷卷积单位元： $\epsilon (n)=\{\begin{array}{ll}1,& n=1\\ 0,& n\ne 1\end{array}$

$4.$ 因子幂和函数： ${\sigma }_k(n)={\sum }_{i|n}{i}^k$

$5.$ 欧拉函数： $\phi (n)={\sum }_{i=1}^{i\le n}1[n\perp i]$

$6.$ 莫比乌斯函数： $\mu (n)=\{\begin{array}{ll}0,& \dots 1\\ (-1{)}^k,& \dots 2\end{array}$

注：$1:n$ 无平方因子。 $2:n$ 有平方因子， $k$ 为 $n$ 的质因子个数。

$7.$ 最大公因数： $gc{d}_k(n)=\gcd (n,k)$

其实大部分今天都没卵用

欧拉函数的求解：

若设 $n=p_1^{k_1}{\cdot }_2^{k_2}\cdot \dots \cdot p_r^{k_r}$

则 $\phi (n)=n\cdot (1-\frac{1}{p_1})\cdot (1-\frac{1}{p_2})\cdot \dots \cdot (1-\frac{1}{p_r})$

习题：

P2303 [SDOI2012] Longge 的问题

P2568 GCD

P2158 [SDOI2008] 仪仗队

7.欧拉定理

内容：

设 $a\perp n$ ，那么 $a^{\phi (n)}\equiv 1(\mathrm{mod}n)$ 。

当 $n$ 为质数时，就是费马小定理。

证明：

设 $1$ 到 $n$ 互质的数为 $x_1,x_2,\dots ,x_{\phi (n)}$ 。

令 $P_i=x_i\ast a\mathrm{mod}n$ ，尝试证 $P_i$ 两两不同。

设 $\exists$ $i,j$ ，使 $P_i\equiv P_j(\mathrm{mod}n)$ ，则 $n|a(x_i-x_j)$ 。

考虑到 $a\perp n$ ， $n|(x_i-x_j)$ ，但又 $x_i-x_j\ne 0$ 且 $|x_i-x_j|<n$ ，显然假设不成立。

设 $\gcd (P_i,n)=d$ ，由于 $a\cdot x_i=P_i+kn$ 且 $d|n,d|P_i$ ， $\therefore d|a\cdot x_i$ 。

又 $a\perp n,x_i\perp n$ ， $\therefore d=1,P_i\perp n$ 。

综上， $P_i$ 是 $x_i$ 的另一种排列，所以

$$a\cdot x_1\cdot a\cdot x_2\cdot \dots \cdot a\cdot x_{\phi (n)}\equiv P_1\cdot P_2\cdot \dots \cdot P_{\phi (n)}(\mathrm{mod}n)$$

$$a\cdot x_1\cdot a\cdot x_2\cdot \dots \cdot a\cdot x_{\phi (n)}\equiv x_1\cdot x_2\cdot \dots \cdot x_{\phi (n)}(\mathrm{mod}n)$$

$$a^{\phi (n)}\equiv 1(\mathrm{mod}n)$$

8.拓展欧拉定理

内容：

$a^{2\phi (n)}\equiv a^{\phi (n)}(\mathrm{mod}n)$

即若 $b>\phi (n)$ ， $a^b\equiv a^{(b\mathrm{mod}\phi (n)+\phi (n))}(\mathrm{mod}n)$

证明：

设 $a=p_1^{k_1}{\cdot }_2^{k_2}\cdot \dots \cdot p_r^{k_r}$

先证 $p^b\equiv p^{(b\mathrm{mod}\phi (n)+\phi (n))}(\mathrm{mod}n)$

$1.$ $p\perp n$ ：欧拉定理本定理。

$2.$ $p,n$ 不互质：

$p|n$ 且 $n\geqq 2p$ 。设 $n=s\cdot p^r$ ，其中 $s\perp p$ 。

那么由欧拉定理， $p^{\phi (n)}\equiv 1(\mathrm{mod}s)$ ， $\phi (n)=\phi (s)\ast \phi (p^r)$ 。

$$p^{\phi (s)}={(p^{\phi (n)})}^{\phi (p^r)}\equiv 1^{\phi (p^r)}\equiv 1(\mathrm{mod}s)$$

$$p^{\phi (n)+r}\equiv p^r(\mathrm{mod}n)$$

$$p^b\equiv p^{b-r}\ast p^r\equiv p^{b-r}\ast p^{\phi (n)+r}\equiv p^{b+\phi (n)}(\mathrm{mod}n)$$

既然 $p^b\equiv p^{b+\phi (n)}(\mathrm{mod}n)$ ，那么当 $b>\phi (n)$ 时， $p^b\equiv p^{b-\phi (n)}(\mathrm{mod}n)$

$\therefore p^b\equiv p^{b\mathrm{mod}\phi (n)+\phi (n)}(\mathrm{mod}n)$

证毕。

习题：

P5091 【模板】扩展欧拉定理