login
The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A192415 G.f. satisfies: A(x) = (1 + x*A(x))*(1 + x^3*A(x)^2). 17
1, 1, 1, 2, 5, 11, 23, 51, 120, 286, 681, 1636, 3985, 9803, 24257, 60338, 150931, 379501, 958360, 2429294, 6179380, 15769380, 40361087, 103579221, 266471500, 687098810, 1775440421, 4596689688, 11922774513, 30977768907, 80615085087, 210103228155, 548352756656, 1433053608502 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,4

LINKS

Vincenzo Librandi, Table of n, a(n) for n = 0..200

FORMULA

G.f. satisfies: A(x) = exp( Sum_{n>=1} x^n/n * Sum_{k=0..n} C(n,k)^2 * x^(2*k) * A(x)^k ).

D-finite with recurrence: 4*(n+1)*(n+2)*(217*n^3 - 1239*n^2 + 1838*n - 336)*a(n) = 6*(n+1)*(434*n^4 - 2261*n^3 + 2339*n^2 + 1792*n - 1344)*a(n-1) - (n-1)*(2821*n^4 - 13286*n^3 + 7829*n^2 + 18464*n - 4032)*a(n-2) + 6*(868*n^5 - 6258*n^4 + 13981*n^3 - 7438*n^2 - 7769*n + 6136)*a(n-3) + 2*(n-3)*(2*n - 5)*(434*n^3 - 1393*n^2 + 211*n + 1048)*a(n-4) + 2*(n-4)*(2*n - 7)*(217*n^3 - 588*n^2 + 11*n + 480)*a(n-5). - Vaclav Kotesovec, Sep 19 2013

a(n) ~ c*d^n/(sqrt(Pi)*n^(3/2)), where d = 2.730683387097269698... is the root of the equation -4 - 8*d - 24*d^2 + 13*d^3 - 12*d^4 + 4*d^5 = 0 and c = 2.078548317061344694159945441842754... is the root of the equation -1 - 67*c^2 - 19811*c^4 + 36463*c^6 - 41664*c^8 + 7936*c^10 = 0. - Vaclav Kotesovec, Sep 19 2013, updated Nov 28 2016

a(n) = Sum_{k=0..n/2+1} C(n-k+2,k-1)*C(n-k+2,2*k-1)/(n-k+2). - Vladimir Kruchinin, Feb 12 2019

EXAMPLE

G.f.: A(x) = 1 + x + x^2 + 2*x^3 + 5*x^4 + 11*x^5 + 23*x^6 + 51*x^7 +...

Related expansions:

A(x)^2 = 1 + 2*x + 3*x^2 + 6*x^3 + 15*x^4 + 36*x^5 + 82*x^6 + 190*x^7 +...

A(x)^3 = 1 + 3*x + 6*x^2 + 13*x^3 + 33*x^4 + 84*x^5 + 205*x^6 + 498*x^7 +...

where A(x) = 1 + x*A(x) + x^3*A(x)^2 + x^4*A(x)^3.

The logarithm of the g.f. equals the series:

log(A(x)) = (1 + x^2*A(x))*x + (1 + 2^2*x^2*A(x) + x^4*A(x)^2)*x^2/2 +

(1 + 3^2*x^2*A(x) + 3^2*x^4*A(x)^2 + x^6*A(x)^3)*x^3/3 +

(1 + 4^2*x^2*A(x) + 6^2*x^4*A(x)^2 + 4^2*x^6*A(x)^3 + x^8*A(x)^4)*x^4/4 +

(1 + 5^2*x^2*A(x) + 10^2*x^4*A(x)^2 + 10^2*x^6*A(x)^3 + 5^2*x^8*A(x)^4 + x^10*A(x)^5)*x^5/5 +...

Explicitly,

log(A(x)) = x + x^2/2 + 4*x^3/3 + 13*x^4/4 + 31*x^5/5 + 70*x^6/6 + 176*x^7/7 + 469*x^8/8 + 1228*x^9/9 + 3161*x^10/10 +...

MATHEMATICA

nmax=20; aa=ConstantArray[0, nmax]; aa[[1]]=1; Do[AGF=1+Sum[aa[[n]]*x^n, {n, 1, j-1}]+koef*x^j; sol=Solve[Coefficient[(1 + x*AGF)*(1 + x^3*AGF^2) - AGF, x, j]==0, koef][[1]]; aa[[j]]=koef/.sol[[1]], {j, 2, nmax}]; Flatten[{1, aa}] (* Vaclav Kotesovec, Sep 19 2013 *)

PROG

(PARI) {a(n)=local(A=1+x); for(i=1, n, A=(1 + x*A)*(1 + x^3*(A+x*O(x^n))^2)); polcoeff(A, n)}

(PARI) {a(n)=local(A=1+x); for(i=1, n, A=exp(sum(m=1, n, sum(j=0, m, binomial(m, j)^2*x^(2*j)*(A+x*O(x^n))^j)*x^m/m))); polcoeff(A, n, x)}

(PARI)

x='x; y='y; Fxy = (1+x*y) * (1 + x^3*y^2) - y;

seq(N) = {

  my(y0 = 1 + O('x^N), y1=0);

  for (k = 1, N,

    y1 = y0 - subst(Fxy, y, y0)/subst(deriv(Fxy, y), y, y0);

    if (y1 == y0, break()); y0 = y1);

  Vec(y0);

};

seq(34)  \\ Gheorghe Coserea, Nov 30 2016

CROSSREFS

Cf. A198951, A198953, A198957, A007863, A036765, A104545.

Sequence in context: A083380 A018112 A274489 * A201779 A067149 A191240

Adjacent sequences:  A192412 A192413 A192414 * A192416 A192417 A192418

KEYWORD

nonn

AUTHOR

Paul D. Hanna, Nov 02 2011

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified July 14 16:21 EDT 2020. Contains 335729 sequences. (Running on oeis4.)