login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A135076
Primes appearing in A001370.
1
2, 7, 5, 11, 13, 7, 19, 19, 29, 31, 41, 37, 29, 43, 41, 37, 47, 61, 59, 67, 71, 61, 73, 79, 89, 109, 103, 89, 109, 107, 107, 113, 139, 151, 127, 137, 107, 113, 167, 173, 167, 181, 191, 173, 193, 223, 233, 211, 199, 229, 251, 239, 281, 251, 277, 281, 239, 241, 239, 269
OFFSET
1,1
FORMULA
a(n)=A007953(2^A076203(n)).
EXAMPLE
a(1)=2 because with s=A076203(1)=1, 2^s=2 and sod(2)=2; sod(x)=sum of digits of x;
a(2)=7 because with s=A076203(2)=4, 2^s=16 and sod(16)=7.
a(7)=9 because with s=A076203(7)=12, 2^s=4096 and sod(4096)=19.
a(8)=9 because with s=A076203(8)=18, 2^s=262144 and sod(262144)=19.
PROG
(PARI) lista(nn) = {for (n=1, nn, sd = sumdigits(2^n); if (isprime(sd), print1(sd, ", ")); ); } \\ Michel Marcus, Oct 13 2013
CROSSREFS
Sequence in context: A220606 A232898 A319903 * A155546 A073132 A294075
KEYWORD
base,nonn
AUTHOR
Zak Seidov, Nov 18 2007
EXTENSIONS
Name corrected by Michel Marcus, Oct 13 2013
STATUS
approved