OFFSET
0,2
COMMENTS
(-1)^n*a(n) is the number of inequivalent elements of norm 8*n-1 in Z[sqrt(2)].
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..1000
Daniel Corson, David Favero, Kate Liesinger, Sarah Zubairy, Characters and q-series in Q(sqrt(2)), J. Number Theory, 107 (2004), 392-405.
Jeremy Lovejoy, Overpartitions and real quadratic fields, J. Number Theory, 106 (2004), 178-186.
MATHEMATICA
With[{m=80}, CoefficientList[Series[1+2*Sum[(-q)^n*QPochhammer[q^4, q^4]*QPochhammer[q^(2*n+1), q]/((1+q^(2*n))*QPochhammer[q^(4*n+4), q^4]*QPochhammer[q, q]), {n, 1, m}], {q, 0, m}], q]] (* G. C. Greubel, Dec 04 2018 *)
PROG
(PARI) m=80; my(q='q+O('q^m)); Vec(1 + 2*sum(n=1, m, ((-q)^n/(1+q^(2*n) ))*prod(k=1, n, (1-q^(4*k))/((1-q^(2*k))*(1-q^(2*k-1)))) )) \\ G. C. Greubel, Dec 04 2018
(Magma) m:=80; R<q>:=PowerSeriesRing(Integers(), m); [1] cat Coefficients(R!( 2*(&+[((-q)^n/(1+q^(2*n)))*(&*[(1-q^(4*k))/((1-q^(2*k))*(1-q^(2*k-1))): k in [1..n]]): n in [1..m]]) )); // G. C. Greubel, Dec 04 2018
(Sage)
from sage.combinat.q_analogues import q_pochhammer
prec = 80
R = PowerSeriesRing(ZZ, 'x')
x = R.gen().O(prec)
s = 1+2*sum( (-x)^n*q_pochhammer(n, x^4, x^4)/((1+x^(2*n))* q_pochhammer(2*n, x, x)) for n in (1..prec))
print(s.coefficients()) # G. C. Greubel, Dec 04 2018
CROSSREFS
KEYWORD
sign
AUTHOR
Jeremy Lovejoy, Jun 12 2009
STATUS
approved