OFFSET
0,3
FORMULA
a(n) = 2*(4^(n-1) - 1)*a(n-1) - (-1)^n for n>0 with a(0)=1.
PROG
(PARI) {a(n)=if(n==0, 1, 2*(4^(n-1)-1)*a(n-1) - (-1)^n)}
(PARI) {a(n)=local(A=x); A=sum(m=0, n, x^m*prod(j=1, m, (2^(2*j-1)-1)/(1+x*2^(2*j-1))+x*O(x^n))); polcoeff(A, n, x)}
for(n=0, 20, print1(a(n), ", "));
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Jul 26 2012
STATUS
approved