login
This site is supported by donations to The OEIS Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A212391 a(n) = A212392(n) / n. 4

%I

%S 1,1,3,14,80,516,3608,26729,206808,1655232,13612512,114466491,

%T 980575020,8533242324,75267759072,671721353474,6056517394512,

%U 55104831724236,505422858053560,4669306663437888,43418090784597696,406109012334694211,3818890067546807794

%N a(n) = A212392(n) / n.

%H Paul D. Hanna, <a href="/A212391/b212391.txt">Table of n, a(n) for n = 1..200</a>

%F Given g.f. A(x), then G(x) = d/dx A(x^3)/3 = Sum_{n>=1} n*a(n)*x^(3*n-1) is the g.f. of A212392 and satisfies: G(x) = (x + G(G(x)))^2.

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

%e G.f.: A(x) = x + x^2 + 3*x^3 + 14*x^4 + 80*x^5 + 516*x^6 + 3608*x^7 + 26729*x^8 +...

%e Let G(x) = d/dx A(x^3)/3, then G(x) = (x + G(G(x)))^2, where

%e G(x) = x^2 + 2*x^5 + 9*x^8 + 56*x^11 + 400*x^14 + 3096*x^17 + 25256*x^20 +...

%e G(G(x)) = x^4 + 4*x^7 + 24*x^10 + 168*x^13 + 1284*x^16 + 10384*x^19 +...

%o (PARI) {a(n)=local(G=x^2+x^3);for(i=1,n,G=(x+subst(G,x,G+O(x^(3*n))))^2);polcoeff(G,3*n-1)/n}

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

%Y Cf. A212392.

%K nonn

%O 1,3

%A _Paul D. Hanna_, May 12 2012

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 April 23 02:43 EDT 2019. Contains 322380 sequences. (Running on oeis4.)