题解：SDOI2015 约数个数和 【莫比乌斯反演】

$$求\sum _{i=1}^n\sum _{j=1}^{m}d(ij)[d(x)表示x的约数个数]$$
$$首先，我们要知道如何去求d(ij)，有这么一个式子$$

$$d(ij)=\sum _{x|i}\sum _{y|j}[gcd(i,j)==1]$$
$$于是，我们的式子就变成了要求\sum _{i=1}^n\sum _{j=1}^m\sum _{x|i}\sum _{y|j}[gcd(i,j)==1]$$
$$改变求和顺序，先枚举因数x和y$$

$$于是我们得到了\sum _{x=1}^n\sum _{y=1}^m\lfloor \frac{n}{x}\rfloor \lfloor \frac{m}{y}\rfloor [gcd(x,y)==1]$$
到了这一步，我们就可以开始上莫比乌斯反演了

$$设f(x)=\sum _{i=1}^n\sum _{j=1}^m\lfloor \frac{n}{i}\rfloor \lfloor \frac{m}{j}\rfloor [gcd(i,j)==x]$$

$$g(n)=\sum _{n|d}f(d)$$
$$于是g(x)=\sum _{i=1}^n\sum _{j=1}^m\lfloor \frac{n}{i}\rfloor \lfloor \frac{m}{j}\rfloor [x|gcd(i,j)]$$

$$同样我们把x提出来消除gcd，于是我们得到$$

$$g(x)=\sum _{i=1}^{\lfloor \frac{n}{x}\rfloor }\sum _{j=1}^{\lfloor \frac{m}{x}\rfloor }\lfloor \frac{n}{ix}\rfloor \lfloor \frac{m}{jx}\rfloor$$
$$把g(x)带回原来的式子$$

$$f(n)=\sum _{n|d}g(d)\mu \left(\frac{d}{n}\right)$$
$$f(1)=\sum _{1|d}g(d)\mu \left(d\right)=\sum _{d=1}^{min(n,m)}\mu (d)\sum _{i=1}^{\lfloor \frac{n}{d}\rfloor }\sum _{j=1}^{\lfloor \frac{m}{d}\rfloor }\lfloor \frac{n}{id}\rfloor \lfloor \frac{m}{jd}\rfloor$$
$$到了这一步之后我们发现，由于\lfloor \frac{n}{id}\rfloor 和枚举y没什么关系，于是又可以把它提前$$

$$\sum _{d=1}^{min(n,m)}\mu (d)\sum _{i=1}^{\lfloor \frac{n}{d}\rfloor }\lfloor \frac{n}{id}\rfloor \sum _{j=1}^{\lfloor \frac{m}{d}\rfloor }\lfloor \frac{m}{jd}\rfloor$$
$$这个复杂的长的很的式子看得十分难受，于是我们就记$$

$$p=\lfloor \frac{n}{id}\rfloor ，q=\lfloor \frac{m}{jd}\rfloor$$

$$于是两个就是很明显的数论分块了，算出来之后再合并起来$$
$$Ans=\sum _{d=1}^{min(n,m)}\mu (d)\ast p\ast q$$

到这里，就完全推完了，整个的复杂度就是$\sqrt{n}\ast \sqrt{n}=n$

于是上一波代码

#include <algorithm>
#include <iostream>
#include <cstring>
#include <cstdio>
using namespace std;
#define re register
#define gc getchar()
#define ll long long
inline int read()
{
	re int x(0); re char ch(gc);
	while(ch>'9'||ch<'0') ch=gc;
	while(ch>='0'&&ch<='9') x=(x*10)+(ch^48),ch=gc;
	return x;
}
const int N=5e4+10;
int mu[N],pri[N],f[N],vis[N],cnt;
void get_mu(int n)	//线性筛求mu函数
{
	mu[1]=1;
	n-=5;
	for(int i=2;i<=n;++i)
	{
		if(!vis[i])
			mu[i]=-1,pri[++cnt]=i;
		for(int j=1;j<=cnt&&i*pri[j]<=n;++j)
		{
			vis[i*pri[j]]=1;
			if(i%pri[j]==0) 
			{
				mu[i*pri[j]]=0;
				break;
			}
			else 
				mu[i*pri[j]]=-mu[i];
		}
	}
	for(int i=1;i<=n;++i)
		mu[i]+=mu[i-1];
	for(int x=1;x<=n;++x)
	{
		ll res=0;
        for(int l=1,r;l<=x;l=r+1) 
			r=x/(x/l),res+=1LL*(r-l+1)*(x/l);
		f[x]=res;
	}
}
void work()
{
	int a=read(),b=read();
	int n=min(a,b);
	ll ans=0;
	for(int l=1,r;l<=n;l=r+1)	//数论分块优化
	{
		r=min(a/(a/l),b/(b/l));
        ans+=1LL*(mu[r]-mu[l-1])*f[a/l]*f[b/l];
	}
	cout<<ans<<endl;
}
int main()
{
	int T=read();
	get_mu(N);
	while(T--) work();
	return 0;
}