login

Reminder: The OEIS is hiring a new managing editor, and the application deadline is January 26.

Sum of digits of Honaker primes (A033548).
1

%I #11 Jun 20 2015 14:48:27

%S 5,11,16,13,14,11,5,11,11,14,14,16,8,7,14,11,17,17,16,19,20,13,16,14,

%T 16,13,22,13,14,13,19,23,14,16,16,14,17,13,14,17,14,20,23,10,11,16,17,

%U 19,20,20,23,23,23,11,13,20,13,20,17,19,10,19,13,14,16,20,10,10,13,10,13,10,13,13,19,17,13,11,14,14,14,22,16,19,20,16,20,19,20,19,20

%N Sum of digits of Honaker primes (A033548).

%H Charles R Greathouse IV, <a href="/A259067/b259067.txt">Table of n, a(n) for n = 1..10000</a>

%F a(n) = A007953(A033548(n)) = A007953(A033549(n)).

%o (PARI) lista(nn) = forprime(p=2, nn, if ((sd=sumdigits(p)) == sumdigits(primepi(p)), print1(sd, ", "));); \\ _Michel Marcus_, Jun 18 2015

%o (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

%Y Cf. A007605, A007953, A033548, A033549.

%K nonn,base

%O 1,1

%A _Zak Seidov_, Jun 18 2015