OFFSET
1,2
COMMENTS
Corresponding primes are listed in A007522. - Altug Alkan, Dec 10 2015
LINKS
Harvey P. Dale, Table of n, a(n) for n = 1..1000
MAPLE
A005122:=n->`if`(isprime(8*n-1), n, NULL): seq(A005122(n), n=1..300); # Wesley Ivan Hurt, Dec 11 2015
MATHEMATICA
Select[Range[160], PrimeQ[8#-1]&] (* Harvey P. Dale, Nov 25 2011 *)
PROG
(PARI) isok(n) = isprime(8*n-1); \\ Michel Marcus, Dec 10 2015
(Magma) [n: n in [1..200] | IsPrime(8*n-1)]; // Vincenzo Librandi, Dec 10 2015
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
STATUS
approved