Eric Weisstein's World of Mathematics, <a href="httphttps://mathworld.wolfram.com/OddPower.html">Odd Powers</a>.
Eric Weisstein's World of Mathematics, <a href="httphttps://mathworld.wolfram.com/OddPower.html">Odd Powers</a>.
reviewed
approved
proposed
reviewed
editing
proposed
return kmax # Chai Wah Wu, Aug 14 2024
(Python)
from sympy import mobius, integer_nthroot
def A070265(n):
def f(x): return int(n-2+x+sum(mobius(k)*(integer_nthroot(x, k)[0]-1) for k in range(3, x.bit_length(), 2)))
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
approved
editing
reviewed
approved
proposed
reviewed
editing
proposed
Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/OddPower.html">Odd Powers</a>.