A075515 -id:A075515 - cp's OEIS frontend

A075513 Triangle read by rows. T(n, m) are the coefficients of Sidi polynomials.

1, -1, 2, 1, -8, 9, -1, 24, -81, 64, 1, -64, 486, -1024, 625, -1, 160, -2430, 10240, -15625, 7776, 1, -384, 10935, -81920, 234375, -279936, 117649, -1, 896, -45927, 573440, -2734375, 5878656, -5764801, 2097152, 1, -2048, 183708, -3670016, 27343750, -94058496, 161414428, -134217728, 43046721

Offset: 1

Wolfdieter Lang, Oct 02 2002

Coefficients of the Sidi polynomials (-1)^(n-1)*D_{n-1,1,n-1}(x), for n >=1, where D_{k,n,m}(z) is given in Theorem 4.2., p. 862, of Sidi [1980].
The row polynomials p(n, x) := Sum_{m=0..n-1} a(n, m)x^m, n >= 1, are obtained from ((Eu(x)^n)*(x-1)^n)/(n*x), where Eu(x) := xd/dx is the Euler-derivative with respect to x.
The row polynomials p(n, y) := Sum_{m=0..n-1} a(n, m)*y^m, n >= 1, are also obtained from ((d^m/dx^m)((exp(x)-1)^m)/m)/exp(x) after replacement of exp(x) by y. Here (d^m/dx^m)f(x), m >= 1, denotes m-fold differentiation of f(x) with respect to x.
b(k,m,n) := (Sum_{p=0..m-1} (a(m, p)*((p+1)*k)^n))/(m-1)!, n >= 0, has g.f. 1/Product_{p=1..m} (1 - k*p*x) for k = 1, 2,... and m = 1, 2,...
The (signed) row sums give A000142(n-1), n >= 1, (factorials) and (unsigned) A074932(n).
The (unsigned) columns give A000012 (powers of 1), 2*A001787(n+1), (3^2)*A027472(n), (4^3)*A038846(n-1), (5^4)*A036071(n-5), (6^5)*A036084(n-6), (7^6)*A036226(n-7), (8^7)*A053107(n-8) for m=0..7.
Right edge of triangle is A000169. - Michel Marcus, May 17 2013

			The triangle T(n, m)  begins:
  n\m 0     1      2        3        4         5         6          7       8
  1:  1
  2: -1     2
  3:  1    -8      9
  4: -1    24    -81       64
  5:  1   -64    486    -1024      625
  6: -1   160  -2430    10240   -15625      7776
  7:  1  -384  10935   -81920   234375   -279936    117649
  8: -1   896 -45927   573440 -2734375   5878656  -5764801    2097152
  9:  1 -2048 183708 -3670016 27343750 -94058496 161414428 -134217728 4304672
  [Reformatted by _Wolfdieter Lang_, Oct 12 2022]
-----------------------------------------------------------------------------
p(2,x) = -1+2*x = (1/(2*x))*x*(d/dx)*x*(d/dx)*(x-1)^2.

A. Sidi, Practical Extrapolation Methods: Theory and Applications, Cambridge University Press, Cambridge, 2003.

Wolfdieter Lang, On a Certain Family of Sidi Polynomials, May 2023.
Harlan J. Brothers, Pascal's triangle, Sidi polynomials, and powers of e, Missouri J. Math. Sci. (2025) Vol. 37, No. 1, 67-78.
Doron S. Lubinsky and Herbert Stahl, Some Explicit Biorthogonal Polynomials, (IN) Approximation Theory XI, (C.K. Chui, M. Neamtu, L. Schumaker, eds.), Nashboro Press, Nashville, 2005, pp. 279-285.
Avram Sidi, Numerical Quadrature and Nonlinear Sequence Transformations; Unified Rules for Efficient Computation of Integrals with Algebraic and Logarithmic Endpoint Singularities, Math. Comp., 35 (1980), 851-874.

Cf. A075510, A075511, A075512, A074932, A075515, A075516, A075906..A075925, A076002..A076013, A258773.

# Assuming offset 0.
seq(seq((-1)^(n-k)*binomial(n, k)*(k+1)^n, k=0..n), n=0..8);
# Alternative:
egf := x -> 1/(exp(LambertW(-exp(-x)*x*y) + x) - x*y):
ser := x -> series(egf(x), x, 12):
row := n -> seq(coeff(n!*coeff(ser(x), x, n), y, k), k=0..n):
seq(print(row(n)), n = 0..8); # Peter Luschny, Oct 21 2022

p[n_, x_] := p[n, x] = Nest[ x*D[#, x]& , (x-1)^n, n]/(n*x); a[n_, m_] := Coefficient[ p[n, x], x, m]; Table[a[n, m], {n, 1, 9}, {m, 0, n-1}] // Flatten (* Jean-François Alcover, Jul 03 2013 *)

tabl(nn) = {for (n=1, nn, for (m=0, n-1, print1((-1)^(n-m-1)*binomial(n-1, m)*(m+1)^(n-1), ", ");); print(););} \\ Michel Marcus, May 17 2013

T(n, m) = ((-1)^(n-m-1)) binomial(n-1, m)*(m+1)^(n-1), n >= m+1 >= 1, else 0.
G.f. for m-th column: ((m+1)^m)(x/(1+(m+1)*x))^(m+1), m >= 0.
E.g.f.: -LambertW(-x*y*exp(-x))/((1+LambertW(-x*y*exp(-x)))*x*y). - Vladeta Jovovic, Feb 13 2008 [corrected for offset 0 <= m <= n. For offset n >= 1 take the integral over x. - Wolfdieter Lang, Oct 12 2022]
T(n, k) = S(n, k+1) / n where S(, ) is triangle in A258773. - Michael Somos, May 13 2018
E.g.f. of column k, with offset n >= 0: exp(-(k + 1)*x)*((k + 1)*x)^k/k!. - Wolfdieter Lang, Oct 20 2022
E.g.f: 1/(exp(LambertW(-exp(-x)*x*y) + x) - x*y) assuming offset = 0. - Peter Luschny, Oct 21 2022

A075498 Stirling2 triangle with scaled diagonals (powers of 3).

Original entry on oeis.org

1, 3, 1, 9, 9, 1, 27, 63, 18, 1, 81, 405, 225, 30, 1, 243, 2511, 2430, 585, 45, 1, 729, 15309, 24381, 9450, 1260, 63, 1, 2187, 92583, 234738, 137781, 28350, 2394, 84, 1, 6561, 557685, 2205225, 1888110, 563031, 71442, 4158, 108, 1

Offset: 1

Wolfdieter Lang, Oct 02 2002

This is a lower triangular infinite matrix of the Jabotinsky type. See the D. E. Knuth reference given in A039692 for exponential convolution arrays.
The row polynomials p(n,x) := Sum_{m=1..n} a(n,m)x^m, n >= 1, have e.g.f. J(x; z)= exp((exp(3*z) - 1)*x/3) - 1.
Subtriangle of the triangle given by (0, 3, 0, 6, 0, 9, 0, 12, 0, 15, 0, ...) DELTA (1, 0, 1, 0, 1, 0, 1, 0, 1, 0, ...) where DELTA is the operator defined in A084938, see example. - Philippe Deléham, Feb 13 2013
Also the Bell transform of A000244. For the definition of the Bell transform see A264428. - Peter Luschny, Jan 26 2016

			[1]; [3,1]; [9,9,1]; ...; p(3,x) = x*(9 + 9*x + x^2).
From _Philippe Deléham_, Feb 13 2013: (Start)
Triangle (0, 3, 0, 6, 0, 9, 0, 12, 0, 15, 0, ...) DELTA (1, 0, 1, 0, 1, 0, 1, 0, ...) begins:
  1;
  0,   1;
  0,   3,   1;
  0,   9,   9,   1;
  0,  27,  63,  18,   1;
  0,  81, 405, 225,  30,   1;
(End)

Andrew Howroyd, Table of n, a(n) for n = 1..1275

Columns 1-7 are A000244, A016137, A017933, A028085, A075515, A075516, A075906. Row sums are A004212.

Cf. A075497, A075499.

# The function BellMatrix is defined in A264428.
# Adds (1, 0, 0, 0, ..) as column 0.
BellMatrix(n -> 3^n, 9); # Peter Luschny, Jan 26 2016

Flatten[Table[3^(n - m) StirlingS2[n, m], {n, 11}, {m, n}]] (* Indranil Ghosh, Mar 25 2017 *)
rows = 9;
t = Table[3^n, {n, 0, rows}];
T[n_, k_] := BellY[n, k, t];
Table[T[n, k], {n, 1, rows}, {k, 1, n}] // Flatten (* Jean-François Alcover, Jun 22 2018, after Peter Luschny *)

for(n=1, 11, for(m=1, n, print1(3^(n - m) * stirling(n, m, 2),", ");); print();) \\ Indranil Ghosh, Mar 25 2017

a(n, m) = (3^(n-m)) * stirling2(n, m).
a(n, m) = (Sum_{p=0..m-1} A075513(m, p)*((p+1)*3)^(n-m))/(m-1)! for n >= m >= 1, else 0.
a(n, m) = 3*m*a(n-1, m) + a(n-1, m-1), n >= m >= 1, else 0, with a(n, 0) := 0 and a(1, 1)=1.
G.f. for m-th column: (x^m)/Product_{k=1..m}(1-3*k*x), m >= 1.
E.g.f. for m-th column: (((exp(3*x)-1)/3)^m)/m!, m >= 1.
From Peter Bala, Jan 13 2018: (Start)
Dobinski-type formulas for row polynomials R(n,x):
R(n,x) = exp(-x/3)*Sum_{i >= 0} (3*i)^n* (x/3)^i/i!;
R(n+1,x) = x*exp(-x/3)*Sum_{i >= 0} (3 + 3*i)^n* (x/3)^i/i!.
R(n+1,x) = x*Sum_{k = 0..n} binomial(n,k)*3^(n-k)*R(k,x).(End)

A099923 Fourth powers of Lucas numbers A000032.

Original entry on oeis.org

16, 1, 81, 256, 2401, 14641, 104976, 707281, 4879681, 33362176, 228886641, 1568239201, 10750371856, 73680216481, 505022001201, 3461445366016, 23725169980801, 162614549665681, 1114577187760656, 7639424429247601

Offset: 0

Ralf Stephan, Nov 01 2004

Arthur T. Benjamin and Jennifer J. Quinn, Proofs that really count: the art of combinatorial proof, M.A.A. 2003, id. 56.

G. C. Greubel, Table of n, a(n) for n = 0..1000
Mohammad K. Azarian, Identities Involving Lucas or Fibonacci and Lucas Numbers as Binomial Sums, International Journal of Contemporary Mathematical Sciences, Vol. 7, No. 45, 2012, pp. 2221-2227.
Pridon Davlianidze, Problem B-1270, Elementary Problems and Solutions, The Fibonacci Quarterly, Vol. 58, No. 2 (2020), p. 179; Four Telescopic Infinite Products, Solution to Problem B-1270 by Jason L. Smith, ibid., Vol. 59, No. 2 (2021), pp. 183-184.
Toufik Mansour, A formula for the generating functions of powers of Horadam's sequence, Australas. J. Combin. 30 (2004) 207-212.
Index entries for linear recurrences with constant coefficients, signature (5,15,-15,-5,1).

Cf. A000032, A000045, A001254, A001622, A075515.

Fourth row of array A103324.

[ Lucas(n)^4 : n in [0..120]]; // Vincenzo Librandi, Apr 14 2011

LucasL[Range[0,20]]^4 (* or *) LinearRecurrence[{5,15,-15,-5,1},{16,1,81,256,2401},21] (* Harvey P. Dale, Jul 04 2014 *)
CoefficientList[Series[(16 - 79 x - 164 x^2 + 76 x^3 + x^4)/((1 - x) (1 + 3*x+x^2)*(1-7*x+x^2)), {x,0,50}], x] (* G. C. Greubel, Dec 21 2017 *)

for(n=0, 30, print1( (fibonacci(n+1) + fibonacci(n-1))^4, ", ")) \\ G. C. Greubel, Dec 21 2017

x='x+O('x^30); Vec((16-79*x-164*x^2+76*x^3+x^4)/((1-x)*(1+3*x+x^2)*(1-7*x+x^2))) \\ G. C. Greubel, Dec 21 2017

a(n) = A000032(n)^4 = A001254(n)^2.
a(n) = L(4*n) + 4*(-1)^n*L(2*n) + 6.
a(n) = L(n-2)*L(n-1)*L(n+1)*L(n+2) + 25, for n >=1.
G.f.: (16-79*x-164*x^2+76*x^3+x^4)/((1-x)*(1+3*x+x^2)*(1-7*x+x^2)). [See Mansour p. 207] - R. J. Mathar, Oct 26 2008
a(0)=16, a(1)=1, a(2)=81, a(3)=256, a(4)=2401, a(n) = 5*a(n-1) + 15*a(n-2) - 15*a(n-3) - 5*a(n-4) + a(n-5). - Harvey P. Dale, Jul 04 2014
Sum_{i=0..n} a(i) = 11 + 6*n + 4*(-1)^n*F(2*n+1) + F(4*n+2), for F = A000045. - Adam Mohamed and Greg Dresden, Jul 02 2024
Product_{n>=2} (1 - 25/a(n)) = phi^5/18, where phi is the golden ratio (A001622) (Davlianidze, 2020). - Amiram Eldar, Dec 04 2024

A075516 Sixth column of triangle A075498.

Original entry on oeis.org

1, 63, 2394, 71442, 1848987, 43615341, 964942308, 20385709344, 416206043253, 8280505692459, 161494678323342, 3101091077181006, 58823743379417199, 1104995938593100617, 20595841868175915096

Offset: 0

Wolfdieter Lang, Oct 02 2002

The e.g.f. given below is Sum_{m=0..5} A075513(6,m)*exp(3*(m+1)*x)/5!.

Cf. A008277, A075498, A075513, A075515, A075906.

a(n) = A075498(n+6, 6) = (3^n)*S2(n+6, 6) with S2(n, m) := A008277(n, m) (Stirling2).
a(n) = Sum_{m=0..5} A075513(6, m)*((m+1)*3)^n/5!.
G.f.: 1/Product_{k=1..6} (1 - 3*k*x).
E.g.f.: (d^6/dx^6)(((exp(3*x)-1)/3)^6)/6! = (-exp(3*x) + 160*exp(6*x) - 2430*exp(9*x) + 10240*exp(12*x) - 15625*exp(15*x) + 7776*exp(18*x))/5!.

A028085 Expansion of 1/((1-3x)(1-6x)(1-9x)(1-12x)).

Original entry on oeis.org

1, 30, 585, 9450, 137781, 1888110, 24862545, 318755250, 4012058061, 49847787990, 613622150505, 7503229474650, 91300979746341, 1106997911204670, 13386607046238465, 161563913916523650

Offset: 0

N. J. A. Sloane

Index entries for linear recurrences with constant coefficients, signature (30,-315,1350,-1944).

Fourth column of triangle A075498.

Cf. A017933, A075515.

CoefficientList[Series[1/((1-3x)(1-6x)(1-9x)(1-12x)),{x,0,30}],x] (* or *) LinearRecurrence[{30,-315,1350,-1944},{1,30,585,9450},30] (* Harvey P. Dale, Feb 06 2015 *)

Vec(1/((1-3*x)*(1-6*x)*(1-9*x)*(1-12*x))+O(x^99)) \\ Charles R Greathouse IV, Sep 26 2012

a(n) = (3^n)*Stirling2(n+4, 4), n >= 0, with Stirling2(n, m) = A008277(n, m).
a(n) = Sum_{m=0..3} (A075513(4, m)*((m+1)*3)^n)/3!.
G.f.: 1/Product_{k=1..4} (1-3*k*x).
E.g.f.: (d^4/dx^4)((((exp(3*x)-1)/3)^4)/4!) = Sum_{m=0..3} (A075513(4, m)*exp(3*(m+1)*x))/3!.
a(n) = (12^(n+3) - 3*9^(n+3) + 3*6^(n+3) - 3^(n+3))/162. - Yahia Kahloune, Jun 10 2013
a(0)=1, a(1)=30, a(2)=585, a(3)=9450, a(n) = 30*a(n-1) - 315*a(n-2) + 1350*a(n-3) - 1944*a(n-4). - Harvey P. Dale, Feb 06 2015

cp's OEIS Frontend

A075513 Triangle read by rows. T(n, m) are the coefficients of Sidi polynomials.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

A075498 Stirling2 triangle with scaled diagonals (powers of 3).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

A099923 Fourth powers of Lucas numbers A000032.

Views

Author

Keywords

References

Links

Crossrefs

Programs

Magma

Mathematica

PARI

PARI

Formula

A075516 Sixth column of triangle A075498.

Views

Author

Keywords

Comments

Crossrefs

Formula

A028085 Expansion of 1/((1-3x)(1-6x)(1-9x)(1-12x)).

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

PARI

Formula