[HAOI2011]Problem b 莫比乌斯函数

最新推荐文章于 2021-08-13 10:48:47 发布

冷酷|射手

最新推荐文章于 2021-08-13 10:48:47 发布

阅读量146

点赞数

分类专栏：数学推导公式

本文链接：https://blog.csdn.net/fangfengwei/article/details/113417673

版权

数学同时被 2 个专栏收录

6 篇文章 0 订阅

订阅专栏

推导公式

5 篇文章 0 订阅

订阅专栏

这篇博客介绍了如何利用莫比乌斯函数解决[HAOI2011]Problem b的竞赛问题。通过推导公式，提出了利用前缀和计算gcd(i, j)=k的和，并给出了线性筛计算莫比乌斯函数及整除分块计算的方法，实现了O(n+Tn)的时间复杂度解决方案。" 105801100,6649953,Redis分布式锁的误解锁问题及解决方案,"['分布式系统', 'Redis', '数据库', '并发控制']

摘要由CSDN通过智能技术生成

题目链接

[HAOI2011]Problem b

题目大意

$\sum_{i=x}^n\sum_{j=y}^m[\gcd(i,j)=k]\quad(1\le T,x,y,n,m,k\le 10^5)$

推柿子

考虑前缀和计算，令 $f(n,m)=\sum_{i=1}^n\sum_{j=1}^m[\gcd(i,j)=k]$ ， $\sum_{i=x}^n\sum_{j=y}^m[\gcd(i,j)=k]=f(n,m)-f(x-1,m)-f(n,y-1)+f(x-1,y-1)$ 。

考虑计算 $f(n,m)\quad(n\le m)$ 。
$\begin{aligned}f(n,m)&=\sum_{i=1}^{n}\sum_{j=1}^{m}[\gcd(i,j)=k]\\&=\sum_{i=1}^{\lfloor\frac{n}{k}\rfloor}\sum_{j=1}^{\lfloor\frac{m}{k}\rfloor}[\gcd(i,j)=1]\\&=\sum_{i=1}^{\lfloor\frac{n}{k}\rfloor}\sum_{j=1}^{\lfloor\frac{m}{k}\rfloor}\sum_{d|\gcd(i,j)}\mu(d)(\sum_{d|n}\mu(d)=[n=1])\\&=\sum_{d=1}^{\lfloor\frac{n}{k}\rfloor}\mu(d)\sum_{i=1}^{\lfloor\frac{n}{k}\rfloor}d|i\sum_{j=1}^{\lfloor\frac{m}{k}\rfloor}d|j\\&=\sum_{d=1}^{\lfloor\frac{n}{k}\rfloor}\mu(d)\lfloor\frac{n}{kd}\rfloor\lfloor\frac{m}{kd}\rfloor\\\end{aligned}$
$O (n)$ 线性筛计算 $\mu$ ，整除分块计算 $\sum_{d=1}^{\lfloor\frac{n}{k}\rfloor}\mu(d)\lfloor\frac{n}{kd}\rfloor\lfloor\frac{m}{kd}\rfloor$ 。时间复杂度 $O(n+T\sqrt{n})$ 。