A248221 Numbers m such that 52*m + 1 is prime.

1, 3, 6, 10, 13, 18, 21, 24, 25, 31, 36, 39, 40, 43, 45, 46, 49, 55, 60, 64, 66, 73, 78, 85, 91, 94, 96, 109, 115, 123, 124, 126, 129, 130, 133, 138, 139, 141, 144, 145, 151, 154, 159, 165, 168, 171, 174, 178, 181, 189, 193, 195, 196, 201, 211, 223, 225, 229

Offset: 1

Author

Zak Seidov, Oct 04 2014

Keywords

• nonn
• easy

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

Crossrefs

Cf. A033210, A142508.

Programs

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 *)

PARI

for(n=1,10^3,if(isprime(52*n+1),print1(n,", "))) \\ Derek Orr, Oct 05 2014

Formula

a(n) = (A142508(n) - 1)/52.

Extensions

More terms from Wesley Ivan Hurt, Oct 05 2014