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

0, 1, 7, 73, 905, 12354, 179305, 2715192, 42414021, 678476755, 11058588574, 182999237590, 3066447596459, 51926183715280, 887204891847960, 15276037569668880, 264797324173666845, 4617195655522976361, 80930337327794271445, 1425171253004955494215, 25202145191953299213490

Offset

0,3

Comments

Reversion of g.f. for 4-dimensional figurate numbers A001296 (with signs).

Links

Robert Israel, Table of n, a(n) for n = 0..781

Eric Weisstein's World of Mathematics, Series Reversion

Formula

a(n) = (1/n) * Sum_{k=0..n-1} binomial(n+k-1,k) * binomial(5*n,n-k-1) * 2^k for n > 0.

a(n) ~ sqrt(32 - 19*sqrt(5/2)) * 3^(4*n - 3/2) * 5^(3*n) / (sqrt(Pi) * n^(3/2) * 2^(2*n + 3/2) * (25 + 34*sqrt(10))^n). - Vaclav Kotesovec, Sep 27 2023

D-finite with recurrence: (-961083984375*n^4 - 1922167968750*n^3 - 1345517578125*n^2 - 384433593750*n - 36905625000)*a(n) + (31514062500*n^4 + 154153125000*n^3 + 183540937500*n^2 - 24300000000*n - 92156535000)*a(n + 1) + (-23749650000*n^4 - 295611300000*n^3 - 1314962730000*n^2 - 2514101580000*n - 1758020837760)*a(n + 2) + (475256000*n^4 + 6490800000*n^3 + 33692228800*n^2 + 78463987200*n + 68917766400)*a(n + 3) + (42909696*n^4 + 643645440*n^3 + 3599646720*n^2 + 8891842560*n + 8181448704)*a(n + 4) = 0. - Robert Israel, Jan 12 2026

Maple

rec := {(-961083984375*n^4 - 1922167968750*n^3 - 1345517578125*n^2 - 384433593750*n - 36905625000)*a(n) + (31514062500*n^4 + 154153125000*n^3 + 183540937500*n^2 - 24300000000*n - 92156535000)*a(n + 1) + (-23749650000*n^4 - 295611300000*n^3 - 1314962730000*n^2 - 2514101580000*n - 1758020837760)*a(n + 2) + (475256000*n^4 + 6490800000*n^3 + 33692228800*n^2 + 78463987200*n + 68917766400)*a(n + 3) + (42909696*n^4 + 643645440*n^3 + 3599646720*n^2 + 8891842560*n + 8181448704)*a(n + 4), a(0) = 0, a(1) = 1, a(2) = 7, a(3) = 73, a(4) = 905}:
f:= gfun:-rectoproc(rec, a(n), remember):
map(f, [$0..30]); # Robert Israel, Jan 12 2026

Mathematica

nmax = 20; A[_] = 0; Do[A[x_] = x (1 + A[x])^5/(1 - 2 A[x]) + O[x]^(nmax + 1) // Normal, nmax + 1]; CoefficientList[A[x], x]
CoefficientList[InverseSeries[Series[x (1 - 2 x)/(1 + x)^5, {x, 0, 20}], x], x]
Join[{0}, Table[1/n Sum[Binomial[n + k - 1, k] Binomial[5 n, n - k - 1] 2^k, {k, 0, n - 1}], {n, 1, 20}]]

Crossrefs

Cf. A001296, A002294, A064063, A365755, A366014, A366035, A366036, A366037.

Sequence in context: A106651 A114429 A365846 * A124547 A084363 A321837

Adjacent sequences: A365665 A365666 A365667 * A365669 A365670 A365671

Keyword

nonn

Author

Ilya Gutkovskiy, Sep 26 2023

Status

approved