Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).
%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