|
EXAMPLE
|
From the table of powers of A(x) (A100224), we see that
2^n-1 = Sum of coefficients [x^0] through [x^n] in A(x)^n:
A^1=[1,0],1,1,0,-2,-3,1,11,...
A^2=[1,0,2],2,1,-2,-5,-2,12,...
A^3=[1,0,3,3],3,0,-5,-6,6,...
A^4=[1,0,4,4,6],4,-2,-8,-3,...
A^5=[1,0,5,5,10,10],5,-5,-10,...
A^6=[1,0,6,6,15,18,17],6,-9,...
A^7=[1,0,7,7,21,28,35,28],7,...
A^8=[1,0,8,8,28,40,60,64,46],...
the main diagonal of which is A001610 = [0,2,3,6,10,17,...],
where Sum_{n>=1} A001610(n-1)/n*x^n = log((1-x)/(1-x-x^2)).
|
|
PROGRAM
|
(PARI) {a(n)=if(n==0, 1, (2^n-1-sum(k=0, n, polcoeff(sum(j=0, min(k, n-1), a(j)*x^j)^n+x*O(x^k), k)))/n)} (PARI) {a(n)=if(n==0, 1, if(n==1, 0, if(n==2, 1, ((2*n-3)*a(n-1)-5*(n-3)*a(n-2))/n)))} (PARI) {a(n)=polcoeff((1+x+sqrt(1-2*x+5*x^2+x^2*O(x^n)))/2, n)}
|