|
MATHEMATICA
|
cols = Tuples[{0, 1, 2, 3}, 4];
tmat = Table[If[Or @@ MapThread[SameQ, cols[[{i, j}]]], 0, 1], {i, 256}, {j, 256}];
vec = vvec = ((x^Count[#, 0] * y^Count[#, 1] * z^Count[#, 2]) & /@ cols);
Prepend[Table[vec = Expand[vvec*(tmat.vec)]; Coefficient[Total[vec], (x*y*z)^n]/24, {n, 2, 10}], 1]
|