Reminder: The OEIS is hiring a new managing editor, and the application deadline is January 26.
%I #22 Dec 12 2023 13:44:34
%S 1,-1,1,2,-2,1,0,2,-3,1,-8,0,4,-4,1,0,-8,0,20,-5,1,32,0,-8,0,10,-6,1,
%T 0,32,0,-56,0,14,-7,1,-128,0,128,0,-112,0,56,-8,1,0,-384,0,128,0,-336,
%U 0,24,-9,1,2560,0,-384,0,320,0,-112,0,30,-10,1
%N Triangle read by rows: a(n,m) = numerator(binomial(n,m)*2^(n-m)*B(n-m)) with B(k) the Bernoulli numbers A027641(k)/A027642(k).
%C The denominator triangle b(n,m) is given in A285865.
%C a(n,m)/b(n,m) = B(2;n,m) is the d = 2 instance of the fractional d-family of triangles B(d;n,m) = binomial(n,m)*d^(n-m)*B(n-m), for d >= 1. They are the coefficient triangles of generalized Bernoulli polynomials PB(d;n,x) = Sum_{m=0..n} B(d;n,m)*x^m for n >= 0.
%C {PB(d;n,x)}_{n>=0} has e.g.f. EB(d;x,z) := Sum_{n>=0} PB(d;n,x)*z^n = d*z*exp(x*z)/(exp(d*z)-1). B(d;n,m) is a Sheffer triangle of the Appell type for each d, denoted by (d*z/(exp(d*z - 1)), z).
%C PB(d;n,x) gives a (trivial) generalization of the Bernoulli polynomials with coefficients given in A196838/A196839 (rising powers of x), and this is PB(1;n,x).
%C The polynomials PB(d;n,x) appear in the generalized Faulhaber formula for sums of powers of arithmetic progressions SP(n,m) := Sum_{j=0..m} (a + d*j)^n, n >= 0, m >= 0, d >= 1, a = 0 for d = 1 and a from the smallest positive restricted residue system modulo d >= 2. For this Faulhaber formula see a comment in A285863, where they are named B(d;n,x).
%C The row sums of the rational triangle B(2;n,m) give A157779(n)/A141459(n). The alternating row sums are given in A285866/A141459(n).
%F a(n,m) = numerator(binomial(n, m)*2^(n-m)*B(n-m)), with the Bernoulli numbers B(k) = A027641(k)/A027642(k).
%F E.g.f.s of the rational column sequences {B(2;n, m)}_{n>=0} are Ecol(m, x) = (2*x/(exp(2*x) - 1))*x^m/m! (Sheffer property). Here the numerators of column m are numerator([x^m/m!] Ecol(m, x)), m >= 0.
%e The triangle a(n,m) begins:
%e n\m 0 1 2 3 4 5 6 7 8 9 10 ...
%e 0: 1
%e 1: -1 1
%e 2: 2 -2 1
%e 3: 0 2 -3 1
%e 4: -8 0 4 -4 1
%e 5: 0 -8 0 20 -5 1
%e 6: 32 0 -8 0 10 -6 1
%e 7: 0 32 0 -56 0 14 -7 1
%e 8: -128 0 128 0 -112 0 56 -8 1
%e 9: 0 -384 0 128 0 -336 0 24 -9 1
%e 10: 2560 0 -384 0 320 0 -112 0 30 -10 1
%e ...
%e The rational triangle B(2;n,m) = a(n,m)/A285865(n,m) begins:
%e n\m 0 1 2 3 4 5 6 7 8 9 10 ...
%e 0: 1
%e 1: -1 1
%e 2: 2/3 -2 1
%e 3: 0 2 -3 1
%e 4: -8/15 0 4 -4 1
%e 5: 0 -8/3 0 20/3 -5 1
%e 6: 32/21 0 -8 0 10 -6 1
%e 7: 0 32/3 0 -56/3 0 14 -7 1
%e 8: -128/15 0 128/3 0 -112/3 0 56/3 -8 1
%e 9: 0 -384/5 0 128 0 -336/5 0 24 -9 1
%e 10: 2560/33 0 -384 0 320 0 -112 0 30 -10 1
%e ...
%p T := d -> (n,m) -> numer(binomial(n, m)*d^(n-m)*bernoulli(n-m)):
%p for n from 0 to 10 do seq(T(2)(n,k),k=0..n) od; # _Peter Luschny_, May 04 2017
%t T[n_, m_]:=Numerator[Binomial[n, m]*2^(n - m)*BernoulliB[n - m]]; Table[T[n, m], {n, 0, 20}, {m, 0, n}] // Flatten (* _Indranil Ghosh_, May 06 2017 *)
%o (PARI) T(n, m) = numerator(binomial(n, m)*2^(n - m)*bernfrac(n - m));
%o for(n=0, 20, for(m=0, n, print1(T(n, m),", ");); print();) \\ _Indranil Ghosh_, May 06 2017
%o (Python)
%o from sympy import binomial, bernoulli
%o def T(n, m): return (binomial(n, m) * 2**(n - m) * bernoulli(n - m)).numerator()
%o for n in range(21): print([T(n, m) for m in range(n + 1)]) # _Indranil Ghosh_, May 06 2017
%Y Cf. A027641/A027642, A157779/A141459, A196838/A196839, A285863, A285865.
%K sign,easy,tabl,frac
%O 0,4
%A _Wolfdieter Lang_, May 03 2017