A021874 -id:A021874 - cp's OEIS frontend

A282629 Sheffer triangle (exp(x), exp(3*x) - 1). Named S2[3,1].

1, 1, 3, 1, 15, 9, 1, 63, 108, 27, 1, 255, 945, 594, 81, 1, 1023, 7380, 8775, 2835, 243, 1, 4095, 54729, 109890, 63180, 12393, 729, 1, 16383, 395388, 1263087, 1151010, 387828, 51030, 2187, 1, 65535, 2816865, 13817034, 18752391, 9658278, 2133054, 201204, 6561, 1, 262143, 19914660, 146620935, 285232185, 210789621, 69502860, 10825650, 767637, 19683

Offset: 0

Wolfdieter Lang, Apr 03 2017

For Sheffer triangles (infinite lower triangular exponential convolution matrices) see the W. Lang link under A006232, with references).
The e.g.f. for the sequence of column m is (Sheffer property) exp(x)*(exp(3*x) - 1)^m/m!.
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].
The a-sequence for this Sheffer triangle has e.g.f. 3*x/log(1+x) and is 3*A006232(n)/ A006233(n) (Cauchy numbers of the first kind).
The z-sequence has e.g.f. (3/(log(1+x)))*(1 - 1/(1+x)^(1/3)) and is A284857(n) / A284858(n).
The main diagonal gives A000244.
The row sums give A284859. The alternating row sums give A284860.
The triangle appears in the o.g.f. G(n, x) of the sequence {(1 + 3*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 + 3*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 A225117 with row reversion.
The first column k sequences divided by 3^k are A000012, A002450 (with a leading 0), A016223, A021874. For the e.g.f.s and o.g.f.s see below. - Wolfdieter Lang, Apr 09 2017
From Wolfdieter Lang, Aug 09 2017: (Start)
The general row polynomials R(d,a;n,x) = Sum_{k=0..n} T(d,a;n,m)*x^m of the Sheffer triangle S2[d,a] satisfy, as special polynomials of the Boas-Buck class, the identity (see the reference, and we use the notation of Rainville, Theorem 50, p. 141, adapted to an exponential generating function)
  (E_x - n*1)*R(d,a;n,x) = - n*a*R(d,a;n-1,x) - Sum_{k=0..n-1} binomial(n, k+1)*(-d)^(k+1)*Bernoulli(k+1)*E_x*R(d,a;n-1-k,x), with E_x = x*d/dx (Euler operator).
This entails a recurrence for the sequence of column m, for n > m:
  T(d,a;n,m) = (1/(n - m))*[(n/2)*(2*a + d*m)*T(d,a;n-1,m) + m*Sum_{p=m..n-2} binomial(n,p)(-d)^(n-p)*Bernoulli(n-p)*T(d,a;p,m)], with input T(d,a;n,n) = d^n. For the present [d,a] = [3,1] case see the formula and example sections below. - Wolfdieter Lang, Aug 09 2017 (End)
The inverse of this triangular Sheffer matrix S2[3,1] is S1[3,1] with rational elements S1[3,1](n, k) = (-1)^(n-k)*A286718(n, k)/3^k. - Wolfdieter Lang, Nov 15 2018
Named after the American mathematician Isador Mitchell Sheffer (1901-1992). - Amiram Eldar, Jun 19 2021

			The triangle T(n, m) begins:
  n\m 0      1        2         3         4         5        6        7      8     9
  0:  1
  1:  1      3
  2:  1     15        9
  3:  1     63      108        27
  4:  1    255      945       594        81
  5:  1   1023     7380      8775      2835       243
  6:  1   4095    54729    109890     63180     12393      729
  7:  1  16383   395388   1263087   1151010    387828    51030     2187
  8:  1  65535  2816865  13817034  18752391   9658278  2133054   201204   6561
  9:  1 262143 19914660 146620935 285232185 210789621 69502860 10825650 767637 19683
  ...
------------------------------------------------------------------------------------
Nontrivial recurrence for m=0 column from z-sequence: T(4,0) = 4*(1*1 + 63*(-1/6) + 108*(11/54) + 27*(-49/108)) = 1.
Recurrence for m=2 column from a-sequence: T(4, 2) = (4/2)*(1*63*3 + 2*108*(3/2) + 3*27*(-3/6)) = 945.
Recurrence for row polynomial R(3, x) (Meixner type): ((3*x + 1) + 3*x*d_x)*(1 + 15*x + 9*x^2) = 1 + 63*x + 108*x^2 + 27*x^3.
E.g.f. and o.g.f. of n = 1 powers {(1 + 3*m)^1}_{m>=0} A016777: E(1, x) = exp(x) * (T(1, 0) + T(1, 1)*x) = exp(x)*(1+3*x). O.g.f.: G(1, x) = T(1, 0)*0!/(1-x) + T(1, 1)*1!*x/(1-x)^2 = (1+2*x)/(1-x)^2.
Boas-Buck recurrence for column m = 2, and n = 4: T(4, 2) = (1/2)*(2*(2 + 3*2)*T(3, 2) + 2*6*(-3)^2*bernoulli(2)*T(2, 2)) = (1/2)*(16*108 + 12*9*(1/6)*9) = 945. - _Wolfdieter Lang_, Aug 09 2017

Ralph P. Boas, Jr. and R. Creighton Buck, Polynomial Expansions of analytic functions, Springer, 1958, pp. 17 - 21, (last sign in eq. (6.11) should be -).
Earl D. Rainville, Special Functions, The Macmillan Company, New York, 1960, ch. 8, sect. 76, 140 - 146.

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:math/1707.04451 [math.NT], July 2017.
Wolfdieter Lang, On Generating functions of Diagonals Sequences of Sheffer and Riordan Number Triangles, arXiv:1708.01421 [math.NT], August 2017.

Cf. A000012, A000244, A006232/A006233, A016777, A024036, A111577, A225117, A225466, A284857, A284858, A284859, A284860, A284861, A286718.

Table[Sum[Binomial[m, k] (-1)^(k - m) (1 + 3 k)^n/m!, {k, 0, m}], {n, 0, 9}, {m, 0, n}] // Flatten (* Michael De Vlieger, Apr 08 2017 *)

T(n, m) = sum(k=0, m, binomial(m, k) * (-1)^(k - m) * (1 + 3*k)^n/m!);
for(n=0, 9, for(m=0, n, print1(T(n, m),", ");); print();) \\ Indranil Ghosh, Apr 08 2017

A nontrivial recurrence for the column m=0 entries T(n, 0) = 1 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) = ((3*x+1) + 3*x*d_x)*R(n-1, x), with differentiation d_x, for n >= 1, with input R(0, x) = 1.
T(n, m) = Sum_{k=0..m} binomial(m,k)*(-1)^(k-m)*(1 + 3*k)^n/m!, 0 <= m <= n.
E.g.f. of triangle: exp(z)*exp(x*(exp(3*z)-1)) (Sheffer type).
E.g.f. for sequence of column m is exp(x)*((exp(3*x) - 1)^m)/m! (Sheffer property).
From Wolfdieter Lang, Apr 09 2017: (Start)
Standard three-term recurrence: T(n, m) = 0 if n < m, T(n,-1) = 0, T(0, 0) = 1, T(n, m) = 3*T(n-1, m-1) + (1+3*m)*T(n-1, m) for n >= 1. From the T(n, m) formula. Compare with the recurrence of S2[3,2] given in A225466.
The o.g.f. for sequence of column m is 3^m*x^m/Product_{j=0..m} (1 - (1+3*j)*x). (End)
In terms of Stirling2 = A048993: T(n, m) = Sum_{k=0..n} binomial(n, k)* 3^k*Stirling2(k, m), 0 <= m <= n. - Wolfdieter Lang, Apr 13 2017
Boas-Buck recurrence for column sequence m: T(n, m) = (1/(n - m))*((n/2)*(2 + 3*m)*T(n-1, m) + m*Sum_{p=m..n-2} binomial(n,p)*(-3)^(n-p)*Bernoulli(n-p)*T(p, m)), for n > m >= 0, with input T(m, m) = 3^m. See a comment above. - Wolfdieter Lang, Aug 09 2017

A016223 Expansion of 1/((1-x) * (1-4x) (1-7*x)).

Original entry on oeis.org

1, 12, 105, 820, 6081, 43932, 312985, 2212740, 15576561, 109385452, 767096265, 5375266260, 37649233441, 263634112572, 1845796701945, 12922008569380, 90459786608721, 633241412753292, 4432781515242025, 31029837110570100, 217210325789494401, 1520478144588475612

Offset: 0

N. J. A. Sloane

Index entries for linear recurrences with constant coefficients, signature (12, -39, 28).

Cf. A002450, A021874, A111577.

Cf. A282629, A383627.

a:=n->sum((7^(n+1-j)-4^(n+1-j))/3, j=0..n+1): seq(a(n), n=0..20); # Zerinvary Lajos, Jan 15 2007

a(n) = (1-2*4^(n+2)+7^(n+2))/18; \\ Seiichi Manyama, May 03 2025

a(n) = (1/18) - (16/9)*4^n + (49/18)*7^n. - Antonio Alberto Olivares, Feb 07 2010 [corrected by Seiichi Manyama, May 03 2025]
a(0)=1, a(1)=12, a(n) = 11*a(n-1) - 28*a(n-2) + 1. - Vincenzo Librandi, Feb 10 2011
E.g.f.: exp(x)*(1 - 32*exp(3*x) + 49*exp(6*x))/(2!*3^2). - This is (d^2/dx^2) (exp(x)*(exp(x) - 1)^2 / (2*3^2)). See also the second column of the Sheffer triangle A282629 divided by 3^2. - Wolfdieter Lang, Apr 08 2017
From Seiichi Manyama, May 03 2025: (Start)
a(n) = Sum_{k=0..n} 3^k * binomial(n+2,k+2) * Stirling2(k+2,2).
G.f.: B(x)^3, where B(x) is the g.f. of A383627. (End)

A284861 Triangle read by rows: T(n, k) = S2[3,1](n, k)k! with the Sheffer triangle S2[3,1] = (exp(x), exp(3x) -1) given in A282629.

Original entry on oeis.org

1, 1, 3, 1, 15, 18, 1, 63, 216, 162, 1, 255, 1890, 3564, 1944, 1, 1023, 14760, 52650, 68040, 29160, 1, 4095, 109458, 659340, 1516320, 1487160, 524880, 1, 16383, 790776, 7578522, 27624240, 46539360, 36741600, 11022480, 1, 65535, 5633730, 82902204, 450057384, 1158993360, 1535798880, 1014068160, 264539520

Offset: 0

Wolfdieter Lang, Apr 09 2017

This is a generalization of triangle A131689(n, k) = Stirling2(n, k)*k!, because S2[3,1] is a generalization of the Stirling2 triangle written as S2[1,0].
This triangle appears in the o.g.f. G(3,1;n,x) of the powers {(1+3*m)^n}{m>=0} as  G(3,1;n,x) = Sum{k>=0..n} T(n, k)*x^k / (1-x)^k.
This triangle is also related to the generalized row reversed Euler triangle rEu[3,1] with row polynomial rEu(3,1;n,x) = Sum_{m=0..n} rEu(3,1;n,m)*x^m with rEu(3,1;n,m) = Sum_{j=0..m} (-1)^(m-j)*binomial(n-j, m-j)*T(n, m). This follows from the above given o.g.f. of powers G(3,1;n,x) = rEu(3,1;n,x)/(1-x)^(n+1). The Euler triangle E[3,1] (row reversed rEu[3,1] is given in A225117. See a formula below.
The e.g.f. of the row polynomials R(3,1;n,x) = Sum_{m=0..n} T(n, m)*x^m follows from the e.g.f. of the row polynomials of the Sheffer triangle A282629. See the formula section.
The diagonal sequence is A032031(k) = k!*3^k.
The row sums give unsigned A151919, and the alternating row sums give A122803.
The first column k sequences divided by A032031(k) are A000012, A002450 (with a leading 0), A016223, A021874. For the e.g.f.s and o.g.f.s see below.

			The triangle T(n, k) begins
n\k 0     1      2       3        4        5        6        7 ...
0:  1
1:  1     3
2:  1    15     18
3:  1    63    216     162
4:  1   255   1890    3564     1944
5:  1  1023  14760   52650    68040    29160
6:  1  4095 109458  659340  1516320  1487160   524880
7:  1 16383 790776 7578522 27624240 46539360 36741600 11022480
...
row n=8: 1 65535 5633730 82902204 450057384 1158993360 1535798880 1014068160 264539520,
row n=9: 1 262143 39829320 879725610 6845572440 25294754520 50042059200 54561276000 30951123840 7142567040,
row n=10: 1 1048575 280378098 9155719980 99549149040 507399658920 1406104706160 2251231315200 2083248720000 1035672220800 214277011200.
------------------------------------------------------------------
T(2, 1) =  -1 + 4^2 = 15 = 2*A225117(2,2) + 1*A225117(2,1) = 2*1 + 1*13.

P. Bala, Deformations of the Hadamard product of power series
M. Dukes, C. D. White, Web Matrices: Structural Properties and Generating Combinatorial Identities, arXiv:1603.01589 [math.CO], 2016.
Wolfdieter Lang, On Sums of Powers of Arithmetic Progressions, and Generalized Stirling, Eulerian and Bernoulli numbers, arXiv:1707.04451 [math.NT], 2017.

Cf. A000012, A002450, A016223, A021874, A032031, A111577, A122803, A131689, |A151919|, A225117, A282629, A019538, A145901.

Table[Sum[Binomial[k, m] (-1)^(k - m) (1 + 3m)^n, {m, 0, k}], {n, 0, 10}, {k, 0, n}]// Flatten (* Indranil Ghosh, Apr 09 2017 *)

for(n=0, 10, for(k=0, n, print1(sum(m=0, k, binomial(k, m) * (-1)^(k - m)*(1 + 3*m)^n),", "); ); print();) \\ Indranil Ghosh, Apr 09 2017

# Indranil Ghosh, Apr 09 2017
from sympy import binomial
for n in range(11):
    print([sum([binomial(k, m)*(-1)**(k - m)*(1 + 3*m)**n for m in range(k + 1)]) for k in range(n + 1)])

E.g.f. of the row polynomials R(n, x) (see a comment above) is exp(z)/(1 - x*(exp(3*z) - 1)). This is the e.g.f. for the triangle.
T(n, k) = Sum_{m=0..k} binomial(k, m)*(-1)^(k-m)*(1+ 3*m)^n, 0 <= k <= n.
T(n, k) = Sum_{m=0..k} binomial(n-m, k-m)*A225117(n,n-m), 0 <= k <= n.
Three term recurrence: T(n, k) = 0 if n < k, T(n,-1) = 0, T(0, 0) = 1, T(n, k) = 3*k*T(n-1, k-1) + (1+3*k)*T(n-1, k) for  n >= 1. See A282629.
The column k sequence has e.g.f. exp(x)*(exp(3*x) - 1)^k (from the Sheffer property of A282629).
  The o.g.f. is A032031(k)*x^k/Product_{j=0..k} (1 - (1+3*j)*x).
From Peter Bala, Jan 12 2018: (Start)
n-th row polynomial R(n,x) = (1 + 3*x) o (1 + 3*x) o ... o (1 + 3*x) (n factors), where o denotes the black diamond multiplication operator of Dukes and White. See example E14 in the Bala link. Cf. A145901.
R(n,x) = Sum_{k = 0..n} binomial(n,k)*3^k*F(k,x) where F(k,x) is the Fubini polynomial of order k, the k-th row polynomial of A019538. (End)

cp's OEIS Frontend

A282629 Sheffer triangle (exp(x), exp(3*x) - 1). Named S2[3,1].

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A016223 Expansion of 1/((1-x) * (1-4x) (1-7*x)).

Views

Author

Keywords

Links

Crossrefs

Programs

Maple

PARI

Formula

A284861 Triangle read by rows: T(n, k) = S2[3,1](n, k)k! with the Sheffer triangle S2[3,1] = (exp(x), exp(3x) -1) given in A282629.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Python

Formula

cp's OEIS Frontend

A282629 Sheffer triangle (exp(x), exp(3*x) - 1). Named S2[3,1].

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A016223 Expansion of 1/((1-x) * (1-4*x) * (1-7*x)).

Views

Author

Keywords

Links

Crossrefs

Programs

Maple

PARI

Formula

A284861 Triangle read by rows: T(n, k) = S2[3,1](n, k)*k! with the Sheffer triangle S2[3,1] = (exp(x), exp(3*x) -1) given in A282629.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Python

Formula

A016223 Expansion of 1/((1-x) * (1-4x) (1-7*x)).

A284861 Triangle read by rows: T(n, k) = S2[3,1](n, k)k! with the Sheffer triangle S2[3,1] = (exp(x), exp(3x) -1) given in A282629.