p[n_] := Product[3^k-1, {k, 1, n}]; t[n_, m_] := p[n]/(p[n-m]*p[m]); Table[t[n, m], {n, 0, 10}, {m, 0, n}] // Flatten (*
Jean-François Alcover, Jan 10 2014, after
R. J. Mathar *)
Table[QBinomial[n, k, 3], {n, 0, 10}, {k, 0, n}] // Flatten
S[n_, k_, q_] /; 1 <= k <= n := S[n-1, k-1, q] + Sum[q^j, {j, 0, k-1}]* S[n-1, k, q]; S[n_, 0, _] := KroneckerDelta[n, 0]; S[0, k_, _] := KroneckerDelta[0, k]; S[_, _, _] = 0;
T[n_, k_] /; n >= k := Sum[Binomial[n, j]*S[n-j, n-k, q]*(q-1)^(k-j) /. q -> 3, {j, 0, k}];