%I #95 Sep 25 2024 09:55:12
%S 1,1,0,0,0,1,0,0,1,0,0,0,0,0,0,0,1,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,1,
%T 0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,1,0,0,
%U 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,1,0
%N a(n) = 0 unless n = 3j^2 + 2j or 3j^2 + 4j + 1 for some j >= 0, in which case a(n) = 1.
%C Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).
%C Also characteristic function of generalized octagonal numbers A001082. - _Omar E. Pol_, Jul 13 2012
%C Number 12 of the 14 primitive eta-products which are holomorphic modular forms of weight 1/2 listed by D. Zagier on page 30 of "The 1-2-3 of Modular Forms". - _Michael Somos_, May 04 2016
%H Seiichi Manyama, <a href="/A089801/b089801.txt">Table of n, a(n) for n = 0..1000</a>
%H George E. Andrews, <a href="https://georgeandrews1.github.io/pdf/320.pdf">The Bhargava-Adiga Summation and Partitions</a>, 2016. See Th. 2.
%H S. Cooper and M. D. Hirschhorn, <a href="http://dx.doi.org/10.1016/S0012-365X(03)00079-7">Results of Hurwitz type for three squares.</a> Discrete Math. 274 (2004), no. 1-3, 9-24. See X(q).
%H James A. Sellers, <a href="https://arxiv.org/abs/2409.12321">Elementary Proofs of Two Congruences for Partitions with Odd Parts Repeated at Most Twice</a>, arXiv:2409.12321 [math.NT], 2024. See p. 4.
%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>.
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/JacobiThetaFunctions.html">Jacobi Theta Functions</a>.
%H I. J. Zucker, <a href="http://dx.doi.org/10.1088/0305-4470/23/2/009">Further Relations Amongst Infinite Series and Products. II. The Evaluation of Three-Dimensional Lattice Sums</a>, J. Phys. A: Math. Gen. 23, 117-132, 1990.
%H <a href="/index/Ch#char_fns">Index entries for characteristic functions</a>.
%F G.f.: Sum_{n=-oo..oo} q^(3n^2+2n).
%F Expansion of Jacobi theta function (theta_3(q^(1/3)) - theta_3(q^3))/(2 q^(1/3)) in powers of q.
%F Euler transform of period 12 sequence [1, -1, 0, 0, 1, -1, 1, 0, 0, -1, 1, -1, ...]. - _Michael Somos_, Apr 12 2005
%F a(n) = b(3*n + 1) where b() is multiplicative with b(3^e) = 0^e, b(p^e) = (1 + (-1)^e) / 2 if p != 3. - _Michael Somos_, Jun 06 2005; b=A033684. - _R. J. Mathar_, Oct 07 2011
%F Expansion of q^(-1/3) * eta(q^2)^2 * eta(q^3) * eta(q^12) / (eta(q) * eta(q^4) * eta(q^6)) in powers of q. - _Michael Somos_, Apr 12 2005
%F Expansion of chi(x) * psi(-x^3) in powers of x where chi(), psi() are Ramanujan theta functions. - _Michael Somos_, Apr 19 2007
%F Expansion of f(x, x^5) in powers of x where f(, ) is Ramanujan's general theta function. - _Michael Somos_, Jun 29 2012
%F G.f. is a period 1 Fourier series which satisfies f(-1 / (36 t)) = 2^(1/2) (t/i)^(1/2) g(t) where q = exp(2 Pi i t) and g() is the g.f. for A089807.
%F a(8*n + 4) = a(4*n + 2) = a(4*n + 3) = 0, a(4*n + 1) = a(n), a(8*n) = A080995(n). - _Michael Somos_, Jan 28 2011
%F a(n) = (-1)^n * A089802(n).
%F For n > 0, a(n) = b(n)-b(n-1) + c(n)-c(n-1), where b(n) = floor(sqrt(n/3+1/9)+2/3) and c(n) = floor(sqrt(n/3+1/9)+4/3). - _Mikael Aaltonen_, Jan 22 2015
%F a(n) = A033684(3*n + 1). - _Michael Somos_, Jan 10 2017
%e G.f. = 1 + x + x^5 + x^8 + x^16 + x^21 + x^33 + x^40 + x^56 + x^65 + x^85 + ...
%e G.f. = q + q^4 + q^16 + q^25 + q^49 + q^64 + q^100 + q^121 + q^169 + q^196 + ...
%p A089801 := proc(n)
%p A033684(3*n+1) ;
%p end proc: # _R. J. Mathar_, Oct 07 2011
%p M:=33;
%p S:=f->series(f,q,500);
%p L:=f->seriestolist(f);
%p X:=add(q^(3*n^2+2*n),n=-M..M);
%p S(%);
%p L(%); # _N. J. A. Sloane_, Jan 31 2012
%p eps:=Array(0..120,0);
%p for j from 0 to 120 do
%p if 3*j^2+2*j <= 120 then eps[3*j^2+2*j] := 1; fi;
%p if 3*j^2+4*j+1 <= 120 then eps[3*j^2+4*j+1] := 1; fi;
%p end do; # _N. J. A. Sloane_, Aug 12 2017
%t a[ n_] := SeriesCoefficient[ (1/2) x^(-1/3) (EllipticTheta[ 3, 0, x^(1/3)] - EllipticTheta[ 3, 0, x^3]), {x, 0, n}]; (* _Michael Somos_, Jun 29 2012 *)
%t a[ n_] := SeriesCoefficient[ 2^(-1/2) x^(-3/8) QPochhammer[ -x, x^2] EllipticTheta[ 2, Pi/4, x^(3/2)], {x, 0, n}]; (* _Michael Somos_, Jun 29 2012 *)
%o (PARI) {a(n) = issquare(3*n + 1)}; /* _Michael Somos_, Apr 12 2005 */
%o (Magma) Basis( ModularForms( Gamma0(36), 1/2), 87) [2]; /* _Michael Somos_, Jul 02 2014 */
%o (Python)
%o from sympy.ntheory.primetest import is_square
%o def A089801(n): return int(is_square(3*n+1)) # _Chai Wah Wu_, May 17 2023
%Y Cf. A001082, A033684, A080995, A089802, A089807, A280739.
%K nonn
%O 0,1
%A _Eric W. Weisstein_, Nov 12 2003
%E Edited with simpler definition by _N. J. A. Sloane_, Jan 31 2012
%E Further edited by _N. J. A. Sloane_, Aug 12 2017