2021年GDCPC广东省大学生程序设计竞赛（热身赛） A.Calculation（贝尔级数）

2021年GDCPC广东省大学生程序设计竞赛（热身赛） A.Calculation（贝尔级数）

传送门：https://ac.nowcoder.com/acm/contest/17796/A?&headNav=acm
题目大意：给定$n$，求$\sum _{i=1}^n\sum _{j=1}^n\sum _{k=1}^n\lambda (i)\lambda (j)[i\mid jandj\mid k]$，$\lambda (n)$为刘维尔函数，$\lambda (n)=(-1{)}^c$，$c$为$n$中所包含的质因数个数且重复计数。
题解：首先考虑改变枚举顺序消去$[i\mid jandj\mid k]$：
$$\begin{array}{rl}& \Rightarrow \sum _{i=1}^n\sum _{j=1}^n\sum _{k=1}^{\lfloor \frac{n}{j}\rfloor }\lambda (i)\lambda (j)[i\mid j]\\ & =\sum _{i=1}^n\sum _{j=1}^{\lfloor \frac{n}{i}\rfloor }\sum _{k=1}^{\lfloor \frac{n}{ij}\rfloor }\lambda (i)\lambda (ij)\\ & =\sum _{i=1}^n\sum _{j=1}^{\lfloor \frac{n}{i}\rfloor }\lfloor \frac{n}{ij}\rfloor \lambda (i)\lambda (ij)\end{array}$$
考虑刘维尔函数的定义，显然这个函数为完全积性函数且对于一个平方数$x^2$，一定有$\lambda (x^2)=1$，那么：
$$\begin{array}{rl}& \Rightarrow \sum _{i=1}^n\sum _{j=1}^{\lfloor \frac{n}{i}\rfloor }\lfloor \frac{n}{ij}\rfloor \lambda (i^{2}j)\\ & =\sum _{i=1}^n\sum _{j=1}^{\lfloor \frac{n}{i}\rfloor }\lfloor \frac{n}{ij}\rfloor \lambda (j)\\ & =\sum _{T=1}^n\lfloor \frac{n}{T}\rfloor \sum _{d\mid T}\lambda (d)\end{array}$$
第三步是枚举$T=ij$时得出来的转变。
现在只需要求解后一部分的$\sum _{d\mid T}\lambda (d)$。
考虑$\lambda (n)$的贝尔级数，不懂贝尔级数的读者可以到这篇文章进行学习。
$$\begin{gathered} {\lambda }_p(x)=1-x+x^2-x^3+\cdots \\ x{\lambda }_p(x)=x-x^2+x^3-x^4+\cdots \\ (1+x){\lambda }_p(x)=1 \\ {\lambda }_p(x)=\frac{1}{1+x} \end{gathered}$$
这里的$\sum _{d\mid T}\lambda (d)$相当于：
$$(\lambda \ast 1)(x)\Rightarrow {\lambda }_p(x)1_p(x)=\frac{1}{1+x}\times \frac{1}{1-x}=\frac{1}{1-x^2}$$
逆推一下这个式子可以得到：
$$\begin{gathered} {f}_p(x)=\frac{1}{1-x^2} \\ {f}_p(x)=1+x^2{f}_p(x) \\ {f}_p(x)=1+0+x^2+0+x^4+\cdots \end{gathered}$$
所以函数$f(x)$的意义是什么？由于两个积性函数的狄利克雷卷积仍是积性函数，那么当且仅当一个数的所有质因数的幂次可以被$2$整除时，$f(x)=1$，否则$f(x)=0$。说明$f(x)=[xisasquarenumber]$。
此时原式等于：
$$\begin{array}{rl}& \Rightarrow \sum _{T=1}^n\lfloor \frac{n}{T}\rfloor f(T)\\ & =\sum _{T=1}^n\lfloor \frac{n}{T}\rfloor [Tisasquarenumber]\\ & =\sum _{T=1}^n\lfloor \frac{n}{T^2}\rfloor \end{array}$$
时间复杂度$O(\sqrt{n})$ 。