OFFSET
1,1
LINKS
Giovanni Resta, Table of n, a(n) for n = 1..39 (terms < 7^29)
EXAMPLE
19865^2 = 394618225 = 12531|12531|1 in base 7.
411772^2 = 169556179984 = 1515152|1515152 in base 7.
MATHEMATICA
psq[b_, nd_] := Block[{sq={}, r, x}, Do[r = Mod[nd, p]; Do[x = (b^( Floor[nd/p] p) - 1)/(b^p - 1) t b^r + Floor[t/b^(p - r)]; If[ IntegerQ@ Sqrt@ x, AppendTo[ sq, x]], {t, b^(p - 1), b^p - 1}], {p, nd/2}]; sq]; Union @@ (psq[7, #] & /@ Range[12]) (* Giovanni Resta, Aug 01 2018 *)
CROSSREFS
KEYWORD
base,nonn
AUTHOR
Erich Friedman, May 28 2002
EXTENSIONS
a(17)-a(21) from Giovanni Resta, Aug 01 2018
STATUS
approved