%I #8 Nov 16 2021 17:29:35
%S 1,1,1,1,1,2,1,1,2,1,2,2,1,1,2,2,2,1,1,1,1,1,2,2,2,2,2,1,2,2,2,1,2,2,
%T 1,1,1,1,1,1,2,2,2,1,1,2,2,1,2,1,1,1,2,2,1,1,2,1,2,1,1,1,1,1,2,1,2,1,
%U 2,2,2,1,2,1,2,2,1,2,2,1,2,1,2,2,2,2,2,2,2,2,2,1,2,1,1,1,2,1,1,2,2,1,1,1,1
%N a(n) = abs(numerator of 2n-th Bernoulli number) modulo 3.
%F a(n) = abs(A000367(n)) (mod 3)
%t Mod[Abs[Numerator[BernoulliB[2*Range[0,120]]]],3] (* _Harvey P. Dale_, Nov 16 2021 *)
%o (PARI) a(n)=abs(numerator(bernfrac(2*n)))%3
%K nonn
%O 0,6
%A _Benoit Cloitre_, Dec 13 2003