POJ 2478 Farey Sequence（欧拉函数前n项和）

#include <iostream>
#include <math.h>
using namespace std;
long long a[1000005]={0};
long long c[1000005]={0};
void enlur()
{
    int i,j;
    for(i=2;i<1000005;i++)
    {
        if(!a[i])
        {
            for(j=i;j<1000005;j=j+i)
            {
                if(!a[j])
                a[j]=j;
                a[j]=a[j]/i*(i-1);
            }
        }
    }
}
int main()
{
    int i,n;
    enlur();
    c[0]=a[0];
    for(i=1;i<1000005;i++)
    c[i]=c[i-1]+a[i];
    while(cin>>n&&n)
    cout<<c[n]<<endl;
    return 0;
}