#include<iostream>
#include<iomanip>
#include<algorithm>
#include<cmath>
#include<cstring>
#include<cstdio>
#define maxn 40010
#define ll long long
using namespace std;
int prime[maxn],primenum,phi[maxn],n,ans;
bool vis[maxn];
void pri(int v)
{
int cnt=0;
memset(vis,0,sizeof(vis));
for(int i=2;i<=v;i++)
{
if(!vis[i])
prime[++primenum]=i,phi[i]=i-1;
for(int j=1;j<=primenum&&i*prime[j]<=v;++j)
{
vis[i*prime[j]]=1;
if(i%prime[j]==0)
{
phi[i*prime[j]]=phi[i]*prime[j];
break;
}
else phi[i*prime[j]]=phi[i]*(prime[j]-1);
}
}
}
int main()
{
scanf("%d",&n);
if(n==1)
{
cout<<0;
return 0;
}
pri(n);
for(int i=2;i<=n-1;++i)
ans+=phi[i];
cout<<ans*2+3<<endl;
return 0;
}
[Luogu2158] 仪仗队 线性筛欧拉函数
最新推荐文章于 2023-04-28 18:06:13 发布