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!)
A133316 G.f.: 1/(1-x) = Sum_{n>=0} a(n)*x^n/(1+x)^(n^2). 7
1, 1, 2, 8, 54, 544, 7508, 133704, 2943194, 77589536, 2391477804, 84582890704, 3382005372970, 151034046369696, 7458091839548356, 403808650013237224, 23801728042233670770, 1517930142778063770304, 104179592763803229618620 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,3

LINKS

Table of n, a(n) for n=0..18.

FORMULA

a(n) = 1 - Sum_{j=0..n-1} a(j) * (-1)^(n-j) * C(j^2 + n-j-1, n-j) for n>0, with a(0)=1.

EXAMPLE

1/(1-x) = 1 + x/(1+x) + 2*x^2/(1+x)^4 + 8*x^3/(1+x)^9 + 54*x^4/(1+x)^16 +...

PROG

(PARI) {a(n)=if(n==0, 1, polcoeff(-(1-x)*sum(m=0, n-1, a(m)*x^m/(1+x +x*O(x^n))^(m^2)), n))}

(PARI) {a(n)=if(n==0, 1, 1 - sum(j=0, n-1, a(j)*(-1)^(n-j)*binomial(j^2+n-j-1, n-j)))}

CROSSREFS

Cf. A141761, A177447, A215241, A215242, A215243.

Sequence in context: A073564 A199576 A005155 * A234301 A005440 A183282

Adjacent sequences:  A133313 A133314 A133315 * A133317 A133318 A133319

KEYWORD

nonn

AUTHOR

Paul D. Hanna, Dec 24 2010

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 April 16 03:56 EDT 2021. Contains 343030 sequences. (Running on oeis4.)