OFFSET
1,2
LINKS
Paul D. Hanna, Table of n, a(n) for n = 1..2050
EXAMPLE
L.g.f.: A(x) = x + 3*x^3/3 + 5*x^5/5 + 7*x^7/7 + 12*x^9/9 + 11*x^11/11 + 26*x^13/13 + 15*x^15/15 + 51*x^17/17 + 19*x^19/19 + 91*x^21/21 + 23*x^23/23 + 155*x^25/25 + 27*x^27/27 + 232*x^29/29 + 62*x^31/31 + 341*x^33/33 + 35*x^35/35 + 592*x^37/37 + 39*x^39/39 + 656*x^41/41 + 344*x^43/43 + 870*x^45/45 + 47*x^47/47 + 1820*x^49/49 + 51*x^51/51 + 1431*x^53/53 + 1441*x^55/55 + 1843*x^57/57 + 59*x^59/59 + 4758*x^61/61 + 63*x^63/63 + 2925*x^65/65 +...
which may be written as
A(x) = x/(1 - x^2) + x^9/(3*(1 - x^4)^3) + x^25/(5*(1 - x^6)^5) + x^49/(7*(1 - x^8)^7) + x^81/(9*(1 - x^10)^9) + x^121/(11*(1 - x^12)^11) + x^169/(13*(1 - x^14)^13) +...+ x^((2*n-1)^2) / ((2*n-1)*(1 - x^(2*n))^(2*n-1)) +...
The coefficient of x^(2^n+1)/(2^n+1) in A(x) for n>=1 begins:
[3, 5, 12, 51, 341, 2925, 169033, 33445209, 21619038033, ..., A293599(n), ...].
MATHEMATICA
nmax = 80; Table[(CoefficientList[Series[Sum[x^((2*k - 1)^2)/((2*k - 1)*(1 - x^(2*k))^(2*k - 1)), {k, 1, 2*nmax + 1}], {x, 0, 2*nmax + 1}], x] * Range[0, 2*nmax + 1])[[2*n]], {n, 1, nmax}] (* Vaclav Kotesovec, Oct 15 2017 *)
PROG
(PARI) {a(n) = my(A, Ox = O(x^(2*n+1)));
A = sum(m=1, sqrtint(n+1), x^((2*m-1)^2) / ( (2*m-1) * (1 - x^(2*m) +Ox)^(2*m-1) ) );
(2*n-1)*polcoeff(A, 2*n-1)}
for(n=1, 80, print1(a(n), ", "))
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Oct 12 2017
STATUS
approved