OFFSET
1,1
COMMENTS
Subsequence of A067841.
LINKS
Colin Barker, Table of n, a(n) for n = 1..350
EXAMPLE
1723 is in the sequence because 1723 % 1000 = 723, prime(1723) = 14723 and 14723 % 1000 = 723.
MATHEMATICA
Select[Prime@ Range@ 15000, Mod[#, 1000] == Mod[Prime@ #, 1000] &] (* Michael De Vlieger, Mar 29 2016 *)
PROG
(PARI) L=List(); forprime(p=2, 300000, if(p%1000==prime(p)%1000, listput(L, p))); Vec(L)
(Python)
from sympy import isprime, prime
for p in range(2, 100000):
if(prime(p)%1000==p%1000 and isprime(p)):print(p)
# Soumil Mandal, Apr 04 2016
CROSSREFS
KEYWORD
nonn,easy,base
AUTHOR
Colin Barker, Mar 29 2016
STATUS
approved