2021SCUACM FallTraining 05 数学 J——中规中矩的莫比乌斯反演
Posted hans774882968
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了2021SCUACM FallTraining 05 数学 J——中规中矩的莫比乌斯反演相关的知识,希望对你有一定的参考价值。
看到一大串式子里面出现了gcd,可以猜测是要设gcd为d,并搞出gcd(...) == 1
的形式,然后莫比乌斯反演了。题目中f
函数的定义实际上就是莫比乌斯函数的平方。
a
n
s
=
∑
d
=
1
n
μ
2
[
d
]
∗
d
∗
(
∑
a
[
1
]
=
1
n
.
.
.
∑
a
[
x
]
=
1
n
∏
j
=
1
x
(
a
[
j
]
)
k
[
d
=
=
g
c
d
(
a
[
1
]
,
.
.
.
,
a
[
x
]
)
]
)
ans = \\sum_{d=1}^n \\mu^2[d]*d * (\\sum_{a[1]=1}^n...\\sum_{a[x]=1}^n \\prod_{j=1}^x (a[j])^k [d == gcd(a[1],...,a[x])])
ans=d=1∑nμ2[d]∗d∗(a[1]=1∑n...a[x]=1∑nj=1∏x(a[j])k[d==gcd(a[1],...,a[x])])
a[j]
都是d的倍数,因此改成枚举倍数系数,得:
a
n
s
=
∑
d
=
1
n
μ
2
[
d
]
∗
d
x
∗
k
+
1
∗
(
∑
a
[
1
]
=
1
n
/
d
.
.
.
∑
a
[
x
]
=
1
n
/
d
∏
j
=
1
x
(
a
[
j
]
)
k
[
1
=
=
g
c
d
(
a
[
1
]
,
.
.
.
,
a
[
x
]
)
]
)
ans = \\sum_{d=1}^n \\mu^2[d]*d^{x*k+1} * (\\sum_{a[1]=1}^{n/d}...\\sum_{a[x]=1}^{n/d} \\prod_{j=1}^x (a[j])^k [1 == gcd(a[1],...,a[x])])
ans=d=1∑nμ2[d]∗dx∗k+1∗(a[1]=1∑n/d...a[x]=1∑n/dj=1∏x(a[j])k[1==gcd(a[1],...,a[x])])
[1 == gcd(a[1],...,a[x])]
可以写成
∑
e
=
1
n
/
d
μ
[
e
]
∗
[
e
∣
g
c
d
(
a
[
1
]
,
.
.
.
,
a
[
x
]
)
]
\\sum_{e=1}^{n/d}\\mu[e]*[e \\ | \\ gcd(a[1],...,a[x])]
e=1∑n/dμ[e]∗[e ∣ gcd(a[1],...,a[x])]
把不受束缚的项拿到外面,再改成枚举倍数系数,然后进行简单化简得:
a
n
s
=
∑
d
=
1
n
μ
2
[
d
]
∗
d
x
∗
k
+
1
∗
∑
e
=
1
n
/
d
μ
[
e
]
∗
e
x
∗
k
∑
a
[
1
]
=
1
n
/
(
d
∗
e
)
.
.
.
∑
a
[
x
]
=
1
n
/
(
d
∗
e
)
∏
j
=
1
x
(
a
[
j
]
)
k
ans = \\sum_{d=1}^n \\mu^2[d]*d^{x*k+1} * \\sum_{e=1}^{n/d} \\mu[e] * e^{x*k} \\sum_{a[1]=1}^{n/(d*e)}...\\sum_{a[x]=1}^{n/(d*e)} \\prod_{j=1}^x (a[j])^k
ans=d=1∑nμ2[d]∗dx∗k+1∗e=1∑n/dμ[e]∗ex∗ka[1]=1∑n/(d∗e)...a[x]=1∑n/(d∗e)j=1∏x(a[j])k
a[t1]
和a[t2]
终于没有绑定关系了,于是右边的式子可以写为
(
∑
i
=
1
n
/
T
i
k
)
x
,
T
=
d
∗
e
(\\sum_{i=1}^{n/T}i^k)^x,T=d*e
(i=1∑n/Tik)x,T=d∗e
按照莫比乌斯反演的套路,此时我们应该引入T = d*e
的限制来进一步处理了。
a
n
s
=
∑
T
=
1
n
∑
d
=
1
n
μ
2
[
d
]
∗
d
x
∗
k
+
1
∗
∑
e
=
1
n
/
d
[
T
=
d
∗
e
]
∗
(
μ
[
e
]
∗
e
x
∗
k
∗
(
∑
i
=
1
n
/
T
i
k
)
x
)
ans = \\sum_{T=1}^n \\sum_{d=1}^n \\mu^2[d]*d^{x*k+1} * \\sum_{e=1}^{n/d} [T=d*e] * (\\mu[e] * e^{x*k} * (\\sum_{i=1}^{n/T}i^k)^x)
ans=T=1∑nd=1∑nμ2[d]∗dx∗k+1∗e=1∑n/d[T=d∗e]∗(μ[e]∗ex∗k∗(i=1∑n/Tik)x)
我们把枚举了d
和e
的式子看成一个整体。因为外层是T
,所以必须保证d
是T
的因子。并且因为枚举e
的时候T
和d
都是固定的,所以只有1个e
满足条件:T/d
。于是e
那层枚举可以消失。
a
n
s
=
∑
T
=
1
n
(
∑
i
=
1
n
/
T
i
k
)
x
∑
d
∣
T
μ
2
[
d
]
∗
d
x
∗
k
+
1
∗
μ
[
T
/
d
]
∗
(
T
/
d
)
x
∗
k
ans = \\sum_{T=1}^n (\\sum_{i=1}^{n/T}i^k)^x \\sum_{d|T} \\mu^2[d]*d^{x*k+1}*\\mu[T/d]*(T/d)^{x*k}
ans=T=1∑n(i=1∑n/Tik)xd∣T∑以上是关于2021SCUACM FallTraining 05 数学 J——中规中矩的莫比乌斯反演的主要内容,如果未能解决你的问题,请参考以下文章
四川大学2021SCUACM新生赛决赛大部分题解(差分dp线段树……)