#include using namespace std; int main(){ int N; cin >> N; int U=pow(N,1.0/3); int L=N/(U+1)/(U+1); vectorlarge(U+2),small(L+1); // small[i]=1以上 i 以下の整数のうち、平方因子を持たないものの個数 // large[i]=1以上N/i/i以下の整数のうち、平方因子を持たないものの個数 int ans=0; for(int i=1;i<=L;i++){ small[i]=i; for(int d=2;d*d<=i;d++)small[i]-=small[i/d/d]; int cnt=sqrt(N/i); ans+=cnt*cnt*(small[i]-small[i-1]); } large[U+1]=small[L]; for(int i=U;i>0;i--){ int R=N/i/i; large[i]=R; for(int d=2;d*d<=R;d++)large[i]-=i*d<=U?large[i*d]:small[R/d/d]; ans+=i*i*(large[i]-large[i+1]); } cout << ans; }