51Nod1238最小公倍数之和V3的另一种做法

前面的讲解【Blog地址】

---

题目意思：求
$$\sum _{i=1}^n\sum _{j=1}^{m}lcm(i,j)$$
对于数据范围$n,m\le 10^{10}$。

---

其实我们前面将式子化简成了：

$$=\sum _{d=1}^{n}d\sum _{w=1}^{\lfloor \frac{n}{d}\rfloor }\mu (w)w^{2}S(\lfloor \frac{n}{dw}\rfloor )S(\lfloor \frac{n}{dw}\rfloor )$$

其中$S(n)=\frac{n\times (n+1)}{2}$

其实这个时候我们用杜教筛筛${\sum }_{i=1}^n\mu (i)i^2$即可在$O(n^{\frac{3}{4}})$时间内做出（分块套分块再加杜教筛）。

对于$f(x)=\mu (x)x^2$的前缀和$F(n)$，我们可以根据杜教筛的套路，令$g(x)=x^2$，那么$h=f⨂g$就等于（$⨂$为狄利克雷卷积）：

$$\begin{gathered} =\sum _{d|n}\mu (d)d^2{\left(\frac{n}{d}\right)}^2 \\ =n^2\sum _{d|n}\mu (d) \\ =n^2[n=1] \\ =[n=1] \end{gathered}$$

此时带入杜教筛公式即可得，其中后面要除以$g(1)$，但因为$g(1)=1$这里略去：

$$\begin{gathered} F(n)=\sum _{i=1}^{n}h(i)-\sum _{i=2}^{n}g(i)F(\lfloor \frac{n}{i}\rfloor ) \\ =1-\sum _{i=2}^n{i}^{2}F(\lfloor \frac{n}{i}\rfloor ) \end{gathered}$$

就可以直接递归筛了。

---

但是我们可以继续看，转而枚举乘积$dw$，令$T=dw$则可以得到原式等于：

$$\begin{gathered} =\sum _{T=1}^{n}S(\lfloor \frac{n}{T}\rfloor )S(\lfloor \frac{m}{T}\rfloor )\sum _{d|T}^n\mu (d)d^2\left(\frac{T}{d}\right) \\ =\sum _{T=1}^{n}S(\lfloor \frac{n}{T}\rfloor )S(\lfloor \frac{m}{T}\rfloor )T\sum _{d|T}\mu (d)d \end{gathered}$$

那么我们只需要杜教筛出后面$T{\sum }_{d|T}\mu (d)d$即可分块加杜教筛，在$O(n^{\frac{2}{3}})$的时间内做出。

但是如何筛后面的那个呢？

我们来看，现在相当于求这个式子：

$$\sum _{d=1}^{n}d\sum _{w|d}\mu (w)w$$

由于前面的式子和现在这个求出来的值是一样的，也就是：

$$\sum _{d=1}^{n}d\sum _{w=1}^{\lfloor \frac{n}{d}\rfloor }\mu (w)w^{2}S(\lfloor \frac{n}{dw}\rfloor )S(\lfloor \frac{n}{dw}\rfloor )=\sum _{T=1}^{n}S(\lfloor \frac{n}{T}\rfloor )S(\lfloor \frac{m}{T}\rfloor )T\sum _{d|T}\mu (d)d$$

所以我们可以将$S$相同的部分减去，那么剩下就是：

$$\sum _{d=1}^{n}d\sum _{w=1}^{\lfloor \frac{n}{d}\rfloor }\mu (w)w^2=\sum _{T=1}^{n}T\sum _{d|T}\mu (d)d$$

所以相当于筛前面的那个东西，我们预处理线性筛出一部分${\sum }_{T=1}^{n}T{\sum }_{d|T}\mu (d)d$的这个，然后用杜教筛加分块继续求前面的那个即可。

另外一种解释：
相当于我们后面需要杜教筛的$d{\sum }_{w|d}\mu (w)w={\sum }_{w|d}\mu (w)w^2\frac{d}{w}$
写成狄利克雷卷积的形式就是$(i{d}^2\cdot \mu )⨂id$，我们令$f=(i{d}^2\cdot \mu )⨂id,g=i{d}^2$，那么$f⨂g=(i{d}^2\cdot \mu )⨂id⨂i{d}^2$，交换一下就是$(i{d}^2\cdot \mu )⨂i{d}^2⨂id$，前两个卷出来就是$ϵ$，然后原式就变成$ϵ⨂id=id$，然后杜教筛即可。
公式即为：
$$F(n)=\sum _{i=1}^{n}i-\sum _{i=2}^n{i}^{2}F(\lfloor \frac{n}{i}\rfloor )$$

我觉得上面这种更好

---

代码，咕咕咕

---

参考资料博客：

• Orz-Dalao%%%
• Orz-Dalao%%%