|
%I #21 May 13 2013 01:54:23
%S 6,42,330,510,690,798,870,1410,1518,1590,1770,1806,2490,3102,3210,
%T 3318,3894,4110,4326,4470,4686,5010,5190,5370,5478,6486,6810,7062,
%U 7890,8070,8142,8646,8790,9366,9510,10410,10770,11022
%N Denominators of Bernoulli numbers which are == 6 (mod 9).
%C The sequence contains the elements of A090801 which are == 6 (mod 9).
%C Conjecture: all first differences 36, 288, 180, 180,... of the sequence are multiples of 36.
%C The conjecture is true, since elements of A090801 are 2 mod 4. - _Charles R Greathouse IV_, Nov 22 2012
%H Charles R Greathouse IV, <a href="/A218755/b218755.txt">Table of n, a(n) for n = 1..10000</a>
%t Take[Union[Select[Denominator[BernoulliB[Range[1000]]],Mod[#,9]==6&]],60] (* _Harvey P. Dale_, Nov 28 2012 *)
%o (PARI) is(n)=if(n%36-6, 0, my(f=factor(n)); if(vecmax(f[, 2])>1, return(0)); fordiv(lcm(apply(k->k-1, f[, 1])), k, if(isprime(k+1) && n%(k+1), return(0))); 1) \\ _Charles R Greathouse IV_, Nov 26 2012
%Y Second subset of the Bernoulli denominators: A090801 which are == 3 (mod 9).
%K nonn,easy
%O 1,1
%A _Paul Curtz_, Nov 05 2012
|
