筛法求欧拉函数之和

题目

#include<iostream>
using namespace std;
const   int N=1000010;
bool st[N];
int prime[N];
int ph[N];
int cnt;
long long ola(int n)
{
    ph[1]=1;
    for(int i=2;i<=n;i++)
    {
        if(!st[i])
        {
            prime[cnt++]=i;
            ph[i]=i-1;
        }
        for(int j=0;prime[j]<=n/i;j++)
        {
            st[i*prime[j]]=1;
            if(i%prime[j]==0)
            {
                ph[i*prime[j]]=ph[i]*prime[j];
                break;
            }
            else
            {
                ph[i*prime[j]]=ph[i]*(prime[j]-1);
            }
        }
    }  
    long long ans=0;
        for(int i=1;i<=n;i++)
        ans+=ph[i];
        return ans;
}
int main()
{
    int n;
    cin>>n;
    long long t=ola(n);
    cout<<t<<endl;
    return 0;
}