A106757 Primes with digit sum = 16.

79, 97, 277, 349, 367, 439, 457, 547, 619, 673, 691, 709, 727, 853, 907, 1069, 1087, 1249, 1429, 1447, 1483, 1609, 1627, 1663, 1753, 1861, 1933, 1951, 2239, 2293, 2347, 2383, 2437, 2473, 2617, 2671, 2707, 2833, 2851, 3049, 3067, 3229, 3319, 3373, 3391

Offset

1,1

Links

Harvey P. Dale, Table of n, a(n) for n = 1..1000

Formula

Intersection of A000040 (primes) and A235227 (digit sum = 16); also equals { p in A000040 | A007953(p) = 16 }. - M. F. Hasler, Mar 09 2022

Mathematica

Reap[Do[If[16==Apply[Plus, IntegerDigits[p=Prime[n]]], Sow[p]], {n, 1000}]][[2, 1]] (* Zak Seidov, Oct 30 2009 *)
Select[Prime[Range[500]], Total[IntegerDigits[#]]==16&] (* Harvey P. Dale, Nov 14 2011 *)

Prog

(Magma) [p: p in PrimesUpTo(4000) | &+Intseq(p) eq 16]; // Vincenzo Librandi, Jul 08 2014
(PARI) select( {is_A106757(n)= sumdigits(n)==16 && isprime(n)}, primes([1, 3333])) \\ M. F. Hasler, Mar 09 2022

Crossrefs

Cf. A000040 (primes), A007953 (sum of digits), A235227 (digit sum = 16).

Cf. A062339 (same for digit sum s = 4), A106756 (s = 14), A106758 (s = 17), and others listed in A244918 (s = 68).

Subsequence of A062342 (primes whose sum of digits is a multiple of 8) and of A107288 (primes with sum of digits a square).

Sequence in context: A043259 A044039 A065211 * A285228 A087537 A117840

Adjacent sequences: A106754 A106755 A106756 * A106758 A106759 A106760

Keyword

nonn,base

Author

Zak Seidov, May 16 2005

Extensions

More terms from Zak Seidov, Oct 30 2009

Status

approved