JavaScript (ES6), 7877 bytes
a=>a.map((_,k)=>a.reduce((t,v,n)=>t+v*(S=k=>ns=n=>n?s(--n+k*S(n,k--)+S+ ++k*s(k-1n)-n++:!k)(kn),0))
Try it online!Try it online!
a => // a[] = coefficients in ascending order
a.map((_, k) => // for each coefficient at index k in a[]:
a.reduce( // for each coefficient v at index n in a[],
(t, v, n) => // with t as the accumulator:
t + // add to t ...
v * ( // ... v multiplied by S(n, k)
// actually implemented as Ss(kn) where Ss is ...
Ss = kn => // ... a recursive function taking only kn explicitly:
n ? // if n is not 0:
s(--n // decrement n
+, k * S(k--) // compute k * S(k)
+ S(kn - 1) // +, S(k - 1)
- n++ + ++k * s(n) // restore n+ (andk cancel* outS(n the- earlier1, --nk)
: // else:
!k // return 1 if k = 0, or 0 otherwise
)(kn), // initial call to S
0 // start with t = 0
) // end of reduce()
) // end of map()