#include <bits/stdc++.h>

using namespace std;
ifstream f("fractii.in");
ofstream g("fractii.out");
int n,phi[1000001];
int64_t sol;
int main()
{
    f>>n;
    for(int i=1;i<=n;i++)
        phi[i]=i;
    for(int i=2;i<=n;i++)
    {
        if(phi[i]==i)
            for(int j=i;j<=n;j+=i)
                phi[j]=phi[j]/i*(i-1);
        sol+=phi[i];
    }
    sol=2*sol+1;
    g<<sol;
    return 0;
}