OFFSET
0,2
LINKS
Paul D. Hanna, Table of n, a(n) for n = 0..180
FORMULA
a(n) = [x^(n*(n+1))] Sum_{k>=0} x^k * (1 + x^k)^k / (1 + x^(k+1))^(k+1).
a(n) = [x^(n*(n+1))] Sum_{k>=0} (-x)^k * (1 - x^k)^k / (1 - x^(k+1))^(k+1).
a(n) = [x^(n*(n+1))] Sum_{m>=0} x^m * Sum_{k=0..m} binomial(m,k) * (x^m - x^k)^(m-k).
a(n) = [x^(n*(n+1))] Sum_{m>=0} x^m * Sum_{k=0..m} binomial(m,k) * (-1)^k * (x^m + x^k)^(m-k).
a(n) = [x^(n*(n+1))] Sum_{m>=0} x^m * Sum_{k=0..m} binomial(m,k) * (-1)^k * Sum_{j=0..m-k} binomial(m-k,j) * x^((m-k)*(m-j)).
EXAMPLE
Given the g.f. of A323557, G(x) = Sum_{n>=0} x^n * (1 + x^n)^n / (1 + x^(n+1))^(n+1), i.e.,
G(x) = 1/(1 + x) + x*(1 + x)/(1 + x^2)^2 + x^2*(1 + x^2)^2/(1 + x^3)^3 + x^3*(1 + x^3)^3/(1 + x^4)^4 + x^4*(1 + x^4)^4/(1 + x^5)^5 + x^5*(1 + x^5)^5/(1 + x^6)^6 + x^6*(1 + x^6)^6/(1 + x^7)^7 + x^7*(1 + x^7)^7/(1 + x^8)^8 + ...
and writing G(x) as a power series in x starting as
G(x) = 1 + 3*x^2 - 2*x^3 + 2*x^4 + 9*x^6 - 14*x^7 + 8*x^8 + 12*x^10 - 12*x^11 + 15*x^12 - 52*x^13 + 76*x^14 - 36*x^15 + 2*x^16 + 50*x^18 - 104*x^19 + 79*x^20 - 140*x^21 + 324*x^22 - 276*x^23 + 128*x^24 - 144*x^25 + 118*x^26 - 28*x^27 + 72*x^28 - 336*x^29 + 657*x^30 + ... + A323557(n)*x^n + ...
then the odd coefficients of x^n in G(x), occurring at n = k*(k+1) for k>=0, form this sequence.
RELATED SEQUENCE.
The terms A323557(n*(n+3)) begin
[1, 2, 12, 50, 72, 142, 5346, 38338, 240902, 1369462, ..., A323678(n), ...].
MATHEMATICA
nmax = 20; s = CoefficientList[Series[Sum[x^m*(1 + x^m)^m/(1 + x^(m + 1))^(m + 1), {m, 0, nmax*(nmax + 1)}], {x, 0, nmax*(nmax + 1)}], x]; Table[s[[n*(n + 1) + 1]], {n, 0, nmax}] (* Vaclav Kotesovec, Jul 08 2019 *)
PROG
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Feb 03 2019
STATUS
approved