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=1nμ2[d]d(a[1]=1n...a[x]=1nj=1x(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=1nμ2[d]dxk+1(a[1]=1n/d...a[x]=1n/dj=1x(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=1n/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=1nμ2[d]dxk+1e=1n/dμ[e]exka[1]=1n/(de)...a[x]=1n/(de)j=1x(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=1n/Tik)x,T=de
按照莫比乌斯反演的套路,此时我们应该引入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=1nd=1nμ2[d]dxk+1e=1n/d[T=de](μ[e]exk(i=1n/Tik)x)
我们把枚举了de的式子看成一个整体。因为外层是T,所以必须保证dT的因子。并且因为枚举e的时候Td都是固定的,所以只有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=1n(i=1n/Tik)xdT以上是关于2021SCUACM FallTraining 05 数学 J——中规中矩的莫比乌斯反演的主要内容,如果未能解决你的问题,请参考以下文章

2021SCUACM集训队冬季选拔2大部分题解

四川大学2021SCUACM新生赛决赛大部分题解(差分dp线段树……)

SCUACM22暑假集训前劝退赛部分题解

(补题 水题 汇总)四川大学第二届SCUACM新生赛

Spring Cloud 2021.0.0-M1(代号 Jubilee)发布

Spring Cloud 2021.0.0-M1(代号 Jubilee)发布