cp's OEIS Frontend

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.

A028834 Numbers whose sum of digits is a prime.

This page as a plain text file.
%I A028834 #44 Jul 05 2024 10:07:38
%S A028834 2,3,5,7,11,12,14,16,20,21,23,25,29,30,32,34,38,41,43,47,49,50,52,56,
%T A028834 58,61,65,67,70,74,76,83,85,89,92,94,98,101,102,104,106,110,111,113,
%U A028834 115,119,120,122,124,128,131,133,137,139,140,142,146,148,151,155,157,160,164,166,173,175,179,182
%N A028834 Numbers whose sum of digits is a prime.
%H A028834 Christian N. K. Anderson, <a href="/A028834/b028834.txt">Table of n, a(n) for n = 1..25000</a> (terms 1..1000 from T. D. Noe)
%H A028834 Christian N. K. Anderson, <a href="/A028834/a028834.png">Ulam Spiral</a> for a(n) <= 100000.
%e A028834 89 included because 8+9 = 17, which is prime.
%p A028834 a:=proc(n) local nn: nn:=convert(n,base,10): if isprime(sum(nn[j],j=1..nops(nn)))=true then n else fi end: seq(a(n),n=1..200); # _Emeric Deutsch_, Mar 17 2007
%t A028834 Select[Range[200],PrimeQ[Total[IntegerDigits[#]]]&]  (* _Harvey P. Dale_, Feb 18 2011 *)
%o A028834 (Haskell)
%o A028834 a028834 n = a028834_list !! (n-1)
%o A028834 a028834_list = filter ((== 1) . a010051 . a007953) [1..]
%o A028834 -- _Reinhard Zumkeller_, Nov 13 2011
%o A028834 (Sage) [x for x in range(200) if (sum(Integer(x).digits(base=10))) in Primes()] # _Bruno Berselli_, May 05 2014
%o A028834 (PARI) is(n)=isprime(sumdigits(n)) \\ _Felix Fröhlich_, Aug 16 2014
%o A028834 (Python)
%o A028834 from sympy import isprime
%o A028834 def ok(n): return isprime(sum(map(int, str(n))))
%o A028834 print(list(filter(ok, range(183)))) # _Michael S. Branicky_, Jun 18 2021
%o A028834 (R) require(gmp); which(sapply(1:1000, function(i) isprime(sum(floor(i/10^(0:(nchar(i)-1)))%%10)))==2) # _Christian N. K. Anderson_, Apr 22 2024
%Y A028834 Cf. A007953, A028835, A028842, A052018, A052019, A052020, A052021, A052022.
%Y A028834 Cf. A010051; A046704 is a subsequence.
%Y A028834 Complement of A104211.
%K A028834 nonn,base,easy,nice
%O A028834 1,1
%A A028834 Armand Turpel (armand(AT)vo.lu, armand_t(AT)geocities.com)
%E A028834 More terms from Scott Lindhurst (ScottL(AT)alumni.princeton.edu)