m = 3; e[n_, 0, m_] := 1;
e[n_, k_, m_] := 0 /; k >= n;
e[n_, k_, 1] := 1 /; k >= n;
e[n_, k_, m_] := (k + m - 1)e[n - 1, k, m] + (m*n - k + 1 - m)e[n - 1, k - 1, m];
Table[Table[e[n, k, m], {k, 0, n - 1}], {n, 1, 10}];
Table[Table[e[n, k, m] + e[n, n - k, m], {k, 0, n}], {n, 0, 10}];
Flatten[%]