cp's OEIS Frontend

A294982 Number of compositions (ordered partitions) of 1 into exactly 3n+1 powers of 1/(n+1).

%I A294982 #17 Sep 20 2019 04:06:15
%S A294982 1,13,217,4245,90376,2019836,46570140,1097525253,26293568950,
%T A294982 638048716305,15643738390215,386826618273420,9633468179090952,
%U A294982 241366000080757480,6078975012187601768,153798067122829610085,3906583987216447704594,99579591801208823965115
%N A294982 Number of compositions (ordered partitions) of 1 into exactly 3n+1 powers of 1/(n+1).
%H A294982 Alois P. Heinz, <a href="/A294982/b294982.txt">Table of n, a(n) for n = 0..700</a>
%F A294982 For n>0, a(n) = binomial(3*n+1, n+1)*binomial(2*n, n) + binomial(3*n+1, n-1). - _Vaclav Kotesovec_, Sep 20 2019
%e A294982 a(0) = 1: [1].
%e A294982 a(1) = 13: [1/4,1/4,1/4,1/4], [1/2,1/4,1/8,1/8], [1/2,1/8,1/4,1/8], [1/2,1/8,1/8,1/4], [1/4,1/2,1/8,1/8], [1/4,1/8,1/2,1/8], [1/4,1/8,1/8,1/2], [1/8,1/2,1/4,1/8], [1/8,1/2,1/8,1/4], [1/8,1/4,1/2,1/8], [1/8,1/4,1/8,1/2], [1/8,1/8,1/2,1/4], [1/8,1/8,1/4,1/2].
%p A294982 a:= proc(n) option remember; `if`(n<2, 12*n+1, (3*n-1)*(3*n+1)*
%p A294982       3*((15*n^3-31*n^2-4*n+8)*n*a(n-1)-3*(3*n-4)*(3*n-2)*
%p A294982       (3*n^2-2*n-2)*a(n-2))/((n+1)*(4*n+2)*(3*n^2-8*n+3)*n^2))
%p A294982     end:
%p A294982 seq(a(n), n=0..20);
%t A294982 a[n_] := a[n] = If[n < 2, 12*n + 1, (3*n - 1)*(3*n + 1)*3*((15*n^3 - 31*n^2 - 4*n + 8)*n*a[n-1] - 3*(3*n - 4)*(3*n - 2)*(3*n^2 - 2*n - 2)*a[n-2])/((n + 1)*(4*n + 2)*(3*n^2 - 8*n + 3)*n^2)];
%t A294982 Table[a[n], {n, 0, 20}] (* _Jean-François Alcover_, May 21 2018, translated from Maple *)
%t A294982 Table[Binomial[3*n + 1, n + 1]*Binomial[2*n, n] + Binomial[3*n + 1, n - 1], {n, 0, 20}] (* _Vaclav Kotesovec_, Sep 20 2019 *)
%Y A294982 Row n=3 of A294746.
%K A294982 nonn
%O A294982 0,2
%A A294982 _Alois P. Heinz_, Nov 12 2017