#define _CRT_SECURE_NO_WARNINGS #pragma GCC target("avx") #pragma GCC optimize("O3") #pragma GCC optimize("unroll-loops") #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #define rep(i, n) for(int i=0;i<(n);i++) #define REP(i, n) for(int i=1;i<=(n);i++) #define all(V) V.begin(),V.end() typedef long long lint; typedef std::pair P; constexpr int INF = INT_MAX / 10; constexpr lint LINF = LLONG_MAX / 10; constexpr double eps = 1e-9; template class prique :public std::priority_queue, std::greater> {}; template inline bool chmax(T& lhs, const U& rhs) { if (lhs < rhs) { lhs = rhs; return 1; } return 0; } template inline bool chmin(T& lhs, const U& rhs) { if (lhs > rhs) { lhs = rhs; return 1; } return 0; } inline lint gcd(lint a, lint b) { while (b) { lint c = a; a = b; b = c % b; } return a; } inline lint lcm(lint a, lint b) { return a / gcd(a, b) * b; } bool isprime(lint n) { if (n == 1)return false; for (int i = 2; i * i <= n; i++) { if (n % i == 0)return false; } return true; } lint mypow(lint a, lint b) { if (!b)return 1; if (b & 1)return mypow(a, b - 1) * a; lint memo = mypow(a, b >> 1); return memo * memo; } lint modpow(lint a, lint b, lint m) { if (!b)return 1; if (b & 1)return modpow(a, b - 1, m) * a % m; lint memo = modpow(a, b >> 1, m); return memo * memo % m; } void printArray(std::vector& vec) { rep(i, vec.size() - 1)std::cout << vec[i] << " "; std::cout << vec.back() << std::endl; } template void printArray(T l, T r) { T rprev = r; rprev--; for (T i = l; i != rprev; i++) { std::cout << *i << " "; } std::cout << *rprev << std::endl; } std::string to_string(std::vector& vec) { std::string res = "["; rep(i, vec.size() - 1)res += std::to_string(vec[i]) + ", "; res += std::to_string(vec.back()) + "]"; return res; } int t; lint n; bool used[100000]; int main() { std::cin >> t; rep(i, t) { std::cin >> n; lint ans = n * n + (n - 1) * n; std::fill(used, used + 100000, false); for (lint j = 2; j * j <= n; j++) { if (used[j])continue; lint m = j; while (m * m <= n) { used[m] = true; m *= j; } for (lint k = 1; mypow(j, k) <= n; k++) { for (lint l = 1; mypow(j, l) <= n; l++) { if (k == l)continue; ans += n / (std::max(k, l) / gcd(k, l)); } } } std::cout << ans << std::endl; } return 0; }