P3327 [SDOI2015] 约数个数和

最新推荐文章于 2024-09-14 22:25:23 发布

C++oj

最新推荐文章于 2024-09-14 22:25:23 发布

阅读量325

点赞数 2

分类专栏： C++数学题，脑袋发烧！！！文章标签：算法数学莫比乌斯反演 c++ P3327 约数个数和

版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。

本文链接：https://blog.csdn.net/2301_78836126/article/details/141037102

版权

C++数学题，脑袋发烧！！！专栏收录该内容

4 篇文章 0 订阅

订阅专栏

需要莫比乌斯函数知识，主要是推式子。挺难的，但不多。

对了以后的数学题不放代码，除非是知识点讲解，或者算法模板。

因为总有几个烧饼喜欢复制代码刷AC率

这题的难点就是这个约数个数函数。因此，我先给出这个函数的一个重要性质：

$d(ij)=\sum_{x\mid i}^{}\sum_{y\mid j }^{}[gcd(x,y)=1]$

至于证明。。。。。我不会。。。。。

知道式子后就可以“愉快”的推式子了!

首先我们很套路的设几个函数，(若不知道为什么那么设的话。。。反正你看就对了）

$f(d)=\sum_{i=1}^{n}\sum_{j=1}^{m}[gcd(i,j)=d]$

$F(n)=\sum_{n\mid d}^{}f(d)$

由莫比乌斯反演可以得到：

$f(n)=\sum_{n\mid d}^{}\mu(\left \lfloor \frac{d}{n} \right \rfloor)F(d)$

我们设ans为求和最终值；

则 $ans=\sum_{i=1}^{n}\sum_{j=1}^{m}d(ij)$

$ans=\sum_{i=1}^{n}\sum_{j=1}^{m}\sum_{x\mid i}^{}\sum_{y\mid j}^{}[gcd(x,y)=1]$

看到这个[gcd(x,y)=1]，我们就可以根据μ的性质把它带进去.

$ans=\sum_{i=1}^{n}\sum_{j=1}^{m}\sum_{x\mid i}^{}\sum_{y\mid j}^{}\sum_{d\mid gcd(x,y)}^{}\mu (d)$

更换枚举项，由枚举gcd(x,y)的约数，改为直接枚举d.

$ans=\sum_{i=1}^{n}\sum_{j=1}^{m}\sum_{x\mid i}^{}\sum_{y\mid j}^{}\sum_{d=1}^{min(n,m)}\mu (d)[d\mid gcd(x,y)]$

将 $\mu (d)$ 提出来：（因为它与i,j无关）

$ans=\sum_{d=1}^{min(n,m)}\mu (d)\sum_{i=1}^{n}\sum_{j=1}^{m}\sum_{x\mid i}^{}\sum_{y\mid j}^{}[d\mid gcd(x,y)]$

我们再一次更换枚举项，将枚举x,y换为枚举dx,dy。这样[d∣gcd(x,y)]这个条件就可以省去。

$ans=\sum_{d=1}^{min(n,m)}\mu (d)\sum_{x=1}^{\left \lfloor \frac{n}{d} \right \rfloor}\sum_{y=1}^{\left \lfloor \frac{m}{d} \right \rfloor}\left \lfloor \frac{n}{dx} \right \rfloor\left \lfloor \frac{m}{dy} \right \rfloor$

我们可以发现 $\left \lfloor \frac{n}{dx} \right \rfloor$ 与y无关，所以可以提前。

$ans=\sum_{d=1}^{min(n,m)}\mu (d)(\sum_{x=1}^{\left \lfloor \frac{n}{d} \right \rfloor})(\sum_{y=1}^{\left \lfloor \frac{m}{d} \right \rfloor})$

当我们将式子化简成这样的时候，我们已经可以看出，这个式子已经可以做到O(n)计算了。

但是，由于存在多组数据，所以我们就可以运用整除分块，将这个式子优化成O(sqrt(n))的时间复杂度。

这样，这道题就可以A了。!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

还是非常有成就感的。

关注

2
点赞
踩
3

收藏

觉得还不错? 一键收藏
0
评论
复制链接

分享到 QQ

分享到新浪微博

扫一扫

专栏目录

评论

被折叠的条评论为什么被折叠?

到【灌水乐园】发言

查看更多评论

添加红包

成就一亿技术人!

hope_wisdom

发出的红包

实付元

使用余额支付

点击重新获取

扫码支付

钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。