A386728 Triangle read by rows: T(n,k) is the denominator of A(n,k), such that A(n,k) satisfies the identity for sums of odd powers: Sum_{k=1..p} k^(2n-1) = 1/(2*n) * Sum_{k=0..n-1} A(n,k) * (p^2+p)^(n-k), for all integers p >= 1.

1, 1, 6, 1, 1, 30, 1, 2, 1, 42, 1, 3, 3, 1, 30, 1, 2, 1, 2, 1, 66, 1, 1, 2, 1, 1, 1, 2730, 1, 6, 15, 3, 15, 30, 1, 6, 1, 1, 3, 3, 3, 1, 1, 1, 510, 1, 2, 1, 1, 5, 2, 5, 10, 1, 798, 1, 3, 2, 7, 1, 3, 42, 21, 21, 1, 330, 1, 2, 3, 2, 1, 6, 15, 3, 5, 10, 1, 138, 1

Offset: 0

Author

Kolosov Petro, Jul 31 2025

Keywords

• nonn
• tabl
• frac
• easy

Comments

The companion triangle with the numerators is A385567.

Extension of A093557 with k in the range 0 <= k <= n.

Examples

			Triangle begins:
  ---------------------------------------------------------
  k =   0  1   2   3    4    5    6   7    8    9    10
  ---------------------------------------------------------
  n=0:  1;
  n=1:  1, 6;
  n=2:  1, 1, 30;
  n=3:  1, 2,  1, 42;
  n=4:  1, 3,  3,  1, 30;
  n=5:  1, 2,  1,  2,  1, 66;
  n=6:  1, 1,  2,  1,  1,  1, 2730;
  n=7:  1, 6, 15,  3, 15, 30,    1,  6;
  n=8:  1, 1,  3,  3,  3,  1,    1,  1, 510;
  n=9:  1, 2,  1,  1,  5,  2,    5, 10,   1, 798;
  n=10: 1, 3,  2,  7,  1,  3,   42, 21,  21,   1, 330;
  ...
		

Links

• Donald E. Knuth, Johann Faulhaber and Sums of Powers, arXiv:9207222 [math.CA], 1992, see page 16.
• Petro Kolosov, Faulhaber's coefficients: Examples, GitHub, 2025.
• Petro Kolosov, Mathematica programs, GitHub, 2025.

Crossrefs

Cf. A385567 (numerators).

Cf. A303675, A304330, A304334, A304336.

Cf. A093558/A093559, A335951/A335952, A093556/A093557.

Programs

Mathematica

FaulhaberCoefficient[n_, k_] := 0;
FaulhaberCoefficient[n_, k_] := (-1)^(n - k) * Sum[Binomial[2 n, n - k - j]* Binomial[n - k + j, j] * (n - k - j)/(n - k + j) * BernoulliB[n + k + j], {j, 0, n - k}] /; 0 <= k < n;
FaulhaberCoefficient[n_, k_] := BernoulliB[2 n] /; k == n;
Flatten[Table[Denominator[FaulhaberCoefficient[n, k]], {n, 0, 10}, {k, 0, n}]]

PARI

T(n,k) = denominator(if (k==n, bernfrac(2*n), if (kMichel Marcus, Aug 03 2025

Formula

A(n,k) = 0 if k>n or n<0;

A(n,k) = (-1)^(n - k) * Sum_{j=0..n-k} binomial(2n, n - k - j) * binomial(n - k + j, j) * (n - k - j)/(n - k + j) * B_{n + k + j}, if 0 <= k < n;

A(n,k) = B_{2n}, if k = n;

T(n,k) = denominator(A(n,k)).