Gym - 101982B Coprime Integers (莫比乌斯反演)

Posted zjl192628928

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了Gym - 101982B Coprime Integers (莫比乌斯反演)相关的知识,希望对你有一定的参考价值。

题目链接:http://codeforces.com/gym/101982/attachments

题目大意:有区间[a,b]和区间[c,d],求gcd(x,y)=1,其中x属于[a,b],y属于[c,d],求这样的x,y有多少对。

解题思路:

第一种反演思路:

把下界变换一下

技术图片

代码:

 

#include<iostream>
#include<cstdio>
using namespace std;
typedef long long ll;
const int maxn=1e7+7;
ll a,b,c,d,mu[maxn],sum[maxn],prime[maxn],tot;
void getMobius(int N){
    for(int i=1;i<=N;i++) prime[i]=1;
    mu[1]=1;
    tot=0;
    for(int i=2;i<=N;i++){
        if(prime[i]){
            prime[tot++]=i;
            mu[i]=-1;
        }
        for(int j=0;j<tot&&prime[j]*i<=N;j++){
            prime[prime[j]*i]=0;
            if(i%prime[j]==0){
                mu[i*prime[j]]=0;
                break;
            }
            mu[i*prime[j]]=-mu[i];
        }
    }
}
int main()
{
    scanf("%lld%lld%lld%lld",&a,&b,&c,&d);
    getMobius(1e7);
    ll ans=0;
    for(int i=1;i<=min(b,d);i++)
        ans+=mu[i]*(b/i-(a-1)/i)*(d/i-(c-1)/i);
    printf("%lld\\n",ans);
    return 0;
}

 

第二种反演思路:

技术图片

右边全部都是已知的,枚举下可取范围内的d(也就是原来n的倍数,这里n是1)

可以利用容斥原理,先求出[1,b]和[1,d],再减去[1,a-1]和[1,d]以及[1,b]和[1,c-1],最后加上多减的部分[1,a-1]和[1,c-1]。

并且很显然,推演最后得到的式子是可以经过整除分块优化的,只需要预处理出莫比乌斯函数的前缀和即可。

代码:

#include<iostream>
#include<cstdio>
using namespace std;
typedef long long ll;
const int maxn=1e7+7;
ll a,b,c,d,mu[maxn],sum[maxn],prime[maxn],tot;
void getMobius(int N){
    for(int i=1;i<=N;i++) prime[i]=1;
    mu[1]=1;
    tot=0;
    for(int i=2;i<=N;i++){
        if(prime[i]){
            prime[tot++]=i;
            mu[i]=-1;
        }
        for(int j=0;j<tot&&prime[j]*i<=N;j++){
            prime[prime[j]*i]=0;
            if(i%prime[j]==0){
                mu[i*prime[j]]=0;
                break;
            }
            mu[i*prime[j]]=-mu[i];
        }
    }
}
ll solve(ll n,ll m){
    ll res=0;
    for(ll l=1,r=0;l<=min(n,m);l=r+1){
        r=min(n/(n/l),m/(m/l));
        res+=(sum[r]-sum[l-1])*(n/l)*(m/l);
    }
    return res;
}
int main()
{
    scanf("%lld%lld%lld%lld",&a,&b,&c,&d);
    getMobius(1e7);
    sum[1]=1;
    for(int i=2;i<=1e7;i++) sum[i]=sum[i-1]+mu[i];
    printf("%lld\\n",solve(b,d)-solve(a-1,d)-solve(c-1,b)+solve(a-1,c-1));
    return 0;
}

 

以上是关于Gym - 101982B Coprime Integers (莫比乌斯反演)的主要内容,如果未能解决你的问题,请参考以下文章

最小圆覆盖 gym-102006 I

HDU 6025 Coprime Sequence

A - Anticlockwise Motion Gym - 101177A

JZOJ1227 Coprime 题解

HDU - 6025 Coprime Sequence(gcd+前缀后缀)

数据结构_coprime_sequence(互质序列)