|
%I #11 Aug 13 2019 08:13:26
%S 0,1,2,3,4,6,8,10,12,16,24,28,30,36,48,60,120
%N Numbers n such that the denominator of the Bernoulli polynomial B(n,x) equals the Clausen number C(n), {n | A144845(n) = A141056(n)}.
%C Is this a finite sequence?
%p # Clausen(n, k) defined in A160014.
%p seq(`if`(denom(bernoulli(i,x))=Clausen(i,1),i,NULL), i=0..120);
%t Clausen[n_, k_] := If[n == 0, 1, Times @@ (Select[Divisors[n], PrimeQ[# + k]&] + k)];
%t Select[Range[0, 120], Denominator[BernoulliB[#, x] // Together] == Clausen[#, 1]&] (* _Jean-François Alcover_, Aug 13 2019 *)
%Y Cf. A141056, A144845, A213621.
%K nonn,more
%O 0,3
%A _Peter Luschny_, Jun 16 2012
|
