HDU1695的题目大意是这样的,给你 a , b , c , d , k 五个值 (题目说明了 你可以认为 a=c=1) x 属于 [1,b] ,y属于[1,d] 让你求有多少对这样的 (x,y)满足gcd(x,y)==k。
这一题和hdu 6053不一样的地方在于,在hdu6053 x,y 是确定的一个数而不是一个范围,所以hdu 6053 不是纯正的莫比乌斯,只是单纯的容斥的时候利用到了莫比乌斯函数而已,如果把hdu 6053 强行归类到 正版的莫比乌斯,像我这样的新手就会很纠结。。hdu 6053应该归类到容斥。
在这里 x和 y是一个定值 那么 因为 gcd(x,y)==k 和 gcd(x/k,y/k)==1 其实是等价的 所以 枚举 x/k 和y/k里面的共有因子 就是 1~min(x/k,y/k) 那么就变成了基础的莫比乌斯入门题了,由于会重复,所以要把两个区间的公共部分/2
#include <bits/stdc++.h>
using namespace std;
const int maxn=1e5+7;
bool vis[maxn];
int prime[maxn],mu[maxn];
int cnt;
void Init(){
int N=maxn;
memset(prime,0,sizeof(prime));
memset(mu,0,sizeof(mu));
memset(vis,0,sizeof(vis));
mu[1] = 1;
cnt = 0;
for(int i=2; i<N; i++){
if(!vis[i]){
prime[cnt++] = i;
mu[i] = -1;
}
for(int j=0; j<cnt&&i*prime[j]<N; j++){
vis[i*prime[j]] = 1;
if(i%prime[j]) mu[i*prime[j]] = -mu[i];
else{
mu[i*prime[j]] = 0;
break;
}
}
}
}
// void getMu(){
// int N=maxn;
// for(int i=1;i<N;++i){
// int target=i==1?1:0;
// int delta=target-mu[i];
// mu[i]=delta;
// for(int j=2*i;j<N;j+=i)
// mu[j]+=delta;
// }
// }
int main(){
int T,cas=0;
int a,b,c,d,k;
scanf("%d",&T);
Init();
while(T--)
{
cout<<"Case "<<++cas<<": ";
scanf("%d%d%d%d%d",&a,&b,&c,&d,&k);
if(k==0)
{
cout<<0<<endl;
continue;
}
b/=k;
d/=k;
int minn=min(b,d);
long long ans1=0,ans2=0;
for(int i=1;i<=minn;i++)
{
ans1+=(long long)mu[i]*(b/i)*(d/i);
}
for(int i=1;i<=minn;i++)
{
ans2+=(long long)mu[i]*(minn/i)*(minn/i);
}
cout<<ans1-ans2/2<<endl;
}
}