This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
%I A028586 #32 Feb 16 2025 08:32:35
%S A028586 1,0,2,4,0,0,0,4,2,0,2,0,4,0,0,4,0,0,6,0,0,0,0,4,0,0,0,8,4,0,0,0,2,0,
%T A028586 0,4,0,0,0,0,2,0,8,4,0,0,0,4,4,0,2,0,0,0,0,0,0,0,4,0,4,0,0,12,0,0,0,4,
%U A028586 0,0,0,0,6,0,0,4,0,0,0,0,0,0,4,4,0,0,0,8,0,0,6,0,4,0,0,0,0,0,6,0,0,0,0,4
%N A028586 Theta series of lattice with Gram matrix [2 1; 1 3].
%C A028586 Ramanujan theta functions: f(q) := Prod_{k>=1} (1-(-q)^k) (see A121373), phi(q) := theta_3(q) := Sum_{k=-oo..oo} q^(k^2) (A000122), psi(q) := Sum_{k=0..oo} q^(k*(k+1)/2) (A010054), chi(q) := Prod_{k>=0} (1+q^(2k+1)) (A000700).
%C A028586 The number of integer solutions (x, y) to 2*x^2 + 2*x*y + 3*y^2 = n, discriminant -20. - _Ray Chandler_, Jul 12 2014
%H A028586 John Cannon, <a href="/A028586/b028586.txt">Table of n, a(n) for n = 0..10000</a>
%H A028586 A. Berkovich and H. Yesilyurt, <a href="https://arxiv.org/abs/math/0611300">Ramanujan's identities and representation of integers by certain binary and quaternary quadratic forms</a>, arXiv:math/0611300 [math.NT], 2006-2007, page 8 equation (3.18).
%H A028586 N. J. A. Sloane et al., <a href="https://oeis.org/wiki/Binary_Quadratic_Forms_and_OEIS">Binary Quadratic Forms and OEIS</a> (Index to related sequences, programs, references)
%H A028586 Michael Somos, <a href="/A010815/a010815.txt">Introduction to Ramanujan theta functions</a>
%H A028586 Eric Weisstein's World of Mathematics, <a href="https://mathworld.wolfram.com/RamanujanThetaFunctions.html">Ramanujan Theta Functions</a>
%F A028586 G.f.: Sum_{n,m} x^(2*n^2 + 2*m*n + 3*m^2). - _Michael Somos_, Jan 31 2011
%F A028586 Expansion of (theta_3(z)*theta_3(5z)+theta_2(z)*theta_2(5z)).
%F A028586 Expansion of phi(q^2) * phi(q^10) + 4 * q^3 * psi(q^4) * psi(q^20) in powers of q where phi(q),psi(q) are Ramanujan theta functions. - _Michael Somos_, Aug 13 2006
%F A028586 If p is prime then a(p) is nonzero iff p is in A106865.
%F A028586 0=a(n)a(2n) and 2*A035170(n)=a(n)+a(2n) if n>0. - _Michael Somos_, Oct 21 2006
%e A028586 1 + 2*q^2 + 4*q^3 + 4*q^7 + 2*q^8 + 2*q^10 + 4*q^12 + 4*q^15 + 6*q^18 + 4*q^23 + 8*q^27 + 4*q^28 + 2*q^32 + 4*q^35 + 2*q^40 + 8*q^42 + 4*q^43 + 4*q^47 + ...
%t A028586 terms = 104; phi[q_] := EllipticTheta[3, 0, q]; chi[q_] := ((1 - InverseEllipticNomeQ[q])*InverseEllipticNomeQ[q]/(16*q))^(-1/24); psi[q_] := (1/2)*q^(-1/8)*EllipticTheta[2, 0, q^(1/2)]; s = phi[q^2]*phi[q^10] + 4*q^3*psi[q^4]*psi[q^20] + O[q]^(terms+1); CoefficientList[s, q] (* _Jean-François Alcover_, Jul 04 2017, after _Michael Somos_ *)
%t A028586 r[n_]:=Reduce[{x, y}.{{2, 1}, {1, 3}}.{x, y}==n, {x, y}, Integers]; Table[rn=r[n]; Which[rn===False, 0, Head[rn]===Or, Length[rn], Head[rn]===And, 1], {n, 0, 105}] (* _Vincenzo Librandi_, Feb 23 2020 *)
%o A028586 (PARI) {a(n) = if( n<1, n==0, qfrep([2, 1; 1, 3], n)[n] * 2)} /* _Michael Somos_, Aug 13 2006 */
%K A028586 nonn
%O A028586 0,3
%A A028586 _N. J. A. Sloane_