login
A088526
a(n) = A061017(n) - pi(n) - 1.
4
0, 0, -1, 0, -1, 0, -1, -1, 0, 0, 0, 0, -1, -1, 0, 0, 0, 0, -1, -1, 0, 0, -1, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 1, 0, 1, 1, 1, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 0, 1, 1, 2, 2, 1, 1, 1, 1, 0, 0, 1, 1, 1, 2, 1, 1, 1, 2, 2, 2, 2, 3, 2, 2, 2, 2, 1, 2, 1, 2, 2
OFFSET
1,77
COMMENTS
a(n) >= 0 for all n >= 24.
The old entry with this sequence number was a duplicate of A074292.
PROG
(Python)
from math import isqrt
from sympy import primepi
from oeis_sequences.OEISsequences import bisection, bsearch
def A088526(n):
def g(x): return x-(s:=isqrt(x))**2+(sum(x//k for k in range(1, s+1))<<1)
def f(x): return n+bsearch(g, x)
return bisection(f, n, n)-n-primepi(n) # Chai Wah Wu, Mar 01 2026
CROSSREFS
KEYWORD
sign
AUTHOR
STATUS
approved