, j n, which may be arranged into an n_n matrix array Z=(zij). These polynomials are indexed by double Gelfand patterns, or equivalently, by pairs of column strict Young tableaux of the same shape. Using the double labeling property, one may define a square matrix D(Z), whose elements are the double-indexed polynomials. These matrices possess the remarkable "group multiplication property" D(XY)=D(X) D(Y) for arbitrary matrices X and Y, even though these matrices may be singular. For Z=U # U(n), these matrices give irreducible unitary representations of U(n). These results are known, but not always fully proved from the extensive physics literature on representation of the unitary groups, where they are often formulated in terms of the boson calculus, and the multiplication property is unrecognized. The generality of the multiplication property is the key to under-standing group representation theory from the purview of combinatorics. The combinatorial structure of the general polynomials is expected to be intricate, and in this paper, we take the first step to explore the combinatorial aspects of a special class which can be defined in terms of the set of integral matrices with given row and column sums. These special polynomials are denoted by LXβ (Z), where ɑ and β are integral vectors representing the row sums and column sums of a class of integral matrices. We present a combinatorial interpretation of the multiplicative properties of these polynomials. We also point out the connections with MacMahon's Master Theorem and Schwinger's inner product formula, which is essentially equiv-alent to MacMahon's Master Theorem. Finally, we give a formula for the double Pfaffian, which is crucial in the studies of the generating function of the 3n-j coef-ficients in angular momentum theory. We also review the background of the general polynomials and give some of their properties.