OFFSET
1,1
COMMENTS
Except for (5,7) and (7,11), all terms are divisible by 12, since safe primes are congruent to 5 modulo 6 except 7 and safe_prime + 6 is not a safe prime. Closest safe primes differ by 12 like (11,23) or (83207,83219).
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000 (terms 1..5001 from G. C. Greubel)
FORMULA
a(n) = A005385(n+1) - A005385(n). [corrected by Harvey P. Dale, Jan 12 2011 and Zak Seidov, Sep 19 2016]
MATHEMATICA
Differences[Select[Prime[Range[500]], PrimeQ[(#-1)/2]&]] (* Harvey P. Dale, Jan 12 2011 *)
PROG
(PARI) list(lim) = {my(p1 = 5); forprime(p2 = 7, lim, if(isprime((p2-1)/2), print1(p2-p1, ", "); p1 = p2)); } \\ Amiram Eldar, Mar 02 2025
CROSSREFS
KEYWORD
nonn,changed
AUTHOR
Labos Elemer, Jan 26 2001
EXTENSIONS
Offset corrected by Amiram Eldar, Mar 02 2025
STATUS
approved