A303558 - OEIS

%I #15 Apr 27 2018 12:38:53

%S 1,1,3,15,93,620,4046,24077,138709,1448110,32196907,697062429,

%T 11628586799,140701747704,1213958140794,21010393015426,

%U 1219119257782371,56771570322826701,1958328935823726853,55307040650361330286,1429761229872782798864,40335112020127634181504,1453061157501054116976825,63044483568580391515906296

%N G.f. A(x) satisfies: 0 = [x^(n-1)] (x*A(x))' / (1 + x*A(x)^n)^n for n>1.

%C Compare to: 0 = [x^(n-1)] (x*F(x))' / (1 + x*F(x)^2)^n for n>1 holds when F(x) = 1 + x*F(x)^2 is a g.f. of the Catalan numbers (A000108).

%C Compare to: 0 = [x^(n-1)] (x*G(x))' / (1 + x*G(x)^k)^n for n>1 holds when G(x) = 1 + x*G(x)^k and k is fixed; this sequence explores the case where k varies directly with n.

%C Related identity: 0 = [x^(n-1)] (x*F(x))' / F(x)^n for n>1 holds when F(0) = 1.

%H Paul D. Hanna, <a href="/A303558/b303558.txt">Table of n, a(n) for n = 0..400</a>

%e G.f.: A(x) = 1 + x + 3*x^2 + 15*x^3 + 93*x^4 + 620*x^5 + 4046*x^6 + 24077*x^7 + 138709*x^8 + 1448110*x^9 + 32196907*x^10 + ...

%e such that 0 = [x^(n-1)] (A(x) + x*A'(x)) / (1 + x*A(x)^n)^n for n>1.

%e RELATED SERIES.

%e (x*A(x))' = 1 + 2*x + 9*x^2 + 60*x^3 + 465*x^4 + 3720*x^5 + 28322*x^6 + 192616*x^7 + 1248381*x^8 + ...

%e A'(x)/A(x) = 1 + 5*x + 37*x^2 + 305*x^3 + 2516*x^4 + 19205*x^5 + 126624*x^6 + 792081*x^7 + 10741078*x^8 + ...

%e ILLUSTRATION OF DEFINITION.

%e The table of coefficients of x^k in (x*A(x))'/(1 + x*A(x)^n)^n begins:

%e n=1: [1, 1, 7, 49, 391, 3151, 23815, 158425, 1013803, ...];

%e n=2: [1, 0, 4, 34, 287, 2354, 17532, 110690, 681487, ...];

%e n=3: [1, -1, 0, 17, 169, 1431, 10215, 55203, 297810, ...];

%e n=4: [1, -2, -5, 0, 62, 554, 3111, 1872, -56195, 5924106, ...];

%e n=5: [1, -3, -11, -15, 0, -31, -2145, -36993, -281119, 4826695, ...];

%e n=6: [1, -4, -18, -26, 26, 0, -3742, -50288, -295155, 5468804, ...];

%e n=7: [1, -5, -26, -31, 192, 1053, 0, -34031, -94648, 7795307, ...];

%e n=8: [1, -6, -35, -28, 559, 3620, 10228, 0, 155958, 10425666, ...];

%e n=9: [1, -7, -45, -15, 1197, 8283, 27059, 12493, 0, 9907761, ...];

%e n=10: [1, -8, -56, 10, 2185, 15718, 48988, -78722, -1474424, 0, ...]; ...

%e in which the main diagonal is all zeros after the initial term, illustrating that: 0 = [x^(n-1)] (x*A(x))' / (1 + x*A(x)^n)^n for n>1.

%e RELATED TABLE.

%e The table of coefficients of x^k in (x*A(x))' / A(x)^n begins:

%e n=1: [1, 1, 5, 37, 305, 2516, 19205, 126624, 792081, 10741078, ...];

%e n=2: [1, 0, 2, 20, 186, 1620, 12495, 79302, 474495, 7998088, ...];

%e n=3: [1, -1, 0, 8, 100, 969, 7680, 46440, 261588, 5988094, ...];

%e n=4: [1, -2, -1, 0, 40, 510, 4337, 24608, 126090, 4504280, ...];

%e n=5: [1, -3, -1, -5, 0, 199, 2120, 11037, 46376, 3386810, ...];

%e n=6: [1, -4, 0, -8, -25, 0, 749, 3514, 5499, 2514128, ...];

%e n=7: [1, -5, 2, -10, -39, -116, 0, 290, -9640, 1795607, ...];

%e n=8: [1, -6, 5, -12, -45, -172, -304, 0, -9000, 1165374, ...];

%e n=9: [1, -7, 9, -15, -45, -186, -301, 1594, 0, 577155, ...];

%e n=10: [1, -8, 14, -20, -40, -172, -97, 4278, 11976, 0, ...]; ...

%e in which the main diagonal is all zeros after the initial term, illustrating the identity: 0 = [x^(n-1)] (x*A(x))' / A(x)^n for n>1 holds since A(0) = 1.

%o (PARI) {a(n) = my(A=[1,1]); for(i=1,n, A = concat(A,0); m=#A; A[m] = -Vec( (Ser(A) + x*Ser(A)')/(1 + x*Ser(A)^m)^m/m )[m]);A[n+1]}

%o for(n=0,30, print1(a(n),", "))

%Y Cf. A303559.

%K nonn

%O 0,3

%A _Paul D. Hanna_, Apr 26 2018