from sage.all import QQ, PolynomialRing, SymmetricFunctions
T = PolynomialRing(QQ, 't')
P = SymmetricFunctions(T.fraction_field()).hall_littlewood().P()
grouped = P([3, 2]).expand(4, alphabet=['x1', 'x2', 'x3', 'x4'])
# Sage groups coefficients in t; this writes the same polynomial
# expanded over QQ[x1, x2, x3, x4, t], as the table does.
R = PolynomialRing(QQ, ['x1', 'x2', 'x3', 'x4', 't'])
print(R(str(grouped)))Every entry is an exact polynomial with integer coefficients. The generator computes Macdonald's symmetrization formula in the fraction field of $\mathbb{Z}[x_1,\dots,x_n,t]$ and refuses an entry unless the denominator cancels. Before the values were written, every entry was compared with Sage's Hall-Littlewood $P$ basis [3].
The same values were checked at $t=0$ against the Schur polynomials, at $t=1$ against the monomial symmetric polynomials, and in the transition identity with the Kostka-Foulkes polynomials.