OFFSET
0,8
LINKS
G. C. Greubel, Rows n = 0..50 of the triangle, flattened
EXAMPLE
Triangle begins as:
0;
0, 0;
1, 1, 1;
1, 2, 3, 4;
1, 4, 9, 16, 25;
2, 9, 28, 65, 126, 217;
3, 20, 87, 264, 635, 1308, 2415;
4, 44, 270, 1072, 3200, 7884, 16954, 32960;
MAPLE
A117716 := proc(n, m)
x^2/(1-(m+1)*x-x^3) ;
if n < 0 then
0;
else
coeftayl(%, x=0, n) ;
end if;
end proc: # R. J. Mathar, May 14 2013
MATHEMATICA
T[n_, k_]:= T[n, k]= Coefficient[Series[x^2/(1-(k+1)*x-x^3), {x, 0, n+ 2}], x, n];
Table[T[n, k], {n, 0, 12}, {k, 0, n}]//Flatten
PROG
(Magma)
m:=12;
R<x>:=PowerSeriesRing(Integers(), m+2);
A117716:= func< n, k | Coefficient(R!( x^2/(1-(k+1)*x-x^3) ), n) >;
[[A117716(n, k): k in [0..n]]: n in [0..m]]; // G. C. Greubel, Jul 23 2023
(SageMath)
def A117716(n, k):
P.<x> = PowerSeriesRing(QQ)
return P( x^2/(1-(k+1)*x-x^3) ).list()[n]
flatten([[A117716(n, k) for k in range(n+1)] for n in range(13)]) # G. C. Greubel, Jul 23 2023
CROSSREFS
KEYWORD
nonn,tabl
AUTHOR
Roger L. Bagula, Apr 13 2006, corrected Apr 15 2006
EXTENSIONS
Edited by G. C. Greubel, Jul 23 2023
STATUS
approved