[国家集训队2012] 和与积

一道数论好题（突然发现我还没写莫比乌斯反演的博客，先挖个坑吧。。）

题目描述

给出N，统计满足下面条件的数对(a,b)的个数：

$1\le a<b\le N1\le a<b\le N$
$a+b$整除$a\times b$

输入格式

一行一个数$N$

输出格式

一行一个数表示答案

输入输出样例

输入 #1

15

输出 #1

4

题目解析

首先，这题我在刷数学题的时候看到过（确信）。。。。

所以这题是数竟题 ，不过当时数学书上N只有50，并且要写出所有解，并且数学书上的计算方法是O($N^2$)的（因为口算时人脑对于一些枚举和为k的互质的一对数的过程计算起来更快，所以比单纯枚举要快）。

好了进入正题，根据数学书上写的：

设$k=gcd(a,b),a=kx,b=ky$
那么$(a+b)|ab⟹(x+y)|kxy$

又$gcd(x,y,k)=1$
易知$gcd(x+y,xy)=1⟹(x+y)!|xy$
所以$(x+y)|k$
然后数学书上就开始枚举$k$了。。。

但是OI很神奇，所以我们接下来继续推式子

$$ans=\sum _{x=1}^n\sum _{y=x+1}^n[gcd(x,y)==1]\lfloor \frac{n}{max(x,y)(x+y)}\rfloor$$

关于$$\lfloor \frac{n}{max(x,y)(x+y)}\rfloor$$感性理解 ，是因为$1\le k\le \frac{n}{max(x,y)(x+y)}$，$k\in N$所以$k$一共有$\lfloor \frac{n}{max(x,y)(x+y)}\rfloor$种取值。

然后因为$y>x，x+y\le k$
所以$y^2<y(x+y)\le ky=b\le n$
所以$y<\sqrt{n}$
$$ans=\sum _{x=1}^{\sqrt{n}}\sum _{y=x+1}^{\sqrt{n}}[gcd(x,y)==1]\lfloor \frac{n}{max(x,y)(x+y)}\rfloor$$
然后我们就得到了优秀的$O(NlogN)$的算法

之后的过程开始加速毒瘤化（莫比乌斯反演），带好护目镜。
因为$$[x==1]=e(x)=\sum _{d|x}\mu (d)$$
所以!!$$ans=\sum _{x=1}^{\sqrt{n}}\sum _{y=x+1}^{\sqrt{n}}\lfloor \frac{n}{max(x,y)(x+y)}\rfloor \sum _{d|gcd(x,y)}\mu (d)$$
把最后一个$\sum$掉换到前面
$$ans=\sum _{\sigma =1}^{\sqrt{n}}\mu (\sigma )\sum _{x=1}^{\frac{\sqrt{n}}{\sigma }}\sum _{y=x+1}^{\frac{\sqrt{n}}{\sigma }}\lfloor \frac{n}{{\sigma }^{2}y(x+y)}\rfloor$$
但是这里的复杂度就很迷模，这里运用调和求和法（详见具体数学）$$\hat{H}=\frac{\sqrt{n}}{1}+\frac{\sqrt{n}}{2}+\frac{\sqrt{n}}{3}+......+\frac{\sqrt{n}}{\sqrt{n}}=lo{g}_{e}n+\gamma +O(\frac{1}{n})$$
这里的$\gamma$是欧拉常数，所以复杂度为$O(\sqrt{n}lo{g}^{2}n)$，简直是快飞了，但是还是可以考虑优化，设$z=x+y$（因为很容易发现对于$x$的枚举并没啥用）
$$ans=\sum _{\sigma =1}^{\sqrt{n}}\mu (\sigma )\sum _{y=1}^{\frac{\sqrt{n}}{\sigma }}\sum _{z=y+1}^{2y-1}\lfloor \frac{n}{{\sigma }^{2}yz}\rfloor$$
这里同样用调和级数计算时间复杂度得到$O(\sqrt{n}logn)$，然后就A掉啦！！