51nod 1227 平均最小公倍数 莫比乌斯反演+杜教筛

题意

Lcm(a,b)表示a和b的最小公倍数，A(n)表示Lcm(n,i)的平均数(1 <= i <= n)，
例如：A(4) = (Lcm(1,4) + Lcm(2,4) + Lcm(3,4) + Lcm(4,4)) / 4 = (4 + 4 + 12 + 4) / 4 = 6。
F(a, b) = A(a) + A(a + 1) + …… A(b)。(F(a,b) = ∑A(k), a <= k <= b)
例如：F(2, 4) = A(2) + A(3) + A(4) = 2 + 4 + 6 = 12。
给出a,b，计算F(a, b)，由于结果可能很大，输出F(a, b) % 1000000007的结果即可。
1 <= a <= b <= 10^9

分析

我们要求的实际就是

$$\sum_{i=1}^n\sum_{j=1}^i\frac{j}{gcd(i,j)}$$
$$=\sum_{d=1}^n\sum_{i=1}^{\lfloor\frac{n}{d}\rfloor}\sum_{j=1}^ij*[gcd(i,j)=1]$$
到这里为止就有两种方法，一种是转化为用 φ(d)