import std.conv, std.functional, std.range, std.stdio, std.string; import std.algorithm, std.array, std.bigint, std.bitmanip, std.complex, std.container, std.math, std.mathspecial, std.numeric, std.regex, std.typecons; import core.bitop; class EOFException : Throwable { this() { super("EOF"); } } string[] tokens; string readToken() { for (; tokens.empty; ) { if (stdin.eof) { throw new EOFException; } tokens = readln.split; } auto token = tokens.front; tokens.popFront; return token; } int readInt() { return readToken.to!int; } long readLong() { return readToken.to!long; } real readReal() { return readToken.to!real; } bool chmin(T)(ref T t, in T f) { if (t > f) { t = f; return true; } else { return false; } } bool chmax(T)(ref T t, in T f) { if (t < f) { t = f; return true; } else { return false; } } int binarySearch(alias pred, T)(in T[] as) { int lo = -1, hi = cast(int)(as.length); for (; lo + 1 < hi; ) { const mid = (lo + hi) >> 1; (unaryFun!pred(as[mid]) ? hi : lo) = mid; } return hi; } int lowerBound(T)(in T[] as, T val) { return as.binarySearch!(a => (a >= val)); } int upperBound(T)(in T[] as, T val) { return as.binarySearch!(a => (a > val)); } // floor(sqrt(a)) long floorSqrt(long a) { import core.bitop : bsr; import std.algorithm : min; long b = a, x = 0, y = 0; for (int e = bsr(a) & ~1; e >= 0; e -= 2) { x <<= 1; y <<= 1; if (b >= (y | 1) << e) { b -= (y | 1) << e; x |= 1; y += 2; } } return x; } long N; long sqrtN; bool[] isPrime; int primesLen; long[] primes; long[] small, large, large1; long ans; long get(long n) { if (n <= sqrtN + 1) return small[n]; if (n <= N && n == N / (N / n)) return large[N / n]; if (n <= N + 1 && n == N / (N / (n - 1)) + 1) return large1[N / (n - 1)]; assert(false); } void init() { sqrtN = floorSqrt(N + 1); isPrime = new bool[sqrtN + 1 + 1]; small = new long[sqrtN + 1 + 1]; large = new long[sqrtN + 1]; large1 = new long[sqrtN + 1]; isPrime[2 .. $] = true; primes = []; foreach (n; 1 .. sqrtN + 1 + 1) small[n] = n; foreach (l; 1 .. sqrtN + 1) large[l] = N / l; foreach (l; 1 .. sqrtN + 1) large1[l] = N / l + 1; foreach (p; 2 .. sqrtN + 1 + 1) { if (isPrime[p]) { primes ~= p; for (long n = p^^2; n <= sqrtN + 1; n += p) isPrime[n] = false; const g1 = get(p - 1); foreach (l; 1 .. sqrtN + 1) { { const n = N / l + 1; if (n < p^^2) break; large1[l] -= (get(n / p) - g1); } { const n = N / l; if (n < p^^2) break; large[l] -= (get(n / p) - g1); } } foreach_reverse (n; 1 .. sqrtN + 1 + 1) { if (n < p^^2) break; small[n] -= (get(n / p) - g1); } } } primesLen = cast(int)(primes.length); small[1 .. $] -= 1; large[1 .. $] -= 1; large1[1 .. $] -= 1; debug { writeln("sqrtN = ", sqrtN); writeln("primes = ", primes); writeln("small = ", small); writeln("large = ", large); writeln("large1 = ", large1); stdout.flush; } } void dfs(int pos, long n) { // primes[i] <= n if (pos >= 0) { debug { writefln("dfs %s %s: +1", pos, n); } ans += 1; } // primes[i] - 1 <= n debug { writefln("dfs %s %s: +(pi(%s) - %s)", pos, n, n + 1, pos + 1); } ans += max(get(n + 1) - (pos + 1), 0); if (pos >= 0) { const nn = n / primes[pos]; if (primes[pos] <= nn) { dfs(pos, nn); } } foreach (i; pos + 1 .. primesLen) { const nn = n / (primes[i] - 1); if (primes[i] > nn) { break; } dfs(i, nn); } } void main() { try { for (; ; ) { N = readLong(); init; ans = 1; dfs(-1, N); writeln(ans); } } catch (EOFException e) { } }