login
A153209
Primes of the form 2*p+1 where p is prime and p+1 is squarefree.
7
5, 11, 59, 83, 227, 347, 563, 1019, 1283, 1307, 1523, 2459, 2579, 2819, 2963, 3803, 3947, 4259, 4547, 5387, 5483, 6779, 6827, 7187, 8147, 9587, 10667, 10883, 11003, 12107, 12227, 12539, 12659, 13043, 13163, 14243, 14387, 15683, 16139, 16187
OFFSET
1,1
COMMENTS
Subsequence of A005385.
LINKS
EXAMPLE
For p = 2 (the only case with p+1 odd), 2*p+1 = 5 is prime and p+1 = 3 is squarefree, so 5 is in the sequence. For p = 3, 2*p+1 = 7 is prime and p+1 = 4 is not squarefree, so 7 is not in the sequence.
MATHEMATICA
lst = {}; Do[p = Prime[n]; If[PrimeQ[Floor[p/2]] && SquareFreeQ[Ceiling[p/2]], AppendTo[lst, p]], {n, 7!}]; lst
Select[2#+1&/@Select[Prime[Range[2000]], SquareFreeQ[#+1]&], PrimeQ] (* Harvey P. Dale, Aug 02 2024 *)
PROG
(Magma) [ q: p in PrimesUpTo(8100) | IsSquarefree(p+1) and IsPrime(q) where q is 2*p+1 ];
CROSSREFS
Cf. A005117 (squarefree numbers), A005385 (safe primes p: (p-1)/2 is also prime), A153207, A153208, A153210.
Sequence in context: A121934 A153812 A269454 * A239026 A106257 A104358
KEYWORD
nonn
AUTHOR
EXTENSIONS
Edited by Klaus Brockhaus, Dec 24 2008
Mathematica updated by Jean-François Alcover, Jul 04 2013
STATUS
approved