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!)
A264905 Expansion of Product_{k>=1} (1 + x^k + x^(3*k)). 7
1, 1, 1, 3, 2, 4, 6, 7, 8, 13, 16, 18, 26, 29, 38, 49, 58, 68, 90, 101, 125, 156, 181, 214, 263, 304, 358, 435, 505, 589, 701, 812, 939, 1115, 1275, 1485, 1736, 1991, 2286, 2667, 3038, 3489, 4028, 4588, 5240, 6036, 6833, 7787, 8904, 10078, 11429, 13020, 14698 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,4

LINKS

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

FORMULA

a(n) ~ c^(1/4) * exp(2*sqrt(c*n)) / (2*sqrt(3*Pi)*n^(3/4)), where c = Integral_{0..infinity} log(1 + exp(-x) + exp(-3*x)) dx = 0.9953865985263189816963357718655148864441174218433250148867... . - Vaclav Kotesovec, Jan 05 2016

MATHEMATICA

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

nmax = 100; p = ConstantArray[0, nmax + 1]; p[[1]] = 1; p[[2]] = 1; p[[4]] = 1; Do[Do[p[[j+1]] = p[[j+1]] + p[[j - k + 1]] + If[j < 3*k, 0, p[[j - 3*k + 1]]], {j, nmax, k, -1}]; , {k, 2, nmax}]; p (* Vaclav Kotesovec, May 10 2018 *)

CROSSREFS

Cf. A000009, A000726, A001935, A100405, A266647, A266648, A266649, A266650, A266686, A275820.

Sequence in context: A239516 A021312 A099258 * A254052 A105746 A166309

Adjacent sequences:  A264902 A264903 A264904 * A264906 A264907 A264908

KEYWORD

nonn

AUTHOR

Vaclav Kotesovec, Nov 28 2015

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 5 06:13 EDT 2020. Contains 333238 sequences. (Running on oeis4.)