[莫比乌斯反演]求sum[gcd(i,j,k)=1]

题目大意

求 $\sum _{i=1}^a\sum _{j=1}^b\sum _{k=1}^c[gcd(i,j,k)=1]$

详细思路

---

先推导这两个式子$\sum _{i=1}^n[gcd(i,m)=1]$,$\sum _{i=1}^n\sum _{j=1}^m[gcd(i,j)=k]$

---

关于$\sum _{i=1}^n[gcd(i,m)=1]$

可以感性思考一下,我们可以通过把$1$到$n$中，把所有关于$m$的质因数的倍数减去,把减多的加回来,再把加多的减回去……

进行了一次容斥操作,可以通过$\mu$实现

$\sum _{i=1}^n[gcd(i,m)=1]=\sum _{d|m}\mu (d)\lfloor \frac{n}{d}\rfloor$

---

关于$\sum _{i=1}^n\sum _{j=1}^m[gcd(i,j)=k]$

设$f(k)=\sum _{i=1}^n\sum _{j=1}^m[gcd(i,j)=k]$

那么它的倍数和$F(k)=\sum _{k|d}f(d)$

把$f(k)=\sum _{i=1}^n\sum _{j=1}^m[gcd(i,j)=k]$代入$F(k)=\sum _{k|d}f(d)$

得到$F(k)=\sum _{k|d}\sum _{i=1}^n\sum _{j=1}^m[gcd(i,j)=d]$

改变一下枚举顺序$\sum _{k|d}\sum _{i=1}^n\sum _{j=1}^m[gcd(i,j)=d]=\sum _{i=1}^n\sum _{j=1}^m\sum _{k|d}[gcd(i,$