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