OFFSET
1,1
LINKS
Indranil Ghosh, Table of n, a(n) for n = 1..6414
FORMULA
MAPLE
isA077510 := proc(n)
isprime(n+numtheory[pi](n)) ;
end proc:
A077510 := proc(n)
local a;
if n = 1 then
return 2;
else
for a from procname(n-1)+1 do
if isA077510(a) then
return a;
end if;
end do:
end if:
end proc:
A076757 := proc(n)
local a10 ;
a10 := A077510(n) ;
a10+numtheory[pi](a10) ;
end proc:
seq(A076757(n), n=1..40) ; # R. J. Mathar, Nov 19 2011
MATHEMATICA
Select[Table[n + PrimePi[n], {n, 500}], PrimeQ] (* T. D. Noe, Nov 19 2011 *)
PROG
(Magma) [a: n in [1..400] | IsPrime(a) where a is (n + #PrimesUpTo(n))]; // Vincenzo Librandi, Jan 29 2017
CROSSREFS
KEYWORD
nonn
AUTHOR
David Garber, Nov 13 2002
EXTENSIONS
Name edited by Michel Marcus, Dec 30 2013
STATUS
approved