|
%I #25 Jul 31 2024 03:16:50
%S 1,1,2,1,7,5,1,15,32,13,1,26,111,123,34,1,40,285,603,429,89,1,57,610,
%T 2094,2748,1408,233,1,77,1155,5845,12170,11196,4437,610,1,100,2002,
%U 14014,42355,60686,42255,13587,1597,1,126,3246,30030,124137,254756,271961,150951,40736,4181
%N Triangle read by rows: T(n, k) = C(n, k)*Sum_{j=0..n} C(n, k-j)*C(n+j, j)/C(2*j, j).
%H G. C. Greubel, <a href="/A339000/b339000.txt">Rows n = 0..50 of the triangle, flattened</a>
%F G.f.: A008459(x,y)/(1-x*y*A008459(x,y)^2).
%F T(n,n) = Fibonacci(2*n+1).
%e Triangle begins as:
%e 1;
%e 1, 2;
%e 1, 7, 5;
%e 1, 15, 32, 13;
%e 1, 26, 111, 123, 34;
%e 1, 40, 285, 603, 429, 89;
%e 1, 57, 610, 2094, 2748, 1408, 233;
%t T[n_, k_]:= With[{B=Binomial}, B[n,k]*Sum[B[n,k-j]*B[n+j,j]/B[2*j,j], {j,0,n}]];
%t Table[T[n,k], {n,0,12}, {k,0,n}]//Flatten (* _G. C. Greubel_, Jul 31 2024 *)
%o (Maxima)
%o T(n,m):=(binomial(n,m))*sum(((binomial(n,m-k))*(binomial(n+k,k)) )/(binomial(2*k,k)),k,0,n);
%o (Magma)
%o b:=Binomial;
%o A339000:= func< n,k | b(n,k)*(&+[b(n,k-j)*b(n+j,j)/b(2*j,j): j in [0..n]]) >;
%o [A339000(n,k): k in [0..n], n in [0..12]]; // _G. C. Greubel_, Jul 31 2024
%o (SageMath)
%o b=binomial
%o def A339000(n,k): return b(n,k)*sum(b(n,k-j)*b(n+j,j)//b(2*j,j) for j in range(n+1))
%o flatten([[A339000(n,k) for k in range(n+1)] for n in range(13)]) # _G. C. Greubel_, Jul 31 2024
%Y Cf. A000045 (Fibonacci), A001519, A008459, A046748 (row sums).
%K nonn,tabl
%O 0,3
%A _Vladimir Kruchinin_, Nov 18 2020
|
