ICode9

精准搜索请尝试: 精确搜索
首页 > 其他分享> 文章详细

杜教筛

2022-01-30 14:33:44  阅读:174  来源: 互联网

标签:phi 积性 sum times mu 杜教 id


认识一下常用数论函数$:$

$\mu$ 无需解释

$I =1\ $恒等函数(恒等于$1$)

$e= \ [n=1]$原函数

然后考虑莫比乌斯反演

性质式$:$

$\sum_{i|n}\mu(i)=[n=1]$

等价于

$\mu*I=e$

$F=\sum_{i|n}f[i]$

$F=(f*I)$

$F*\mu=(f*I)*\mu$

$F*\mu=f*(I*\mu)$

$F*\mu=f*e$

$F*\mu=f$

$f=F*\mu$

证毕

$n=\sum_{i|n}\phi(i)$

继续推导

$id=\phi*I$

$id*\mu=\phi*I*\mu$

$id*\mu=\phi*e$

$id*\mu=\phi$

$\phi(n)=\sum_{i|n}n/i\times\mu(i)$

$\phi(n)/n=\sum_{i|n}\mu(i)/i$

这貌似就是$\mu$和$\phi$的互相转化

那么杜教筛是什么$:$

考虑求$\sum f(i)$

$f(i)$为积性函数

积性函数$*$积性函数$=$积性函数

设$h=f*g$

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

我们先求一个$h$的前缀和

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

$\sum_{i=1}^{n}h(i)=\sum_{d=1}^{n}g(d)\sum_{i=1}^{n/d}f(i)$

$\sum_{i=1}^{n}h(i)=\sum_{d=1}^{n}g(d)S(n/d)$

$\sum_{i=1}^{n}h(i)=g(1)S(n)+\sum_{d=2}^{n}g(d)S(n/d)$

$g(1)S(n)=\sum_{i=1}^{n}h(i)-\sum_{d=2}^{n}g(d)S(n/d)$

$g(1)S(n)=\sum_{i=1}^{n}(f*g)(i)-\sum_{d=2}^{n}g(d)S(n/d)$

当$h$的前缀和很好求,可以在$O(n^{2/3})$得到$S(n)$

感觉不会的话都尝试卷一下就可以了.

例$1.$

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

如果是我,我会猜$g=I,h=e$

$S(n)=\sum_{i=1}^{n}e(i)-\sum_{d=2}^{n}S(n/d)$

$S(n)=1-\sum_{d=2}^{n}S(n/d)$

后面的显然可以整除分块,然后我发现这个可以递归?不用整除分块?

例$2.$

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

$g=I,h=id$

$S(n)=\sum_{i=1}^{n}id(i)-\sum_{d=2}^{n}I(d)S(n/d)$

$S(n)=\sum_{i=1}^{n}i-\sum_{d=2}^{n}S(n/d)$

这也很显然了

例$3.$

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

考虑如果卷一下

$f*g=\sum_{i|n}(i\times\phi(i))*g(n/i)$

显然把后面配成$id$会很好,我一开始也这么想的...

$f*g=\sum_{i|n}n\times\phi(i)$

$f*g=n\times\sum_{i|n}\phi(i)$

$f*g=n^{2}$

$S(n)=\sum_{i=1}^{n}i^{2}-\sum_{d=2}^{n}id(d)\times S(n/d)$

我只能说很强了.

标签:phi,积性,sum,times,mu,杜教,id
来源: https://www.cnblogs.com/Force-A/p/15856799.html

本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有