|
|
A006256
|
|
a(n) = Sum_{k=0..n} binomial(3k,k)*binomial(3n-3k,n-k).
(Formerly M4229)
|
|
11
|
|
|
1, 6, 39, 258, 1719, 11496, 77052, 517194, 3475071, 23366598, 157206519, 1058119992, 7124428836, 47983020624, 323240752272, 2177956129818, 14677216121871, 98923498131762, 666819212874501, 4495342330033938, 30308036621747679, 204356509814519712
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,2
|
|
COMMENTS
|
The right-hand sides of several of the "Ruehr identities". - N. J. A. Sloane, Feb 20 2020
|
|
REFERENCES
|
Allouche, J-P. "Two binomial identities of Ruehr Revisited." The American Mathematical Monthly 126.3 (2019): 217-225.
Alzer, Horst, and Helmut Prodinger. "On Ruehr's Identities." Ars Comb. 139 (2018): 247-254.
Bai, Mei, and Wenchang Chu. "Seven equivalent binomial sums." Discrete Mathematics 343.2 (2020): 111691.
M. Petkovsek et al., A=B, Peters, 1996, p. 165.
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
|
|
LINKS
|
D. Merlini, R. Sprugnoli and M. C. Verri, The tennis ball problem, J. Combin. Theory, A 99 (2002), 307-344 (Y_n for s=3).
|
|
FORMULA
|
a(n) = (3/4)*(27/4)^n*(1+c/sqrt(n)+o(n^(-1/2))) where c = (2/3)*sqrt(1/(3*Pi)) = 0.217156671956853298... More generally, a(n, m) = sum(k=0, n, C(m*k,k) *C(m*(n-k),n-k)) is asymptotic to (1/2)*m/(m-1)*(m^m/(m-1)^(m-1))^n. See A000302, A078995 for cases m=2 and 4. - Benoit Cloitre, Jan 26 2003, extended by Vaclav Kotesovec, Nov 06 2012
G.f.: 1/(1-3*z*g^2)^2, where g=g(z) is given by g=1+z*g^3, g(0)=1, i.e. (in Maple command) g := 2*sin(arcsin(3*sqrt(3*z)/2)/3)/sqrt(3*z). - Emeric Deutsch, May 22 2003
D-finite with recurrence: 6*(36*n^2-45*n+16)*a(n-1) - 81*(3*n-4)*(3*n-2)*a(n-2) - 8*n*(2*n-1)*a(n) = 0. - Vaclav Kotesovec, Oct 05 2012
From Rui Duarte and António Guedes de Oliveira, Feb 17 2013: (Start)
a(n) = sum(k=0, n, C(3*k+l,k)*C(3*(n-k)-l,n-k)) for every real number l.
a(n) = sum(k=0, n, 2^(n-k)*C(3n+1,k)).
a(n) = sum(k=0, n, 3^(n-k)*c(2n+k,k)). (End)
From Akalu Tefera, Sean Meehan, Michael Weselcouch, and Aklilu Zeleke, May 11 2013: (Start)
a(n) = sum(k=0, 2n, (-3)^k*C(3n - k, n)).
a(n) = sum(k=0, 2n, (-4)^k*C(3n + 1, 2n - k)).
a(n) = sum(k=0, n, 3^k*C(3n - k, 2n)).
a(n) = sum(k=0, n, 2^k*C(3n + 1, n - k)). (End)
a(n) = C(3*n+1,n)*Hyper2F1(1,-n,2*n+2,-2). - Peter Luschny, May 19 2015
|
|
MAPLE
|
a:= proc(n) option remember; `if`(n<2, 5*n+1,
((216*n^2-270*n+96) *a(n-1)
-81*(3*n-2)*(3*n-4) *a(n-2)) /(n*(16*n-8)))
end:
|
|
MATHEMATICA
|
a[n_] := HypergeometricPFQ[{1/3, 2/3, 1/2-n, -n}, {1/2, 1/3-n, 2/3-n}, 1]*(3n)!/(n!*(2n)!); Table[a[n], {n, 0, 20}] (* Jean-François Alcover, Jun 20 2012 *)
Table[Sum[Binomial[3k, k]Binomial[3n-3k, n-k], {k, 0, n}], {n, 0, 30}] (* Harvey P. Dale, Oct 23 2013 *)
|
|
PROG
|
(Haskell)
a006256 n = a006256_list !! n
a006256_list = f (tail a005809_list) [1] where
f (x:xs) zs = (sum $ zipWith (*) zs a005809_list) : f xs (x : zs)
(Sage)
a = lambda n: binomial(3*n+1, n)*hypergeometric([1, -n], [2*n+2], -2)
(Magma) [&+[Binomial(3*k, k) *Binomial(3*n-3*k, n-k): k in [0..n]]:n in [0..22]]; // Vincenzo Librandi, Feb 21 2020
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,easy,nice
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|