BZOJ 2301 Problem b (莫比乌斯反演+区间容斥+分块)

思路:

这个题算是bzoj 1695的加强版吧,同样是给你a,b,c,d,k五个数,只不过这次a和c都不一定是1了。那么我们需要进行区间的容斥定理(注意端点值)。同时由于数据量很大,我们还需要分块。

容斥:

fanyan(b,d)fanyan(a1,d)fanyan(b,c1)+fanyan(a1,c1)

分块方法:
对于区间(1,a)和(1,b),求gcd(i,j)==k的数对的个数。

int minab = min(a,b),l;
a/=k,b/=k;
forint i = 1;i <= minab;i=l+1){
    l = min(a/(a/i),b/(b/i));//注意,这个方法必须是在i连续的情况下,所以才有了a/=k,b/=k。
    ans += (sum[l]-sum[i-1]) * (a/i) * (b/i);//sum[i] 为μ(i)的前缀和。
}

所以本题也就迎刃而解了。

#include<stdio.h>
#include <iostream>
#include<string.h>
#include<math.h>
#include<algorithm>
#define eps 1e-8
typedef long long int lli;
using namespace std;

const int maxn = 5e4+20;
bool isprime[maxn];
//int phi[maxn];
int prime[maxn],miu[maxn];
//int cnt[maxn];质因数的个数 8:2*2*2 3个
void moblus(){
    int num = 0;miu[1] = 1;
    for(lli i = 2;i < maxn;i++){
        if(!isprime[i]){
            prime[num++] = i,miu[i] = -1;//phi[i] = i-1;
            //cnt[i] = 1;
        }
        for(lli j = 0;j < num && i*prime[j] < maxn;j++){
            lli x = prime[j];
            isprime[i*x] = 1;//cnt[i*x] = cnt[i]+1;
            if(i%x){
                miu[i*x] = -miu[i];
                //phi[i*x] = phi[i] * phi[x];
            }
            else{
                miu[i*x] = 0;
                //phi[i*x] = phi[i] * x;
                break;
            }
        }
    }
}
int sum[maxn];
lli fy(lli n,lli m,lli k){
    lli ans = 0,l;
    if(n>m) swap(n,m);
    n/=k;m/=k;
    for(lli i = 1;i <= n;i=l+1){
        l = min(n/(n/i),m/(m/i));
        ans += (lli)(sum[l]-sum[i-1])*(n/i)*(m/i);
    }
    return ans;
}

int main(){
    moblus();
    for(int i = 1;i <= maxn;i++){
        sum[i] = sum[i-1]+miu[i];
    }
    int q,n,a,b,c,d,k;
    scanf("%d",&q);
    while(q--){
        scanf("%d%d%d%d%d",&a,&b,&c,&d,&k);
        lli ans = 0;
        ans = fy(b,d,k) + fy(a-1,c-1,k) - fy(b,c-1,k) - fy(a-1,d,k);
        printf("%lld\n",ans);
    }

    return 0;
}
发布了203 篇原创文章 · 获赞 69 · 访问量 5万+
展开阅读全文

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

©️2019 CSDN 皮肤主题: 大白 设计师: CSDN官方博客

分享到微信朋友圈

×

扫一扫,手机浏览