OFFSET
0,2
MATHEMATICA
Floor[(16^(1/5))^Range[0, 40]+1/2] (* Harvey P. Dale, May 23 2012 *)
PROG
(Python)
from gmpy2 import iroot_rem
def A018160(n):
i, j = iroot_rem(1<<(n<<2), 5)
return int(i)+int(j<<5>=10*i*((i*((i*(i+1)<<1)+1)<<2)+1)+1) # Chai Wah Wu, Jun 20 2024
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved