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”).
%I #35 Aug 11 2024 14:41:30
%S 0,1,2,3,4,5,6,7,8,9,11,323,464,717,858,999,1111,39593,59095,420024,
%T 546645,9046409,9578759,9813189,535505535,564303465,595121595,
%U 5736116375,6758008576,10476867401,11652825611,14203330241
%N Numbers that are palindromic in bases 10 and 14.
%H Robert G. Wilson v, <a href="/A029969/b029969.txt">Table of n, a(n) for n = 1..70</a> (first 63 terms from Ray Chandler)
%H P. De Geest, <a href="https://www.worldofnumbers.com/nobase10.htm">Palindromic numbers beyond base 10</a>
%t NextPalindrome[n_] := Block[{l = Floor[ Log[10, n] + 1], idn = IntegerDigits[n]}, If[ Union[idn] == {9}, Return[n + 2], If[l < 2, Return[n + 1], If[ FromDigits[ Reverse[ Take[idn, Ceiling[l/2]] ]] FromDigits[ Take[idn, -Ceiling[l/2]]], FromDigits[ Join[ Take[idn, Ceiling[l/2]], Reverse[ Take[idn, Floor[l/2]] ]]], idfhn = FromDigits[ Take[idn, Ceiling[l/2]]] + 1; idp = FromDigits[ Join[ IntegerDigits[idfhn], Drop[ Reverse[ IntegerDigits[idfhn]], Mod[l, 2]] ]]] ]]]; palQ[n_Integer, base_Integer] := Block[{idn = IntegerDigits[n, base]}, idn == Reverse[idn]]; l = {0}; a = 0; Do[a = NextPalindrome[a]; If[ palQ[a, 14], AppendTo[l, a]], {n, 300000}]; l (* _Robert G. Wilson v_, Sep 03 2004 *)
%t b1=10; b2=14; lst={}; Do[d1=IntegerDigits[n, b1]; d2=IntegerDigits[n, b2]; If[d1==Reverse[d1]&&d2==Reverse[d2], AppendTo[lst, n]], {n, 0, 1000000}]; lst (* _Vincenzo Librandi_, Nov 23 2014 *)
%t palQ[n_]:=Module[{idn14=IntegerDigits[n,14]},n==IntegerReverse[n]&&idn14==Reverse[idn14]]; Select[Range[10^7],palQ] (* The program uses the IntegerReverse function from Mathematica version 10 *) (* _Harvey P. Dale_, Apr 23 2016 *)
%t Select[Range[0, 10^5],
%t PalindromeQ[#] && # == IntegerReverse[#, 14] &] (* _Robert Price_, Nov 09 2019 *)
%o (Magma) [n: n in [0..10000000] | Intseq(n, 10) eq Reverse(Intseq(n, 10))and Intseq(n, 14) eq Reverse(Intseq(n, 14))]; // _Vincenzo Librandi_, Nov 23 2014
%Y Cf. A007632, A007633, A029961, A029962, A029963, A029964, A029804, A029965, A029966, A029967, A029968, A029970, A029731, A097855, A099165.
%K nonn,base
%O 1,3
%A _Patrick De Geest_