login
G.f. satisfies: [x^n] A(x)^(n+1) = [x^n] A(x)^n for n>1 with A(0)=A'(0)=1.
5

%I #15 Feb 18 2020 11:28:31

%S 1,1,-2,9,-56,425,-3726,36652,-397440,4695489,-59941550,821711605,

%T -12037503384,187689245588,-3104186515976,54295661153700,

%U -1001685184237056,19444296845046033,-396260414466644574,8460628832978195683,-188898511962856879400

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

%H Alois P. Heinz, <a href="/A158883/b158883.txt">Table of n, a(n) for n = 0..448</a>

%F G.f. satisfies: A(x) = 1 + x*(d/dx)(x/A(x)) so that x^2*A'(x) = x*A(x) + A(x)^2 - A(x)^3.

%F a(n) = (-1)^(n-1)*n*A088716(n-1) for n >= 1.

%F G.f.: A(x) = 1/(Sum_{n>=0} (-1)^n*A088716(n)*x^n), where g.f. F(x) of A088716 satisfies: F(x) = 1 + x*F(x)*(d/dx)(x*F(x)).

%F G.f. satisfies: [x^n] A(x)^(n+1) = (n+1)*A158884(n) for n > 1.

%e G.f.: A(x) = 1 + x - 2*x^2 + 9*x^3 - 56*x^4 + 425*x^5 - 3726*x^6 + ...

%e (d/dx) (x/A(x)) = 1 - 2*x + 9*x^2 - 56*x^3 + 425*x^4 - 3726*x^5 + ...

%e 1/A(x) = 1 - x + 3*x^2 - 14*x^3 + 85*x^4 + ... + (-1)^n*A088716(n)*x^n + ...

%e where a(n) = (-1)^(n-1)*n*A088716(n-1) for n >= 1.

%e ...

%e Coefficients of powers of g.f. A(x) begin:

%e A^1: 1,1,-2,9,-56,425,-3726,36652,-397440,4695489,...;

%e A^2: 1,2,(-3),14,-90,702,-6297,63144,-695886,8334822,...;

%e A^3: 1,3,(-3),(16),-108,870,-7997,81774,-915798,11116902,...;

%e A^4: 1,4,-2,(16),(-115),960,-9050,94368,-1073658,13204560,...;

%e A^5: 1,5,0,15,(-115),(996),-9630,102365,-1182690,14730890,...;

%e A^6: 1,6,3,14,-111,(996),(-9870),106890,-1253466,15804548,...;

%e A^7: 1,7,7,14,-105,973,(-9870),(108816),-1294412,16514162,...;

%e A^8: 1,8,12,16,-98,936,-9704,(108816),(-1312227),16931984,...;

%e A^9: 1,9,18,21,-90,891,-9426,107406,(-1312227),(17116900),...;

%e A^10:1,10,25,30,-80,842,-9075,104980,-1298625,(17116900),...; ...

%e where coefficients [x^n] A(x)^(n+1) and [x^n] A(x)^n are

%e enclosed in parenthesis and equal (n+1)*A158884(n) for n > 1:

%e [ -3,16,-115,996,-9870,108816,-1312227,17116900,...];

%e compare to A158884:

%e [1,1,-1,4,-23,166,-1410,13602,-145803,1711690,-21785618,...]

%e and also to the logarithmic derivative of A158884:

%e [1,-3,16,-115,996,-9870,108816,-1312227,17116900,...].

%p b:= proc(n) option remember; `if`(n=0, 1, add(

%p b(j)*b(n-j-1)*(j+1), j=0..n-1))

%p end:

%p a:= n-> `if`(n=0, 1, -(-1)^n*n*b(n-1)):

%p seq(a(n), n=0..25); # _Alois P. Heinz_, Feb 18 2020

%t m = 19; A[_] = 1;

%t Do[A[x_] = 1 + x*D[x/A[x], x] + O[x]^m // Normal, {m}];

%t CoefficientList[A[x], x] (* _Jean-François Alcover_, Feb 18 2020 *)

%o (PARI) {a(n)=local(A=[1,1]);for(i=2,n,A=concat(A,0);A[ #A]=(Vec(Ser(A)^(#A-1))-Vec(Ser(A)^(#A)))[ #A]);A[n+1]}

%o (Maxima)

%o Composita(n,k,F):=if k=1 then F(n) else sum(F(i+1)*Composita(n-i-1,k-1,F),i,0,n-k);

%o array(a, 10);

%o a[1]:1;

%o af(n):=a[n];

%o for n:2 thru 10 do a[n]:n*sum(Composita(n-1, k, af)*(-1)^k , k, 1, n-1);

%o makelist(af(n),n,1,10); /* _Vladimir Kruchinin_, Dec 01 2011 */

%Y Cf. A088716, A158884, variant: A158882.

%K sign

%O 0,3

%A _Paul D. Hanna_, Apr 30 2009