Consider a recent arXiv preprint 1805.11445. The author of 1805.11445 has done an overview of classical problem of simplifying of power sum

$$\sum_{1\leq k\leq n}k^m, \ (n,m)\geq 0, \ m=\mathrm{const}$$

and proposed "A Fast Algorithm to Calculate Power Sum of Natural Numbers". The follow identites are found in 1805.11445 for $\Sigma_{1\leq k\leq n} k^4$ and $\Sigma_{1\leq k\leq n} k^{12}$,

$$\sum_{1\leq k\leq n} k^4=\binom{n+4}{n-1}+11\binom{n+3}{n-2}+11\binom{n+2}{n-3}+\binom{n+1}{n-4}$$

And, respectively

\begin{eqnarray*} % \nonumber to remove numbering (before each equation) \sum_{1\leq k\leq n} k^{12} &=& \binom{n+12}{n-1}+4083\binom{n+11}{n-2}+478271\binom{n+10}{n-3} \\ &+&10187685\binom{n+9}{n-4}+66318474\binom{n+8}{n-5}+162512286\binom{n+7}{n-6}\\ &+&162512286\binom{n+6}{n-7}+66318474\binom{n+5}{n-8}+10187685\binom{n+4}{n-9}\\ &+&478271\binom{n+3}{n-10}+4083\binom{n+2}{n-11}+\binom{n+1}{n-12} \end{eqnarray*}

Now we have to mention that the coefficients in corresponding sums $\Sigma_{1\leq k\leq n} k^4$ and $\Sigma_{1\leq k\leq n} k^{12}$ are terms of forth and twelfth rows of Eulerian Triangle. Therefore, we have right to generalize above identities as follows

$$(\star)\quad\sum_{1\leq k\leq n} k^{m}=\sum_{k\geq0}E_{m,k}\binom{n+m-k}{n-1-k}$$

that is analog of Faulhaber's formula for every positive integers $(n,m)$.

**Question**: Is there any reference to equation $(\star)$ ?