A142968 -id:A142968 - cp's OEIS frontend

A156919 Table of coefficients of polynomials related to the Dirichlet eta function.

1, 2, 1, 4, 10, 1, 8, 60, 36, 1, 16, 296, 516, 116, 1, 32, 1328, 5168, 3508, 358, 1, 64, 5664, 42960, 64240, 21120, 1086, 1, 128, 23488, 320064, 900560, 660880, 118632, 3272, 1, 256, 95872, 2225728

Offset: 0

Johannes W. Meijer, Feb 20 2009, Jun 24 2009

Essentially the same as A185411. Row reverse of A185410. - Peter Bala, Jul 24 2012
The SF(z; n) formulas, see below, were discovered while studying certain properties of the Dirichlet eta function.
From Peter Bala, Apr 03 2011: (Start)
Let D be the differential operator 2*x*d/dx. The row polynomials of this table come from repeated application of the operator D to the function g(x) = 1/sqrt(1 - x). For example,
  D(g) = x*g^3
  D^2(g) = x*(2 + x)*g^5
  D^3(g) = x*(4 + 10*x + x^2)*g^7
  D^4(g) = x*(8 + 60*x + 36*x^2 + x^3)*g^9.
Thus this triangle is analogous to the triangle of Eulerian numbers A008292, whose row polynomials come from the  repeated application of the operator x*d/dx to the function 1/(1 - x). (End)

			The first few rows of the triangle are:
  [1]
  [2, 1]
  [4, 10, 1]
  [8, 60, 36, 1]
  [16, 296, 516, 116, 1]
The first few P(z;n) are:
  P(z; n=0) = 1
  P(z; n=1) = 2 + z
  P(z; n=2) = 4 + 10*z + z^2
  P(z; n=3) = 8 + 60*z + 36*z^2 + z^3
The first few SF(z;n) are:
  SF(z; n=0) = (1/2)*(1)/(1-z)^(3/2);
  SF(z; n=1) = (1/4)*(2+z)/(1-z)^(5/2);
  SF(z; n=2) = (1/8)*(4+10*z+z^2)/(1-z)^(7/2);
  SF(z; n=3) = (1/16)*(8+60*z+36*z^2+z^3)/(1-z)^(9/2);
In the Savage-Viswanathan paper, the coefficients appear as
  1;
  1,    2;
  1,   10,     4;
  1,   36,    60,     8;
  1,  116,   516,   296,    16;
  1,  358,  3508,  5168,  1328,   32;
  1, 1086, 21120, 64240, 42960, 5664, 64;
  ...

D. H. Lehmer, Interesting Series Involving the Central Binomial Coefficient, Am. Math. Monthly 92 (1985) 449-457, Polynomial V in eq (17). [R. J. Mathar, Feb 24 2009]
Shi-Mei Ma, A family of two-variable derivative polynomials for tangent and secant, arXiv: 1204.4963v3 [math.CO], 2012.
Shi-Mei Ma, A family of two-variable derivative polynomials for tangent and secant, El J. Combinat. 20 (1) (2013) P11
S.-M. Ma and T. Mansour, The 1/k-Eulerian polynomials and k-Stirling permutations, arXiv preprint arXiv:1409.6525 [math.CO], 2014.
S.-M. Ma and Y.-N. Yeh, Stirling permutations, cycle structures of permutations and perfect matchings, arXiv preprint arXiv:1503.06601 [math.CO], 2015.
Carla D. Savage and Gopal Viswanathan, The 1/k-Eulerian polynomials, Elec. J. of Comb., Vol. 19, Issue 1, #P9 (2012).
Eric Weisstein's World of Mathematics, Dirichlet Eta Function

A142963 and this sequence can be mapped onto the A156920 triangle.
FP1 sequences A000340, A156922, A156923, A156924.
FP2 sequences A050488, A142965, A142966, A142968.
Appears in A162005, A000182, A162006 and A162007.
Cf. A186695, A187075.
Cf. A185410 (row reverse), A185411.

A156919 := proc(n,m) if n=m then 1; elif m=0 then 2^n ; elif m<0 or m>n then 0; else 2*(m+1)*procname(n-1,m)+(2*n-2*m+1)*procname(n-1,m-1) ; end if; end proc: seq(seq(A156919(n,m), m=0..n), n=0..7); # R. J. Mathar, Feb 03 2011

g[0] = 1/Sqrt[1-x]; g[n_] := g[n] = 2x*D[g[n-1], x]; p[n_] := g[n] / g[0]^(2n+1) // Cancel; row[n_] := CoefficientList[p[n], x] // Rest; Table[row[n], {n, 0, 9}] // Flatten (* Jean-François Alcover, Aug 09 2012, after Peter Bala *)
Flatten[Table[Rest[CoefficientList[Nest[2 x D[#, x] &, (1 - x)^(-1/2), k] (1 - x)^(k + 1/2), x]], {k, 10}]] (* Jan Mangaldan, Mar 15 2013 *)

SF(z; n) = Sum_{m >= 1} m^(n-1)*4^(-m)*z^(m-1)*Gamma(2*m+1)/(Gamma(m)^2) = P(z;n) / (2^(n+1)*(1-z)^((2*n+3)/2)) for n >= 0. The polynomials P(z;n) = Sum_{k = 0..n} a(k)*z^k generate the a(n) sequence.
If we write the sequence as a triangle the following relation holds: T(n,m) = (2*m+2)*T(n-1,m) + (2*n-2*m+1)*T(n-1,m-1) with T(n,m=0) = 2^n and T(n,n) = 1, n >= 0 and 0 <= m <= n.
G.f.: 1/(1-xy-2x/(1-3xy/(1-4x/(1-5xy/(1-6x/(1-7xy/(1-8x/(1-... (continued fraction). - Paul Barry, Jan 26 2011
From Peter Bala, Apr 03 2011 (Start)
E.g.f.: exp(z*(x + 2)) * (1 - x)/(exp(2*x*z) - x*exp(2*z))^(3/2) = Sum_{n >= 0} P(x,n)*z^n/n! = 1 + (2 + x)*z + (4 + 10*x + x^2)*z^2/2! + (8 + 60*x + 36*x^2 + x^3)*z^3/3! + ... .
Explicit formula for the row polynomials:
P(x,n-1) = Sum_{k = 1..n} 2^(n-2*k)*binomial(2k,k)*k!*Stirling2(n,k)*x^(k-1)*(1 - x)^(n-k).
The polynomials x*(1 + x)^n * P(x/(x + 1),n) are the row polynomials of A187075.
The polynomials x^(n+1) * P((x + 1)/x,n) are the row polynomials of A186695.
Row sums are A001147(n+1). (End)
Sum_{k = 0..n} (-1)^k*T(n,k) = (-1)^binomial(n,2)*A012259(n+1). - Johannes W. Meijer, Sep 27 2011

Minor edits from Johannes W. Meijer, Sep 27 2011

A156925 FP2 polynomials related to the generating functions of the left hand columns of the A156920 triangle.

Original entry on oeis.org

1, 1, 1, 1, 8, -11, -6, 1, 38, -108, -242, 839, -444, -180, 1, 144, -425, -7382, 48451, -96764, -2559, 257002, -312444, 88344, 30240, 1, 487, 720, -130472, 1277794, -4193514, -6504496

Offset: 0

Johannes W. Meijer, Feb 20 2009

The FP2 polynomials appear in the numerators of the GF2 o.g.f.s. of the left hand columns of A156920. The FP2 can be calculated with the formula of the LHC sequence, see A156920, and the formula for the general structure of the generating function GF2, see below.
An appropriate name for the FP2 polynomials seems to be the flower polynomials of the second kind because the zero patterns of these polynomials look like flowers. The zero patterns of the FP2 and the FP1, see A156921, resemble each other closely.
A Maple program that generates for a left hand column with a certain LHCnr its GF2 and FP2 can be found below. LHCnr stands for left hand column number and starts from 1.

			The first few rows of the "triangle" of the coefficients of the FP2 polynomials.
In the columns the coefficients of the powers of z^m, m=0,1,2,..., appear.
  [1]
  [1, 1]
  [1, 8, -11, -6]
  [1, 38, -108, -242, 839, -444, -180]
  [1, 144, -425, -7382, 48451, -96764, -2559, 257002, -312444, 88344, 30240]
Matrix of the coefficients of the FP2 polynomials. The coefficients in the columns of this matrix are the powers of z^m, m=0,1,2,...
  [1, 0, 0, 0, 0, 0, 0]
  [1, 1, 0, 0, 0, 0, 0]
  [1, 8 , -11, -6, 0, 0, 0]
  [1, 38, -108, -242, 839, -444, -180]
The first few FP2 polynomials are:
  FP2(z; LHCnr = 1) = 1
  FP2(z; LHCnr = 2) = (1+z)
  FP2(z; LHCnr = 3) = 1+8*z-11*z^2-6*z^3
Some GF2(z;LHCnr) are:
  GF2(z; LHCnr = 3) = (1+8*z-11*z^2-6*z^3)/((1-z)^3*(1-2*z)^2*(1-3*z))
  GF2(z; LHCnr = 4) = (1+38*z-108*z^2-242*z^3+839*z^4-444*z^5-180*z^6)/((1-z)^4*(1-2*z)^3*(1-3*z)^2*(1-4*z))

Cf. A156920, A156921, A156927, A156933.
For the first few GF2's see A050488, A142965, A142966 and A142968.
Row sums(n) = A156926(n).
The number of FP2 terms follow the 'Lazy Caterer's sequence' A000124.
For the polynomials in the denominators of the GF2(z;LHCnr) see A157703.

LHCnr:=5; LHCmax:=(LHCnr)*(LHCnr-1)/2: RHCend:=LHCnr+LHCmax: for k from LHCnr to RHCend do for n from 0 to k do S2[k,n]:=sum((-1)^(n+i)*binomial(n,i)*i^k/n!,i=0..n) end do: G(k,x):= sum(S2[k,p]*((2*p)!/p!)*x^p/(1-4*x)^(p+1),p=0..k)/ (((-1)^(k+1)*2*x)/(-1+4*x)^(k+1)): fx:=simplify(G(k,x)): nmax:=degree(fx); for n from 0 to nmax do d[n]:= coeff(fx,x,n)/2^n end do: LHC[n]:=d[LHCnr-1] end do: a:=n-> LHC[n]: seq(a(n), n=LHCnr..RHCend); for nx from 0 to LHCmax do num:=sort(sum(A[t]*z^t,t=0..LHCmax)): nom:=product((1-u*z)^(LHCnr-u+1),u=1..LHCnr); LHCb:=series(num/nom,z,nx+1); y:=coeff(LHCb,z,nx)-A[nx]; x:=LHC[LHCnr+nx]; A[nx]:=x-y; end do: FP2[LHCnr]:=sort(num,z, ascending); GenFun[LHCnr]:= FP2[LHCnr]/ product((1-m*z)^(LHCnr-m+1), m=1..LHCnr);

G.f.: GF2(z; LHCnr) = FP2(z; LHCnr)/Product_{m=1..LHCnr} (1-m*z)^(LHCnr-m+1).

Row sum(n+1) = (-1)^(n)*2*(n+1)!*Row sum(n); Row sum(n=0) = 1.

A156920 Triangle of the normalized A142963 and A156919 sequences.

Original entry on oeis.org

1, 1, 1, 1, 5, 1, 1, 15, 18, 1, 1, 37, 129, 58, 1, 1, 83, 646, 877, 179, 1, 1, 177, 2685, 8030, 5280, 543, 1, 1, 367, 10002, 56285, 82610, 29658, 1636, 1, 1, 749, 34777, 335162, 919615, 756218, 159742, 4916, 1

Offset: 0

Johannes W. Meijer, Feb 20 2009

The originator sequences are A142963 and A156919.
The Flower Triangle seems to be an appropriate name for the triangular array of this sequence. The zero patterns of the Flower Polynomials of the first, see A156921, the second, see A156925, the third, see A156927, and the fourth kind, see A156933, look like flowers.
The first Maple program generates the Flower Triangle sequence.
The second program generates the Right Hand Columns sequences and the third one generates the Left Hand Column sequences. For an explanation of these two algorithms see A142963.

			The first few rows of the triangle are:
  [1]
  [1, 1]
  [1, 5, 1]
  [1, 15, 18, 1]
  [1, 37, 129, 58, 1]
  [1, 83, 646, 877, 179, 1]

Shi-Mei Ma, Some combinatorial sequences associated with context-free grammars, arXiv:1208.3104v2 [math.CO]. - From N. J. A. Sloane, Aug 21 2012

Originator sequences A142963, A156919.
Related sequences A156921, A156925, A156927, A156933.
Left hand column sequences A050488, A142965, A142966, A142968.
Right hand column sequences A000340, A156922, A156923, A156924.
Row sums A014307(n+1).

A156920 := proc(n,m): if n=m then 1; elif m=0 then 1 ; elif m<0 or m>n then 0; else (m+1)*procname(n-1, m)+(2*n-2*m+1)*procname(n-1, m-1) ; end if; end proc: seq(seq(A156920(n, m), m=0..n), n=0..8);
RHCnr:=5; RHCmax:=10; RHCend:=RHCnr+RHCmax: for k from RHCnr to RHCend do for n from 0 to k do S2[k,n]:=sum((-1)^(n+i)*binomial(n,i)*i^k/n!,i=0..n) end do: G(k,x):= sum(S2[k,p]*((2*p)!/p!) *x^p/(1-4*x)^(p+1),p=0..k)/(((-1)^(k+1)*2*x)/(-1+4*x)^(k+1)): fx:=simplify(G(k,x)): nmax:=degree(fx); RHC[k-RHCnr+1]:= coeff(fx,x,k-RHCnr)/2^(k-RHCnr) end do: a:=n-> RHC[n]: seq(a(n), n=1..RHCend-RHCnr);
LHCnr:=5; LHCmax:=10: LHCend:=LHCnr+LHCmax: for k from LHCnr to LHCend do for n from 0 to k do S2[k,n]:=sum((-1)^(n+i)*binomial(n,i)*i^k/n!,i=0..n) end do: G(k,x):= sum(S2[k,p]*((2*p)!/p!)*x^p/(1-4*x)^(p+1),p=0..k)/ (((-1)^(k+1)*2*x)/(-1+4*x)^(k+1)): fx:=simplify(G(k,x)): nmax:=degree(fx); for n from 0 to nmax do d[n]:= coeff(fx,x,n)/2^n end do: LHC[n]:=d[LHCnr-1] end do: a:=n-> LHC[n]: seq(a(n), n=LHCnr..LHCend-1);

T[, 0] = 1; T[n, n_] = 1; T[n_, m_] := T[n, m] = (m + 1)*T[n - 1, m] + (2*n - 2*m + 1)*T[n - 1, m - 1];
Table[T[n, m], {n, 0, 8}, {m, 0, n}] // Flatten (* Jean-François Alcover, Nov 14 2017 *)

T(n,m) = (m+1)*T(n-1,m) + (2*n-2*m+1)*T(n-1,m-1) with T(n,m=0) = 1 and T(n,n) = 1, n>=0 and 0 <= m <= n.
From Peter Bala, Jul 22 2012: (Start)
T(n,k) = 1/(2^(n-k))*A156919(n,k).
E.g.f.: 1 + t*x + (t+t^2)*x^2/2! + (t+5*t^2+t^3)*x^3/3! + ... = sqrt(E(x,2*t)), where E(x,t) = (1-t)*exp(x*t)/(exp(x*t)-t*exp(x)) = 1 + t*x + (t+t^2)*x^2/2! + (t+4*t^2+t^3)*x^3/3! + ... is the e.g.f. for the Eulerian numbers A008292.
The row polynomials R(n,x) satisfy 1/sqrt(1-2*x)*(x*d/dx)^n(1/sqrt(1-2*x)) = R(n,x)/(1-2*x)^(n+1). (End)

Minor edits by Johannes W. Meijer, Sep 28 2011

A142963 Triangle read by rows, coefficients of the polynomials P(k, x) = (1/2) Sum_{p=0..k-1} Stirling2(k, p+1)x^p(1-4x)^(k-1-p)(2*p+2)!/(p+1)!.

Original entry on oeis.org

1, 1, 2, 1, 10, 4, 1, 30, 72, 8, 1, 74, 516, 464, 16, 1, 166, 2584, 7016, 2864, 32, 1, 354, 10740, 64240, 84480, 17376, 64, 1, 734, 40008, 450280, 1321760, 949056, 104704, 128, 1, 1498, 139108, 2681296, 14713840, 24198976, 10223488, 629248, 256, 1, 3030, 462264, 14341992

Offset: 1

Wolfdieter Lang, Sep 15 2008

Previous name: Table of coefficients of row polynomials of certain o.g.f.s.
The o.g.f.s G(k, x) for the k-family of sequences S(k, n):= Sum_{p=0..n} p^k*binomial(2*p, p)*binomial(2*(n-p), n-p), k=0,1,... (convolution of two sequences involving the central binomial coefficients) are 1/(1-4*x) for k=0 and 2*x*P(k, x)/(1-4*x)^(k+1) for k=1,2,..., with the row polynomials P(k, x) = Sum_{m=0..k-1} a(n,m)*x^m).
The author was led to compute the sums S(k, n) by a question asked by M. Greiter, Jun 27 2008.
In order to keep the index k>=1 of Sigma(k, n) also for the polynomials P(k, x), their degree is then k-1.

			Triangle starts:
[1]
[1,   2]
[1,  10,     4]
[1,  30,    72,      8]
[1,  74,   516,    464,      16]
[1, 166,  2584,   7016,    2864,     32]
[1, 354, 10740,  64240,   84480,  17376,     64]
[1, 734, 40008, 450280, 1321760, 949056, 104704, 128]
...
P(3,x) = 1+10*x+4*x^2.
G(3,x) = 2*x*(1+10*x+4*x^2)/(1-4*x)^4.

Vincenzo Librandi, Table of n, a(n) for n = 1..210
Wolfdieter Lang, First 10 rows and more.
L. Liu, Y. Wang, A unified approach to polynomial sequences with only real zeros, arXiv:math/0509207v5 [math.CO], 2005-2006.

Left hand column sequences 2*A142964, 4*A142965, 8*A142966, 16*A142968.
Row sums A142967.
From Johannes W. Meijer, Feb 20 2009: (Start)
A156919 and this sequence can be mapped onto A156920.
Cf. A156921, A156925, A156927, A156933.
Right hand column sequences 2^n*A000340, 2^n*A156922, 2^n*A156923, 2^n*A156924. (End)
Cf. A142961, A142962.

A142963 := proc(n,m): if n=m+1 then 2^(n-1); elif m=0 then 1 ; elif m<0 or m>n-1 then 0; else (m+1)*procname(n-1, m)+(4*n-4*m-2)*procname(n-1, m-1); end if; end proc: seq(seq(A142963(n,m), m=0..n-1), n=1..9); # Johannes W. Meijer, Sep 28 2011
# Alternatively (assumes offset 0):
p := (n,x) -> (1/2)*add(Stirling2(n+1,k+1)*x^k*(1-4*x)^(n-k)*(2*k+2)!/(k+1)!, k=0..n): for n from 0 to 7 do [n], PolynomialTools:-CoefficientList(p(n,x), x) od;
# Peter Luschny, Jun 18 2017

t[, 0] = 1; t[n, m_] /; m == n-1 := 2^m; t[n_, m_] := (m+1)*t[n-1, m] + (4*n-4*m-2)*t[n-1, m-1]; Table[t[n, m], {n, 1, 10}, {m, 0, n-1}] // Flatten (* Jean-François Alcover, Jun 21 2013, after Johannes W. Meijer *)

G(k, x) = Sum_{p=0..k} S2(k, p)*((2*p)!/p!)*x^p/(1-4*x)^(p+1), k >= 0 (here k >= 1), with the Stirling2 triangle S2(k, p):=A048993(k, p). (Proof from the product of the o.g.f.s of the two convoluted sequences and the normal ordering (x^d_x)^k = Sum_{p=0..k} S2(k, p)*x^p*d_x^p, with the derivative operator d_x.)
a(k,m) = [x^m]P(k, x) = [x^m] ((1-4*x)^(k+1))*G(k,x)/(2*x), k>=1, m=0,1,...,k-1.
For the triangle coefficients the following relation holds: T(n,m) = (m+1)*T(n-1,m) + (4*n-4*m-2)*T(n-1,m-1) with T(n,m=0) = 1 and T(n,m=n-1) = 2^(n-1), n >= 1 and 0 <= m <= n-1. - Johannes W. Meijer, Feb 20 2009
From Peter Bala, Jan 18 2018: (Start)
(x*d/dx)^n (1/(sqrt(1 - 4*x)) = 2*x*P(n,x)/sqrt(1 - 4*x)^(n+1/2) for n >= 1.
x*P(n,x)/(1 - 4*x)^(n+1/2) = (1/2)*Sum_{k >= 1} binomial(2*k,k)* k^n*x^k for n >= 1.
P(n+1,x) = ((4*n - 2)*x + 1)*P(n,x) - x*(4*x - 1)*d/dx(P(n,x)).
Hence the polynomial P(n,x) has all real zeros by Liu et al., Theorem 1.1, Corollary 1.2. (End)

Minor edits by Johannes W. Meijer, Sep 28 2011
A more precise name by Peter Luschny, Jun 18 2017
Name reformulated with offset corrected, edited by Wolfdieter Lang, Aug 23 2019

A211399 Triangle T(n,k), 0 <= k <= n, given by (0, 1, 0, 2, 0, 3, 0, 4, 0, 5, 0, ...) DELTA (1, 0, 3, 0, 5, 0, 7, 0, 9, 0, ...) where DELTA is the operator defined in A084938.

Original entry on oeis.org

1, 0, 1, 0, 1, 1, 0, 1, 5, 1, 0, 1, 15, 18, 1, 0, 1, 37, 129, 58, 1, 0, 1, 83, 646, 877, 179, 1, 0, 1, 177, 2685, 8030, 5280, 543, 1, 0, 1, 367, 10002, 56285, 82610, 29658, 1636, 1, 0, 1, 749, 34777, 335162

Offset: 0

Philippe Deléham, Feb 08 2013

Contains A156920 as submatrix.

Row-reversal of A102365. - Philippe Deléham, Feb 12 2013

			Triangle begins :
1
0, 1
0, 1, 1
0, 1, 5, 1
0, 1, 15, 18, 1
0, 1, 37, 129, 58, 1
0, 1, 83, 646, 877, 179, 1

Left hand column sequences: A000007, A000012, A050488, A142965, A142966, A142968.
Right hand column sequences: A000340, A156922, A156923, A156924.
Row sums A014307(n).

T(n,k) = k*T(n-1,k) + (2n-2k+1)*T(n-1,k-1) , T(n,n) = 1, T(n,k) = 0 if k<0 or if k>n.
T(n,k) = A185411(n,k)/(2^(n-k)).
Sum_{k, 0<=k<=n} T(n,k)*x^(n-k) = A000012(n), A014307(n), A001147(n) for x = 0, 1, 2 respectively .
G.f.: 1/(1-xy/(1-x/(1-3xy/(1-2x/(1-5xy/(1-3x/(1-7xy/(1- ...(continued fraction).

cp's OEIS Frontend

A156919 Table of coefficients of polynomials related to the Dirichlet eta function.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

Extensions

A156925 FP2 polynomials related to the generating functions of the left hand columns of the A156920 triangle.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Maple

Formula

A156920 Triangle of the normalized A142963 and A156919 sequences.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

Extensions

A142963 Triangle read by rows, coefficients of the polynomials P(k, x) = (1/2) Sum_{p=0..k-1} Stirling2(k, p+1)*x^p*(1-4*x)^(k-1-p)*(2*p+2)!/(p+1)!.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

Extensions

A211399 Triangle T(n,k), 0 <= k <= n, given by (0, 1, 0, 2, 0, 3, 0, 4, 0, 5, 0, ...) DELTA (1, 0, 3, 0, 5, 0, 7, 0, 9, 0, ...) where DELTA is the operator defined in A084938.

Views

Author

Keywords

Comments

Examples

Crossrefs

Formula

A142963 Triangle read by rows, coefficients of the polynomials P(k, x) = (1/2) Sum_{p=0..k-1} Stirling2(k, p+1)x^p(1-4x)^(k-1-p)(2*p+2)!/(p+1)!.