OFFSET
1,1
COMMENTS
All terms are odd. - Karl-Heinz Hofmann, Sep 15 2025
LINKS
Paolo Xausa, Table of n, a(n) for n = 1..5000
MATHEMATICA
With[{upto=10^4}, Flatten[Position[Accumulate[Prime[Range[upto]]], _?(PrimeQ[#-1]&&PrimeQ[#+1]&)]]]
PROG
(PARI) isok(k) = my(s=vecsum(primes(k))); isprime(s-1) && isprime(s+1); \\ Michel Marcus, Oct 04 2023
(Python)
from sympy import isprime, sieve
print([k for k in range(1, 10402, 2) if isprime((s:=sum(sieve[1:k+1]))-1) and isprime(s+1)])
# Karl-Heinz Hofmann, Sep 15 2025
CROSSREFS
KEYWORD
nonn
AUTHOR
Paolo Xausa, Oct 04 2023
STATUS
approved
