OFFSET
1,1
LINKS
Harvey P. Dale, Table of n, a(n) for n = 1..1000
MATHEMATICA
s = 0; Do[s = s + Mod[ Prime[n], 9]; If[ PrimeQ[s], Print[s]], {n, 1, 350} ]
Select[Accumulate[Mod[Prime[Range[1000]], 9]], PrimeQ] (* Harvey P. Dale, May 04 2012 *)
PROG
(PARI) sumdigits2(n)=local(c); while (n>9, c=0; while (n>0, c=c+n%10; n=n-n%10; n=n/10); n=c); n
c=0; forprime (p=2, 1000, c=c+sumdigits2(p); if (isprime(c), print1(c, ", ")))
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Jon Perry, Oct 11 2002
EXTENSIONS
Edited by Robert G. Wilson v, Oct 14 2002
Definition corrected by M. F. Hasler, Apr 04 2007
Corrected by Harvey P. Dale, May 04 2012
STATUS
approved