#pragma region Macros #include using namespace std; using namespace __gnu_pbds; // using namespace __gnu_cxx; // #include // namespace mp = boost::multiprecision; // using Bint = mp::cpp_int; #define TO_STRING(var) # var #define pb emplace_back #define int ll #define endl '\n' using ll = long long; using ld = long double; const ld PI = acos(-1); const ld EPS = 1e-10; const int INF = 1 << 30; const ll INFL = 1LL << 61; // const int MOD = 998244353; const int MOD = 1000000007; __attribute__((constructor)) void constructor() { ios::sync_with_stdio(false); cin.tie(nullptr); // ifstream in("input.txt"); // cin.rdbuf(in.rdbuf()); cout << fixed << setprecision(15); } template class modint{ public: int val = 0; modint(int x = 0) { while (x < 0) x += mod; val = x % mod; } modint(const modint &r) { val = r.val; } modint operator -() { return modint(-val); } modint operator +(const modint &r) { return modint(*this) += r; } modint operator -(const modint &r) { return modint(*this) -= r; } modint operator *(const modint &r) { return modint(*this) *= r; } modint operator /(const modint &r) { return modint(*this) /= r; } modint &operator +=(const modint &r) { val += r.val; if (val >= mod) val -= mod; return *this; } modint &operator -=(const modint &r) { if (val < r.val) val += mod; val -= r.val; return *this; } modint &operator *=(const modint &r) { val = val * r.val % mod; return *this; } modint &operator /=(const modint &r) { int a = r.val, b = mod, u = 1, v = 0; while (b) { int t = a / b; a -= t * b; swap(a, b); u -= t * v; swap(u, v); } val = val * u % mod; if (val < 0) val += mod; return *this; } bool operator ==(const modint& r) { return this -> val == r.val; } bool operator <(const modint& r) { return this -> val < r.val; } bool operator !=(const modint& r) { return this -> val != r.val; } }; using mint = modint; istream &operator >>(istream &is, mint& x) { int t; is >> t; x = t; return (is); } ostream &operator <<(ostream &os, const mint& x) { return os << x.val; } mint modpow(const mint &a, int n) { if (n == 0) return 1; mint t = modpow(a, n / 2); t = t * t; if (n & 1) t = t * a; return t; } int modpow(int x, int n, int mod) { int ret = 1; while (n > 0) { if (n % 2 == 1) ret = ret * x % mod; x = x * x % mod; n /= 2; } return ret; } int ceil(int x, int y) { return (x > 0 ? (x + y - 1) / y : x / y); } #pragma endregion vector erat(int n) { vector prime(n + 1, true); prime[0] = false; if (n >= 1) prime[1] = false; for (int i = 2; i * i <= n; i++) { if (prime[i]) { for (int j = i + i; j <= n; j += i) { prime[j] = false; } } } return prime; } vector prime_fact(int x) { vector ret; if (x == 1) { ret.push_back(1); return ret; } for (int i = 2; i * i <= x; i++) { while (x % i == 0) { ret.push_back(i); x /= i; } } if (x > 1) ret.push_back(x); return ret; } signed main() { int L, R; cin >> L >> R; vector E = erat(2e6); vector P; for (int i = 2; i < 2e6; i++) { if (E[i]) P.pb(i); } vector P2; for (auto &x : P) { P2.pb(x * x); } int sz = P.size(); int sz2 = P2.size(); vector V; for (int i = L; i <= R; i++) { bool flag = false; for (int j = 0; j < sz2; j++) { if (i % P2[j] == 0) { V.pb(i); flag = true; break; } } if (flag) continue; for (int j = 0; j < sz; j++) { if (i != P[j] && i % P[j] == 0) { int p = round(sqrt(i / P[j])); if (p * p * P[j] == i) { V.pb(i); break; } } } } sort(V.begin(),V.end()); V.erase(unique(V.begin(),V.end()),V.end()); cout << R - L + 1 - V.size() << endl; // for (auto &x : mp) { // cout << x.first << " " << x.second << endl; // } }