GCD - Extreme (II) UVA - 11426

#include<bits/stdc++.h>
using namespace std;

#define rep(i,a,b) for(int i=a;i<b;i++)
#define per(i,a,b) for(int i=b-1;i>=a;i--)
#define bug(x) printf("x:%d****\n",x)
typedef long long ll;

const int maxn=4e6+10;

int phi[maxn];
ll f[maxn];
/*
gcd(a,b)=d 的充要条件： gcd(a/d,b/d)=1 

*/


ll ans[maxn];
void init(){
    phi[1]=1;
    for(int i=2;i<maxn;i++){
        if(phi[i])continue;
        for(int j=i;j<maxn;j+=i){
            if(!phi[j])phi[j]=j;
            phi[j]-=phi[j]/i;
        }
    }

    rep(i,1,maxn){
        for(int j=i*2;j<maxn;j+=i)
            f[j]+=(ll)i*phi[j/i];
    }
    rep(i,1,maxn)ans[i]=ans[i-1]+f[i];
}


int main(){

   // double t1=clock();
    init();
   // printf("%.2f\n",(clock()-t1)/CLOCKS_PER_SEC);

    int n;
    while(scanf("%d",&n)==1&&n){
        printf("%lld\n",ans[n]);
    }
    return 0;
}