莫比乌斯反演经典例题（2）

链接：P3327 [SDOI2015]约数个数和 - 洛谷 | 计算机科学教育新生态 (luogu.com.cn)
题意：多组数据给定 n,m，求 ${\sum }_{i=1}^n{\sum }_{j=1}^{m}d(i\ast j)$。其中 d 函数为约数和个数函数。
题解：

1. 首先要了解约数和个数的另一表达形式 $d(i\ast j)={\sum }_{x|i}{\sum }_{y|j}[gcd(x,y)==1]$。

   证明：对于任意质数因子 p，我们有 $p^c$ 按先在 i 里取出 $p^c$ ，如果 i 中最多只能取出 $p^a$ ，那么在 j 中需取出 $p^{c-a}$ ，对于每一个 c ，其取法是唯一的。对于每一种取法，其得到的 c 也是唯一的。那么可以这么表示，在 i 中取出 $p^a$ 时表示取出 $p^a$ ，而在 j 中取出 b 时表示 $p^{a+b}$，a 为 i 中 p 因子最高次幂。因此对于任意质数因子 p ，只会从 i 或者 j 中取出，而不会重叠。因此 $gcd(x,y)$ 为 1 。这样通过枚举 x,y ，可以不重不漏的计算出 $d(i\ast j)$ 的所有因数个数。（取自某位大佬，稍加修改）

2. 从 1 可将原始化为 $d(i\ast j)={\sum }_{i=1}^n{\sum }_{j=1}^m{\sum }_{x|i}{\sum }_{y|j}[gcd(x,y)==1]$。

3. 对于 2 式，转化为直接枚举 x,y ，有 $d(i\ast j)={\sum }_{i=1}^n{\sum }_{x|i}{\sum }_{j=1}^m{\sum }_{y|j}[gcd(x,y)==1]={\sum }_{x=1}^n{\sum }_{y=1}^m[gcd(x,y)==1]\lfloor \frac{n}{x}\rfloor \lfloor \frac{m}{y}\rfloor$ 。这么化简的原因是对于每个 x 或者 y ，其都会在它的倍数出现时贡献次数。

4. 此时可以令 $f(z)={\sum }_{x=1}^n{\sum }_{y=1}^m[gcd(x,y)==d]\lfloor \frac{n}{x}\rfloor \lfloor \frac{m}{y}\rfloor$，令 $F(z)={\sum }_{z|d}f(d)$ 。

5. 对于 F(z) 而言是容易快速求出答案的，因为 F(z) 求的是其 $gcd$ 的倍数，那么 $F(z)={\sum }_{i=1}^n{\sum }_{j=1}^m[gcd(i,j)==k\ast z]\lfloor \frac{n}{i}\rfloor \lfloor \frac{m}{j}\rfloor ={\sum }_{i=1}^{n/z}{\sum }_{j=1}^{m/z}[gcd(i,j)==k]\lfloor \frac{n}{i}\rfloor \lfloor \frac{m}{j}\rfloor$ ，其中 k 为任意数，则 $gcd$ 为任意数可直接化去。有 $F(z)={\sum }_{i=1}^{n/z}{\sum }_{j=1}^{m/z}\lfloor \frac{n}{i}\rfloor \lfloor \frac{m}{j}\rfloor$。

6. 对于 5 ，其实还可以化简，由于 $\lfloor \frac{n}{i}\rfloor$ 与 j 无关，则两个除法满足分配率有 $F(z)=({\sum }_{i=1}^{n/z}\lfloor \frac{n}{i}\rfloor )({\sum }_{j=1}^{m/z}\lfloor \frac{m}{j}\rfloor )$ 。

7. 对于 $F(z)$ 可以预处理出来，枚举每一个 n ，然后整出分块求出其答案。

8. 此时回到 4，有反演形式 $f(z)={\sum }_{z|d}\mu (\frac{d}{z})F(d)$ 。则所求为 $f(1)={\sum }_{i=1}^{min(n,m)}\mu (i)F(i)$ 。对于某些 i 的取值 ，F 函数值是一样的，则可以整出分块来求答案。总复杂度为 $O(n\sqrt{n}+T\sqrt{(}n))$ 。

#pragma GCC optimize("Ofast")
#pragma GCC optimize("unroll-loops")
#include<iostream>
#include<algorithm>
#include<vector>
#include<cstring>
#include<functional>
#include<queue>
#include<unordered_map>
#include<map>
#include<set>

using namespace std;
using ll=long long;
using P=pair<int,int>;
const ll inf=1e18;

struct MobiusInversion{
	static constexpr int maxn=5e4+5;
	int cnt;
	vector<int>pri,vs,sm,mu;
	MobiusInversion(int x=maxn):pri(x+5),vs(x+5),sm(x+5),mu(x+5){
		vs[0]=vs[1]=1,mu[1]=1;
		for(int i=2;i<=x;i++)
		{
			if(!vs[i])pri[++cnt]=i,mu[i]=-1;
			for(int j=1;j<=cnt&&i*pri[j]<=x;j++)
			{
				vs[i*pri[j]]=1;
				if(i%pri[j]==0)break;
				mu[i*pri[j]]=-mu[i];
			}
			mu[i]+=mu[i-1];
		}

		for(int i=1;i<=x;i++)
		{
			for(int l=1,r;l<=i;l=r+1)
			{
				r=i/(i/l);
				sm[i]+=(i/l)*(r-l+1);
			}
		}

	}

	ll getans(int n,int m){
		ll ans=0;
		if(n>m)swap(n,m);
		for(int l=1,r;l<=n;l=r+1)
		{
			r=min(n/(n/l),m/(m/l));
			ans+=1ll*(mu[r]-mu[l-1])*sm[n/l]*sm[m/l];
		}
		return ans;
	}

}tt;

void solve()
{
	int n,m; cin>>n>>m;
	cout<<tt.getans(n,m)<<"\n";
}

int main()
{
	ios::sync_with_stdio(false);
	cin.tie(0); cout.tie(0);
	int t; cin>>t;
	while(t--)solve();
	return 0;
}