OFFSET
0,2
LINKS
Paolo Xausa, Table of n, a(n) for n = 0..1000
MATHEMATICA
Round[(8^(1/5))^Range[0, 50]] (* Paolo Xausa, Jul 01 2024 *)
PROG
(Python)
from gmpy2 import iroot_rem
def A018136(n):
i, j = iroot_rem(1<<3*n, 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