|
|
A269588
|
|
Numbers n such that n^2 ends with the digits of n reversed (A004086(n)).
|
|
2
|
|
|
1, 5, 6, 963, 9867, 65766, 69714, 6317056, 90899553, 169605719, 4270981082, 96528287587, 465454256742, 692153612536, 182921919071841, 655785969669834, 650700037578750084, 125631041500927357539, 673774165549097456624, 16719041449406813636569
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
COMMENTS
|
a(29)>10^32 (if it exists)
|
|
LINKS
|
|
|
EXAMPLE
|
6317056^2 = 39905196507136 which ends with 6507136, so 6317056 is a term.
|
|
MATHEMATICA
|
Select[Range[10^7], Function[k, Take[IntegerDigits[#^2], -Length@ k] == Reverse@ k]@ IntegerDigits@ # &] (* Michael De Vlieger, Mar 04 2016 *)
|
|
PROG
|
(PARI) isA269588(n)=dn = digits(n); rn = subst(Polrev(dn), x, 10); nbd = #dn; (n^2 - rn) % 10^nbd == 0; \\ Michel Marcus, Mar 01 2016
(PARI) \\ printA269588len(d) prints all terms of the sequence with d digits
rev(n) = eval(concat(Vecrev(Str(n))));
{ printA269588len(d) = my(l, u, n); l=ceil(d/2); u=floor(d/2); for(y=0, 10^l-1, n=rev(y^2 % 10^u)*10^l+y; if(#Str(n)==d && Mod(n, 10^d)^2==rev(n), print(n)); ); }
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,base
|
|
AUTHOR
|
|
|
EXTENSIONS
|
a(28) from Dieter Beckerle, Jun 09 2016
|
|
STATUS
|
approved
|
|
|
|