P3455 [POI2007]ZAP-Queries
2020-10-19
题意就是
∑
i
=
1
a
∑
j
=
1
b
[
g
c
d
(
i
,
j
)
=
x
]
\sum_{i=1}^{a}{\sum_{j=1}^{b}{[gcd(i,j)=x]}}
i=1∑aj=1∑b[gcd(i,j)=x]
就是这个看这个简单的式子,看着简单,他不一定简单。看题分析一下吧。
我们先把这个式子化一下,化成
g
c
d
(
i
,
j
)
=
1
gcd(i,j)=1
gcd(i,j)=1因为在数学中我们比较喜欢质数
∑
i
=
1
a
x
∑
j
=
1
b
x
[
g
c
d
(
i
,
j
)
=
1
]
\sum_{i=1}^{\frac{a}{x}}{\sum_{j=1}^{\frac{b}{x}}{[gcd(i,j)=1]}}
i=1∑xaj=1∑xb[gcd(i,j)=1]
简单点之后的
a
x
\frac{a}{x}
xa我们用
l
l
l代替,
b
x
\frac{b}{x}
xb我们用
r
r
r代替,这个式子就可以化成:
∑
i
=
1
l
∑
j
=
1
r
[
g
c
d
(
i
,
j
)
=
1
]
\sum_{i=1}^{l}{\sum_{j=1}^{r}{[gcd(i,j)=1]}}
i=1∑lj=1∑r[gcd(i,j)=1]
之后我们根据莫比乌斯函数
[
x
=
1
]
=
∑
d
∣
x
μ
(
d
)
[x=1] = \sum_{d|x}{\mu(d)}
[x=1]=∑d∣xμ(d) 我们就可以得出
∑
i
=
1
a
x
∑
j
=
1
b
x
∑
a
∣
g
c
d
(
i
,
j
)
μ
(
d
)
\sum_{i=1}^{\frac{a}{x}}{\sum_{j=1}^{\frac{b}{x}}{\sum_{a|gcd(i,j)}{\mu(d)}}}
i=1∑xaj=1∑xba∣gcd(i,j)∑μ(d)