#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#define MAXN 50000
#define qc std::ios::sync_with_stdio(false),std::cin.tie(0)
using namespace std;
int n,q,t,prime[MAXN+5],phi[MAXN+5],mo[MAXN+5];
bool isf[MAXN+5];
void iniz_prime(int n){
memset(isf,true,sizeof(isf));
phi[1]=1,isf[0]=isf[1]=false,mo[1]=1;
for (int i=2;i<=n;i++){
if(isf[i]){
phi[i]=i-1;
prime[++prime[0]]=i;
mo[i]=-1;
}
for(int j=1;j<=prime[0]&&i*prime[j]<=n;j++){
isf[i*prime[j]]=false;
if(i%prime[j]==0){
phi[i*prime[j]]=phi[i]*prime[j];
mo[i*prime[j]]=0;
break;
}else{
phi[i*prime[j]]=phi[i]*(prime[j]-1);
mo[i*prime[j]]=-mo[i];
}
}
}
for(int i=2;i<=n;i++){
mo[i]+=mo[i-1];
}
}
int main(){
qc;
iniz_prime(MAXN);
int a,b,d;
int T;
cin>>T;
while (T--)
{
cin>>a>>b>>d;
a/=d,b/=d;
if(a>b)swap(a,b);
int ans=0,pos;
for(int i=1;i<=a;i=pos+1){
pos=min(a/(a/i),b/(b/i));
ans+=(mo[pos]-mo[i-1])*(a/i)*(b/i);
}
cout<<ans<<"\n";
}
}
Neko and function HDU-6537(莫比乌斯反演）
©著作权归作者所有,转载或内容合作请联系作者
- 文/潘晓璐 我一进店门，熙熙楼的掌柜王于贵愁眉苦脸地迎上来，“玉大人，你说我怎么就摊上这事。” “怎么了？”我有些...
- 文/花漫 我一把揭开白布。 她就那样静静地躺着，像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上，一...
- 文/苍兰香墨 我猛地睁开眼，长吁一口气：“原来是场噩梦啊……” “哼！你这毒妇竟也来了？” 一声冷哼从身侧响起，我...
推荐阅读更多精彩内容
- The programs we have looked at thus far have dealt with t...
- flow control if: The condition judegement allow one varia...
- 题目：http://www.lydsy.com/JudgeOnline/problem.php?id=2820 这...
