A066651 Primes of the form 2*s + 1, where s is a squarefree number (A005117).
3, 5, 7, 11, 13, 23, 29, 31, 43, 47, 53, 59, 61, 67, 71, 79, 83, 103, 107, 131, 139, 149, 157, 167, 173, 179, 191, 211, 223, 227, 229, 239, 263, 269, 277, 283, 293, 311, 317, 331, 347, 349, 359, 367, 373, 383, 389, 419, 421, 431
Offset: 1
Keywords
Examples
a(13) = A000040(18) = 61 = 2*30+1 = 2*A005117(19)+1.
Links
- Amiram Eldar, Table of n, a(n) for n = 1..10000
Programs
-
Mathematica
Select[2 * Select[Range[200], SquareFreeQ] + 1, PrimeQ] (* Amiram Eldar, Feb 22 2021 *)
-
PARI
isok(p) = isprime(p) && (p>2) && issquarefree((p-1)/2); \\ Michel Marcus, Feb 22 2021
Comments