#include using namespace std; #define DEBUG(x) cerr<<#x<<": "< #define vl vector #define vii vector< vector > #define vll vector< vector > #define vs vector #define pii pair #define pis pair #define psi pair #define pll pair template pair operator+(const pair &s, const pair &t) { return pair(s.first + t.first, s.second + t.second); } template pair operator-(const pair &s, const pair &t) { return pair(s.first - t.first, s.second - t.second); } template ostream& operator<<(ostream& os, pair p) { os << "(" << p.first << ", " << p.second << ")"; return os; } #define X first #define Y second #define rep(i,n) for(int i=0;i<(int)(n);i++) #define rep1(i,n) for(int i=1;i<=(int)(n);i++) #define rrep(i,n) for(int i=(int)(n)-1;i>=0;i--) #define rrep1(i,n) for(int i=(int)(n);i>0;i--) #define REP(i,a,b) for(int i=a;i bool chmax(T &a, const T &b) { if (a bool chmin(T &a, const T &b) { if (a>b) { a = b; return 1; } return 0; } #define UNIQUE(v) v.erase(std::unique(v.begin(), v.end()), v.end()); const ll inf = 1000000001; const ll INF = (ll)1e18 + 1; const long double pi = 3.1415926535897932384626433832795028841971L; #define Sp(p) cout< is_prime(33); ll nibu(ll p, ll n) { ll ok = 1, ng = 1e8; while (ok + 1 < ng) { ll x = (ok + ng) / 2; ll mul = 1; rep (_, p) { mul *= x; if (mul > n) break; } if (mul > n) ng = x; else ok = x; } return ok; } ll gcd(ll a, ll b) { if (b > a) { swap(a, b); } ll r = a%b; while (r != 0) { a = b; b = r; r = a%b; } return b; } ll lcm(ll a, ll b) { return (a / gcd(a, b))*b; } void solve() { ll n; cin >> n; ll ans = n * n; // a = c = 1; ans += (n - 1) * n; // a = c != 1; for (int p = 2; p <= 31; p++) { //if (not is_prime[p]) continue; for (int q = 1; q < p; q++) { //if (q != 1 and not is_prime[q]) continue; if (gcd(p, q) != 1) continue; ll y = nibu(p, n); ll g = gcd(p, q); ll p2 = p / g; ll d = n / p2; if (y != 1) { //DEBUG_VEC(vl({p, q, y, d})); ans += d * (y - 1) * 2; } } } cout << ans << endl; } signed main() { fio(); for (int i = 2; i <= 32; i++) { int x = i; is_prime[i] = true; for (int j = 2; j < x; j++) { if (x % j == 0) is_prime[i] = false; } } int t; cin >> t; while (t--) { solve(); } }