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

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A027999 Expansion of Product(1+q^m)^(m(m-1)/2); m=1..inf. 9
1, 0, 1, 3, 6, 13, 24, 49, 91, 181, 334, 632, 1163, 2138, 3880, 7006, 12531, 22279, 39369, 69078, 120597, 209282, 361405, 620829, 1061687, 1807014, 3062642, 5168784, 8688820, 14549659, 24274226, 40353748, 66854518, 110391391, 181695436, 298129605, 487706902 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,4

LINKS

Alois P. Heinz, Table of n, a(n) for n = 0..1000

Vaclav Kotesovec, Graph - the asymptotic ratio

FORMULA

a(n) ~ 7^(1/8) / (2^(47/24) * 15^(1/8) * n^(5/8)) * exp(-2025 * Zeta(3)^3 / (98*Pi^8) - 135*(15/7)^(1/4) * Zeta(3)^2 / (28*Pi^5) * n^(1/4) - 3*sqrt(15/7) * Zeta(3) / (2*Pi^2) * sqrt(n) + 2*(7/15)^(1/4) * Pi/3 * n^(3/4)), where Zeta(3) = A002117. - Vaclav Kotesovec, May 27 2015

MAPLE

b:= proc(n, i) option remember; `if`(n=0, 1, `if`(i<1, 0,

      add(binomial(binomial(i, 2), j)*b(n-i*j, i-1), j=0..n/i)))

    end:

a:= n-> b(n$2):

seq(a(n), n=0..50);  # Alois P. Heinz, Aug 03 2013

MATHEMATICA

b[n_, i_] := b[n, i] = If[n == 0, 1, If[i < 1, 0, Sum[Binomial[Binomial[i, 2], j]*b[n-i*j, i-1], {j, 0, n/i}]]]; a[n_] := b[n, n]; Table[a[n], {n, 0, 50}] (* Jean-Fran├žois Alcover, Oct 13 2014, after Alois P. Heinz *)

CROSSREFS

Cf. A027998, A028377, A258349, A258341, A258344.

Sequence in context: A225199 A000219 A191782 * A005196 A320286 A032287

Adjacent sequences:  A027996 A027997 A027998 * A028000 A028001 A028002

KEYWORD

nonn

AUTHOR

N. J. A. Sloane.

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 June 24 05:21 EDT 2019. Contains 324318 sequences. (Running on oeis4.)