|
| |
|
|
A092833
|
|
Expansion of eta(q^2)eta(q^46)/(eta(q)eta(q^23)) in powers of q.
|
|
0
| |
|
|
0, 1, 1, 1, 2, 2, 3, 4, 5, 6, 8, 10, 12, 15, 18, 22, 27, 32, 38, 46, 54, 64, 76, 89, 105, 123, 143, 167, 194, 225, 260, 301, 346, 398, 458, 524, 600, 686, 782, 891, 1014, 1151, 1306, 1480, 1674, 1892, 2137, 2409, 2713, 3053, 3431, 3852, 4322, 4842, 5421, 6064, 6776
(list; graph; refs; listen; history; internal format)
|
|
|
|
OFFSET
| 0,5
|
|
|
COMMENTS
| Euler transform of period 46 sequence with g.f. x/(1-x^2)+x^23/(1-x^46).
G.f. A(x) satisfies 0=f(A(x),A(x^2)) where f(u,v)=u^2-v-2uv(1+v).
|
|
|
FORMULA
| G.f.: x(Product_{k>0} (1+x^k)(1+x^(23k))).
|
|
|
PROG
| (PARI) {a(n)=local(A, m); if(n<0, 0, A=x+O(x^2); m=1; while(m<=n, m*=2; A=subst(A, x, x^2); A=A+A^2+sqrt(A+(A+A^2)^2)); polcoeff(A, n))}
(PARI) {a(n)=local(A); if(n<1, 0, n--; A=x*O(x^n); polcoeff(eta(x^2+A)*eta(x^46+A)/eta(x+A)/eta(x^23+A), n))}
|
|
|
CROSSREFS
| Sequence in context: A000009 A081360 A117409 * A100926 A179241 A157046
Adjacent sequences: A092830 A092831 A092832 * A092834 A092835 A092836
|
|
|
KEYWORD
| nonn
|
|
|
AUTHOR
| Michael Somos, Mar 06 2004
|
| |
|
|
