OFFSET
0,2
COMMENTS
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..1000
Michael Somos, Introduction to Ramanujan theta functions, 2019.
Eric Weisstein's World of Mathematics, Ramanujan Theta Functions.
FORMULA
Expansion of phi(-x^3)^4 / (chi(-x) * chi(-x^3))^3 in powers of x where phi(), chi() are Ramanujan theta functions.
Euler transform of period 6 sequence [ 3, 0, -2, 0, 3, -4, ...].
a(n) = b(2*n + 1) where b(n) is multiplicative with b(2^e) = 0^e, b(3^e) = 3^e, b(p^e) = (p^(e+1) - 1) / (p - 1) if p>3.
G.f. is a period 1 Fourier series which satisfies f(-1 / (12 t)) = 2 (t/i)^2 g(t) where q = exp(2 Pi i t) and g() is g.f. for A118271.
Sum_{k=1..n} a(k) ~ c * n^2, where c = Pi^2/9 = 1.0966227... (A100044). - Amiram Eldar, Dec 28 2023
EXAMPLE
1 + 3*x + 6*x^2 + 8*x^3 + 9*x^4 + 12*x^5 + 14*x^6 + 18*x^7 + 18*x^8 + ...
q + 3*q^3 + 6*q^5 + 8*q^7 + 9*q^9 + 12*q^11 + 14*q^13 + 18*q^15 + ...
MATHEMATICA
A185717[n_] := SeriesCoefficient[(QPochhammer[q^3, q^3]/QPochhammer[-q^3, q^3])^4*(1/(QPochhammer[q, q^2]*QPochhammer[q^3, q^6])^3), {q, 0, n}];
Table[A185717[n], {n, 0, 50}] (* G. C. Greubel, Jul 10 2017 *)
PROG
(PARI) {a(n) = if( n<0, 0, n = 2*n + 1; sumdiv( n, d, if( (n/d) % 3, 1, 0) * d))}
(PARI) {a(n) = local(A); if( n<0, 0, A = x * O(x^n); polcoeff( eta(x^2 + A)^3 * eta(x^3 + A)^5 / (eta(x + A)^3 * eta(x^6 + A)), n))}
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Michael Somos, Feb 10 2011
STATUS
approved