#include<bits/stdc++.h>
#define LL long long
using namespace std;
const int maxn=1e7+;
bool vis[maxn];
int prime[maxn];
int mu[maxn];
int sum1[maxn];
int sum2[maxn];
int tot=;
void get_mu()// mo bi su si han shu
{
mu[]=; vis[]=;
for(int i=;i<maxn;i++) // prime = 0; other = 1;
{
if(!vis[i]){ prime[++tot]=i; mu[i]=-;}
for(int j=;j<=tot&& prime[j]*i<maxn;j++)
{
vis[prime[j]*i]=;
if(i%prime[j]==)break;
mu[i*prime[j]]=-mu[i];
}
}
for(int i=;i<=tot;i++)
{
for(int j=prime[i];j<maxn;j+=prime[i])
{
sum1[j]+=mu[j/prime[i]];
}
}
//for(int i=1;i<maxn;i++) sum2[i]=sum2[i-1]+sum1[i];
}
int main()
{
get_mu();
int n; cin>>n;
LL ans=;
for(int i=;i<=n;i++)
{
ans+=1LL*(n/i)*(n/i)*sum1[i];
}
cout<<ans<<endl;
}
过度代码
整除分块 (看起来更麻烦)
#include<bits/stdc++.h>
#define LL long long
using namespace std;
const int maxn=1e7+;
bool vis[maxn];
int prime[maxn];
int mu[maxn];
int sum1[maxn];
int sum2[maxn];
int tot=;
void get_mu()// mo bi su si han shu
{
mu[]=; vis[]=;
for(int i=;i<maxn;i++) // prime = 0; other = 1;
{
if(!vis[i]){ prime[++tot]=i; mu[i]=-;}
for(int j=;j<=tot&& prime[j]*i<maxn;j++)
{
vis[prime[j]*i]=;
if(i%prime[j]==)break;
mu[i*prime[j]]=-mu[i];
}
}
for(int i=;i<=tot;i++)
{
for(int j=prime[i];j<maxn;j+=prime[i])
{
sum1[j]+=mu[j/prime[i]];
}
}
for(int i=;i<maxn;i++) sum2[i]=sum2[i-]+sum1[i];
}
int main()
{
get_mu();
int n; cin>>n;
LL ans=;
//for(int i=1;i<=n;i++) ans+=1LL*(n/i)*(n/i)*sum1[i];
for(int l=,r;l<=n;l=r+)
{
r=n/(n/l); // l-r 区间相同值 区间值n/l
ans+=1LL*(n/l)*(n/l)*(sum2[r]-sum2[l-]);
}
cout<<ans<<endl;
}