有两种做法一个是欧拉函数的前缀和另一个是莫比乌斯反演
貌似数据水能让莫比乌斯反演过去
我用的是莫比乌斯反演
#include<cstdio>
#include<iostream>
#include<cstring>
using namespace std;
const int maxn=40001;
int mu[maxn],pre[maxn],prime[maxn];
bool check[maxn];
int tot;
int main()
{
int i,j,k;
mu[1]=1;
for(i=2;i<maxn;i++)
{
if(!check[i])
prime[++tot]=i,mu[i]=-1;
for(j=1;j<=tot;j++)
{
k=prime[j]*i;
if(k>=maxn)break;
if(i%prime[j]==0){check[k]=true,mu[k]=0;break;}
mu[k]=-mu[i];check[k]=true;
}
}
for(i=1;i<maxn;i++)pre[i]=pre[i-1]+mu[i];
int n;
scanf("%d",&n);
n--;
int T=1,Div1,Div2,l;
int ans=0;
while(T<=n)
{
Div1=n/T;
l=n/Div1+1;
ans+=Div1*Div1*(pre[l-1]-pre[T-1]);
T=l;
}
ans+=2;
printf("%d\n",ans);
return 0;
}