SPOJ 3871 GCD Extreme

题目链接：http://www.spoj.com/problems/GCDEX/

题意：

思路：

 

int phi[N];

i64 f[N];



void init()

{

    phi[1]=1;

    int i,j;

    for(i=2;i<N;i++) if(!phi[i]) for(j=i;j<N;j+=i)

    {

        if(!phi[j]) phi[j]=j;

        phi[j]-=phi[j]/i;

    }

    for(i=2;i<N;i++) f[i]=phi[i];

    for(i=2;i*i<N;i++)

    {

        f[i*i]+=(i64)phi[i]*i;

        for(j=i+1;j*i<N;j++) f[i*j]+=(i64)phi[i]*j+(i64)phi[j]*i;

    }

    for(i=2;i<N;i++) f[i]+=f[i-1];

}



int n;



int main()

{

    init();

    Rush(n)

    {

        if(n==0) break;

        printf("%lld\n",f[n]);

    }

}