#include<cstdio>
using namespace std;
const int NMAX = 1000005;
int n,i,j;
long long phi[NMAX],sol;
int main()
{
    freopen("fractii.in","r",stdin);
    freopen("fractii.out","w",stdout);
    scanf("%d",&n);
    for(i=2;i<=n;i++) phi[i]=i;
    for(i=2;i<=n;i++)
    {
        sol+=2*phi[i];
        if(phi[i]==i)
            for(j=i;j<=n;j+=i) phi[j]=phi[j]*(i-1)/i;
    }
    printf("%lld\n",sol+1);
    return 0;
}