OFFSET
1,1
LINKS
Charles R Greathouse IV, Table of n, a(n) for n = 1..10000
PROG
(PARI) lista(nn) = forprime(p=2, nn, if ((sd=sumdigits(p)) == sumdigits(primepi(p)), print1(sd, ", ")); ); \\ Michel Marcus, Jun 18 2015
(PARI) go(lim)=my(v=List(), n, s); forprime(p=2, lim, s=sumdigits(n++); if(sumdigits(p)==s, listput(v, s))); Vec(v) \\ Charles R Greathouse IV, Jun 18 2015
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Zak Seidov, Jun 18 2015
STATUS
approved