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