AcWing 874. 筛法求欧拉函数
原题链接
简单
作者:
术
,
2021-03-01 16:09:42
,
所有人可见
,
阅读 113
#include <iostream>
using namespace std;
const int N=1000005;
long long phi[N];
int primes[N];
int cnt;
bool st[N];
long long get_euler(int n){
phi[1]=1;
for(int i=2;i<=n;i++){
if(!st[i]){
primes[cnt++]=i;
phi[i]=i-1;
}
for(int j=0;primes[j]<=n/i;j++){
st[primes[j]*i]=true;
if(i%primes[j]==0){
phi[primes[j]*i]=phi[i]*primes[j];
break;
}
phi[primes[j]*i]=phi[i]*(primes[j]-1);
}
}
long long res=0;
for(int i=1;i<=n;i++)
res+=phi[i];
return res;
}
int main()
{
int n;
cin>>n;
cout<<get_euler(n);
//cout << "Hello world!" << endl;
return 0;
}