This site is supported by donations to The OEIS Foundation.

 Hints (Greetings from The On-Line Encyclopedia of Integer Sequences!)
 A103402 Palindromes p such that pi(p) is a palindromic prime. 5
 3, 4, 5, 6, 11, 33, 555, 878, 5775, 6116, 919919, 58633685, 129707921, 16958285961, 867275572768 (list; graph; refs; listen; history; text; internal format)
 OFFSET 1,1 COMMENTS From a suggestion from Zak Seidov, Feb 02 2005. a(16) > 32*10^12. [From Donovan Johnson, Dec 03 2009] LINKS MATHEMATICA 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]]]] ]]]]; p = 0; a = {}; Do[p = NextPalindrome[p]; q = PrimePi[p]; If[PrimeQ[q], r = IntegerDigits[q]; If[Reverse[r] == r, Print[{p, q}]; AppendTo[a, p]]], {n, 10^6}]; a palQ[n_] := Reverse[x = IntegerDigits[n]] == x; t = {}; Do[If[palQ[n] && PrimeQ[x = PrimePi[n]] && palQ[x], AppendTo[t, n]], {n, 10^6}]; t (* Jayanta Basu, Jun 24 2013 *) CROSSREFS Cf. A046941, A046942, A103357, A103358, A103403. Sequence in context: A047425 A048989 A299299 * A154664 A191750 A034448 Adjacent sequences:  A103399 A103400 A103401 * A103403 A103404 A103405 KEYWORD base,nonn AUTHOR Robert G. Wilson v, Feb 03 2005 EXTENSIONS a(15) from Donovan Johnson, Dec 03 2009 STATUS approved

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

Last modified January 23 21:17 EST 2019. Contains 319404 sequences. (Running on oeis4.)