cp's OEIS Frontend

A349332 G.f. A(x) satisfies A(x) = 1 + x * A(x)^5 / (1 - x).

%I A349332 #25 Nov 23 2024 09:09:24
%S A349332 1,1,6,46,406,3901,39627,418592,4551672,50610692,572807157,6577068383,
%T A349332 76426719408,897078662538,10620634999318,126676885170703,
%U A349332 1520759193166329,18361269213121164,222814883564042704,2716125963857227904,33244557641365865109
%N A349332 G.f. A(x) satisfies A(x) = 1 + x * A(x)^5 / (1 - x).
%H A349332 Seiichi Manyama, <a href="/A349332/b349332.txt">Table of n, a(n) for n = 0..897</a>
%F A349332 a(n) = Sum_{k=0..n} binomial(n-1,k-1) * binomial(5*k,k) / (4*k+1).
%F A349332 a(n) ~ 3381^(n + 1/2) / (25 * sqrt(Pi) * n^(3/2) * 2^(8*n + 7/2)). - _Vaclav Kotesovec_, Nov 15 2021
%F A349332 Recurrence: 8*n*(2*n - 1)*(4*n - 1)*(4*n + 1)*a(n) = (4405*n^4 - 10346*n^3 + 9575*n^2 - 4354*n + 840)*a(n-1) - 12*(n-2)*(1255*n^3 - 3957*n^2 + 4492*n - 1820)*a(n-2) + 2*(n-3)*(n-2)*(10655*n^2 - 32733*n + 26908)*a(n-3) - 4*(n-4)*(n-3)*(n-2)*(3445*n - 6986)*a(n-4) + 3381*(n-5)*(n-4)*(n-3)*(n-2)*a(n-5). - _Vaclav Kotesovec_, Nov 17 2021
%p A349332 a:= n-> coeff(series(RootOf(1+x*A^5/(1-x)-A, A), x, n+1), x, n):
%p A349332 seq(a(n), n=0..20);  # _Alois P. Heinz_, Nov 15 2021
%t A349332 nmax = 20; A[_] = 0; Do[A[x_] = 1 + x A[x]^5/(1 - x) + O[x]^(nmax + 1) // Normal, nmax + 1]; CoefficientList[A[x], x]
%t A349332 Table[Sum[Binomial[n - 1, k - 1] Binomial[5 k, k]/(4 k + 1), {k, 0, n}], {n, 0, 20}]
%o A349332 (PARI) {a(n) = my(A=[1]); for(m=1, n, A=concat(A, 0);
%o A349332 A[#A] = 1 + sum(k=1, m-1, (polcoeff(Ser(A)^5, k)) )); A[n+1]}
%o A349332 for(n=0, 30, print1(a(n), ", ")) \\ _Vaclav Kotesovec_, Nov 23 2024, after _Paul D. Hanna_
%Y A349332 Cf. A002212, A307678, A349331, A349333, A349334, A349335.
%Y A349332 Cf. A002294, A346647 (partial sums), A349361, A364748.
%K A349332 nonn
%O A349332 0,3
%A A349332 _Ilya Gutkovskiy_, Nov 15 2021