标签：frac log 复杂度 powerful number sqrt 逆元 sum

心血来潮跑来实现以下这个东西

我们应该知道杜教筛的理论是 \(f * g=h\)，那么问题在于如何找 \(g\)。

之前的blog应该提到过可以令 \(g(p)=-f(p)\)，这样一来 \(h\) 就只会在PN处有值。于是可以大力爆搜 \(h\)，而 \(g\) 的块筛又很好处理。

但是这样复杂度会有一个下限为 \(O(n^{\frac 2 3})\)，有没有办法去除呢？

办法是有的，反过来，设 \(h * g=f\)。

此时我们构造 \(g(p)=f(p)\) 即可得到和上面相同的结论，但此时只需处理 \(g\) 的块筛即可，复杂度下降至 \(O(\frac{n^{\frac 3 4}}{\log n})\) 或者更低。

问题来了，当 \(k>2\) 时，\(g(p^k)\) 应该是多少？

实际上是多少否无所谓，因为有 \(\sum_{i=0}^k h(p^i)g(p^{k-i})=f(p^k)\)，一般情况令 \(g(p^k)=0\)。

但是我在实现的时候 推 错 了 \(g\) 的 块 筛 柿 子，懒得重新推。又发现我的柿子实际上是令 \(g(p^k)=f(p)^k\)，如果直接暴力做卷积那么复杂度会变成 \(O(\sqrt n\log n)\)，于是来优化一下：

\[\sum_{i=0}^kh(p^i)f(p)^{k-i}=f(p^k) \]

\[f(p)^k\sum_{i=0}^kh(p^i)f(p)^{-i}=f(p^k) \]

于是处理出 \(\frac{f(p^k)}{f(p)^k}\)，然后来个差分求逆元就好，复杂度变回了 \(O(\sqrt n)\)，不过需要存在逆元才行。

但是？

\[h(p^k)f(p)^{-k}=\frac{f(p^k)}{f(p)^k}-\frac{f(p^{k-1})}{f(p)^{k-1}} \]

\[h(p^k)=f(p^k)-f(p^{k-1})f(p) \]

不需要逆元也可以。

代码暂时咕了（

标签：frac,log,复杂度,powerful,number,sqrt,逆元,sum
来源： https://www.cnblogs.com/lmpp/p/15827671.html

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