Luogu P2257 YY的GCD

题目大意

给定 $n,m$，求$$\sum _{i=1}^n\sum _{j=1}^m[(i,j)\in prime]$$每个测试点有 $T$ 组数据。
数据范围　$1⩽n,m⩽10^7,1⩽T⩽10^4$

题解

先把 $(i,j)\in prime$ 提出来
$$\sum _{p\in prime}^{\min (n,m)}\sum _{i=1}^n\sum _{j=1}^m[(i,j)=p]$$

然后常规套路后面两个 $\text{sigma}$。
$$\begin{array}{rl}原式& =\sum _{p\in prime}^{\min (n,m)}\sum _{i=1}^{\lfloor \frac{n}{p}\rfloor }\sum _{j=1}^{\lfloor \frac{m}{p}\rfloor }[(i,j)=1]\\ & =\sum _{p\in prime}^{\min (n,m)}\sum _{i=1}^{\lfloor \frac{n}{p}\rfloor }\sum _{j=1}^{\lfloor \frac{m}{p}\rfloor }\sum _{d|i,d|j}\mu (d)\\ & =\sum _{p\in prime}^{\min (n,m)}\sum _{d=1}^{\min (\lfloor \frac{n}{p}\rfloor ,\lfloor \frac{m}{p}\rfloor )}\mu (d)(\sum _{i=1}^{\lfloor \frac{n}{p}\rfloor }[d|i]\left)\right(\sum _{j=1}^{\lfloor \frac{m}{p}\rfloor }[d|j])\\ & =\sum _{p\in prime}^{\min (n,m)}\sum _{d=1}^{\min (\lfloor \frac{n}{p}\rfloor ,\lfloor \frac{m}{p}\rfloor )}\mu (d)\lfloor \frac{n}{pd}\rfloor \lfloor \frac{m}{pd}\rfloor \end{array}$$

如此设 $n$ 以内有 $\pi (n)$ 个质数，则运行一次的时间复杂度为 $$T(n)=\Theta (\pi (n)\times \sqrt{\min (n,m)})$$根据素数定理，有 $\pi (n)\sim \frac{1}{\ln n}$，因而程序的时间复杂度为 $$T(n)=\Omega (n+T\times \frac{\sqrt{\min (n,m)}}{\ln n})$$

按照题目的数据范围，计算得结果约为 $6214206884.3>5\times 10^7$，而且素数定理的估计相对于我们的复杂度计算是乐观的，也就是说，结果可能会更大。因而我们只能考虑优化。

可以发现，第一个 $\text{sigma}$ 的枚举顺序是没有前途的，是没有前途的，是没有前途的……因此我们要更换枚举顺序，至少不能让第一个 $\text{sigma}$ 那么没有前途……尝试枚举 $k=pd$，则
$$原式=\sum _{k=1}^{\min (n,m)}\sum _{p\in prime,p|k}^{\min (n,m)}\mu \left(\frac{k}{p}\right)\lfloor \frac{n}{pd}\rfloor \lfloor \frac{m}{pd}\rfloor$$

令 $$f(n)=\sum _{p\in prime,p|k}^n\mu \left(\frac{k}{p}\right)$$则 $$原式=\sum _{k=1}^{\min (n,m)}f(\text{min}(n,m))\lfloor \frac{n}{pd}\rfloor \lfloor \frac{m}{pd}\rfloor$$

现在如果我们能够得到 $f$ 的前缀和，那么就可以用除法分块了。