[杜教筛+莫比乌斯反演] 51Nod1238: 最小公倍数之和 V3

题意

求$\sum_{i=1}^n \sum_{j=1}^n lcm(i,j)$
$n \le 10^{10}$

题解

自己推到后面就不会了，然后去找题解……
大概就是

$$\sum_{i=1}^n \sum_{j=1}^n \frac{ij}{gcd(i,j)}=\sum_{d=1}^n \sum_{i=1}^{\lfloor \frac n d \rfloor} \sum_{j=1}^{\lfloor \frac n d \rfloor} [gcd(i,j)=1]\frac {i*j*d^2} d$$
反演一下：
$$=\sum_{d=1}^n\sum_{i=1}^{\lfloor \frac n d \rfloor} \sum_{j=1}^{\lfloor \frac n d \rfloor} i*j*d \sum_{k|gcd(i,j)}\mu(k)=\sum_{k=1}^n \sum_{d=1}^n \sum_{k|i}^{\lfloor \frac n d \rfloor} \sum_{k|j}^{\lfloor \frac n d \rfloor} i*j*d*\mu(k)$$
$$=\sum_{k=1}^n \sum_{d=1}^n \sum_{i=1}^{\lfloor \frac n {dk}\rfloor} \sum_{j=1}^{\lfloor \frac n {dk} \rfloor} ik*jk*d*\mu(k)$$
设 $T=dk$ :
$$=\sum_{T=1}^n \sum_{d|T} (\frac T d)^2*d*\mu(\frac T d) \sum_{i=1}^{\lfloor \frac n T \rfloor} \sum_{j=1}^{\lfloor \frac n T \rfloor} i*j=\sum_{T=1}^n(\frac{(\lfloor \frac n T \rfloor+1)*\lfloor \frac n T \rfloor}{2})^2 T \sum_{d|T} \mu(d)*d$$
现在考虑如何求 $f(T)= T \sum_{d|T} \mu(d)*d$ 的前缀和：
$$\sum_{i=1}^n f(i)=\sum_{i=1}^n i \sum_{d|i} \mu(d)*d=\sum_{d=1}^n \sum_{i=1}^{\lfloor \frac n d \rfloor} \mu(d)*d*i*d=\sum_{d=1}^n \mu(d)*d^2\sum_{i=1}^{{\lfloor \frac n d \rfloor} }i=\sum_{d=1}^n \frac{(1+\lfloor \frac n d \rfloor)*\lfloor \frac n d \rfloor}{2} \mu(d)*d^2$$
在考虑求 $\mu(d)*d^2$ 的前缀和: (我好像用了最复杂的方法......)
杜教筛就好了，卷上 $id^2$,
$$S(n)=1-\sum_{i=2}^n i^2*S(\lfloor \frac n i \rfloor)$$