#include using namespace std; #pragma region Math Osa_k struct osa_k { private: std::vector spf, pr; public: osa_k() = default; osa_k(int MAX) : spf(MAX + 1) { for(int i = 2; i <= MAX; i++) { if(spf[i] == 0) { spf[i] = i; pr.push_back(i); } for(int j = 0; j < pr.size() and pr[j] <= spf[i] and i * pr[j] <= MAX; j++) spf[i * pr[j]] = pr[j]; } } std::vector> prime_factorization(int n) { std::vector> divisor; if(n == 1) return divisor; int before = spf[n], cnt = 0; while(n > 1) { if(spf[n] == before) { cnt++; n /= spf[n]; } else { divisor.emplace_back(before, cnt); before = spf[n]; cnt = 1; n /= spf[n]; } } divisor.emplace_back(before, cnt); return divisor; } int smallestprimefactor(const int n) const { return spf[n]; } bool is_prime(const int n) const { return n == spf[n]; } }; #pragma endregion int embed[6] = {0, 99371160, 207162592, 318140852, 431159896, 545761220}; int n, res; int main() { scanf("%d", &n); osa_k pf(n); /* 愚直 for(int a = 1; a <= n; a++) { auto v = pf.prime_factorization(a); int b = 1; for(auto [val, cnt] : v) if(cnt & 1) b *= val; res += sqrt(n / b); } printf("%d\n", res); */ int base = n / 10000000 * 10000000; res = embed[n / 10000000]; for(int a = base + 1; a <= n; a++) { auto v = pf.prime_factorization(a); int b = 1; for(auto [val, cnt] : v) if(cnt & 1) b *= val; int d = sqrt(n / b); res += d + min(d, (int)sqrt(base / b)); } printf("%d\n", res); }