LibreOJ #6207. 米缇(杜教筛 + 拉格朗日插值)
#6207. 米緹
推式子
∑i=1n∑j=1ndK(ij)∑i=1n∑j=1n∑x∣i∑y∣j[gcd(x,y)=1]ixkyk∑i=1n∑j=1n∑x∣i∑y∣j∑d∣gcd(x,y)μ(d)ixkyk∑d=1nμ(d)dk∑i=1nd∑x∣iixk∑j=1nd∑y∣iyk∑d=1nμ(d)dk∑i=1nd∑x∣ixk∑j=1nd∑y∣iyk∑d=1nμ(d)dk(∑i=1nd∑x∣ixk)2∑d=1nμ(d)dk(∑x=1ndxk∑x∣i)2∑d=1nμ(d)dk(∑x=1ndxk∑x=1nxd)2∑d=1nμ(d)dk(∑i=1ndiknid)2\sum_{i = 1} ^{n} \sum_{j = 1} ^{n} d_K(ij)\\ \sum_{i = 1} ^{n} \sum_{j = 1} ^{n} \sum_{x \mid i} \sum_{y \mid j} [gcd(x, y) = 1]\frac{i}{x} ^ ky ^ k\\ \sum_{i = 1} ^{n} \sum_{j = 1} ^{n} \sum_{x \mid i} \sum_{y \mid j} \sum_{d \mid gcd(x, y)} \mu(d) \frac{i}{x} ^ ky ^ k\\ \sum_{d = 1} ^{n} \mu(d) d ^ k \sum_{i = 1} ^{\frac{n}ozvdkddzhkzd} \sum_{x \mid i} \frac{i}{x} ^ k \sum_{j = 1} ^{\frac{n}ozvdkddzhkzd} \sum_{y \mid i} y ^ k\\ \sum_{d = 1} ^{n} \mu(d) d ^ k \sum_{i = 1} ^{\frac{n}ozvdkddzhkzd} \sum_{x \mid i} x ^ k \sum_{j = 1} ^{\frac{n}ozvdkddzhkzd} \sum_{y \mid i} y ^ k\\ \sum_{d = 1} ^{n} \mu(d) d ^ k (\sum_{i = 1} ^{\frac{n}ozvdkddzhkzd} \sum_{x \mid i} x ^ k) ^ 2\\ \sum_{d = 1} ^{n} \mu(d) d ^ k (\sum_{x = 1} ^{\frac{n}ozvdkddzhkzd}x ^ k \sum_{x \mid i}) ^ 2\\ \sum_{d = 1} ^{n} \mu(d) d ^ k (\sum_{x = 1} ^{\frac{n}ozvdkddzhkzd}x ^ k \sum_{x = 1} ^{\frac{n}{xd}}) ^ 2\\ \sum_{d = 1} ^{n} \mu(d) d ^ k (\sum_{i = 1} ^{\frac{n}ozvdkddzhkzd}i ^ k \frac{n}{id}) ^ 2\\ i=1∑n?j=1∑n?dK?(ij)i=1∑n?j=1∑n?x∣i∑?y∣j∑?[gcd(x,y)=1]xi?kyki=1∑n?j=1∑n?x∣i∑?y∣j∑?d∣gcd(x,y)∑?μ(d)xi?kykd=1∑n?μ(d)dki=1∑dn??x∣i∑?xi?kj=1∑dn??y∣i∑?ykd=1∑n?μ(d)dki=1∑dn??x∣i∑?xkj=1∑dn??y∣i∑?ykd=1∑n?μ(d)dk(i=1∑dn??x∣i∑?xk)2d=1∑n?μ(d)dk(x=1∑dn??xkx∣i∑?)2d=1∑n?μ(d)dk(x=1∑dn??xkx=1∑xdn??)2d=1∑n?μ(d)dk(i=1∑dn??ikidn?)2
考慮求∑i=1nμ(i)ik\sum\limits_{i = 1} ^{n} \mu(i) i ^ ki=1∑n?μ(i)ik我們卷上IidkIid ^kIidk得到S(n)=1?∑i=2nikS(ni)S(n) = 1 - \sum\limits_{i = 2} ^{n} i ^ k S(\frac{n}{i})S(n)=1?i=2∑n?ikS(in?)這一項(xiàng)可以通過拉格朗日插值跟杜教篩得到。
考慮求∑i=1nikni\sum\limits_{i = 1} ^{n} i ^ k \frac{n}{i}i=1∑n?ikin?,我們可以通過數(shù)論分塊加拉格朗日插值來求。
注意這題拉格朗日插值求值同樣也上記憶化。
代碼
待補(bǔ)......總結(jié)
以上是生活随笔為你收集整理的LibreOJ #6207. 米缇(杜教筛 + 拉格朗日插值)的全部?jī)?nèi)容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: E. Jamie and Tree(树链
- 下一篇: dft + fft(模板)