P3327 [SDOI2015]约数个数和（思维+莫比乌斯反演）

题目传送门
这是一个很有意思的莫反
这道题要你求
$$\sum _{i=1}^n\sum _{j=1}^{m}d(ij)$$
$d()$表示约数的个数

这里有一个结论需要记住就是：
$$d(i,j)=\sum _{x|i}\sum _{y|j}[gcd(x,y)=1]$$
具体证明就是：按情况分类，因数个数就是分情况考虑
如果$ij$中包含$a^k$那么在这两个$\sum$中就是枚举了$k+1$次$a^0$(也就是一次都没有枚举过，也算一次)，因为毕竟是因数和函数，性质我们都不好搞，所以搞一个比较熟悉的$gcd$

那么这个式子就是：
$$\sum _{i=1}^n\sum _{j=1}^m\sum _{x|i}\sum _{y|j}[gcd(x,y)=1]$$通过莫比乌斯反演我们可以知道

$$\sum _{i=1}^n\sum _{j=1}^m\sum _{x|i}\sum _{y|j}\sum _{d|gcd(x,y)}\mu (d)$$
交换枚举顺序
$$\sum _{d=1}^{min(n,m)}\mu (d)\sum _{i=1}^{n/d}\sum _{j=1}^{m/d}\sum _{x|i}\sum _{y|j}$$
继续交换枚举顺序
$$\sum _{d=1}^{min(n,m)}\mu (d)\sum _{x=1}^{n/d}\sum _{y=1}^{m/d}\sum _{i=1}^{n/(d\ast x)}\sum _{j=1}^{m/(d\ast y)}$$
然后得到
$$\sum _{d=1}^{min(n,m)}\mu (d)\sum _{x=1}^{n/d}\sum _{y=1}^{m/d}[\frac{n}{d\ast x}][\frac{m}{d\ast y}]$$
再进行分离，为啥？因为你前后无瓜！
$$\sum _{d=1}^{min(n,m)}\mu (d)(\sum _{x=1}^{n/d}[\frac{n}{d\ast x}])(\sum _{y=1}^{m/d}[\frac{m}{d\ast y}])$$
看这个东西：
$$(\sum _{y=1}^{m/d}[\frac{m}{d\ast y}])$$
让$t=m/d$那不就成了
$$(\sum _{y=1}^t[\frac{t}{y}])$$
对这个进行数论分块，预处理一下

for(ll i=1;i<=N;i++)
{
	for(ll l=1,r;l<=i;l=r+1)
	{
	  	r=i/(i/l);
	  	check[i]+=(r-l+1)*(i/l);
	}
}

因为到最后枚举的是$$\frac{n}{d\ast x}$$
有前面的一个$d$这个问题
所以还要再进行一个数论分块（因为数组访问的，你每回计算的角标有大量的相等的情况）

#include<iostream>
using namespace std;
typedef long long ll;
const ll N = 50005;
ll check[N+5];
ll mu[N + 5], p[N + 5];
bool flg[N + 5];
void mobius() {
	  ll tot = 0;
	  mu[1] = 1;
	  for (ll i = 2; i <= N; ++i) {
	    if (!flg[i]) {
	      p[++tot] = i;
	      mu[i] = -1;
	    }
	    for (ll j = 1; j <= tot && i * p[j] <= N; ++j) {
	      flg[i * p[j]] = 1;
	      if (i % p[j] == 0) {
	        mu[i * p[j]] = 0;
	        break;
	      }
	      mu[i * p[j]] = -mu[i];
	    }
	  }
	  //数论分块要用到的前缀和 
	  for (ll i = 1; i <= N; ++i) mu[i] += mu[i - 1];
	  for(ll i=1;i<=N;i++){
	  	for(ll l=1,r;l<=i;l=r+1){
	  		r=i/(i/l);
	  		check[i]+=(r-l+1)*(i/l);
	  		
		  }//cout<<check[i]<<endl;
	  		//if(i==50)exit(0);
	  }
}
int main(){
	mobius();
	int T;
	cin>>T;
	while(T--){
		ll n,m,ans=0;
		cin>>n>>m;
		for(ll l=1,r;l<=min(m,n);l=r+1){
		//	cout<<check[l]<<" "<<l<<endl;
			r=min(m/(m/l),n/(n/l));
			ans+=(mu[r]-mu[l-1])*(check[m/l])*check[n/l];
		//	cout<<ans<<endl;
		}
		cout<<ans<<endl;
	} 
}