A208546 - OEIS

%I #17 Mar 12 2021 22:24:46

%S 1,1,0,-1,0,0,0,1,1,0,0,0,0,0,1,0,0,0,0,0,-1,0,0,0,0,0,0,0,0,-1,0,1,0,

%T 0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,-1,0,0,0,0,0,0,0,0,0,0,0,0,0,

%U -1,0,0,-1,0,0,0,0,0,0,0,0,0,0,0,0,0,0

%N Expansion of f(-x^29, x^31) + x * f(-x^19, x^41) - x^3 * f(-x^11, x^49) + x^7 * f(x, -x^59) in powers of x where f() is Ramanujan's two-variable theta function.

%H G. C. Greubel, <a href="/A208546/b208546.txt">Table of n, a(n) for n = 0..2500</a>

%H Michael Somos, <a href="/A010815/a010815.txt">Introduction to Ramanujan theta functions</a>

%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/RamanujanThetaFunctions.html">Ramanujan Theta Functions</a>

%F |a(n)| is the characteristic function of A093722.

%F The exponents in the q-series q * A(q^120) are the squares of the numbers in A057538.

%F Euler transform of a period 80 sequence.

%F G.f.: Sum_{k} (-1)^floor(k/4) * x^(3*k * (5*k + 1)/2) * (x^(4*k + 1) + x^(-16*k + 7)).

%e G.f. = 1 + x - x^3 + x^7 + x^8 + x^14 - x^20 - x^29 + x^31 + x^42 - x^52 - x^66 + ...

%e G.f. = q + q^121 - q^361 + q^841 + q^961 + q^1681 - q^2401 - q^3481 + q^3721 + ...

%t f[x_, y_]:= QPochhammer[-x, x*y]*QPochhammer[-y, x*y]*QPochhammer[x*y, x*y]; CoefficientList[Series[f[-x^29, x^31] + x*f[-x^19, x^41] - x^3*f[-x^11, x^49] + x^7*f[x, -x^59], {x, 0, 50}], x] (* _G. C. Greubel_, Aug 11 2018 *)

%o (PARI) {a(n) = local(m); if( issquare( 120*n + 1, &m), (-1)^(m \ 40 + m \ 12))}

%Y Cf. A057538, A093722.

%K sign

%O 0,1

%A _Michael Somos_, Feb 28 2012