|
PROG
|
(PARI) {a(n)=polcoeff(serreverse(x-sum(k=1, n, x^prime(k))+x*O(x^n)), n)}
for(n=1, 25, print1(a(n), ", "))
(PARI) {Dx(n, F)=local(G=F); for(i=1, n, G=deriv(G)); G}
{a(n)=local(A=x, P=sum(m=1, n, x^prime(m))+x*O(x^n)); A=x+sum(m=1, n, Dx(m-1, P^m/m!)+x*O(x^n)); polcoeff(A, n)}
for(n=1, 25, print1(a(n), ", "))
(PARI) {Dx(n, F)=local(G=F); for(i=1, n, G=deriv(G)); G}
{a(n)=local(A=x, P=sum(m=1, n, x^prime(m))+x*O(x^n)); A=x*exp(sum(m=1, n, Dx(m-1, P^m/x/m!)+x*O(x^n))); polcoeff(A, n)}
for(n=1, 25, print1(a(n), ", "))
|