A068807 Primes whose sum of digits is a power of 2.
2, 11, 13, 17, 31, 53, 71, 79, 97, 101, 103, 107, 211, 233, 251, 277, 349, 367, 431, 439, 457, 503, 521, 547, 619, 673, 691, 701, 709, 727, 853, 907, 1021, 1061, 1069, 1087, 1151, 1201, 1223, 1249, 1429, 1447, 1483, 1511, 1601, 1609, 1627, 1663, 1753, 1861, 1933, 1951, 2011, 2141
Offset: 1
Links
- Aaron Toponce, Table of n, a(n) for n = 1..10000
Crossrefs
Different from A178796.
Programs
-
Magma
[p: p in PrimesUpTo(2200) | PrimeDivisors(s) eq [2] where s is &+Intseq(p)]; // Bruno Berselli, Dec 26 2012
-
Mathematica
Select[ Prime@Range[ 300 ],IntegerQ[ Log[ 2,Plus@@IntegerDigits[ # ] ] ]& ] (* Ray Chandler, Dec 03 2009 *)
-
PARI
lista(nn) = forprime (n=1, nn, sd = sumdigits(n); if (2^valuation(sd,2) == sd, print1(n, ", "))); \\ Michel Marcus, Aug 19 2015
Extensions
Extended by Ray Chandler, Dec 03 2009
Comments