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

Author

Zak Seidov, May 16 2005

Keywords

• nonn
• base

Links

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

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).

Programs

Magma

[p: p in PrimesUpTo(4000) | &+Intseq(p) eq 16]; // Vincenzo Librandi, Jul 08 2014

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 *)

PARI

select( {is_A106757(n)= sumdigits(n)==16 && isprime(n)}, primes([1, 3333])) \\ M. F. Hasler, Mar 09 2022

Formula

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

Extensions

More terms from Zak Seidov, Oct 30 2009