OFFSET
1,2
COMMENTS
These are the squares in A080151.
x^2 is a term iff x mod 81 is in {0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 18, 21, 24, 27, 30, 33, 36, 45, 48, 51, 54, 57, 60, 63, 72, 73, 74, 75, 76, 77, 78, 79, 80}. - Robert Israel, Jul 30 2025
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
Eric Weisstein's World of Mathematics, Demlo Number
MAPLE
b:=n->sum(convert(((10^(n+1)-1)/9)^2, base, 10)[j], j=1..2*n+1): a:=proc(n) if type(sqrt(b(n)), integer)=true then b(n) else fi end: seq(a(n), n=0..2000); # Emeric Deutsch, Jun 19 2005
MATHEMATICA
A080151[n_] := (9^2)*(n/9 - FractionalPart[n/9] + FractionalPart[n/9]^2)
(* Enrique Pérez Herrero, Nov 05 2022 *)
CROSSREFS
KEYWORD
nonn,base
AUTHOR
N. J. A. Sloane, Jun 19 2005
EXTENSIONS
More terms from Emeric Deutsch, Jun 19 2005
STATUS
approved
