login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A099165 Palindromic in bases 10 and 32. 35

%I #30 Mar 22 2020 17:19:19

%S 0,1,2,3,4,5,6,7,8,9,11,22,33,66,99,363,858,1441,2882,5445,6886,9449,

%T 15951,19891,21012,29692,32223,54945,369963,477774,564465,585585,

%U 609906,672276,717717,780087,804408,912219,1251521,2639362,3825283

%N Palindromic in bases 10 and 32.

%H Ray Chandler and Robert G. Wilson v, <a href="/A099165/b099165.txt">Table of n, a(n) for n = 1..115</a>, terms a(88)-a(111) from Ray Chandler.

%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, 32], AppendTo[l, a]], {n, 10000}]; l

%t Select[Range[0, 10^5],

%t PalindromeQ[#] && # == IntegerReverse[#, 32] &] (* _Robert Price_, Nov 09 2019 *)

%o (Python)

%o from gmpy2 import digits

%o def palQ(n,b): # check if n is a palindrome in base b

%o s = digits(n,b)

%o return s == s[::-1]

%o def palQgen10(l): # unordered generator of palindromes of length <= 2*l

%o if l > 0:

%o yield 0

%o for x in range(1,10**l):

%o s = str(x)

%o yield int(s+s[-2::-1])

%o yield int(s+s[::-1])

%o A099165_list = sorted([n for n in palQgen10(6) if palQ(n,32)])

%o # _Chai Wah Wu_, Nov 25 2014

%Y Cf. A007632, A007633, A029961, A029962, A029963, A029964, A029804, A029965, A029966, A029967, A029968, A029969, A029970, A029731, A097855.

%K base,nonn

%O 1,3

%A _Robert G. Wilson v_, Sep 30 2004

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 28 11:59 EDT 2024. Contains 371254 sequences. (Running on oeis4.)