[Luogu2158] 儀仗隊 線性篩歐拉函數

Luogu2158

#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;
 }