标签：25 前缀 Min 质数 最小 sqrt 因子 Sieve sum

这个东西不是人想的。

解决问题：积性函数前缀和。

适用条件：可以快速计算 \(f(p)\) 的前缀和，\(f(p^k)\) 可以被表示成若干完全积性函数的线性组合（指对应项可以快速组合出来）。

时空复杂度：就当是 \(O(\dfrac{n^\frac{3}{4}}{\log n}+n^{1-\epsilon})-O(\sqrt n)\)。

以下默认 \(f(p)\) 为关于 \(p\) 的单项式，\(p_i\) 为第 \(i\) 个质数。

求质数处的前缀和

我们需要对于每个能被表示成 $\lfloor \dfrac{n}{x}\rfloor $ 的数，求所有在其前的质数 \(p\) 的 \(\sum f(p)\)。注意到一个数的最小质因子小等于 \(\sqrt n\)，我们设计一个 dp：\(g_{i,j}\) 表示最小质因子大于 \(p_i\)，\([1,j]\) 的所有数的 \(f\) 之和。转移时把最小质因子为 的 \(p_i\) 数的 \(f\) 全部删去，即得转移，注意容斥掉前面的质数的贡献：

\[g_{i,j} = g_{i-1,j} - f(p_i)\big{(}g_{i-1,j/p_i} - \sum_{k\le i}f(p_k)\big{)} \]

这就体现出完全积性的重要性。求和号部分可以在线筛质数的时候预处理出来。第一维显然可以滚动数组优化掉，而对于第二维，$\lfloor \dfrac{n}{x}\rfloor $ 只有 \(O(\sqrt n)\) 个，预处理出来，大于 \(\sqrt n\) 的编号为 \(n/\sqrt n\) 即可。

我们把筛去所有最小质因子后的得到的 dp 数组记为 \(g\)。

求原函数前缀和

设 \(s_{n,i}\) 表示 \([1,n]\) 最小质因子大于 \(p_i\) 的所有数的 \(f\) 之和。枚举最小质因子及其次数，得到转移

\[s_{n,i} = g_n - \sum_{j\le i}f(p_j) + \sum\limits_{p_k^t \le n} f(p_k^t)\bigg{(}s_{n/p_k^t},k+1)+[t\gt 1]\bigg{)} \]

根据神奇的复杂度分析，可以直接递归下去做。

最后加上 \(1\) 的贡献。

标签：25,前缀,Min,质数,最小,sqrt,因子,Sieve,sum
来源： https://www.cnblogs.com/defK/p/16097216.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。