欧拉函数2:最大公约数
作者:
总打瞌睡的天天啊
,
2024-08-06 22:46:16
,
所有人可见
,
阅读 1
//看思路
#include<iostream>
#include<algorithm>
#include<cstring>
using namespace std;
typedef long long LL;
const int N=1e7+10;
int primes[N],cnt;
bool st[N];
LL phi[N],s[N];
void init(int n)
{
for(int i=2;i<=n;i++)
{
if(!st[i])
{
primes[cnt++]=i;
phi[i]=i-1;
}
for(int j=0;primes[j]*i<=n;j++)
{
st[primes[j]*i]=true;
if(i%primes[j]==0)
{
phi[i*primes[j]]=primes[j]*phi[i];
break;
}
phi[i*primes[j]]=phi[i]*(primes[j]-1);
}
}
for(int i=1;i<=n;i++)s[i]=s[i-1]+phi[i];
}
int main()
{
int n;
cin>>n;
init(n);
LL res=0;
for(int i=0;i<cnt;i++)
{
int p=primes[i];
res+=s[n/p]*2+1;
}
cout<<res<<endl;
return 0;
}