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

 Hints (Greetings from The On-Line Encyclopedia of Integer Sequences!)
 A232188 Primes p with same last four digits as pi(p). 2
 99551, 165103, 208697, 263071, 284833, 588229, 663853, 728681, 1073819, 1181617, 1213909, 1407647, 1515259, 1590487, 1676497, 2116991, 2170681, 2202857, 2417477, 2481833, 2589143, 2664229 (list; graph; refs; listen; history; text; internal format)
 OFFSET 1,1 LINKS Robert Israel and Charles R Greathouse IV, Table of n, a(n) for n = 1..10000 (first 607 terms from Israel) FORMULA a(n) = prime(A232189(n)). EXAMPLE 15103 and prime(15103)=165103, both end with 5103. MAPLE Primes:= select(isprime, [2, seq(2*i+1, i=1..10^6)]): A232189:= select(t -> Primes[t]-t mod 10^4=0, [\$1..nops(Primes)]): Primes[A232189]; # Robert Israel, Jul 02 2015 MATHEMATICA Select[Prime[Range[200000]], Mod[#, 10000]==Mod[PrimePi[#], 10000]&] (* Harvey P. Dale, Jul 16 2021 *) PROG (PARI) p=10007; n=1230; while(n<10^6, p=nextprime(p+1); n=n+1; if(p%10^4==n%10^4, print1(p, ", "))) (MATLAB) P = primes(10^7); R = mod(P - [1:size(P, 2)], 10000); A232189 = find(R==0); P(A232189) % Robert Israel, Jul 02 2015 (PARI) n=0; forprime(p=2, 1e7, if(Mod(n++, 10000)==p, print1(p", "))) \\ Charles R Greathouse IV, Jul 05 2015 CROSSREFS Cf. A085598, A232102, A232104, A232189. Sequence in context: A138732 A108224 A268417 * A032428 A204523 A182698 Adjacent sequences: A232185 A232186 A232187 * A232189 A232190 A232191 KEYWORD nonn,base,less AUTHOR Antonio Roldán, Nov 20 2013 STATUS approved

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.

Last modified May 29 01:47 EDT 2023. Contains 363029 sequences. (Running on oeis4.)