|
|
A122135
|
|
Expansion of f(x, -x^4) / phi(-x^2) in powers of x where f(, ) and phi() are Ramanujan theta functions.
|
|
21
|
|
|
1, 1, 2, 2, 3, 4, 6, 7, 10, 12, 16, 20, 26, 31, 40, 48, 60, 72, 89, 106, 130, 154, 186, 220, 264, 310, 370, 433, 512, 598, 704, 818, 958, 1110, 1293, 1494, 1734, 1996, 2308, 2650, 3052, 3496, 4014, 4584, 5248, 5980, 6825, 7760, 8834, 10020, 11380, 12882, 14594
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,3
|
|
COMMENTS
|
Generating function arises naturally in Rodney Baxter's solution of the Hard Hexagon Model according to George Andrews.
Conjecture: Also the number of integer partitions y of n such that y_i > y_{i+1} for all even i. For example, the a(1) = 1 through a(9) = 12 partitions are:
(1) (2) (3) (4) (5) (6) (7) (8) (9)
(11) (21) (22) (32) (33) (43) (44) (54)
(31) (41) (42) (52) (53) (63)
(221) (51) (61) (62) (72)
(321) (331) (71) (81)
(2211) (421) (332) (432)
(3211) (431) (441)
(521) (531)
(3311) (621)
(4211) (3321)
(4311)
(5211)
The even-length case appears to be A122134.
(End)
|
|
REFERENCES
|
G. E. Andrews, q-series, CBMS Regional Conference Series in Mathematics, 66, Amer. Math. Soc. 1986, see p. 8, Eq. (1.5). MR0858826 (88b:11063)
G. E. Andrews, R. Askey and R. Roy, Special Functions, Cambridge University Press, 1999; Exercise 6(d), p. 591.
|
|
LINKS
|
|
|
FORMULA
|
Expansion of f(x^2, x^8) / f(-x, -x^4) in powers of x where f(, ) is Ramanujan's general theta function. - Michael Somos, Nov 12 2016
Expansion of f(-x^3, -x^7) * f(-x^4, -x^16) / ( f(-x) * f(-x^20) ) in powers of x where f(, ) is Ramanujan's general theta function.
Euler transform of period 20 sequence [ 1, 1, 0, 0, 1, 1, 0, 1, 1, 0, 1, 1, 0, 1, 1, 0, 0, 1, 1, 0, ...].
G.f.: Sum_{k>=0} x^(k^2 + k) / ((1 - x) * (1 - x^2) * ... * (1 - x^(2*k+1))).
Let f(n) = 1/Product_{k >= 0} (1-q^(20k+n)). Then g.f. is f(1)*f(2)*f(5)*f(6)*f(8)*f(9)*f(11)*f(12)*f(14)*f(15)*f(18)*f(19); - N. J. A. Sloane, Mar 19 2012.
a(n) ~ (3 + sqrt(5))^(1/4) * exp(Pi*sqrt(2*n/5)) / (4*sqrt(5)*n^(3/4)). - Vaclav Kotesovec, Nov 12 2016
|
|
EXAMPLE
|
G.f. = 1 + x + 2*x^2 + 2*x^3 + 3*x^4 + 4*x^5 + 6*x^6 + 7*x^7 + 10*x^8 + ...
G.f. = q^9 + q^49 + 2*q^89 + 2*q^129 + 3*q^169 + 4*q^209 + 6*q^249 + ...
|
|
MAPLE
|
f:=n->1/mul(1-q^(20*k+n), k=0..20);
f(1)*f(2)*f(5)*f(6)*f(8)*f(9)*f(11)*f(12)*f(14)*f(15)*f(18)*f(19);
|
|
MATHEMATICA
|
a[ n_] := SeriesCoefficient[ QPochhammer[ -x, -x^5] QPochhammer[ x^4, -x^5] QPochhammer[-x^5] / EllipticTheta[ 4, 0, x^2], {x, 0, n}]; (* Michael Somos, Nov 12 2016 *)
nmax = 50; CoefficientList[Series[Product[1/((1 - x^(20*k+1))*(1 - x^(20*k+2))*(1 - x^(20*k+5))*(1 - x^(20*k+6))*(1 - x^(20*k+8))*(1 - x^(20*k+9))*(1 - x^(20*k+11))*(1 - x^(20*k+12))*(1 - x^(20*k+14))*(1 - x^(20*k+15))*(1 - x^(20*k+18))*(1 - x^(20*k+19)) ), {k, 0, nmax}], {x, 0, nmax}], x] (* Vaclav Kotesovec, Nov 12 2016 *)
|
|
PROG
|
(PARI) {a(n) = if( n<0, 0, polcoeff( sum(k=0, (sqrtint(4*n+1) - 1) \2, x^(k^2 + k) / prod(i=1, 2*k+1, 1 - x^i, 1 + x * O(x^(n-k^2-k)))), n))};
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|