|
Search: id:A129390
|
|
|
| A129390 |
|
Expansion of phi(q)* phi(-q^5)/ (chi(-q^2)* chi(-q^10)) in powers of q. |
|
+0 2
|
|
| 1, 2, 1, 2, 3, 0, 0, 2, 0, 0, 4, 2, 1, 4, 2, 0, 0, 2, 0, 0, 2, 2, 3, 2, 3, 0, 0, 0, 0, 0, 2, 6, 0, 2, 4, 0, 0, 2, 0, 0, 5, 2, 0, 4, 2, 0, 0, 0, 0, 0, 2, 2, 4, 2, 2, 0, 0, 2, 0, 0, 1, 4, 1, 2, 4, 0, 0, 4, 0, 0, 4, 0, 2, 6, 2, 0, 0, 0, 0, 0, 4, 2, 0, 2, 1, 0, 0, 2, 0, 0, 2, 4, 0, 0, 8, 0, 0, 0, 0, 0, 4, 4, 2, 6, 0
(list; graph; listen)
|
|
|
OFFSET
|
0,2
|
|
|
FORMULA
|
Expansion of q^(-1/2)* eta(q^2)^4* eta(q^5)^2* eta(q^20)/ (eta(q)^2* eta(q^4)* eta(q^10)^2) in powers of q.
Euler transform of period 20 sequence [ 2, -2, 2, -1, 0, -2, 2, -1, 2, -2, 2, -1, 2, -2, 0, -1, 2, -2, 2, -2, ...].
a(n)= b(2*n+1) where b(n) is multiplicative with b(2^e) = 0, b(5^e) = 1, b(p^e) = e+1 if p == 1, 3, 7, 9 (mod 20), b(p^e) = (1+(-1)^e)/2 if p == 11, 13, 17, 19 (mod 20).
G.f.: Sum_{k>0} a(k)* x^(2*k-1) = Sum_{k>0} f(x^(2*k-1)) where f(x)= x* (1+x^2)* (1+x^6)/ (1+x^10).
|
|
EXAMPLE
|
q + 2*q^3 + q^5 + 2*q^7 + 3*q^9 + 2*q^15 + 4*q^21 + 2*q^23 + q^25 + ...
|
|
PROGRAM
|
(PARI) {a(n)= if(n<0, 0, n= 2*n+1; sumdiv(n, d, kronecker(-20, d)))}
(PARI) {a(n)= local(A, p, e); if(n<0, 0, n= 2*n+1; A= factor(n); prod(k=1, matsize(A)[1], if(p= A[k, 1], e= A[k, 2]; if(p==5, 1, if(p%20 <10, e+1, !(e%2))))))}
(PARI) {a(n)= local(A); if(n<0, 0, A= x*O(x^n); polcoeff( eta(x^2+A)^4* eta(x^5+A)^2* eta(x^20+A)/ (eta(x+A)^2* eta(x^4+A)* eta(x^10+A)^2), n))}
|
|
CROSSREFS
|
A035710(2*n+1)= a(n). A129391(n)= (-1)^n* a(n).
Sequence in context: A138262 A127510 A158810 this_sequence A129391 A123590 A092872
Adjacent sequences: A129387 A129388 A129389 this_sequence A129391 A129392 A129393
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
Michael Somos, Apr 13 2007
|
|
|
Search completed in 0.002 seconds
|