|
Search: id:A094362
|
|
|
| A094362 |
|
Expansion of (eta(q^3)eta(q^13))/(eta(q)eta(q^39)) in powers of q. |
|
+0
1
|
|
| 1, 1, 2, 2, 4, 5, 7, 9, 13, 16, 22, 27, 36, 43, 56, 68, 87, 104, 130, 156, 193, 230, 281, 333, 404, 477, 572, 673, 802, 940, 1113, 1299, 1531, 1780, 2085, 2418, 2820, 3259, 3784, 4362, 5047, 5799, 6685, 7662, 8806, 10066, 11532, 13152, 15026, 17098, 19482
(list; graph; listen)
|
|
|
OFFSET
|
-1,3
|
|
|
COMMENT
|
Euler transform of period 39 sequence [1,1,0,1,1,0,1,1,0,1,1,0,0,1,0,1,1,0,1,1,0,1,1,0,1,0,0,1,1,0,1,1,0,1,1,0,1,1,0,...].
G.f. A(x) satisfies 0=f(A(x),A(x^2))=f(1/A(x),1/A(x^2)) where f(u,v)=u^3+v^3+2uv(u+v)-u^2v^2-uv.
|
|
FORMULA
|
G.f.: x^-1 Product_{k>0} (1-x^(3k))(1-x^(13k))/((1-x^k)(1-x^(39k))).
|
|
PROGRAM
|
(PARI) a(n)=local(A); if(n<-1, 0, n++; A=x*O(x^n); polcoeff(eta(x^3+A)*eta(x^13+A)/eta(x+A)/eta(x^39+A), n))
(PARI) a(n)=local(A, u, v); if(n<0, 0, A=1/x; for(k=2, n, u=A+x*O(x^k); v=subst(u, x, x^2); A-=x^k*polcoeff(u^3+v^3+2*u*v*(u+v)-u^2*v^2-u*v, k+2)/2); polcoeff(A, n))
|
|
CROSSREFS
|
Sequence in context: A027592 A007209 A058661 this_sequence A000726 A128663 A135833
Adjacent sequences: A094359 A094360 A094361 this_sequence A094363 A094364 A094365
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
Michael Somos, May 03 2004
|
|
|
Search completed in 0.002 seconds
|
