AmazingCounters.com
随笔- 138  评论- 7  文章- 0 

[BZOJ]2301: [HAOI2011]Problem b

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

思路:

我们只要能求出1<=x<=n,1<=y<=m时的数对个数,容斥一下即可,求法如下:

复制黄学长复制的鏼爷的题解

推导:




用莫比乌斯函数的性质把求和的式子换掉,

其中,更换求和指标,

容易知道单调不上升,且最多有种不同的取值。所以按取值分成个段分别处理,一个连续段内的和可以用预处理出的莫比乌斯函数前缀和求出

#include<cstdio>
#include<algorithm>
using namespace std;
inline int read()
{
    int x;char c;
    while((c=getchar())<'0'||c>'9');
    for(x=c-'0';(c=getchar())>='0'&&c<='9';)x=(x<<3)+(x<<1)+c-'0';
    return x;
}
#define MN 50000
int u[MN+5],f[MN+5],p[MN+5],pn;
int cal(int n,int m)
{
    if(n>m)swap(n,m);
    int res=0,ls,i;
    for(i=1;i<=n;i=ls+1)
    {
        ls=min(n/(n/i),m/(m/i));
        res+=(u[ls]-u[i-1])*(n/i)*(m/i);
    }
    return res;
}
int main()
{
    int n=read(),i,j,a,b,c,d,k;
    for(u[1]=1,i=2;i<=MN;++i)
    {
        if(!f[i])p[++pn]=i,u[i]=-1;
        for(j=1;i*p[j]<=MN&&(f[i*p[j]]=1);++j)
            if(i%p[j])u[i*p[j]]=-u[i];else break;
        u[i]+=u[i-1];
    }
    while(n--)a=read()-1,b=read(),c=read()-1,d=read(),k=read(),
        printf("%d\n",cal(b/k,d/k)-cal(a/k,d/k)-cal(b/k,c/k)+cal(a/k,c/k));
}

 

posted on 2017-04-01 14:58  ditoly  阅读(...)  评论(...编辑  收藏