OFFSET
0,3
COMMENTS
See A193722 for the definition of fusion of two sequences of polynomials or triangular arrays.
Triangle T(n,k), read by rows, given by (1,1,0,0,0,0,0,0,0,...) DELTA (2,3,0,0,0,0,0,0,0,...) where DELTA is the operator defined in A084938. - Philippe Deléham, Oct 05 2011
LINKS
G. C. Greubel, Rows n = 0..50 of the triangle, flattened
FORMULA
T(n,k) = 5*T(n-1,k-1) + 2*T(n-1,k) with T(0,0)=T(1,0)=1 and T(1,1)=2. - Philippe Deléham, Oct 05 2011
G.f.: (1-x-3*x*y)/(1-2*x-5*x*y). - R. J. Mathar, Aug 11 2015
From G. C. Greubel, Dec 02 2023: (Start)
T(n, 0) = A011782(n).
T(n, n) = A020699(n).
T(n, n-1) = A081040(n-1).
Sum_{k=0..n} T(n, k) = A169634(n-1) + (4/7)*[n=0].
EXAMPLE
First six rows:
1;
1, 2;
2, 9, 10;
4, 28, 65, 50;
8, 76, 270, 425, 250;
16, 192, 920, 2200, 2625, 1250;
MATHEMATICA
(* First program *)
z = 8; a = 1; b = 2; c = 1; d = 2;
p[n_, x_] := (a*x + b)^n ; q[n_, x_] := (c*x + d)^n
t[n_, k_] := Coefficient[p[n, x], x^k]; t[n_, 0] := p[n, x] /. x -> 0;
w[n_, x_] := Sum[t[n, k]*q[n + 1 - k, x], {k, 0, n}]; w[-1, x_] := 1
g[n_] := CoefficientList[w[n, x], {x}]
TableForm[Table[Reverse[g[n]], {n, -1, z}]]
Flatten[Table[Reverse[g[n]], {n, -1, z}]] (* A193726 *)
TableForm[Table[g[n], {n, -1, z}]]
Flatten[Table[g[n], {n, -1, z}]] (* A193727 *)
(* Second program *)
T[n_, k_]:= T[n, k]= If[k<0 || k>n, 0, If[n<2, k+1, 2*T[n-1, k] + 5*T[n -1, k-1]]];
Table[T[n, k], {n, 0, 12}, {k, 0, n}]//Flatten (* G. C. Greubel, Dec 02 2023 *)
PROG
(Magma)
function T(n, k) // T = A193726
if k lt 0 or k gt n then return 0;
elif n lt 2 then return k+1;
else return 2*T(n-1, k) + 5*T(n-1, k-1);
end if;
end function;
[T(n, k): k in [0..n], n in [0..12]]; // G. C. Greubel, Dec 02 2023
(SageMath)
def T(n, k): # T = A193726
if (k<0 or k>n): return 0
elif (n<2): return k+1
else: return 2*T(n-1, k) + 5*T(n-1, k-1)
flatten([[T(n, k) for k in range(n+1)] for n in range(13)]) # G. C. Greubel, Dec 02 2023
CROSSREFS
KEYWORD
nonn,tabl
AUTHOR
Clark Kimberling, Aug 04 2011
STATUS
approved