|
%I #12 Jul 24 2023 01:24:49
%S 0,0,0,1,1,1,0,0,0,0,1,2,3,4,5,1,1,1,1,1,1,1,4,9,16,25,36,49,2,4,6,8,
%T 10,12,14,16,2,9,28,65,126,217,344,513,730,3,12,27,48,75,108,147,192,
%U 243,300,4,22,90,268,640,1314,2422,4120,6588,10030,14674
%N Triangle T(n,k) = coefficient [x^n] of x^2/(1-(k+1)*x^2-x^3) for row n, and columns k = 0..n, read by rows.
%H Nathaniel Johnston, <a href="/A117724/b117724.txt">Rows n = 0..50, flattened</a>
%F T(n,k) = coefficient [x^n] ( x^2/(1-(k+1)*x^2-x^3) ).
%F T(n, 0) = A000931(n+1).
%F T(n, 1) = A008346(n-2) = (-1)^(n-1)*A119282(n-1).
%F T(n, 2) = A052931(n-2).
%e The table starts:
%e 0;
%e 0, 0;
%e 1, 1, 1;
%e 0, 0, 0, 0;
%e 1, 2, 3, 4, 5;
%e 1, 1, 1, 1, 1, 1;
%e 1, 4, 9, 16, 25, 36, 49;
%e 2, 4, 6, 8, 10, 12, 14, 16;
%e 2, 9, 28, 65, 126, 217, 344, 513, 730;
%e 3, 12, 27, 48, 75, 108, 147, 192, 243, 300;
%p t:=taylor(x^2/(1-(k+1)*x^2-x^3), x, 15):
%p seq(seq(coeff(t,x,n), k=0..n),n=0..12); # _Nathaniel Johnston_, Apr 27 2011
%t T[n_, k_]:= T[n, k]= Coefficient[Series[x^2/(1-(k+1)*x^2-x^3), {x,0,n+ 2}], x, n];
%t Table[T[n,k], {n,0,12}, {k,0,n}]//Flatten
%o (Magma)
%o m:=12;
%o R<x>:=PowerSeriesRing(Integers(), m+2);
%o A117724:= func< n, k | Coefficient(R!( x^2/(1-(k+1)*x^2-x^3) ), n) >;
%o [A117724(n, k): k in [0..n], n in [0..m]]; // _G. C. Greubel_, Jul 23 2023
%o (SageMath)
%o def A117724(n, k):
%o P.<x> = PowerSeriesRing(QQ)
%o return P( x^2/(1-(k+1)*x^2-x^3) ).list()[n]
%o flatten([[A117724(n, k) for k in range(n+1)] for n in range(13)]) # _G. C. Greubel_, Jul 23 2023
%Y Cf. A000931, A008346, A052931, A117716, A119282.
%K nonn,tabl,easy
%O 0,12
%A _Roger L. Bagula_, Apr 13 2006
%E Sign in definition corrected, offset set to -1 by Assoc. Eds. of the OEIS, Jun 15 2010
%E Edited by _G. C. Greubel_, Jul 23 2023
|
