洛谷 P3327 [SDOI2015]约数个数和 (莫比乌斯反演)

题目描述

设$d(x)$为$x$的约数个数，给定$N$、$M$，求 $\sum^N_{i=1}\sum^M_{j=1}d(ij)$

---

输入输出格式

输入格式：
输入文件包含多组测试数据。第一行，一个整数T，表示测试数据的组数。接下来的T行，每行两个整数N、M。

输出格式：
T行，每行一个整数，表示你所求的答案。

---

输入输出样例

输入样例#1：

2
7 4
5 6

输出样例#1：

110
121

---

说明

1<=N, M<=50000

1<=T<=50000

---

题目分析

数论题对我这种智障来说就是要命的。这题我想了一会还是没有想出来。

某Cla口中的水题对蒟蒻我来说都是难题，所谓的好题就是我根本不会做的题。。。

本题必须用到的结论：

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

证明很简单：

将一个数唯一分解为$p_1^{k_1}p_2^{k_2}p_3^{k_3}...p_n^{k_n}$，则其约数个数为$(k_1+1)(k_2+1)(k_3+1)...(k_n+1)$。因为对于任意$i\neq j$，$(p_i,p_j)=1$。

然后考虑一个质数$p$对$d(ij)$的贡献。假设$i$的质因数分解中有$k$个$p$，$j$的质因数分解中有$q$个$p$，那么对 d(