众所周知,\(d\)和\(\sigma\)都是积性函数,今天我们来对这俩动刀子。
我打算只讲稍难\(\sigma\),剩下的那个\(d\),留给大家自己思考。
根据定义得到:
\(\begin{aligned}\sigma(p^k)=1+p+p^2+...+p^k\end{aligned}\)
然后我们就能推导出\(\sigma\)的求值式:
\(\begin{aligned}n=\prod_{i=1}p_i^{c_i}\end{aligned},\)
\(\begin{aligned}\sigma(n)=\prod_{i=1}1+p_i+p_i^2...+p_i^{c_i}\end{aligned}\)
尝试用线性筛构造这个函数:
设\(x=i\times p_1\),其中\(p\)为\(x\)的最小质因子,\(\begin{aligned}x=\prod_{i=1}p_i^{c_i}\end{aligned}\)。
\(x\)为质数:\(\sigma(x)=1+x\)
当\(i\)与\(p_1\)互质时:\(\sigma(x)=\sigma(i)\times\sigma(p_1)\)
当\(i\%p_1=0\)时:
我们具体来讨论一下这个情况。
\(\begin{aligned}\sigma(x)=\prod_{i=1}1+p_i+p_i^2+...+p_i^{c_i}\end{aligned}\)
\(\begin{aligned}\sigma(i)=(1+p_1+p_2^2+..+p_1^{c_1-1})\prod_{i=2}1+p_i+p_i^2+...+p_i^{c_i-1}\end{aligned}\)
我们想要直接通过\(\sigma(i)\)来求\(\sigma(x)\)显然是不现实的,因为我们不知道\(c_1\)是多少。
那么我们应该怎么做呢?
我们尝试把这个包含绝对关系的式子转换成包含相对关系的式子。
因为我们\(x\)和\(i\),从第二项开始就是一样的了,所以我们可以设:\(\begin{aligned}T=\prod_{i=2}1+p_i+p_i^2+p_i^{c_i}\end{aligned}\)
然后式子就转换为了:
\(\begin{aligned}\sigma(x)=\sigma(i)+p_1^{c_1}\times T\end{aligned}\)
我们来考虑一下怎么把这个讨厌的\(p_1^{c_1}\)给处理了。
注意我们一开始的目标:把绝对关系转化成相对关系。
既然我们不知道\(p_1^{c_1}\times T\)的值,那么我们能不能求出\(p_1^{c_1-1}\times T\)的值呢?
注意到\(\begin{aligned}\sigma(i)=\sigma(\frac i {p_1})+p_1^{c_1-1}\times T\end{aligned}\)。
然后怎么做应该就不用我多说了吧——盘他!
\(\begin{aligned}\sigma(x)&=\sigma(i)+p\times(\sigma(i)-\sigma(\frac i {p_1}))\\&=(1+p)\times\sigma(i)-p\times\sigma(\frac i{p_1})\end{aligned}\)
然后我们就能直接用线性筛来求约数函数了。
最后,我在稍微写一下\(d\)的式子吧:
\(\begin{aligned}d(x)=2\times d(i)-d(\frac i {p_1})\end{aligned}\)