BZOJ 2301 [HAOI2011]Problem b - 莫比乌斯反演+容斥

x_1023

于 2017-08-17 10:50:11 发布

阅读量255

点赞数

分类专栏：莫比乌斯反演数论

本文链接：https://blog.csdn.net/x_1023/article/details/77314751

版权

莫比乌斯反演同时被 2 个专栏收录

2 篇文章

订阅专栏

数论

2 篇文章

订阅专栏

本文介绍了一种使用线筛法快速求解莫比乌斯函数的方法，并通过一个具体的C++实现示例展示了如何计算区间内的莫比乌斯函数值，进而解决涉及数论中莫比乌斯反演的问题。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

这里写图片描述

#include<iostream>
#include<cstring>
#include<cstdio>
#include<cstdlib>
#include<algorithm>

using namespace std;

const int maxn=100005;

int n,cnt,k;
bool is[maxn];
long long pre[maxn];
long long prime[maxn];
long long miu[maxn];

void linear_shaker()
{
    miu[1]=1;
    int t;
    for(int i=2;i<maxn;i++)
    {
        if(!is[i])
        {
            miu[i]=-1;
            prime[++cnt]=i;
        }
        for(int j=1;j<=cnt&&(t=prime[j]*i)<maxn;j++)
        {
            is[t]=true;
            if(i%prime[j]==0)
            {
                miu[t]=0;
                break;
            }
            miu[t]=-miu[i];
        }
    }
    for(int i=1;i<maxn;i++)
        pre[i]=pre[i-1]+miu[i];
}
long long cal(int a,int b)
{
    long long res=0,last;
    a/=k;b/=k;
    for(int i=1;i<=min(a,b);i=last+1)
    {
        last=min(b/(b/i),a/(a/i));
        res+=(pre[last]-pre[i-1])*(a/last)*(b/last);
    }
    return res;
}
int main()
{
    linear_shaker();
    scanf("%d",&n);
    while(n--)
    {
        int a,b,c,d;
        scanf("%d%d%d%d%d",&a,&b,&c,&d,&k);
        printf("%lld\n",cal(b,d)-cal(a-1,d)-cal(b,c-1)+cal(a-1,c-1));
    }
}