A243919 G.f.: Sum_{n>=0} (x^n - 1)^n * x^n / (1-x)^(n+1).

1, 0, 1, 2, 1, -2, -3, 1, 9, 13, 6, -13, -33, -39, -18, 33, 102, 159, 158, 46, -199, -525, -781, -744, -221, 808, 2089, 3101, 3212, 1933, -831, -4650, -8635, -11669, -12652, -10615, -4635, 6314, 23181, 45855, 71241, 91000, 90501, 51044, -43113, -193400, -374830, -529580, -573509

Offset

0,4

Comments

This sequence has a complex pattern of signs.

Links

Paul D. Hanna, Table of n, a(n) for n = 0..1000

Formula

G.f.: Sum_{n>=0} x^(n*(n+1)) / (1-x + x^(n+1))^(n+1).

G.f.: Sum_{n>=0} x^n * Sum_{k=0..n} binomial(n,k) * (x^k - 1)^k.

G.f.: Sum_{n>=0} x^n * Sum_{k=0..n} binomial(n,k) * (1 - x^k)^(n-k) * x^(k^2).

Example

G.f.: A(x) = 1 + x^2 + 2*x^3 + x^4 - 2*x^5 - 3*x^6 + x^7 + 9*x^8 + 13*x^9 +...
where we have the series identity:
A(x) = 1/(1-x) + (x-1)*x/(1-x)^2 + (x^2-1)^2*x^2/(1-x)^3 + (x^3-1)^3*x^3/(1-x)^4 + (x^4-1)^4*x^4/(1-x)^5 +...+ (x^n-1)^n * x^n / (1-x)^(n+1) +...
A(x) = 1 + x^2/(1-x+x^2)^2 + x^6/(1-x+x^3)^3 + x^12/(1-x+x^4)^4 + x^20/(1-x+x^5)^5 + x^30/(1-x+x^6)^6 +...+ x^(n*(n+1)) / (1-x + x^(n+1))^(n+1) +...
as well as the binomial identity:
A(x) = 1 + x*(1 + (x-1)) + x^2*(1 + 2*(x-1) + (x^2-1)^2) + x^3*(1 + 3*(x-1) + 3*(x^2-1)^2 + (x^3-1)^3) + x^4*(1 + 4*(x-1) + 6*(x^2-1)^2 + 4*(x^3-1)^3 + (x^4-1)^4) + x^5*(1 + 5*(x-1) + 10*(x^2-1)^2 + 10*(x^3-1)^3 + 5*(x^4-1)^4 + (x^5-1)^5) +...+ x^n * Sum_{k=0..n} binomial(n,k) * (x^k-1)^k +...
A(x) = 1 + x*(0 + x) + x^2*(0 + 2*(1-x)*x + x^4) + x^3*(0 + 3*(1-x)^2*x + 3*(1-x^2)*x^4 + x^9) + x^4*(0 + 4*(1-x)^3*x + 6*(1-x^2)^2*x^4 + 4*(1-x^3)*x^9 + x^16) + x^5*(0 + 5*(1-x)^4*x + 10*(1-x^2)^3*x^4 + 10*(1-x^3)^2*x^9 + 5*(1-x^4)*x^16 + x^25) +...+ x^n * Sum_{k=0..n} binomial(n,k) * (1-x^k)^(n-k) * x^(k^2) +...

Prog

(PARI) {a(n)=local(A); A=sum(m=0, n, (x^m - 1)^m * x^m / (1-x +x*O(x^n) )^(m+1) ); polcoeff(A, n)}
for(n=0, 60, print1(a(n), ", "))
(PARI) {a(n)=local(A); A=sum(m=0, sqrtint(n+1), x^(m*(m+1)) / (1-x + x^(m+1) +x*O(x^n) )^(m+1) ); polcoeff(A, n)}
for(n=0, 60, print1(a(n), ", "))
(PARI) {a(n)=polcoeff(sum(m=0, n, x^m*sum(k=0, m, binomial(m, k)*(x^k-1)^k) +x*O(x^n)), n)}
for(n=0, 40, print1(a(n), ", "))
(PARI) {a(n)=polcoeff(sum(m=0, n, x^m*sum(k=0, m, binomial(m, k)*(1-x^k)^(m-k)*x^(k^2)) +x*O(x^n)), n)}
for(n=0, 40, print1(a(n), ", "))

Crossrefs

Cf. A243988.

Sequence in context: A232177 A111377 A014046 * A231775 A372840 A128065

Adjacent sequences: A243916 A243917 A243918 * A243920 A243921 A243922

Keyword

sign

Author

Paul D. Hanna, Jun 18 2014

Extensions

Name changed by Paul D. Hanna, Jul 02 2014

Status

approved