login

Reminder: The OEIS is hiring a new managing editor, and the application deadline is January 26.

Safe primes p such that p + 24 is also a safe prime.
1

%I #40 Dec 31 2019 19:19:19

%S 23,59,83,359,479,563,839,863,1283,2039,2879,2999,3779,4259,4679,5483,

%T 7703,10079,12203,13103,23603,26903,27803,30323,31583,33623,35339,

%U 41519,43403,44519,44939,53759,55079,57119,57899,58043,65123,66359,70139,70199,76379,77723,79943

%N Safe primes p such that p + 24 is also a safe prime.

%H John Cerkan, <a href="/A274381/b274381.txt">Table of n, a(n) for n = 1..10000</a>

%F { x | both x and x + 24 are safe primes }.

%e 83 is a safe prime and so is 83 + 24 = 107.

%e 4679 is a safe prime and so is 4679 + 24 = 4703.

%p A274381:=n->`if`(type((n-1)/2, prime) and isprime(n+24) and isprime(n) and type((n+23)/2, prime), n, NULL): seq(A274381(n), n=1..10^5); # _Wesley Ivan Hurt_, Jun 25 2016

%t Select[Prime@ Range[10^4], PrimeQ[(# - 1)/2] && PrimeQ[# + 24] && PrimeQ[(23 + #)/2] &] (* _Giovanni Resta_, Jun 23 2016 *)

%t Select[Prime[Range[8000]],AllTrue[{(#-1)/2,(#+23)/2,#+24},PrimeQ]&] (* The program uses the AllTrue function from Mathematica version 10 *) (* _Harvey P. Dale_, Dec 31 2019 *)

%o (PARI) lista(nn) = forprime(p=3, nn, if (isprime((p-1)/2) && isprime(q=p+24) && isprime((q-1)/2), print1(p, ", "))); \\ _Michel Marcus_, Jun 23 2016

%Y Cf. A005385, A033560, A059323.

%K nonn

%O 1,1

%A _Debapriyay Mukhopadhyay_, Jun 23 2016