|
%I #25 Sep 08 2022 08:46:13
%S 3,5,7,19,23,29,37,41,43,47,59,61,67,73,83,89,109,113,127,131,137,139,
%T 149,151,157,163,167,173,179,181,191,193,197,199,223,227,229,239,241,
%U 257,263,269,271,281,283,293,307,311,313,317,331,337,347,353,359,373
%N Prime numbers whose sum of digits is not a power of 2.
%C Complement of A068807 with respect to primes.
%H Felix Fröhlich, <a href="/A261456/b261456.txt">Table of n, a(n) for n = 1..10000</a>
%t Select[Prime@Range[80], !IntegerQ[Log[2, Plus@@IntegerDigits[#]]] &] (* _Vincenzo Librandi_, Aug 19 2015 *)
%o (PARI) lista(nn) = forprime (n=1, nn, sd = sumdigits(n); if (2^valuation(sd,2) != sd, print1(n, ", "))); \\ _Michel Marcus_, Aug 19 2015
%o (Magma) [p: p in PrimesUpTo(400) | not PrimeDivisors(s) eq [2] where s is &+Intseq(p)]; // _Vincenzo Librandi_, Aug 19 2015
%Y Cf. A007605, A068807.
%K nonn,base
%O 1,1
%A _Altug Alkan_, Aug 19 2015
|
