|
%I #17 Sep 08 2022 08:45:03
%S 13,17,31,53,71,79,97,103,107,211,233,251,277,349,367,389,431,439,457,
%T 479,503,521,547,569,587,619,659,673,677,691,701,709,727,839,853,857,
%U 907,929,947,983,1021,1061,1069,1087,1151,1201,1223,1249,1289,1429
%N Primes whose sum of digits is a multiple of 4.
%H Robert Israel, <a href="/A062338/b062338.txt">Table of n, a(n) for n = 1..10000</a>
%e 349 is a prime with sum of digits =16=4*4, hence belongs to the sequence.
%p filter:= x -> (convert(convert(x,base,10),`+`) mod 4 = 0) and isprime(x);
%p A062338:= select(filter, [seq(2*i+1,i=0..1000)]); # _Robert Israel_, Apr 20 2014
%o (Magma) [ p: p in PrimesUpTo(10000) | &+Intseq(p) mod 4 eq 0 ]; /* _Vincenzo Librandi_, Apr 02 2011 */
%o (PARI) is(n)=isprime(n) && sumdigits(n)%4==0 \\ _Charles R Greathouse IV_, Mar 09 2022
%Y Subsequence of A119449.
%Y Cf. A062336, A062337.
%K nonn,base,easy
%O 1,1
%A _Amarnath Murthy_, Jun 21 2001
%E Corrected and extended by Larry Reeves (larryr(AT)acm.org), Jul 06 2001
|
