login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A156052 Triangle read by rows: T(n, k) = binomial(n, k)/Beta(n+1, n-k+1) + binomial(n, n-k)/Beta(n+1, k+1). 1

%I #7 Sep 08 2022 08:45:41

%S 2,8,8,33,48,33,144,240,240,144,635,1240,1260,1240,635,2778,6510,6720,

%T 6720,6510,2778,12019,33600,38430,33600,38430,33600,12019,51488,

%U 168672,223776,184800,184800,223776,168672,51488,218799,824400,1275120,1119888,900900,1119888,1275120,824400,218799

%N Triangle read by rows: T(n, k) = binomial(n, k)/Beta(n+1, n-k+1) + binomial(n, n-k)/Beta(n+1, k+1).

%C Row sums are 2*A108666(n+1): {2, 16, 114, 768, 5010, 32016, 201698, 1257472, 7777314, 47800080, 292292946, ...}.

%H G. C. Greubel, <a href="/A156052/b156052.txt">Rows n = 0..100 of triangle, flattened</a>

%F T(n, k) = binomial(n, k)/Beta(n+1, n-k+1) + binomial(n, n-k)/Beta(n+1, k+1).

%F T(n, k) = (k+1)*binomial(n+1, k+1)*( binomial(2*n-k+1, n+1) + binomial(n+k+1, n+1) ). - _G. C. Greubel_, Dec 01 2019

%e Triangle begins as:

%e 2;

%e 8, 8;

%e 33, 48, 33;

%e 144, 240, 240, 144;

%e 635, 1240, 1260, 1240, 635;

%e 2778, 6510, 6720, 6720, 6510, 2778;

%e 12019, 33600, 38430, 33600, 38430, 33600, 12019;

%e 51488, 168672, 223776, 184800, 184800, 223776, 168672, 51488;

%p b:=binomial; seq(seq( (k+1)*b(n+1, k+1)*( b(2*n-k+1, n+1) + b(n+k+1, n+1) ), k=0..n), n=0..10); # _G. C. Greubel_, Dec 01 2019

%t Table[Binomial[n, k]/Beta[n+1, n-k+1] + Binomial[n, n-k]/Beta[n+1, k+1], {n, 0, 10}, {k, 0, n}]//FlattenTable[(k+1)*Binomial[n+1, k+1]*(Binomial[n+k+1, n+1] + Binomial[2*n-k+1, n+1]), {n, 0, 10}, {k, 0, n}]//Flatten (* _G. C. Greubel_, Dec 01 2019 *)

%o (PARI) T(n, k) = my(b=binomial); (k+1)*b(n+1, k+1)*( b(2*n-k+1, n+1) + b(n+k+1, n+1) ); \\ _G. C. Greubel_, Dec 01 2019

%o (Magma) B:=Binomial; [(k+1)*B(n+1, k+1)*( B(2*n-k+1, n+1) + B(n+k+1, n+1) ): k in [0..n], n in [0..10]]; // _G. C. Greubel_, Dec 01 2019

%o (Sage) b=binomial; [[(k+1)*b(n+1, k+1)*( b(2*n-k+1, n+1) + b(n+k+1, n+1) ) for k in (0..n)] for n in (0..10)] # _G. C. Greubel_, Dec 01 2019

%o (GAP) B:=Binomial;; Flat(List([0..10], n-> List([0..n], k-> (k+1)*B(n+1, k+1)*( B(2*n-k+1, n+1) + B(n+k+1, n+1) ) ))); # _G. C. Greubel_, Dec 01 2019

%K nonn,tabl

%O 0,1

%A _Roger L. Bagula_, Feb 02 2009

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified April 19 08:08 EDT 2024. Contains 371782 sequences. (Running on oeis4.)