A335952 - OEIS

%I #10 May 17 2022 16:12:09

%S 1,1,1,3,3,5,3,105,3,45,105,165,45,273,7,15,231,1785,9,25935,105,3465,

%T 1155,2415,315,16575,429,1485,399,435,45,465465,2145,8415,41055,525,

%U 1485,10555545,57057,585,105,268345,19635,18723705,15015

%N Denominators of the Faulhaber polynomials.

%C See A335951 for comments, programs and references.

%o (Python)

%o from math import lcm

%o from sympy import simplify, sqrt, bernoulli

%o from sympy.abc import x

%o def A335952(n): return lcm(*(d.q for d in simplify((bernoulli(2*n,(sqrt(8*x+1)+1)/2)-bernoulli(2*n,1))/(2*n)).as_poly().all_coeffs())) # _Chai Wah Wu_, May 17 2022

%Y Cf. A335951 (numerators).

%K nonn,frac

%O 0,4

%A _Peter Luschny_, Jul 16 2020