OFFSET
0,11
COMMENTS
LINKS
Seiichi Manyama, Table of n, a(n) for n = 0..1000 (first 251 terms from Reinhard Zumkeller)
Michael Somos, Introduction to Ramanujan theta functions
Eric Weisstein's World of Mathematics, Ramanujan Theta Functions
FORMULA
Expansion of f( x, x^4) / f(-x^5, -x^10) in powers of x where f() is the Ramanujan two-variable theta function. - Michael Somos, Mar 23 2013
Expansion of (-x; x^5)_oo (-x^4; x^5)_oo in powers of x where (x; q)_oo is the q-Pochhammer symbol. - Michael Somos, Mar 23 2013
Euler transform of period 10 sequence [ 1, -1, 0, 1, 0, 1, 0, -1, 1, 0, ...]. - Michael Somos, Mar 23 2013
G.f.: Product_{k>0} (1 + x^(5*k - 1)) * (1 + x^(5*k - 4)). - Michael Somos, Mar 23 2013
a(n) ~ exp(sqrt(2*n/15)*Pi) / (2*30^(1/4)*n^(3/4)) * (1 + (Pi/(60*sqrt(30)) - 3*sqrt(15/2)/(8*Pi)) / sqrt(n)). - Vaclav Kotesovec, Jan 18 2017, extended Jan 24 2017
EXAMPLE
a(10) = #{9+1, 6+4} = 2;
a(20) = #{19+1, 16+4, 14+6, 11+9, 9+6+4+1} = 5.
1 + x + x^4 + x^5 + x^6 + x^7 + x^9 + 2*x^10 + 2*x^11 + x^12 + x^13 + 2*x^14 + ...
q + q^61 + q^241 + q^301 + q^361 + q^421 + q^541 + 2*q^601 + 2*q^661 + q^721 + ...
MATHEMATICA
a[ n_] := SeriesCoefficient[ Product[ (1 + x^(5 k - 1)) (1 + x^(5 k - 4)), {k, Ceiling[ n / 5]}], {x, 0, n}] (* Michael Somos, Mar 23 2013 *)
a[ n_] := SeriesCoefficient[ QPochhammer[ -x, x^5] QPochhammer[ -x^4, x^5], {x, 0, n}] (* Michael Somos, Mar 23 2013 *)
PROG
(Haskell)
a203776 = p a047209_list where
p _ 0 = 1
p (k:ks) m = if m < k then 0 else p ks (m - k) + p ks m
(PARI) {a(n) = polcoeff( prod( k=1, ceil(n / 5), (1 + x^(5*k - 1)) * (1 + x^(5*k - 4)), 1 + x * O(x^n)), n)} /* Michael Somos, Mar 23 2013 */
CROSSREFS
KEYWORD
nonn
AUTHOR
Reinhard Zumkeller, Jan 05 2012
STATUS
approved