OFFSET
1,1
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 1..1000
FORMULA
a(n)=8*A023228(n)+1. - R. J. Mathar, Mar 15 2009
a(n) >> n log^2 n. - Charles R Greathouse IV, Oct 21 2021
MATHEMATICA
Flatten[Table[If[PrimeQ[n] && PrimeQ[(n - 1)/8], n, {}], {n, 1, 10000}]]
Select[Prime[Range[1500]], PrimeQ[(# - 1) / 8]&] (* Vincenzo Librandi, Apr 14 2013 *)
PROG
(PARI) list(lim)=my(v=List()); forprime(p=2, (lim-1)\8, if(isprime(8*p+1), listput(v, 8*p+1))); Vec(v) \\ Charles R Greathouse IV, Oct 20 2021
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Roger L. Bagula, Mar 11 2009
EXTENSIONS
Edited by the Associate Editors of the OEIS, Apr 22 2009
STATUS
approved