SPOJ 3871 GCD Extreme

//数学题，对于a,b满足gcd(a,b)=1,则ka和kb的最大公约数为k，所以其实每对数的最大公约数都是由两个互质的数分别乘以相同倍数得到的，所以只要找出1到n中互质的数对就好做了。这里可以用欧拉函数求出小于 i 的并且与 i 互质的数的个数，那么形如(kx,ki)(x<i)的最大公约数

之和就是k*phi(i)了。

#include<iostream>
#include<cstring>
#include<cmath>
#include<cstdio>
#define N 1000001
#define ll long long
using namespace std;
int p[N];
ll a[N];
void solve()
{
         for(int i=2;i<=N;i++)p[i]=i;
         for(ll i=2;i<=N;i++)
         {
                 if(p[i]==i)
                 {
                            for(int j=i;j<=N;j+=i)
                            p[j]=p[j]/i*(i-1);
                 }
         }
         for(int i=2;i<=N;i++)a[i]=p[i];
         int t=sqrt(N);
         for(int i=2;i<=t;i++)
         {
            a[i*i]+=p[i]*i;
            for(int j=i*i+i;j<=N;j+=i)
           {
             a[j]+=p[i]*j/i+p[j/i]*i;
           }
        }
         for(int i=1;i<=N;i++)a[i]+=a[i-1];
}
int main()
{
   solve();
   int n;
   while(scanf("%d",&n),n)
   {
             printf("%lld\n",a[n]);
   }
   return 0;
}