SDOI 2015 约数个数和 题解

题目传送门

题目大意： 求${\sum }_{i=1}^n{\sum }_{j=1}^{m}d(ij)$，$d(x)$ 是 $x$ 的因数个数。

题解

首先要知道 $d$ 的一个性质，有了这个性质就可以愉快地推柿子了！
$$\begin{gathered} \sum _{i=1}^n\sum _{j=1}^{m}d(ij) \\ =\sum _{i=1}^n\sum _{j=1}^m\sum _{x|i}\sum _{y|j}[\gcd (x,y)=1] \end{gathered}$$

更改枚举对象，于是有
$$=\sum _{x=1}^n\sum _{y=1}^m[\gcd (x,y)=1]\lfloor \frac{n}{x}\rfloor \lfloor \frac{m}{y}\rfloor$$

套个莫比乌斯反演（强行安利博客qwq），于是变成
$$\begin{gathered} =\sum _{x=1}^n\sum _{y=1}^m\lfloor \frac{n}{x}\rfloor \lfloor \frac{m}{y}\rfloor \sum _{d|gcd(x,y)}\mu (d) \\ =\sum _{x=1}^n\sum _{y=1}^m\lfloor \frac{n}{x}\rfloor \lfloor \frac{m}{y}\rfloor \sum _{d=1}^{\min (n,m)}\mu (d)[d|\gcd (x,y)] \\ =\sum _{d=1}^{\min (n,m)}\mu (d)\sum _{x=1}^n\sum _{y=1}^m\lfloor \frac{n}{x}\rfloor \lfloor \frac{m}{y}\rfloor [d|\gcd (x,y)] \\ =\sum _{d=1}^{\min (n,m)}\mu (d)\sum _{x=1}^{\lfloor \frac{n}{d}\rfloor }\sum _{y=1}^{\lfloor \frac{m}{d}\rfloor }\lfloor \frac{n}{dx}\rfloor \lfloor \frac{m}{dy}\rfloor \\ =\sum _{d=1}^{\min (n,m)}\mu (d)\sum _{x=1}^{\lfloor \frac{n}{d}\rfloor }\lfloor \frac{n}{dx}\rfloor \sum _{y=1}^{\lfloor \frac{m}{d}\rfloor }\lfloor \frac{m}{dy}\rfloor \end{gathered}$$

设 $g(x)={\sum }_{i=1}^x\lfloor \frac{x}{i}\rfloor$，那么柿子变成
$$=\sum _{d=1}^{\min (n,m)}\mu (d)\times g(\lfloor \frac{n}{d}\rfloor )\times g(\lfloor \frac{m}{d}\rfloor )$$

先用除法分块预处理 $g$，然后线性筛出 $\mu$ ，在求个前缀和，最后上面这条式子用除法分块搞一搞即可。

代码如下：

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
#define maxn 50010
#define ll long long

int t,n,m;
ll ans;
ll mu[maxn],g[maxn];
int prime[maxn],tt=0;
bool v[maxn];
void work()
{
	mu[1]=1;
	for(int i=2;i<=maxn-10;i++)//线性筛mu
	{
		if(!v[i])prime[++tt]=i,mu[i]=-1;
		for(int j=1;j<=tt&&i*prime[j]<=maxn-10;j++)
		{
			v[i*prime[j]]=true;
			if(i%prime[j]==0)break;
			mu[i*prime[j]]=mu[i]*mu[prime[j]];
		}
	}
	for(int i=2;i<=maxn-10;i++)//求前缀和
	mu[i]+=mu[i-1];
	int l,r;
	for(int i=1;i<=maxn-10;i++)//求g
	{
		l=1;
		while(l<=i)
		{
			r=i/(i/l);
			g[i]+=(ll)(r-l+1)*(i/l);
			l=r+1;
		}
	}
}

int main()
{
	scanf("%d",&t);
	work();
	while(t--)
	{
		scanf("%d %d",&n,&m);
		ans=0;
		int l=1,r;
		while(l<=min(n,m))//除法分块求解
		{
			r=min(n/(n/l),m/(m/l));
			ans+=(mu[r]-mu[l-1])*g[n/l]*g[m/l];
			l=r+1;
		}
		printf("%lld\n",ans);
	}
}