A062338 - 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