|
%I #15 Sep 08 2022 08:45:54
%S 1,1,6,41,362,3542,37692,424377,4990722,60704138,758665388,9694652838,
%T 126203947828,1668947978908,22370427181624,303383342784729,
%U 4156846359584754,57473870722327874,801081711581734764
%N a(n) = Sum_{k=0..n} binomial(n,k)^4/(n+1).
%H G. C. Greubel, <a href="/A178824/b178824.txt">Table of n, a(n) for n = 0..830</a>
%F a(n) = A005260(n)/(n+1).
%p a:=n->add(binomial(n,k)^4/(n+1),k=0..n): seq(a(n),n=0..20); # _Muniru A Asiru_, Jan 22 2019
%t Table[Sum[Binomial[n,k]^4/(n+1), {k,0,n}], {n,0,20}] (* _G. C. Greubel_, Jan 22 2019 *)
%o (PARI) {a(n)=sum(k=0, n, binomial(n, k)^4)/(n+1)}
%o (Magma) [(&+[Binomial(n,k)^4/(n+1): k in [0..n]]): n in [0..20]]; // _G. C. Greubel_, Jan 22 2019
%o (Sage) [sum(binomial(n,k)^4/(n+1) for k in (0..n)) for n in (0..20)] # _G. C. Greubel_, Jan 22 2019
%o (GAP) List([0..20], n-> Sum([0..n], k-> Binomial(n,k)^4/(n+1) )); # _G. C. Greubel_, Jan 22 2019
%Y Cf. A005260.
%K nonn
%O 0,3
%A _Paul D. Hanna_, Dec 27 2010
|
