定义 Q(n)=∑i=1⌊n⌋μ(i)⌊ni2⌋Q(n) = \sum\limits_{i=1}^{\lfloor \sqrt n \rfloor} \mu(i) \left\lfloor \dfrac n {i^2} \right\rfloorQ(n)=i=1∑⌊n⌋μ(i)⌊i2n⌋,即 [1,n][1,n][1,n] 中无平方因子数的个数。
// summu 为莫比乌斯函数前缀和 ll Q(ll n) { // [1, n] 中无平方因子数的个数 ll ans = 0; for (ll l = 1, r; l <= sqrt(n); l = r + 1) { ll val = n / (l * l); r = sqrt(n / val); ans += val * (summu[r] - summu[l-1]); } return ans; }
(浅色调 的题解提到了这一内容,但是给出了错误的时间复杂度)
所以此题答案约为输入的 π26\dfrac {\pi^2} 66π2 倍。
这个数约为 1.6449340671.6449340671.644934067,第一篇题解中的二分上限 164493408116449340811644934081 就差不多是这个数。
这也是 P8883 幻想中成为原神 的结论。
参考资料:https://zhuanlan.zhihu.com/p/590842929