#include #define REP_(i, a_, b_, a, b, ...) \ for (int i = (a), END_##i = (b); i < END_##i; ++i) #define REP(i, ...) REP_(i, __VA_ARGS__, __VA_ARGS__, 0, __VA_ARGS__) #define ALL(x) std::begin(x), std::end(x) using i64 = long long; template inline bool chmax(T &a, U b) { return a < b and ((a = std::move(b)), true); } template inline bool chmin(T &a, U b) { return a > b and ((a = std::move(b)), true); } template inline int ssize(const T &a) { return (int)std::size(a); } template std::istream &operator>>(std::istream &is, std::vector &a) { for (auto &x : a) is >> x; return is; } template std::ostream &operator<<(std::ostream &os, const std::pair &a) { return os << "(" << a.first << ", " << a.second << ")"; } template std::ostream &print_seq(const Container &a, std::string_view sep = " ", std::string_view ends = "\n", std::ostream &os = std::cout) { auto b = std::begin(a), e = std::end(a); for (auto it = std::begin(a); it != e; ++it) { if (it != b) os << sep; os << *it; } return os << ends; } template struct is_iterable : std::false_type {}; template struct is_iterable())), decltype(std::end(std::declval()))>> : std::true_type {}; template ::value && !std::is_same::value && !std::is_same::value>> std::ostream &operator<<(std::ostream &os, const T &a) { return print_seq(a, ", ", "", (os << "{")) << "}"; } #ifdef ENABLE_DEBUG #include "debug_dump.hpp" #else #define DUMP(...) #endif using namespace std; // Returns a vector that maps x to x's smallest prime factor (> 1). // O(n log log n) std::vector sieve_smallest_prime_factors(int n) { std::vector res(n + 1); for (int i = 1; i <= n; ++i) res[i] = i; for (int i = 2; i <= n; i += 2) res[i] = 2; for (int i = 3; i * i <= n; i += 2) { if (res[i] != i) continue; for (int j = i * i; j <= n; j += i) { if (res[j] == j) res[j] = i; } } return res; } // Quick factorization. i64 quick_factorize(int n, const std::vector &spf) { // assert(0 < n and n < int(spf.size())); // std::vector> res; i64 b = 1; for (;;) { const int p = spf[n]; if (p == 1) break; int count = 0; do { n /= p; ++count; } while (n % p == 0); if (count & 1) { b *= p; } } return b; } auto solve() { i64 n; cin >> n; auto spf = sieve_smallest_prime_factors(n); i64 ans = 0; for (i64 x = 1; x <= n; ++x) { i64 b = quick_factorize(x, spf); i64 ub = n / b; i64 tv = 1, fv = n + 1; while (fv - tv > 1) { i64 mid = (tv + fv) / 2; if (mid * mid <= ub) { tv = mid; } else { fv = mid; } } ans += tv; } return ans; } int main() { ios_base::sync_with_stdio(false), cin.tie(nullptr); cout << solve() << "\n"; }