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
David Applegate and N. J. A. Sloane, Oct 22 2008
STATUS
approved
