$$\sum\limits_{n = 1}^{N - 1} {{q^n}} \left( {_M^{n + K}} \right) = {q^N}\sum\limits_{g = 0}^M {{{( - 1)}^{^g}}} \frac{{\left( {_{M - g}^{N + K - 1 - g}} \right)}}{{{{(q - 1)}^{g + 1}}}} + \frac{{{q^{M - k}}}}{{{{(1 - q)}^{M + 1}}}}$$
Expand nM
$${n^M} = \sum\limits_{g = 0}^M {g!\left\{ {_g^M} \right\}} \left( {_g^n} \right) = \sum\limits_{g = 0}^M {{{( - 1)}^{M - g}}g!\left\{ {_g^M} \right\}} \left( {_g^{n + g - 1}} \right) = \sum\limits_{g = 0}^M {\left\langle {_g^M} \right\rangle } \left( {_M^{n + g}} \right)$$
By combining the two, the formulas can be obtained:
$$A_q^M = \sum\limits_{k = 0}^M {{{(1 - q)}^{M - k}}} {q^k}\left\{ {_k^M} \right\}k! = q\sum\limits_{k = 0}^M {{{(q - 1)}^{M - k}}} \left\{ {_k^M} \right\}k! = \sum\limits_{k = 0}^M {\left\langle {_k^M} \right\rangle } {q^{1 + k}}$$
$$\sum\limits_{k = 0}^M {{q^n}} {n^M} = \frac{{{q^N}}}{{{{(q - 1)}^{M + 1}}}}\sum\limits_{g = 0}^M {{{(q - 1)}^g}} {( - 1)^{M - g}}\left( {_g^M} \right)A_q^{M - g}{N^g} + \frac{{A_q^M}}{{{{(1 - q)}^{M + 1}}}}$$
$$|q|<1,\mathop {\lim }\limits_{N \to \infty } \sum\limits_{n = 0}^{N - 1} {{q^n}} {n^M} = \frac{{A_q^M}}{{{{(1 - q)}^{M + 1}}}} \to A_t^M = tA(M,t)$$
A(M,t) is Eulerian polynomials.See Further Study of the Shape of the Numbers and More Calculation Formulas
More progress at https://www.preprints.org/manuscript/202305.0311/v1