OFFSET
1,2
COMMENTS
All terms are == {0,1} mod 3, because 52*(3k+2) + 1 is divisible by 3. - Zak Seidov, Oct 05 2014
LINKS
Zak Seidov, Table of n, a(n) for n = 1..2000
FORMULA
a(n) = (A142508(n) - 1)/52.
MAPLE
A248221:=n->`if`(isprime(52*n+1), n, NULL): seq(A248221(n), n=1..500); # Wesley Ivan Hurt, Oct 05 2014
MATHEMATICA
Select[Range[250], PrimeQ[52# + 1] &] (* Alonso del Arte, Oct 04 2014 *)
PROG
(PARI)
for(n=1, 10^3, if(isprime(52*n+1), print1(n, ", "))) \\ Derek Orr, Oct 05 2014
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Zak Seidov, Oct 04 2014
EXTENSIONS
More terms from Wesley Ivan Hurt, Oct 05 2014
STATUS
approved