OFFSET
0,3
COMMENTS
See A193842 for the definition of fission of two sequences of polynomials or triangular arrays.
FORMULA
From Peter Bala, Jul 16 2013: (Start)
T(n,k) = sum {i = 0..k} (-1)^k*binomial(n+1,k-i)*(-3)^(k-i) for 0 <= k <= n.
O.g.f.: 1/( (1 - 2*x*t)*(1 - (3*x + 1)*t) )= 1 + (1 + 5*x)*t + (1 + 8*x + 19*x^2)*t^2 + .... Cf. A193860.
The n-th row polynomial R(n,x) = 1/(x + 1)*( (3*x + 1)^(n+1) - (2*x)^(n+1) ). (End)
T(n, k) = 3^k*binomial(n+1, k)*hypergeom([1, -k], [n-k+2], 1/3). - Peter Luschny, Nov 19 2018
EXAMPLE
First six rows:
1
1...5
1...8....19
1...11...43....65
1...14...76....194...211
1...17...118...422...793...665
MAPLE
T := (n, k) -> 3^k*binomial(n+1, k)*hypergeom([1, -k], [n-k+2], 1/3):
for n from 0 to 6 do seq(simplify(T(n, k)), k=0..n) od; # Peter Luschny, Nov 19 2018
MATHEMATICA
z = 10;
p[n_, x_] := (2 x + 1)^n;
q[n_, x_] := (x + 1)^n;
p1[n_, k_] := Coefficient[p[n, x], x^k];
p1[n_, 0] := p[n, x] /. x -> 0;
d[n_, x_] := Sum[p1[n, k]*q[n - 1 - k, x], {k, 0, n - 1}]
h[n_] := CoefficientList[d[n, x], {x}]
TableForm[Table[Reverse[h[n]], {n, 0, z}]]
Flatten[Table[Reverse[h[n]], {n, -1, z}]] (* A193856 *)
TableForm[Table[h[n], {n, 0, z}]]
Flatten[Table[h[n], {n, -1, z}]] (* A193857 *)
CROSSREFS
KEYWORD
nonn,tabl
AUTHOR
Clark Kimberling, Aug 07 2011
STATUS
approved