Let
- $B_n$ be the $n$-th Bernoulli number with $B_1 = -\frac{1}{2}$.
- $T_m(n,k)$ be the family of integer coefficients such that $$T_m(n,k) = \sum\limits_{j=k}^{n} {n \brace j} {j \brack k} j^m. $$
- $S_1(n,m)$ be an integer coefficients such that $$ S_1(n,m) = \sum\limits_{k=1}^{n} T_m(n,k). $$
- $S_2(n,m)$ be an integer coefficients such that $$ S_2(n,m) = \sum\limits_{k=1}^{n} T_m(n,k) B_{k-1}. $$
I conjecture that $S_1(n,m) = S_2(n,m+2)$ for all $n \in \mathbb{N}, m \in \mathbb{Z}$.
Here is the PARI/GP program to check it numerically:
T(n,k,m) = sum(j=k, n, stirling(n, j, 2)*abs(stirling(j, k, 1))*j^m)
S1(n,m) = sum(k=1, n, T(n,k,m))
S2(n,m) = sum(k=1, n, T(n,k,m)*bernfrac(k-1))
test(n,m) = S1(n,m) == S2(n,m+2)
Is there a way to prove it?