OFFSET
1,1
LINKS
Robert Israel, Table of n, a(n) for n = 1..5839
FORMULA
a(n) = prime(A067841(n)).
EXAMPLE
1723, and prime(1723)= 14723, both end with 723.
MAPLE
P:= select(isprime, [2, seq(i, i=3..10^6, 2)]):
P[select(t -> P[t]-t mod 1000 = 0, [seq(i, i=3..nops(P), 2)])]; # Robert Israel, Oct 03 2017
MATHEMATICA
Select[Prime[Range[34000]], Mod[#, 1000]==Mod[PrimePi[#], 1000]&] (* Harvey P. Dale, Dec 30 2017 *)
PROG
(PARI)
cutdigit(a, p, q)=(a%10^q)\10^(p-1)
{for(n=1, 40000, p=prime(n); if(cutdigit(p, 1, 3)==cutdigit(n, 1, 3), print(p)))}
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Antonio Roldán, Nov 18 2013
STATUS
approved