#include<cstdio>
using namespace std;
const int NMAX = 1000005;
int n,i,j,phi[NMAX];
long long 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++)
    {
        if(phi[i]==i)
            for(j=i;j<=n;j+=i) phi[j]=phi[j]/i*(i-1);
        sol+=2*phi[i];
    }
    printf("%lld\n",sol+1);
    return 0;
}