数论笔记 (20200205)

类欧几里得算法：

求${\sum }_{x=0}^n\lfloor \frac{ax+b}{c}\rfloor$
类欧几里得算法的推导

欧拉定理：

对于$(a,m)=1$，有 $a^{\phi (m)}\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}m)$

扩展欧拉定理：

对于$(a,m)\ne 1$，有$a^b\equiv a^{\min (b,b\mathrm{m}\mathrm{o}\mathrm{d}\phi (m)+\phi (m))}(\mathrm{m}\mathrm{o}\mathrm{d}m)$
等价的表述为对于$b\ge \phi (m)$，有$a^b\equiv a^{b+\phi (m)}(\mathrm{m}\mathrm{o}\mathrm{d}m)$
证明：
首先证明对于质数$p$成立，令$m=s\ast p^r$，且$(s,p)=1$，即将$m$中的$p$提出来。
因为$(s,p)=1$，所以$p^{\phi (s)}\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}s)$
因为$\phi (s)|\phi (m)$，所以$p^{\phi (m)}\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}s)$
根据同余性质，左右以及模数同时乘以$p^r$，有$p^{\phi (m)+r}\equiv p^r(\mathrm{m}\mathrm{o}\mathrm{d}m)$
$\phi (m)=\phi (s)(p-1)p^{r-1}\ge p^{r-1}\ge 2^{r-1}\ge r$，即$\phi (m)\ge r$
所以当$b\ge r时$，$p^b\equiv p^{\phi (m)+b}(\mathrm{m}\mathrm{o}\mathrm{d}m)$，而$r\le \phi (m)$，得证。

然后证明对于$p^k$成立：
当$b\ge \phi (m)$时，$(p^k{)}^b=p^{kb}\equiv p^{kb+k\phi (m)}\equiv (p^k{)}^{b+\phi (m)}$，得证。

对于任意的$a$，将其分解为$p^k$后即可得证。

---

例题：上帝与集合的正确用法

求 $2^{2^{2^{...}}}\mathrm{m}\mathrm{o}\mathrm{d}p$，$p\le 10^9$（不一定为质数）

结论：$\phi (\phi (n))\le \frac{1}{2}n$
证明：$\phi (n)=n\ast \prod \frac{p_i-1}{p_i}$，如果$n$为奇数，那么$\phi (n)$一定为偶数；如果$n$为偶数，那么$\phi (n)\le \frac{1}{2}n$。

应用拓展欧拉定理：$2^{2^{2^{...}}}\equiv 2^{2^{2^{2^{...}}}\mathrm{m}\mathrm{o}\mathrm{d}\phi (p)+\phi (p)}(\mathrm{m}\mathrm{o}\mathrm{d}p)$，递归，由上面的结论递归的层数为$O(logp)$层。

---

Pollard_Rho和Miller_Rabin

中国剩余定理：

有$k$个方程组：$X\equiv x_i(\mathrm{m}\mathrm{o}\mathrm{d}{m}_i)$，$m_i$两两互质，求$X$
令$M=\prod m_i$，$M_i=\frac{M}{m_i}$，$M_i^{-1}$为$M_i$模$m_i$的逆元。$$X=\sum _{i=1}^k{x}_i\ast M_i\ast M_i^{-1}(\mathrm{m}\mathrm{o}\mathrm{d}M)$$

扩展中国剩余定理：

模数不一定互质，将方程依次合并：
$\{\begin{array}{l}x\equiv a(\mathrm{m}\mathrm{o}\mathrm{d}b)\\ x\equiv c(\mathrm{m}\mathrm{o}\mathrm{d}d)\end{array}$
可知 $b\cdot t+a\equiv c(\mathrm{m}\mathrm{o}\mathrm{d}d)$
用$\text{exgcd}$求出$t\equiv \frac{c-a}{(b,d)}\ast (\frac{b}{(b,d)}{)}^{-1}(\mathrm{m}\mathrm{o}\mathrm{d}\frac{d}{(b,d)})$
代回得$x\equiv b\cdot t+a(\mathrm{m}\mathrm{o}\mathrm{d}[b,d])$
易验证$x$满足原方程，且方程的通解为$x+k[b,d]$

两个定理都旨在求出一个特解，易发现在模数范围内只有一个解，通解形式就是特解±模数的倍数。

例题：NOI2018 屠龙勇士

例题：CF#338D GCD Table

---

阶：

BSGS：

BSGS算法

多项式复合：

直接暴算$g(x{)}^i$是$O(nmlogm)$的，乘系数加起来是$O(nm)$的，复杂度瓶颈在算$g(x{)}^i$。
可以预处理$g(x),g(x{)}^2,g(x{)}^3...g(x{)}^S$，以及$g(x{)}^S,g(x{)}^{2S},...$
那么$g(x{)}^{i+kS}=g(x{)}^i\ast g(x{)}^{kS}$
所以可以先对每块求出${\sum }_{i=0}^{S-1}{f}_{i+kS}\ast g(x{)}^i$，再乘上$g(x{)}^{kS}$
复杂度$O(nm+\sqrt{n}\cdot mlogm)$。

---

Lucas定理：

$\left(\genfrac{}{}{0ex}{}{n}{m}\right)$中质因子$p$的个数为$\sum _i\lfloor \frac{n}{p^i}\rfloor -\lfloor \frac{m}{p^i}\rfloor -\lfloor \frac{n-m}{p^i}\rfloor =p进制下m+(n-m)的进位次数$

例题：


例题2：

---

$r=1\ast 2...\ast (p-1)\ast (p+1)\ast (p+2)...\ast (2p-1)\ast (2p+1)...\ast (\lfloor \frac{n}{p}\rfloor \ast p-1)\ast 剩下部分$
记$F$为$r$在$\le p^k$部分的乘积，那么$r=F^{\lfloor \frac{n}{p^k}\rfloor }\ast 剩下部分$
$\left(\genfrac{}{}{0ex}{}{n}{m}\right)=r_n\ast (r_m{)}^{-1}\ast (r_{n-m}{)}^{-1}\ast p^{t_n-t_m-t_{n-m}}$

下面这两个比较懵，存疑

---

Stern-Brocot Tree

待更。

---

Min25筛学习笔记

---

例题：

$[(x,n)=d]v=[(\frac{x}{d},\frac{n}{d})=1]v={\sum }_{k|\frac{x}{d},k|\frac{n}{d}}\mu (k)v$
直接修改$a$不好修改，设辅助数组$a_i={\sum }_{d|i}f(d)$，则$f(i)={\sum }_{x|i}a(x)\ast \mu (\frac{i}{x})$。考虑对$f$进行修改，枚举$k$（下面的$d$是定值，$x$是枚举的）：
$$f(i)+=\sum _{x|i}a(x)\mu (\frac{i}{x})=\sum _{kd|x|i}\mu (k)v\mu (\frac{i}{x})=\mu (k)v\sum _{kd|x|i}\mu (\frac{i}{x})$$
令$i^{\prime }=\frac{i}{kd},x^{\prime }=\frac{x}{kd}$：
$$\begin{gathered} =\mu (k)v\sum _{x^{\prime }|i^{\prime }}\mu (\frac{i^{\prime }}{x^{\prime }}) \\ =\mu (k)v[i==kd] \end{gathered}$$
所以只需要枚举$\frac{n}{d}$的约数$k$，令$f[kd]+=\mu (k)v$即可。

${\sum }_{i=1}^{n}a(i)={\sum }_{i=1}^n{\sum }_{d|i}f(d)={\sum }_{i=1}^n{\sum }_{d=1}^{\frac{n}{i}}f(d)$
维护$f$的前缀和然后整除分块回答询问即可。

例题：

求${\sum }_{i=1}^n{\mu }^2(i)，n\le 10^{18}$，即求$n$以内无平方因子数的个数。
${\mu }^2(i)={\sum }_{d^2|i}\mu (d)$
${\sum }_{i=1}^n{\mu }^2(i)={\sum }_{i=1}^n{\sum }_{d^2|i}\mu (d)=\sum _{d=1}^{\sqrt{n}}\mu (d)\lfloor \frac{n}{d^2}\rfloor$
整除分块的复杂度是$O(n^{\frac{1}{3}})$的，当$d\le n^{\frac{1}{3}}$时只有$n^{\frac{1}{3}}$种，当$d>n^{\frac{1}{3}}$时$\lfloor \frac{n}{d^2}\rfloor <n^{\frac{1}{3}}$，而杜教筛筛$\mu (d)$的前缀和的复杂度按照每次都跑满$O(n^{\frac{2}{3}})$算复杂度为$O(n^{\frac{3}{7}})$（实际比这个小）

例题三连：

51nod1222 最小公倍数计数
51Nod1220 约数之和
SPOJ DIVCNT2 - Counting Divisors