BZOJ2301[HAOI2011]problem b

Description

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

Input

第一行一个整数n,接下来n行每行五个整数,分别表示a、b、c、d、k

Output

共n行,每行一个整数表示满足要求的数对(x,y)的个数

Sample Input
2

2 5 1 5 1

1 5 1 5 2

Sample Output

14

3

HINT

100%的数据满足:1≤n≤50000,1≤a≤b≤50000,1≤c≤d≤50000,1≤k≤50000


算是写的第一道反演题。题目要求的是一段区间,通常的做法是求出1到右边界,然后再用容斥来求解区间。
我们考虑如何求1到右边界,设右边界为nm,则要求的式子是i=1nj=1m[gcd(i,j)=k],我们来推下式子。
i=1nj=1m[gcd(i,j)=k]=i=1nkj=1mk[gcd(i,j)=1]=i=1nkj=1mkx|gcd(i,j)μ(x)(根据性质1)
接下来要做最重要的一步转换,大家要好好理解:i=1nkj=1mkx|gcd(i,j)μ(x)=x=1min(n,m)kμ(x)[nkx][mkx]
这步大家可能不太理解,我刚看的时候也是一脸懵逼,简单的说就是将x提到外面,而右边[nkx][mkx]表示的是x的倍数的个数。
我们发现[nkx][mkx]可以用数论分块,所以直接分块枚举x即可。

这道题还有另外一种推法,需要构造。
我们令f(x)=i=1nj=1m[gcd(i,j)=x],F(x)=i=1nj=1m[x|gcd(i,j)]
那么显然F(x)=x|df(d),那么根据莫比乌斯反演f(x)=x|dμ(dx)F(d)
我们考虑F(x)是什么,F(x)必须要求i,j都是x的倍数,所以F(x)=[nx][mx]
所以答案就是f(k)=k|xμ(xk)[nx][mx]=x=1min(n,m)kμ(x)[nkx][mkx]
与上面的做法殊途同归。
#include<bits/stdc++.h>
using namespace std;
int read(){
    char c;int x;while(c=getchar(),c<'0'||c>'9');x=c-'0';
    while(c=getchar(),c>='0'&&c<='9') x=x*10+c-'0';return x;
}
void print(int x){
    if(x/10) print(x/10);
    putchar(x%10+'0');
}
int T,a,b,c,d,k,top,ans,mu[50005],sum[50005],pri[50005],vis[50005];
int calc(int n,int m){
    if(n>m) swap(n,m);
    int l=1,r=0,res=0;n=n/k;m=m/k;
    while(l<=n){
        r=min(n/(n/l),m/(m/l));r=min(r,n);
        res+=(n/l)*(m/l)*(sum[r]-sum[l-1]);
        l=r+1;
    }
    return res;
}
int main()
{
    T=read();mu[1]=sum[1]=1;
    for(int i=2;i<=50000;i++){
        if(!vis[i]) pri[++top]=i,mu[i]=-1;
        for(int j=1;j<=top&&pri[j]*i<=50000;j++){
            vis[i*pri[j]]=1;
            if(i%pri[j]==0) break;
            mu[i*pri[j]]=-mu[i];
        }
        sum[i]=sum[i-1]+mu[i];
    }
    while(T--){
        a=read();b=read();c=read();d=read();k=read();
        ans=calc(b,d)-calc(a-1,d)-calc(b,c-1)+calc(a-1,c-1);
        print(ans);puts("");
    }
    return 0;
}
发布了184 篇原创文章 · 获赞 40 · 访问量 5万+
展开阅读全文

没有更多推荐了,返回首页

©️2019 CSDN 皮肤主题: 创作都市 设计师: CSDN官方博客

分享到微信朋友圈

×

扫一扫,手机浏览