GCD - Extreme (II) UVA - 11426 (欧拉函数）

题意：
G=0;
for(i=1;i<N;i++)
for(j=i+1;j<=N;j++)
{
G+=gcd(i,j);
}
给出n，求出想应的G
思路：
进行打表。G[i]表示要求的值，a[i]表示从1–(i-1)与i的GCD之和。可以推出G[i]=G[i-1]+a[i]。
所以我们的重点在于求出a[i]。这里就需要对欧拉函数进行利用了，很巧妙。
依次暴力来求肯定超时，，i与（1–i-1）肯定存在相同的最大公因子，a[i]可以用i与1–(i-1)的最大公因子与公因子个数乘积之和来表示，gcd(i,j)==x那么gcd(i/x,j/x)==1,所以与i的最大公因子为x的个数为phi(i/x);
当然求a[i]时，不能从1-n依次求，这样需要判断因子，可以让因子从1-n依次累加到a[i]上。

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<set>
#include<map>
#include<vector>
#include<queue>
using namespace std;
#define maxn 4000009
long long a[maxn+4],dp[maxn+4];
int phi[maxn+4];
void phi_()
{
    for(int i=1;i<=maxn;i++)
        phi[i]=i;
    for(int i=2;i<=maxn;i+=2)
        phi[i]/=2;
    for(int i=3;i<=maxn;i+=2)
    {
        if(phi[i]==i)
        {
            for(int j=i;j<=maxn;j+=i)
                phi[j]=phi[j]-phi[j]/i;
        }
    }
}
void init()
{
    memset(a,0,sizeof(a));
    for(int i=1;i<=maxn;i++)
    {
        for(int j=i+i;j<=maxn;j+=i)
        {
            a[j]+=i*phi[j/i];
        }
    }
    dp[1]=0;
    for(int i=2;i<=maxn;i++)
    {
        dp[i]=dp[i-1]+a[i];
    }
}
int  main()
{
    int n;
    phi_();
    init();
   while(scanf("%d",&n)!=EOF)
         {
             if(n==0)
             break;
             printf("%lld\n",dp[n]);
         }
    return 0;
}