1,3

In A160464 we defined the ES1 matrix by ES1[2*m-1,n=1] and in A094665 it was shown that the nth term of the coefficients of matrix row ES1[1-2*m,n] for m = 1, 2, 3, .. can be generated with the RES1(1-2*m,n) polynomials.

We define the o.g.f.s. of these polynomials by GFRES1(z,1-2*m) = sum(RES1(1-2*m,n)*z^(n-1), n=1..infinity) for m = 1, 2, .. . The general expression of the o.g.f.s. is GFRES1(z,1-2*m) = (-1)*RE(z,1-2*m)/(2*p(m-1)*(z-1)^(m)). The p(m-1) sequence is Gould's sequence A001316.

The coefficients of the RE(z,1-2*m) polynomials lead to the triangle given above.

The E(z,n) = numer(sum((-1)^(n+1)*k^n*z^(k-1),k=1..infinity)) polynomials with n = 1, 2, .. , see the Maple algorithm, lead to the Eulerian numbers A008292.

Some of our results are conjectures based on numerical evidence.

The first few rows are:

[1]

[1]

[2, 1]

[17, 26, 2]

[62, 192, 60, 1]

The first few polynomials RE(z,m) are:

RE(z,-1) = 1

RE(z,-3) = 1

RE(z,-5) = 2+z

RE(z,-7) = 17+26*z+2*z^2

The first few GFRES1(z,m) are:

GFRES1(z,-1) = -(1/1)*(1)/(2*(z-1)^1)

GFRES1(z,-3) = -(1/2)*(1)/(2*(z-1)^2)

GFRES1(z,-5) = -(1/2)*(2+z)/(2*(z-1)^3)

GFRES1(z,-7) = -(1/4)*(17+26*z+2*z^2)/(2*(z-1)^4)

restart; nmax:=8; mmax:=nmax: T(0, x):=1: for i from 1 to nmax do dgr:=degree(T(i-1, x), x): for na from 0 to dgr do c(na):=coeff(T(i-1, x), x, na) od: T(i-1, x+1):=0: for nb from 0 to dgr do T(i-1, x+1):=T(i-1, x+1)+c(nb)*(x+1)^nb od: for nc from 0 to dgr do ECGP(i-1, nc+1):=coeff(T(i-1, x), x, nc) od: T(i, x):=expand((2*x+1)*(x+1)*T(i-1, x+1)- 2*x^2*T(i-1, x)) od: dgr:= degree (T(nmax, x), x): kmax:=nmax: for k from 1 to kmax do p:=k: for m from 1 to k do E(m, k):= sum((-1)^(m-q)*(q^k)*binomial(k+1, m-q), q=1..m) od: fx(p):=(-1)^(p+1)* (sum(E(r, k)*z^(k-r), r=1..k))/(z-1)^(p+1):

GF(-(2*p+1)):= sort(simplify(((-1)^p* 1/2^(p+1)) * sum(ECGP(k-1, k-s)*fx(k-s), s=0..k-1)), z, ascending): numer(-GF(-2*k-1)); NUMGF(-(2*p+1)):=-numer(GF(-(2*p+1))): for n from 1 to mmax+1 do A(k+1, n):= coeff(NUMGF(-(2*p+1)), z, n-1) od: od: for m from 2 to mmax do A(1, m):=0 od: A(1, 1):=1: FT(1):=1: for n from 1 to nmax do for m from 1 to n do FT((n)*(n-1)/2+m+1):=A(n+1, m) end do end do: a:=n-> FT(n):seq(a(n), n = 1..(nmax+1)*(nmax)/2+1);

Cf. A160464, A094665 and A083061.

For the Eulerian numbers E(n, k) see A008292.

The p(n) sequence equals Gould's sequence A001316.

The first right hand column of the triangle equals A048896.

The first left hand column equals A160469.

The row sums equal the absolute values of A117972.

Sequence in context: A012963 A013120 A012968 this_sequence A012889 A013072 A012895

Adjacent sequences: A160465 A160466 A160467 this_sequence A160469 A160470 A160471

easy,nonn,tabf

Johannes W. Meijer (meijgia(AT)hotmail.com), May 24 2009

Maple program fixed by Johannes W. Meijer (meijgia(AT)hotmail.com), Jun 21 2009