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
KEYWORD
nonn,base,less
AUTHOR
Antonio Roldán, Nov 20 2013
STATUS
approved