洛谷P3327：[SDOI2015]约数个数和（莫比乌斯反演）

枚举倍数的一种灵活的变形：$g(d)={\sum }_{d|i}^{n}f(i)={\sum }_{i=1}^{\lfloor \frac{n}{d}\rfloor }f(i\cdot d)$
很显然，但有时能发挥大作用。
其实本质还是要理解西格玛究竟是在算什么

解析

一个重要的前置结论：
$$d(i\cdot j)=\sum _{x|i,y|j}[\gcd (x,y)=1]$$
证明：
尝试把所有互质的数对 $(x,y)$ 和 $i\cdot j$ 的所有因子建立双射关系。首先，对于只有 $i,j$ 一方有的质因子，直接归给有该项质因子的那一方对应的 $x$ 或 $y$。对于双方都有的质因子 $p$，若 $i$ 的对应次数为 $a$，$j$ 的对应次数为 $b$，如果 $i\cdot j$ 的某个因子 $p$ 的次数 $c$ 满足 $c\le a$，就当成 $x$ 含有 $p^c$ 的因子；若 $c>a$，就当成 $y$ 含有 $p^{c-a}$ 次方的因子。这样我们就建立起了双射关系，且所有的 $x,y$ 都是互质的，证毕。

我们再化一下这个 $d$ 函数：
$$\begin{gathered} d(i,j)=\sum _{x|i}\sum _{y|j}[\gcd (i,j)=1] \\ =\sum _{x|i}\sum _{y|j}\sum _{k|\gcd (x,y)}\mu (k) \\ =\sum _{k|i,k|j}\mu (k)\sum _{k|x,x|i}\sum _{k|y,y|j}1 \\ =\sum _{k|i,k|j}\mu (k)\sum _{x|\frac{i}{k}}^{\frac{i}{k}}\sum _{y|\frac{j}{k}}^{\frac{j}{k}}1 \\ =\sum _{k|i,k|j}\mu (k)d(\frac{i}{k})d(\frac{j}{k}) \end{gathered}$$
带回原式：
$$\begin{gathered} \sum _{i=1}^n\sum _{j=1}^{m}d(ij)=\sum _{i=1}^n\sum _{j=1}^m\sum _{k|i,k|j}\mu (k)d(\frac{i}{k})d(\frac{j}{k}) \\ =\sum _k\mu (k)\sum _{k|i}^n\sum _{k|j}^{m}d(\frac{i}{k})d(\frac{j}{k}) \\ =\sum _k\mu (k)\sum _{i=1}^{\lfloor \frac{n}{k}\rfloor }\sum _{j=1}^{\lfloor \frac{m}{k}\rfloor }d(i)d(j) \\ =\sum _k\mu (k)(\sum _{i=1}^{\lfloor \frac{n}{k}\rfloor }d(i))\cdot (\sum _{j=1}^{\lfloor \frac{m}{k}\rfloor }d(j)) \end{gathered}$$
线性筛出 $\mu ,d$，预处理出前缀和，整除分块即可。

代码

#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define ull unsigned long long
#define debug(...) fprintf(stderr,__VA_ARGS__)
inline ll read(){
	ll x(0),f(1);char c=getchar();
	while(!isdigit(c)) {if(c=='-')f=-1;c=getchar();}
	while(isdigit(c)) {x=(x<<1)+(x<<3)+c-'0';c=getchar();}
	return x*f;
}
const int N=1e5+100;
int n,m;

int p[N],tot,vis[N],mu[N],msum[N],d[N];
ll dsum[N];
void init(int n){
	mu[1]=1;d[1]=1;
	for(int i=2;i<=n;i++){
		if(!vis[i]) p[++tot]=i,mu[i]=-1,d[i]=2;
		for(int j=1;j<=tot&&p[j]<=n/i;j++){
			vis[i*p[j]]=1;
			if(i%p[j]==0){
				mu[i*p[j]]=0;
				d[i*p[j]]=d[i]+d[i]-d[i/p[j]];
				break;
			}
			mu[i*p[j]]=-mu[i];
			d[i*p[j]]=d[i]*2;
		}
	}
	for(int i=1;i<=n;i++){
		msum[i]=msum[i-1]+mu[i];
		dsum[i]=dsum[i-1]+d[i];
	}
	return;
}

signed main() {
#ifndef ONLINE_JUDGE
	//freopen("a.in","r",stdin);
	//freopen("a.out","w",stdout);
#endif
	init(5e4);
	int T=read();
	while(T--){
		n=read();m=read();
		ll ans(0);
		for(int l=1,r;l<=min(n,m);l=r+1){
			r=min(n/(n/l),m/(m/l));
			ans+=(msum[r]-msum[l-1])*dsum[n/l]*dsum[m/l];
		}
		printf("%lld\n",ans);
	}
	return 0;
}
/*
5
5
8
1
46
3564
*/