(PARI) {a(n)=polcoeff(sum(m=0, n, x^m/(1-(m+1)*x+x*O(x^n))), n)} /*
Paul D. Hanna, Sep 13 2011 */
(PARI) {INTEGRATE(n, F)=local(G=F); for(i=1, n, G=intformal(G)); G}
{a(n)=local(A=1+x); A=sum(k=0, n, INTEGRATE(k, exp((k+1)*x+x*O(x^n)))); n!*polcoeff(A, n)} \\
Paul D. Hanna, Dec 28 2013
for(n=0, 30, print1(a(n), ", "))
(PARI)
{a(n)=polcoeff( sum(m=0, n, m!*x^m/(1-x +x*O(x^n))^(m+1)/prod(k=1, m, 1+k*x +x*O(x^n))), n)} /* From o.g.f. (
Paul D. Hanna, Jul 20 2014) */
for(n=0, 25, print1(a(n), ", "))
(Haskell)
a026898 n = sum $ zipWith (^) [n + 1, n .. 1] [0 ..]
(Magma) [(&+[(n-k+1)^k: k in [0..n]]): n in [0..50]]; //
Stefano Spezia, Jan 09 2019
(SageMath) [sum((n-j+1)^j for j in (0..n)) for n in (0..30)] #
G. C. Greubel, Jun 15 2021