洛谷P2568 GCD（莫比乌斯反演）

 

这题和一样……不过是n和m相同而已……

所以虽然正解是欧拉函数然而直接改改就行了所以懒得再码一遍了2333

不过这题卡空间，记得mu开short，vis开bool

//minamoto
#include<cstdio>
#define ll long long
const int N=1e7+5;
int p[1000005],n,m;short mu[N];bool vis[N];ll sum[N],ans;
void init(int n){
    mu[1]=1;
    for(int i=2;i<=n;++i){
        if(!vis[i]) mu[i]=-1,p[++m]=i;
        for(int j=1;j<=m&&p[j]*i<=n;++j){
            vis[i*p[j]]=1;
            if(i%p[j]==0) break;
            mu[i*p[j]]=-mu[i];
        }
    }
    for(int j=1;j<=m;++j)
    for(int i=1;i*p[j]<=n;++i)
    sum[i*p[j]]+=mu[i];
    for(int i=1;i<=n;++i) sum[i]+=sum[i-1];
}
int main(){
    scanf("%d",&n),init(n);
    for(int l=1,r;l<=n;l=r+1){
        r=n/(n/l);
        ans+=(sum[r]-sum[l-1])*(n/l)*(n/l);
    }
    printf("%lld\n",ans);
    return 0;
}