[HAOI2011]Problem b

题目描述：

给出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的前缀和，用分块优化即可.

题目链接：

Luogu 2522
BZOJ 2301

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;
}