Suppose I have two independent integer random variables $X_1$, $X_2$ (with constraint that $X_1+X_2\le N,0\le X_1\le N,0\le X_2\le N$), with probability generating functions $g_1(z)$, $g_2(z)$. Now I have a joint-distribution $P(X_1-X_2,2X_1+X_2-N)$, whose probability generating function is $G(z_1,z_2)$. $N$ is a constant.
What is the relation between $G(z_1,z_2)$ and $g_1(z)$, $g_2(z)$?
The following question is related but not the same: multivariate probability generating function