bzoj2301 [HAOI2011]Problem b

题目链接：bzoj2301
题目大意：
对于给出的n个询问，每次求有多少个数对(x,y)，满足a≤x≤b，c≤y≤d，且gcd(x,y) = k，gcd(x,y)函数为x和y的最大公约数。

题解：
分块+莫比乌斯反演
原式即

$$\sum_{a≤x≤b}\sum_{c≤y≤d}\sum_{gcd(x,y)=k}1$$
同样的，化式子: $$\sum_{\frac{a}{k}≤ \frac{x}{k} ≤ \frac{b}{k}}\sum_{ \frac{c}{k} ≤ \frac{y}{k} ≤ \frac{d}{k}}\sum_{ gcd(\frac{x}{k},\frac{y}{k})=1 }1$$
设 $\frac{x}{k}$为 $i$，$\frac{y}{k}$为 $j$，把区间$[\frac{a}{k}, \frac{b}{k}]$、 $[\frac{c}{k}, \frac{d}{k}]$容斥一下，就只需计算 $[1,n]$这样的范围了。
于是式子就变成了： $$\sum_{1≤i≤n}\sum_{1≤j≤m}\sum_{gcd(i,j)=1}1$$
即 $$\sum_{1≤i≤n}\sum_{1≤j≤m}\epsilon(gcd(i,j))$$
因为有 $1*\mu=\epsilon$，即 $\sum_{d|n}\mu(d)=\epsilon(n)$
所以式子变成： $$\sum_{1≤i≤n}\sum_{1≤j≤m}\sum_{d|(i,j)}\mu(d)$$
即 $$\sum_{1≤i≤n}\sum_{1≤j≤m}\sum_{d|i , d|j}\mu(d)$$
于是套路变形： $$\sum_{1≤d≤min(n,m)}\mu(d)\times \lfloor\frac{n}{d}\rfloor\times\lfloor\frac{m}{d}\rfloor$$

然后就可以做了。
还是套路..对$\lfloor\frac{i}{d}\rfloor\times\lfloor\frac{j}{d}\rfloor$相等的进行分块，每次直接计算这段区间的和。但是因为$d$是会变的，所以计算完$mobius$之后要做个前缀和。

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
#define maxn 50100

bool ispri[maxn];
int k,cnt,pri[maxn],sum[maxn],mu[maxn];
int mymin(int x,int y){return (x<y)?x:y;}
void mobius(int lim)
{
    cnt=0;sum[0]=0;mu[1]=1;
    for (int i=2;i<=lim;i++)
    {
        if (!ispri[i])
        {
            pri[++cnt]=i;
            mu[i]=-1;
        }
        for (int j=1;j<=cnt && pri[j]*i<=lim;j++)
        {
            ispri[pri[j]*i]=true;
            if (i%pri[j]==0)
            {
                mu[i*pri[j]]=0;
                break;
            }
            mu[i*pri[j]]=-mu[i];
        }
    }
    for (int i=1;i<=50000;i++) sum[i]=sum[i-1]+mu[i];
}
int get(int n,int m)
{
    int lim,ret,x,y,lx,ly,r;
    ret=0;n/=k;m/=k;
    lim=mymin(n,m);
    for (int i=1;i<=lim;i=r+1)
    {
        x=n/i;y=m/i;
        lx=mymin(lim,n/x);
        ly=mymin(lim,m/y);
        r=mymin(lx,ly);
        ret+=(sum[r]-sum[i-1])*x*y;
    }
    return ret;
}
int main()
{
    //freopen("a.in","r",stdin);
    //freopen("a.out","w",stdout);
    int T,a,b,c,d;
    scanf("%d",&T);
    mobius(50000);
    while (T--)
    {
        scanf("%d%d%d%d%d",&a,&b,&c,&d,&k);
        int ans=get(b,d)-get(a-1,d)-get(b,c-1)+get(a-1,c-1);
        printf("%d\n",ans);
    }
    return 0;
}