login
A076757
Primes of the form n + pi(n), that is, generated in A077510.
5
3, 5, 11, 13, 17, 19, 29, 37, 43, 47, 59, 67, 71, 73, 79, 83, 89, 97, 101, 103, 107, 109, 113, 127, 131, 139, 149, 151, 163, 167, 173, 179, 181, 191, 197, 199, 211, 223, 227, 229, 239, 251, 263, 269, 271, 277, 281, 307, 313, 317, 331, 337, 347
OFFSET
1,1
LINKS
FORMULA
a(n) = k+A000720(k) where k=A077510(n). - R. J. Mathar, Nov 19 2011
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
Sequence in context: A059349 A045316 A040100 * A178942 A045404 A154500
KEYWORD
nonn
AUTHOR
David Garber, Nov 13 2002
EXTENSIONS
Name edited by Michel Marcus, Dec 30 2013
STATUS
approved