杜教筛

最新推荐文章于 2023-09-29 19:51:34 发布

我是一只计算鸡

最新推荐文章于 2023-09-29 19:51:34 发布

阅读量289

点赞数

分类专栏：数论与数论函数文章标签：杜教筛

版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。

本文链接：https://blog.csdn.net/giftedpanda/article/details/100049469

版权

数论与数论函数专栏收录该内容

16 篇文章 0 订阅

订阅专栏

杜教筛用于求数论函数前缀和。

比如常见的莫比乌斯函数、欧拉函数。

当 n = 10^4 , 我会 for 循环

当n = 10^6 ，我会线性筛

当 n = 10^9 ，我。。。

杜教筛可以低于线性时间复杂度求解数论函数前缀和， $\Theta(n^{\frac{2}{3}})$ 。

对于数论函数，我们要求 $S(n) = \sum_{i = 1}^nf(i)$

对于任意一个数论函数，

$\sum_{i = 1}^n(f*g)(i)$

$=\sum_{i = 1}^n\sum_{d|i}f(d)g(\frac{i}{d})$

$=\sum_{i=1}^n\sum_{j=1}^{\left \lfloor \frac{n}{i} \right \rfloor}g(i)f(j)$

$=\sum_{i=1}^ng(i)\sum_{j=1}^{\left \lfloor \frac{n}{i} \right \rfloor}f(i)$

$=\sum_{i=1}^ng(i)S(\left \lfloor \frac{n}{i} \right \rfloor)$

$=g(1)S(n) + \sum_{i=2}^ng(i)S(\left \lfloor \frac{n}{i} \right \rfloor)$

$\Rightarrow g(1)S(n) = \sum_{i=1}^n(f*g)(i)-\sum_{i=2}^ng(i)S(\left \lfloor \frac{n}{i} \right \rfloor)$

如果我们可以很快的算出 $\sum_{i=1}^n(f*g)(i)$ ，然后用数论分块算出 $\sum_{i=2}^n g(i)S(\left \lfloor \frac{n}{i} \right \rfloor)$ ，我们就能很快得到 g(1)S(n)

莫比乌斯函数

$S(n) = \sum_{i = 1} ^n \mu(i)$

$f(i) = \mu(i)$

$\sum_{i = 1} ^n(1*\mu)(i) = \sum_{i = 1} ^nS(\left \lfloor \frac{n}{i} \right \rfloor)$ , 由于 $1 * \mu = \varepsilon$

$1 = S(n) + \sum_{i = 2} ^ nS(\left \lfloor \frac{n}{i} \right \rfloor)$

$S(n) = 1 - \sum_{i = 2} ^nS(\left \lfloor \frac{n}{i} \right \rfloor)$

欧拉函数

$\sum_{i = 1} ^ n \varphi(i)$

$f(i) = \varphi(i)$

$\sum_{i=1}^n(1*\varphi)(i) = \sum_{i=1}^nS(\left \lfloor \frac{n}{i}\right \rfloor)$ ，由于 $1*\varphi = ID$

$\frac{1}{2}n(n+1) =S(n)+ \sum_{i=2}^nS(\left \lfloor \frac{n}{i} \right \rfloor)$

$S(n) = \frac{1}{2}n(n+1) - \sum_{i=2}^nS(\left \lfloor \frac{n}{i} \right \rfloor)$

当 n = 10^4 , 我会 for 循环

当 n = 10^6 ，我会线性筛

当 n = 10^9 ，我会杜教筛

我是一只计算鸡

关注

0
点赞
踩
2

收藏

觉得还不错? 一键收藏
0
评论
杜教筛

杜教筛用于求数论函数前缀和。比如常见的莫比乌斯函数、欧拉函数。当, 我会循环当n =，我会线性筛当，我。。。杜教筛可以低于线性时间复杂度求解数论函数前缀和，。对于数论函数，我们要求对于任意一个数论函数，如果我们可以很快的算出，然后用数论分块算出，我们就能很快得到莫比乌斯函数, 由于欧拉函数，由...
复制链接

扫一扫

专栏目录

评论

被折叠的条评论为什么被折叠?

到【灌水乐园】发言

查看更多评论

添加红包

成就一亿技术人!

hope_wisdom

发出的红包

实付元

使用余额支付

点击重新获取

扫码支付

钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。