OFFSET
0
LINKS
Seiichi Manyama, Table of n, a(n) for n = 0..10000
Michael Somos, Introduction to Ramanujan theta functions
Eric Weisstein's World of Mathematics, Ramanujan Theta Functions
FORMULA
G.f.: f(x, x^6) = Sum_{k in Z} x^((7k^2 + 5k)/2).
G.f.: Product_{k>0} (1 + x^(7*k-6)) * (1 + x^(7*k-1)) * (1 - x^(7*k)). - Michael Somos, Jan 28 2017
a(3*n + 2) = a(5*n + 2) = a(5*n + 3) = 0. - Michael Somos, Jan 28 2017
a(n) = abs(A232714(n)). - Michael Somos, Jan 28 2017
EXAMPLE
G.f. = 1 + x + x^6 + x^9 + x^19 + x^24 + x^39 + x^46 + x^66 + x^75 + x^100 + ...
G.f. = q^25 + q^61 + q^241 + q^349 + q^709 + q^889 + q^1429 + q^1681 + q^2401 + ...
MATHEMATICA
a[ n_] := If[ n < 0, 0, Boole @ IntegerQ @ Sqrt @ n]; (* Michael Somos, Jan 28 2017 *)
a[ n_] := SeriesCoefficient[ QPochhammer[ -x, x^7] QPochhammer[ -x^6, x^7] QPochhammer[ x^7], {x, 0, n}]; (* Michael Somos, Jan 28 2017 *)
a[ n_] := SeriesCoefficient[ Product[ (1 - x^k)^{-1, 1, 0, 0, 0, -1, 1, -1, 0, 0, 0, 1, -1, 1}[[Mod[k, 14, 1]]], {k, n}], {x, 0, n}]; (* Michael Somos, Jan 28 2017 *)
PROG
(PARI) {a(n) = issquare(56*n + 25)}; /* Michael Somos, Jan 28 2017 */
(PARI) {a(n) = if( n<0, 0, polcoeff( prod(k=1, n, (1 - x^k)^[1, -1, 1, 0, 0, 0, -1, 1, -1, 0, 0, 0, 1, -1][k%14 + 1], 1 + x * O(x^n)), n))}; /* Michael Somos, Jan 28 2017 */
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, Jun 12 2016
STATUS
approved