A156933 -id:A156933 - cp's OEIS frontend

A157705 G.f.s of the z^p coefficients of the polynomials in the GF4 denominators of A156933.

1, 1, 3, 2, 18, 128, 171, 42, 1, 22, 1348, 11738, 26293, 17693, 3271, 115, 13, 6122, 228986, 2070813, 6324083, 7397855, 3361536, 544247, 24590, 155, 3, 17248, 2413434, 67035224, 612026240, 2274148882

Offset: 0

Johannes W. Meijer, Mar 07 2009

The formula for the PDGF4(z;n) polynomials in the GF4 denominators of A156933 can be found below.
The general structure of the GFKT4(z;p) that generate the z^p coefficients of the PDGF4(z;n) polynomials can also be found below. The KT4(z;p) polynomials in the numerators of the GFKT4(z;p) have a nice symmetrical structure.
The sequence of the number of terms of the first few KT4(z;p) polynomials is 1, 3, 5, 7, 10, 13, 15, 18, 20, 23, 26, 29, 32, 34, 37, 40, 42. The differences of this sequence and that of the number of terms of the KT3(z;p), see A157704, follow a simple pattern.
A Maple algorithm that generates relevant GFKT4(z;p) information can be found below.

			Some PDGF4 (z;n) are:
  PDGF4(z; n=3) = (1-7*z)*(1-5*z)^4*(1-3*z)^7*(1-z)^10
  PDGF4(z; n=4) = (1-9*z)*(1-7*z)^4*(1-5*z)^7*(1-3*z)^10*(1-z)^13
The first few GFKT4's are:
  GFKT4(z;p=0) = 1/(1-z)
  GFKT4(z;p=1) = -(1+3*z+2*z^2)/(1-z)^4
  GFKT4(z;p=2) = z*(18+128*z+171*z^2+42*z^3+z^4)/(1-z)^7
Some KT4(z,p) polynomials are:
  KT4(z;p=2) = 18+128*z+171*z^2+42*z^3+z^4
  KT4(z;p=3) = 22+1348*z+11738*z^2+26293*z^3+17693*z^4+3271*z^5+115*z^6

Originator sequence A156933.
See A081436 for the z^1 coefficients and A157708 for the z^2 coefficients.
Row sums equal A064350 and those of A157704.
Cf. A157702, A157703, A157704.

p:=2; fn:=sum((-1)^(n1+1)*binomial(3*p+1,n1) *a(n-n1),n1=1..3*p+1): fk:=rsolve(a(n) = fn,a(k)): for n2 from 0 to 3*p+1 do fz(n2):=product((1-(2*n2+1-(2*k))*z)^(3*k+1), k=0..n2): a(n2):= coeff(fz(n2),z,p): end do: b:=n-> a(n): seq(b(n), n=0..3*p+1); a(n)=fn; a(k)= sort (simplify(fk)); GFKT4(p):=sum((fk)*z^k,k=0..infinity); q4:=ldegree((numer (GFKT4(p)))): KT4(p):=sort((-1)^(p)*simplify((GFKT4(p)*(1-z)^(3*p+1))/z^q4),z, ascending);

PDGF4(z;n) = Product_{k=0..n} (1-(2*n+1-2*k)*z)^(3*k+1) with n = 1, 2, 3, ...
GFKT4(z;p) = (-1)^(p)*(z^q4)*KT4(z, p)/(1-z)^(3*p+1) with p = 0, 1, 2, ...
The recurrence relation for the z^p coefficients a(n) is a(n) = Sum_{k=1..3*p+1} (-1)^(k+1)*binomial(3*p + 1, k)*a(n-k) with p = 0, 1, 2, ... .

A156938 Row sums of the FP4 polynomials of A156933.

Original entry on oeis.org

1, 2, -640, -55050240, 2678771102515200, 138029104604220676374528000, -12755118515091679234002154772496384000000, -3319735985081076411622842555075338617937492705280000000000

Offset: 0

Johannes W. Meijer, Feb 20 2009

The row sums of the FP4 have simple prime factors.

Row sums of A156933.

A157708 The z^2 coefficients of the polynomials in the GF4 denominators of A156933.

Original entry on oeis.org

18, 254, 1571, 6335, 19615, 50743, 115234, 237066, 451320, 807180, 1371293, 2231489, 3500861, 5322205, 7872820, 11369668, 16074894, 22301706, 30420615, 40866035, 54143243, 70835699, 91612726

Offset: 1

Johannes W. Meijer, Mar 07 2009

See A157705 for background information.

Cf. A156933, A157705

nmax:=23; for n from 0 to nmax do fz(n):=product((1-(2*n+1-2*k)*z)^(3*k+1), k=0..n); c(n):= coeff(fz(n),z,2); end do: a:=n-> c(n): seq(a(n), n=1..nmax);

a(n) = 7*a(n-1)-21*a(n-2)+35*a(n-3)-35*a(n-4)+21*a(n-5)-7*a(n-6)+a(n-7)
a(n) = 1/2*n^6+5/2*n^5+41/8*n^4+67/12*n^3+27/8*n^2+11/12*n
G.f.: (18 + 128*z + 171*z^2 + 42*z^3 + z^4)/(1-z)^7

A000340 a(0)=1, a(n) = 3*a(n-1) + n + 1.

Original entry on oeis.org

1, 5, 18, 58, 179, 543, 1636, 4916, 14757, 44281, 132854, 398574, 1195735, 3587219, 10761672, 32285032, 96855113, 290565357, 871696090, 2615088290, 7845264891, 23535794695, 70607384108, 211822152348, 635466457069

Offset: 0

N. J. A. Sloane, Simon Plouffe

From Johannes W. Meijer, Feb 20 2009: (Start)
Second right hand column (n-m=1) of the A156920 triangle.
The generating function of this sequence enabled the analysis of the polynomials A156921 and A156925.
(End)
Partial sums of A003462, and thus the second partial sums of A000244 (3^n). Also column k=2 of A106516. - John Keith, Jan 04 2022

			G.f. = 1 + 5*x + 18*x^2 + 58*x^3 + 179*x^4 + 543*x^5 + 1636*x^6 + ...

F. N. David, M. G. Kendall and D. E. Barton, Symmetric Function and Allied Tables, Cambridge, 1966, p. 260.
N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Vincenzo Librandi, Table of n, a(n) for n = 0..1000
Dillan Agrawal, Selena Ge, Jate Greene, Tanya Khovanova, Dohun Kim, Rajarshi Mandal, Tanish Parida, Anirudh Pulugurtha, Gordon Redwine, Soham Samanta, and Albert Xu, Chip-Firing on Infinite k-ary Trees, arXiv:2501.06675 [math.CO], 2025. See pp. 9, 18.
Shaoshi Chen, Hanqian Fang, Sergey Kitaev, and Candice X.T. Zhang, Patterns in Multi-dimensional Permutations, arXiv:2411.02897 [math.CO], 2024. See p. 7.
INRIA Algorithms Project, Encyclopedia of Combinatorial Structures 389
Simon Plouffe, Approximations de séries génératrices et quelques conjectures, Dissertation, Université du Québec à Montréal, 1992; arXiv:0911.4975 [math.NT], 2009.
Simon Plouffe, 1031 Generating Functions, Appendix to Thesis, Montreal, 1992
László Tóth, On Schizophrenic Patterns in b-ary Expansions of Some Irrational Numbers, arXiv:2002.06584 [math.NT], 2020. Mentions this sequence. See also Proc. Amer. Math. Soc. 148 (2020), 461-469.
Index entries for linear recurrences with constant coefficients, signature (5,-7,3).

From Johannes W. Meijer, Feb 20 2009: (Start)
Cf. A156921, A156925, A156927, A156933. Other columns A156922, A156923, A156924.
Equals A156920 second right hand column.
Equals A142963 second right hand column divided by 2^n.
Equals A156919 second right hand column divided by 2.
(End)
Cf. A014915.
Equals column k=1 of A008971 (shifted). - Jeremy Dover, Jul 11 2021
Cf. A000340, A003462 (first differences), A106516.

[(3^(n+2)-2*n-5)/4: n in [0..30]]; // Vincenzo Librandi, Aug 15 2011

a[ -1]:=0:a[0]:=1:for n from 1 to 50 do a[n]:=4*a[n-1]-3*a[n-2]+1 od: seq(a[n],n=0..50); # Miklos Kristof, Mar 09 2005
A000340:=-1/(3*z-1)/(z-1)**2; # conjectured by Simon Plouffe in his 1992 dissertation

a[ n_] := MatrixPower[ {{1, 0, 0}, {1, 1, 0}, {1, 1, 3}}, n + 1][[3, 1]]; (* Michael Somos, May 28 2014 *)
RecurrenceTable[{a[0]==1,a[n]==3a[n-1]+n+1},a,{n,30}] (* or *) LinearRecurrence[{5,-7,3},{1,5,18},30] (* Harvey P. Dale, Jan 31 2017 *)

G.f.: 1/((1-3*x)*(1-x)^2).
a(n) = (3^(n+2) - 2*n - 5)/4.
a(n) = Sum_{k=0..n+1} (n-k+1)*3^k = Sum_{k=0..n+1} k*3^(n-k+1). - Paul Barry, Jul 30 2004
a(n) = Sum_{k=0..n} binomial(n+2, k+2)*2^k. - Paul Barry, Jul 30 2004
a(-1)=0, a(0)=1, a(n) = 4*a(n-1) - 3*a(n-2) + 1. - Miklos Kristof, Mar 09 2005
a(n) = 5*a(n-1) - 7*a(n-2) + 3*a(n-3). - Johannes W. Meijer, Feb 20 2009
a(-2 - n) = 3^-n * A014915(n). - Michael Somos, May 28 2014
E.g.f.: exp(x)*(9*exp(2*x) - 2*x - 5)/4. - Stefano Spezia, Nov 09 2024

A081436 Fifth subdiagonal in array of n-gonal numbers A081422.

Original entry on oeis.org

1, 7, 24, 58, 115, 201, 322, 484, 693, 955, 1276, 1662, 2119, 2653, 3270, 3976, 4777, 5679, 6688, 7810, 9051, 10417, 11914, 13548, 15325, 17251, 19332, 21574, 23983, 26565, 29326, 32272, 35409, 38743, 42280, 46026, 49987, 54169, 58578, 63220

Offset: 0

Paul Barry, Mar 21 2003

One of a family of sequences with palindromic generators.
Also as A(n) = (1/6)*(6*n^3 - 3*n^2 + 3*n), n>0: structured pentagonal diamond numbers (vertex structure 5). (Cf. A004068 = alternate vertex; A000447 = structured diamonds; A100145 for more on structured numbers.) - James A. Record (james.record(AT)gmail.com), Nov 07 2004
Sequence of the absolute values of the z^1 coefficients of the polynomials in the GF4 denominators of A156933. See A157705 for background information. - Johannes W. Meijer, Mar 07 2009
Row 1 of the convolution arrays A213831 and A213833. - Clark Kimberling, Jul 04 2012
Partial sums of A056109. - J. M. Bergot, Jun 22 2013
Number of ordered pairs of intersecting multisets of size 2, each chosen with repetition from {1,...,n}. - Robin Whitty, Feb 12 2014
Row sums of A244418. - L. Edson Jeffery, Jan 10 2015

Vincenzo Librandi, Table of n, a(n) for n = 0..5000
Christian Barrientos, The number of spanning trees of cyclic snakes, Indones. J. Comb. (2025) Vol. 9, No. 1, 21-30. See p. 29.
J. A. Dias da Silva and Pedro J. Freitas, Counting Spectral Radii of Matrices with Positive Entries, arXiv:1305.1139 [math.CO], 2013.
Theorem of the Day, Lovász Local Lemma example involving intersecting pairs of multisets
Index entries for linear recurrences with constant coefficients, signature (4,-6,4,-1).

Cf. A081434, A081435, A081437, A156933, A157705, A244418.

List([0..40], n-> (n+1)*(2*(n+1)^2-n)/2); # G. C. Greubel, Aug 14 2019

[(2*n^3+5*n^2+5*n+2)/2: n in [0..40]]; // Vincenzo Librandi, Jul 19 2011

A081436 := proc(n)
    (n+1)*(2*n^2+3*n+2)/2 ;
end proc:
seq(A081436(n),n=0..60) ; # R. J. Mathar, Jun 26 2013

LinearRecurrence[{4, -6, 4, -1}, {1, 7, 24, 58}, 40] (* Jean-François Alcover, Sep 21 2017 *)

a(n)=n^3+5/2*n*(n+1)+1 \\ Charles R Greathouse IV, Jun 20 2013

[(n+1)*(2*(n+1)^2-n)/2 for n in (0..40)] # G. C. Greubel, Aug 14 2019

a(n) = (n+1)*(2*n^2 + 3*n + 2)/2.
G.f.: (1+x)*(1+2*x)/(1-x)^4. (Convolution of A005408 and A016777.)
a(n) = A110449(n, n-1), for n>1.
a(n) = (n+1)*T(n+1) + n*T(n), where T( ) are triangular numbers. Binomial transform of [1, 6, 11, 6, 0, 0, 0, ...]. - Gary W. Adamson, Dec 28 2007
E.g.f.: exp(x)*(2 + 12*x + 11*x^2 + 2*x^3)/2. - Stefano Spezia, Apr 13 2021
a(n) = 4*a(n-1) - 6*a(n-2) + 4*a(n-3) - a(n-4). - Wesley Ivan Hurt, Apr 14 2021

G.f. simplified and crossrefs added by Johannes W. Meijer, Mar 07 2009

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

A156921 FP1 polynomials related to the generating functions of the right hand columns of the A156920 triangle.

Original entry on oeis.org

1, 1, 1, 1, -6, 1, 7, -79, 119, 126, -270, 1, 28, -515, 1654, 8689, -65864, 142371, -82242, -99090, 113400, 1, 86, -2255, 5784, 300930, -3904584, 20663714, -41517272, -80232259, 657717054

Offset: 0

Johannes W. Meijer, Feb 20 2009

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

			The first few rows of the "triangle" of the coefficients of the FP1 polynomials.
In the columns the coefficients of the powers of z^m, m=0,1,2,... , appear.
  [1]
  [1]
  [1, 1, -6]
  [1, 7, -79, 119, 126, -270]
  [1, 28, -515, 1654, 8689, -65864, 142371, -82242, -99090, 113400]
Matrix of the coefficients of the FP1 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, 0, 0, 0]
  [1, 0 ,0, 0, 0, 0, 0, 0, 0, 0]
  [1, 1, -6, 0 ,0, 0, 0, 0, 0, 0]
  [1, 7, -79, 119, 126, -270, 0, 0, 0, 0]
  [1, 28, -515, 1654, 8689, -65864, 142371, -82242, -99090, 113400]
The first few FP1 polynomials are:
  FP1(z; RHCnr=1) = 1
  FP1(z; RHCnr=2) = 1
  FP1(z; RHCnr =3) = 1+z-6*z^2
Some GF1(z;RHCnr) are:
  GF1(z;RHCnr= 3) = (1+z-6*z^2)/((1-5*z)*(1-3*z)^2*(1-z)^3)
  GF1(z;RHCnr= 4) = (1+7*z-79*z^2+119*z^3+126*z^4-270*z^5)/((1-7*z)*(1-5*z)^2*(1-3*z)^3*(1-z)^4)

Cf. A156920, A156925, A156927, A156933.
For the first few GF1's see A000340, A156922, A156923, A156924.
The number of FP1 terms follow the triangular numbers A000217, with quite surprisingly one exception here a(0)=1.
Abs(Row sums (n)) = A098695(n).
For the polynomials in the denominators of the GF1(z;RHCnr) see A157702.

RHCnr:=4: if RHCnr=1 then RHCmax :=1; else RHCmax:=(RHCnr-1)*(RHCnr)/2 end if: 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+1); for nx from 0 to RHCmax do num:=sort(sum(A[t]*z^t, t=0..RHCmax)); nom:=Product((1-(2*u-1)*z)^(RHCnr-u+1),u=1..RHCnr): RHCa:= series(num/nom,z,nx+1); y:=coeff(RHCa,z,nx)-A[nx]; x:=RHC[nx+1]; A[nx]:=x-y; end do: FP1[RHCnr]:=sort(num,z, ascending); GenFun[RHCnr] :=FP1[RHCnr]/product((1-(2*m-1)*z)^(RHCnr-m+1),m=1..RHCnr);

G.f.: GF1(z;RHCnr) := FP1(z;RHCnr)/product((1-(2*m-1)*z)^(RHCnr+1-m),m=1..RHCnr)

Row sums (n) = (-1)^(1+(n+1)*(n+2)/2)*A098695(n).

A156927 FP3 polynomials related to the generating functions of the columns of the A156921 matrix.

Original entry on oeis.org

1, 1, 1, -6, 29, 31, -283, 245, 298, -286, -108, 119, -3106, 29469, -104585, -220481, 3601363, -15487305, 34949165, -39821950, 4356011, 46881744, -51274736, 9005908, 14663472, -5205168, -1456704, -20736

Offset: 0

Johannes W. Meijer, Feb 20 2009

For the matrix of the FP1 polynomials see A156921. The coefficients in the columns of this matrix are the powers of z^m, m=0, 1, 2, ... . The columns are numbered 1, 2, 3... .
The GF3(z;m) generate the sequences of the z^m coefficients. The general structure of the GF3(z;m) is given below.
The FP3(z,m) in the numerator of the GF3(z;m) is a polynomial of a certain degree, let's say k3. The (k3+1) coefficients of this polynomial can be determined one by one by comparing the series expansion of the FP3(z,m) with the values of the powers of z^m in column (m+1). These values can be generated with the GF1 formulas, see A156921.
An appropriate name for the polynomials FP3(z;m) in the numerators of the GF(3;m) seems to be the flower polynomials of the third kind, the FP3, because the zero patterns of these polynomials look like flowers. The zero patterns of the FP3 and the FP4, see A156933, resemble each other closely and look like the zero patterns of the FP1 and FP2.
The sequence of the (k3+1) number of terms of the FP3(z;m) polynomials for m from 0 to 11 is 1, 2, 8, 17, 29, 45, 63, 84, 109, 137, 167, 200.

			The first few rows of the "triangle" of the FP3(z,m) coefficients are:
  [1]
  [1, 1]
  [-6, 29, 31, -283, 245, 298, -286, -108]
The first few FP3 polynomials are:
  FP3(z; m=0) = 1
  FP3(z; m=1) = (1+z)
  FP3(z; m=2) = (-6+29*z+31*z^2-283*z^3+245*z^4+298*z^5-286*z^6-108*z^7)
Some GF3(z;m) are:
  GF3(z;m=1) = z^2*(1+z)/((1-z)^4*(1-2*z))
  GF3(z;m=2) = z^2*(-6+29*z+31*z^2-283*z^3+245*z^4+298*z^5-286*z^6-108*z^7)/((1-z)^7*(1-2*z)^4*(1-3*z))

Cf. A156920, A156921, A156925, A156933.
For the first few GF3's see A156928, A156929, A156930, A156931.
Row sums A156932.
For the polynomials in the denominators of the GF3(z;m) see A157704.

G.f.: GF3(z;m):= z^p*FP3(z;m)/Product_{k=0..m} (1-(k+1)*z)^(1+3*k).

cp's OEIS Frontend

A157705 G.f.s of the z^p coefficients of the polynomials in the GF4 denominators of A156933.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Maple

Formula

A156938 Row sums of the FP4 polynomials of A156933.

Views

Author

Keywords

Comments

Crossrefs

A157708 The z^2 coefficients of the polynomials in the GF4 denominators of A156933.

Views

Author

Keywords

Comments

Crossrefs

Programs

Maple

Formula

A000340 a(0)=1, a(n) = 3*a(n-1) + n + 1.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Magma

Maple

Mathematica

Formula

A081436 Fifth subdiagonal in array of n-gonal numbers A081422.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

GAP

Magma

Maple

Mathematica

PARI

Sage

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)!.

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)!.