OFFSET
0,6
COMMENTS
FORMULA
G.f. A(x) satisfies: A(x)^4 (mod 8) = g.f. of A083954.
EXAMPLE
A(x) = 1 + x + x^2 + x^3 - x^4 + 2*x^5 + x^7 - 2*x^8 + 5*x^9 +...
A(x)^2 = 1 + 2*x + 3*x^2 + 4*x^3 + x^4 + 4*x^5 + 3*x^6 +...
A(x)^4 = 1 + 4*x + 10*x^2 + 20*x^3 + 27*x^4 + 36*x^5 +...
A(x)^4 (mod 8) = 1 + 4*x + 2*x^2 + 4*x^3 + 3*x^4 + 4*x^5 +...
G(x) = 1 + 4*x + 2*x^2 + 4*x^3 + 3*x^4 + 4*x^5 + 4*x^6 +...
where G(x) is the g.f. of A083954.
PROG
(PARI) {a(n)=local(d=2, m=4, A=1+m*x); for(j=2, d*n, for(k=1, m, t=polcoeff((A+k*x^j+x*O(x^j))^(1/m), j); if(denominator(t)==1, A=A+k*x^j; break))); polcoeff(Ser(vector(n+1, i, polcoeff(A, d*(i-1))))^(1/2), n)}
CROSSREFS
KEYWORD
sign
AUTHOR
Paul D. Hanna, Sep 14 2005
STATUS
approved