A285066 -id:A285066 - cp's OEIS frontend

A285061 Sheffer triangle S2[4,1] = (exp(x), exp(4*x) - 1).

1, 1, 4, 1, 24, 16, 1, 124, 240, 64, 1, 624, 2656, 1792, 256, 1, 3124, 26400, 33920, 11520, 1024, 1, 15624, 250096, 546560, 331520, 67584, 4096, 1, 78124, 2313360, 8105664, 7822080, 2745344, 372736, 16384, 1, 390624, 21132736, 114627072, 165398016, 88940544, 20299776, 1966080, 65536, 1, 1953124, 191757120, 1574682880, 3270274560, 2529343488, 863256576, 138215424, 10027008, 262144

Offset: 0

Wolfdieter Lang, Apr 13 2017

For Sheffer triangles (infinite lower triangular exponential convolution matrices) see the W. Lang link under A006232, with references.
This is a generalization of the Sheffer triangle Stirling2(n, m) = A048993(n, m) denoted by (exp(x), exp(x)-1), which could be named S2[1,0].
For the Sheffer triangle (exp(x), (1/4)*(exp(4*x) - 1)) see A111578, also the P. Bala link where this triangle, T(n, m)/4^m, is named S_{(4,0,1)}. The triangle T(n, m)*m! is given in A285066.
The a-sequence for this Sheffer triangle has e.g.f. 4*x/log(1+x) and is 4*A006232(n)/A006233(n) (Cauchy numbers of the first kind).
The z-sequence has e.g.f. (4/(log(1+x)))*(1 - 1/(1+x)^(1/4)) and is A285062(n)/A285063(n).
The main diagonal gives A000302.
The row sums give A285064. The alternating row sums give A285065.
The first column sequences are A000012, 4*A003463, 4^2*A016234. For the e.g.f.s and o.g.f.s see the formula section.
This triangle appears in the o.g.f. G(n, x) of the sequence {(1 + 4*m)^n}{m>=0}, as G(n, x) = Sum{m=0..n} T(n, m)*m!*x^m/(1-x)^(m+1), n >= 0. Hence the corresponding e.g.f. is, by the linear inverse Laplace transform, E(n, t) = Sum_{m >=0}(1 + 4*m)^n t^m/m! = exp(t)*Sum_{m=0..n} T(n, m)*t^m.
The corresponding Euler triangle with reversed rows is rEu(n, k) = Sum_{m=0..k} (-1)^(k-m)*binomial(n-m, k-m)*T(n, k)*k!, 0 <= k <= n. This is A225118 with row reversion.
In general the Sheffer triangle S2[d,a] appears in the reordering of the operator (a*1 + d*E_x) = Sum_{m=0..n} S2(d,a;n,m) x^m*(d_x)^m with the derivative d_x and the Euler operator E_x = x*d_x. For [d,a] = [1,0] this becomes the standard Stirling2 property.

			The triangle T(n,m) begins:
  n\m  0      1        2         3         4        6        7       8     9
  0:   1
  1:   1      4
  2:   1     24       16
  3:   1    124      240        64
  4:   1    624     2656      1792       256
  5:   1   3124    26400     33920     11520     1024
  6:   1  15624   250096    546560    331520    67584     4096
  7:   1  78124  2313360   8105664   7822080  2745344   372736   16384
  8:   1 390624 21132736 114627072 165398016 88940544 20299776 1966080 65536
  ...
Three term recurrence: T(4, 1) = 4*1 + (1 + 4*1)*124 = 624.
Recurrence for row polynomial R(3, x) (Meixner type): ((1 + 4*x) + 4*x*d_x)*(1 + 24*x + 16*x^2) = 1 + 124*x + 240*x^2 + 64*x^3.
Boas-Buck recurrence for column m = 2, and n = 4: T(4, 2) = (1/2)*(2*(2 + 4*2)*T(3, 2) + 2*6*(-4)^2*Bernoulli(2)*T(2, 2)) = (1/2)*(20*240 + 12*16*(1/6)*16) = 2656. - _Wolfdieter Lang_, Aug 11 2017

Michael De Vlieger, Table of n, a(n) for n = 0..11475, rows n = 0..150, flattened.
Paweł Hitczenko, A class of polynomial recurrences resulting in (n/log n, n/log^2 n)-asymptotic normality, arXiv:2403.03422 [math.CO], 2024. See p. 9.
Wolfdieter Lang, On Sums of Powers of Arithmetic Progressions, and Generalized Stirling, Eulerian and Bernoulli numbers, arXiv:1707.04451 [math.NT], 2017.

Cf. A000012, A000302, 4*A003463, A111578, A006232/A006233, 16*A016234, A282629, A285062/A285063, A285064, A285065, A285066.

Table[Sum[Binomial[n, k]*4^k*StirlingS2[k, m], {k, 0, n}], {n, 0, 20}, {m, 0, n}] // Flatten (* Indranil Ghosh, May 06 2017 *)

T(n, m) = sum(k=0, n, binomial(n, k)*4^k*stirling(k, m, 2));
for(n=0, 20, for(m=0, n, print1(T(n, m),", ");); print();) \\ Indranil Ghosh, May 06 2017

Three term recurrence (from the conversion property mentioned above, with [d,a] =[4,1]): T(n, -1) = 0, T(0, 0) = 1, T(n, m) = 0 if n < m. T(n, m) = 4*T(n-1, m-1) + (1 + 4*m)*T(n-1, m) for n >= 1, m = 0..n.
T(n, m) = Sum_{k=0..m} binomial(m,k)*(-1)^(k-m)*(1 + 4*k)^n/m!, 0 <= m <= n, satisfying the recurrence.
E.g.f. of the row polynomials R(n, x) = Sum_{m=0..n} T(n, m)*x^m: exp(z)*exp(x*(exp(4*z) - 1)). This is the e.g.f. of the triangle.
E.g.f. for the sequence of column m: exp(x)*((exp(3*x) - 1)^m)/m! (Sheffer property).
O.g.f. for sequence of column m: (4*x)^m/Product_{j=0..m} (1 - (1 + 4*j)*x) (by Laplace transform of the e.g.f.).
T(n, m) = Sum_{k=0..n} binomial(n, k)* 4^k*Stirling2(k, m), 0 <= m <= n, where Stirling2 is given in A048993.
A nontrivial recurrence for the column m=0 entries T(n, 0) = 4^n from the z-sequence given above: T(n,0) = n*Sum_{j=0..n-1} z(j)*T(n-1,j), n >= 1, T(0, 0) = 1.
Recurrence for column m >= 1 entries from the a-sequence given above: T(n, m) = (n/m)* Sum_{j=0..n-m} binomial(m-1+j, m-1)*a(j)*T(n-1, m-1+j), m >= 1.
Recurrence for row polynomials R(n, x) (Meixner type): R(n, x) = ((1 + 4*x) + 4*x*d_x)*R(n-1, x), with differentiation d_x, for n >= 1, with input R(0, x) = 1.
Boas-Buck recurrence for column sequence m: T(n, m) = (1/(n - m))*((n/2)*(2 + 4*m)*T(n-1, m) + m*Sum_{p=m..n-2} binomial(n, p)*(-4)^(n-p)*Bernoulli(n-p)*T(p, m)), for n > m >= 0, with input T(m,m) = 4^m. See a comment and references in A282629. An example is given below. - Wolfdieter Lang, Aug 11 2017

A141413 Inverse binomial transform of A140962.

Original entry on oeis.org

0, -1, 1, -3, 9, -27, 81, -243, 729, -2187, 6561, -19683, 59049, -177147, 531441, -1594323, 4782969, -14348907, 43046721, -129140163, 387420489, -1162261467, 3486784401, -10460353203, 31381059609, -94143178827, 282429536481, -847288609443

Offset: 0

Paul Curtz, Aug 04 2008

Vincenzo Librandi, Table of n, a(n) for n = 0..1000
Index entries for linear recurrences with constant coefficients, signature (-3).

Cf. A285066 (alternating row sums, starting here with n >= 2). - Wolfdieter Lang, Apr 23 2017

Cf. A000244, A133494.

[0,-1] cat [(-3)^(n-2): n in [2..30]]; // G. C. Greubel, Mar 30 2021

A141413:= n-> `if`(n<2, -n, (-3)^(n-2)); seq(A141413(n), n=0..30); # G. C. Greubel, Mar 30 2021

CoefficientList[Series[-x(1+2x)/(1+3x), {x, 0, 30}], x] (* Vincenzo Librandi, Feb 24 2013 *)
Join[{0, -1},LinearRecurrence[{-3},{1},26]] (* Ray Chandler, Aug 12 2015 *)

[-n if n<2 else (-3)^(n-2) for n in (0..30)] # G. C. Greubel, Mar 30 2021

a(n) = (-1)^n*A133494(n-1), n>0.
G.f.: (-x)*(1 + 2*x)/(1+3*x). - R. J. Mathar, Nov 11 2008
G.f.: x^2/( Q(0)+ 2*x)- x where Q(k) =  1 - x/(x*(k+1) - 1 )/Q(k+1) ; (recursively defined continued fraction). - Sergei N. Gladkovskii, Feb 24 2013
E.g.f.: (exp(-3*x) - (1 + 6*x))/9. - Wolfdieter Lang, Apr 19 2017
a(n) = (-3)^(n-2) for n >= 2, with a(0) = 0 and a(1) = -1. - G. C. Greubel, Mar 30 2021

Edited and extended by R. J. Mathar, Nov 11 2008

cp's OEIS Frontend

A285067 Row sums of triangle A285066.

Views

Author

Keywords

Comments

Crossrefs

Formula

A285061 Sheffer triangle S2[4,1] = (exp(x), exp(4*x) - 1).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A141413 Inverse binomial transform of A140962.

Views

Author

Keywords

Links

Crossrefs

Programs

Magma

Maple

Mathematica

Sage

Formula

Extensions