gym101982 B.Coprime Integers（莫比乌斯反演，容斥）

题意：

$给定a,b,c,d,计算{\sum }_{x=a}^b{\sum }_{y=c}^d[gcd(x,y)=1]$

数据范围：a,b,c,d<=1e7

解法：

$给定a,b,c,d,计算{\sum }_{x=a}^b{\sum }_{y=c}^d[gcd(x,y)=1]$

$写过一道更难的题，是计算{\sum }_{x=a}^b{\sum }_{y=c}^d[gcd(x,y)=k]$

$令cal(a,b)={\sum }_{i=1}^a{\sum }_{j=1}^b[gcd(i,j)=k]$

$根据容斥,答案为cal(b,d)-cal(a-1,d)-cal(b,c-1)+cal(a-1,c-1)$

$那么问题变为如何计算cal(a,b)$

$---$

$设f(k)={\sum }_{i=1}^a{\sum }_{j=1}^b[gcd(i,j)=k]$

$f(k)表示gcd为k的方案数$

$F(n)={\sum }_{n|k}f(k)=\lfloor \frac{a}{n}\rfloor \lfloor \frac{b}{n}\rfloor$

$F(n)表示gcd为n的倍数的方案数$

$反演得:f(k)={\sum }_{n|k}\mu (\frac{k}{n})F(k)={\sum }_{n|k}\mu (\frac{k}{n})\lfloor \frac{a}{k}\rfloor \lfloor \frac{b}{k}\rfloor$

$到这一步发现可以O(n)计算了,我就不继续优化了(可以优化到O(sq))$

code：

#include<bits/stdc++.h>
using namespace std;
#define ll long long
const int maxm=1e7+5;
int np[maxm];
int p[maxm],cnt;
int mu[maxm];
int a,b,c,d;
void init(){
    mu[1]=1;
    for(int i=2;i<maxm;i++){
        if(!np[i]){
            p[cnt++]=i;
            mu[i]=-1;
        }
        for(int j=0;j<cnt;j++){
            if(1ll*p[j]*i>=maxm)break;
            np[p[j]*i]=1;
            mu[p[j]*i]=i%p[j]?-mu[i]:0;
            if(i%p[j]==0)break;
        }
    }
}
ll cal(int a,int b){
    ll ans=0;
    for(int k=1;k<=min(a,b);k++){
        ans+=1ll*mu[k]*(a/k)*(b/k);
    }
    return ans;
}
signed main(){
    init();
    scanf("%d%d%d%d",&a,&b,&c,&d);
    ll ans=cal(b,d)-cal(a-1,d)-cal(b,c-1)+cal(a-1,c-1);
    printf("%lld\n",ans);
    return 0;
}

---