UVA11426 欧拉函数

大白书P125

 

#include <iostream>
#include <cstring>
using namespace std;
#define MMX 4000010
#define LL long long
int phi[MMX],f[MMX];
LL S[MMX];

void calc_phi(int n)        //求1--n的欧拉函数，phi[i]=φ(i)
{
    for (int i=2;i<=n;i++)
        phi[i]=0;
    phi[1]=1;
    for (int i=2;i<=n;i++)
        if (!phi[i])
            for (int j=i;j<=n;j+=i)
            {
                if (!phi[j])    phi[j]=j;
                phi[j]=phi[j]/i*(i-1);
            }
}

int main()
{
    calc_phi(MMX);
    memset(f,0,sizeof(f));
    for (int i=1;i<=MMX;i++)
        for (int j=2*i;j<=MMX;j+=i)
            f[j]+=i*phi[j/i];

    S[2]=f[2];
    for (int i=3;i<=MMX;i++)
        S[i]=S[i-1]+f[i];

    int n;
    while (cin>>n)
    {
        if (n==0) break;
        cout<<S[n]<<endl;
    }
}