|
PROG
|
(PARI) {a(n)=local(k=0, A=vector(n+1), p); A[1]=1; for(j=1, n-k-1, p=(n-1)^2-(n-j-1)^2; A=Vec((Polrev(A)+x*O(x^p))/(1-x))); A=Vec((Polrev(A)+x*O(x^p))/(1-x)); A[p+1]}
for(n=0, 25, print1(a(n), ", "))
(PARI) {a(n)=local(A=[1]); for(i=1, n, A=Vec(Ser(A)/(1-x)^(2*(#A)-3)); A=concat(A, A[#A])); A[#A]}
for(n=0, 25, print1(a(n), ", "))
|