874. 筛法求欧拉函数
题目链接https://www.acwing.com/problem/content/876/
题目:
欧拉函数
思路:线性筛法的基础上修改一些就可以得出答案
图片来着该佬佬的截图https://www.acwing.com/solution/content/3952/
#include<iostream>
#include<cstdio>
using namespace std;
const int N=1e6+10;
bool vis[N];
int primes[N];
int phi[N];
int cnt=0;
void get(int u){
phi[1]=1;
for(int i=2;i<=u;i++){
if(!vis[i]){
primes[cnt++]=i;
phi[i]=i-1;
}
for(int j=0;primes[j]<=u/i;j++){
vis[primes[j]*i]=1;
if(i%primes[j]==0){
phi[i*primes[j]]=phi[i]*primes[j];
break;
}
phi[i*primes[j]]=phi[i]*(primes[j]-1);
}
}
}
int main(){
int n;
cin>>n;
get(n);
long long ans=0;
for(int i=1;i<=n;i++){
ans+=phi[i];
}
printf("%lld",ans);
return 0;
}