OFFSET
1,3
COMMENTS
Only digits {0, 1, 2, 3} seem to arise.
Numbers (other than 0) that end in zero are excluded. - N. J. A. Sloane, Mar 20 2010
REFERENCES
David Wells, The Dictionary of Curious and Interesting Numbers. London: Penguin Books (1997): p. 124.
LINKS
FORMULA
Solutions to rev(x^2) = rev(x)^2.
EXAMPLE
n = 13 is a term because 31^2 = 961 = rev(169) = rev(13^2) = rev(rev(31)^2).
MATHEMATICA
rt[x_] := tn[Reverse[IntegerDigits[x]]] Do[s = rt[n^2]; s1=rt[n]^2; If[Equal[s, s1]&&!Equal[Mod[n, 10], 0], Print[{n, s, rt[s1]}]], {n, 0, 1000000}]
(* Second program: *)
Select[Range[0, 1999], Mod[#, 10] != 0 && FromDigits[Reverse[IntegerDigits[#^2]]] == FromDigits[Reverse[IntegerDigits[#]]]^2 &] (* Alonso del Arte, Oct 08 2012; corrected by Jean-François Alcover, Jan 11 2021 *)
PROG
a085305 n = a085305_list !! (n-1)
a085305_list = 0 : filter (\x -> x `mod` 10 > 0
&& a004086 (x^2) == (a004086 x)^2) [1..]
-- Reinhard Zumkeller, Jul 08 2011
(Magma) [0] cat [ m: n in [1..1810] | Reverse(Intseq(m^2)) eq Intseq(Seqint(Reverse(Intseq(m)))^2) where m is n+Floor((n-1)/9) ]; // Bruno Berselli, Jul 08 2011
(PARI) isok(x) = (x==0) || ((x%10) && fromdigits(Vecrev(digits(x^2))) == fromdigits(Vecrev(digits(x)))^2); \\ Michel Marcus, Jan 11 2021
CROSSREFS
KEYWORD
base,nonn
AUTHOR
Labos Elemer, Jun 27 2003
STATUS
approved