A291452 - OEIS

%I #12 Jul 23 2019 04:23:52

%S 1,0,1,0,1,35,0,1,495,5775,0,1,8255,450450,2627625,0,1,130815,

%T 35586525,727476750,2546168625,0,1,2098175,2941884000,181262956875,

%U 1932541986375,4509264634875

%N Triangle read by rows, expansion of e.g.f. exp(x*(cos(z) + cosh(z) - 2)/2), nonzero coefficients of z.

%e Triangle starts:

%e [1]

%e [0, 1]

%e [0, 1,      35]

%e [0, 1,     495,       5775]

%e [0, 1,    8255,     450450,      2627625]

%e [0, 1,  130815,   35586525,    727476750,    2546168625]

%e [0, 1, 2098175, 2941884000, 181262956875, 1932541986375, 4509264634875]

%p CL := (f,x) -> PolynomialTools:-CoefficientList(f,x):

%p A291452_row := proc(n) exp(x*(cos(z)+cosh(z)-2)/2):

%p series(%, z, 88): CL((4*n)!*coeff(series(%,z,4*(n+1)),z,4*n),x) end:

%p for n from 0 to 7 do A291452_row(n) od;

%p # Alternative:

%p A291452row := proc(n) local P; P := proc(m, n) option remember;

%p if n = 0 then 1 else add(binomial(m*n, m*k)*P(m, n-k)*x, k=1..n) fi end:

%p CL(P(4, n), x); seq(%[k+1]/k!, k=0..n) end: # _Peter Luschny_, Sep 03 2018

%t P[m_, n_] := P[m, n] = If[n == 0, 1, Sum[Binomial[m*n, m*k]*P[m, n - k]*x, {k, 1, n}]];

%t row[n_] := Module[{cl = CoefficientList[P[4, n], x]}, Table[cl[[k + 1]]/k!, {k, 0, n}]];

%t Table[row[n], {n, 0, 6}] // Flatten (* _Jean-François Alcover_, Jul 23 2019, after _Peter Luschny_ *)

%Y Cf. A048993 (m=1), A156289 (m=2), A291451 (m=3), this seq. (m=4).

%Y Diagonal: A000012 (m=1), A001147 (m=2), A025035 (m=3), A025036 (m=4).

%Y Row sums: A000110 (m=1), A005046 (m=2), A291973 (m=3), A291975 (m=4).

%Y Alternating row sums: A000587 (m=1), A260884 (m=2), A291974 (m=3), A291976 (m=4).

%K nonn,tabl

%O 0,6

%A _Peter Luschny_, Sep 07 2017