题意:
给定n,m,求公式,d(x)为x的约数个数
思路:
这题比较适合反演的学习,因为要用到反演经常用的技巧(公式)
首先需要将d(ij)进行变化,这里有公式:
( 这个公式还可以推到两个,三个的乘积 )
那么就可以推到:
然后就出了gcd()==1的式子,这时候就可以用莫比乌斯反演经常会用到的公式之一:
将公式变化整理一下为:
下一步也要用到莫比乌斯反演经常要用到的技巧——转换枚举的变量:我们可以将后三个求和式转换成枚举gcd(n,m)的因子d,然后统计d对应的贡献:
然后继续转换枚举变量:枚举所有d,统计它对应的贡献:
到这一步就能发现后面的两个求和是一个东西,而且都可以在O( n*sqrt(n) )内预处理出来
设这个式子为,可以发现g(x+1) - g(x) = d(x+1),利用这个递推式可以利用类似于埃氏筛的方法更快的在O(n log(n) )内完成g(x)的预处理。
最后就可以对每个询问在O( sqrt(n) ) 内得出结果。
代码:
#include<bits/stdc++.h>
using namespace std;
const int MAXN = 5e4+10;
bool Np[MAXN];
vector<int>pr;
int mu[MAXN];
int sum[MAXN];
long long F[MAXN];
long long D[MAXN];
void init(){
mu[1] = sum[1] = 1;
for(int i=2;i<MAXN;i++){
if(!Np[i]){
pr.push_back(i);
mu[i] = -1;
}for(int j=0;j<pr.size();j++){
int t = pr[j] * i;
if(t>=MAXN)break;
Np[t] = true;
if(i%pr[j]==0){
mu[t] = 0;
break;
}mu[t] = -mu[i];
}sum[i] = sum[i-1] + mu[i];
}
for(int n=1;n<MAXN;n++){
for(int j=n;j<MAXN;j+=n){
D[j]++;
}F[n] = F[n-1] + D[n];
}
}
int main()
{
init();
int T;
scanf("%d",&T);
int n,m;
while(T--){
scanf("%d%d",&n,&m);
long long ans = 0;
if(n>m)swap(n,m);
for(int g=1,last;g<=n;g=last+1){
last = min( n/(n/g) , m/(m/g) );
int a = n/g , b = m/g;
ans += 1LL * (sum[last]-sum[g-1]) * F[a] * F[b];
}printf("%lld\n",ans);
}return 0;
}