数论综合（进阶）

最新推荐文章于 2024-07-25 20:22:48 发布

Znloye

最新推荐文章于 2024-07-25 20:22:48 发布

阅读量107

点赞数

分类专栏：数论文章标签： c++

本文链接：https://blog.csdn.net/Znloye/article/details/120226878

版权

数论专栏收录该内容

2 篇文章 0 订阅

订阅专栏

前言

之前自行整理的数论讲义存在公式不规范、条理不清楚等问题，这里结合一下各个老师各个书本的讲法，试图系统化这部分内容。

一、莫比乌斯反演

1. O(n) 预处理 $\mu$

根据莫比乌斯函数 $\mu$ 的定义，有：

当没有平方项因子时： $\mu(n)=(-1)^{k}$ ，否则 $\mu(n)=0$ 。

2. trick：

$\mathbf 1 \times \varphi=\mathbf{id}$

证明 $\Leftarrow$

3. $O(\sqrt{n})$ 计算单个的 $\varphi(n)$ 值

根据公式

4. O(n) 预处理 $\varphi$ （ $1\sim n$ 中与互质的个数）

根据公式

5. trick

$\sigma_0(nm)=\sum_{i|n}\sum_{j|m}[\gcd(i,j)=1]$

证明 $\Leftarrow$

6. O(n) 预处理 $\sigma_0$ （ $1\sim n$ 中的约数个数）

设是的最小素因子

$\sigma_0(x)=2\sigma_0\left(\dfrac{x}{p}\right)-\sigma_0\left(\dfrac{x}{p^2}\right)$

证明 $\Leftarrow$

7. O(n) 预处理 $\sigma_1$

若

$n=\prod\limits_{i=1}^{k}p_i^{\alpha_i}$

则

$\begin{aligned}\sigma_1(n)&=\prod\limits_{i=1}^k\sum\limits_{j=0}^{\alpha_i}p_i^j\\&=\color{red}\sum_{j=0}^{\alpha_1}p_1^{j}\color{black}\times\sum_{j=0}^{\alpha_2}p_2^{j}\times\cdots\times\sum_{j=0}^{\alpha_k}p_k^{j}\\&=\color{red}(p_1^0+p_1^1+\cdots+p_1^{\alpha_1})\color{black}(p_2^0+p_2^1+\cdots+p_2^{\alpha_2})\cdots(p_k^{0}+p_k^{1}+\cdots+p_k^{\alpha_k})\end{aligned}$

设：

$\begin{aligned}tp(n)&=\color{red}\sum\limits_{j=0}^{\alpha}p_1^j\\&=\color{red}(p_1^0+p_1^1+\cdots+p_1^{\alpha_1})\end{aligned}$

则有递推式：（ p_1 即的最小素因子）

$\begin{aligned}tp(n)&=\begin{cases}tp\left(\dfrac{n}{p_1}\right)\times p_1+1~~~\left(\dfrac{n}{p_1}\mod p_1=0\right)\\p_1+1~~~~~~~~~~~~~~~~~\left(\dfrac{n}{p_1}\mod p_1\ne0\right)\end{cases}\\\sigma_1(n)&=\begin{cases}\dfrac{\sigma_1\left(\dfrac{n}{p_1}\right)}{tp\left(\dfrac{n}{p_1}\right)}\times tp(n)~~~~~~~~~\left(\dfrac{n}{p_1}\mod p_1=0\right)\\\sigma_1\left(\dfrac{n}{p_1}\right)\times(1+p_1)~~~~~\left(\dfrac{n}{p_1}\mod p_1\ne0\right)\end{cases}\end{aligned}$

对于素数有：

Znloye

关注

0
点赞
踩
0

收藏

觉得还不错? 一键收藏
0
评论
数论综合（进阶）

前言之前自行整理的数论讲义存在公式不规范、条理不清楚等问题，这里结合一下各个老师各个书本的讲法，试图系统化这部分内容。一、莫比乌斯反演1. 预处理根据莫比乌斯函数的定义，有：当没有平方项因子时：，否则。for(int j = 1; j < N; ++j) mu[j] = 1;for(int i = 2; i < N; ++i)if(!b[i]){ mu[i] = -1, p[++p[0]] = i; for(int j = i <&...
复制链接

扫一扫