题目描述:
给出a,b,c,d,k 求 a<=x<=b c<=y<=d gcd(x,y)==k的个数
题目分析:
令g(a,b,k)为 1<=x<=a 1<=y<=b gcd(x,y)==k的个数
根据容斥原理 ans=g(b,d,k)-g(a-1,d,k)-g(b,c-1,k)+g(a-1,c-1,k)
满足gcd(x,y)是k的(x,y)的对数也等价于1<=x<=n/k,1<=y<=m/k,(x,y)互质的对数
即 g(a,b,k)==g(a/k,b/k,1)
设f[i]=g(a,b,i)的个数
令F[i]=g(a,b,k*i) (k为正整数)
显然 F[i]=floor(n/i)*floor(m/i)
根据莫比乌斯反演可得
直接枚举是O(n)的.
然而会T!
观察一下 floor(n/i)至多只有2*sqrt(n)个取值,也就是说很长一段区间内值是不变的
于是我们处理mu的前缀和,用分块优化即可.
题目链接:
Ac code:
#include <cstdio>
#include <iostream>
#include <cmath>
#define sqr(x) x*x
#define ll long long
const int maxm=51000;
int mu[maxm],prime[maxm],sum[maxm],cnt,n;
bool vis[maxm];
void getmu()
{
mu[1]=1,vis[1]=1;
for(int i=2;i<=maxm;i++)
{
if(!vis[i])
{
prime[++cnt]=i;
//printf("%d\n",i);
mu[i]=-1;
}
for(int j=1;j<=cnt&&prime[j]*i<=maxm;j++)
{
vis[prime[j]*i]=1;
if(i%prime[j]==0)
{
mu[prime[j]*i]=0;
break;
}
mu[prime[j]*i]=-mu[i];
}
}
for(int i=1;i<=maxm;i++) sum[i]=mu[i]+sum[i-1];
}
int cal(int x,int y)
{
int ans=0,last;
for(int i=1;i<=std::min(x,y);i=last+1)
{
last=std::min(x/(x/i),y/(y/i));
ans+=(sum[last]-sum[i-1])*(x/i)*(y/i);
}
return ans;
}
void work()
{
int a,b,c,d,k;
scanf("%d%d%d%d%d",&a,&b,&c,&d,&k);
a--,c--;
a/=k,b/=k,c/=k,d/=k;
printf("%d\n",cal(b,d)-cal(b,c)-cal(a,d)+cal(a,c));
}
int main()
{
getmu();
int q;
scanf("%d",&q);
while(q--)
work();
return 0;
}