OFFSET
1,1
COMMENTS
MATHEMATICA
lim=2025; Select[Union[Range[Floor[lim^(1/2)]]^2, Range[Floor[lim^(1/3)]]^3], !IntegerQ[#^(1/6)]&] (* James C. McMahon, Oct 16 2024 *)
PROG
(Python)
from math import isqrt
from sympy import integer_nthroot
def A377025(n):
def bisection(f, kmin=0, kmax=1):
while f(kmax) > kmax: kmax <<= 1
while kmax-kmin > 1:
kmid = kmax+kmin>>1
if f(kmid) <= kmid:
kmax = kmid
else:
kmin = kmid
return kmax
def f(x): return n+x+(integer_nthroot(x, 6)[0]<<1)-integer_nthroot(x, 3)[0]-isqrt(x)
return bisection(f, n, n)
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Chai Wah Wu, Oct 13 2024
STATUS
approved