OFFSET
1,2
MATHEMATICA
t = Union@ Flatten@ Join[{1}, Table[ n^i, {n, 2, Sqrt[5775]}, {i, 2, Log[n, 5775]}]]; t[[2# - 1]] & /@ Range@(Length@t/2) (* Robert G. Wilson v *)
PROG
(Python)
from sympy import mobius, integer_nthroot
def A099997(n):
def f(x): return int((n<<1)-3+x+sum(mobius(k)*(integer_nthroot(x, k)[0]-1) for k in range(2, x.bit_length())))
kmin, kmax = 1, 2
while f(kmax) >= kmax:
kmax <<= 1
while True:
kmid = kmax+kmin>>1
if f(kmid) < kmid:
kmax = kmid
else:
kmin = kmid
if kmax-kmin <= 1:
break
return kmax # Chai Wah Wu, Aug 14 2024
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
N. J. A. Sloane, Nov 20 2004
EXTENSIONS
More terms from Robert G. Wilson v, Dec 14 2005
STATUS
approved