#include //ios::sync_with_stdio(false);cin.tie(0);cout.tie(0); //clock_t start=clock();clock_t end=clock();cout<<(double)(end-start)/CLOCKS_PER_SEC< pii; typedef pair ppii; typedef pair pipi; typedef pair pll; typedef pair ppll; typedef pair plpl; typedef pair P; typedef tuple tl; typedef pair pdd; typedef vector> mat; ll mod=1000000007; ll mod2=998244353; ll mod3=1000003; ll mod4=998244853; ll mod5=1000000009; ll inf=1LL<<60; double pi=3.141592653589793238462643383279L; double eps=1e-14; #define rep(i,m,n) for(ll i=m;i=m;i--) #define srep(itr,st) for(auto itr=st.begin();itr!=st.end();itr++) #define mrep(itr,mp) for(auto itr=mp.begin();itr!=mp.end();itr++) #define Max(a,b) a=max(a,b) #define Min(a,b) a=min(a,b) int dh[4]={1,-1,0,0}; int dw[4]={0,0,1,-1}; int ddh[8]={-1,-1,-1,0,0,1,1,1}; int ddw[8]={-1,0,1,-1,1,-1,0,1}; struct custom_hash { static uint64_t splitmix64(uint64_t x) { // http://xorshift.di.unimi.it/splitmix64.c x += 0x9e3779b97f4a7c15; x = (x ^ (x >> 30)) * 0xbf58476d1ce4e5b9; x = (x ^ (x >> 27)) * 0x94d049bb133111eb; return x ^ (x >> 31); } size_t operator()(uint64_t x) const { static const uint64_t FIXED_RANDOM = chrono::steady_clock::now().time_since_epoch().count(); return splitmix64(x + FIXED_RANDOM); } }; #define umh unordered_map ll gcd(ll a,ll b){ if(a0){ if(k&1)ret*=now; now*=now; k/=2; } return ret; } ll beki(ll n,ll k,ll md){ ll ret=1; ll now=n; now%=md; while(k>0){ if(k%2==1){ ret*=now; ret%=md; } now*=now; now%=md; k/=2; } return ret; } ll gyaku(ll n,ll md){ return beki(n,md-2,md); } ll popcount(ll n){ ll ret=0; ll u=n; while(u>0){ ret+=u%2; u/=2; } return ret; } int main(){ ios::sync_with_stdio(false);cin.tie(0); int t;cin>>t; int used[100000]; fill(used,used+100000,1); bool used2[100000]; fill(used2,used2+100000,false); vector v,w; rep(i,2,100000){ if(used2[i])continue; rep(j,1,100000){ ll now=i*j; if(now>100000)break; used2[now]=true; if(gcd(i,j)==1)used[now]*=1; else used[now]*=0; } } rep(i,2,100000){ if(i*i>1000000000)break; if(used[i])v.push_back(i); } //rep(i,0,20)cout<>n; ull ans=n*n; ll cnt=1; rep(i,2,100000){ if(i*i>n)break; if(used[i])continue; ll ma=1; ll pp=i; used[i]=true; while(pp<=n){ pp*=i; if(pp<100000)used[pp]=true; ma++; } cnt+=ma-1; rep(j,1,ma){ rep(k,1,ma){ ll f=max(j,k); f/=gcd(j,k); ans+=n/f; } } } //rep(i,2,n+1)cout<