OFFSET
1,2
COMMENTS
What is intriguing about this sequence is that the number of primes less than 10^m is of the same magnitude as A006880. Here they begin 7, 25, 122, 934.
LINKS
Z. Stankova-Frenkel and J. West, Explicit enumeration of 321,hexagon-avoiding permutations, arXiv:math/0106073 [math.CO], 2001.
EXAMPLE
5 is a term because sum of digits of first 5 primes, 2+3+5+7+(1+1)=19, is prime.
a(5) = 6 because in A051351(6) = 2 + 3 + 5 + 7 + 2 (sum of eleven's digits) + 4 (sum of thirteen's digits) which equals the sum of the digits through the sixth prime = 23 which itself is a prime.
MATHEMATICA
s = 0; Do[ s = s + Apply[ Plus, RealDigits[ Prime[ n ]] [[1]] ]; If[ PrimeQ[ s ], Print[ n ] ], {n, 0, 1000} ].
PROG
(PARI) isok(n) = isprime(sum(k=1, n, sumdigits(prime(k)))); \\ Michel Marcus, Mar 11 2017
(Python)
from sympy import isprime, nextprime
def sd(n): return sum(map(int, str(n)))
def aupto(limit):
alst, k, p, s = [], 1, 2, 2
while k <= limit:
if isprime(s): alst.append(k)
k += 1; p = nextprime(p); s += sd(p)
return alst
print(aupto(306)) # Michael S. Branicky, Jul 18 2021
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Robert G. Wilson v, Nov 18 2000
EXTENSIONS
Edited by R. J. Mathar, Aug 04 2008
STATUS
approved