题目大意:求
看到gcd(i,j)==d 首先想到把i,j同时除d,变成统计互质的数的个数。
令题目转化为求
想统计众多式子中等于1的个数,可以想到莫比乌斯函数的一个性质
(这里的d不同于题目中的d,题目中的d已经没用了)
令上式中的,则代入原式有
将条件转换一下,得到
将式子由枚举i,j转换成枚举d,得到
而(a’/d)与(b’/d)的取值在一定区间内不会发生变化,分块预处理,算法时间复杂度为O(n*sqrt(n)),其中n=min(a’,b’)
#include<cstdio>
#include<algorithm>
#define N 50000
using namespace std;
int mu[N+5],n,ans,prime[N],tot,sum[N];
bool np[N+5];
void init(){
mu[1]=1;
for(int i=2;i<=N;i++) {
if(!np[i]) prime[++tot]=i, mu[i]=-1;
for(int j=1;prime[j]*i<=N;j++){
np[i*prime[j]]=true;
if(!(i%prime[j])) break;
mu[i*prime[j]]=mu[i]*mu[prime[j]];
}
}
return ;
}
int main() {
init();
for(int i=1;i<=N;i++)
sum[i]=sum[i-1]+mu[i];
int T;
scanf("%d",&T);
while(T--) {
int a,b,d;
scanf("%d%d%d",&a,&b,&d);
if(a>b) swap(a,b);
a/=d, b/=d, ans=0;
int l=1,r;
while(l<=a) {
int r=min(a/(a/l), b/(b/l));
ans+=(a/l)*(b/l)*(sum[r]-sum[l-1]);
l=r+1;
}
printf("%d\n",ans);
}
return 0;
}