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

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A280664 G.f.: Product_{k>=1, j>=1} (1 + x^(j*k^4)). 4
1, 1, 1, 2, 2, 3, 4, 5, 6, 8, 10, 12, 15, 18, 22, 27, 33, 39, 47, 56, 66, 79, 93, 109, 128, 150, 175, 204, 237, 274, 318, 367, 423, 487, 559, 641, 734, 839, 957, 1091, 1241, 1410, 1601, 1814, 2053, 2322, 2622, 2957, 3334, 3752, 4218, 4740, 5318, 5962, 6679 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,4

COMMENTS

In general, if m>=3 and g.f. = Product_{k>=1, j>=1} (1+x^(j*k^m)), then a(n, m) ~ exp(Pi*sqrt(Zeta(m)*n/3) + (2^(1/m)-1) * Pi^(-1/m) * Gamma(1+1/m) * Zeta(1+1/m) * Zeta(1/m) * (3*n/Zeta(m))^(1/(2*m))) * Zeta(m)^(1/4) / (2^(5/4) * 3^(1/4) * n^(3/4)).

LINKS

Vaclav Kotesovec, Table of n, a(n) for n = 0..10000

FORMULA

a(n) ~ exp(Pi^3 * sqrt(n/30)/3 + 2^(-15/8) * 3^(3/8) * 5^(1/8) * (2^(1/4)-1) * Pi^(-3/4) * Gamma(1/4) * Zeta(5/4) * Zeta(1/4) * n^(1/8)) * Pi / (2^(3/2) * 3^(3/4) * 5^(1/4) * n^(3/4)).

MATHEMATICA

nmax = 100; CoefficientList[Series[Product[1+x^(j*k^4), {k, 1, Floor[nmax^(1/4)]+1}, {j, 1, Floor[nmax/k^4]+1}], {x, 0, nmax}], x]

CROSSREFS

Cf. A107742 (m=1), A280451 (m=2), A280663 (m=3).

Cf. A280662.

Sequence in context: A081360 A117409 A092833 * A100926 A258875 A179241

Adjacent sequences:  A280661 A280662 A280663 * A280665 A280666 A280667

KEYWORD

nonn

AUTHOR

Vaclav Kotesovec, Jan 07 2017

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 November 22 13:47 EST 2019. Contains 329393 sequences. (Running on oeis4.)