OFFSET
1,2
COMMENTS
Also: 1 plus numbers k such that prime(k) - 6 is prime. - Michael S. Branicky, Aug 16 2025
LINKS
Charles R Greathouse IV, Table of n, a(n) for n = 1..10000
EXAMPLE
a(1) = 1 because prime(1)-5/2-7/2 = 2-5/2-7/2 = -4 (nonprime) and prime(1)-5/2+7/2 = 2-5/2+7/2 = 3 (prime).
MATHEMATICA
Join[{1}, Select[Range[5, 150], PrimeQ[Prime[#] - 6] &]] (* Paolo Xausa, Aug 17 2025 *)
PROG
(PARI) is(n, p=prime(n))=isprime(p-6) || isprime(p+1) \\ Charles R Greathouse IV, Jun 13 2017
(Python)
from sympy import isprime, prime
def ok(n): return n == 1 or isprime(prime(n)-6)
print([k for k in range(1, 145) if ok(k)]) # Michael S. Branicky, Aug 16 2025
CROSSREFS
KEYWORD
nonn
AUTHOR
Juri-Stepan Gerasimov, Jan 25 2010
EXTENSIONS
Corrected (57 inserted) by R. J. Mathar, May 02 2010
STATUS
approved
