|
Search: id:A125095
|
|
|
| A125095 |
|
Expansion of phi(-q)psi(q^4) in powers of q where psi(), phi() are Ramanujan theta functions. |
|
+0 2
|
|
| 1, -2, 0, 0, 3, -2, 0, 0, 2, -2, 0, 0, 1, -4, 0, 0, 4, 0, 0, 0, 2, -2, 0, 0, 1, -4, 0, 0, 4, -2, 0, 0, 0, -2, 0, 0, 2, -2, 0, 0, 5, -2, 0, 0, 2, 0, 0, 0, 2, -6, 0, 0, 0, -2, 0, 0, 2, 0, 0, 0, 3, -4, 0, 0, 4, -2, 0, 0, 2, -2, 0, 0, 0, -2, 0, 0, 6, 0, 0, 0, 0, -2, 0, 0, 1, -6, 0, 0, 4, -2, 0, 0, 0, -4, 0, 0, 2, 0, 0, 0, 4, 0, 0, 0, 4
(list; graph; listen)
|
|
|
OFFSET
|
0,2
|
|
|
FORMULA
|
Expansion of q^(-1/2)(eta(q)^2*eta(q^8)^2)/(eta(q^2)*eta(q^4)) in powers of q.
Given g.f. A(x), then B(x)=x*A(x^2) satisfies 0=f(B(x), B(x^2), B(x^3), B(x^6)) where f(u1, u2, u3, u6)= +u1^2*u6*(u1+3*u3) +2*u2^2*u3*(u2+3*u6) -3*u3^2*u2*(u1+u3) -6*u6^2*u1*(u2+u6).
a(n) = b(2n+1) where b(n) is multiplicative and b(2^e) = 0^e, b(p^e) = (e+1)(-1)^e if p == 1, 3 (mod 8), b(p^e) = (1+(-1)^e)/2 if p == 5, 7 (mod 8).
Euler transform of period 8 sequence [ -2, -1, -2, 0, -2, -1, -2, -2, ...].
a(4n+2)=a(4n+3)=0.
G.f.: (Sum_{k} (-1)^k*x^k^2)(Sum_{k>=0} x^(2k^2+2k)).
|
|
EXAMPLE
|
q - 2*q^3 + 3*q^9 - 2*q^11 + 2*q^17 - 2*q^19 + q^25 - 4*q^27 +...
|
|
PROGRAM
|
(PARI) {a(n)= if(n<0, 0, (-1)^n*sumdiv(2*n+1, d, (-1)^(d%8>3)))}
(PARI) {a(n)= if(n<0, 0, n=2*n+1; qfrep([1, 0; 0, 8], n)[n] -qfrep([3, 1; 1, 3], n)[n])}
(PARI) {a(n)= local(A); if(n<0, 0, A=x*O(x^n); polcoeff( eta(x+A)^2*eta(x^8+A)^2/ eta(x^2+A)/eta(x^4+A), n))}
(PARI) {a(n) = if( n<0, 0, n = 2*n + 1; sumdiv(n, d, kronecker(2, d) * kronecker(-4, n/d)))}
|
|
CROSSREFS
|
A113411(n)*(-1)^n=a(n).
Sequence in context: A077888 A167634 A113411 this_sequence A143161 A142886 A099026
Adjacent sequences: A125092 A125093 A125094 this_sequence A125096 A125097 A125098
|
|
KEYWORD
|
sign
|
|
AUTHOR
|
Michael Somos, Nov 20 2006
|
|
|
Search completed in 0.002 seconds
|