显然枚举
a
i
,
a
j
a_i,a_j
ai,aj 很不好看,我们改成枚举数值,直接莫比乌斯反演
A
n
s
=
∑
d
n
d
φ
(
d
)
∑
i
n
/
d
∑
j
n
/
d
i
j
[
(
i
,
j
)
=
1
]
c
n
t
i
d
c
n
t
j
d
=
∑
d
d
φ
(
d
)
∑
l
n
/
d
μ
(
l
)
(
∑
l
∣
i
i
c
n
t
i
d
)
(
∑
l
∣
j
j
c
n
t
j
d
)
=
∑
d
d
φ
(
d
)
∑
l
μ
(
l
)
l
2
(
∑
i
n
/
d
l
i
c
n
t
i
d
l
)
2
=
∑
T
(
∑
i
n
/
T
i
c
n
t
i
T
)
2
(
∑
d
∣
T
d
φ
(
d
)
μ
(
T
d
)
(
T
d
)
2
)
Ans=\sum_d^n d\varphi(d)\sum_i^{n/d}\sum_j^{n/d}ij[(i,j)=1]cnt_{id}cnt_{jd}\\=\sum_d d\varphi(d)\sum_{l}^{n/d}\mu(l)(\sum_{l|i}icnt_{id})(\sum_{l|j}jcnt_{jd})\\=\sum_dd\varphi(d)\sum_l\mu(l)l^2(\sum_{i}^{n/dl}icnt_{idl})^2\\=\sum_T(\sum_i^{n/T}icnt_{iT})^2(\sum_{d|T}d\varphi(d)\mu(\frac{T}{d})(\frac{T}{d})^2)
Ans=d∑ndφ(d)i∑n/dj∑n/dij[(i,j)=1]cntidcntjd=d∑dφ(d)l∑n/dμ(l)(l∣i∑icntid)(l∣j∑jcntjd)=d∑dφ(d)l∑μ(l)l2(i∑n/dlicntidl)2=T∑(i∑n/TicntiT)2(d∣T∑dφ(d)μ(dT)(dT)2) 后面提一个
T
T
T 出来,是
(
μ
⋅
I
d
)
∗
φ
(\mu \cdot Id)*\varphi
(μ⋅Id)∗φ,而
(
μ
⋅
I
d
)
∗
φ
∗
I
=
(
μ
⋅
I
d
)
∗
I
d
=
e
(\mu \cdot Id)*\varphi*I=(\mu \cdot Id)*Id=e
(μ⋅Id)∗φ∗I=(μ⋅Id)∗Id=e,故
(
μ
⋅
I
d
)
∗
φ
=
μ
(\mu \cdot Id)*\varphi=\mu
(μ⋅Id)∗φ=μ,故后面一坨就是
T
μ
(
T
)
T\mu(T)
Tμ(T),莫队即可