One has the nice identities
$${xy\choose 1}={x\choose 1}{y\choose 1},$$
$${xy+1\choose 2}={x+1\choose 2}{y+1\choose 2}+{x\choose 2}{y\choose 2}$$
and
$${xy+2\choose 3}={x+2\choose 3}{y+2\choose 3}+4{x+1\choose 3}{y+1\choose 3}+{x\choose 3}{y\choose 3}.$$


(The proof is essentially trivial by interpreting ${z\choose k}$ as a polynomial of degree $k$.)


This sequence of identities stops: There seems to be no nice expression of ${xy+k-1\choose k}$ as a linear combination of
${x+k-i\choose k}{y+k-i\choose k},i=1,\ldots,k$ for $k\geq 4$.


*Is there a good reason for this breakdown?* (Probably a better question is: Is there a reason for these identities to hold for $k=2$ and, especially, for $k=3$?)



Reason for breakdown of a nice binomial identity
$\def\des{\operatorname{des}}$Let $\des(\pi)$ be the number of descents of the permutation $\pi$. Then for any permutation $\pi$ in $S\_k$, we have
\begin{equation\*}\binom{xy+k-\des(\pi)-1}{k} =\sum\_{\sigma\tau=\pi}\binom{x+k-\des(\tau)-1}{k}
\binom{y+k-\des(\sigma)-1}{k}.\tag{$\*$}\label{star}
\end{equation\*}


If $\pi$ is the identity then $\des(\pi)=0$, so the left side of \eqref{star} is $\binom{xy+k-1}{k}$ and on the right, $\tau=\sigma^{-1}$. The OP's identities correspond to the fact that if $k\le 3$ then the number of descents of $\pi$ is the same as the number of descents of $\pi^{-1}$, but this is not true for $k>3$. This also explains the occurrence of the Eulerian numbers as coefficients for $k\le3$.


As far as I know, the earliest proof of \eqref{star} is in
Bogdan Mielnik and Jerzy Plebański,
*[Combinatorial approach to Baker–Campbell–Hausdorff exponents](http://www.numdam.org/item/?id=AIHPA_1970__12_3_215_0)*,
Annales de l’I. H. P., section A, tome 12, no 3 (1970), pp. 215–254 (see equation (11.10)); further references can be found in Jason Fulman and T. Kyle Petersen,
*[Card shuffling and P-partitions](https://doi.org/10.48550/arXiv.2004.01659)*, arXiv:2004.01659 [math.CO].


