OFFSET
0,2
COMMENTS
Compare to the curious identity: Sum_{n=-oo..+oo} x^n * (1 - x^n)^n = 0.
More generally, for all k we have the identity:
Sum_{n=-oo..+oo} x^n * (1 - k^n*x^n)^n = (-1) * Sum_{n=-oo..+oo} k*(k*x)^n * (1 - k*(k*x)^n)^n. - Paul D. Hanna, Dec 25 2015
LINKS
Paul D. Hanna, Table of n, a(n) for n = 0..500
FORMULA
G.f.: (-1) * Sum_{n=-oo..+oo} 2*(2*x)^n * (1 - 2*(2*x)^n)^n. - Paul D. Hanna, Dec 25 2015
It appears that for prime p >= 3, a(p) = 1 - 2^(p+1). - Peter Bala, Aug 06 2023
EXAMPLE
G.f.: A(x) = -1 - 3*x + 7*x^2 - 15*x^3 + 89*x^4 - 63*x^5 + 121*x^6 - 255*x^7 + 3521*x^8 - 13119*x^9 + 18273*x^10 - 4095*x^11 - 40319*x^12 + ...
PROG
(PARI) {a(n) = local(A=1); A = sum(k=-sqrtint(n)-1, n+1, x^k*(1 - 2^k*x^k + x*O(x^n))^k ); polcoeff(A, n)}
for(n=0, 60, print1(a(n), ", "))
CROSSREFS
KEYWORD
sign
AUTHOR
Paul D. Hanna, Nov 06 2015
STATUS
approved