cp's OEIS Frontend

A303060 G.f. A(x) satisfies: [x^(n-1)] (1+x)^(n^2) / A(x)^n = 0 for n>1.

%I A303060 #14 Oct 06 2020 07:22:19
%S A303060 1,2,2,8,94,1754,44638,1430912,55127738,2474943058,126686496842,
%T A303060 7274993163208,462851928538166,32302608870561034,2452988579883365414,
%U A303060 201318961408580839648,17755213859921917042610,1674529624268435698338210,168166630871054819006126610,17916303403273387932439315944,2018304444655346718894816153102
%N A303060 G.f. A(x) satisfies: [x^(n-1)] (1+x)^(n^2) / A(x)^n = 0 for n>1.
%H A303060 Paul D. Hanna, <a href="/A303060/b303060.txt">Table of n, a(n) for n = 0..300</a>
%F A303060 a(n) = A303061(n) + A303061(n-1) for n>=0.
%F A303060 a(n) ~ c * d^n * n! / n^2, where d = -4 / (LambertW(-2*exp(-2)) * (2 + LambertW(-2*exp(-2)))) = 6.17655460948348035823168... and c = 0.06049920104... - _Vaclav Kotesovec_, Oct 06 2020
%e A303060 G.f.: A(x) = 1 + 2*x + 2*x^2 + 8*x^3 + 94*x^4 + 1754*x^5 + 44638*x^6 + 1430912*x^7 + 55127738*x^8 + 2474943058*x^9 + 126686496842*x^10 + ...
%e A303060 ILLUSTRATION OF DEFINITION.
%e A303060 The table of coefficients in (1+x)^(n^2) / A(x)^n begins:
%e A303060 n=1: [1, -1, 0, -6, -74, -1500, -39688, -1302742, ...];
%e A303060 n=2: [1, 0, -2, -12, -159, -3136, -82180, -2680752, ...];
%e A303060 n=3: [1, 3, 0, -26, -300, -5454, -137764, -4398210, ...];
%e A303060 n=4: [1, 8, 24, 0, -548, -9576, -223760, -6847536, ...];
%e A303060 n=5: [1, 15, 100, 350, 0, -16022, -376660, -10771830, ...];
%e A303060 n=6: [1, 24, 270, 1844, 7641, 0, -596908, -17643792, ...];
%e A303060 n=7: [1, 35, 588, 6258, 46186, 224196, 0, -26940146, ...];
%e A303060 n=8: [1, 48, 1120, 16864, 182640, 1478160, 8281968, 0, ...]; ...
%e A303060 in which the main diagonal equals all zeros after the initial term, illustrating that [x^(n-1)] (1+x)^(n^2) / A(x)^n = 0 for n>1.
%e A303060 RELATED SERIES.
%e A303060 Define B(x) by [x^(n-1)] A(x)^n / B(x)^(n^2) = 0, then B(x) is an integer series that begins
%e A303060 B(x) = 1 + x + x^2 + x^3 + 33*x^4 + 161*x^5 + 11061*x^6 + 97837*x^7 + 10690245*x^8 + 130880301*x^9 + 20319470589*x^10 + 308187480741*x^11 + 63533886021453*x^12 + 1128112773827877*x^13 + 294884291428281145*x^14 + 5933604178561639785*x^15 + ...
%o A303060 (PARI) {a(n) = my(A=[1]); for(m=1,n+1, A=concat(A,0); A[m] = Vec( (1+x +x*O(x^n))^(m^2)/Ser(A)^m )[m]/m ); A[n+1]}
%o A303060 for(n=0,30, print1(a(n),", "))
%Y A303060 Cf. A303061.
%K A303060 nonn
%O A303060 0,2
%A A303060 _Paul D. Hanna_, Apr 17 2018