#include <iostream>
using namespace std;
const int N=1000010;
typedef long long LL;
int primes[N],cnt,phi[N];
bool str[N];
int get_eulers(int x)
{
phi[1]=1;
LL res=0;
for(int i=2;i<=x;i++)
{
if(!str[i])
{
primes[cnt++]=i;
phi[i]=i-1;
}
for(int j=0;primes[j]<=x/i;j++)
{
str[i*primes[j]]=true;
if(i%primes[j]==0)
{
phi[i*primes[j]]=phi[i]*primes[j];
break;
}
phi[i*primes[j]]=phi[i]*(primes[j]-1);
}
}
for(int i=1;i<=x;i++)
{
res+=phi[i];
}
return res;
}
int main()
{
int n;
cin>>n;
cout<<get_eulers(n)<<endl;
return 0;
}
筛质法求欧拉函数
最新推荐文章于 2025-08-05 20:04:56 发布