ABC206 E - Divide Both(容斥)

Posted live4m

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了ABC206 E - Divide Both(容斥)相关的知识,希望对你有一定的参考价值。

题意:

在这里插入图片描述

解法:

考虑枚举gcd=i,计算gcd=i的数对数量(先不考虑g!=x和g!=y的这个条件),
然后减去不合法条件(即g=x或者g=y的情况),

令d[i]为gcd=i的数对数量,可以容斥计算:
d[i]=[r/i-(l-1)/i]^2-d[i*2]-d[i*3]....

先ans+=sum{d[2,r]}.

由于需要减去g=x或者g=y的情况,
枚举g=[max(2,l),r],
那么ans-={数对(g,kg)的数量}+{数对(kg,g)的数量}+{(g,g)的数量},
其中k>=1,设cnt=r/i-(l-1)/i,容易得出上式可以转化为ans-=cnt*2+1.

code:

#include<bits/stdc++.h>
#define int long long
using namespace std;
const int maxm=4e6+5;
int d[maxm];
void solve(){
    int l,r;cin>>l>>r;
    for(int i=r;i>=1;i--){
        int cnt=r/i-(l-1)/i;
        d[i]=cnt*cnt;
        for(int j=i+i;j<=r;j+=i){
            d[i]-=d[j];
        }
    }
    int ans=0;
    for(int i=2;i<=r;i++){//枚举g
        ans+=d[i];
        if(i>=l&&i<=r){
            int cnt=r/i-(l-1)/i;
            ans-=cnt*2-1;
        }
    }
    cout<<ans<<endl;
}
signed main(){
    ios::sync_with_stdio(0);cin.tie(0);
    solve();
    return 0;
}

以上是关于ABC206 E - Divide Both(容斥)的主要内容,如果未能解决你的问题,请参考以下文章

ABC206 C - Swappable(map计数)

ABC172

Divide by Zero 2021 and Codeforces Round #714 (Div. 2) ABC题解

ABC206 D - KAIBUNsyo(思维,并查集求连通块大小)

异常访问HDFS权限异常:Permission denied: user=СÎ÷?Îè, access=WRITE, inode="/"

Divide by Zero 2017 D&E&F