|
%I #16 Mar 10 2022 03:21:48
%S 67,139,157,193,229,283,337,373,409,463,571,607,643,661,733,751,823,
%T 1039,1093,1129,1237,1291,1327,1381,1453,1471,1543,1723,1741,1831,
%U 2029,2083,2137,2281,2371,2551,2713,2731,2803,3019,3037,3109,3163,3181,3217,3253
%N Primes p with digital sum equal to 13.
%H Vincenzo Librandi, <a href="/A106755/b106755.txt">Table of n, a(n) for n = 1..4000</a>
%F Intersection of A000040 (primes) and A143164 (digit sum = 13); also equals { p in A000040 | A007953(p) = 13 }. - _M. F. Hasler_, Mar 09 2022
%t Select[Prime[Range[100000]], Total[IntegerDigits[#]]==13 &] (* _Vincenzo Librandi_, Jul 08 2014 *)
%o (Magma) [p: p in PrimesUpTo(10000) | &+Intseq(p) eq 13]; // _Vincenzo Librandi_, Jul 08 2014
%o (PARI) select( {is_A106755(n)=sumdigits(n)==13&&isprime(n)}, primes([1, 3333])) \\ _M. F. Hasler_, Mar 09 2022
%Y Cf. A000040 (primes), A007953 (sum of digits), A143164 (digit sum = 13).
%Y Cf. A062339 (same for digit sum s = 4), ..., A106754 (s = 11), A106756 (s = 14), and others listed in A244918 (s = 68).
%K nonn,base
%O 1,1
%A _Zak Seidov_, May 16 2005
|
