OFFSET
0,2
FORMULA
a(n)=b(12n+1) where b(n) is multiplicative and b(2^e) = b(3^e) = 0^e, b(p^e) = (1+(-1)^e)/2*p^e if p == 7,11 (mod 12), b(p^e) = (-1)^(e/2)(1+(-1)^e)/2*p^e if p == 5 (mod 12), b(p^e) = (e+1)*((-1)^y*p)^e where p == 1 (mod 12) and p = x^2+9y^2.
G.f.: (1 -24*Sum_{k>0} x^k/(1-x^k)^2)(Product_{k>0} 1-x^k)^2.
PROG
(PARI) {a(n)=local(A); if(n<0, 0, A=x*O(x^n); polcoeff( eta(x+A)^2*sum(k=1, n, -24*sigma(k)*x^k, 1+A), n))}
(PARI) {a(n)= local(A, p, e, y); if(n<0, 0, n=12*n+1; A=factor(n); prod(k=1, matsize(A)[1], if(p=A[k, 1], e=A[k, 2]; if(p<5, 0, if(p%12>1, if(e%2, 0, ((-1)^(p%12==5)*p^2)^(e/2)), for(i=1, sqrtint(p\9), if(issquare(p-9*i^2), y=i; break)); (e+1)*((-1)^y*p)^e)))))}
CROSSREFS
KEYWORD
sign
AUTHOR
Michael Somos, Aug 22 2006
STATUS
approved