gcd(1,n)+gcd(2,n)....gcd(n-1,n); Uva11426

#include<bits/stdc++.h>
#define int long long
using namespace std;
const int maxn=4e6+;
int phi[maxn];
int prime[maxn];
int visit[maxn];
int tot=;
int num[maxn];
int ans[maxn];
void build_phi(int n)
{
    phi[]=;
    for(int i=;i<=n;i++)
    {
        if(!visit[i])
        {
          prime[++tot]=i;
           phi[i]=i-;
        }
        for(int j=;j<=tot && i*prime[j]<=n;j++)
        {
            visit[i*prime[j]]=;
            if(i%prime[j]==)
            {
                phi[i*prime[j]]=phi[i]*prime[j]; break;
            }
            else
            {
                phi[i*prime[j]]=phi[i]*(prime[j]-);
            }
        }
    }
}
void solve(int n)
{
        for(int i=;i<=n;i++)
          for(int j=i*;j<=n;j+=i)
            num[j] += phi[j/i]*i;
        for(int i=;i<=n;i++)
        {
            ans[i]=ans[i-]+num[i];
        }

}
int32_t main()
{
    build_phi(maxn);
    solve(maxn);
    //cout<<phi[1000000]<<endl;
    int n;
    while(cin>>n&&n)
    {
        cout<<ans[n]<<endl;
    }
    return ;
}