Ra

Problem

求$\sum_{a=1}^n\sum_{b=1}^n[lcm(a,b)>n]$
$1\le n\le 10^{10}$

Solution

莫比乌斯反演。
首先要求$lcm(a,b)>n$，那么就可以用$n^2$减去$lcm(a,b)\le n$，即：
$\sum_{a=1}^n\sum_{b=1}^n[lcm(a,b)\le n]$
即求满足(a,b)=1且$abc\le n$的三元组(a,b,c)的个数
即$\sum_{(a,b)=1}[\frac{n}{ab}]$
于是我们设：
$f(d)=\sum_{(a,b)=d}[\frac{n}{ab}]$
$g(d)=\sum_{d|(a,b)}[\frac{n}{ab}]=\sum_{a=1}^{[\frac{n}{d}]}\sum_{b=1}^{[\frac{n}{d}]}[\frac{[\frac{n}{d^2}]}{ab}]$
发现g(d)的式子中a和b的取值不可能超过$[\frac{n}{d^2}]$于是：
$g(d)=\sum_{d|(a,b)}[\frac{n}{ab}]=\sum_{a=1}^{[\frac{n}{d^2}]}\sum_{b=1}^{[\frac{n}{d^2}]}[\frac{[\frac{n}{d^2}]}{ab}]$
由g(d)和f(d)的定义可得：
$g(d)=\sum_{i=1}^{\frac{n}{d}}f(i\times d)$
反演得到：
$f(d)=\sum_{i=1}^{\frac{n}{d}}\mu(i)\times g(i\times d)$
于是：
$Ans=f(1)=\sum_{i=1}^{\frac{n}{d}}\mu(i)\times g(i)$
$=\sum_{i=1}^{\frac{n}{d}}\mu(i)\times[abc\le n]$
按照张俊的方法，可以假定$a\le b\le c$，然后特殊情况另外讨论，发现a,b的范围是有个区间的，然后c是可以根据a,b算出来的，于是就可以做了。