//#pragma GCC optimize("Ofast") //#pragma GCC optimize("unroll-loops") #include using namespace std; using ll = long long; using ull = unsigned long long; using pii = pair; template using V = vector; template using VV = V>; template V make_vec(size_t a) { return V(a); } template auto make_vec(size_t a, Ts... ts) { return V(ts...))>(a, make_vec(ts...)); } #define pb push_back #define eb emplace_back #define mp make_pair #define fi first #define se second #define rep(i, n) rep2(i, 0, n) #define rep2(i, m, n) for (int i = m; i < (n); i++) #define per(i, b) per2(i, 0, b) #define per2(i, a, b) for (int i = int(b) - 1; i >= int(a); i--) #define ALL(c) (c).begin(), (c).end() #define SZ(x) ((int)(x).size()) constexpr ll TEN(int n) { return (n == 0) ? 1 : 10 * TEN(n - 1); } template void chmin(T& t, const U& u) { if (t > u) t = u; } template void chmax(T& t, const U& u) { if (t < u) t = u; } template ostream& operator<<(ostream& os, const pair& p) { os << "(" << p.first << "," << p.second << ")"; return os; } template ostream& operator<<(ostream& os, const vector& v) { os << "{"; rep(i, v.size()) { if (i) os << ","; os << v[i]; } os << "}"; return os; } #ifdef LOCAL void debug_out() { cerr << endl; } template void debug_out(Head H, Tail... T) { cerr << " " << H; debug_out(T...); } #define debug(...) \ cerr << __LINE__ << " [" << #__VA_ARGS__ << "]:", debug_out(__VA_ARGS__) #define dump(x) cerr << __LINE__ << " " << #x << " = " << (x) << endl #else #define debug(...) (void(0)) #define dump(x) (void(0)) #endif // https://cp-algorithms.com/algebra/prime-sieve-linear.html const int maxv = 1000010; // lowest prime factor int lp[maxv]; bool pri[maxv]; vector pr; void linear_sieve() { for (int i = 2; i < maxv; ++i) pri[i] = true; for (int i = 2; i < maxv; ++i) { if (pri[i]) { lp[i] = i; pr.pb(i); } for (int j = 0; j < (int)pr.size() && pr[j] <= lp[i] && i * pr[j] < maxv; ++j) { pri[i * pr[j]] = false; lp[i * pr[j]] = pr[j]; if (i % pr[j] == 0) break; } } } constexpr ll floor_sqrt(ll n) { if (n == 0) return 0; return llround(floor(sqrt(n))); } int main() { cin.tie(nullptr); ios::sync_with_stdio(false); ll L, R; cin >> L >> R; linear_sieve(); V f(R - L + 1); V ng(R - L + 1); for (ll i = L; i <= R; ++i) { f[i - L] = i; } for (auto p : pr) { ll l = (L + p - 1) / p * p; for (ll i = l; i <= R; i += p) { ll& x = f[i - L]; x /= p; while (x % p == 0) { x /= p; ng[i - L] = 1; } } } int ans = 0; for (ll i = L; i <= R; ++i) { if (ng[i - L]) continue; ll x = floor_sqrt(f[i - L]); if (x == 1 || x * x != f[i - L]) { // debug(i); ans++; } } cout << ans << endl; return 0; }