莫比乌斯反演 hdu6390

最新推荐文章于 2020-12-30 10:00:46 发布

animalcoder

最新推荐文章于 2020-12-30 10:00:46 发布

阅读量1k

点赞数

分类专栏： ACM--莫比乌斯反演文章标签： hdu6390 莫比乌斯反演 hdu多校第7场1005 GuGuFishtion

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

本文链接：https://blog.csdn.net/animalcoder/article/details/81639129

版权

ACM--莫比乌斯反演专栏收录该内容

10 篇文章 0 订阅

订阅专栏

hdu多校第7场1005 GuGuFishtion

好久没做莫比乌斯反演了，推柿子走起ooo

$ans=\sum_{i=1}^{m}\sum_{j=1}^{n}\frac{\phi (ij)}{\phi (i)\phi(j)},n,m<=1e6$

$=\sum_{i=1}^{m}\sum_{j=1}^{n}\frac{(i,j)}{\phi ( (i,j) )}$ ,解释：用phi[x]=x*（1-质因子倒数的积）来证

$=\sum_{i=1}^{m}\sum_{j=1}^{n}\sum_{k=1}^{min(n,m)}\frac{k}{\phi(k) }[(i,j)==k]$ ,令i=ik,j=jk

$= \sum_{k=1}^{min(n,m)}\frac{k}{\phi(k) } \sum_{i=1}^{m/k} \sum_{j=1}^{n/k} [(i,j)==1]$

$= \sum_{k=1}^{min(n,m)}\frac{k}{\phi(k) } \sum_{i=1}^{m/k} \sum_{j=1}^{n/k} \sum_{d|(i,j)}\mu (d)$

$= \sum_{k=1}^{min(n,m)}\frac{k}{\phi(k) } \sum_{i=1}^{m/k} \sum_{j=1}^{n/k} \sum_{d=1}^{min(n/k,m/k)} \mu (d)[d|i][d|j]$ ,d求和提前

$= \sum_{k=1}^{min(n,m)}\frac{k}{\phi(k) } \sum_{d=1}^{min(n/k,m/k)} \mu (d) \sum_{i=1}^{m/k} [d|i] \sum_{j=1}^{n/k} [d|j]$ ,令i=id,j=jd,并假设n<m

$= \sum_{k=1}^{n}\frac{k}{\phi(k) } \sum_{d=1}^{n/k} \mu (d) \left \lfloor \frac{n}{kd} \right \rfloor * \left \lfloor \frac{m}{kd} \right \rfloor$ ，令kd=T;交换求和顺序,则k为T的约数

$= \sum_{T=1}^{n} \left \lfloor \frac{n}{T} \right \rfloor \left \lfloor \frac{m}{T} \right \rfloor \sum_{k|T} \mu (T/k) \frac{k }{\phi(k) }$ ，令 $Fuck(T)=\sum_{k|T} \mu (T/k) \frac{k }{\phi(k) }$

$ans= \sum_{T=1}^{n} \left \lfloor \frac{n}{T} \right \rfloor \left \lfloor \frac{m}{T} \right \rfloor Fuck(T)$ ，考虑到Fuck是积性函数

用线性筛预处理Fuck前缀和加上分块即可O(Tn(sqrt(n)+sqrt(m)))解决

$Fuck(p)=\frac{p}{p-1}-1=\frac{1}{p-1}$ ,若p为质数

$Fuck(ap)=Fuck(a)*Fuck(p)$ ,若a%p!=0

$Fuck(ap)=0$ ,若a%p==0,考虑用Fuck(a)&&G(p)构造之

考虑a的质因子f(p^k),跟f(p^k+1)均为0，所以G(p)=0

速度rank2,能吹一天

关注

0
点赞
踩
2

收藏

觉得还不错? 一键收藏
0
评论
莫比乌斯反演 hdu6390

hdu多校第7场1005 GuGuFishtion好久没做莫比乌斯反演了，推柿子走起ooo,解释：用phi[x]=x*（1-质因子倒数的积）来证,令i=ik,j=jk,d求和提前,令i=id,j=jd,并假设n&lt;m，令kd=T;交换求和顺序,则k为T的约数，令，考虑到Fuck是积性函数用线性筛预处理Fuck前缀和加上分块即可O(T*n*(sqrt(...
复制链接

扫一扫

专栏目录

评论

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

到【灌水乐园】发言

查看更多评论

添加红包

成就一亿技术人!

hope_wisdom

发出的红包

实付元

使用余额支付

点击重新获取

扫码支付

钱包余额 0

抵扣说明：

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