OFFSET
1,1
LINKS
Alois P. Heinz, Table of n, a(n) for n = 1..10000
EXAMPLE
184 = 20211_3 and 184^2 = 1201102221_3.
MAPLE
p:= n-> add(x^i, i=convert(n, base, 3)):
a:= proc(n) option remember; local k;
for k from 1+`if`(n=1, 0, a(n-1))
while p(k)*2<>p(k^2) do od; k
end:
seq(a(n), n=1..50); # Alois P. Heinz, May 10 2015
CROSSREFS
KEYWORD
base,easy,nonn
AUTHOR
Erich Friedman, Jun 16 2001
EXTENSIONS
Offset changed to 1 by Alois P. Heinz, May 10 2015
STATUS
approved