login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A177950
Numbers k that divide the digit reversal of k^2.
3
1, 2, 3, 9, 11, 22, 26, 33, 48, 66, 87, 99, 101, 111, 117, 121, 202, 212, 216, 264, 273, 288, 297, 307, 333, 484, 513, 528, 666, 783, 819, 836, 999, 1001, 1111, 1323, 1331, 1452, 1474, 1602, 2002, 2178, 2285, 2622, 2636, 3333, 4884, 4961, 6666, 7161, 7575
OFFSET
1,2
LINKS
Harvey P. Dale, Table of n, a(n) for n = 1..171 (all terms up to k = 10 million)
EXAMPLE
273 is in the sequence because 273^2 = 74529 and 273|92547.
MAPLE
with(numtheory):T:=array(1..94):k:=1:for n from 1 to 200000 do: n0:=n^2:l:=length(n0) :s:=0:for m from l by -1 to 1 do:q:=n0:u:=irem(q, 10):v:=iquo(q, 10): n0:=v :s:=s+ u*10^(m-1):od:if irem(s, n)=0 then T[k]:=n:k:=k+1:else fi: od:print(T):
MATHEMATICA
Select[Range[1000], Mod[FromDigits[Reverse[IntegerDigits[ (#)^2 ]]], # ]==0 & ]
Select[Range[8000], Divisible[IntegerReverse[#^2], #]&] (* Requires Mathematica version 10 or later *) (* Harvey P. Dale, Dec 12 2016 *)
PROG
(PARI) isok(k) = !(fromdigits(Vecrev(digits(k^2))) % k); \\ Michel Marcus, Aug 20 2021
CROSSREFS
Sequence in context: A369338 A110350 A057569 * A049618 A057292 A098016
KEYWORD
nonn,base
AUTHOR
Michel Lagneau, May 15 2010
EXTENSIONS
Edited by Charles R Greathouse IV, Aug 02 2010
STATUS
approved