|
%I #25 Sep 08 2022 08:45:03
%S 5,19,23,37,41,73,109,113,127,131,163,181,271,307,311,389,401,433,479,
%T 523,541,569,587,613,631,659,677,811,839,857,929,947,983,997,1009,
%U 1013,1031,1063,1103,1117,1153,1171,1289,1301,1423,1487,1531,1559,1621,1667
%N Primes whose sum of digits is a multiple of 5.
%H Harvey P. Dale, <a href="/A062340/b062340.txt">Table of n, a(n) for n = 1..1000</a>
%F Intersection of A000040 (primes) and A227793 (sum of digits in 5Z). - _M. F. Hasler_, Mar 10 2022
%e 569 is a prime with sum of digits = 20, hence belongs to the sequence.
%t Select[Prime[Range[300]],Divisible[Total[IntegerDigits[#]],5]&] (* _Harvey P. Dale_, Jul 06 2020 *)
%o (Magma) [ p: p in PrimesUpTo(10000) | &+Intseq(p) mod 5 eq 0 ]; // _Vincenzo Librandi_, Apr 02 2011
%o (Python)
%o from sympy import primerange as primes
%o def ok(p): return sum(map(int, str(p)))%5 == 0
%o print(list(filter(ok, primes(1, 1668)))) # _Michael S. Branicky_, May 19 2021
%o (PARI) select( {is_A062340(n)=sumdigits(n)%5==0&&isprime(n)}, primes([1,2000])) \\ _M. F. Hasler_, Mar 10 2022
%Y Cf. A007953 (sum of digits), A227793 (sum of digits divisible by 5).
%Y Has as subsequence A062341 (primes with sum of digits s = 5), A107579 (s = 10), A106760 (s = 20), A106763 (s = 25), A106770 (s = 35), A106773 (s = 40), A106780 (s = 50), A106783 (s = 55), A107619 (s = 65) and A181321 (s = 70).
%Y Cf. A062340 (equivalent for 8).
%K nonn,base,easy
%O 1,1
%A _Amarnath Murthy_, Jun 21 2001
%E Corrected and extended by _Harvey P. Dale_ and Larry Reeves (larryr(AT)acm.org), Jul 04 2001
|
