OFFSET
0,3
LINKS
Robert Israel, Table of n, a(n) for n = 0..2093
EXAMPLE
a(1) = ( 1)_3 = 1,
a(2) = ( 21)_3 = 7,
a(3) = ( 121)_3 = 16,
a(4) = (2121)_3 = 70.
MAPLE
N:= 30: # for a(0) to a(N)
with(padic):
A:= rootp(x^2+1/2, 3, N):
if ratvaluep(A[1], 1) = 1 then A:= A[1] else A:= A[2] fi:
seq(ratvaluep(A, i), i=0..N); # Robert Israel, May 11 2020
PROG
(PARI) {a(n) = truncate(sqrt(-1/2+O(3^n)))}
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, Aug 04 2019
STATUS
approved