I got this general formula for $ n\in N$ (I showed it here)
$$\int_0^1 \left(\frac{x}{1-x} \ln x \right)^n dx=n \sum_{p=0}^{n-1}a(n,p+1) (-1)^{n-p} \zeta(p+2)+n! $$ where $a(n,k)$ is the coefficient of $x^{k-1}$ in the expression $$\prod_{p=2}^n (1+xp)$$
I found in A145324 this formula using Stirling numbers of first kind $$ a(n,k)=\sum_{p=0}^{k-1} (-1)^p |s_1(n+1,n-p+1)|$$
MY QUESTION
how can I prove the last formula that I found in OEIS and is there a closed form for it without using series ?