#define _USE_MATH_DEFINES #include using namespace std; #define FOR(i,m,n) for(int i=(m);i<(n);++i) #define REP(i,n) FOR(i,0,n) #define ALL(v) (v).begin(),(v).end() using ll = long long; constexpr int INF = 0x3f3f3f3f; constexpr long long LINF = 0x3f3f3f3f3f3f3f3fLL; constexpr double EPS = 1e-8; constexpr int MOD = 1000000007; // constexpr int MOD = 998244353; constexpr int dy[] = {1, 0, -1, 0}, dx[] = {0, -1, 0, 1}; constexpr int dy8[] = {1, 1, 0, -1, -1, -1, 0, 1}, dx8[] = {0, -1, -1, -1, 0, 1, 1, 1}; template inline bool chmax(T &a, U b) { return a < b ? (a = b, true) : false; } template inline bool chmin(T &a, U b) { return a > b ? (a = b, true) : false; } struct IOSetup { IOSetup() { std::cin.tie(nullptr); std::ios_base::sync_with_stdio(false); std::cout << fixed << setprecision(20); } } iosetup; std::vector prime_sieve(int n, bool get_only_prime) { std::vector prime, smallest_prime_factor(n + 1); std::iota(smallest_prime_factor.begin(), smallest_prime_factor.end(), 0); for (int i = 2; i <= n; ++i) { if (smallest_prime_factor[i] == i) prime.emplace_back(i); for (int p : prime) { if (i * p > n || p > smallest_prime_factor[i]) break; smallest_prime_factor[i * p] = p; } } return get_only_prime ? prime : smallest_prime_factor; } struct osa_k { std::vector smallest_prime_factor; osa_k(int n = 10000000) : smallest_prime_factor(prime_sieve(n, false)) {} std::vector> query(int n) const { std::vector> res; while (n > 1) { int prime = smallest_prime_factor[n], exponent = 0; while (smallest_prime_factor[n] == prime) { ++exponent; n /= prime; } res.emplace_back(prime, exponent); } return res; } }; int main() { constexpr int N = 50000000; int n; cin >> n; bitset is_sq = 0; vector sq; for (int i = 1; i * i <= n; ++i) { is_sq.set(i * i); sq.emplace_back(i * i); } vector spf = prime_sieve(n, false); int p = int(sq.size()) - 1; ll ans = 0; for (int i = 1; i <= n; ++i) { if (i == 1 || !is_sq[i]) { bool is_valid = true; int tmp = i; while (tmp > 1) { int prime = spf[tmp], ex = 0; while (spf[tmp] == prime && ex <= 1) { ++ex; tmp /= prime; } if (ex > 1) { is_valid = false; break; } } if (is_valid) { while (p >= 0 && 1LL * i * sq[p] > n) --p; ans += 1LL * (p + 1) * (p + 1); } } } cout << ans << '\n'; return 0; }