似乎只要会了莫比乌斯反演就不用会phi了、
感觉反演更顺手、
码:
#include<iostream>
#include<cstdio>
using namespace std;
#define N 40005
int mu[N],i,j,n,tot,ans,su[N];
bool he[N];
void eular()
{
mu[1]=1;
for(i=2;i<=n;i++)
{
if(!he[i])
{
su[++tot]=i;
mu[i]=-1;
}
for(j=1;j<=tot&&su[j]*i<=n;j++)
{
he[su[j]*i]=1;
if(i%su[j]==0)
{
mu[i*su[j]]=0;
break;
}else mu[i*su[j]]=-mu[i];
}
}
}
int main()
{
scanf("%d",&n);
if(n==1)
{
printf("0");
return 0;
}
n--;
eular();
for(i=1;i<=n;i++)
ans+=mu[i]*(n/i)*(n/i);
//cout<<ans;
ans+=2;
printf("%d",ans);
}