%I #21 May 20 2024 10:16:37
%S 0,1,11,121,10201,11111,112211,122221,1002001,1120211,11022011,
%T 100020001,101212101,122111221,1012112101,1100220011,10000200001,
%U 10111011101,110002200011,111221122111,1000002000001,1001221221001,1012200022101,1101202021011,1221221221221,10101111110101
%N Palindromes which are base-3 representations of squares.
%H Robert Israel, <a href="/A263608/b263608.txt">Table of n, a(n) for n = 1..143</a>
%H G. J. Simmons, <a href="/A002778/a002778.pdf">On palindromic squares of non-palindromic numbers</a>, J. Rec. Math., 5 (No. 1, 1972), 11-19. [Annotated scanned copy]
%p rev3:= proc(n) local L,i; L:= convert(n,base,3); add(L[-i]*3^(i-1),i=1..nops(L)) end proc:
%p c3:= proc(n) local L,i; L:= convert(n,base,3); add(L[i]*10^(i-1),i=1..nops(L)) end proc:
%p R:= 0,1: count:= 2:
%p for d from 2 while count < 100 do
%p if d::odd then
%p V:= select(issqr, [seq(seq(a*3^((d+1)/2) + b*3^((d-1)/2)+rev3(a),b=0..2),a=3^((d-3)/2) .. 3^((d-1)/2)-1)])
%p else
%p V:= select(issqr, [seq(a*3^(d/2) + rev3(a), a=3^(d/2-1) .. 3^(d/2)-1)]);
%p fi;
%p count:= count+nops(V);
%p R:= R, op(map(c3,V));
%p od:
%p R; # _Robert Israel_, May 19 2024
%Y Cf. A002778, A003166, A029984, A262607, A029985.
%Y Intersection of A001738 and A118594.
%K nonn,base,changed
%O 1,3
%A _N. J. A. Sloane_, Oct 22 2015
%E Name edited by _Robert Israel_, May 19 2024
|