This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
%I A076203 #23 Feb 21 2024 10:53:23 %S A076203 1,4,5,7,8,10,12,18,19,20,23,24,25,28,29,30,31,34,35,38,39,40,48,52, %T A076203 63,66,68,69,72,81,87,89,92,94,96,97,99,101,119,121,125,132,133,145, %U A076203 146,148,153,158,162,164,165,167,175,177,184,187,191 %N A076203 Numbers n such that the sum of the digits of 2^n is prime. %H A076203 Zak Seidov, <a href="/A076203/b076203.txt">Table of n, a(n) for n = 1..1000.</a> %e A076203 For n=7, 2^7=128 and 1+2+8=11, a prime. For n=23, 2^23=8388608 having the sum of its digits equaling 41. %t A076203 a={};For[n=1,n<200, n++, If[PrimeQ[Plus @@ IntegerDigits[2^n]],AppendTo[a, n]]];a (* _Stefan Steinerberger_, Jun 08 2007 *) %t A076203 Select[Range[200],PrimeQ[Total[IntegerDigits[2^#]]]&] (* _Harvey P. Dale_, Jan 30 2018 *) %o A076203 (R) %o A076203 library(gmp); digsum<-function(x) sum(as.numeric(unlist(strsplit(as.character(x),split="")))) %o A076203 which(isprime(sapply(as.character(as.bigz(2)^(1:100)),digsum))>0) #_Christian N. K. Anderson_, Apr 29 2013 %K A076203 nonn,base %O A076203 1,2 %A A076203 _Zak Seidov_, Nov 02 2002 %E A076203 More terms from _Stefan Steinerberger_, Jun 08 2007